简介:地震数据中的面波是严重降低地震资料信噪比的干扰波,它的存在影响了后续地震资料的处理与解释。本文根据地震记录中面波与反射波信号形态结构的差异,采用基于二维字典形态成分分析方法对面波噪声与反射波进行分离。根据面波信号的低频、低视速度和频散的特性,选择二维非抽样离散小波变换作为面波的稀疏表示字典,根据反射波局部相关性较强的特点,选择二维局部离散余弦变换作为反射波的稀疏表示字典,构建地震记录在联合二维字典下的稀疏表示模型并采用块协调松弛算法进行求解,将地震记录分解为反射波部分和面波部分。对合成地震信号以及实际地震资料的处理结果表明本文方法不仅能有效压制强能量的面波干扰,而且还能很好保护反射波信号的波形。
简介:太行山山前断裂位于太行山脉与华北平原过渡地带,是华北及我国东部地区一条重要的构造带,有关该断裂带的活动性有以下观点:1)它属深大断裂带;2)它是活动断裂带和地震构造带;3)它不是地震构造带。为研究和评价该断裂的活动性,利用最近获得的地震资料和钻探联合剖面研究了太行山山前断裂的深部构造背景和第四系以来的活动特征。钻探验证以及研究结果与表明:太行山山前断裂不是岩石圈深大断裂,该断裂在第四系仅错断了中更新世QP2早期,中更新世QP2晚期和晚更新世地层没有发生错断。由此得出,太行山山前断裂不是活动断裂,也不是岩石圈深大断裂。该研究结果为跨越太行山山前断裂的城市进行规划和建设提供了重要的地质和地球物理资料,对于首都圈地震危险性评价也具有十分重要的意义。
简介:传统的地层压力预测方法大都建立在地层欠压实机理之上,对于由其它因素引起的异常高压现象未必适用。Yassir等(1999)发现异常地层压力通常与异常水平应力有密切联系,异常高压带的最大主应力往往很大。基于此,在分析重力场和构造应力场的双重应力场特征基础上,本次研究首先通过常规测井资料建立构造稳定地区的最大主应力的计算模型,将利用该模型计算构造挤压强烈的非稳定地区的主应力结果定义为虚拟最大主应力,再结合相对构造应力的贡献值得到最大主应力,最后在一定的超压范围之内,根据最大主应力与地层压力的拟合关系预测地层压力。通过对A气田的实际资料处理表明,利用该方法得到的地层压力预测值与实际测量值吻合较好,预测精度较传统的等效深度法有明显的提高。
简介:本文提出了一种矿山微地震震源函数反演的新方法。以球坐标系波动方程为基础,对时空域中的微地震观测信号进行层析成像投影到慢度时间域,我们可以获得微地震波成像能量最大值及其对应的震源位置、发震时间和传播速度以及慢度时间信号。通过分析能量最大值处的慢度时间信号与震源函数之间的关系,推导出了计算矿山微地震震源函数的层析成像公式和利用震源函数计算微地震有效辐射能量的公式。进一步利用最小二乘法,将震源函数的振幅谱拟合成模型震源谱,确定了微地震波的零频极限值和拐角频率,最终计算出该微震事件的震源破裂半径。利用这一方法对理论模型数据和实际资料处理,结果表明:本方法具有既能求出任意一个地震事件的震源位置、发震时间和传播速度,又能同时求出震源函数及其频率特征值和震源破裂半径等震源参数的优点,因此,本方法是一种实时快速计算方法。
简介:传统的AVO正演只考虑了单。界面的反射系数对地震波波场振幅的影响,忽略了地震波在介质中传播的各种传播效应。通过引入地震波在介质中传播的几何扩散、吸收衰减以及透射损失等传播效应,提出了基于射线理论的水平层状介质多波保幅AVO正演方法。推导了水平层状介质多波几何扩散校正公式,来描述多波在介质中传播的几何扩散效应。通过直接引入复旅行时,而无需借助复速度,建立了复旅行时与品质因子的关系,来描述粘弹介质的吸收衰减。直接求解Zoeppritz方程计算多波的透射系数,用于描述多波在介质中传播时的透射损失。数值计算表明,几何扩散、吸收衰减以及透射损失对多波振幅的影响是随偏移距变化而变化的,多波保幅AVO正演需要考虑波传播效应对反射波振幅的改造。
简介:储层物性参数的改变伴随着弹性参数某种程度的改变,而这种改变并非简单的线性关系。加上观测信息的缺失、重叠、噪声破坏及模型理想化等原因,致使这种改变关系的获取具有较大的不确定性。针对传统统计岩石物理物性参数反演方法的不足,以利用弹性参数反演储层物性参数为目的,依据贝叶斯反演框架,我们建立了新的储层物性参数目标反演函数。首先,采用兼具确定性与随机性特点的统计岩石物理模型,考虑到不同弹性参数间的精度存在差异,引入权重系数,建立起储层物性参数与弹性参数间的加权统计关系。其次,基于这种加权统计关系,结合马尔科夫链蒙特卡洛随机模拟技术产生储层物性参数、弹性参数随机联合样本空间作为目标函数求解样本空间。最后,建立解的快速求解准则,求取最大后验概率密度对应的储层物性参数取值作为最终解。实际应用表明,该方法具有较高的反演效率,应用前景较好。
简介:本文提出的储层物性参数同步反演是一种高分辨率的非线性反演方法,该方法综合利用岩石物理和地质统计先验信息,在贝叶斯理论框架下,首先通过变差结构分析得到合理的变差函数,进而利用快速傅里叶滑动平均模拟算法(FastFourierTransformMovingAverage,FFT-MA)和逐渐变形算法(GradualDeformationMethod,GDM)得到基于地质统计学的储层物性参数先验信息,然后根据统计岩石物理模型建立弹性参数与储层物性参数之间的关系,构建似然函数,最终利用Metropolis算法实现后验概率密度的抽样,得到物性参数反演结果。并将此方法处理了中国陆上探区的一块实际资料,本方法的反演结果具有较高的分辨率,与测井数据吻合度较高;由于可以直接反演储层物性参数,避免了误差的累积,大大减少了不确定性的传递,且计算效率较高。
简介:常规多波联合反演采用Zoeppritz方程的近似式构建正演方程,反演过程中需要假定背景纵横波速度比为常数,其反演精度不高,稳定性不好。本文提出了一种基于精确Zoeppritz方程的多波联合反演方法,结合贝叶斯方法进行广义线性反演。本方法基于精确Zoeppritz方程构建正演方程,避免了近似式反演在大角度时引起的误差;利用贝叶斯方法引入模型参数的先验分布信息,作为反演的正则化项,降低了反演的不适定性;反演目标函数中引入低频软约束,稳定了反演低频结果,提高了反演的鲁棒性;在求解反演目标函数时,利用快速算法,降低了反演的运算量。经过模型试算,证明了该方法的优越性和抗噪性;并在实际资料的应用中证明了该方法的实用性和有效性。
简介:本文以中观孔隙结构的White模型为基础,构建了部分饱和孔隙介质模型,利用Biot方程的建立思路和Johnson推导的体变模量,推导了部分饱和孔隙介质中的纵波方程,并以平面波为例,求取了方程的衰减系数,分析了地震勘探频带范围内地震波的衰减特性。结果表明:在部分饱和孔隙介质中,地震波在低频段也会发生明显的衰减和频散现象,频率越大,衰减越大;且第二纵波的衰减比第一纵波更为明显;这一结论弥补了Biot理论在描述地震勘探频带范围内波的衰减现象的不足。文中还研究了孔隙度、饱和度和模型内径尺寸对第纵波衰减特性的影响机理,主要表现在在地震勘探频带范围内,波的衰减随孔隙度的增大而增大,随含油气饱和度的增大而减小,当孔隙内径尺寸小于二分之一外径尺寸时,波的衰减随内径尺寸的增大而增大,当内径尺寸大于二分之一外径尺寸时,波的衰减随内径尺寸增大而减小。
简介:频率-波数域单程波算子能高效地模拟地震波在复杂介质中的传播,但是在描述波的大角度传播和速度横向扰动变化较大介质中传播的问题时仍然存在一定误差。这类误差是由于对单平方根算子使用Taylor展开式的近似程度不足所造成。为了进一步提高泰勒展开式的精确性,本文提出一种利用粒子群智能算法优化级数展开系数的高阶广义屏算子对单平方根算子的展开级数进行优化处理。新的偏移算法能在保持单程波偏移算法高效的前提下进一步提高偏移算子在大角度的成像精度和对强横向速度变化介质的适应性。通过脉冲响应实验,验证了基于粒子群算法优化级数的高阶广义屏算子能够提高常规的高阶广义屏算子的成像精度和成像角度。根据对二维SEG/EAGE盐丘模型的成像处理,基于粒子群算法优化级数的高阶广义屏算子对盐丘下面的断层取得了更高质量的成像,说明粒子群优化级数的高阶广义屏算子比常规的高阶广义屏算子具有更好的横向速度适应性。为了检验本文所提算法对实际资料的处理能力,我们利用常规的偏移处理技术和本文所提算法对一条海上二维数据进行了偏移成像处理,对比分析成像剖面发现本文所提算法描述了更加清晰的层位信息和更高质量的偏移剖面。本文所提算法能有效提高高阶广义屏偏移在广角度成像的能力,具有一定实际应用价值。
简介:噪声衰减是探地雷达信号处理中的关键问题之一。当探测目标埋藏深度比较浅时,其反射信号与直耦信号和地面回波信号相互重叠,直接影响目标反射波到达时刻的检测及目标的正确定位。针对这个问题,本文提出了一种基于Curvelet变换的噪声衰减方法。通过对理论数值模拟数据和实测数据的处理,以及与平均消去法和二维连续小波该方法处理结果的对比,验证了该方法的可行性和有效性。处理结果显示,该方法不仅可以去除背景噪声、同时可以衰减倾斜相关的相干干扰和数据中的随机噪声。与二维连续小波变换方法相比有更高的计算效率。