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

  • 标签: 波动方程 延拓 地形起伏 变换
  • 简介:重构反演是一种改进的全波形反演理论。该反演方法通过将波动方程引入目标函数中拓宽了解的寻找空间,通过重构真实来计算模型梯度,大大提高了计算效率的同时还减弱了局部极小值的影响。但目前该理论基本在频率域进行,而频率域反演对计算内存的需求太高,并且很难应用到实际生产中。因此,本文将重构反演拓展到时间域,推导了时间域重构的增广方程,结合模型试算结果对重构的模型梯度进行了修改。数值实验表明,时间域重构反演准确性较高并且对低频信息具有良好的重建能力。

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

  • 标签: 弹性波 逆时偏移 矢量波场 角道集 保幅偏移
  • 简介:便于不同观测系统的统一,本文定义了广义炮检距概念,给出了空间平界面广义炮检距不同阶地震菲涅耳带表达式。基于波动理论,推导出了广义炮检距地震菲涅耳带横向叠加公式。以不同阶几何菲涅耳带形状为参考,分为零炮检距和非零炮检距情况,进行了相似菲涅耳带不同绕射面元大小的横向叠加振幅分析。结果表明:①绕射面元对观测点的贡献与炮检距、地表起伏程度、界面倾斜程度、激发点到界面深度、观测方式和干涉叠加区大小等因素有关;②第一菲涅耳带是主要绕射干涉叠加区,并且该区域对观测点的半振幅贡献约小于所有阶菲涅耳带的振幅贡献;③当绕射面元小于第一菲涅耳带时,即使采用非自激自收观测方式,仍有绕射面元越大观测点振幅越大的结论。

  • 标签: 炮检距 菲涅耳带 波动理论 叠加 振幅
  • 简介:针对有限差分数值模拟的频散问题,本文将交错网格技术和紧致差分格式相结合,推导了横向各向同性介质一阶速度一应力波动方程的紧致交错网格差分格式;对比分析了紧致交错网格差分格式、交错网格差分格式以及紧致差分格式的截断误差主项,并利用Fourier误差分析方法分析了上述三种差分格式的近似精度;在此基础上,分别采用上述三种差分格式进行了数值模拟。结果表明,当差分方程阶数相同时,紧致交错网格差分格式截断误差最小,数值频散最弱,差分精度最高,证实了该方法的有效性。

  • 标签: 横向各向同性介质 紧致交错网格 一阶速度-应力波动方程 数值频散 波场模拟
  • 简介:应用多分量地震资料进行成像时通常需要先做分离,然后再对分离的型进行成像。其中,分离可以在空间域或波数域实现。然而,由于用交错网格有限差分进行弹性数值模拟时,用来进行波数域分离的质点振动速度分量定义在不同网格节点上,本文提出了利用波数域插值方法来估算同一网格节点所需质点振动速度值;进而给出了先进行波数域插值后进行波分离的波数域保幅分离方案。数值实验结果表明波数域插值方法具有较高的插值精度且保幅分离方法具有较好的保幅性,将本文方法进一步应用于弹性逆时偏移可以获得保幅性较好的成像结果且对存在一定程度速度误差情况具有较好的适应性。

  • 标签: 波场分离 保幅性 交错网格有限差分 波数域插值 逆时偏移
  • 简介:基于Aki-Richards公式和贝叶斯原理,本文发展了利用叠前PP和PS资料联合反演P波速度比、S波速度比和密度比的方法。该方法假设参数之间满足正态分布,引入参数协方差矩阵来描述反演参数之间的相关性以提高反演过程的稳定性,并同时使反演的参数序列服从Cauchy分布,引入矩阵Q来描述参数序列的稀疏性以提高反演结果的分辨率。采用本文提出的方法对模型数据和实际多资料进行反演,结果表明:本文方法正确有效;与传统的单一PP反演相比,PP和PSAVO联合反演具有稳定性更好和反演精度更高等优点。

  • 标签: PP波 PS波 联合反演 贝叶斯原理
  • 简介:在文献中,Kirchhoff型反偏移的稳相分析主要是在下列两个条件下进行的:(1)等时面和目标反射面相切;(2)深度偏移像信号的长度接近于零。对于与目标反射面不相切的等时面和长度远大于零的深度偏移像场子波,已有的结果将不再成立。为了在等时面和目标反射面不相切和深度偏移像的子波长度远大于零的条件(一般条件)下对Kirchhoff型反偏移进行稳相分析,我推导了出现在二维稳相分析公式中的诸因子的计算公式,并从中发现:(1)对于不同的等时面,距离差函数的稳相点具有不同的水平坐标;(2)Kirchhoff型真振幅反偏移的输出场由两部分(真振幅反偏移信号与振幅畸变因子)的乘积组成。由此得到下列两个结论:(1)一个给定的反偏移信号由多个深度偏移信号上的采样点组装而成,反偏移信号上的采样点个数等于对于这种组装有贡献的偏移信号的个数。(2)振幅畸变效应是Kirchhoff型反偏移中的固有效应,靠反偏移本身无法消除。如果一定要消除这种振幅畸变效应,必须对反偏移结果进行振幅校正。

  • 标签: Kirchhoff型反偏移 稳相分析 振幅畸变
  • 简介:常规地震勘探大多基于反射理论,对尺度小、倾角陡、构造复杂的灾害性地质体,很难得到精确的成像结果。本文建立2类典型矿井地质模型,采用有限差分正演模拟算法,研究了矿井地震散射的特征,将等效偏移距偏移(EOM)和基于等效偏移距干涉叠加偏移的散射成像方法用于矿井地震勘探中,重点分析了散射成像相比常规反射地震成像的技术优势。研究表明:1)散射成像方法能够提高覆盖次数,充分利用有效的弱散射信号,对不均匀复杂地质体的散射的成像效果明显优于反射叠后偏移;2)对于巷道超前探测,它弥补了基于射线理论的叠加偏移的不足,为提高矿井地震勘探分辨率提供了有效的成像方法。

  • 标签: 矿井 地震勘探 散射波 地震成像 数值模拟
  • 简介:裂缝识别对碳酸盐岩储层测井评价至关重要,传统的测井仪器由于探测深度浅(小于3米)而无法对井壁外围(大于3米)裂缝发育情况进行评价,远探测声波测井仪器采用相控阵发射、同时加大源距,有效提升了测井仪器的探测深度。但由于缺少正演模拟研究,目前对于裂缝的解释往往是基于经验而缺乏理论依据,导致很多现象难以解释。本文利用高阶有限差分方法对远探测声反射测井裂缝识别进行了正演模拟及叠前逆时偏移成像研究,首先在理论研究的基础上构建正演模型,重点研究不同裂缝离井壁距离、裂缝张开度和倾角的响应特征;其次在单因素变化基础上提取成像区域的能量强度,分析确定出在实际地层速度有波动变化时远探测声反射测井方法能够识别裂缝的条件;最后通过对裂缝识别的影响因素定量化分析,确定裂缝识别的最大距离、最小张开度和最小倾角,降低了裂缝识别与评价中的多解性。研究成果对远探测声反射测井仪器的发展、数据处理方法的改进,以及后续的测井解释工作都有一定的借鉴意义。

  • 标签: 远探测声波 逆时偏移 高阶有限差分 波动方程 反射波成像
  • 简介:传统叠前深度偏移只能够提供地下的构造信息,但工业界在需要构造信息的同时还要与地下界面反射系数成比例的振幅信息。最近几年,基于单程方程的保幅叠前深度偏移算法有了一定的发展,但是,基于炮域、单程的保幅型叠前深度偏移必须应用反褶积型的成像条件,这种成像条件在构造复杂、速度变化剧烈的地区会出现不稳定现象。基于角度域的保幅深度偏移克服了这一不稳定性缺点的同时,还域的保幅深度偏移,模型和实际资料的试算分析验证该思路方法的正确性和有效性。

  • 标签: 保幅偏移 叠前深度偏移 角度域 单程波方程 成像条件
  • 简介:震源定位是微地震监测关键技术之一。本文提出用于微地震定位的弹性和加权弹性(WEW)干涉成像方法。该方法在保证定位精度的同时,还可避免震源假象。通过各向同性水平层介质状模型的数值试验,初步表明该方法可适应低信噪比微震信号、速度随机扰动、较稀疏的检波器分布等情况,并在速度模型存在一定的系统误差时也仍保持较高的定位精度。由于干涉成像方法不需要进行初至拾取,定位效率相对传统走时方法也得到了提高。采用二维断层模型试算Nnumericalresultsofusingatwo-dimensionalfaultmodel,表明方法还能实现多震源定位,且比逆时成像有更高的定位精度。

  • 标签: 微地震监测 震源定位 弹性波 干涉偏移 逆时成像
  • 简介:为了探明介质的非完全弹性对瑞利型槽传播的影响以及利用瑞利型槽品质因子QR反演煤层碳化程度、裂隙、煤层厚度等煤层属性,在本文中,我们将复速度引入瑞利型槽频散方程,计算了煤层瑞利型槽的品质因子QR曲线;采用控制变量法,分析了瑞利型槽品质因子QR随煤层厚度、煤层纵横波品质因子以及围岩纵横波品质因子的变化规律。研究表明:瑞利型槽品质因子QR曲线与群速度曲线变化趋势一致;埃里相(Airy-phase)频率附近瑞利型槽的品质因子QR最小,且埃里相频率随煤层厚度增加而减小;瑞利型槽品质因子QR随煤层横波品质因子QS2的增大而增大。利用瑞利型槽品质因子QR曲线可以实现对瑞利型槽的吸收衰减进行补偿、对煤层厚度进行预测以及岩性成像。更多还原

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

  • 标签: 地震波场模拟 完全匹配层PML 衰减因子
  • 简介:电磁全息测量数据包含两种模态,因此重建流动图像也需采用“双模”融合的敏感。首先,通过电磁全息探测物理分析,结合层析成像数学理论Radon反变换,从定解问题推导了全息测量敏感函数表达式。其次,将有限元仿真计算得到的全息测量敏感应用于模拟流动试验和全息成像。结果表明,基于复电位φ关于极径r的偏导数的全息测量敏感契合了Radon反变换的数学表达,且充分体现了幅度、相位测量敏感性;反演所得流动图像更加吻合实际流型。全息测量敏感构建有效克服了传统单模敏感在计算精度和计算效率等方面的局限性。

  • 标签: 全息测量敏感场 全息探测物理场 有限元 全息成像
  • 简介:地震波数值模拟方法对理解和分析地震波的传播规律具有着重要的意义。弹性波动方程能够模拟地下介质的实际情况,为偏移和成像提供有效的依据。在弹性波波数值模拟中,旋转交错网格数值模拟(RSM)修改了标准交错网格数值模拟(SSM)方法,将同类的参数定义在同样的节点上,拓宽了稳定性条件的约束,但在低速区会出现较严重的频散。变阶数差分方法是自适应空间算子长度方法的一种变化和推广。它以理论频散误差研究为基础,结合实际传播的情况进行误差计算,对不同速度匹配不同的差分阶数。本文研究了变阶数旋转交错网格数值模拟(VRSM),即是籽变阶数方法应用到RSM中,它可以很好地解决RSM在低速区域的数值频散问题,以及减少不必要的时间损耗;同时讨论了旋转交错网格的理论频散特性,并基于分离的方法分析了实际传播的频散误差,将原方法的应用范围由声波推广到剪切,由理论值推广到时变值。在数值模拟试验中,VRSM将被应用于水平层模型和Overthrust模型。通过阶数分配以及相应传播效果和计算时间的分析,验证了该方法应用于复杂介质模拟中的实用性和有效性。实验的结果表明VRSM能够合理分配不同速度所对应的差分阶数,能保证计算的精确性,并合理控制计算的时间。

  • 标签: 弹性波 变阶数 旋转交错网格 横波 时变
  • 简介:传统的AVO正演只考虑了单。界面的反射系数对地震波波振幅的影响,忽略了地震波在介质中传播的各种传播效应。通过引入地震波在介质中传播的几何扩散、吸收衰减以及透射损失等传播效应,提出了基于射线理论的水平层状介质多保幅AVO正演方法。推导了水平层状介质多几何扩散校正公式,来描述多在介质中传播的几何扩散效应。通过直接引入复旅行时,而无需借助复速度,建立了复旅行时与品质因子的关系,来描述粘弹介质的吸收衰减。直接求解Zoeppritz方程计算多的透射系数,用于描述多在介质中传播时的透射损失。数值计算表明,几何扩散、吸收衰减以及透射损失对多振幅的影响是随偏移距变化而变化的,多保幅AVO正演需要考虑传播效应对反射振幅的改造。

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

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

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

  • 标签: 转换波 AVO 反射系数 联合反演 多分量