复杂介质中地震多次反射波快速正演模拟

多次波的压制方法研究及多次波对一次波振幅随炮检距变化（AVO）影响的分析，需要实用的多次波正演模拟方法，提出适用复杂层状介质的多次波射线码的计算方法，使得多次波的射线追踪易于快速实现，同时，提出了减少多次波计算量的量板表法及插值法，模拟的多次波计算结果表明，多次波振幅的计算精度高，计算速度快，对文中设计的模型，振幅的最大相对误差为1．8％，而计算时间减少了0．6倍。     

碳酸盐岩喀斯特溶洞和裂缝系统的地震模拟与预测

提出了有助于确定碳酸盐岩裂缝带和喀斯特溶洞系统空间分布的地震正演模拟和预测方法。地震正演方法包含两个关键步骤，一是利用弹性波动方程计算含缝洞介质的有效弹性参数，另一步骤是根据有效弹性参数构成的等效地质模型，用单程声学波动方程进行波场延拓。数值模拟和实际地震数据的例子说明，将地震正演剖面与地震偏移剖面相结合进行对比分析有可能识别和预测碳酸盐岩地层中的缝洞发育带。     

用伪谱法合成零偏移距地震记录

零偏移距地震记录是叠后偏移的基础,常用射线追踪法来制作,但当介质复杂时,射线会出现交叉现象,对三维情况计算量更大。从二维标量波动方程出发,用伪谱法研究了零偏移距地震记录的合成,给出了不同情况下的数值计算公式,并对边界吸收问题作了有效处理。数值计算结果表明了方法的正确性和有效性。该方法可用于偏移中正演模拟的快速计算,且易推广到三维情形。     

叠前正演模拟与偏移的网络并行计算

在叠前深度偏移和非零炮检距声波方程正演计算过程中包含了大量的可并行计算的成分。作者在本文中提出叠前正演模拟与偏移的网络并行计算算法，并基于TCP／IP协议，将该算法设计成网络并行处理程序，极大地提高了计算效率。实际运算结果证明，本文提出的并行算法和技术路线是切实可行的。     

复杂模型高精度数值模拟研究

通过对声波方程采用二阶有限差分与四阶有限差分近似,进行数值模拟及波场分析,边界条件利用吸收边界和增加衰减带混合使用的方法。计算结果表明,四阶有限差分的精度明显高于二阶有限差分,并且边界吸收效果好。对于单炮记录难于模拟复杂模型的情况下,实现了单次叠加剖面和爆炸反射面的零炮检距剖面,并进行了偏移处理。地震剖面更加清晰地反映了复杂构造的信息,更好地实现了对复杂地质模型的正演模拟。     

转换横波共炮记录单程声波方程正演模拟

本文分别把分离的纵波和横波方程分裂为各自的两个单向旁轴方程（Ｊ．ＦＣｌｅａｒ－ｂｏｕｔ）,然后用单程纵波和横波方程,结合相位移加有限差分混合波场延拓法（Ｄ．Ｒｉｓｔｏｗ）,在共炮记录叠前深度偏移原理的基础上,借鉴几何地震学的射线法两点问题原理,提出合成转换Ｐ—ＳＶ波共炮记录的时间一致性原理和逆过程方法,从而推导出精确的转换Ｐ—ＳＶ波共炮记录正演模拟算法,由于单程延拓波场,可完全解决多次反射波问题。     

波动方程垂向插值斜缆采集正演模拟

为了解决海洋斜缆地震数据正演模拟时大多数检波器不在网格点上,不能直接记录波场值的问题,在使用有限差分法求解波动方程时,增加了垂向插值计算,得到了实际检波点的波场值。通过建立精细模型对正演模拟结果进行了误差分析,验证了该垂向插值方法的可靠性。分析对比了正演模拟数据和实际生产数据的波场特征,进一步证明了该波动方程与垂向插值正演模拟方法的正确性。该方法在倾斜电缆地震采集参数设计方面,具有较高的实用价值。     

煤田陷落柱特殊波对陷落柱解释的影响

对典型陷落柱建立圆锥状陷落柱数学模型,采用波动方程法对陷落柱模型进行地震射线追踪和波场模拟。模拟结果说明:由于陷落柱构造形态的特殊性,在陷落柱附近形成了正常煤层反射波、断陷点绕射波、延迟绕射波、延迟反射波及柱顶绕射波组成的纵波波场。在此基础上,模拟野外放炮,然后对生成的单炮记录进行常规处理,形成叠加、偏移剖面。通过对正演剖面和实际揭露资料分析,认为陷落柱特殊波的存在,是判别陷落柱构造的一个重要特征,但又造成了实际揭露陷落柱规模往往小于地震资料解释结果。将研究结果用于实际勘查,效果较好。     

AVO简化方程的物理意义及其在油气识别中的应用

地震反射振幅随炮检距变化研究(AVO)的基本原理是反射系数随入射角的变化而变化。这里从精确的Zoeppritz方程出发，根据Shuey及Hilterman提出的P—P波反射系数计算公式，设计多层地质模型进行AVO正演模拟并合成地震记录，利用地震记录对比讨论了Shuey及Hilterman的简化方程各项的物理意义及其在不同入射角范围内对地震记录的贡献。     

地震波场正演模拟中的有限单元方法

本文介绍了地震波场正演模拟中的有限单元方法。首先从粘弹性波动理论出发,推导了满足粘弹性介质的有限元方程,其次用不含阻尼矩阵的有限元方程模拟了地震记录。地表实例是自由界面的瑞雷面波。VSP记录有水平层、倾斜层;有零井源距和非零井源距;有线状和对称震源。 在程序中采用了超压缩存贮和分块计算措施,对降低有限元法要求计算机内存和运算速度的苛刻条件上做了有益的尝试。     

球坐标系下三维射线追踪的方法和应用

作者在文中详细地论述了适用于任意曲线坐标系的射线追踪系统的建立和重要公式的推导,并给出了该追踪系统在球坐标系下的具体应用,讨论了理论模型和实际三维模型下的射线追踪效果,球坐标三维射球追踪结果表明,该方法和弯曲法相比能更精确地反映介质的速度分布情况以及地震波在介质中的传播情况。     

频率域单程方程三维叠前正演中的波场替换

利用定位原理进行三维叠前正演,仅靠单程声波方程的延拓即可形成复杂地质条件下的正演记录,运算效率高无多次反射出现。但速度模型中对于过大或过小的速度常使波场在延拓时有频率和波数成分的损失,这种损失引起波场畸变并严重影响正演记录的质量。本文采用波场替换技术,在每次延拓时用零相位子波形成一个同时保留了原波场的动力学特征与空间位置信息的新波场,然后用新波场取代原畸变波场并参与正演记录合成。在频率域单程方程三维叠前正演中特别适合做波场替换,因为震源脉冲与接收响应(均可视为点脉冲)的延拓波场在形态上相似,而且在延拓波场中每道记录之振幅包络只有一个极值点。理论试算结果表明,采用波场替换技术后最终正演记录质量比替换前有了很大改善。      

基于P-SV波分离的VTI介质射线追踪方法

近年来,VTI介质参数建模和成像技术发展迅速,其中VTI介质射线追踪是基础,其精度和稳定性直接影响建模和成像的效果。VTI介质射线追踪大多通过前人提出的频散关系推导出程函方程和射线方程,该方法主要基于纵横波相互独立的假设,采用令横波速度等于零的手段,在某些特殊情况下精度较低、稳定性较差。针对其缺点提出一种新的射线追踪方法,该方法基于P-SV波分离得到的q P波频散关系,从而推导出程函方程和射线方程,提高了精度和稳定性。最后通过模型试算验证该方法的正确性和稳定性。     

复杂海底各种地震波的射线追踪与运动学特征

为了实现适应崎岖海底、大陡坡、海底火山等复杂海底地质条件且灵活、稳定、精度高的射线追踪方法,并基于该方法详细分析复杂海底条件下各种地震波型的运动学特征,综合多种算法的优势,实现了一种快速推进迎风插值射线追踪方法。首先,采用混合网格法剖分复杂海底地质模型;其次,通过融入迎风差分思想的线性插值策略来构建精度高、无条件稳定且灵活的局部走时和射线路径计算公式;然后,综合应用这些公式和多级次快速推进法,灵活计算整个模型各种类型地震波的走时,并基于逆向追踪方法计算射线路径;最后,对该方法进行了精度分析,发现其能够获得相对高的走时和射线路径计算精度,且反射波的计算精度远高于入射波。此外,计算实例分析显示,初至波中富含折射波和陡倾构造的反射波在很大偏移距还能被接收,崎岖海底各种波型传播路径复杂;基于此提出加大采集排列长度和采用直达波走时可提高复杂构造成像质量等复杂海底地震数据采集与处理方面的思考与建议。     

三维定位原理与三维反射波场模拟

在三维正演模拟定位原理基础上 ,以单程声波方程作为波场外推手段 ,实现三维叠前反射波场模拟 ,完全克服了全波方程模拟的 3个问题。在波场外推中 ,使用波动方程的混合解法 ,即以波动方程解析解 (相位移法 )为主 ,在介质速度横向变化的地方 ,用数值解 (45°有限差分 )作补充 ,使得整个算法精确、快速 ,既能处理复杂构造 ,又能适应速度的任意变化。     

基于ADIPI的单程波三维VSP正演模拟方法

三维VSP正演模拟在三维VSP观测系统设计、三维VSP波场分析及三维VSP偏移成像方法研究中均有重要的作用。以三维单程波动方程为基础,采用基于正交方向加插值(ADIPI)的傅里叶有限差分(FFD)算子,实现了三维VSP正演模拟。克服了对三维算子的双向正交分裂所导致的方位各向异性误差,既保证了模拟精度,又提高了计算效率,还可根据实际需要单独模拟上行波或下行波。模型试验结果证明了方法的有效性。      

基于模型参数化的地震波走时与射线路径计算

模型的建立与描述是地震资料分析与处理的基础,合理的模型参数化方式对于地震正反演的各个方面都有很好的效果,其中射线追踪在正演模拟、层析成像及偏移等研究领域中都占有很重要的地位。本文针对地震波走时与射线路径计算,设计了一种将离散模型连续化的最小二乘模型参数化方法,分别对速度模型和走时模型进行模型参数化,然后进行梯度速度模型验算及误差分析。计算结果表明,该模型参数化方法使走时计算精度从10^-2提高到10^-3,使射线路径的计算精度提高了14.33%。最后通过经典Marmousi模型和Sigsbee 2A模型实例验证,充分证明了该模型参数化方法的有效性及普遍适用性。     

基于保幅波动方程的广义高阶屏地震偏移成像方法

原始地震数据中饱含丰富的反映相位的走时信息和反映反射率的振幅信息。基于波场延拓理论的波动方程保幅地震偏移成像,是在给出正确构造成像位置的同时也给出真实反射振幅的有效完善。基于全声波方程,利用严格的解耦理论进行单程波动保幅分解,得到一个由波场传播项与振幅补偿项构成的,在走时与振幅上满足全声波方程对应的程函方程与输运方程的保幅单程波动方程;利用摄动理论进行单平方根算子渐进展开,推导出基于保幅波动方程的广义高阶屏地震偏移算子方程。模型测试和实际资料处理表明,该方法不但可以凭借更准确的相位归位和散射能量聚焦提高构造成像精度,而且输出了能更正确反映地下反射属性的能量信息,从而可以为更深层次的勘探开发,提供地球物理技术支撑。     

偏移速度分析与偏移成像

偏移速度分析和偏移成像是地震资料处理的两个重要组成部分.目前时间偏移技术比较成熟,而深度偏移技术也在逐步完善,时间域偏移成像主要推崇叠前时间偏移法.采用沿层叠加速度分析技术可获得层面上准确的叠加速度,并通过倾角校正、叠前时间偏移和CRP反偏移速度分析逐步优化速度,得到一个符合地质规律、准确的均方根速度场;通过深度偏移方法的研究,总结了建立精确偏移速度场的方法,并提出了一种地震资料处理的思路,即基于射线追踪的Kirchhoff偏移和基于波场延拓的波动方程偏移的结合,使偏移速度分析和偏移成像在应用效果和效率上得到了很大的提高.