简介：常规多波联合反演采用Zoeppritz方程的近似式构建正演方程,反演过程中需要假定背景纵横波速度比为常数,其反演精度不高,稳定性不好。本文提出了一种基于精确Zoeppritz方程的多波联合反演方法,结合贝叶斯方法进行广义线性反演。本方法基于精确Zoeppritz方程构建正演方程,避免了近似式反演在大角度时引起的误差;利用贝叶斯方法引入模型参数的先验分布信息,作为反演的正则化项,降低了反演的不适定性;反演目标函数中引入低频软约束,稳定了反演低频结果,提高了反演的鲁棒性;在求解反演目标函数时,利用快速算法,降低了反演的运算量。经过模型试算,证明了该方法的优越性和抗噪性;并在实际资料的应用中证明了该方法的实用性和有效性。

简介：本文以中观孔隙结构的White模型为基础，构建了部分饱和孔隙介质模型，利用Biot方程的建立思路和Johnson推导的体变模量，推导了部分饱和孔隙介质中的纵波方程，并以平面波为例，求取了方程的衰减系数，分析了地震勘探频带范围内地震波的衰减特性。结果表明：在部分饱和孔隙介质中，地震波在低频段也会发生明显的衰减和频散现象，频率越大，衰减越大；且第二纵波的衰减比第一纵波更为明显；这一结论弥补了Biot理论在描述地震勘探频带范围内波的衰减现象的不足。文中还研究了孔隙度、饱和度和模型内径尺寸对第纵波衰减特性的影响机理，主要表现在在地震勘探频带范围内，波的衰减随孔隙度的增大而增大，随含油气饱和度的增大而减小，当孔隙内径尺寸小于二分之一外径尺寸时，波的衰减随内径尺寸的增大而增大，当内径尺寸大于二分之一外径尺寸时，波的衰减随内径尺寸增大而减小。

简介：作者基于数学检波器和等时叠加原理,实现了复杂地表的单程波动方程地震叠前正演模拟。该方法采用虚拟的数学检波器接收地下的反射地震信号,灵活地将接收点布置在地表的任何地方,从而满足地表起伏的要求。此外,根据等时叠加原理,该方法采用单程波动方程进行波场延拓和成像,计算简单快速。通过复杂正断层的数值模拟,得到了高信噪比的共炮集地震记录,并采用适用于起伏地形的深度偏移方法对该共炮集地震记录进行了叠前深度偏移,较好地实现了地震波的偏移归位,从而证明了这里提出的起伏地表的单程波动方程地震叠前正演方法是正确和有效的。

简介：用基于张量格林函数的体积分方程法对三维异常体进行瞬变电磁响应的正演模拟,首先在频率域内计算电磁场分量的频率域响应,然后利用快速数字滤波技术将计算结果转换到时间域。设计和计算了水平电偶极子源激发下层状水平地层模型背景下的常见地形如山谷、山峰地形的模型,并考察分别把源和接收器放于这些地形中的瞬变电磁场响应,详细分析了这些地形对长偏移距瞬变电磁测深（LOTEM）的影响。结果表明,山谷和山峰地形对LOTEM的结果均有不同程度的影响。当电偶极子源放在山谷谷底时,地形对观测异常场的畸变非常严重;当接收器放在山谷中时,接收器处地形的影响强烈但该影响在空间和时间上只是局部的。总体来讲,不论山峰地形位于何处,其对LOTEM的影响相对较小。当地形处于发射源与接收器之间时,地形对LOTEM的影响非常小,表明在进行LOTEM勘探时,选择发射源的放置比接收器的位置更加重要,野外勘探是尽量把发射源选择在开阔的平坦位置。

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

简介：在地震波场数值模拟中，为了消除由人为边界产生的边界反射，需要引进边界吸收条件。本文从声波方程完全匹配层吸收边界的经典方法出发，基于高斯函数任意阶光滑的特点，提出了一种高斯型衰减因子，分析比较该因子与一般衰减因子的性质，并基于均匀与层状速度模型分别进行了数值模拟计算。数值结果显示，当选择相同的PML边界吸收层层数时，高斯型衰减因子的吸收效果明显优于一般的衰减因子，边界反射更少，信噪比更高；对比最近提出的正弦型衰减因子，在信噪比接近的情况下，高斯型衰减因子所需的PML吸收层层数更少。

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

简介：由于地震数据中包含的噪声在不同频率或者频带数据中的分布强度存在差异,使得全频带数据上进行的噪声衰减处理改变了地震反射波信号的动力学特征,干扰后期的地震资料解释、储层预测、油气检测等问题,提出边界和振幅特性保持自适应噪声衰减方法。首先应用小波包变换对全频带地震数据进行多频段划分,然后对分频段数据进行非线性各向异性倾角导向边界保持自适应滤波处理。在该方法中,由结构张量计算的扩散张量实现自适应地确定平滑滤波方向,加入的不连续结构置信度量和不连续性算子自适应地控制不连续结构特征的保持程度,引入的去相关滤波迭代停止准则自适应地确定滤波迭代次数。这些参数的引入具有减少处理人员的干预和人为的主观性,且执行简单的特点。对合成地震记录和实际地震记录处理结果表明,提议的方法能够自适应地衰减地震数据中噪声,同时既能保持地震反射波中有效的不连续性信息,也能有效地保持有效信号的频率分布规律。能够为后期的地震资料解释和分析提供高品质的基础数据。

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

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

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

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

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

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

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

简介：渗透率评价是储层评价中的一个复杂问题，传统测井方法难以对储层渗透率参数进行直接和有效的评估。储层渗透率对钻井过程中的泥浆滤液侵入有着较大的影响，因此本文提出一种利用泥浆侵入效应来评价储层渗透率的数学模型和方法。首先构造含泥饼的泥浆侵入数值模型，然后以达西渗流理论为基础导出泥浆侵人深度和储层渗透率的近似数学关系，以此可以利用侵入深度测量值来估算储层渗透率。对孔隙度、渗透率和含水饱和度三个主要储层参数的敏感性分析，发现该方法适用于受到泥浆侵入影响的低孔隙度、低渗透率的油层或油水同层。采用现场测井和取芯数据进行数值模拟计算，结果表明估算出来的渗透率曲线和预设的渗透率曲线吻合较好，证明该方法具有一定的可行性。

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

简介：由于地震数据中包含的噪声在不同频率或者频带数据中的分布强度存在差异,使得全频带数据上进行的噪声衰减处理改变了地震反射波信号的动力学特征,干扰后期的地震资料解释、储层预测、油气检测等问题,提出边界和振幅特性保持自适应噪声衰减方法。首先应用小波包变换对全频带地震数据进行多频段划分,然后对分频段数据进行非线性各向异性倾角导向边界保持自适应滤波处理。在该方法中,由结构张量计算的扩散张量实现自适应地确定平滑滤波方向,加入的不连续结构置信度量和不连续性算子自适应地控制不连续结构特征的保持程度,引入的去相关滤波迭代停止准则自适应地确定滤波迭代次数。这些参数的引入具有减少处理人员的干预和人为的主观性,且执行简单的特点。对合成地震记录和实际地震记录处理结果表明,提议的方法能够自适应地衰减地震数据中噪声,同时既能保持地震反射波中有效的不连续性信息,也能有效地保持有效信号的频率分布规律。能够为后期的地震资料解释和分析提供高品质的基础数据。

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