CN104977608B - 利用固定网格声波波场模拟的时间域全波形反演方法 - Google Patents
利用固定网格声波波场模拟的时间域全波形反演方法 Download PDFInfo
- Publication number
- CN104977608B CN104977608B CN201410140367.3A CN201410140367A CN104977608B CN 104977608 B CN104977608 B CN 104977608B CN 201410140367 A CN201410140367 A CN 201410140367A CN 104977608 B CN104977608 B CN 104977608B
- Authority
- CN
- China
- Prior art keywords
- wave
- time
- speed
- field
- wave field
- 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
本发明是用固定网格声波波场模拟的时间域全波形反演方法,利用规则的固定网格进行剖分建立初始速度场,在时间‑波数域得到精确的波场随时间演化关系式,根据空间点的波场传播速度和时间延拓步长,用最小二乘寻优的方法确定星形网格各点的褶积系数,使星形网格各点对应的波数域时间演化关系与精确的波场随时间演化关系式中演化关系间的误差最小化,用最小二乘寻优的方法,分别求取对应速度的褶积算子,利用褶积算子和速度索引进行波场正演模拟并记录模拟的地震波场,使得野外记录的地震波场和正演模拟的地震波场误差减小,迭代更新速度模型完成反演。本发明解决了现有方法时间导数近似计算精度不足问题,运算量小,效率高。
Description
技术领域
本发明涉及石油勘探技术,是一种应用于地震数据深度建模中的利用固定网格声波波场模拟的时间域地震波形反演方法。
背景技术
地震数据处理的重要目的之一是通过深度偏移(Depth Migration)得到高质量的深度域图像。深度域图像的质量依赖于地震速度模型的质量,因此建模技术在整个地震数据处理起着举足轻重的作用。全波形反演(Full Waveform Inversion)是一种数据域(DataDomain)的地震建模技术,其基本原理是通过比较正演得到的模拟数据和实际观测数据的差别来确定模型改变量,反复迭代此过程直到模拟数据和观测数据的差别达到最小。和传统的基于图像域(Image Domain)的建模方法(如层析反演)相比,全波形反演方法的优点在于(1)分辨率高;(2)需要的人工干预少。
目前业内全波形反演方法的实现方法主要有两个方向:(1)在时间域实现;(2)在频率域实现。其中时间域算法中,高精度高效的波场模拟算法是其关键因素。在众多模拟算法中,有限差分方法是应用最为普遍的一个。在有限差分算法中,受计算机内存和计算量的限制,尽管空间导数会使用高阶逼近,但时间导数一般都使用二阶逼近,由于时间方向引入的误差偏大,会导致所谓的数字频散现象。这种误差在正演模拟的结果中会表现为走时误差和相位旋转。为提高时间方向的模拟精度,有学者提出了伪谱法(Tal-Ezer et al,1987,Etgen,1989,Zhang and Zhang 2009),但这些方法要在毎个时间步中多次运用髙维付氏变换,比有限差分方法运算量大很多。
发明内容
本发明的目的是提供一种在解决现有时间域全波形反演实现中,由于时间方向模拟精度不足导致的数值频散,影响反演精度的问题的利用固定网格声波波场模拟的时间域全波形反演方法。
本发明实现步骤如下:
1)在工区沿炮线方向人工激发地震波,由检波器记录地震波波场;
2)对工区利用规则的固定网格进行剖分,根据工区地质信息等建立初始速度场;
3)以速度为常速,在时间-波数域求出声波方程的解析解,得到精确的波场随时间演化关系式:
式中:
为波场,v为当地速度,Δt为时间步长,k为波数,t为时间;
4)根据空间点的波场传播速度和时间延拓步长,用最小二乘寻优的方法确定星形网格各点的褶积系数,使星形网格各点对应的波数域时间演化关系与步骤3)中精确的波场随时间演化关系式中演化关系间的误差最小化;
5)在模拟空间范围内,分析速度变化范围,并建立速度索引,利用步骤4)中最小二乘寻优的方法,分别求取对应速度的褶积算子;
6)在时间空间域,利用得到的褶积算子和速度索引进行波场正演模拟并记录模拟的地震波场;
7)利用高斯-牛顿优化方法使得野外记录的地震波场和正演模拟的地震波场误差减小,迭代更新速度模型完成反演。
本发明相对于现有技术具有以下优点和技术效果:
解决了现有方法时间导数近似计算精度不足问题。利用声波方程在时间波数域的表达形式,得到精确的波场时间演化关系。可以证明波数域中二阶时间差分的有限差分格式是这个时间演化关系的二阶截断。因此时间演化关系提高了时间导数近似计算的精度。解决了高维fourier变换运算量太大的问题。在空间域求取褶积算子,利用最小二乘法,使之在波数域的有效波数带内拟合时间演化算子,避免多次运用高维fourier变换,减少运算量。解决了波场模拟运算量大计算效率低的问题。尽管在每一个时间步中,时间演化算子的计算量比相应阶数的有限差分算子略大,但要达到相同的波场模拟精度,有限差分法必须釆用小于数倍的时间步长。因此时间演化算子在精度上和效率上都有很强的优势。
附图说明
图1是规则网格显格式时间演化算子的星形格点分布图;
图2是时间演化算法与差分方法分别在近场和远场的单道模拟波形对比图;
图3是时间演化算法与差分方法远道模拟效率对比图。
具体实施方式
以下结合附图和实例详细描述本发明。
本发明所要解决的技术问题是一种应用于全波场反演的高精度高效率声波波场模拟技术。具体的实施方案如下:
(1)在工区沿炮线方向人工激发地震波,由检波器接收并记录地震波场,得到地震观测数据,其中xi和xj,分别代表震源和接收点的空间位置,t代表观测时间。
(2)对工区利用规则的固定网格进行剖分,根据工区地质信息等建立初始速度场V。图1是规则网格显格式时间演化算子的星形格点分布图。
(3)以速度为常速,在时间-波数域求出声波方程的解析解,得到精确的波场随时间演化关系式:
式中:
为波场,v为当地速度,Δt为时间步长,k为波数,t为时间;
基于初始速度场V,应用本发明的精确的波场随时间演化关系式进行数字正演,模拟地震观测数据,得到地震模拟数据;
(4)根据空间点的波场传播速度和时间延拓步长,用最小二乘寻优的方法确定星形网格各点的褶积系数,使星形网格各点对应的波数域时间演化关系与步骤(3)中精确的波场随时间演化关系式中演化关系间的误差最小化;
图2是时间演化算法与差分方法分别在近场和远场的单道模拟波形对比图;在近场位置,时间演化算法与差分方法模拟波形与参考值接近,误差大致相等;随时间推移,远场位置时间演化算法的波形与参考值吻合很好,而差分方法发生频散。
图3是时间演化算法与差分方法远道模拟效率对比图;如要达到相同的传播精度,差分算法的时间传播步长需要减少至时间演化算法的1/5。
(5)在模拟空间范围内,分析速度变化范围,并建立速度索引,利用步骤(4)中最小二乘寻优的方法,分别求取对应速度的褶积算子;
(6)在时间空间域,利用得到的褶积算子和速度索引进行波场正演模拟并记录模拟的地震波场,对观测数据和模拟数据求取残差rij(t)=dij(t)-uij(t),利用残差反演校正方向g(xi,xj,xk),反演过程中应用本发明的波场模拟方法进行炮点正向波场模拟和残差反向传播的波场模拟;
(7)利用求取的校正方向和合理的步长,求取速度模型修正量,对初始模型进行修改,利用高斯-牛顿优化方法使得野外记录的地震波场和正演模拟的地震波场误差rij(t)减小,设立目标函数,多次迭代反演过程,使目标函数最小化,改进速度模型的精确度。
综上所述,本发明提供了一种时间域高效率高精度的波场模拟方法,在时间域全波形反演技术的实现过程中,可以用此方法代替原来的有限差分波场模拟方法,提高波场模拟的精度,有利于波形的匹配。
Claims (1)
1.一利用固定网格声波波场模拟的时间域全波形反演方法,特点是采用以下步骤实现:
1)在工区沿炮线方向人工激发地震波,由检波器记录地震波波场;
2)对工区利用规则的固定网格进行剖分,根据工区地质信息建立初始速度场;
3)以速度为常速,在时间-波数域求出声波方程的解析解,得到精确的波场随时间演化关系式:
式中:
为波场,v为当地速度,Δt为时间步长,k为波数,t为时间;
4)根据空间点的波场传播速度和时间延拓步长,用最小二乘寻优的方法确定星形网格各点的褶积系数,使星形网格各点对应的波数域时间演化关系与步骤3)中精确的波场随时间演化关系式中演化关系间的误差最小化;
5)在模拟空间范围内,分析速度变化范围,并建立速度索引,利用步骤4)中最小二乘寻优的方法,分别求取对应速度的褶积系数;
6)在时间空间域,利用得到的褶积系数和速度索引进行波场正演模拟并记录模拟的地震波场;
7)利用高斯-牛顿优化方法使得野外记录的地震波场和正演模拟的地震波场误差减小,迭代更新速度模型完成反演。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410140367.3A CN104977608B (zh) | 2014-04-09 | 2014-04-09 | 利用固定网格声波波场模拟的时间域全波形反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410140367.3A CN104977608B (zh) | 2014-04-09 | 2014-04-09 | 利用固定网格声波波场模拟的时间域全波形反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104977608A CN104977608A (zh) | 2015-10-14 |
CN104977608B true CN104977608B (zh) | 2017-07-07 |
Family
ID=54274283
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410140367.3A Active CN104977608B (zh) | 2014-04-09 | 2014-04-09 | 利用固定网格声波波场模拟的时间域全波形反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104977608B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102183790A (zh) * | 2011-02-12 | 2011-09-14 | 中国石油大学(华东) | 基于时空双变网格的弹性波正演模拟技术 |
CN103091711A (zh) * | 2013-01-24 | 2013-05-08 | 中国石油天然气集团公司 | 全波形反演方法及装置 |
CN103119552A (zh) * | 2010-09-27 | 2013-05-22 | 埃克森美孚上游研究公司 | 同时源编码和源分离作为全波场反演的实际解决方案 |
CN103207409A (zh) * | 2013-04-17 | 2013-07-17 | 中国海洋石油总公司 | 一种频率域全波形反演地震速度建模方法 |
CN103499835A (zh) * | 2013-10-13 | 2014-01-08 | 中国石油集团西北地质研究所 | 初至波波形反演近地表速度模型方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6999880B2 (en) * | 2003-03-18 | 2006-02-14 | The Regents Of The University Of California | Source-independent full waveform inversion of seismic data |
US9291734B2 (en) * | 2011-06-10 | 2016-03-22 | International Business Machines Corporation | Full waveform inversion using combined shot data and no scratch disk |
-
2014
- 2014-04-09 CN CN201410140367.3A patent/CN104977608B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103119552A (zh) * | 2010-09-27 | 2013-05-22 | 埃克森美孚上游研究公司 | 同时源编码和源分离作为全波场反演的实际解决方案 |
CN102183790A (zh) * | 2011-02-12 | 2011-09-14 | 中国石油大学(华东) | 基于时空双变网格的弹性波正演模拟技术 |
CN103091711A (zh) * | 2013-01-24 | 2013-05-08 | 中国石油天然气集团公司 | 全波形反演方法及装置 |
CN103207409A (zh) * | 2013-04-17 | 2013-07-17 | 中国海洋石油总公司 | 一种频率域全波形反演地震速度建模方法 |
CN103499835A (zh) * | 2013-10-13 | 2014-01-08 | 中国石油集团西北地质研究所 | 初至波波形反演近地表速度模型方法 |
Non-Patent Citations (7)
Title |
---|
"Fourier finite-difference wave propagation";Sergey Formel 等;《Geophysics》;20111031;第76卷(第5期);第123-129页 * |
"Multiscale seismic waveform inversion";Carey Bunks 等;《Geophysics》;19951031;第60卷(第5期);第1457-1473页 * |
"波动方程变网格步长有限差分数值模拟";孙成禹 等;《石油勘探》;20080331;第47卷(第2期);第123-127页 * |
"波动方程数值模拟的隐式差分法";朱生旺 等;《石油物探》;20060331;第45卷(第2期);第151-156页 * |
"用于声波方程数值模拟的时间-空间域有限差分系数确定新方法";梁文全 等;《地球物理学报》;20131024;第56卷(第10期);第3497-3506页 * |
"错格傅里叶伪谱微分算子在波场模拟中的应用";龙桂华 等;《地球物理学报》;20090131;第52卷(第1期);第193-199页 * |
"高维波动方程数值模拟的隐式***有限差分格式";张厚柱 等;《石油物探》;20071130;第46卷(第6期);第594-597页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104977608A (zh) | 2015-10-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104977607B (zh) | 利用变步长网格声波波场模拟的时间域全波形反演方法 | |
CN103238158B (zh) | 利用互相关目标函数进行的海洋拖缆数据同时源反演 | |
CN107843925B (zh) | 一种基于修正相位的反射波波形反演方法 | |
Rawlinson et al. | Seismic ray tracing and wavefront tracking in laterally heterogeneous media | |
CN106908835B (zh) | 带限格林函数滤波多尺度全波形反演方法 | |
CN103703391B (zh) | 使用频谱整形的全波场反演的***和计算机实施的方法 | |
CN103630933B (zh) | 基于非线性优化的时空域交错网格有限差分方法和装置 | |
CN107505654B (zh) | 基于地震记录积分的全波形反演方法 | |
CN105467444B (zh) | 一种弹性波全波形反演方法及装置 | |
CN109407151B (zh) | 基于波场局部相关时移的时间域全波形反演方法 | |
CN108549100B (zh) | 基于非线性高次拓频的时间域多尺度全波形反演方法 | |
KR20130060231A (ko) | 지구 물리학적 데이터의 반복 반전의 아티팩트 감소 | |
AU2015383134B2 (en) | Multistage full wavefield inversion process that generates a multiple free data set | |
Shragge | Solving the 3D acoustic wave equation on generalized structured meshes: A finite-difference time-domain approach | |
CN107894618B (zh) | 一种基于模型平滑算法的全波形反演梯度预处理方法 | |
CN109459789B (zh) | 基于振幅衰减与线性插值的时间域全波形反演方法 | |
CN111239806B (zh) | 基于振幅增量编码的时间域全波形反演方法 | |
CN112464520A (zh) | 局部重力异常深度反演方法及装置 | |
CN108680968B (zh) | 复杂构造区地震勘探数据采集观测***评价方法及装置 | |
CN109239776B (zh) | 一种地震波传播正演模拟方法和装置 | |
CN109164488B (zh) | 一种梯形网格有限差分地震波场模拟方法 | |
CN109490954B (zh) | 波场正演模拟方法及装置 | |
CN104977608B (zh) | 利用固定网格声波波场模拟的时间域全波形反演方法 | |
CN115270579A (zh) | 二阶声波方程有限差分数值模拟参数选取方法 | |
CN108680957A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |