三维混合延拓一步法波动方程正演模拟

提出三维相位移加有限差分混合延拓法,以此为基础,实现速度纵横向变化,构造任意复杂的三维介质一步法正演模拟,为三维复杂波场的地震,地质解释提供了精确,实用的方法。对于二维介质,相位移延拓方法对陡倾斜地层成像精度高,稳定性好,计算快,但实现的速度的横向变化困难,而频率、空间域的有限差分法算法简单,稳定性好,能适应速度的纵横向任意变化,但偏移陡倾地层存在很多问题,用快速的４５°有限差分对相位移延拓作补充     

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

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

基于傅氏变换有限差分的三维波动方程正演模拟并行算法

傅氏变换有限差分(即FFD)混合波场延拓法波动方程正演模拟精度高,但三维FFD因其计算复杂、数据存储巨大而难以推广。这里提出三维FFD一步法正演模拟并行算法,把一个串行的大规模频率循环分割成并行的多个小规模频率循环,并发执行多个波场延拓来提高运算速度,并通过由给定的资源设定并行数,以及适时重排三维数据,大大增强了算法的适应性和时效性。通过对比并行和非并行正演实验,验证了算法的有效性。     

波动方程显式短算子三维正演模拟

通过对三维波动方程的显式差分法的研究,并对模小于1.0的稳定性约束作线性化处理时引入加权因子,提出了一个基于二次规划方法的一维显式短算子设计方法.与以往的显式算子设计方法的不同之处,在于对传播区内算子的模也引入了最大值限制,这样可更严格地保证深层波场延拓的稳定性.同时,给出了基于一维显式算子及McClellan变换的三维波场深度延拓方法,大大减少了计算量,提高了计算速度,并能更好地处理介质横向非均匀性.     

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

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

三维高密度电法高分辨数值模拟与分析

随着三维高密度电阻率法在工程中的推广与应用, 提高资料的处理与解译的分辨率的研究显得十分重要。本文基于三维高密度电法的理论, 提出建立不同地电模型、利用有限差分法和最小二乘法进行资料处理。论文进一步对模拟结果进行空间定位分析和形态效应分析, 进而得出异常体的变化规律。模拟与分析表明, 三维高密度电法在探测异常体的空间定位上具备较好的准确度。     

某铬铁矿区地震波场正演模拟研究

依据某铬铁矿区地质、物性资料,分别对矿区铬铁矿体、含矿超基性岩体、地质剖面进行地震正演模拟研究,分析对比模拟剖面和实际地震剖面的波场特征,发现实际剖面上存在与模拟剖面含矿位置相似的波场特征,结合地质资料,推测含矿岩体的反射区域,辅助金属矿区地震勘探对矿体的识别。     

非均匀弹性介质地震波场伪谱法正演模拟

基于一阶速度—应力弹性波波动方程,利用傅里叶变换将空间导数变换到频率域求解,对介质波场进行正演模拟。本文结合实际勘探中可能存在的地质问题,设计了多种地质模型,包括背斜模型、向斜模型、透镜体模型和断层--侵入体模型等,进行了伪谱法地震波场正演模拟、合成地震记录的实现,并分析了其波场动力学变化特征。结果表明,应用伪谱法对复杂模型的模拟有较高的精度和可靠性,可直接用于辅助地震野外工作设计、反演解释和资料解释等工作。     

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

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

地震波混合阶褶积算法模拟

基于正反傅立叶变换,提出了地震波模拟的混合阶褶积算法。该方法原理简单、易于实现,结合了有限差分与伪谱法的优点,具有较高的精度和计算效率,适用于地震波场的正演计算,同时还给出了二维地震波场的理论计算实例。计算结果表明,此算法模拟结果正确、精度高、速度快、能适应较为复杂地质模型,并且易于推广到各向异性介质中去。     

复杂地质构造波动方程数值模拟

利用全声波方程进行地震波场的数值模拟，虽然可以计算全声波波场值，但是也常常会出现层间多次波，给资料解释造成一定的困难。作者在本文中根据无反射声波方程，采用有限差分数值解法，实现了自激自收时间剖面的制作，避免了层间多次波的产生。本方法不但允许速度有纵向变化，而且允许速度有横向变化，这使得模拟复杂地质构造的地震响应成为可能。实际计算结果表明，该方法是一种行之有效的自激自收时间剖面的计算机数值模拟方法。     

二阶弹性波动方程高精度交错网格波场分离数值模拟

给出了一种等价的二阶弹性波动方程,以解决弹性波场中完全弹性波动方程不能完全分离耦合的纵、横波波场问题.应用高阶交错网格有限差分法求解该波动方程,并使用通量校正技术(FCT)进一步压制频散,采用均匀介质模型和层状介质模型进行波场分离数值试验,精确得到了混合波场、完全分离的纯纵波及纯横波波场.数值结果分析表明,本文方法在均匀介质情况下准确可靠,在分离后的纯纵、横波波场中可观察到较为丰富的能量转换信息,这对认识复杂弹性波的传播规律及弹性波理论具有重要意义.     

三维任意形体高密度电阻率法数值模拟方法

为研究高密度电阻率法对各种三维地电体的反映规律，本文介绍了正演程序编制的方法并给出了计算实例。     

槽波地震勘探的数学模型研究

为研究槽波的形成及地震时间剖面上的表现规律，利用非均匀各向同性介质的波动方程，应用有限差分法进行求解并编制了计算程度，给出了一个3层不对称模型的两分量计算结果实例。     

基于散度和旋度的弹性波波场分离数值模拟方法

弹性波波场分离数值模拟方法是研究混合波场中纯纵波和纯横波波场传播规律的一种重要手段,通常采用算子分裂来实现,这不但增加了计算量,而且还浪费了大量计算时间.因此提出了基于散度和旋度的弹性波波场分离数值模拟新方法,仅依赖于交错网格来计算出散度场(纯纵波)和旋度场(纯横波),从而提高计算效率.在人工截断边界处采用完全匹配层吸收边界条件(PML).数值模拟结果表明,该方法能够从混合波场中准确地分离出纯纵横波波场,且精度高,边界吸收效果好,计算速度快,计算过程稳定,可用于指导实际多波多分量地震资料处理,同时指出本文算法对相互耦合的波型(如面波)无法实现准确分离.     

直流电阻率三维正演的预条件共轭梯度法研究

三维电阻率法对反演的精度和速度的要求越来越高,而正演是反演的基础,因此直流电阻率三维正演计算的速度和精度是三维电阻率反演实用化的关键。这里利用对称超松弛预条件共轭梯度法(SSOR-PCG),求解有限差分法离散生成的大型稀疏线性方程组,预条件矩阵的选择大大降低了系数矩阵的条件数,结合矩阵的一维非零元素压缩存储模式,使得正演计算速度得以提高,而内存占用量明显减小。在直流电阻率三维正演中采用异常场法,提高了电源点附近的解的精度。利用编制的有限差分正演程序,对两层模型、垂直接触带模型和低阻异常体模型进行了数值模拟,计算结果表明该算法是可行的,且可以明显提高正演计算的速度和精度。