共查询到20条相似文献，搜索用时 15 毫秒
1.
地球介质相对于地震波波长尺度的定向非均匀性会导致波速的各向异性,进而影响地震波场的运动学与动力学特征．各向异性弹性波动方程是描述该类介质波场传播的基本工具,在正演模拟、偏移成像与参数反演中起着关键作用．为了面向实际应用构建灵活、简便的各向异性波场传播算子,人们一直在寻求简化的各向异性波动方程．本文借鉴各向异性弹性波波型分离思想,通过对平面波形式的弹性波方程(即Christoffel方程)实施一种代表向波矢量方向投影的相似变换,推导出了一种适应任意各向异性介质、运动学上与原始弹性波方程完全等价,在动力学上突出qP波的新方程,即qP波伪纯模式波动方程．文中以横向各向同性(TI)介质为例,给出了相应的qP波伪纯模式波动方程及其声学与各向同性近似,并在此基础上开展了正演模拟和逆时偏移试验,展示了这种描述各向异性波场传播的新方程的特点与优势． 相似文献
2.
三、基于波动方程反演的散射层析成象方法地震波的传播是服从弹性波方程的,但是要严格的从弹性波方程反演来研究散射层析成象是远远超出当前数学发展水平的,其初步探讨可见[13]。熟知弹性介质中存在两种波:压缩波(纵波)与剪切波(横波),如介质参数的变化率较小,则压缩波场u近似地满足变系数的声波方程如下 相似文献
3.
声波方程逆散射反演的近似方法 总被引：7，自引：0，他引：7
我们在文献[1]里研究了介质参考波速沿某个方向线性变化时的三维声散射理论,导出了声波方程逆散射反演问题解的解析表达式.考虑到应用时的实际条件,本文根据上述反演方法导出2.5维模型的声波方程逆散射反演的波速扰动计算公式,给出该方法在“高频”近似条件下的波速扰动反演计算公式,从而使我们提出的“参考波速线性变化时的声波方程逆散射反演”理论更接近实际应用条件.本文给出的这些反演公式仍然具有原方法的优点,即不但可以使Born近似的假定在大多数情况下能得以满足,而且可以利用快速Fourier变换来快速实现介质波速扰动的反演成象. 相似文献
4.
我们在文献[1]里研究了介质参考波速沿某个方向线性变化时的三维声散射理论,导出了声波方程逆散射反演问题解的解析表达式.考虑到应用时的实际条件,本文根据上述反演方法导出2.5维模型的声波方程逆散射反演的波速扰动计算公式,给出该方法在“高频”近似条件下的波速扰动反演计算公式,从而使我们提出的“参考波速线性变化时的声波方程逆散射反演”理论更接近实际应用条件.本文给出的这些反演公式仍然具有原方法的优点,即不但可以使Born近似的假定在大多数情况下能得以满足,而且可以利用快速Fourier变换来快速实现介质波速扰动的反演成象. 相似文献
5.
土是由一定尺寸大小颗粒所构成的多孔介质,具有明显的颗粒特性,当土颗粒间的孔隙被流体（如水或油）充满时则成为饱和土.利用微极理论和Biot波动理论的研究成果,把饱和土中多孔固体骨架部分近似地视为微极介质,孔隙中的流体部分视为质点介质,获得饱和多孔微极介质的弹性波动方程.借鉴Greetsma理论,建立了饱和多孔微极介质弹性本构方程力学参数与相应单相介质弹性参数的相互关系,使饱和多孔微极介质弹性波动方程中的物理参数具有明确的物理意义,易于在试验中确定.运用场论理论把饱和多孔微极介质的波动方程简化为势函数方程,建立了饱和多孔微极介质中五种弹性波的弥散方程,数值分析了五种简谐体波在无限饱和多孔微极介质中的传播特性. 结果表明,P1波、P2波和剪切S1波的波速弥散曲线与经典饱和多孔介质基本相同,当频率小于临界频率ω0时旋转纵波θ波和横波S2波不存在,当频率大于临界频率ω0时,θ波和S2波的传播速度随频率增加而减小. 相似文献
6.
声波方程的逆散射反演乃是求解双曲型偏微分方程系数项反问题的一种解析方法,一般利用Born近似把这一非线性反问题线性化,并给出了恒参考波速介质中反问题解的解析表达式.由于Born近似假定波速扰动为一级无穷小,因此,在大多数情况下,恒参考波速介质模型的反问题的解无法得以应用.本文研究介质参考波速沿某个方向线性变化时的声散射理论,导出了声波方程逆散射问题解的解析表达式,从而既可使Born近似的假定在大多数情况下能得以满足,又可利用快速Fourier变换快速实现介质波速扰动的反演成象. 相似文献
7.
当地震波在地层中传播时,能量的衰减是传播介质非弹性性质的反映.这种介质所固有的衰减特性通常用品质因子(Quality Factor, Q)来描述.本文基于地震波在传播过程中主频降低的衰减特性,利用质心频率偏移法(Centroid Frequency Shift,简称 CFS)对零偏VSP正演模型进行Q值计算,分析薄层、上下行波场、频带宽度、反演波速、震源位置等因素的影响.模型结果显示:CFS法比频谱比法、振幅衰减法能更加准确识别薄层界面,可准确反演出厚度10m的地层;上行波场的加入降低反演结果准确度,尤其对高Q值层;当有效波频带增宽时,高频部分衰减明显,Q值反演结果接近理论值;反演波速误差对反演准确性影响很大.最后,利用CFS对实际井旁地面地震资料进行Q值反演,反Q滤波后其主频由25Hz提高到35Hz,分辨率得到改善,证明了CFS方法的实用性. 相似文献
8.
借助Christoffel方程可求解出各向异性介质弹性波精确频散关系.利用近似方法进行处理,再通过傅里叶逆变换将频率波数域算子变换为时空域算子,可导出解耦的 qP波或 qS波波动方程.本文在 TTI介质弹性波精确频散关系的基础上,利用近似配方法推导了 qP波和 qSV波近似频散关系,通过傅里叶逆变换推导了 TTI介质 qP波和 qSV波解耦的波动方程.为了验证近似频散关系的有效性,利用两组模型参数对其进行数值计算,分析了相对误差在不同传播方向上的分布.随后使用有限差分方法分别对均匀、层状及复杂 TTI介质弹性波近似解耦波动方程进行数值模拟,结果显示 qP波和 qSV波完全解耦,并且在各向异性参数η＜0 以及介质对称轴倾角变化较大的情况下,纯 qP波和纯 qSV波近似波动方程依然可以保持稳定. 相似文献
9.
当前全波形反演方法研究大多针对弹性介质,忽略了真实介质的吸收衰减作用。现有针对粘弹介质的波形反演方法也主要在Q模型给定的假设前提下聚焦于速度参数反演,其Q模型的获取通常利用层析反演得到,分辨率低,不能精确匹配传播过程中的振幅衰减和相位畸变作用,一定程度上影响反演精度。本文提出了一种波动方程Q值波形反演新方法。从描述衰减机制的标准线性体理论出发,本文首先推导得到明确表征振幅衰减和相位畸变作用的简化粘滞声波方程。相比传统方程,本文中方程形式不涉及记忆变量,计算过程中内存需求少,且衰减补偿过程更易于实现。本文进一步得到该方程对应的伴随方程形式和目标函数关于模型参数的梯度表达式,并针对性地给出了克服伴随方程传播不稳定的规则化方案。在速度参数已知的条件下,利用L-BFGS方法实现了Q值波形反演。为缓解波形反演对初始模型的依赖性,一定程度上避免"周波跳跃"问题,本文采用多尺度分析策略进行复杂模型Q值波形反演。此外本文还进一步讨论了目标函数的抗噪性及速度和品质因子双参数波形反演问题,并结合模型试算给出结论与认识。 相似文献
10.
目前转换波AVO(Amplitude Versus Offset)反演或纵波、转换波联合AVO反演是基于Aki等的反射系数近似公式,通常反演结果为上、下层速度比值和速度平均比值.本文在Shuey和郑晓东近似公式的基础上,运用最小平方算法,提出一种p-p波和P-SV波联合AVO反演下层介质纵、横波速度及泊松比的方案,并应用于双层含气砂岩模型上,误差在10%以内.此方案要求上、下层物性差别不大且入射角小于30°,但不需要速度比值为常数2的假设. 相似文献
11.
地下岩石的速度各向异性影响地震波的传播与成像.横向各向同性(TI)介质为最普遍的等效各向异性模型.引入TI介质拟声波方程可以避免复杂的弹性波方程求解以及各向异性介质波场分离,以满足对纵波成像的实际需要.本文从垂直横向各向同性(VTI)介质弹性波方程出发,推导出正应力表达的拟声波方程以及相应的纵波分量的表达式,进而分析从频散关系得到的拟声波方程的物理意义,而后将拟声波方程扩展到更一般的倾斜横向各向同性(TTI)介质中.波前快照与群速度平面的对比验证了拟声波方程可以很好地近似描述qP波的运动学特征.在此基础上,将拟声波方程应用在逆时偏移中并与其特例声波近似方程进行对比,讨论了计算效率、稳定性等实际问题.数值试验表明VTI介质情况下采用声波近似方程可以提高计算效率,而TTI介质qP-qSV波方程则在效率相当的情况下可以保证稳定性.SEG/HESS模型和逆冲模型逆时偏移试验验证了本文TI介质拟声波方程的实用性. 相似文献
12.
常见弹性波动理论的建立是基于介质均匀这一基本假设,实际介质的非均匀性非常普遍.为研究连续介质中波的传播特征,本文从弹性力学中建立弹性波动方程的三个基本方程出发,考虑连续介质弹性参数的空变特征,建立非均匀介质的弹性波动方程,利用Alkhalifah声学近似思想建立位移表征的纵波波动方程,利用本征值问题求解方法建立标量波频率-波数域传播算子,从而建立描述纵波传播的标量波方程,其中波函数为纵波位移的散度,不同于均匀介质标量波方程的波函数为位移势.随后推导含PML边界波动方程差分格式并建立不同模型数值模拟进行数值试算,与均匀假设标量波方程和变密度方程对比证明本方法的准确性和稳定性. 相似文献
13.
频率域波动方程的双参数识别方法 总被引：1，自引：0，他引：1
利用反射波地震记录恢复波动方程的两个系数以研究地球内部构造，是当前应用地球物理研究的热点之一，本文在波动方程反演方法的理论基础上，导出了频率域识别介质的密度波速两个重要参数的新算法，从而把当前的波动方程单参数反演双参数反演方法，通过理论模型及实际测井资料的反演试算对比，均藜得良好效果，在算法中还省去了Ｇｒｅｅｎ函数的繁杂计算，节省了计算量，增加实用性。 相似文献
14.
随着地震勘探和开发的不断深入,面向地质目标的精细储层预测技术变得越来越重要.由于透射损失、层间多次波、波模式转换以及随机噪声等的影响,观测地震数据和待反演的地下介质属性之间呈现出很强的非线性.考虑到这些非线性,本文基于积分波动方程开展叠前地震反演,从观测地震数据中恢复出介质属性和整体波场,其中反演参数是波动方程中的压缩系数、剪切柔度和密度的对比度,相比于常规线性AVO反演的波阻抗弹性参数,它们对流体指示有更强的敏感性.在反演过程中,从平滑的低频背景场出发,交替迭代求解数据方程和目标方程.采用乘性正则化方法于共轭梯度框架下求解反演参数,采用优化的散射级数Neumann序列获得整体波场,这种方法不易陷入局部极值,能收敛到正确解.测井资料和典型山前带模型测试表明,利用上述反演方法能获得高分辨率的深度域地下介质属性,可直接进行储层预测和解释. 相似文献
15.
基于地震波反射系数近似公式的叠前反演是油气勘探的重要工具.本文在已有研究的基础上,推导了一个改进的射线参数域地震纵波反射系数近似方程.该方程建立了地震纵波反射系数与纵波阻抗和横波阻抗的非线性关系,在中、小角度的范围内较现有的反射系数线性近似公式精度更高.另外,由于该方程仅包含纵波和横波阻抗反射系数项,因此基于新方程的反演能够有效地降低同步反演纵波速度、横波速度、密度三个参数的不适定性.在此基础上,结合广义线性反演法(GLI)理论和贝叶斯理论,相应地发展了一种叠前地震同步反演方法.模型测试和实际资料的应用表明,基于新方程的反演方法能够利用有限角度(偏移距)的数据稳定地反演纵波和横波阻抗,由于在反演过程中,不需要假设纵横波速度为常数,因此该方法还能有效地提高反演结果的精度. 相似文献
16.
基于Zoeppritz方程对介质密度偏导数所建立的偏导方程的精确解,构造了多角度反演地层介质密度的反演方程,在偏导数求解过程中考虑了介质密度对波速度的影响因素,并由此实现了利用反射系数梯度精确解计算地层密度的多角度联合反演.通过数值算例考察了计算方法,结果显示:反演方法对层状地层模型不论反射波是否存在相干现象均获得了较好的反演结果,反演迭代10次后计算结果的最大相对误差能够收敛到1%之内;随着反演角度的增加地层介质密度反演的精度逐步提高,反演具有自动校正能力,有快的计算速度.本方法克服了传统AVO(Amplitude Versus Offset)基于Zoeppritz方程近似所遇到的困难,不受反演角度大小及反射界面对波反射强弱的限制,为地层介质密度的多角度包括大角度反演提供了一种新的快速有效的计算方法. 相似文献
17.
Multi-angle inversion of formation densities based on the accurate solutions of Zoeppritz's partial derivative equations 
下载免费PDF全文
下载免费PDF全文 基于Zoeppritz方程对介质密度偏导数所建立的偏导方程的精确解,构造了多角度反演地层介质密度的反演方程,在偏导数求解过程中考虑了介质密度对波速度的影响因素,并由此实现了利用反射系数梯度精确解计算地层密度的多角度联合反演.通过数值算例考察了计算方法,结果显示:反演方法对层状地层模型不论反射波是否存在相干现象均获得了较好的反演结果,反演迭代10次后计算结果的最大相对误差能够收敛到1％之内;随着反演角度的增加地层介质密度反演的精度逐步提高,反演具有自动校正能力,有快的计算速度.本方法克服了传统AVO(Amplitude Versus Offset)基于Zoeppritz方程近似所遇到的困难,不受反演角度大小及反射界面对波反射强弱的限制,为地层介质密度的多角度包括大角度反演提供了一种新的快速有效的计算方法. 相似文献
18.
传统的伪谱(PS)方法,采用傅里叶变换(FT)计算空间导数具有很高的精度,每个波长仅需要两个采样点,而时间导数采用有限差分(FD)近似因而精度较低.当采用大时间步长时,由于时空精度不平衡,PS法存在不稳定性问题.原始的k-space方法可以有效地克服这些问题但是却无法适用于非均匀介质.为了提高原始k-space方法模拟非均匀介质波动方程的精度,我们提出了一种新的k-space算子族.它是用非均匀介质的变速度代替原k-space算子中的常数补偿速度构造得到,引入低秩近似可以高效求解.我们将构造的新的k-space算子应用于耦合的二阶位移波动方程,而不是交错网格一阶速度应力波动方程,使模拟弹性波的计算存储量减少.我们从数学上证明了基于二阶波动方程的k-space方法与基于一阶波动方程的k-space方法是等价的.数值模拟实验表明,与传统的PS、交错网格PS和原始的k-space方法相比,我们的新方法可以在时间和空间步长较大的均匀和非均匀介质中,为弹性波的传播提供更精确的数值解.在保持稳定性和精度的同时,采用较大的时空采样间隔,可以大大降低数值模拟的计算成本. 相似文献
19.
正浅层地表介质的剪切波速是工程抗震设计、场地地震危险性分析、场地地震效应研究等研究的重要基础。面波勘探方法是一种利用人工震源激发产生的瑞雷面波来反演地下介质弹性参数特别是剪切波速的现场原位测试方法,工作效率高且经济成本低。本文针对面波勘探中相对较新的多道面波分析方法(MASW)在地下介质剪切波速探测中的应用进行了研究。主要研究内容及取得的成果如下:(1)利用传递矩阵法(Dunkin)计算了水平层位条件下的瑞雷波相速度频散曲线,并总结 相似文献
20.
本文以波动理论为基础, 半解析化求解地震勘探中常用的SH波方程. 获得的主要结果包括: 给出了二维均匀介质中SH波方程的解析解; 利用Cagniard-de Hoop方法详细推导了二维双层介质中SH波方程的解析解, 获得了透射波的解析解表达式. 同时, 基于SH波方程的解析表达式, 给出了包含各种波(如直达波、反射波、首波以及透射波)的解析解和波形图. 对于比较复杂的积分型解析解, 利用数值积分方法给出了数值结果, 并与优化的近似解析离散化方法(ONADM)和4阶Lax-Wendroff修正方法(LWC)的数值结果进行了比较, 以验证解析解的正确性. 本文的研究成果有望在检验波动方程数值新方法的有效性、波传播理论分析等方面得到应用. 相似文献