简介：裂缝识别对碳酸盐岩储层测井评价至关重要，传统的测井仪器由于探测深度浅（小于3米）而无法对井壁外围（大于3米）裂缝发育情况进行评价，远探测声波测井仪器采用相控阵发射、同时加大源距，有效提升了测井仪器的探测深度。但由于缺少正演模拟研究，目前对于裂缝的解释往往是基于经验而缺乏理论依据，导致很多现象难以解释。本文利用高阶有限差分方法对远探测声反射波测井裂缝识别进行了正演模拟及叠前逆时偏移成像研究，首先在理论研究的基础上构建正演模型，重点研究不同裂缝离井壁距离、裂缝张开度和倾角的响应特征；其次在单因素变化基础上提取成像区域的能量强度，分析确定出在实际地层速度有波动变化时远探测声反射波测井方法能够识别裂缝的条件；最后通过对裂缝识别的影响因素定量化分析，确定裂缝识别的最大距离、最小张开度和最小倾角，降低了裂缝识别与评价中的多解性。研究成果对远探测声反射波测井仪器的发展、数据处理方法的改进，以及后续的测井解释工作都有一定的借鉴意义。

简介：有限差分方法广泛应用于求解许多科技领域所涉及的偏微分方程，高阶显式有限差分方法通常用来提高求解精度，已经提出的高阶隐式有限差分方法和截断高阶显式有限差分方法可用来进一步提高模拟精度而不增加计算量。本文首先计算了针对常规网格上的一阶导数和二阶导数、交错网格上的一阶导数的有限差分系数，发现高阶隐式有限差分系数中存在一些小的系数。频散分析结果表明：忽略这些小的差分系数能够近似维持有限差分的精度，但是显著减小了计算量。然后，引入镜像对称边界条件来提高隐式有限差分方法的精度和稳定性，采用混合吸收边界条件来减小来自模型边界所不需要的反射。最后，给出了针对均匀和非均匀介质模型的弹性波模拟例子，表明了本文方法的优点。

简介：复杂地层曲面构建可广泛应用于石油勘探、地质建模、地质构造分析等领域中,也是这些领域数据可视化和可视分析的重要基础。现有方法在多种类型断层处理、曲面光滑性等方面存在不足。为解决地层数据分布不规则性带来的曲面光滑性问题,本文引入几何偏微分方程的曲面造型方法开展地层曲面构建;针对多种类型复杂断层存在情况下的曲面构建问题,提出三维空间与二维平面的投影互换算法,给出了基于几何偏微分方程的多种形式断层的统一处理方法;针对复杂地层曲面构建问题,建立了相应几何偏微分方程,给出演化求解算法。实际地层数据的空间曲面构建验证了本文方法计算效率高、能处理不规则的数据分布,同时能对存在断层尤其是逆掩断层的地层面进行曲面恢复。

简介：在复杂地表条件的区域,地震数据的采集和处理是一项极大的挑战。虽然可以通过静校正来消除起伏地表的影响,然而当地表高程以及近地表速度剧烈变化时,简单的垂直时移对地震波场造成的畸变会严重降低偏移成像的质量。基于射线的偏移方法可以直接在起伏地表面进行波场的延拓成像,是解决上述问题的有效手段。本文针对复杂地表条件下的高斯波束叠前深度偏移进行研究,对倾斜叠加公式进行修改,使之包含地表高程以及速度的信息,通过直接在复杂地表面进行平面波的合成,得到了一种具有更高成像精度的改进方法。首先简单介绍常规高斯波束偏移的基本原理和计算流程,并以此为基础,给出复杂地表条件下高斯波束偏移原有的实现方法以及本文的改进方法,最后通过模型和实际资料的试算验证本文方法的有效性。

简介：地震波场数值模拟中不可避免地会出现边界反射,一般采用吸收边界条件以压制人工边界反射。目前常用的分裂式完全匹配层（PML）边界条件需要在边界处进行特殊处理,尤其是在三维情况下需要将变量分裂为三个分量,增加了数值模拟的计算时间和内存占用量。与分裂式PML吸收边界条件相比,混合吸收边界条件（HABC）具有易于实现、计算量小和吸收效果好等优点,可以提高三维波动方程数值模拟的计算效率。本文将基于一阶Higdon单程波方程的混合吸收边界条件从二维计算域发展到三维,提出了适用于三维弹性波数值模拟的混合吸收边界条件。均匀模型以及复杂模型的三维数值模拟结果表明,混合吸收边界条件与传统的完全匹配层边界条件相比,具有效率高、吸收效果好的优势。

简介：初至拾取是起伏地表地震资料处理的关键环节,直接影响近地表建模的精度及静校正效果,但起伏地表资料极低的信噪比使传统的自动初至拾取算法几乎失效,而手动拾取要耗费巨大的人工和时间成本。本文研究了改进超虚干涉法,结合多道多域初至质量监控技术,实现起伏地表资料初至波自动拾取。改进超虚干涉法首次将近地表散射波纳入干涉法提高信噪比的范围,并通过折射波和散射波的线性组合使干涉法适用于起伏地表条件下任意复杂射线路径的初至波类型;提出反向干涉和多域干涉的概念,显著增强了所估计的虚源信号;使用波形反褶积滤波器较好地抑制了干涉导致的＂假事件＂的形成;采用多道多域初至质量监控技术,实现错误＂假事件＂初至的自动归位,提高了拾取初至的稳定性。本文研究的初至波自动拾取理论与技术具有突出的鲁棒性和稳定批量处理大量实际三维地震数据的能力,在中国西部山区三维地震勘探数据处理的应用效果显著,质量优于某常用商业软件。

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

简介：基于双相介质理论的AVO正演技术是储层性质描述和流体预测的有效技术手段之一,但是输入参数中基质矿物模量的准确性和双相介质模型的的合理性极大地影响双相介质AVO正演效果。因此,本文采用基于流体因子的基质矿物模量反演方法,自适应反演基质矿物体积模量。引入具有岩石物理意义的多约束条件,使得流体替换技术制作的双相介质模型具有岩石物理意义。保证获得的双相介质AVO特征反映实际地层响应,真实可靠。通过不同岩性岩样的对比分析,说明反演方法的优越性和准确性。同时LH地区实际资料应用,获得孔隙度和流体饱和度等重要岩性参数变化时双相介质AVO特征,特别是不同储层孔隙度在同一入射角对应快纵波和横波反射系数幅值的大小差异和突变角差异是分辨储层孔隙度大小的依据。

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

简介：已有的均匀磁化长方体磁场计算方法大都假设观测点位于上半无源空间，对于起伏地形条件而言，这些方法可能存在解析“奇点”。为此，本文基于地磁场基本理论，采用变量替换的积分方法，导出了改进的磁场表达式，详细讨论并有效解决了整个无源空间的所有奇点问题。相比前人的方法，其积分过程更自然、简单，最后的积分结果形式更加统一，并且不需要坐标变换即可求出无源区任意点处的磁场值，从而简化了正演过程。对比模型试验表明，新导出的磁场无解析“奇点”理论表达式是正确的，并能适应地形起伏的情况。

简介：提高地震资料的分辨率是地震数据处理流程中的重要环节,对后续的精细构造解释起到重要作用。传统的提高分辨率方法大都假设地震资料是稳态的并且噪声水平不随空间发生变化,而实际情况不满足这一假设,导致提高分辨率处理后的效果达不到预期要求。针对这一问题,本文提出了一种基于时频二次谱的提高地震资料分辨率方法。首先,文中提出了基于S变换的时频二次谱,并结合模型论述了时变子波和反射系数在时频二次谱中的特征及其可分离性;其次,依据时变子波和反射系数在时频二次谱中的特征差异,构建了二维滤波器在地震记录的时频二次谱中提取时变子波的振幅谱;再次,文中研究了噪声环境中时变提高分辨率算子设计方法,并提出了依据时频谱能量强弱相对关系自适应确定频带拓宽范围的时变提高分辨率算子设计,进行提高分辨率处理;最后,文中对该方法进行了模型和实际数据的试处理,并与传统谱模拟方法和Q补偿方法的处理结果进行了对比分析,对比结果表明：本方法不需要估计Q值,提高分辨率能力不受震源子波频带的限制,在兼顾信噪比的前提下能够充分提高不同时间局部的地震数据的分辨率。

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

简介：Gabor变换和S变换是常用的时频分析工具。根据测不准原理,它们的时频分解结果无法在时间域和频率域同时具有很高的分辨率。为了提高非平稳信号时频分解结果的分辨率,本文提出瞬时频率分布函数（IFDF）并利用它表达非平稳信号。当非平稳信号时频成分的分布满足测不准原理对信号可分辨的要求时,瞬时频率分布函数的支集和短时Fourier变换的小波脊支集是同一个集合。利用IFDF的该特征,本文提出一种迭代算法（Sparse-STFT）实现了信号的稀疏时频分解。该算法在每次迭代过程中利用残留信号的短时Fourier变换结果的脊支集更新信号的时频成分,每次迭代得到的时频成分的叠加结果即为最终的稀疏时频分解结果。文中的数值实验证明了Sparse-STFT可以有效地提高非平稳信号时频分解结果的分辨率。最后,本文将该方法应用于地震数据面波的压制中,取得了理想的处理结果。

简介：本文通过地震正演模型和面波振幅数据对比分析了天然地震、爆炸,与核试验数据,以探索改进辨别天然地震与爆炸的方法。所提方法是基于波数域和频率域的双积分变换解。用于研究的地震事件都发生在北朝鲜,其中包括2001年6月26日（39.212°N,125.383°E）和2001年10月30日（38.748°N,125.267°E）所记录的爆炸数据,2006年10月9日（41.275°N,129.095°E）所记录的一次核试验数据、以及2002年4月14日（39.207°N,125.686°E）和2002年6月7日（38.703°N,125.638°E）所记录的两次天然地震。基于这些数据通过地震波分析与理论模拟表明爆炸引起的地震与天然地震有着不同的波型特征。爆炸引起的信号特征是P波的能量较S波强,在爆炸记录上0.05-0.5Hz频率之间Rg波清晰呈现,而在然地震记录上没有。这是由于爆炸地震记录上P波是优势波,与SH成分发生了耦合。

简介：常规地震观测系统设计方法基于地下水平层状介质的假设,通常不能适应复杂构造情况。我们从控制照明的思想出发,提出了一种面向目标成像的地震观测系统设计方法,该方法需要一个由初步地震解释得到的速度模型。利用单程傅立叶有限差分波场传播算子将目标层的平面源延拓到地表,通过分析从目标层延拓到地表的波场能量的分布,可以确定目标层成像所需要的炮点或者检波点的位置。利用SEG-EAGE盐丘模型数值试算结果表明,该方法用于设计面向目标成像的特定地震采集系统。

简介：本文通过人工变换T2分布和建立管-球模型模拟法研究含水合物地层渗透率与水合物含量之间的关系。首先,在渗透率的模拟试验中,我们改变了束缚水与可动水的比例、总孔隙度以及与之关联的T2分布。试验结果表明,相对渗透率与水合物含量之间的关系受到这些因素的制约。随后,我们用管-球模型表示水合物生长的孔隙空间,并把水合物的生长过程看成是向孔隙空间随机扔小球的过程。在此过程中,采用两种方法计算渗透率,一是Schlumberger＇sT2公式（即SDR模型）,二是Darcy定律与Poiseuille流动方程相结合的方法。前人的实验研究表明,在一定的水合物含量范围内,渗透率基本保持不变。以此为参考,我们将计算结果与之进行比较。我们发现,采用SDR模型时,渗透率的数值模拟曲线与Masuda模型N=15时的结果相近。而采用Darcy定律时,渗透率模拟值较高,但与实验结果的趋势相一致,都会出现渗透率的平直阶段。尤其,当水合物晶体在孔隙体内优先生成时,优先的概率越高,渗透率的平直范围越大。

简介：海底电缆双检采集是利用水检和陆检对海水层鸣震的响应极性相反,通过水、陆检数据的合并处理来消除海水层的鸣震,进而拓展了数据的频带范围,提高了资料的保真性和分辨率,而得到了广泛应用。但海底电缆双检资料中的陆检数据,由于受洋流、海底地形和施工工艺等的影响,很难保证其与海底的良好耦合,因而受耦合系统传输函数的影响比较大,造成水、陆检数据在振幅、频率、相位等方面的差异,降低了陆检数据的信噪比,影响了双检地震数据的合并处理效果。而海底电缆双检中的水检资料检测的是海水压力的变化,与海底不存在耦合问题,因而水检数据信噪比往往比较高。本文首先给出了陆检与海底耦合系统传输函数的数学表达式,然后利用水检数据作为约束条件估算出了陆检与海底耦合系统的传输函数,并利用估算出的传输函数对陆检数据进行了耦合校正处理,解决了陆检与海底的耦合效应对陆检数据的振幅、相位等的影响,提高了陆检数据的信噪比,取得了较好的应用效果。

简介：在地震频段内,中观尺度局域流是引起孔隙岩层速度频散与衰减的主要原因,研究中观尺度孔隙中局域流的过渡频率及其对应的衰减（逆品质因子）峰值频率之间的关系,有助于深入分析孔隙介质内部构造的细节。基于Biot孔隙弹性理论得出时间域内局域流流体通量,引入傅立叶变换得到频域上的流体通量。为避免不同单元体之间局域流的相互影响,建立双相孔隙介质单重单元体模型,并且在流体通量基础上选择含水层的过渡频率来等价局域流的过渡频率。此外,针对不同孔隙介质参数,分析了局域流整体过渡频率及相应的衰减峰值频率的变化特征。为了了解边界条件对计算结果造成的影响,建立了多重单元体模型,引出局部局域流过渡频率及相应的衰减峰值频率。数值模拟结果表明,在低饱和度状态下,两种频率的变化趋势相同,但变化幅度存在较大差异;在高饱和度状态下,两种频率的变化趋势与变化幅度都很接近。

简介：3D地震波走时计算是偏移、反偏移、层析等诸多地震勘探技术中的重要中间步骤。快速推进法计算3D地震波走时具有高效率、稳定性及适应能力强的特点，但快速推进法在震源附近区域的计算精度不高，降低了整个走时算法的计算精度。本文提出了一种联合3D走时计算方法来解决这一问题。该方法在震源附近小范围内使用计算精度较高的波前构建法计算走时，在剩余区域使用快速推进法计算走时，由于模型中绝大多数网格节点走时是通过快速推进法计算的，故新方法保留了快速推进法高效的特点，同时由于震源附近网格节点走时精度的提高，整个新算法的计算精度相对于快速推进法而言有了较大的改善。文中通过数值分析对上述结论进行了验证并使用三维岩丘模型验证了新方法的稳定性和适应能力。

简介：致密储层与围岩波阻抗差异小，常规储层反演方法难以有效预测出有利区。中国西部致密储层非均质性强，常规储层反演预测难度大，本文提出了基于Xu-White模型的叠前储层预测方法。首先修正Xu-White模型的试验参数，直至该模型计算出来的纵、横波速度及密度与实际测井得到的数值基本吻合，然后用修正的模型计算不同储层物性、岩性、含油性下的纵、横波速度、密度及泊松比，进而建立不同的砂泥岩互层地质模型。其次，采用Zoeppritz方程近似式进行叠前正演，分析不同储层物性、岩性、含油性时的叠前道集响应特征，在对中国西部Z区实际叠前道集优化处理的基础上，选择其中1种对岩性、含油性、物性最敏感的弹性参数，对其有利储层分布区进行预测，其预测结果具有较高的精度。