【摘 要】
：
本文针对地震波速度应力方程，构造了求解方法，并通过理论分析和大量的数值实验表明，该方法具有较好的时空稳定性以及较高的波动方程求解精度效率。本文将优化的辛格式交错网格褶积微分算子作为地震波逆时偏移算子，对叠前资料作振幅补偿处理，以地面速度分量记录作为边值条件进行波场延拓，震源波场与接收波场做零延迟互相关，最后对成像结果运用Laplace滤波过滤低频假象。本方法相比于传统方法的优势在于继承了传统方法在
【机 构】
：
中国科学院地质与地球物理研究所地球深部重点实验室
论文部分内容阅读
本文针对地震波速度应力方程，构造了求解方法，并通过理论分析和大量的数值实验表明，该方法具有较好的时空稳定性以及较高的波动方程求解精度效率。本文将优化的辛格式交错网格褶积微分算子作为地震波逆时偏移算子，对叠前资料作振幅补偿处理，以地面速度分量记录作为边值条件进行波场延拓，震源波场与接收波场做零延迟互相关，最后对成像结果运用Laplace滤波过滤低频假象。本方法相比于传统方法的优势在于继承了传统方法在计算效率，内存要求、易于并行等方面的优点，同时在精度与长时间计算能力方面有明显的改进，在成像中表现为在不显著增加内存要求和计算量的前提下，运用较少的震源成像即能获得准确的构造剖面，这大大提高了成像效率。
其他文献
本文用水玻璃胶结石英砂，压制成型然后烧结。实验过程中进行了多角度的研究，摸索出了一套制备可控裂隙人造砂岩的方法，并且各向参数都达到了理想的效果。分别检测纵横波在含气和含水条件下的波速，然后分析得到含水含气条件下的纵波各向异性系数和横波分裂系数。
本文针对长偏移距地震资料，研究了适用于VTI介质的时差方程，并且比较了各时差方程的逼近精度，在对前人工作加以改进的基础上，得到了适用于水平多层VTI介质和各向同性层状的非均匀介质的时差方程。通过将该改进的时差方程与实际的走时方程进行误差分析，对比分析可以得出结论：该方法所得到的时距曲线能更好的逼近实际的走时曲线，并且该方法基本上不改变目前动校正公式的形式，只是对常规双曲时差方程中的速度平方项进行修
本文介绍qP波逆时偏移具有伪横波噪声压制、数值稳定性的特点。本文基于二阶耦合qP-qSV波方程，在逆时延拓过程中采用时间二阶、空间十阶的高阶有限差分方法，并结合PML完全匹配层吸收边界进行实现。滤波前后的脉冲响应显示了该方法的有效性。TTI介质叠前逆时偏移相对于传统的成像方法而言，可以获得更加精确的成像效果。本文给出的空间域压制伪横波噪声的方法简单易行、几乎不增加计算量，是一种高效稳定的偏移方法。
本文介绍无论均匀裂隙各向异性地壳还是强裂缝上地壳模型，SKS分裂T分量振幅和快、慢波时差均随方位变化，但各方位的整体振幅水平和时差均随裂隙密度的增大而增强。直立裂缝裂隙为HTI。对于HTI模型，在对称面两侧的对称方位，T分量上的振幅强度显示出方位对称性和极性反向，而ATI(任意空间取向TI)时这种记录特征消失。以往SKS分裂的资料处理和各向异性解释多采用HTI模型，T分量上振幅和极性的方位对称性可
本文利用zhang(2007)发展的三维横波分裂成像方法,采用1999年Izmit和Duzce地震长达6个月的4602个余震事件对土耳其北安纳托利亚断层进行各向异性成像。结果显示该地区具有强烈的各向异性,并且在断层附近出现明显的各向异性条带,覆盖断层附近3～4km,深达4km,并且结合之前的研究分析了各向异性的可能成因。
本文首次将ONADM方法应用于TTI介质中的声波方程模拟，研究TTI介质中的qP波传播规律。数值实验表明，在达到相同计算精度的前提下，ONADM方法的计算速度约为4阶Lax-Wendroff修正格式的16倍：存贮量仅需4阶Lax-Wendroff修正格式的28%。ONADM方法比Lax-Wendroff修正格式具有更快的计算速度和更小的存贮量需求，并能够更加有效地抑制数值频散。除此之外，本文还对二
本文对这两次气枪激发试验产生的地震信号特别是岸上地震台站记录到的气枪信号进行了分析，发现：①相同激发条件下岸上地震台站记录到的气枪信号的相关系数由信噪比决定，虽然信噪比较高时，气枪信号的相关系数随信噪比变化不大，但当信噪比低于一定程度时，气枪信号的相关系数受信噪比影响明显；②在不同激发压力下(8~14 MPa)，岸上地震台站记录到的气枪信号的振幅的变化规律与oBs水听器分量记录到的气枪信号振幅的变
本文针对地震波方程的二阶位移波方程，推导了卷积PML的时间域方程及二阶交错网格差分格式，给出了两种复合卷积迭代方法，并针对卷积PML的特点，给出了吸收域角点处不同于经典PML的处理方式。
本文利用IRIS上获取的LHZ频段数据计算出功率谱密度，画出各个台站功率谱密度图(spectrograms )再与相应的海冰年际变化趋势相比较。可以看出：当海冰面积减少时，格林兰岛DAG和ALE两处台站背景噪声信号强度减弱，但变化趋势却不是非常明显，而且另外三个台站的背景噪声强度似乎却和海冰年纪变化趋势无相关性。为了进一步分析结果，本文将各个台站的频率从小到大依次以05Hz为间隔划分为不同的频率段
为了检验球坐标系中的最短路径算法的计算精度和有效性，本文用Kennett等人(1995)得到的AK135地球速度模型计算全球49种不同的地壳、地幔、地核震相，并与AKI35全球主要震相理论走时表进行对比研究，结果表明：该算法能快速、有效、准确的追踪计算全球常见的主要震相，最大绝对走时误差约为0.1 s，而平均走时误差一般小于0.06 s。因此，该算法可作为区域或全球地震走时层析成像中多震相地震射线