弹性介质中垂直地震剖面有限差分正演

垂直地震剖面法可以通过认识反射和透射过程,深入了解地震波传播特性,又可以改善地表地震资料关于构造、地层和岩性的解释.本文是在基于垂直地震剖面的基本原理基础上,采用高阶交错网格有限差分数值求解二维弹性波方程,然后对常见的层状介质和岩溶介质模型分别进行模拟,根据模型合成了相关的垂直地震剖面记录,期间采用的是完全匹配层法吸收...     

各向同性介质中射线理论与有限差分地震波场模拟方法比较研究

地震波场数值模拟不仅是研究复杂地区地震资料采集、处理和解释的有效辅助手段,而且是研究地球深部精细构造和地球深部探测的有效工具.射线理论和波动方程理论是地震波场数值模拟的理论基础.射线理论主要刻画地震波在介质中传播的走时场、地震射线等运动学属性;波动方程理论通过求解波动方程来描述地震波在介质中传播的弹性动力学响应(能量衰减、相位特征、偏振属性、以及全波形等).基于波动方程理论的波场数值模拟由于能够引入丰富的波场信息,使得人们对不同介质中地震波的传播过程有了较全面的了解.本文以二维层状均匀介质模型为例,通过射线追踪法和交错网格有限差分法模拟得到的波场快照图、单炮地震记录剖面、合成理论地震图的分析比较,不但对地震波在各向同性层状均匀介质中的传播规律和特点有了深刻的理解和认识,同时又可以相互验证两种不同方法的正确性和有效性.     

地震波数值模拟中优化的通量校正传输方法

为了有效提高地震波正演数值模拟精度,提出了一种优化的通量校正传输方法,即在每一步时间递推过程中,沿坐标轴方向和对角线方向均引入通量校正参数进行波场校正处理,以高频散的均匀介质地震模拟波场为例,研究对比了不同中心网格有限差分近似阶数、传统的FCT方法和优化的FCT方法校正得到的地震波场效果.数值计算结果表明,传统的FCT方法容易产生虚假波场,而优化的FCT方法可以削弱或消除这些假波动,并指出结合高差分阶数和应用优化的FCT方法可以更好地提高数值模拟记录的信噪比.     

层状弹性介质全波射线追踪技术研究

射线追踪法是一种有效的地震波场数值模拟方法,在层析成像、叠前深度偏移及正演模拟等研究领域均占据重要地位.但相较于波动方程数值模拟方法,射线追踪法在获取较高运算速度的同时损失掉了地震波场的动力学信息,对于整个地震波场的地震波衰减以及转换波等问题都没有进行考虑.为了能够得到较为完整的合成地震记录并对地震资料做出更合理解释,本文在使用射线追踪法进行波场数值模拟的同时考虑了衰减信息,追踪转换波以弥补射线追踪法在动力学上的缺陷.在全波射线追踪过程中,通过追踪同类波和转换波,并计算该段射线的能量衰减,在保证高速计算的同时获得了较为完整的合成地震记录,取得较好的效果.     

地震波有限差分模拟综述

本文从有限差分法数值模拟技术的各个方面对地震波有限差分模拟的发展和现状进行了论述.波场的数值模拟技术是认识地震波传播规律,检验各种处理方法正确性的重要工具,地震波的数值模拟是地震波传播规律研究的必要手段,贯穿于地震资料的采集、处理、解释的整个过程中.有限差分法数值模拟技术相对于射线方法具有更高的精度,同时比有限元方法计算量小,因此在实际应用中占很重要的地位.     

VTI介质准P波旋转交错有限差分数值模拟

本文采用旋转交错网格差分格式对VTI（垂直对称轴的横向各向同性）介质准P波一阶应力-速度方程进行数值模拟。并在PML边界条件和稳定性条件下得出Marmousi等复杂模型的高精度波场快照和地震记录,分析了各向异性对地震波的影响。数值结果表明：旋转交错网格有限差分能获得高精度的地震模拟数据,PML边界有较好的吸收效果。     

黏弹TTI介质中旋转交错网格高阶有限差分数值模拟

以Carcione黏弹各向异性理论为基础,给出了适用于黏弹性具有任意倾斜对称轴横向各向同性介质(黏弹TTI介质)的二维三分量一阶速度-应力方程,采用旋转交错网格任意偶数阶精度有限差分格式求解该方程,并推导出了二维黏弹TTI介质完全匹配层(PML)吸收边界条件公式和相应的旋转交错网格任意偶数阶精度有限差分格式,实现了该类介质的地震波场数值模拟.数值模拟结果表明:该方法模拟精度高,边界吸收效果好,可以得到高精度的波场快照和合成记录;并且波场快照和合成记录能较好地反映地下介质的各向异性特征和黏弹性特征.     

VTI介质qP波方程高精度有限差分算子

波动方程有限差分法是一种使用广泛的地震波数值模拟方法.但是有限差分法本身固有存在着数值频散问题,会降低地震波场模拟的精度与分辨率.为了克服常规有限差分算子的数值频散,本文针对VTI介质地震波数值模拟问题,构造了频率-空间域qP波波动方程高精度有限差分优化算子,根据最优化理论中高斯-牛顿法确定了高精度有限差分算子的优化系数.利用常规差分算子和高精度优化差分算子对归一化相速度的频散关系精度进行了对比分析,并对均匀各向同性介质和均匀VTI介质中的qP波地震波场进行了有限差分数值模拟,通过频散关系精度分析和波场数值模拟结果表明：有限差分优化算子具有较高的波场数值模拟精度,有效压制了传统有限差分算子数值模拟中的数值频散现象,提高了有限差分算子精度,为VTI介质频率-空间域qP波正演模拟奠定了基础.     

基于接收阵列的时域地震波束形成方法

地震勘探中,针对目的层地震数据信噪比较低情况,为改善信噪比,本文提出一种基于接收阵列的时域地震波束形成方法(TSBBRA).该方法的理论依据是波束形成原理,文中通过对相控震源激发定向地震波场的分析,提出在反射地震勘探中,在信号接收端同样能生成定向地震波的思路.具体地,论文给出在地震数据集上先对单炮记录分组,再对组内记录时域内延时、叠加,生成对应定向地震波场的地震记录的方法.TSBBRA方法适于杂波干扰强,且杂波与目的层反射波方向不一致的地震条件.本文先通过波场模拟,展示了TSBBRA方法具有形成定向地震波场的能力;而后通过对比相控震源与TSBBBRA时间剖面,说明了TSBBRA方法的有效性;最后在吉林某油页岩矿区,用TSBBRA方法,分别针对浅层目标与深层目标,生成定向地震波数据.与原始数据相比,目标层信噪比提高了3~35 dB,在时间剖面上观测到的深部目标信号从200 ms扩展到1000 ms,TSBBRA剖面资料与测区附近区域的地质与地震资料一致.由此得出结论,TSBBRA方法有助于改善目的层地震资料信噪比,当关注深部目的层时,由于TSBBRA方法能有效压制噪声,因此可用于能够实现更大的勘探深度.实验中,我们给出了浅层勘探结果,由于本方法对震源及观测系统没有限制,实际应用中TSBBRA方法可作为深部资源勘探领域提高深部地震资料信噪比的有效方法.     

起伏层状介质中曲线网格有限差分法与射线法波场模拟对比研究

针对处理起伏地表(或含地下不规则波阻抗界面)条件下发展起来的地震波场数值模拟算法的模拟结果与解析解(大多数情形下无法得到)无法进行对比,且其有效性和正确性难以验证的情况,本文提出了一种可以相互验证波场数值模拟结果与射线追踪数值模拟结果的正确性和有效性的佐证方法,验证了参考射线追踪法.其中,波场数值模拟中采用曲线网格DRP/opt MacCormack有限差分法,射线追踪模拟则采用分区多步三角网格最短路径算法.通过系统对比上述两种方法得到的波场快照、单炮地震记录,以及合成理论地震图的结果显示,本方法相互作证了两种方法所得结果的正确性和有效性.双层和三层起伏层状模型的对比分析结果表明,这种方法不但可以加深理解地震波在复杂介质中的传播规律,同时射线法的引入为清晰识别和标定地震波场数值模拟中各种不同震相提供了一种便捷的途径.      

井间电磁场时域有限差分数值模拟

为进一步了解井间地下介质结构及其电性特征,本文用时域有限差分方法(FDTD),实现了井间电磁场响应的数值模拟.激发源是探测深度相对磁偶极子更大一些的电偶极子,首先选择一个位于井轴上的垂直电偶极子作为场源,并假设参与计算的介质相对于发射井井轴是轴向对称的,这样可将研究区域作为二维问题处理.推导了二维井间电磁波传播时域差分公式.该方法适用于任何方向入射的电偶极子源,尤其善于解决频域差分方法所难以描述的宽频脉冲.给出了井间金属圆柱和矿体圆柱数值模拟2个例子,结果表明,时域有限差分方法能有效的模拟井间地下介质中电磁波的传播,揭示电磁波传播规律.该方法速度快、精度高、结果稳定,适合用于井间电磁场反演成像的正演响应计算.     

基于非正交网格的带地形三维瞬变电磁场模拟

本文采用非正交网格时间域有限差分（FDTD）算法开展带地形瞬变电磁场的正演研究,并提出了一种采用非等权投影的改进非正交网格算法,提高了带地形瞬变电磁场三维正演算法的精度.对不同地形模型的数值模拟结果表明,地形的深度、宽度以及地形与源的相对位置等都可能对瞬变电磁场产生影响.     

高精度频率域弹性波方程有限差分方法及波场模拟

有限差分方法是波场数值模拟的一个重要方法,但常规的有限差分法本身存在着数值频散问题,会降低波场模拟的精度与分辨率,为了克服常规差分算子的数值频散,本文采用25点优化差分算子,再根据最优化理论求取的优化系数,建立了频率空间域中弹性波波动方程的差分格式;为了消除边界反射,引入最佳匹配层,构造了各向同性介质中弹性波方程在不同边界和角点处的边界条件. 最后由弹性波波动方程和边界条件,通过频率域有限差分法,分别利用不同震源对弹性波在均匀各向同性介质、层状介质及凹陷模型中的传播过程进行了数值正演模拟,得到了单频波波场、时间切片和共炮点道集,为下一步的研究工作（如成像、反演）提供了研究基础.     

基于卷积完全匹配层的非规则网格时域有限元探地雷达数值模拟

复频移完全匹配层(Complex Frequency-Shifted PML,CFS-PML)在长时间时域计算中对凋落波、倏失波具有好的吸收效果,并被广泛应用于时域有限差分模拟中.而本文采用卷积方法将CFS-PML应用于时域有限元求解GPR波动方程的数值模拟中.论文以TM波为例,推导了基于CPML(Convolutional PML)边界的时域有限元GPR波动方程求解公式,采用Newmark-β方法对时间导数进行离散,有效改善了时域有限元GPR数值计算程序的稳定性.并以狭长模型为例,开展了CPML边界中关键参数m、R和κ的选取实验,通过对比反射误差大小确定了综合最优参数组合.相同时刻UPML与CPML波场快照、3个检测点的反射误差比较,说明CPML较UPML具有更好的吸收效果.最后,采用非规则四边形网格对1个复杂GPR模型进行剖分,应用加载CPML边界条件的FETD程序对该模型进行了正演,得到了二维剖面法、宽角法正演GPR剖面图,说明非规则四边形对复杂模型的良好适应性,基于CPML边界条件的FETD可有效减少边界反射误差,能实现对任意复杂不规则模型的正演模拟.     

基于变换光学有限差分探地雷达数值模拟研究

在采用有限差分方法开展探地雷达复杂目标体精细结构模拟时,为了提高计算精度,常采用非均匀网格对目标区域划分小尺寸的网格,以压制离散网格频散现象和保证有限差分方法的稳定性.常规非均匀网格和自适应亚网格技术在网格剖分数量和粗细网格边界处理上难以达到计算效率和计算精度的均衡.本文根据隐形斗篷(invisible cloak)理论,将基于变换光学(Transformation optics)理论应用于有限差分探地雷达数值计算中.该理论的主要思想是基于目标参数变化而保持电磁场的传播不变性,在坐标变换后,Maxwell方程的形式可以维持不变,而使得相对介电常数与磁导率的表达式变得复杂.通过这种方式可以虚拟地扩大目标体所占的网格节点数,减少背景介质区域的网格数,不增加模型空间的网格总数.另外,这种网格划分方式不但提高了计算效率,同时也可以克服亚网格技术边界反射误差的影响.本文推导实现了基于变换光学的二维有限差分方法,通过典型探地雷达模型测试,对比分析了该方法与常规有限差分、变网格有限差分和自适应亚网格有限差分的优缺点.计算结果验证了基于变换光学的有限差分可用于探地雷达目标精细结构模拟,具有较高的计算精度和计算效率.     

分段光滑曲线边界波动方程数值模拟研究

矩形网格有限差分法在地震波传播数值模拟方面具有计算速度快的显著优势,但该方法在处理复杂边界问题上存在着效率低的严重缺陷.本文针对分段光滑曲线边界定义了尖点处的一种正则导数,给出了矩形网格情形分段光滑曲线网格边界点法向导数的一种插值计算方法.采用矩形网格有限差分法对复杂边界地球介质模型进行地震波场数值模拟,并采用波场系列快照技术揭示地震波在起伏地表和复杂介质中的传播规律.模拟结果表明：法向导数插值计算方法为矩形网格有限差分法处理复杂边界提供了有效途径,采用波场系列快照技术可以清晰地展现地震波在反射界面的反射和透射规律、在尖点的绕射规律以及在自由表面的直达波和多次反射规律.     

黏弹性VTI介质频率空间域准P波正演模拟

有限差分方法是波场数值模拟的一个重要方法,时间域有限差分计算方法因按时间片递推计算,每个时间片的舍入误差会累积到下一片中,当时间片较多,最终会导致累积误差太大.而频率域计算是按频率片对空间网格进行整体求解方程组,其计算误差分配到了每个网格点上,并且各个频率片之间是独立计算的,因此不存在累计误差,而且在频率-空间域更易于...     

适于声波方程数值模拟的时间-空间域有限差分算子改进线性算法

有限差分法广泛应用于地震波场的数值延拓,确定合适的有限差分算子以减小数值频散是有限差分法的一个重要研究内容。近年来为了进一步抑制数值频散和增加时间步长,新的有限差分模板得到了应用,对于此,前人使用泰勒展开方法和最小二乘方法确定有限差分算子系数。本文在以前工作的基础上,使用改进的线性方法确定新模板的有限差分系数,并与传统模板线性方法进行对比;通过频散分析和正演模拟验证出新模板线性方法能够更好地保持频散关系,在相同的精度下效率提高了一倍,从而说明了改进的线性方法的有效性。