简介：坐标变换法通过将物理空间的曲网格映射为计算空间的矩形网格,将起伏地表转化为水平地表,同时将物理空间的波动方程转化为计算空间的波动方程,在计算空间完成数值模拟,坐标变换的方法对处理起伏自由边界具有较好的适应性和应用效果。本文在传统坐标变换方法的基础上,根据计算区域速度差异采用不同的网格大小和采样时间步长,提出了一种基于时空双变网格的起伏地表坐标变换正演模拟方法。在编程实现算法的基础上,通过典型模型波场模拟试算结果分析可知：（1）变网格方法与常规方法波场模拟误差在0.5%左右;（2）变网格方法计算效率视不同的变网格区域面积及变网格大小可提高几倍量级,在本文模型和计算参数下提高约5倍。（3）在满足模拟精度及频散条件要求下,变网格方法较全局细网格算法能显著节约计算内存。为此,针对起伏地表数值模拟,本文方法具有较高的模拟计算精度和一定的适应性。

简介：频率-波数域单程波算子能高效地模拟地震波在复杂介质中的传播，但是在描述波的大角度传播和速度横向扰动变化较大介质中传播的问题时仍然存在一定误差。这类误差是由于对单平方根算子使用Taylor展开式的近似程度不足所造成。为了进一步提高泰勒展开式的精确性，本文提出一种利用粒子群智能算法优化级数展开系数的高阶广义屏算子对单平方根算子的展开级数进行优化处理。新的偏移算法能在保持单程波偏移算法高效的前提下进一步提高偏移算子在大角度的成像精度和对强横向速度变化介质的适应性。通过脉冲响应实验，验证了基于粒子群算法优化级数的高阶广义屏算子能够提高常规的高阶广义屏算子的成像精度和成像角度。根据对二维SEG/EAGE盐丘模型的成像处理，基于粒子群算法优化级数的高阶广义屏算子对盐丘下面的断层取得了更高质量的成像，说明粒子群优化级数的高阶广义屏算子比常规的高阶广义屏算子具有更好的横向速度适应性。为了检验本文所提算法对实际资料的处理能力，我们利用常规的偏移处理技术和本文所提算法对一条海上二维数据进行了偏移成像处理，对比分析成像剖面发现本文所提算法描述了更加清晰的层位信息和更高质量的偏移剖面。本文所提算法能有效提高高阶广义屏偏移在广角度成像的能力，具有一定实际应用价值。

简介：针对转换波各向异性叠前时间偏移实践中的四参数速度模型估计中参数间相互影响、难于准确的确定,本文将简化的转换波两参数动校正方程从叠加速度分析拓展到叠前时间偏移速度模型修正,形成了基于两参数动校正方程的新的四参数法转换波叠前时间偏移速度修正新方法和流程。在该方法和流程中,先进行转换波两参数叠加速度分析,获得转换波叠前时间偏移初始的速度和各向异性参数,然后通过分析偏移后的共成像点道集中的剩余时差修正速度和各向异性参数。其中叠前时间偏移速度模型的垂直速度比仍是需要利用纵波和转换波叠加剖面,层位标定后,利用相关的方法获得。有效速度比需结合纵波速度分析,利用Thomsen（1999）公式计算得到初始模型,最终依靠百分比扫描偏移处理获得最终的有效速度比模型。该方法简化了转换波高质量成像的速度估计方法,减小了多个参数估计的不确定性,在实际应用中也取得了较好的成像效果。

简介：多分量勘探相比常规纵波勘探而言有许多的优势，PP／PS联合AVO分析和反演是强而有效的储层识别方法。在本文中，我们推导了更精确的转换波AVO反射系数公式；并在入射角小于30度时进行了进一步简化，其对于转换波AVO分析和反演而言更加简单有效。基于该近似，我们进行了PP／PS联合AVO反演。实际资料的实例表明，反演得到的纵横波速度比可很好地识别岩性及油气。反演得到的流体因子和泊松比等其他属性在储层处也显示出明显异常，效果显著。

简介：随着并行计算技术的发展,非线性反演计算效率在不断提高,但对于基于单点搜索的非线性反演方法,其并行算法的实现则是一个难题。本文将群体搜索的思想引入到基于单点搜索的非线性反演方法,构建了并行算法,以量子蒙特卡罗方法为例进行了二维地震波速度反演及实际资料波阻抗反演,并测试了使用不同节点数进行计算的效率。计算结果表明：该并行算法在理论和实际资料反演中是可行的和有效的,具有很好的通用性;算法计算效率随着使用节点数的增加而提高,但算法计算效率的提高幅度随着使用节点数的增加逐渐减小。

简介：为了进一步提高大地电磁三维快速松弛反演的计算效率，在深入分析大地电磁三维快速松弛反演算法的基础上，结合MPI自身的优越性，确定了并行计算的思路，实现了三维快速松弛反演的并行计算。通过理论模型合成数据和实测数据对实现的三维快速松弛反演并行程序进行了试算，分析对比了在多种情况下程序的执行效率。测试结果表明，所实现的三维快速松弛反演并行程序运行结果正确，效率提高明显。此成果在普通微机上实现，推动了大地电磁三维反演技术的实用化，可为其它地球物理三维正反演研究所借鉴。

简介：多波勘探是提高致密、低渗透等复杂油气藏勘探开发精度的有效手段,但转换波属于典型的低信噪比、低分辨率资料,如何最大程度提高转换波纵向分辨率已成为多波处理的一个难点。谱反演技术能有效解决常规反褶积提频技术受频带限制影响,提高分辨率能力有限的问题,能最大程度提高数据分辨率,便于识别薄层,但该技术的难点是如何通过稳定的反演算法得到高精度反射系数,如何把反射系数恢复为宽频的地震数据？本文在前人研究的基础上利用全变差作为先验信息有助于求解欠定问题的优势,提出了一种基于全变差约束的最小二乘反演算法,提高反演的精度和稳定性,并利用高斯拟合振幅谱模拟得到宽频子波数据,通过在蓬莱地区的高分辨率数据恢复处理,得到分辨率更高的转换波。理论试验和实际资料应用证明了该方法能较大幅度地提高转换波资料分辨率,为后续速度反演、储层反射信息提取提供更准确的数据。

简介：针对矩形网射线追踪存在的模型剖分灵活性差、速度界面描述精度差等问题,研究了复杂结构三角网最小走时射线追踪全局算法。（1）根据剖分区域点、线、面的结构关系,遵循Delaunay三角剖分的优化准则进行三角网格剖分;（2）定义三角单元射线追踪的拓扑关系;（3）波源点及某一时刻波到达的每一个节点点构成波行面,在波行面扩展过程中计算节点的最小走时和次级源位置,实际次级源检索采用双曲线近似算法;（4）利用各节点走时和次级源方向信息,通过最小走时搜索,拾取从接收点到源点的射线路径。数值模拟结果表明,三角网射线追踪方法模型剖分时灵活性强、速度间断面的描述精度高,追踪结果准确。

简介：正则化因子的选择方法决定了正则化方法的稳定性和准确度。在分析改进的奇异值分解法与正则化方法的关系的基础上，给出一个正则因子计算公式。数值模拟试验表明改进的奇异值分解法和正则化方法适应低信噪比。当信噪比低于30时正则化方法比奇异值分解法结果更好，当信噪比高于30时奇异值分解法更优。采用本文提出的正则因子的正则化方法可以适用于实际核磁测井的低信噪比（5〈SNR）情况。数值模拟及实际资料的计算实验表明，改进的SVD反演算法与正则化方法都具有适应低信噪比、求解速度快等优点，可以满足核磁共振测井工作的需求。

简介：纵波和转换波联合的多波地震勘探技术是解决复杂油气勘探的有效技术，提高转换波的分辨率是其关键问题之一。影响转换波分辨率的主要原因是地层对地震波的吸收，有效地计算地层Q值、消除地层吸收对转换波传播的影响，是提高转换波分辨率的关键。本文提出了从叠前转换波道集中估算横波Q值的方法，并利用沿射线路径的波场延拓，将一种稳定有效的反Q滤波方法应用到叠前共炮点纵波和转换波道集的衰减补偿中。模型资料结果表明，本文提出的估算转换横波Q值的方法精度较高；模型资料和实际资料的吸收补偿结果表明，此稳定全反Q滤波能有效地提高叠前纵波和转换波资料的分辨率。

简介：转换波技术是解决气云区成像问题的有效手段，而二极化现象则是制约转换波构造成像精度提高的一个重要因素。解决该问题的一个有效办法是二极化校正：将转换波速度分解为基本速度和速度扰动，分两步对转换波地震数据进行叠加速度分析和叠前时间偏移速度分析。柴达木盆地三湖地区天然气资源丰富，气云广泛发育，通过对三湖地区二维三分量地震数据的研究分析，发现该地区转换波二极化现象明显，对转换波的准确成像有较大影响。通过对转换波二极化校正方法的应用，不论是在叠加成像，还是叠前偏移成像都取得了较好的效果，提高了气云区转换波构造成像的精度。二极化校正对解决气云区转换波剩余静校正有重要促进作用，二极化校正与转换波剩余静校正的结合是陆上气云区转换波高精度成像流程中必不可少的环节。

简介：基于曾新吾和MacBeth提出的时域内双源累积旋转方法,建立了频率相关介质中分析多分量VSP数据的横波分裂参数提取算法（DCTF）。该算法可以在频域中针对单个频率提取横波各向异性参数（快横波的极化方向以及快、慢横波间的时间差）,从而避免了目前常用方法中使用窄带通滤波可能带来的误差。通过对地震合成记录的数值分析,确定了该算法的可行性和正确性,并与目前常用方法的应用结果进行了比较。结果显示,频率相关横波分裂参数可以利用DCTF从地震四分量数据中直接提取。在地震频率范围内,含较大尺度裂缝时各向异性参数将表现出频率相关性,这意味着在地震频率范围会出现频散。随着频率的增加,各向异性有降低的趋势。

简介：在具有垂直对称轴横向各向同性介质中,利用四种参数来确定中间至远偏移距转换波(C-波)动校正.它们是C-波叠加速度VC2,垂直速度比和有效速度比γ0和γeff,以及各向异性参数χeff.我们将这四种参数作为C波叠加速度模型.C-波速度分析的目的就是确定这种叠加速度模型.C-波叠加速度模型VC2,γ0,γeff,和χeff可以由P-波和C-波反射动校正资料获得.然而错误的传播是C-波反射动校正反演中的严重问题.当前短排列叠加速度由于是从双曲线动校正推算而得,因而其精度不足以为各向异性参数提供有意义的反演值.中间偏移非双曲线动校正不再被人们所勿略,而是可以用一个背景γ加以量化.非双曲线分析通过中间偏移距的γ校正量可以产生VC2,若数据不含燥音,其误差小于1%.方法稳健,允许γ启始假定值的误差达20%.该方法也适用垂直非均匀各向异性介质.精度的提高使能够用4分量地震资料计算各向异性参数.为此提出了两种工作流程:双扫描和单扫描流程.理论数据和实际数据的应用表明这两种流程得出的结果其精度相似,但是单扫描流程比双扫描更有效.

简介：随着全张量重力梯度（FTG）测量技术的不断发展,重力梯度数据的三维反演技术在油气和矿产勘探中日益受到重视与关注。为了快速处理和解释大规模的高精度数据,图形处理器GPU（GraphicsProcessingUnit）和预处理分解技术（Preconditioningmethods）在地球物理反演中的使用变得十分重要。本文结合对称逐次超松弛（SSOR）技术与不完全乔列斯基分解共轭梯度算法（ICCG）提出改进的预处理共轭梯度法,并考虑到方法预处理分解占用额外的时间,开发该算法的GPU并行算法来提高加速效果。然后通过含噪的模型数据反演来证明改进的并行预处理方法在三维全张量重力梯度数据反演中的适应性。由此,基于NVIDIATeslaC2050GPU的并行SSOR-ICCG算法和在2.0GHzCPU上的串行程序比较,达到了大约25倍的加速比。最后,我们将该算法应用于美国路易斯安那州南方Vinton盐丘的实测航空重力梯度数据反演中,反演出良好的反演结果,验证了该方法在三维全张量重力梯度数据快速反演中的优势和可行性。