简介：地震正演模拟技术是研究地震波在复杂介质中传播规律的有效途经，尤其在地质构造及其复杂的中国西部地区，其意义更是重大。本文介绍了两种新的正演模拟技术：有限元有限差分方法（FE—FDM）和任意精细积分方法（ADPI），并结合实例分析来验证FE—FDM和ADPI算法的实际效果，结果证明这两种方法能够有效地适用于复杂介质下的地震波传播性质的研究。

简介：噪声衰减是探地雷达信号处理中的关键问题之一。当探测目标埋藏深度比较浅时,其反射信号与直耦信号和地面回波信号相互重叠,直接影响目标反射波到达时刻的检测及目标的正确定位。针对这个问题,本文提出了一种基于Curvelet变换的噪声衰减方法。通过对理论数值模拟数据和实测数据的处理,以及与平均消去法和二维连续小波该方法处理结果的对比,验证了该方法的可行性和有效性。处理结果显示,该方法不仅可以去除背景噪声、同时可以衰减倾斜相关的相干干扰和数据中的随机噪声。与二维连续小波变换方法相比有更高的计算效率。

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

简介：提出一种对整幅高光谱图像的稀疏表示结果进行直接显示的方法,图中不仅包含了稀疏表示中保留的光谱信息,还可显示整幅图像的空间信息。稀疏表示后,将字典中的各有效原子根据光谱特性选择颜色标签,之后根据稀疏系数进行混合颜色显示,此时的图像能够同时满足可分性及距离保持特性。针对局部地物时,提出的单像素混合阵列表示法及改进的裂片纹理技术能够直观且完整的显示出每个像元的具体组成情况,还能够根据所生成图像中的信息对原始HSI进行重建,进而提高数据的利用率。该模型不仅能够良好地显示地物的空间特性,同时能够显示稀疏系数的组成,同时单像素混合阵列表示法及裂片纹理技术弥补了混合像素彩色显示中颜色表达混乱的弊端。对真实地物数据进行实验,结果证明该模型产生的彩色图像具有良好的视觉效果及可分性,满足距离保持特性。

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

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

简介：海洋地震勘探过程中，由于采集设备的老化或电源的不稳定而造成的漏电，在地震记录表现为强噪音干扰，利用常规噪音衰减方法处理此类强噪音效果并不理想。鉴于强噪音在统计学上具有相同的特性，本文在基于峰度的盲分离（blindsourceseparation，BSS）算法研究基础上，推导出一种基于多用户峰度（multiuserkurtosis，MUK）准则的噪音衰减算法来估计地震记录中具有相同统计特性的强噪音，并将其从地震记录中分离，从而达到衰减强噪音的目的。模型试验与实际资料的处理表明：该方法能够在好的衰减海洋地震勘探记录中的强噪音，保留了更多的有效信息，提高海洋地震数据的信噪比，具有可行性和应用前景。

简介：噪声衰减是探地雷达信号处理中的关键问题之一。当探测目标埋藏深度比较浅时,其反射信号与直耦信号和地面回波信号相互重叠,直接影响目标反射波到达时刻的检测及目标的正确定位。针对这个问题,本文提出了一种基于Curvelet变换的噪声衰减方法。通过对理论数值模拟数据和实测数据的处理,以及与平均消去法和二维连续小波该方法处理结果的对比,验证了该方法的可行性和有效性。处理结果显示,该方法不仅可以去除背景噪声、同时可以衰减倾斜相关的相干干扰和数据中的随机噪声。与二维连续小波变换方法相比有更高的计算效率。更多还原

简介：为了解决海底起伏变化对地震波场的影响问题,本文提出将（x-z）域中的曲网格映射成（ξ,η）域中的矩形网格,推导出（ξ,η）域中的二维标量声波方程,根据推导出来的波动方程采用逆时有限差分法将海面上采集到的地震波场在（ξ,η）域中向下延拓至海底面,延拓时采用海水的速度,然后采用顺时有限差分法将延拓后的地震波场再反延拓到海面上,延拓时采用海底面以下地层的速度,从而消除了海底起伏带来的负面影响。模型及实际地震资料的计算分析表明该方法不但能够校正由于海底起伏所引起的海底面下地层反射波场的不连续性还能够校正由于海底起伏所引起的地震波的动力学特征的变化。对延拓前后的地震波场进行速度反演,延拓后反演的地层速度比延拓前反演的地层速度的精度提高很多,延拓前后地震波场的叠加剖面对比表明该延拓方法能够明显提高地震波场的成像质量。

简介：波动方程有限差分法是地震数值模拟中的一种重要的方法，对理解和分析地震传播规律、分析地震属性和解释地震资料有着非常重要的意义。但是有限差分法由于其离散化的思想，产生了不稳定性。精细积分法在有限差分法的基础上，在时间域采用解析解的表达形式，在空间域保留任意差分格式，发展成为半解析的数值方法。本文结合并发展了以往学者的成果，推导了任意精细积分法的三维弹性波正演模拟计算公式，并对其稳定性进行了数值分析。在计算实例中，实现了精细积分法二维和三维弹性波模型的地震正演模拟，对计算结果的分析表明，精细积分法反射信号走时准确，稳定性好，弹性波场相较于声波波场，弹性波波场成分更为丰富，包含了更多波型成分（PP-和PS-反射波、透射波和绕射波），这对实际地震资料的解释和储层分析有重要的意义。实践证明，该方法可直接应用到弹性波的地质模型的数值模拟中。

简介：Inthispaper,wetheoreticallyandnumericallystudyacombinedmonopole–dipolemeasurementmodetoshowitscapabilitytoovercometheissuesencounteredinconventionalsingle-wellimaging,i.e.,thelowsignal-to-noiseratioofthereflectionsandazimuthambiguity.First,theazimuthambiguity,whichexistsextensivelyinconventionalsingle-wellimaging,issolvedwithanimprovedimagingprocedureusingcombinedmonopole–dipoleloggingdatainadditiontoconventionalloggingdata.Furthermore,wedemonstratethatthedirectwavespropagatingalongtheboreholeswithstrongenergy,canbeeffectivelyeliminatedwiththeproposedcombinedmonopole–dipolemeasurementmode.Thereflectionsarethereforepredominantinthecombinedmonopole–dipoledataevenbeforethesignalsarefiltered;thus,thereflections’arrivaltimesineachreceiverareidentified,whichmayhelpminimizethedifficultiesinfilteringconventionalloggingdata.Theoptimizedprocessingflowofthecombinedmeasurementmodeloggingimageisgiveninthispaper.Theproposedcombinedmonopole–dipolemeasurementmodemayimprovetheaccuracyofsingle-wellimaging.

简介：相干属性是检测地下介质不连续性特征的主要参数之一，但倾斜地层会对相干属性产生干扰，故在计算相干属性时需要消除地层倾角的影响。传统的相干属性计算一般使用时域倾角扫描法来消除地层倾角的影响，但这种方法的精度受到采样间隔等因素的限制，并不能很好地消除地层倾角的影响。根据振幅谱具有不受地层倾角影响的特点，本文提出将预先定义的子数据体内的地震数据变换为振幅谱，再利用振幅谱构造协方差矩阵，协方差矩阵最大特征根与矩阵迹的比值作为子数据体中心点的相干值。利用振幅谱计算相干属性的方法能够较好地消除地层倾角的影响，而且可以使用不同频段的振幅谱计算相干体，具备多尺度相干体算法的特性。模型数据和实际地震数据的应用结果表明，振幅谱相干属性可以更好地消除地层倾角对相干属性的影响。

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

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

简介：尾波是持续时间比主要波型长得多的多重散射波，它包含了直达波之外的部分有用信号。在油气田的开采过程中，由储层物性参数的微小变化而导致的储层速度的微小变化对初至的影响很小，无法通过初至变化直接观测，但是因为尾波为多重散射波，故储层速度的微小变化会在尾波的传播过程中被放大，因此利用尾波可以观测到这种变化。本文通过实验与数值模拟研究了尾波干涉方法在储层微小时移差异监测方面的有效性。合成地震数据是基于部分MarmousiII模型采用有限差分波动方程正演方法计算获得，为了模拟时移地震中的储层波速变化，在模型中选取波速发生微小变化的目标区，计算速度扰动前后的合成地震数据，进而观测尾波变化。研究表明利用尾波干涉可检测储层物性的微小变化，为时移地震信号监测提供新的方法和手段，从而提高时移监测的准确度，也为开发生产工作提供可靠的指导。

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

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

简介：本文提出利用在不同深度处倾斜放置气枪子阵,对不同深度处单枪或相干枪延时激发的方法压制震源端鬼波,从而改善震源子波分辨率。采用延时激发的方法可以构造同相叠加的首波以及被相对压制掉的鬼波。为了合理评价倾斜震源对于鬼波压制的效果,提出了以实际子波和不受鬼波影响的期望子波频谱之间的归一化平方误差为标准,以深度间隔、气枪组合次序、子阵数目等为变量的最优化问题。针对680cu.in经典枪阵研究表明,在综合考虑去鬼波效果和野外施工可操作性,相邻单枪之间的深度间隔为1m-1.5为较为合理的选择。多个子阵在进行组合时,需要对组合效果进行前期模拟确定最佳深度组合,过多的子阵数目由于会引人大量陷波频率,反而会影响鬼波压制效果,在本次实例研究中,2-3个子阵是较为合理的选择。气枪组合次序会一定程度上影响鬼波压制效果,但不同次序结果差异可以通过简单匹配滤波消除。方向性对比表明多深度倾斜气枪震源能够有效地补偿消除陷波效应的影响,同时改善能量传播的稳定性。

简介：瑞利波勘探主要是建立在弹性介质分层半空间模型基础上的。当实际地层中包含孔隙介质层时，需要将孔隙介质简化为弹性介质来进行分析和处理。这种简化处理究竟对瑞利波勘探造成怎样的影响是本文所研究的问题。本文基于弹性介质与孔隙介质共同构成的分层半空间模型，首先推导了分层半空间中两种介质处于不同相对位置时的瑞利波频散方程，解决了不同阶数矩阵之间的变量传递问题然后，针对传递矩阵法在求解频散函数时可能出现的溢出问题，给出了一种可以有效提高计算范围的解决方案；同时，提出了一套新的数值算法用于复频散方程的快速求解。数值计算结果表明：当孔隙介质位于半空间表面时对低频瑞利波频散特性的影响最为显著，而在其它情况下对瑞利波频散特性的影响相对较小。

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