基于高斯束传播算子的成像域层析方法.docxVIP

基于高斯束传播算子的成像域层析方法 1 研究内容介绍 石油和天然气勘探的最终目标是存储层。对于存储层的描述有两种重要方法：一类是偏移速度估算、真气候分析和角度提取、ava（振幅随角度变化率）反演和存储层描述，这是当前方法和技术相对成熟的第一种方法。另一方面，将重叠前的地震数据直接反演为fwi（全波形反演）宽度波数带的弹性参数场直接用于存储层描述。上述两种路线的核心是地震波的反演成像。最终目标是在全频带中建立速度领域，精确测量和精细速度领域的建立可以达到储层的目的。第一条路是石油地震勘探学家发展的实用路线。主要采用nmo（nmo动态校正）梯度分析、psdm（初始速度位移分析）、psdm（初始速度位移分析）、psdm（初始速度位移分析）、psdm（初始速度位移分析）、psdm（初始速度位移分析）、折射波反射波层析和背景速度分析。高波数的变化矩阵（高波数组成的空间位置）。现在，它仍然是石油行业的主要技术路线。一些成熟的技术已经在石油行业很流行。第二条路线是基于数学反演理论框架的路径，主要是基于fwi和lsrtm（最小二乘时间位移位移）。预测源速度，lsrtm估计速度场的高波数变化成分。这是理论上接近全波数带反演的最好方式。目前，它仍处于理论研究和探索和应用阶段。这项工作的研究内容属于第一条路线的反射波层析和估计背景速度。 对现有主流的层析成像方法进行分类,根据执行域可以分为成像域和数据域层析成像;根据正问题可以分为基于波动理论和基于射线理论的层析成像;根据利用数据属性可以分为波形和走时层析成像;根据利用数据类型主要可以分为透射波和反射波层析成像.根据正问题及执行域,分别总结和评价基于波动理论和基于射线理论的、在成像域和数据域实现的速度层析反演方法特点,整体可以概括为:(1)基于射线理论的层析偏移速度反演方法(成像域)发展较为成熟,在实际应用中已经实现自动化层析(层析控制点的自动拾取、RMO(Residue Move Out,剩余时差)的自动计算、利用某些相似性准则自动筛选数据等).然而,射线理论的固有问题(阴影区、多路径问题等)限定了该方法的反演精度和在复杂速度模型中的适用性,因而反演效果仍然有待改善(这也是本文考虑的方向);(2)波动方程偏移速度分析(WEMVA)方法以波动方程作为正演算子,克服了射线理论的高频近似假设;其提取的成像道集(扩展成像条件)可实现自动化反演,可适用复杂模型的速度估计.其反演精度理论上优于射线类层析方法.但是,相较于射线类层析,WEMVA方法的计算量巨大,更存在地下照明不均衡及泛函梯度假象等问题( 上述前两类方法属于成像域速度反演方法,主要利用反射数据的走时信息,对初始模型的依赖性较低.该方法仅能反演速度模型的低波数成分,因而主要用于背景速度建模,为偏移成像服务及为更高精度的反演方法(如FWI)提供初始模型.后两类方法属于数据域速度反演方法,可选择利用走时、相位、振幅或波形属性,理论反演精度较高.但受到叠前数据复杂波场及低信噪比的影响,难以准确高效地实现全自动拾取,实际应用较为困难.从实用化的角度,本文关注成像域走时层析速度反演方法. 常规的成像域走时层析速度反演的基本过程是首先自动拾取成像道集ADCIGs(Angle Domain Common Image Gathers,角度域共成像点道集)或ODCIGs(Offset Domain Common Image Gathers,偏移距域共成像点道集)的剩余时差(RMO),然后通过射线追踪计算射线路径并构造相应的层析线性方程,求解方程实现剩余时差的反投影从而估计速度更新量.通常而言,层析速度反演一般在偏移距域或角度域成像道集中实现.相对于偏移距域,角度域层析更有优势(射线出发角确定,无需费时且不稳定的两点射线追踪匹配地表偏移距,且避免了多路径等问题).基于射线理论的层析速度反演方法技术相对较为成熟,在实际应用中已经基本实现自动化(包括剖面上自动判断/拾取层析控制点、自动计算界面倾角、成像道集上自动筛选数据并计算RMO等).然而,射线理论的固有问题限定了该方法的反演精度、稳定性及适用范围,因而有必要进一步发展更加适合复杂速度模型的射线束层析反演方法. 与成像方法类似,在层析速度反演中也有介于常规射线理论和波动理论之间的方法,如胖射线层析( 2 成像域走时扰动与速度扰动的关系式 从角度域高斯束偏移(GBM)成像条件出发: (1)式是频率域高斯束成像条件,其中pS和pR分别代表从炮点和检波点传播到成像点x的慢度矢量,表示波场在该点的传播方向.I 在波动方程的一阶Born近似下,波场U可以分解为背景波场U 因此从成像条件(1)式可以近似得到扰动像: 其中,ΔS,ΔR分别为一阶Born近似散射场,S 进一步地,根据 (4)式中的V 如果将(4)式代入(3)式,形式