简介:常规多波联合反演采用Zoeppritz方程的近似式构建正演方程,反演过程中需要假定背景纵横波速度比为常数,其反演精度不高,稳定性不好。本文提出了一种基于精确Zoeppritz方程的多波联合反演方法,结合贝叶斯方法进行广义线性反演。本方法基于精确Zoeppritz方程构建正演方程,避免了近似式反演在大角度时引起的误差;利用贝叶斯方法引入模型参数的先验分布信息,作为反演的正则化项,降低了反演的不适定性;反演目标函数中引入低频软约束,稳定了反演低频结果,提高了反演的鲁棒性;在求解反演目标函数时,利用快速算法,降低了反演的运算量。经过模型试算,证明了该方法的优越性和抗噪性;并在实际资料的应用中证明了该方法的实用性和有效性。
简介:频率-波数域单程波算子能高效地模拟地震波在复杂介质中的传播,但是在描述波的大角度传播和速度横向扰动变化较大介质中传播的问题时仍然存在一定误差。这类误差是由于对单平方根算子使用Taylor展开式的近似程度不足所造成。为了进一步提高泰勒展开式的精确性,本文提出一种利用粒子群智能算法优化级数展开系数的高阶广义屏算子对单平方根算子的展开级数进行优化处理。新的偏移算法能在保持单程波偏移算法高效的前提下进一步提高偏移算子在大角度的成像精度和对强横向速度变化介质的适应性。通过脉冲响应实验,验证了基于粒子群算法优化级数的高阶广义屏算子能够提高常规的高阶广义屏算子的成像精度和成像角度。根据对二维SEG/EAGE盐丘模型的成像处理,基于粒子群算法优化级数的高阶广义屏算子对盐丘下面的断层取得了更高质量的成像,说明粒子群优化级数的高阶广义屏算子比常规的高阶广义屏算子具有更好的横向速度适应性。为了检验本文所提算法对实际资料的处理能力,我们利用常规的偏移处理技术和本文所提算法对一条海上二维数据进行了偏移成像处理,对比分析成像剖面发现本文所提算法描述了更加清晰的层位信息和更高质量的偏移剖面。本文所提算法能有效提高高阶广义屏偏移在广角度成像的能力,具有一定实际应用价值。
简介:从大量的地震属性中提取最能反映地质特征的综合属性是储层预测技术的关键,通常选用降维方法来优选属性。目前应用最为广泛的线性降维方法。但是,由于地震属性与地质特征的关系通常是非线性的,基于线性变换的地震属性降维优化方法不能充分地反映这种非线性关系,降低了储层预测的精度。流形学习是一种新的非线性学习方法,它是通过保持数据局部结构的方式将高维数据投影到低维空间,挖掘和发现隐藏在数据中的内在特征与规律性,开拓了地震属性降维优化研究的新领域。本文首次实现了3D地震数据的层问属性特征提取,讨论了LLE方法及其关键技术,并以奥陶系礁滩相储层实例说明LLE和PCA两种方法降维及聚类的不同效果。理论模型分析和实例应用表明:LLE较好地保持了数据本身的原始结构;提取的综合属性和聚类相图较好地刻画了沉积相带、储层和流体的特征。这说明流形学习具有更好的特征提取性能。
简介:由于地震数据中包含的噪声在不同频率或者频带数据中的分布强度存在差异,使得全频带数据上进行的噪声衰减处理改变了地震反射波信号的动力学特征,干扰后期的地震资料解释、储层预测、油气检测等问题,提出边界和振幅特性保持自适应噪声衰减方法。首先应用小波包变换对全频带地震数据进行多频段划分,然后对分频段数据进行非线性各向异性倾角导向边界保持自适应滤波处理。在该方法中,由结构张量计算的扩散张量实现自适应地确定平滑滤波方向,加入的不连续结构置信度量和不连续性算子自适应地控制不连续结构特征的保持程度,引入的去相关滤波迭代停止准则自适应地确定滤波迭代次数。这些参数的引入具有减少处理人员的干预和人为的主观性,且执行简单的特点。对合成地震记录和实际地震记录处理结果表明,提议的方法能够自适应地衰减地震数据中噪声,同时既能保持地震反射波中有效的不连续性信息,也能有效地保持有效信号的频率分布规律。能够为后期的地震资料解释和分析提供高品质的基础数据。
简介:噪声衰减是探地雷达信号处理中的关键问题之一。当探测目标埋藏深度比较浅时,其反射信号与直耦信号和地面回波信号相互重叠,直接影响目标反射波到达时刻的检测及目标的正确定位。针对这个问题,本文提出了一种基于Curvelet变换的噪声衰减方法。通过对理论数值模拟数据和实测数据的处理,以及与平均消去法和二维连续小波该方法处理结果的对比,验证了该方法的可行性和有效性。处理结果显示,该方法不仅可以去除背景噪声、同时可以衰减倾斜相关的相干干扰和数据中的随机噪声。与二维连续小波变换方法相比有更高的计算效率。
简介:提出一种对整幅高光谱图像的稀疏表示结果进行直接显示的方法,图中不仅包含了稀疏表示中保留的光谱信息,还可显示整幅图像的空间信息。稀疏表示后,将字典中的各有效原子根据光谱特性选择颜色标签,之后根据稀疏系数进行混合颜色显示,此时的图像能够同时满足可分性及距离保持特性。针对局部地物时,提出的单像素混合阵列表示法及改进的裂片纹理技术能够直观且完整的显示出每个像元的具体组成情况,还能够根据所生成图像中的信息对原始HSI进行重建,进而提高数据的利用率。该模型不仅能够良好地显示地物的空间特性,同时能够显示稀疏系数的组成,同时单像素混合阵列表示法及裂片纹理技术弥补了混合像素彩色显示中颜色表达混乱的弊端。对真实地物数据进行实验,结果证明该模型产生的彩色图像具有良好的视觉效果及可分性,满足距离保持特性。
简介:针对转换波各向异性叠前时间偏移实践中的四参数速度模型估计中参数间相互影响、难于准确的确定,本文将简化的转换波两参数动校正方程从叠加速度分析拓展到叠前时间偏移速度模型修正,形成了基于两参数动校正方程的新的四参数法转换波叠前时间偏移速度修正新方法和流程。在该方法和流程中,先进行转换波两参数叠加速度分析,获得转换波叠前时间偏移初始的速度和各向异性参数,然后通过分析偏移后的共成像点道集中的剩余时差修正速度和各向异性参数。其中叠前时间偏移速度模型的垂直速度比仍是需要利用纵波和转换波叠加剖面,层位标定后,利用相关的方法获得。有效速度比需结合纵波速度分析,利用Thomsen(1999)公式计算得到初始模型,最终依靠百分比扫描偏移处理获得最终的有效速度比模型。该方法简化了转换波高质量成像的速度估计方法,减小了多个参数估计的不确定性,在实际应用中也取得了较好的成像效果。
简介:海洋地震勘探过程中,由于采集设备的老化或电源的不稳定而造成的漏电,在地震记录表现为强噪音干扰,利用常规噪音衰减方法处理此类强噪音效果并不理想。鉴于强噪音在统计学上具有相同的特性,本文在基于峰度的盲分离(blindsourceseparation,BSS)算法研究基础上,推导出一种基于多用户峰度(multiuserkurtosis,MUK)准则的噪音衰减算法来估计地震记录中具有相同统计特性的强噪音,并将其从地震记录中分离,从而达到衰减强噪音的目的。模型试验与实际资料的处理表明:该方法能够在好的衰减海洋地震勘探记录中的强噪音,保留了更多的有效信息,提高海洋地震数据的信噪比,具有可行性和应用前景。
简介:噪声衰减是探地雷达信号处理中的关键问题之一。当探测目标埋藏深度比较浅时,其反射信号与直耦信号和地面回波信号相互重叠,直接影响目标反射波到达时刻的检测及目标的正确定位。针对这个问题,本文提出了一种基于Curvelet变换的噪声衰减方法。通过对理论数值模拟数据和实测数据的处理,以及与平均消去法和二维连续小波该方法处理结果的对比,验证了该方法的可行性和有效性。处理结果显示,该方法不仅可以去除背景噪声、同时可以衰减倾斜相关的相干干扰和数据中的随机噪声。与二维连续小波变换方法相比有更高的计算效率。更多还原
简介:渗透率评价是储层评价中的一个复杂问题,传统测井方法难以对储层渗透率参数进行直接和有效的评估。储层渗透率对钻井过程中的泥浆滤液侵入有着较大的影响,因此本文提出一种利用泥浆侵入效应来评价储层渗透率的数学模型和方法。首先构造含泥饼的泥浆侵入数值模型,然后以达西渗流理论为基础导出泥浆侵人深度和储层渗透率的近似数学关系,以此可以利用侵入深度测量值来估算储层渗透率。对孔隙度、渗透率和含水饱和度三个主要储层参数的敏感性分析,发现该方法适用于受到泥浆侵入影响的低孔隙度、低渗透率的油层或油水同层。采用现场测井和取芯数据进行数值模拟计算,结果表明估算出来的渗透率曲线和预设的渗透率曲线吻合较好,证明该方法具有一定的可行性。
简介:波动方程有限差分法是地震数值模拟中的一种重要的方法,对理解和分析地震传播规律、分析地震属性和解释地震资料有着非常重要的意义。但是有限差分法由于其离散化的思想,产生了不稳定性。精细积分法在有限差分法的基础上,在时间域采用解析解的表达形式,在空间域保留任意差分格式,发展成为半解析的数值方法。本文结合并发展了以往学者的成果,推导了任意精细积分法的三维弹性波正演模拟计算公式,并对其稳定性进行了数值分析。在计算实例中,实现了精细积分法二维和三维弹性波模型的地震正演模拟,对计算结果的分析表明,精细积分法反射信号走时准确,稳定性好,弹性波场相较于声波波场,弹性波波场成分更为丰富,包含了更多波型成分(PP-和PS-反射波、透射波和绕射波),这对实际地震资料的解释和储层分析有重要的意义。实践证明,该方法可直接应用到弹性波的地质模型的数值模拟中。
简介:由于地震数据中包含的噪声在不同频率或者频带数据中的分布强度存在差异,使得全频带数据上进行的噪声衰减处理改变了地震反射波信号的动力学特征,干扰后期的地震资料解释、储层预测、油气检测等问题,提出边界和振幅特性保持自适应噪声衰减方法。首先应用小波包变换对全频带地震数据进行多频段划分,然后对分频段数据进行非线性各向异性倾角导向边界保持自适应滤波处理。在该方法中,由结构张量计算的扩散张量实现自适应地确定平滑滤波方向,加入的不连续结构置信度量和不连续性算子自适应地控制不连续结构特征的保持程度,引入的去相关滤波迭代停止准则自适应地确定滤波迭代次数。这些参数的引入具有减少处理人员的干预和人为的主观性,且执行简单的特点。对合成地震记录和实际地震记录处理结果表明,提议的方法能够自适应地衰减地震数据中噪声,同时既能保持地震反射波中有效的不连续性信息,也能有效地保持有效信号的频率分布规律。能够为后期的地震资料解释和分析提供高品质的基础数据。
简介:相干属性是检测地下介质不连续性特征的主要参数之一,但倾斜地层会对相干属性产生干扰,故在计算相干属性时需要消除地层倾角的影响。传统的相干属性计算一般使用时域倾角扫描法来消除地层倾角的影响,但这种方法的精度受到采样间隔等因素的限制,并不能很好地消除地层倾角的影响。根据振幅谱具有不受地层倾角影响的特点,本文提出将预先定义的子数据体内的地震数据变换为振幅谱,再利用振幅谱构造协方差矩阵,协方差矩阵最大特征根与矩阵迹的比值作为子数据体中心点的相干值。利用振幅谱计算相干属性的方法能够较好地消除地层倾角的影响,而且可以使用不同频段的振幅谱计算相干体,具备多尺度相干体算法的特性。模型数据和实际地震数据的应用结果表明,振幅谱相干属性可以更好地消除地层倾角对相干属性的影响。
简介:复杂地层曲面构建可广泛应用于石油勘探、地质建模、地质构造分析等领域中,也是这些领域数据可视化和可视分析的重要基础。现有方法在多种类型断层处理、曲面光滑性等方面存在不足。为解决地层数据分布不规则性带来的曲面光滑性问题,本文引入几何偏微分方程的曲面造型方法开展地层曲面构建;针对多种类型复杂断层存在情况下的曲面构建问题,提出三维空间与二维平面的投影互换算法,给出了基于几何偏微分方程的多种形式断层的统一处理方法;针对复杂地层曲面构建问题,建立了相应几何偏微分方程,给出演化求解算法。实际地层数据的空间曲面构建验证了本文方法计算效率高、能处理不规则的数据分布,同时能对存在断层尤其是逆掩断层的地层面进行曲面恢复。
简介:本文提出基于原始含表层多次波数据实现叠前共炮集地震数据插值。相对于利用相邻道的信息变换或外推插值用于缺失的地震数据重建,本文方法利用表层多次波数据互相关构建准一次波,将蕴含在表层多次波数据中的,而在采集记录中表现为缺失的近炮检距信息提取出来,并在滑动时间空间窗内采用最小二乘匹配滤波和均方根振幅校正方法进行准一次波校正而后用于数据插值重建。本文方法适用于表层多次波比较发育,同时又存在数据缺失尤其是近炮检距数据缺失情况。方法易于实现,不需多次波和一次波的提取,利用多次波中蕴含的信息实现缺失的地震数据弥补,为含有表层多次波的数据进行近炮检距地震信息的插值重建提供了一个很好的思路。
简介:本文提出利用在不同深度处倾斜放置气枪子阵,对不同深度处单枪或相干枪延时激发的方法压制震源端鬼波,从而改善震源子波分辨率。采用延时激发的方法可以构造同相叠加的首波以及被相对压制掉的鬼波。为了合理评价倾斜震源对于鬼波压制的效果,提出了以实际子波和不受鬼波影响的期望子波频谱之间的归一化平方误差为标准,以深度间隔、气枪组合次序、子阵数目等为变量的最优化问题。针对680cu.in经典枪阵研究表明,在综合考虑去鬼波效果和野外施工可操作性,相邻单枪之间的深度间隔为1m-1.5为较为合理的选择。多个子阵在进行组合时,需要对组合效果进行前期模拟确定最佳深度组合,过多的子阵数目由于会引人大量陷波频率,反而会影响鬼波压制效果,在本次实例研究中,2-3个子阵是较为合理的选择。气枪组合次序会一定程度上影响鬼波压制效果,但不同次序结果差异可以通过简单匹配滤波消除。方向性对比表明多深度倾斜气枪震源能够有效地补偿消除陷波效应的影响,同时改善能量传播的稳定性。
简介:随着煤矿的深部开采,多层积水采空区探测成为煤田水文物探工作的新内容。中心回线瞬变电磁(简称TEM)法因场结构与地层的耦合关系,对低阻层探测有利。但低阻层的屏蔽作用不仅使得探测同样的深度需要更长的观测时间,而且还会减弱下伏地层的异常响应。本文通过直接时域数值模拟和水平分层大地的模型正演,估算了探测目的层所需要的时间长度,根据噪声对观测数据造成的影响,给出了多层积水采空区可分辨的标准。山西大同达子沟煤矿水文勘探中获得的TEM实测曲线,表现了对多层积水采空区的探测能力。在实测曲线的定量反演解释中,利用电测井资料作为初始参数进行定量反演,约束了煤系薄层的等值性。所推断的三层和两层积水采空区,为钻孔所证实。研究结果表明,当观测时间有足够的长度、下伏地层的异常显示大于观测误差,中心回线TEM法探测多层采空区积水是可行性的。
简介:将偏移后的炮域偏移距道集转换为角度域共成像道集(ADCIGs)可为偏移速度分析(MVA)和叠前反演提供输入道集,并且ADCIGs是理论上没有假象的叠前反演道集,也是目前公认的精度最高的叠前反演道集。本文研究了基于矢量波场逆时偏移的弹性波保幅ADCIGs的提取方法,以保幅弹性波逆时偏移方程为基础,其核心是求取不同震源位置的纵、横波场在地下各成像点的入射角,对于转换波勘探,二者共享一个入射角,即震源纵波入射角。根据几何关系,震源纵波波场的传播角、构造的局部地层倾角之差为震源纵波入射角,震源纵波波场的传播角利用解耦后纵波场的极化向量得到,构造的局部地层倾角利用偏移叠加剖面的复波数得到。对纵、横波的共炮点偏移道集按入射角重新排列即可得到各自的ADCIGs。文中利用水平层状介质模型、倾斜层状介质模型、Marmousi—II弹性波部分模型和实测资料验证了算法的有效性,计算结果表明,本文方法计算的纵、横波角度具有较高的精度,提取的角道集具有较好的振幅保真性,能够为MVA和叠前反演提供可靠的输入道集。