CN107843925B - 一种基于修正相位的反射波波形反演方法 - Google Patents
一种基于修正相位的反射波波形反演方法 Download PDFInfo
- Publication number
- CN107843925B CN107843925B CN201710916080.9A CN201710916080A CN107843925B CN 107843925 B CN107843925 B CN 107843925B CN 201710916080 A CN201710916080 A CN 201710916080A CN 107843925 B CN107843925 B CN 107843925B
- Authority
- CN
- China
- Prior art keywords
- wave
- field
- cost functional
- model
- phase place
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 93
- 238000013508 migration Methods 0.000 claims description 19
- 230000005012 migration Effects 0.000 claims description 19
- 238000002474 experimental method Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- 238000000354 decomposition reaction Methods 0.000 claims description 4
- 230000009189 diving Effects 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 3
- 230000009471 action Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 9
- 238000004804 winding Methods 0.000 abstract description 8
- 238000000926 separation method Methods 0.000 abstract description 3
- 238000000614 phase inversion technique Methods 0.000 abstract description 2
- 238000003672 processing method Methods 0.000 abstract description 2
- 238000002939 conjugate gradient method Methods 0.000 description 4
- 238000007689 inspection Methods 0.000 description 4
- 238000005457 optimization Methods 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000004587 chromatography analysis Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 239000000729 antidote Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 238000002945 steepest descent method Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及地震资料的处理方法,具体为全波形反演速度场建模方法。提出了一种基于修正相位目标泛函的反射波波形反演方法。通过修改传统的相位计算方法得到新的目标泛函,利用伴随状态法求得相应的伴随震源;通过最优化方法对目标泛函寻优以获得有效的背景速度场;最终与传统的波形反演结合得到高精度的速度场建模结果。本发明针对传统相位反演方法不稳定问题,避免了传统的相位缠绕问题,通过相位与振幅信息的分离,所对应的目标泛函具有更高的全局收敛性。能够有效更新模型深层的中低波数组分,为全波形反演提供准确的初始模型,有效降低全波形反演对初始模型的依赖性。
Description
技术领域
本发明涉及地震资料处理方法,具体是地震资料的波形反演速度建模方法。
背景技术
由于模型参数众多和正演问题的复杂性,FWI通常采用梯度类优化方法,即从足够准确的初始模型开始,利用伴随状态法求取目标泛函在当前模型的梯度信息并对模型进行更新直至目标泛函收敛。目标泛函通常用以衡量观测数据与正演数据的误差,其选择合适与否直接决定的波形反演问题能否成功。Tejero等(2015)详细对比了几种时间域全波形反演目标泛函性态,指出基于相位信息目标泛函的反演方法对低频信息缺失和噪音具有更高的稳定性。基于瞬时相位目标泛函的反演方法分别在时频域和时间域应用于区域地震数据。然而,传统的相位反演存在固有的缠绕问题。
准确的中低波数成分是波形反演成功的重要保障。众所周知,受采集能力和处理技术的限制,地震数据主要以近偏移距的反射波为主。为了能够有效提取由反射波所携带的地下介质中低波数组分,Xu等(2012)提出了反射波波形反演(RFWI)方法。由于数据振幅、走时信息与模型之间的非线性关系,以L2范数为目标泛函的反射波波形反演同样是极强的病态问题。振幅与走时信息的分离可以有效提高RFWI线性程度。为了提高反射波波形反演的适用性,Chi等(2015)提出了基于时间和空间互相关的目标泛函,在实际资料处理上取得了较好的效果。Wang等人(2015)提出了利用加窗互相关拾取走时的方法以提高RFWI稳定性。将动态图像矫正方法应用于反射波旅行时差提取上,Ma&Hale(2013)利用反射波走时信息可以获得较为准确的中低波数模型组分。
发明内容
本发明的目的在于,提供一种基于修正相位的反射波波形反演方法,一方面通过定义新的相位求取方法有效避免传统定义方法的相位缠绕问题;另一方面通过相位信息定义反射波波形反演以提高其目标泛函的线性程度,提高反演稳定性。
为了解决上述问题,本发明所采用的详细技术方案如下。
一种基于修正相位的反射波波形反演方法,该方法包括如下步骤:
步骤1:利用地震数据的包络信息修正相位的表达式,求取地震数据修正后的相位信息;
步骤2:定义基于修正相位信息的反射波波形反演的目标泛函表达式,利用伴随状态法获得目标泛函关于模型的梯度表达式;
步骤3:利用获得的梯度对模型背景场进行更新,获得走时信息的准确的模型背景场,以保证走时信息的准确;
步骤4:最后利用获得走时信息的准确的模型背景场做常规的全波形反演,获得高精度的速度场建模结果。
进一步,在步骤1中,利用地震数据的包络信息修正相位的表达式的具体方法是:
将相位的表达式修改为原信号与信号包络比值的反正切值
其中atan表示反正切,u为地震数据,E(u)为地震信号的包络,θ(u)为相位。
进一步,步骤2中,所述目标泛函为基于L2范数的目标泛函。。
进一步,步骤2中,所述目标泛函为为将修正的相位反演方法与反射波波形反演(RFWI)方法相结合,定义相应的目标泛函。
进一步,步骤2中,所述基于修正相位信息的反射波波形反演的目标泛函为:
其中χ表示目标泛函,θ(dref),θ(δu)分别表示t时刻在震源点s和检波点r处的反射波观测数据和正演数据的瞬时相位信息。将模型分解为背景速度场m和扰动速度场δm,则对应的波场可以分解为背景波场u和扰动波场δu。
进一步,步骤2中,通过对正传震源波场和反传伴随波场的互相关求取目标泛函关于模型的梯度。
进一步,目标泛函的梯度表示为:
其中χ表示目标泛函,B表示正演算子,表示微分算子,m和δm分别表示背景速度场扰动速度场,u和δu分别表示背景波场扰动波场,则对应的波场可以分解为背景波场u和扰动波场δu表示模型的扰动,值得说明的是与前面的δm含义是不同的,在实际应用中,δm往往由偏移或近偏移距的FWI获得,代表模型中的高波数组分,具有明确的物理意义;而表示与目标泛函梯度的点积项,仅仅是数学上的含义,没有明确的物理意义;和表示伴随波场。伴随方程可以表示为:
其中表示伴随正演算子。
进一步,对应于伴随方程的伴随震源可以表示为:
其中,Δθ(t)表示相位数据残差,E(δu)表示反射波数据包络,H表示希尔伯特变换,ξ为保证稳定性的极小值。
进一步,步骤3中,所述利用获得的梯度对模型背景场进行更新具体包括如下步骤。
步骤A::获得观测数据,在数值实验中,观测数据通过在真实模型中正演生成,利用声波有限差分做正演获得地震记录作为观测记录,或者在实际应用中,观测记录由野外采集得到。
步骤B:选取初始速度场,波形反演方法通过对数据匹配以更新速度场,因此除数据外,还需要输入初始速度场,在数值实验中,初始模型通过对真实模型平滑得到,或者在实际数据应用中,初始速度场通过走时层析处理得到;在本方法中窗口的位置由直达波的走时确定,直达波的走时通过炮检点距离以及模型浅层速度确定;
步骤C:由于本发明中反演方法只对反射波进行反演,需要提取观测数据中反射波数据,利用加窗方法去除观测记录中的直达波、回折波和首波,只保留反射波成分;在本方法中窗口的位置由直达波的走时确定,直达波的走时通过炮检点距离以及模型浅层速度确定;
步骤D:对反射波做近偏移距的偏移,获得反射系数场δm;然后利用生成的反射系数场做反偏移获得正演反射波数据;将正演数据与观测数据代入目标泛函中获得目标泛函值。
步骤E:利用伴随震源公式求取伴随源,分别利用两个伴随波动方程求得伴随波场,伴随波场是由伴随状态法定义的,其物理意义是通过对伴随震源的逆时反传获得的,除了与传播的时间方向不同外,伴随波场利用的是终值条件,而传统的波动方程利用的初值条件;根据伴随波场与正传波场互相关得到目标泛函关于模型的梯度,这里的互相关是由公式3定义的,可见目标泛函关于模型参数的梯度是由正传波场与反传伴随波场的互相关获得的。
步骤F:对模型背景场进行优化更新,判断目标泛函是否收敛,如果目标泛函收敛,流程进入步骤G,否则重复步骤D-E,直至目标泛函收敛为止。
步骤G:将目标泛函收敛时所获得的背景速度场作为偏移速度场,与利用真实速度场得到的偏移结果相比较,如果符合要求,即认为该背景速度场为走时信息的准确的模型背景场;否则,重复步骤D-F,直至符合要求为止。
进一步,步骤F中,对模型背景场进行优化更新采用梯度类方法,梯度类方法包括共轭梯度法、BFGS法以及牛顿法等,受计算条件限制,通常采用共轭梯度法和L-BFGS法,对模型的背景场进行更新以保证走时信息的准确性。
进一步,步骤F中,对模型背景场进行优化更新采用共轭梯度法,共轭梯度法是传统最速下降法的修正,其利用当前梯度方向与上一步的更新方向加权组成当前的搜索方向,能够在不增加计算量的同时明显提升计算效率。
进一步,步骤F中,对模型背景场进行优化更新采用BFGS法,BFGS法是对牛顿法的近似,其利用多次梯度方向组成对汉森矩阵(目标泛函的二阶导数)的近似,计算中仅涉及矩阵的相乘和相加运算,并且避免了对汉森矩阵的直接存储,在多次迭代的情况下能够得到与牛顿法接近的迭代效率,有效减少了计算存储量,提升了计算效率。
进一步,步骤G中,当真实速度场得到的偏移结果与目标泛函收敛时所获得的背景速度场之间的差距小于给定阈值时,判定为符合要求。
进一步,步骤4中,将获得的背景速度场作为常规全波形反演的输入初始速度场,对模型的细节进行精细刻画,获得高精度的速度场建模结果。
本发明的的有益效果可以总结为:①相比于传统的相位计算方法,本发明中提出的修正的相位计算方法避免了传统的缠绕问题,所对应的目标泛函具有更高的全局收敛性。②RFWI能够有效更新模型深层的中低波数组分,为FWI提供准确的初始模型,从而有效降低全波形反演对初始模型的依赖性。③将修正的相位目标泛函与传统的RFWI相结合,能够分离地震数据中振幅与相位信息,有效提高RFWI线性程度。避免了采用传统的基于L2范数的RFWI时每次循环做最小二乘偏移的要求。
附图说明
图1为修正后的相位计算方法避免相位缠绕的试算,其中(a)为两个时间域信号,(b)为传统相位计算方法获得的相位数据(c)为本发明中相位计算方法获得的相位数据;
图2为采用的真实模型,为国际标准Sigsbee2A模型中的一部分;
图3为反演采用的初始模型;
图4为第一次迭代中形成的反射系数场;
图5为更新后的背景速度场;
图6为利用更新后的速度场得到的偏移结果;
图7为利用真实速度场得到的偏移结果;
图8为利用更新后的速度场做传统的全波形反演得到的反演结果;
图9为直接利用传统全波形反演方法得到的反演结果;
图10为传统全波形反演方法(FWI)和利用Terejo定义的相位函数(RFWI(RIAR)+FWI)与本文中定义的相位函数(RFWI(RRE)+FWI)更新速度场后的全波形反演的模型误差的收敛曲线;
图11整体反演流程图。
具体实施方式
为了使本发明的上述和其他目的、特点和优点能够更加明显易懂,下文特举出较佳的实施实例,并结合附图和相应理论,做详细说明如下。
实施例1。一种基于修正相位的反射波波形反演方法,该方法包括如下步骤:
步骤1:利用地震数据的包络信息修正相位的表达式,求取地震数据修正后的相位信息;
步骤2:定义基于修正相位信息的反射波波形反演的目标泛函表达式,利用伴随状态法获得目标泛函关于模型的梯度表达式;
步骤3:利用获得的梯度对模型背景场进行更新,获得走时信息的准确的模型背景场,以保证走时信息的准确;
步骤4:最后利用获得走时信息的准确的模型背景场做常规的全波形反演,获得高精度的速度场建模结果。
在步骤1中,利用地震数据的包络信息修正相位的表达式的具体方法是:
将相位的表达式修改为原信号与信号包络比值的反正切值
其中atan表示反正切,u为地震数据,E(u)为地震信号的包络,θ(u)为相位。
步骤2中,所述目标泛函为基于L2范数的目标泛函。所述目标泛函为:
其中χ表示目标泛函,θ(dref),θ(δu)分别表示t时刻在震源点s和检波点r处的反射波观测数据和正演数据的瞬时相位信息。将模型分解为背景速度场m和扰动速度场δm,则对应的波场可以分解为背景波场u和扰动波场δu。
目标泛函的梯度表示为:
其中χ表示目标泛函,B表示正演算子,表示微分算子,m和δm分别表示背景速度场扰动速度场,u和δu分别表示背景波场扰动波场,则对应的波场可以分解为背景波场u和扰动波场δu表示模型的扰动,值得说明的是与前面的δm含义是不同的,在实际应用中,δm往往由偏移或近偏移距的FWI获得,代表模型中的高波数组分,具有明确的物理意义;而表示与目标泛函梯度的点积项,仅仅是数学上的含义,没有明确的物理意义;和表示伴随波场;伴随方程可以表示为:
其中表示伴随正演算子。
进一步,对应于伴随方程的伴随震源可以表示为:
其中,Δθ(t)表示相位数据残差,E(δu)表示反射波数据包络,H表示希尔伯特变换,ξ为保证稳定性的极小值。。
步骤3中,所述利用获得的梯度对模型背景场进行更新具体包括如下步骤:
步骤A:采用如图2所示的速度场作为真实速度场,利用声波有限差分做正演获得地震记录作为观测记录。
步骤B:采用如图3所示的速度场作为初始速度场。在本数值试验中初始模型由对真实模型平滑得到,更加符合实际应用情况,在实际应用中初始模型一般由走时层析等常规处理方法得到,精度较低;
步骤C:利用加窗方法去除观测记录中的直达波、回折波和首波,只保留反射波成分。在本方法中窗口的位置由直达波的走时确定,直达波的走时通过炮检点距离以及模型浅层速度确定。
步骤D:对反射波做近偏移距的偏移,获得反射系数场δm;然后利用生成的反射系数场做反偏移获得正演反射波数据;将正演数据与观测数据代入目标泛函中获得目标泛函值;
步骤E:利用伴随震源公式求取伴随源,分别利用两个伴随波动方程求得伴随波场;根据伴随波场与正传波场互相关得到目标泛函关于模型的梯度;
步骤F:对模型背景场进行优化更新,判断目标泛函是否收敛,如果目标泛函收敛,流程进入步骤G,否则重复步骤D-E,直至目标泛函收敛为止;
步骤G:将目标泛函收敛时所获得的背景速度场作为偏移速度场,与利用真实速度场得到的偏移结果相比较,如果符合要求,即认为该背景速度场为走时信息的准确的模型背景场;否则,重复步骤D-F,直至符合要求为止。
步骤F中,对模型背景场进行优化更新采用梯度类方法,对模型的背景场进行更新以保证走时信息的准确性。
本实施例提出了一种基于修正相位目标泛函的反射波波形反演方法。通过修改传统的相位计算方法得到新的目标泛函,利用伴随状态法求得相应的伴随震源;通过最优化方法对目标泛函寻优以获得有效的背景速度场;最终与传统的波形反演结合得到高精度的速度场建模结果。本发明针对传统相位反演方法不稳定问题,避免了传统的相位缠绕问题,通过相位与振幅信息的分离,所对应的目标泛函具有更高的全局收敛性。能够有效更新模型深层的中低波数组分,为全波形反演提供准确的初始模型,有效降低全波形反演对初始模型的依赖性。
实施例2。一种基于修正相位的反射波波形反演方法,该方法首先提出了一种修正的相位反演方法,利用地震数据的包络信息来有效避免传统相位反演的缠绕问题;然后将修正的相位反演方法与RFWI方法相结合,定义相应的目标泛函,并且利用伴随状态法求得该目标泛函的伴随方程形式;通过对正传震源波场和反传伴随波场的互相关求取目标泛函关于模型的梯度;利用优化方法,例如梯度类方法等对模型的背景场进行更新以保证走时信息的准确性;最后将获得的背景速度场作为传统波形反演的输入初始速度场,对模型的细节进行精细刻画,以获得高精度的地下速度模型。
具体包括如下步骤。
(1)首先将相位公式修改为原信号与信号包络比值的反正切值
其中atan表示反正切,u为地震数据,E(u)为地震信号的包络,如图1所示该定义能够有效避免相位的缠绕现象。
(2)定义基于修正相位信息的反射波波形反演的目标泛函:
其中χ表示目标泛函,θ(dref),θ(δu)分别表示t时刻在震源点s和检波点r处的反射波观测数据和正演数据的瞬时相位信息。将模型分解为背景速度场m和扰动速度场δm,则对应的波场可以分解为背景波场u和扰动波场δu;
(3)则对应于目标泛函的梯度可以表示为:
其中χ表示目标泛函,B表示正演算子,表示微分算子,m和δm分别表示背景速度场扰动速度场,u和δu分别表示背景波场扰动波场,则对应的波场可以分解为背景波场u和扰动波场δu表示模型的扰动,值得说明的是与前面的δm含义是不同的,在实际应用中,δm往往由偏移或近偏移距的FWI获得,代表模型中的高波数组分,具有明确的物理意义;而表示与目标泛函梯度的点积项,仅仅是数学上的含义,没有明确的物理意义;和表示伴随波场;伴随方程可以表示为:
其中表示伴随正演算子。
(4)对应于伴随方程的伴随震源可以表示为:
(5)采用如图2所示的速度场作为真实速度场,利用声波有限差分做正演,获得地震记录作为观测记录。由于是数值实验,因此仅利用真实模型正演获得的数据作为观测数据作为反演的输入,真实模型不参与反演过程,仅产生正演数据。真实模型为Sigsbee2A模型的一部分,由SMAART JV公开,其中包含岩丘,平层等,是测试算法的标准模型之一。
(6)采用如图3所示的速度场作为初始速度场,在本数值试验中初始模型由对真实模型平滑得到,更加符合实际应用情况(如图3所示),由于平滑得到的初始速度场是对真实模型的近似,并且由于平滑的窗口较大,速度场的精度较低,因此与实际应用中利用走时层析或者速度分析等方法获得的初始速度场相似。在实际应用中初始模型一般由走时层析等常规处理方法得到,精度较低。
(7)利用加窗方法去除观测记录中的直达波、回折波和首波,只保留反射波成分,在本方法中窗口的位置由直达波的走时确定,直达波的走时通过炮检点距离以及模型浅层速度确定。
(8)对反射波做近偏移距的偏移,获得反射系数场δm,如图4所示;然后利用生成的反射系数场做反偏移获得正演反射波数据;将正演数据与观测数据代入目标泛函中获得目标泛函值。
(9)利用伴随震源公式求取伴随源,分别利用两个伴随波动方程求得伴随波场;伴随波场与正传波场互相关可以得到目标泛函关于模型的梯度。
(10)利用优化方法,例如共轭状态法或BFGS方法对模型背景场进行更新;
(11)重复步骤(8)-(10)直至目标泛函收敛。目标泛函收敛时所获得的背景速度场如图5所示。该背景速度场可以作为偏移速度场,偏移结果如图6所示,与利用真实速度场得到的偏移结果(图7)相比较说明得到的背景场较为准确。
(12)将获得的背景速度场作为常规全波形反演的输入初始速度场,进行精细模型的刻画。更新后的速度场经过全波形反演后的结果如图8所示,相比于直接利用传统的全波形反演进行速度场建模(图9),该方法效果有了明显改善。对比两者的模型残差收敛曲线(图10)可见本发明中方法有效地恢复了速度场,整体的反演流程如图11所示。
本实施例利用地震数据的包络信息,提出了一种新的相位定义方法,可以避免传统相位求取中的缠绕问题。通过提取地震数据的相位信息,可以避免振幅信息的复杂性,增强反射波波形反演的稳定性。通过伴随状态法求得了以修正相位信息所定义的波形反演伴随震源,形成了相应的反演理论基础。
Claims (13)
1.一种基于修正相位的反射波波形反演方法,其特征在于,该方法包括如下步骤:
步骤1:利用地震数据的包络信息修正相位的表达式,求取地震数据修正后的相位信息;
步骤2:定义基于修正相位信息的反射波波形反演的目标泛函表达式,利用伴随状态法获得目标泛函关于模型的梯度表达式;
步骤3:利用获得的梯度对模型背景场进行更新,获得走时信息的准确的模型背景场,以保证走时信息的准确;
步骤4:最后利用获得走时信息的准确的模型背景场做常规的全波形反演,获得高精度的速度场建模结果;
在步骤1中,利用地震数据的包络信息修正相位的表达式的具体方法是:
将相位的表达式修改为地震数据与地震信号的包络比值的反正切值
其中atan表示反正切,p为地震数据,E(p)为地震信号的包络,θ(p)为相位。
2.如权利要求1所述的基于修正相位的反射波波形反演方法,其特征在于:步骤2中,所述目标泛函为基于L2范数的目标泛函。
3.如权利要求1所述的基于修正相位的反射波波形反演方法,其特征在于:步骤2中,所述目标泛函为将修正的相位反演方法与反射波波形反演方法相结合,定义相应的目标泛函。
4.如权利要求1所述的基于修正相位的反射波波形反演方法,其特征在于:步骤2中,所述基于修正相位信息的反射波波形反演的目标泛函为:
其中χ表示目标泛函,θ(dref),θ(δu)分别表示t时刻在震源点s和检波点r处的反射波观测数据和正演数据的瞬时相位信息;将模型分解为背景速度场m和扰动速度场δm,则对应的波场可以分解为背景波场u和扰动波场δu,扰动波场作为反射波数据。
5.如权利要求1所述的基于修正相位的反射波波形反演方法,其特征在于:
步骤2中,通过对正传震源波场和反传伴随波场的互相关求取目标泛函关于模型的梯度。
6.如权利要求5所述的基于修正相位的反射波波形反演方法,其特征在于:目标泛函的梯度表示为:
其中χ表示目标泛函,B表示正演算子,表示微分算子,m和δm分别表示背景速度场扰动速度场,u和δu分别表示背景波场和扰动波场,则对应的波场可以分解为背景波场u和扰动波场δu,表示模型的扰动,值得说明的是与前面的δm含义是不同的,在实际应用中,δm往往由近偏移距的FWI获得,代表模型中的高波数组分,具有明确的物理意义;而表示与目标泛函梯度的点积项,仅仅是数学上的含义,没有明确的物理意义;和表示伴随波场;伴随方程可以表示为:
其中表示伴随正演算子,Sadj表示伴随震源。
7.如权利要求4所述的基于修正相位的反射波波形反演方法,其特征在于:
对应于公式(1)所定义的相位信息和公式(2)所定义的目标泛函的伴随震源可以表示为:
其中,Δθ(t)表示相位数据残差,E(δu)表示反射波数据包络,H表示希尔伯特变换,ξ为保证稳定性的极小值。
8.如权利要求2所述的基于修正相位的反射波波形反演方法,其特征在于,步骤3中,所述利用获得的梯度对模型背景场进行更新具体包括如下步骤:
步骤A:获得观测记录,在数值实验中,观测记录通过在真实模型中正演生成,利用声波有限差分做正演获得地震记录作为观测记录,或者在实际应用中,观测记录由野外采集得到;
步骤B:选取初始速度场,波形反演方法通过对数据匹配以更新速度场,因此除数据外,还需要输入初始速度场,在数值实验中,初始速度场通过对真实模型平滑得到,或者在实际数据应用中,初始速度场通过走时层析处理得到;
步骤C:由于反演方法只对反射波进行反演,需要提取观测数据中反射波数据,利用加窗方法去除观测记录中的直达波、回折波和首波,只保留反射波成分;
步骤D:对反射波做近偏移距的偏移,获得反射系数场;然后利用生成的反射系数场做反偏移获得正演反射波数据;将正演反射波数据与观测数据中截取得到的反射波代入目标泛函中获得目标泛函值;
步骤E:利用伴随震源公式求取伴随源,分别利用两个伴随波动方程求得伴随波场;根据伴随波场与正传波场互相关得到目标泛函关于模型的梯度;
步骤F:对模型背景场进行优化更新,判断目标泛函是否收敛,如果目标泛函收敛,流程进入步骤G,否则重复步骤D-E,直至目标泛函收敛为止;
步骤G:将目标泛函收敛时所获得的背景速度场作为偏移速度场,与利用真实速度场得到的偏移结果相比较,如果符合要求,即认为该背景速度场为走时信息的准确的模型背景场。
9.如权利要求8所述的基于修正相位的反射波波形反演方法,其特征在于:步骤F中,对模型背景场进行优化更新采用梯度类方法,对模型背景场进行更新以保证走时信息的准确性。
10.如权利要求8所述的基于修正相位的反射波波形反演方法,其特征在于:步骤F中,对模型背景场进行优化更新采用共轭状态法。
11.如权利要求8所述的基于修正相位的反射波波形反演方法,其特征在于:步骤F中,对模型背景场进行优化更新采用BFGS法。
12.如权利要求8所述的基于修正相位的反射波波形反演方法,其特征在于:步骤G中,当真实速度场得到的偏移结果与目标泛函收敛时所获得的背景速度场之间的差距小于给定阈值时,判定为符合要求。
13.如权利要求8所述的基于修正相位的反射波波形反演方法,其特征在于:步骤4中,将获得的背景速度场作为常规全波形反演的输入初始速度场,对模型的细节进行那个精细刻画,获得高精度的速度场建模结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710916080.9A CN107843925B (zh) | 2017-09-29 | 2017-09-29 | 一种基于修正相位的反射波波形反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710916080.9A CN107843925B (zh) | 2017-09-29 | 2017-09-29 | 一种基于修正相位的反射波波形反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107843925A CN107843925A (zh) | 2018-03-27 |
CN107843925B true CN107843925B (zh) | 2019-03-08 |
Family
ID=61661542
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710916080.9A Active CN107843925B (zh) | 2017-09-29 | 2017-09-29 | 一种基于修正相位的反射波波形反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107843925B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108680957B (zh) * | 2018-05-21 | 2019-08-13 | 吉林大学 | 基于加权的局部互相关时频域相位反演方法 |
CN108873066B (zh) * | 2018-06-26 | 2020-03-10 | 中国石油大学(华东) | 弹性介质波动方程反射波旅行时反演方法 |
CN109655910B (zh) * | 2019-01-18 | 2019-12-24 | 吉林大学 | 基于相位校正的探地雷达双参数全波形反演方法 |
CN111665551B (zh) * | 2019-03-07 | 2023-05-02 | 中普宝信(北京)科技有限公司 | 用于桥梁基底探测的声学参数获取方法 |
CN112444848B (zh) * | 2019-08-29 | 2024-01-23 | 中国石油化工股份有限公司 | 一种全波形反演的方法及系统 |
CN112630830B (zh) * | 2019-10-08 | 2024-04-09 | 中国石油化工股份有限公司 | 一种基于高斯加权的反射波全波形反演方法及系统 |
CN110687594B (zh) * | 2019-10-16 | 2020-11-24 | 中国石油大学(华东) | 一种瞬时相位解缠绕方法、全波形反演方法及计算机设备 |
CN110764146B (zh) * | 2019-10-24 | 2020-05-19 | 南京信息工程大学 | 一种基于声波算子的空间互相关弹性波反射波形反演方法 |
CN111060960B (zh) * | 2019-12-27 | 2022-03-18 | 恒泰艾普(北京)能源科技研究院有限公司 | 一种基于合成炮记录的fwi建模方法 |
CN113156493B (zh) * | 2021-05-06 | 2022-02-18 | 中国矿业大学 | 一种使用归一化震源的时频域全波形反演方法及装置 |
CN113050163B (zh) * | 2021-05-06 | 2022-05-03 | 中国矿业大学 | 一种振幅、相位信息可调节的全波形反演方法及装置 |
US11867857B2 (en) | 2021-07-13 | 2024-01-09 | Saudi Arabian Oil Company | Method and system for updating a seismic velocity model |
CN116660981B (zh) * | 2023-07-25 | 2023-10-24 | 北京中矿大地地球探测工程技术有限公司 | 基于包络的各向异性参数反演方法、装置及存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104391323A (zh) * | 2014-11-21 | 2015-03-04 | 中国石油大学(华东) | 一种利用反射波信息反演速度场中低波数成分的方法 |
CN104977614A (zh) * | 2014-04-09 | 2015-10-14 | 中国石油集团东方地球物理勘探有限责任公司 | 一种基于相邻频率相位差目标函数的频率域全波形反演方法 |
CN105467444A (zh) * | 2015-12-10 | 2016-04-06 | 中国石油天然气集团公司 | 一种弹性波全波形反演方法及装置 |
-
2017
- 2017-09-29 CN CN201710916080.9A patent/CN107843925B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104977614A (zh) * | 2014-04-09 | 2015-10-14 | 中国石油集团东方地球物理勘探有限责任公司 | 一种基于相邻频率相位差目标函数的频率域全波形反演方法 |
CN104391323A (zh) * | 2014-11-21 | 2015-03-04 | 中国石油大学(华东) | 一种利用反射波信息反演速度场中低波数成分的方法 |
CN105467444A (zh) * | 2015-12-10 | 2016-04-06 | 中国石油天然气集团公司 | 一种弹性波全波形反演方法及装置 |
Non-Patent Citations (5)
Title |
---|
Comparative study of objective functions to overcome noise and bandwidth limitations in full waveform inversion;C.E. Jimenez Tejero et al.;《Geophysical Journal International》;20151231;第634-635页 |
Inversion on Reflected Seismic Wave;Sheng Xu et al.;《Inversion on Reflected Seismic Wave》;20121231;第1-7页 |
Misfit functions for full waveform inversion based on instantaneous phase and envelope measurements;Ebru Bozdag et al.;《Geophysical Journal International》;20111231;第848-849页 |
RFWI with Source Independent Inversion Strategy;Jianping Huang et al.;《2017 CGS/SEG International Geophysical Conference》;20170420;第303-305页,图2、5-6 |
反射波全波形反演;姚刚等;《中国科学:地球科学》;20170919;第47卷(第10期);第1220-1232页 |
Also Published As
Publication number | Publication date |
---|---|
CN107843925A (zh) | 2018-03-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107843925B (zh) | 一种基于修正相位的反射波波形反演方法 | |
CN103703391B (zh) | 使用频谱整形的全波场反演的系统和计算机实施的方法 | |
KR101948509B1 (ko) | 지구 물리학적 데이터의 반복 반전의 아티팩트 감소 | |
Alkhalifah | Full-model wavenumber inversion: An emphasis on the appropriate wavenumber continuation | |
CN103238158B (zh) | 利用互相关目标函数进行的海洋拖缆数据同时源反演 | |
Guasch et al. | Adaptive waveform inversion: Practice | |
CN102884447A (zh) | Q层析成像方法 | |
US20130258810A1 (en) | Method and System for Tomographic Inversion | |
CN113740901B (zh) | 基于复杂起伏地表的陆上地震数据全波形反演方法及装置 | |
CN113156493B (zh) | 一种使用归一化震源的时频域全波形反演方法及装置 | |
Araujo et al. | Symplectic scheme and the Poynting vector in reverse-time migration | |
CN110187382A (zh) | 一种回折波和反射波波动方程旅行时反演方法 | |
Tognarelli et al. | Two-grid stochastic full waveform inversion of 2D marine seismic data | |
Calderón Agudo et al. | Addressing viscous effects in acoustic full-waveform inversion | |
Jia et al. | Superwide-angle one-way wave propagator and its application in imaging steep salt flanks | |
CN111999764B (zh) | 基于时频域目标函数的盐下构造最小二乘逆时偏移方法 | |
CN112630830B (zh) | 一种基于高斯加权的反射波全波形反演方法及系统 | |
Yu et al. | Application of weighted early-arrival waveform inversion to shallow land data | |
CN108680957B (zh) | 基于加权的局部互相关时频域相位反演方法 | |
Feng et al. | Automatic traveltime inversion via sparse decomposition of seismic data | |
CN111208568B (zh) | 一种时间域多尺度全波形反演方法及系统 | |
CN115170428A (zh) | 一种声波远探测成像图的降噪方法 | |
CN110376642B (zh) | 一种基于锥面波的三维地震速度反演方法 | |
Qiu et al. | Mitigating the cycle-skipping of full-waveform inversion: An optimal transport approach with exponential encoding | |
Díaz et al. | Data domain wavefield tomography using local correlation functions |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |