简介：中国西部复杂山地是地震勘探的难点区域,地表和地下都复杂,从而导致地表激发困难,地震资料干扰严重、信噪比低。地下构造形态复杂,导致反射点分散严重,前人尝试过弯线采集,但很少考虑地下实际构造特征。本文章提出基于CRP的不规则采集设计方法,使地表避开难以激发的区域,地下考虑构造特征及CRP点分布位置,优化地表接收线位置,通过地表有目的布置接收点位置,最大程度的保证地下CRP点不分散,从而提高地震数据的信噪比.并以四川盆地某工地的实际资料通过模型数据验证了方法的合理性与可行性。表明该方法能够解决复杂地区地震采集存在的部分难题,为复杂地区的地震勘探采集提供一个新的思路。

简介：本文在Oristaglio等(1984)和Adhidjaja等(1985)工作基础上,给出线源二维时间域瞬变电磁二次场的DuFort-Frankel有限差分数值解,有效避免了在总场求解法中场源附近的奇异问题,并对地-空边界电导率的处理、归一化感应电动势偏导数的计算、推进时间步的确定,提出了改进方法;吸取前人成就中二次场地-空边界向上延拓和零值边界处理技术,从而简化了计算方法;通过对均匀大地、水平层状大地模型的计算,二次场求解法与解析法的最大相对误差小于0.01%,计算速度比总场求解法提高了约3倍;模拟计算不同时刻瞬变电磁场在地下的分布形态,描绘出感应涡流向下向外的传播特征,以及与地下异常体相互作用的物理过程。

简介：中国西部复杂山地是地震勘探的难点区域,地表和地下都复杂,从而导致地表激发困难,地震资料干扰严重、信噪比低。地下构造形态复杂,导致反射点分散严重,前人尝试过弯线采集,但很少考虑地下实际构造特征。本文章提出基于CRP的不规则采集设计方法,使地表避开难以激发的区域,地下考虑构造特征及CRP点分布位置,优化地表接收线位置,通过地表有目的布置接收点位置,最大程度的保证地下CRP点不分散,从而提高地震数据的信噪比.并以四川盆地某工地的实际资料通过模型数据验证了方法的合理性与可行性。表明该方法能够解决复杂地区地震采集存在的部分难题,为复杂地区的地震勘探采集提供一个新的思路。

简介：为了有效处理人工源的影响,本文开发了带源的CSAMT二维正反演算法,可用于全区（近区、过渡区和远区）资料的反演。引入正则化因子完成磁法二维反演,并且将模型参数调整为磁化率的对数,保证反演过程中磁化率始终为正值。本文基于交叉梯度原理,将CSAMT和磁法进行联合反演,通过搜索交叉梯度项权重的方法,避免了不同异常源引起的两种异常相互干扰的问题。理论模型算例表明基于交叉梯度的联合反演方法优于单独反演。本文开发的带源CSAMT二维正反演算法,有效处理了人工源的影响,保证了最终联合反演算法的可靠性。

简介：频率域全波形反演充分利用全波场的振幅、相位以及频率信息，采用较少的频率便能反演得到精度很高的速度模型。本文以有限单元法为基础，对起伏地形条件下二维声波频率域全波形反演进行了研究。在正演算法中，针对截断边界问题，并考虑多频率联合反演中计算区域采用同一套剖分网格的需求，提出了一种适用于起伏地形的衰减边界条件算法。该算法的核心思想是在控制方程波数项中引入衰减因子，通过一定方式调节衰减因子使得声波在衰减层中充分衰减，达到压制截断边界影响的目的。根据指数衰减规律，文中推导出了一种新的衰减因子计算公式，并给出了不同频率条件下衰减层厚度计算公式；在反演算法中，采用共轭梯度法求解高斯牛顿反演迭代方程组，避免直接求解雅克比矩阵和Hessian矩阵带来的巨额计算量，并采用相同的反演模型，对比分析了不同初始模型和频率组合对全波形反演结果的影响。起伏地形模型数值模拟和全波形反演数值试验表明，本文提出的指数衰减边界条件算法和基于该算法的全波形反演算法具有很好的应用效果。

简介：地震数据中的面波是严重降低地震资料信噪比的干扰波,它的存在影响了后续地震资料的处理与解释。本文根据地震记录中面波与反射波信号形态结构的差异,采用基于二维字典形态成分分析方法对面波噪声与反射波进行分离。根据面波信号的低频、低视速度和频散的特性,选择二维非抽样离散小波变换作为面波的稀疏表示字典,根据反射波局部相关性较强的特点,选择二维局部离散余弦变换作为反射波的稀疏表示字典,构建地震记录在联合二维字典下的稀疏表示模型并采用块协调松弛算法进行求解,将地震记录分解为反射波部分和面波部分。对合成地震信号以及实际地震资料的处理结果表明本文方法不仅能有效压制强能量的面波干扰,而且还能很好保护反射波信号的波形。

简介：无网格法形函数构造不依赖预定义的单元，具有计算精度高、处理复杂模型便利等优点。本文介绍了无单元Galerkin法（EFGM）、点插值法（PIM）与径向基点插值法（RPIM）三种全域弱式无网格法的近似原理及特点；以二维泊松方程为例研究了支持域无量纲尺寸、场节点与背景网格设置对无网格法计算精度的影响。将RPIM与EFGM应用于频率域线源二维正演，给出了RPIM形状参数的推荐值；分析了均匀介质模型大地电磁（MT）二维正演无网格法边界条件直接加载与罚函数法加载的精度差异，结合PIM与RPIM边界条件加载便利及EFGM计算复杂模型精度高的优势，提出了EFG—PIM及EFG．RPIM耦合算法，数值计算结果验证了耦合算法的有效性。研究发现：无网格法及其耦合方法适用于电磁法数值模拟；支持域无量纲尺寸取1．0时无网格法精度与效率高，场节点与背景网格重合时计算效果佳；泊松方程求解PIM及RPIM精度较EFGM低，计算均匀介质MT响应精度较EFGM高；RPIM改善了PIM计算涉及的奇异性问题，对应支持域无量纲尺寸选择空间大。

简介：频率空间域地震波数值模拟具有独特的优势:可以同时模拟多源的波传播、每个频率之间独立并行地计算、计算频带选择灵活、不存在累计误差、容易模拟粘弹性介质中地震波传播。但是该方法的最大瓶颈是对于计算机内存的巨大需求。我们使用压缩存储系数矩阵的方法,极大地减少了计算机内存的需求量。同时为了减少短差分算子的数值频散,引用了频率空间域25点弹性波波动方程的差分格式,并使用了最小二乘意义下求出的优化差分系数。为了克服边界反射,采用了最佳匹配层吸收边界条件。数值模拟试验证明:用压缩存储系数矩阵及优化差分系数的频率空间域25点差分格式进行弹性波正演模拟,可以减少数值频散,提高计算精度。使用较大的网格间距,降低计算机内存需求,并保持较高的计算效率。该正演方法为后续弹性波偏移和弹性参数反演提供较好的基础。

简介：在地震勘探数据采集中,随机噪声严重影响地震资料质量,给后期解释工作带来很大困难。如何在不损失剖面有效信息的前提下压制随机噪声,有效地提高地震资料的信噪比和保真度,是本文的研究目标。构造导向滤波技术的核心是构造方向表征的求取以及如何实现非平稳滤波,来达到提高地震数据信噪比和保真度的目的。本文首先通过分析函数二维导数与希尔伯特变换的频率响应关系,推导出了基于二维希尔伯特变换的非迭代地震同相轴倾角求取算子,进而达到了构造方向表征的求取;其次选取多项式拟合作为构造导向滤波中的非平稳滤波方法,扩展了非平稳多项式拟合的应用范围;最后沿构造倾角方向进行变振幅同相轴的非平稳多项式拟合,实现和构建了新的自适应构造导向滤波方法。理论模型和实际地震资料处理的结果表明,所提出的方法实现了既保护构造信息又有效地压制了随机噪声的目的。

简介：在地震勘探数据采集中,随机噪声严重影响地震资料质量,给后期解释工作带来很大困难。如何在不损失剖面有效信息的前提下压制随机噪声,有效地提高地震资料的信噪比和保真度,是本文的研究目标。构造导向滤波技术的核心是构造方向表征的求取以及如何实现非平稳滤波,来达到提高地震数据信噪比和保真度的目的。本文首先通过分析函数二维导数与希尔伯特变换的频率响应关系,推导出了基于二维希尔伯特变换的非迭代地震同相轴倾角求取算子,进而达到了构造方向表征的求取;其次选取多项式拟合作为构造导向滤波中的非平稳滤波方法,扩展了非平稳多项式拟合的应用范围;最后沿构造倾角方向进行变振幅同相轴的非平稳多项式拟合,实现和构建了新的自适应构造导向滤波方法。理论模型和实际地震资料处理的结果表明,所提出的方法实现了既保护构造信息又有效地压制了随机噪声的目的。

简介：工区位于在大民屯凹陷荣胜堡-前进地区工业繁华地带,地表水系发育,表层岩性复杂,深层各种油气藏圈闭类型众多,勘探难度大,以往采集资料质量差.本文详细介绍了为改进资料质量所采取的一系列措施,其中包括科学的观测系统设计、认真细致的表层调查、高精度的测量,选择最佳岩性组合激发、低噪声接收、大型障碍区特观施工、定量化质量控制等工作.成功地进行了二次三维采集,实现了二次采集资料质量的飞跃,减少浅层资料缺失,保证了深层资料品质,保证了资料解释的准确性.文中展示了了实际效果.

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

简介：本文提出基于原始含表层多次波数据实现叠前共炮集地震数据插值。相对于利用相邻道的信息变换或外推插值用于缺失的地震数据重建，本文方法利用表层多次波数据互相关构建准一次波，将蕴含在表层多次波数据中的，而在采集记录中表现为缺失的近炮检距信息提取出来，并在滑动时间空间窗内采用最小二乘匹配滤波和均方根振幅校正方法进行准一次波校正而后用于数据插值重建。本文方法适用于表层多次波比较发育，同时又存在数据缺失尤其是近炮检距数据缺失情况。方法易于实现，不需多次波和一次波的提取，利用多次波中蕴含的信息实现缺失的地震数据弥补，为含有表层多次波的数据进行近炮检距地震信息的插值重建提供了一个很好的思路。

简介：将偏移后的炮域偏移距道集转换为角度域共成像道集（ADCIGs）可为偏移速度分析（MVA）和叠前反演提供输入道集，并且ADCIGs是理论上没有假象的叠前反演道集，也是目前公认的精度最高的叠前反演道集。本文研究了基于矢量波场逆时偏移的弹性波保幅ADCIGs的提取方法，以保幅弹性波逆时偏移方程为基础，其核心是求取不同震源位置的纵、横波场在地下各成像点的入射角，对于转换波勘探，二者共享一个入射角，即震源纵波入射角。根据几何关系，震源纵波波场的传播角、构造的局部地层倾角之差为震源纵波入射角，震源纵波波场的传播角利用解耦后纵波场的极化向量得到，构造的局部地层倾角利用偏移叠加剖面的复波数得到。对纵、横波的共炮点偏移道集按入射角重新排列即可得到各自的ADCIGs。文中利用水平层状介质模型、倾斜层状介质模型、Marmousi—II弹性波部分模型和实测资料验证了算法的有效性，计算结果表明，本文方法计算的纵、横波角度具有较高的精度，提取的角道集具有较好的振幅保真性，能够为MVA和叠前反演提供可靠的输入道集。

简介：Theattenuationfactororqualityfactor(Q-factororQ)hasbeenusedtomeasuretheenergyattenuationofseismicwavespropagatinginundergroundmedia.ManymethodsareusedtoestimatetheQ-factor.WeproposeamethodtocalculatetheQ-factorbasedontheprestackQ-factorinversionandthegeneralizedS-transform.TheproposedmethodspecifiesastandardprimarywaveletandcalculatesthecumulativeQ-factors;then,itfindstheinterlaminarQ-factorsusingtherelationbetweenQandoffset(QVO)andtheDixformula.TheproposedmethodisalternativetomethodsthatcalculateinterlaminarQ-factorsafterhorizonpicking.Becausethefrequencyspectrumofeachhorizoncanbeextractedcontinuouslyona2Dtime–frequencyspectrum,themethodiscalledthecontinuousspectralratioslope(CSRS)method.ComparedwiththeotherQ-inversionmethods,themethodoffersnearlyeffortlesscomputationsandstability,andhasmathematicalandphysicalsignificance.WeusenumericalmodelingtoverifythefeasibilityofthemethodandapplyittorealdatafromanoilfieldinAhdeb,Iraq.TheresultssuggestthattheresolutionandspatialstabilityoftheQ-profileareoptimalandcontainabundantinterlaminarinformationthatisextremelyhelpfulinmakinglithologyandfluidpredictions.

简介：角度域共成像点道集（ADCIGS）是偏移速度分析和振幅随角度变化分析（AVA）的基础数据。传统Kirchhoff叠前深度偏移（KPSDM）按偏移距组织数据,能方便的输出偏移距域共成像点道集（ODCIGS）,其高效的角度道集输出是有挑战的。本文提出基于旅行时梯度场的KPSDM角道集输出方法。其核心步骤为：（1）利用任意介质中的动态规划法旅行时计算方法提供炮点和检波点的旅行时场;（2）根据旅行时场的梯度方向计算反射张角;（3）在偏移过程中抽取ADCIGS。由于本文旅行时计算方法没有射线阴影区,也没有对速度光滑性的要求,其角度道集输出在阴影区比传统射线追踪更有优势。基于该角度道集输出方法,本文发展了一种适合大规模三维地震数据的KPSDM及角道集输出的并行实现方案。其基本思想是：（1）按照炮数据来组织输入数据;（2）旅行时场的输入与单炮覆盖范围相联系以节省内存;（3）多炮数据间采用MPI并行处理,单炮深度切片之间采用OpenMp并行处理,可进一步提高内存利用率和并行力度。数值试验结果证明本文角度道集生成方法的优越性和本文实现方案的有效性。

简介：在储层发育带，通常地震波波形会变复杂，所以关联维会变小。但由于地震波是有一定带宽的信号，波形受这个带宽内所有频率成分的影响，如果直接进行关联维的计算，其结果往往不能很好地反映储层。鉴于这种情况，提出了经验模态分解（EMD）与关联维相结合的方法，利用EMD对地震波进行分解，对分解后的每一个平稳的固有模函数（IntrinsicModeFunction，IMF）分量进行关联维计算，将计算结果与有效的IMF分量进行储层预测。利用该方法对xx区三叠系中油组进行了试算，取得了较好的效果。

简介：角度域共成像点道集(ADCIGS)是偏移速度分析和振幅随角度变化分析(AVA)的基础数据。传统Kirchhoff叠前深度偏移(KPSDM)按偏移距组织数据,能方便的输出偏移距域共成像点道集(ODCIGS),其高效的角度道集输出是有挑战的。本文提出基于旅行时梯度场的KPSDM角道集输出方法。其核心步骤为:(1)利用任意介质中的动态规划法旅行时计算方法提供炮点和检波点的旅行时场;(2)根据旅行时场的梯度方向计算反射张角;(3)在偏移过程中抽取ADCIGS。由于本文旅行时计算方法没有射线阴影区,也没有对速度光滑性的要求,其角度道集输出在阴影区比传统射线追踪更有优势。基于该角度道集输出方法,本文发展了一种适合大规模三维地震数据的KPSDM及角道集输出的并行实现方案。其基本思想是:(1)按照炮数据来组织输入数据;(2)旅行时场的输入与单炮覆盖范围相联系以节省内存;(3)多炮数据间采用MPI并行处理,单炮深度切片之间采用OpenMp并行处理,可进一步提高内存利用率和并行力度。数值试验结果证明本文角度道集生成方法的优越性和本文实现方案的有效性。

简介：本文将相干成像用于WVSP数据的一次反射波和下行地表多次波的成像中,通过对叠加次数分布的分析,发现一次反射波相干成像的成像范围小但在井位置处的叠加次数高;下行地表多次波相干成像的成像范围大而在井位置处的叠加次数少。基于此,本文尝试利用一种WVSP的联合相干成像方法来解决这一问题。该方法利用WVSP的直达波与下行地表多次波和直达波与一次反射波来进行联合相干成像。这种联合相干成像方法在有效拓宽成像面积的同时增加了覆盖次数(特别是在井旁),使干涉记录中的虚假同相轴得到了有效地压制,提高了成像剖面的信噪比,使成像结果更加可靠。正演模拟得到的WVSP合成记录和野外实际采集到的WVSP数据处理结果均表明这种方法的有效性。更多还原

简介：野外空间采样密度的提高将增加室内数据分析的工作量，常规的基于地震数据的点与线分析方法有一定的局限性。本文简要说明空间子集的抽取方法，列举了正交子集和斜交子集的特点，并通过三维可视化展示了子集数据的空间特性。提出在数据处理中利用子集的时间切片：（1）分析地震道空间分布的均匀性和规则性；（2）研究面波及规则干扰的空间分布特点；（3）检测叠前数据中的异常信息；（4）监控叠前去噪的效果。实际数据的应用结果表明基于空间子集的分析方法是一种独特有效的地震数据分析方法，从另一种视角观察地震数据，可以发现其中某些新的特征，以提高处理人员对数据的洞察能力。