人工源宽角地震实验的谱元法模拟和逆时偏移成像方法研究

人工源宽角反射(折射)地震资料具有偏移距较大、信噪比较低等特点,通常用于地震波走时反演重建地壳速度结构。逆时偏移成像方法作为勘探地震学领域获取地下构造形态的有效手段之一,可以有效弥补走时反演方法的不足。本文针对大偏移距宽角反射(折射)地震实验,利用四边形网格谱元法进行波场模拟,结合了有限元法的灵活性和谱方法的指数收敛性,高效且高精度获取模型合成地震记录,后采用逆时偏移成像方法将合成地震记录偏移归位,获取地壳几何结构,验证了逆时偏移成像方法在宽角地震资料处理及结果解释中的适用性,为后期实际地震资料的偏移成像提供了理论依据和支持。     

不同产状断层错动的地表重力变化和形变

根据点源位错引起的地表重力变化和形变的理论公式,应用数值算法求得半无限空间任意产状断层位错在地表的重力变化和形变．计算结果表明,产状和位错方式是决定断层的地表重力和形变效应的基本因素．铲式断层活动的地表重力和形变效应与矩形断层相比有显著差异．断层位错引起的重力变化、垂直形变、视垂直形变的图象特征相似．根据重力变化和形变的特殊图象特征,有可能从观测资料中提取断层活动引起的重力变化和形变信息,判识和估计发震断层的几何学和运动学特征．      

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

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

钦杭—武夷山成矿带上地壳速度结构与基底特征:万载—惠安宽角反射/折射地震剖面约束

钦杭成矿带和武夷山成矿带是华南大陆两个重要的成矿带,成矿作用主要发生于中生代陆内造山时期.地质研究表明,基底和地表断裂的特征对成矿过程有重要的控制作用,研究上地壳结构特征对成矿差异性特征的认识有重要的参考价值.为此,本文基于跨越钦杭、武夷山成矿带江西万载至福建惠安的NW-SE向深地震测深剖面初至波数据,利用有限差分走时反演方法,获得了钦杭、武夷山成矿带8 km深度范围内的上地壳P波速度结构,其主要特征为：（1）钦杭、武夷山成矿带上地壳P波速度横向非均匀特征明显,以5.8 km·s^-1速度等值线作为基底参考面,发现剖面基底埋深较浅,约1.0~3.0 km;钦杭成矿带的基底埋深总体小于武夷山成矿带,分别为0.5~2.0 km和1.5~3.0 km;（2）P波高速区（速度正异常区）与地表出露的岩浆岩对应较一致,P波低速区（速度负异常区）与主要的断裂位置或沉积盆地对应较一致,绍兴—江山—萍乡断裂和政和—大浦断裂下方的低速特征显示两条断裂至少向深部延伸8 km以上,暗示两条断裂具备深大断裂的性质,推测绍兴—江山—萍乡断裂可能是扬子块体和华夏块体的边界;（3）综合已有的地质、地球物理资料,我们推测钦杭成矿带和武夷山成矿带上地壳P波速度的不同,反映了深部岩浆作用过程的差异,基底深度及断裂性质是造成两个成矿带成矿差异的重要因素.     

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

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

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

为了理解长江中下游地区在中生代成矿的深部动力学过程,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面上隆、下地壳及上地幔的低速异常等壳幔结构特征,预示下地壳不以榴辉岩残体为主,支持燕山期地幔岩浆的上涌和侵入并成矿,是热上涌物质的源地。     

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…     

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

地震属性分析技术在地球物理勘探领域的广泛应用,启发研究人员将其应用于人工源宽角反射/折射深地震测深剖面的资料预处理和震相识别。采用札达-泉水沟深地震测深资料,提取振幅、信噪比、主频、瞬时带宽、瞬时高频能量等地震属性参数,分析不同参数的物理含义,挑选其中对界面变化敏感的参数,对深地震测深资料进行预处理,并利用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°,为一次高倾角走滑型地震.两种方法获得的余震震源机制解比较一致,验证了结果的可靠性.研究结果显示,鲁甸地震的高倾角破裂特征,会导致其能量释放快速且完整,这可能是导致鲁甸地震地面破坏加重及缺乏较大震级余震的主要原因.