简介:频率域全波形反演充分利用全波场的振幅、相位以及频率信息,采用较少的频率便能反演得到精度很高的速度模型。本文以有限单元法为基础,对起伏地形条件下二维声波频率域全波形反演进行了研究。在正演算法中,针对截断边界问题,并考虑多频率联合反演中计算区域采用同一套剖分网格的需求,提出了一种适用于起伏地形的衰减边界条件算法。该算法的核心思想是在控制方程波数项中引入衰减因子,通过一定方式调节衰减因子使得声波在衰减层中充分衰减,达到压制截断边界影响的目的。根据指数衰减规律,文中推导出了一种新的衰减因子计算公式,并给出了不同频率条件下衰减层厚度计算公式;在反演算法中,采用共轭梯度法求解高斯牛顿反演迭代方程组,避免直接求解雅克比矩阵和Hessian矩阵带来的巨额计算量,并采用相同的反演模型,对比分析了不同初始模型和频率组合对全波形反演结果的影响。起伏地形模型数值模拟和全波形反演数值试验表明,本文提出的指数衰减边界条件算法和基于该算法的全波形反演算法具有很好的应用效果。
简介:CO2地质封存四维地震监测的主要目标是监测地质封存的安全性和CO2驱油效果.注入CO2过程中,储层中流体饱和度和储层压力变化导致储层弹性参数及地震响应发生相应变化.从理论上讲,储层开发的基本地质条件是不变的,所以两次四维地震属性相减得到的成像显示,消去了油气藏静态性质(如构造、岩性等),得到的是油气藏的流体动态成像(如压力、饱和度等).本文采用的加拿大Weyburn油田CO2地质封存项目四维地震数据进行分析和研究,两次三维勘探数据采集时间分别是1999年(Baseline)和2002年(Monitor),CO2封存量为280×104t.首先,以测井资料为基础,通过Gassmann方程和人工合成地震记录,计算替换前后其速度和振幅等属性的变化,以此约束实际地震属性的筛选;然后,在对两次观测的四维地震数据进行四维地震匹配处理的基础上,进行两次观测地震数据的地震属性筛选与对比分析,利用两次地震属性的差异,识别注入CO2后储层内流体的分布特征;最后,以盖层重复性为指标,标定储层,提取储层的不同属性,同时结合人工合成地震记录的结果,得到一个能更好表现其流体变化的属性.
简介:标量CSAMT只适合一维及测量方向与构造方向垂直的二维情况,对于复杂的三维地电结构,CSAMT需采用张量测量。本文试图采用矢量有限元法实现三维张量CSAMT的正演模拟。为了验证算法的正确性,本文在层状介质中计算了三维CSAMT远区的电场,磁场及阻抗张量,并且与层状介质中的理论解进行了比较,接着还模拟了均匀半空间中含有三维异常体的模型,并且分析了四个阻抗张量、视电阻率及阻抗相位的响应特征。得出如下结论:采用矢量有限元法来模拟三维张量CSAMT,其电磁场及阻抗张量的实虚部计算精度都比较高,并且该方法本身满足电场法向不连续,不用进行散度校正。
简介:标量CSAMT只适合一维及测量方向与构造方向垂直的二维情况,对于复杂的三维地电结构,CSAMT需采用张量测量。本文试图采用矢量有限元法实现三维张量CSAMT的正演模拟。为了验证算法的正确性,本文在层状介质中计算了三维CSAMT远区的电场,磁场及阻抗张量,并且与层状介质中的理论解进行了比较,接着还模拟了均匀半空间中含有三维异常体的模型,并且分析了四个阻抗张量、视电阻率及阻抗相位的响应特征。得出如下结论:采用矢量有限元法来模拟三维张量CSAMT,其电磁场及阻抗张量的实虚部计算精度都比较高,并且该方法本身满足电场法向不连续,不用进行散度校正。
简介:无网格法形函数构造不依赖预定义的单元,具有计算精度高、处理复杂模型便利等优点。本文介绍了无单元Galerkin法(EFGM)、点插值法(PIM)与径向基点插值法(RPIM)三种全域弱式无网格法的近似原理及特点;以二维泊松方程为例研究了支持域无量纲尺寸、场节点与背景网格设置对无网格法计算精度的影响。将RPIM与EFGM应用于频率域线源二维正演,给出了RPIM形状参数的推荐值;分析了均匀介质模型大地电磁(MT)二维正演无网格法边界条件直接加载与罚函数法加载的精度差异,结合PIM与RPIM边界条件加载便利及EFGM计算复杂模型精度高的优势,提出了EFG—PIM及EFG.RPIM耦合算法,数值计算结果验证了耦合算法的有效性。研究发现:无网格法及其耦合方法适用于电磁法数值模拟;支持域无量纲尺寸取1.0时无网格法精度与效率高,场节点与背景网格重合时计算效果佳;泊松方程求解PIM及RPIM精度较EFGM低,计算均匀介质MT响应精度较EFGM高;RPIM改善了PIM计算涉及的奇异性问题,对应支持域无量纲尺寸选择空间大。
简介:在地震勘探数据采集中,随机噪声严重影响地震资料质量,给后期解释工作带来很大困难。如何在不损失剖面有效信息的前提下压制随机噪声,有效地提高地震资料的信噪比和保真度,是本文的研究目标。构造导向滤波技术的核心是构造方向表征的求取以及如何实现非平稳滤波,来达到提高地震数据信噪比和保真度的目的。本文首先通过分析函数二维导数与希尔伯特变换的频率响应关系,推导出了基于二维希尔伯特变换的非迭代地震同相轴倾角求取算子,进而达到了构造方向表征的求取;其次选取多项式拟合作为构造导向滤波中的非平稳滤波方法,扩展了非平稳多项式拟合的应用范围;最后沿构造倾角方向进行变振幅同相轴的非平稳多项式拟合,实现和构建了新的自适应构造导向滤波方法。理论模型和实际地震资料处理的结果表明,所提出的方法实现了既保护构造信息又有效地压制了随机噪声的目的。
简介:在地震勘探数据采集中,随机噪声严重影响地震资料质量,给后期解释工作带来很大困难。如何在不损失剖面有效信息的前提下压制随机噪声,有效地提高地震资料的信噪比和保真度,是本文的研究目标。构造导向滤波技术的核心是构造方向表征的求取以及如何实现非平稳滤波,来达到提高地震数据信噪比和保真度的目的。本文首先通过分析函数二维导数与希尔伯特变换的频率响应关系,推导出了基于二维希尔伯特变换的非迭代地震同相轴倾角求取算子,进而达到了构造方向表征的求取;其次选取多项式拟合作为构造导向滤波中的非平稳滤波方法,扩展了非平稳多项式拟合的应用范围;最后沿构造倾角方向进行变振幅同相轴的非平稳多项式拟合,实现和构建了新的自适应构造导向滤波方法。理论模型和实际地震资料处理的结果表明,所提出的方法实现了既保护构造信息又有效地压制了随机噪声的目的。