井中震源的远场波场特征研究

井中震源在逆VSP、随钻地震和采矿地球物理研究中都有广泛应用.满足"小井孔"(井孔半径远小于特征波长)及"远场"(炮检距大于特征波长)假设时,井中震源的远场波场存在解析解.为了检验解析解在不同情况下的适用性,本文使用最速下降积分计算了不满足上述假设时井中震源远场波场的合成地震记录,即半解析解.模型试验表明,解析解只能在同时满足"小井孔"和"远场"假设时使用;当这两个假设条件不满足时,解析解的振幅和波形相对于半解析解会有明显的偏差.随着假设不满足程度的增加,偏差会逐渐增加,并会逐渐影响走时的准确拾取;这种条件下,采用半解析解才能获得准确的井中震源波场.     

贴体网格各向异性对坐标变换法求解起伏地表下地震初至波走时的影响

笛卡尔坐标系中的经典程函方程在静校正、叠前偏移、走时反演、地震定位、层析成像等很多地球物理工作中都有应用,然而用其计算起伏地表的地震波走时却比较困难.本文通过把曲线坐标系中的矩形网格映射到笛卡尔坐标系的贴体网格,推导出曲线坐标中的程函方程,而后,用Lax-Friedrichs快速扫描算法求解曲线坐标系的程函方程.研究表明本文方法能有效处理地表起伏的情况,得到准确稳定的计算结果.由于地表起伏,导致与之拟合的贴体网格在空间上的展布呈各向异性,且这种各向异性的强弱对坐标变换法求解地震初至波的走时具有重要影响.本文研究表明,随着贴体网格的各向异性增强,用坐标变换法求解地表起伏区域的走时计算误差增大,且计算效率降低,这在实际应用具有指导意义.     

地震属性参数在札达—泉水沟深地震测深资料预处理与震相识别中的应用

地震属性分析技术在地球物理勘探领域的广泛应用，启发研究人员将其应用于人工源宽角反射/折射深地震测深剖面的资料预处理和震相识别。采用札达-泉水沟深地震测深资料，提取振幅、信噪比、主频、瞬时带宽、瞬时高频能量等地震属性参数，分析不同参数的物理含义，挑选其中对界面变化敏感的参数，对深地震测深资料进行预处理，并利用P波和S波的联合扫描，提高震相识别的准确性。走时互换结果显示，采用地震属性参数可有效提高震相拾取的准确性，进而提高后续地壳速度结构反演结果的精度。     

青藏高原深部地球物理探测70年

青藏高原是全球海拔最高、规模最大、时代最新的陆-陆碰撞造山带.几十个百万年以来高原隆升、喜马拉雅山系崛起是地球演化史上最为壮观的构造事件之一.青藏高原壳幔结构和深部过程备受国内外地球科学家关注.近七十年来的地球物理研究与探索表明:(1)青藏高原地壳巨厚,岩石圈较薄;(2)壳内存在软弱层,但厚度和联通性有限;(3)高原下地壳及Moho面广泛发育叠瓦状反射特征,存在明显的脆性变形;(4)喜马拉雅和拉萨块体南部存在双Moho现象/迹象;(5)印度大陆岩石圈向高原下方俯冲的形态存在显著的东西向差异;(6)高原主体上地幔各向异性以NEE向为主;(7)青藏高原布格重力异常四周高、中间低,异常场边界与地形梯度变化密切相关;(8)高原内部磁异常较弱,周边地区较强,其分界与区域构造边界基本一致;(9)青藏高原水热活动强烈,大地热流值高,主要来自加厚地壳的贡献.但是,有关青藏高原深部过程,诸如是否存在中/下地壳流、印度与欧亚大陆岩石圈的俯冲模式等重大科学问题目前仍存在争议.青藏高原地球物理和动力学研究是一个复杂的系统工程,以科学问题为导向,结合国家重大需求,在关键区域组织实施综合地球物理探测,可望在地学领域取得创新与突破.     

基于分频波形拟合方法的鲁甸M_S6.5地震序列震源机制解研究

2014年鲁甸M_S6.5地震造成了极大的人员伤亡和经济损失,其孕震机理和动力学成因为国内外学者广泛关注.由于目前的反演主要采用单一的反演方法进行处理,导致结果的不确定性较大,使得该地震的发震构造及孕震机理等问题仍存在较大争议.针对目前研究存在的不足,本文基于近震波形求解震源机制解的gCAP(general Cut And Paste)方法和全波形拟合方法,分别反演了鲁甸地震主震及其9次余震的震源机制解,以此评价近震数据反演震源机制解结果的可靠性与稳定性.结果表明数据方位角覆盖对gCAP解的稳定性有较大影响.针对全波形拟合方法反演过程中,低频信号稳定、高频信号解析度高的特点,采用0.01～0.05Hz和0.01～0.2Hz分频段波形拟合思路:在低频段剔除拟合差的数据,进一步在高频段进行高解析度波形拟合,从而获得主震可靠稳定的震源机制解.研究结果表明鲁甸地震主震倾角为76°～83°、滑动角为-157°～-164°,为一次高倾角走滑型地震.两种方法获得的余震震源机制解比较一致,验证了结果的可靠性.研究结果显示,鲁甸地震的高倾角破裂特征,会导致其能量释放快速且完整,这可能是导致鲁甸地震地面破坏加重及缺乏较大震级余震的主要原因.     

阿坝—遂宁宽角地震剖面重建藏东缘龙门山地区地壳速度结构

龙门山断裂带位于青藏高原东缘,在中生代和晚新生代经历强烈的构造变形,急剧抬升,是研究青藏高原隆升和扩展动力学过程的重要窗口.本文利用起伏地形下的高精度成像方法,对"阿坝一龙门山一遂宁"宽角反射/折射地震数据重新处理,通过走时反演重建研究区地壳速度结构.剖面自西向东跨越松潘一甘孜块体、龙门山断裂带和四川盆地,不同块体速度结构表现了显著的差异.松潘甘孜块体地表复理石沉积层内有高速岩体侵入,低速层低界面起伏不平反映了该区的逆冲推覆构造.中下地壳速度横向上连续变化,平均速度较低(约6.26 km·s~(-1)).四川盆地沉积层西厚东薄,并在西侧出现与挤压和剥蚀作用相关的压扭形态.中下地壳西薄东厚,平均速度较高(约6.39 km.s~(-1)).龙门山断裂带是地壳速度和厚度的陡变带,Moho面自西向东抬升约13 km.在整个剖面上Moho面表现为韧性挠曲,中下地壳横向上连续变化,推测古扬子块体已到达松潘甘孜块体下方.松潘甘孜块体下方中下地壳韧性变形,并在底部拖曳着被断裂切割的脆性上地壳,应力在不同断裂上积累和释放,诱发大量地震.     

弹性波逆时偏移中的稳定激发振幅成像条件

本文针对弹性波逆时偏移,提出稳定的激发振幅成像条件.在震源波场的正向传播过程中,计算每个网格点的能量,并保存最大能量密度的时刻和相应的波场值;在检波器波场的逆时传播过程中,在每个网格点提取最大能量密度时刻的检波器波场值,并利用保存的最大能量震源波场做归一化,获得角度依赖的反射系数成像剖面.相比于归一化互相关成像条件,该成像条件在震源波场的正向传播过程中无需存储波场快照,节省大量磁盘空间和I/O吞吐任务,提高了计算效率;相比于弹性波的激发时间成像条件,该成像条件自动校正了水平分量在震源两侧的极性反转,在多炮叠加时避免振幅损失.数值试验表明,与归一化成像条件相比,稳定激发振幅成像条件具有更小的计算量,偏移剖面的低频假象更弱,水平分量的成像能力更优,具有更高的空间分辨率.     

Gravity changes and surface deformations due to faults with different geometry

IntroductionNowadaystheprecisionofgravitymeasurementisupto10xlo--8InjsZ(Okubo,1992),theprecisionofaltimetryisoftheorderofmm,andtheprecisionofrelativemeasurementinhorizontaldirectionusingGPStechniquesreaches(0.01--0.001)x10ed(Liu,etal,1998).Thegravitychangesinducedbyalargeearthquake(e.g.,Tangshanearthquakeof1976)canbetensorevenhundredsof10--8m/s',andthemaximumdisplacementinhorizontalandverticaldirectioncreatedbyitalwaysexceedI.0m(Chen,etal,1979,1980,Zhang,etal,1981).Therefore,structuresrela…     

VTI介质起伏地表地震波场模拟

起伏地表下地震波场模拟有助于解释主动源和被动源地震探测中穿过山脉和盆地的测线所获得的资料.然而传统的有限差分法处理起伏的自由边界比较困难,为了克服这一困难,我们将笛卡尔坐标系的各向异性介质弹性波方程和自由边界条件变换到曲线坐标系中,采用一种稳定的、显式的二阶精度有限差分方法离散(曲线坐标系)VTI介质中的弹性波方程;对...     

长江中下游成矿带及邻区地壳速度结构：来自利辛-宜兴宽角地震资料的约束

为了理解长江中下游地区在中生代成矿的深部动力学过程,Sinoprobe-03-02项目于2011年9月至10月,在跨宁芜矿集区和郯庐断裂带实施了从安徽利辛至江苏宜兴450km长的宽角反射/折射地震剖面。速度剖面结果显示,Moho面深度和地壳速度结构在郯庐断裂两侧东西方向存在明显的差异:(1)在东部扬子块体内部,地壳覆盖层厚3~5km,西部的合肥盆地下方,则达到4~7km。(2)剖面平均Moho面深度为30~32km左右,在郯庐断裂下方,Moho面深度在35km左右;在宁芜矿集区下方,Moho面整体深度偏浅,达30~31km左右,但局部范围内,Moho面深度至34km左右。(3)剖面的下地壳平均速度在6.5~6.6km/s左右,在宁芜矿集区下方,下地壳速度偏低,为6.4~6.5km/s左右。剖面上地幔顶部的速度结构平均在8.0~8.2km/s。在宁芜矿集区下方,速度偏低,为7.9~8.1km/s左右。(4)郯庐断裂带的下方,从地表开始,还存在20多千米长的低速异常带,一直延伸到Moho面附近。剖面的宁芜矿集区下方Moho面上隆、下地壳及上地幔的低速异常等壳幔结构特征,预示下地壳不以榴辉岩残体为主,支持燕山期地幔岩浆的上涌和侵入并成矿,是热上涌物质的源地。