CN104123440B - 一种在振幅不匹配时使时间域地震波形反演收敛的方法 - Google Patents
一种在振幅不匹配时使时间域地震波形反演收敛的方法 Download PDFInfo
- Publication number
- CN104123440B CN104123440B CN201410321860.5A CN201410321860A CN104123440B CN 104123440 B CN104123440 B CN 104123440B CN 201410321860 A CN201410321860 A CN 201410321860A CN 104123440 B CN104123440 B CN 104123440B
- Authority
- CN
- China
- Prior art keywords
- data
- amplitude
- seismic
- amplitude curve
- observation
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
一种在振幅不匹配时使时间域地震波形反演收敛的方法,涉及石油勘探技术,通过以下步骤实现:1)野外激发人工地震波,由检波器接收并记录地震波场,得到地震观测数据dij(t)=d(xi,xj,t),其中:xi和xj,分别代表震源和接收点的空间位置,t代表观测时间;2)基于初始速度模型V,用有限差分算法模拟相应的地震观测数据,得到地震模拟数据uij(t)=u(xi,xj,t),其中:xi和xj,分别代表震源和接收点的空间位置,t代表观测时间;3)基于反演使用频段,应用以下公式求取时间窗口长度N,其中Δt是时间采样率,fd是当前反演频段的主频;4)对观测数据dij(t)求取时间窗口长度N的振幅曲线aij(t);5)对模拟数据uij(t)求取时间窗口长度N的振幅曲线bij(t);本发明反演过程收敛并取得可靠结果,计算简单,可操作性强。
Description
技术领域
本发明涉及石油勘探技术,用于地震数据深度建模中的时间域地震波形反演,是一种在振幅不匹配时使时间域地震波形反演收敛的方法
背景技术
地震数据处理的重要目的之一是通过深度偏移(Depth Migration)得到高质量的深度域图像。深度域图像的质量依赖于地震速度模型的质量,因此建模技术在整个地震数据处理起着举足轻重的作用。全波形反演(Full Waveform Inversion)是一种数据域(DataDomain)的地震建模技术,其基本原理是通过比较正演得到的模拟数据和实际观测数据的差别来确定模型改变量,反复迭代此过程直到模拟数据和观测数据的差别达到最小。和传统的基于图像域(Image Domain)的建模方法(如层析反演)相比,全波形反演方法的优点在于(1)分辨率高;
(2)需要的人工干预少。
如前所述,全波形反演的目标是通过反复迭代最小化模拟数据和观测数据之间的差别。用数学的语言说,全波形反演是一个最小化问题,其目标函数是模拟数据和观测数据的差。在此,目标函数可以理解为一个“判断迭代收敛的准则”,每一轮迭代后,都通过这个准则判断此轮迭代后,模型的误差在增大或是减小。
最简单的判断准则是采用以下的目标函数L2范数:
其中:d代表模拟数据,u代表观测数据,脚标i表示地震道,j表示样点序号。通过该目标函数比较正演得到的模拟数据和实际观测数据的差别来确定模型改变量,反复迭代此过程,目标函数值减小则说明收敛。
但在实际应用中,判断迭代收敛的准则采用这类简单的目标函数很难获得成功,其中一个重要的问题是模拟数据和观测数据在振幅上的不匹配,其原因包括以下几种:
(1)检波器通过与地表介质的耦合把当地的机械振动转换成电信号,经过数字化后记录下来,所记录的幅度与检波器所处的耦合条件及其变化关系很大,存在不可靠性。
(2)地震振幅直接与震源能量与震源周边介质物理特性有关,野外震源条件与环境的变化会导致观测数据中振幅的变化。
(3)波场的数字模拟过程是基于简化的物理模型,与实测数据所在的物理环境有一定的差异,可以导致观测数据与模拟数据在振动幅度不同。
(4)观测数据中包含一定量的测量,当噪音中包含较低频率成份时,可导致振幅偏离理论值。
观测数据与模拟数据在振幅上的不匹配使原有目标函数中的残差不是直接与速度模型的误差有关,其后果是导致反演过程不收敛或失败。
为了克服这些困难,业界研究人员设计了以下类型的目标函数,
(1)Luo等人提出了旅行时反演方法,其目标是最小化地震波event之间的旅行时之差,这个方法理论上解决了所有问题,在实际应用中其困难之处在于实际数据中存在很强的噪声,导致地震波Event难以分辨拾取。
(2)Shin等人提出的Laplace-Fourier域的全波形反演配合Phase Only目标函数。它的优点是:(1)通过Laplace-Fourier变换通过引入时间域的指数衰减,能起到很强的去随机噪声效果;(2)phase only目标函数能去除振幅的影响;(3)计算简单。该方法的局限在于它只便于用在频率域反演中。
发明目的
本发明的目的是提供一种保证反演过程收敛并取得可靠结果,计算简单,可操作性强的一种在振幅不匹配时使时间域地震波形反演收敛的方法。
发明内容
本发明通过以下步骤实现:
1)野外激发人工地震波,由检波器接收并记录地震波场,得到地震观测数据dij(t)=d(xi,xj,t),
其中:xi和xj,分别代表震源和接收点的空间位置,t代表观测时间;
2)基于初始速度模型V,用有限差分算法模拟相应的地震观测数据,得到地震模拟数据uij(t)=u(xi,xj,t),
其中:xi和xj,分别代表震源和接收点的空间位置,t代表观测时间;
3)基于反演使用频段,应用以下公式求取时间窗口长度N,其中Δt是时间采样率,fd是当前反演频段的主频;
4)对观测数据dij(t)求取时间窗口长度N的振幅曲线aij(t);
5)对模拟数据uij(t)求取时间窗口长度N的振幅曲线bij(t);
所述的步骤4)和步骤5)求取振幅曲线采用均方根振幅曲线公式或高斯窗振幅曲线公式。
所述的均方根振幅曲线公式是:
其中:N为时窗大小,Δt为时间采样率,n为求和公式中的变量取值从到t代表观测时间。
所述的高斯窗振幅曲线公式是:
其中:N为时窗大小,Δt为时间采样率,σ是一个小于等于0.5的常数,n为求和公式中的变量取值从到t代表观测时间。
6)对观测数据作归一化处理:Dij(t)=dij(t)/aij(t);
其中:dij(t)为观测数据,aij(t)为步骤4)所求振幅曲线;
7)对模拟数据作归一化处理:Uij(t)=uij(t)/bij(t);
其中:uij(t)为模拟数据,bij(t)为步骤5)所求振幅曲线;
8)设立目标函数,通过高斯-牛顿优化方法使得野外记录的地震波场和正演模拟的地震波场误差减小,多次迭代反演使目标函数最小化,达到时间域地震波形反演收敛。
步骤8)所述的目标函数是:
其中:Dij(t)为归一化后的观测数据,Uij(t)为归一化后的模拟数据,bij(t)为步骤5)所求振幅曲线。
本发明所设计的工作流程和目标函数有以下优点:
(1)当地震观测数据由于震源或检波器的耦合问题或地层介质吸收等原因与正演模拟数据相比振幅不匹配的情况下,使反演收敛于合理结果。
(2)使地震波形反演可以只反演速度模型而忽略密度的变化,因为在地下介质的质面通常会同时存在速度和密度的变化,而密度的变化主要影响地震观测信号的幅度。
附图说明
图1是当振幅匹配情况下应用本发明的一个单道地震时间记录实例图,a)观测地震道d;b)模拟地震道u;c)直接相减得到的残差;d)应用本发明做振幅匹配后相减的残差;
在本例中,由于d和u的振幅是匹配的,因此振幅匹配后的残差的结果和直接相减的残差的结果相同。
图2是当振幅不匹配情况下应用本发明的另一个单道地震时间记录实例图,a)观测地震道d,振幅沿时间方向被衰减;b)模拟地震道u;c)基于算术减法得到的残差;d)振幅匹配后的残差;
本例中的观测数据d延时间方向的衰减较快,其振幅无法与模拟数据u匹配。直接相减导致了错误的残差,而应用本发明做振幅匹配后的残差则在“波形”上更为正确(与上例相同)。
图3是观测数据图;
图4是模拟数据图;
图5是观测数据归一化结果图;
图6是模拟数据归一化结果图;
图7是归一化残差结果图;
图8是归一化残差振幅恢复结果图。
具体实施方式
以下结合附图和实例详细说明本发明。
本发明通过以下步骤实现:
1)野外激发人工地震波,由检波器接收并记录地震波场,得到地震观测数据dij(t)=d(xi,xj,t),单道示例如图3观测数据;
其中:xi和xj,分别代表震源和接收点的空间位置,t代表观测时间;
2)基于初始速度模型V,用有限差分算法模拟相应的地震观测数据,得到地震模拟数据uij(t)=u(xi,xj,t),单道示例如图4模拟数据;
其中:xi和xj,分别代表震源和接收点的空间位置,t代表观测时间;
3)基于反演使用频段,应用以下公式求取时间窗口长度N,其中Δt是时间采样率,fd是当前反演频段的主频;
4)按前文所述均方根振幅曲线公式或高斯窗振幅曲线公式对观测数据dij(t)求取时间窗口长度N的振幅曲线aij(t);
5)按前文所述均方根振幅曲线公式或高斯窗振幅曲线公式对模拟数据uij(t)求取时间窗口长度N的振幅曲线bij(t);
6)对观测数据作归一化处理:Dij(t)=dij(t)/aij(t),结果如图5观测数据归一化结果;
其中:dij(t)为观测数据,aij(t)为步骤4)所求振幅曲线;
7)对模拟数据作归一化处理:Uij(t)=uij(t)/bij(t),结果如图6模拟数据归一化结果;
其中:uij(t)为模拟数据,bij(t)为步骤5)所求振幅曲线;
8)设立目标函数,通过高斯-牛顿优化方法使得野外记录的地震波场和正演模拟的地震波场误差减小,多次迭代反演使目标函数最小化,达到时间域地震波形反演收敛。
步骤8)所述的目标函数是:
其中:Dij(t)为归一化后的观测数据,Uij(t)为归一化后的模拟数据,bij(t)为步骤5)所求振幅曲线;
Dij(t)-Uij(t)表示观测数据Dij(t)与模拟数据Uij(t)归一化残差结果,如图7归一化残差结果;
[Dij(t)-Uij(t)]*bij(t)表示归一化残差振幅恢复结果,如图8归一化残差振幅恢复结果。
Claims (4)
1.一种在振幅不匹配时使时间域地震波形反演收敛的方法,特点是通过以下步骤实现:
1)野外激发人工地震波,由检波器接收并记录地震波场,得到地震观测数据dij(t)=d(xi,xj,t),
其中:xi和xj,分别代表震源和接收点的空间位置,t代表观测时间;
2)基于初始速度模型V,用有限差分算法模拟相应的地震观测数据,得到地震模拟数据uij(t)=u(xi,xj,t),
其中:xi和xj,分别代表震源和接收点的空间位置,t代表观测时间;
3)基于反演使用频段,应用以下公式求取时间窗口长度N,
其中Δt是时间采样率,fd是当前反演频段的主频;
4)对观测数据dij(t)求取时间窗口长度N的振幅曲线aij(t);
5)对模拟数据uij(t)求取时间窗口长度N的振幅曲线bij(t);
6)对观测数据作归一化处理:Dij(t)=dij(t)/aij(t);
其中:dij(t)为观测数据,aij(t)为步骤4)所求振幅曲线;
7)对模拟数据作归一化处理:Uij(t)=uij(t)/bij(t);
其中:uij(t)为模拟数据,bij(t)为步骤5)所求振幅曲线;
8)设立目标函数,通过高斯-牛顿优化方法使得野外记录的地震波场和正演模拟的地震波场误差减小,多次迭代反演使目标函数最小化,达到时间域地震波形反演收敛;
所述的步骤8)所述的目标函数是:
其中:Dij(t)为归一化后的观测数据,Uij(t)为归一化后的模拟数据,bij(t)为步骤5)所求振幅曲线。
2.根据权利要求1的方法,特点是所述的步骤4)和步骤5)求取振幅曲线采用均方根振幅曲线公式或高斯窗振幅曲线公式。
3.根据权利要求2的方法,特点是所述的均方根振幅曲线公式是:
其中:N为时窗大小,Δt为时间采样率,n为求和公式中的变量取值从到
4.根据权利要求2的方法,特点是所述的高斯窗振幅曲线公式是:
其中:N为时窗大小,Δt为时间采样率,σ是一个小于等于0.5的常数,n为求和公式中的变量取值从到
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410321860.5A CN104123440B (zh) | 2014-07-08 | 2014-07-08 | 一种在振幅不匹配时使时间域地震波形反演收敛的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410321860.5A CN104123440B (zh) | 2014-07-08 | 2014-07-08 | 一种在振幅不匹配时使时间域地震波形反演收敛的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104123440A CN104123440A (zh) | 2014-10-29 |
CN104123440B true CN104123440B (zh) | 2017-03-15 |
Family
ID=51768849
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410321860.5A Active CN104123440B (zh) | 2014-07-08 | 2014-07-08 | 一种在振幅不匹配时使时间域地震波形反演收敛的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104123440B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105929446B (zh) * | 2016-04-19 | 2018-06-01 | 中国石油天然气集团公司 | 一种全波形反演中的数据处理方法及装置 |
CN107462922B (zh) * | 2016-06-06 | 2019-08-27 | 中国石油化工股份有限公司 | 全波形反演方法及装置 |
-
2014
- 2014-07-08 CN CN201410321860.5A patent/CN104123440B/zh active Active
Non-Patent Citations (4)
Title |
---|
"基于探地雷达的混凝土无损检测反演成像方法";丁亮等;《地球物理学报》;20120131;第55卷(第1期);全文 * |
"基于极快速模拟退火算法的地层横波各向异性反演";王才志等;《石油学报》;20070131;第28卷(第1期);全文 * |
"煤田地震多属性分析技术的应用";石瑛等;《煤炭学报》;20081231;第33卷(第12期);全文 * |
"频率域叠前波动方程反演及其应用";李国平等;《石油地球物理勘探》;20110630;第46卷(第3期);第411-416页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104123440A (zh) | 2014-10-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Song et al. | Application of particle swarm optimization to interpret Rayleigh wave dispersion curves | |
Chen et al. | Detecting a known near-surface target through application of frequency-dependent traveltime tomography and full-waveform inversion to P-and SH-wave seismic refraction data | |
KR101861060B1 (ko) | 오프셋 의존 탄성 fwi를 통한 다중 파라미터 역산 | |
US20170335675A1 (en) | Method To Predict Pore Pressure And Seal Integrity Using Full Wavefield Inversion | |
CN103460074A (zh) | 全波场反演中小波估计和多次波预测的方法 | |
KR20110057124A (ko) | 지진 표면파들의 파형들을 사용하는 토양 특성들의 추정 | |
CN102937721A (zh) | 利用初至波走时的有限频层析成像方法 | |
CN103913768A (zh) | 基于地震波资料对地表中浅层进行建模的方法及装置 | |
CN106324662A (zh) | 一种针对目标层的全波形反演方法及系统 | |
Lee et al. | Analysis of nonlinear site response using the LSST downhole accelerometer array data | |
CN104123440B (zh) | 一种在振幅不匹配时使时间域地震波形反演收敛的方法 | |
Tran et al. | Inversion of first-arrival time using simulated annealing | |
CN104536046A (zh) | 基于地震记录的震源激发信号一致性评价方法 | |
CN117031539A (zh) | 一种自监督深度学习地震数据低频重建方法及系统 | |
CN107918152B (zh) | 一种地震相干层析成像方法 | |
Gao et al. | Acquisition and processing pitfall with clipped traces in surface-wave analysis | |
CN105445795A (zh) | 倾斜界面的速度分析方法及其系统 | |
CN113267809B (zh) | I类页岩储层预测方法及装置 | |
CN112415601A (zh) | 表层品质因子q值的确定方法及装置 | |
CN116626751B (zh) | 基于多目标函数的黏弹性参数同步反演方法、装置和设备 | |
CN108260359A (zh) | 一种地震信号处理方法、装置和系统 | |
KR101578894B1 (ko) | 지하 매질구조 추정장치 및 그 방법 | |
CN112526611A (zh) | 表层地震波品质因子的提取方法及装置 | |
Gao et al. | Multimode wave-equation-based Scholte wave inversion | |
CN103675898B (zh) | 一种快速精确估算地层吸收所致时差的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |