利用背景噪声成像研究合肥-金华剖面地壳速度结构及径向各向异性的东西差异

背景噪声成像能够为地壳内部结构提供地震波速的观测证据,能够增进对地壳结构与物性的认识.本文利用背景噪声成像方法对合肥-金华流动地震台阵55个台站数据进行了分析反演,获得了华北克拉通与华南块体东部构造边界附近区域的地壳S波速度与径向各向异性结构.成像结果表明,以东经118°-118.5°附近为界,西北部表现为地壳低速异常和中下地壳正的径向各向异性,东南部表现为地壳高速异常,下部地壳存在负的径向各向异性.由此,本文推测地震波速结构差异反映了地壳内部物质与温度的差异,这种差异与两个地区自中生代晚期以来经历了截然不同的岩浆活动相关.     

来自海底高速层径向波的理论地震图研究

本文利用各向异性反射率技术计算理论地震图,提出海底高速薄层会产生沿高速层水平传播的波（简称径向波）,这种波在水层中作为P波,在固液界面激发下行横波,该均匀横波以临界角入射高速薄层,在层内作为超临界角的非均匀横波水平传播,再以临界角转换为上行传播的均匀横波,最终在固液界面上行透射转换为水层中P波.高速薄层传播的径向波不同于界面折射波,也不同于具有频散的面波和通道波.理论地震图的研究表明,径向波具有线性时距,能与海底强反射具有同等振幅水平;径向波有其振幅、时距位置和斜率这些观测记录参数,分别对应高速层的厚度、深度和近似的横波速度;径向波可以克服折射波解释中遇到的振幅强弱和高速层速度等困难.径向波可作为探测海底高速薄层的有力工具,对于研究高速层屏蔽、海底反射类型的多样性和相应的资料处理解释有重要意义.     

基于交错网格有限差分法的断面波叠前时间偏移成像影响因素正演分析

断层和断裂带的有效识别是地震资料解释中的重要环节,断层在地震信号响应中以断面波的形式体现,因此断面波成像的质量关系到断层的精细识别与刻画.本文利用精度较高的交错网格有限差分正演模拟方法对断面波成像的影响因素进行了正演研究,主要正演分析的参数包括采集因素中的电缆长度和采集方向,地质因素中的断层倾角、断距、反射系数,以及处理因素中的偏移方法等几个方面.通过正演论证得出：采用合理的采集参数能够提高断面波的照明度;有效结合地质因素能够提高断面波的解释精度;利用合理的偏移方法能够使断层归位更加准确,断面波有效成像.基于以上结论,对于断面波的精确识别与刻画,应综合采集因素,处理因素及地质因素,只有这样才能提高断层的解释精度,有效减小解释误差.     

含孔隙介质的分层半空间表面瑞利波的衰减特性

以分层半空间内部含有一层孔隙介质为物理模型进行数值计算,研究半空间表面瑞利波的传播和衰减特性.为更加接近实际,结合瑞利波的激发特性,确定了瑞利波的主衰减曲线,并主要以此进行规律分析.针对速度递增和含低速层这两种典型的地质模型,讨论了瑞利波的传播衰减在不同地质模型下的特性,并分析了各自的规律.结果表明,在这两种模型下瑞利波的主衰减曲线都受孔隙介质所处空间位置影响产生比较明显的变化,但衰减系数极大值对应的波长与模型的表层厚度存在较明显的线性对应关系,利用这一关系,可以在实际勘探中快速得到表层介质厚度.另外,通过对比分析还可以看到,瑞利波主衰减曲线随孔隙介质的孔隙度和渗透率的变化都强于主频散曲线的变化,表明衰减曲线对孔隙度和渗透率的变化更加敏感,理论上更加适合进行介质参数反演工作.综合对比结果,我们认为瑞利波主衰减曲线中包含了更丰富的介质参数信息,如果能够有效利用,将可以提高瑞利波勘探的准确性和应用范围.     

利用P波接收函数研究青藏高原东南缘地幔转换带结构

研究青藏高原东南缘的深部结构对于理解印度-欧亚板块的碰撞机理和青藏高原的形成演化具有重要的科学意义.本研究对布设在研究区域内566个固定和流动地震台站的波形资料进行了处理,获得77853条高质量P波接收函数,应用接收函数共转换点（CCP）叠加技术获得了研究区域下方精细的地幔转换带间断面起伏形态及转换带厚度变化图像.结果表明：研究区域南北方向上具有两个明显的转换带增厚异常区,南侧异常区位于滇中次级块体与印支块体下方,可能是新特提斯洋板片与上部印度板块间断离并部分滞留在转换带底部的结果;北侧川西地区异常增厚可能与上方岩石圈拆沉并降至转换带有关;腾冲火山起源可能是板块俯冲过程中发生断离造成软流圈物质部分熔融,湿热物质上涌所致.     

基于广义旋转法双九点格式标量波方程数值模拟

前人虽然基于传统旋转法提出了四阶精度最优17点差分格式,用于提高频率域地震波场数值模拟精度,但其仅适用于等网格间距的情形.这大大限制了该格式的使用范围.为了进一步提高17点有限差分格式的数值精度并将其推广到矩形网格,本文基于二阶精度有限差分算子利用广义旋转法提出了双九点格式,由于其差分格点与前人17点格式在分布上一样,所以也可称为二阶精度最优17点格式,但由差分格式构造原理上来讲,称其为双九点格式更妥.前人四阶精度17点格式仅为本文等网格间距情形时的特殊情况,本文方法单位波长网格点数仅需要2.2个即可.本文格式在继承传统旋转法良好几何旋转性质的同时,拥有平均导数方法适用于矩形网格的特点,和平均导数法所得的广义17点格式相比,本文格式数值精度更高,数值频散抑制性能和差分算子对称性更好.同时,本文双九点格式方法和思想对于后人借助传统九点格式的构造方法将其扩展到17点格式求解各类波动方程具有十分重要的意义.     

自适应虚同相轴方法压制地震层间多次波

地震数据中发育的层间多次波是影响速度分析和偏移成像的精度和可靠性的关键.通常情况下,层间多次波的动校正量、叠加速度和频率与一次波并无明显差异,从而对识别、预测和压制多次波带来了极大挑战.传统虚同相轴方法基于物理图像和定性公式,其预测的层间多次波振幅和相位精度难以满足实际需求,造成了其对匹配算法的过度依赖.本文针对传统虚同相轴方法的理论缺陷和计算精度问题,通过理论推导得到了新的自适应虚同相轴方法.相比于传统方法,自适应虚同相轴方法能够显著提高压制多次波能力,同时减少对匹配算法的依赖.本文给出了自适应虚同相轴方法的推导过程,并运用一维和二维模型算例验证了方法相较于传统虚同相轴方法的多次波预测精度优势.通过在PLUTO模型和实际陆地地震数据上的应用实例,证明了本文新研究的自适应虚同相轴方法对去除层间多次波,恢复并突出目标储层同相轴,提高地震成像分辨率的显著作用.     

基于背景噪声研究青藏高原东北缘瑞利波相速度和方位各向异性

本研究收集了甘肃、青海、宁夏等118个宽频带数字地震台站的连续波形资料,利用噪声互相关,经过计算和筛选,在5~38 s范围内,共得到5773条瑞利波相速度频散曲线.然后采用1°×1°的网格划分,反演获得青藏高原东北缘相速度和方位各向异性分布.结果表明：短周期8~12 s内,鄂尔多斯从低速异常变为高速异常;该周期范围内各向异性结果与区域断裂走向有很好的一致性.18~25 s周期内,祁连地块、松潘-甘孜地块、羌塘地块低速异常范围逐渐变大,随周期增加地壳低速异常与人工探测结果相符;鄂尔多斯表现为速度随周期增加逐渐变大,说明其中下地壳速度相对偏高,不存在低速异常;该周期范围内的各向异性特征表现为,祁连地块和松潘甘孜地块大致呈NW-SE方向,而青藏高原内部快波方向显示了顺时针旋转的形态.在30~35 s范围内面波速度主要受莫霍面深度和莫霍面附近介质速度的影响,与地壳厚度分布有非常好的吻合.综合不同方法获得的各向异性研究结果,支持印度-欧亚板块的碰撞使青藏高原东北缘地壳发生缩短和逐渐隆升的观点,认为整个岩石圈的垂直缩短变形是青藏高原东北缘的主要形成机制.     

各向异性介质的矢量波场分离与应用

针对裂缝各向异性介质,本文提出一种非正交假设下的矢量波场分离方法.本文首先对多分量地震勘探中常见的波型泄漏现象进行了数学描述,提出在纵、横波波场分离的同时应该考虑恢复纵、横波的矢量振幅.为了对裂缝方位角与各向异性系数进行定量预测,本文将矢量波场分离拆分成三个步骤来实施：第一步,用Z、R两分量的仿射坐标系变换分离ZR平面内的P波投影与SV波;第二步,用ZR平面内的P波投影与T分量的仿射坐标系变换分离P波与SH波;第三步,用纯净的SV波与SH波的成像剖面分离快慢横波,并预测裂缝发育参数.模型数据与实际数据的试验结果表明,本文提出的纵、横波波场分离方法能够获得完整的矢量振幅信息,并提供裂缝预测的精度.     

地震槽波的数学-物理模拟初探

针对地震槽波在低速层的传播特性,开展了煤层内地震槽波勘探的数值模拟和物理模拟研究的初探工作.在数值模拟研究方面,采用交错网格有限差分法对煤层中的地震槽波进行三分量全波场模拟.基于波场快照和人工合成地震记录研究了不同模型中的波场特征和各种波型的传播规律.在物理模拟方面,通过选用不同配比的环氧树脂和硅橡胶类材料构建地震槽波物理模型,利用透射法和反射法观测系统获得了清晰的地震槽波记录以研究槽波的地震学特征.研究表明,在煤层内槽波的地震波场中,Love型槽波的能量小于Rayleigh型槽波的SV分量,大于Rayleigh型槽波的SH分量.相对于Love型槽波和Rayleigh型槽波的SH分量,Rayleigh型槽波的SV分量在围岩中的泄露能量较强.在煤层界面附近的围岩中,地震波仍以槽波形式传播,随着距离的增加能量逐渐衰减.随着煤层变薄,煤层槽波主频向高频方向移动,频散现象增强,传播速度增大.