CN103852789B - 用于地震数据的非线性层析方法及其装置 - Google Patents

用于地震数据的非线性层析方法及其装置 Download PDF

Info

Publication number
CN103852789B
CN103852789B CN201410091200.2A CN201410091200A CN103852789B CN 103852789 B CN103852789 B CN 103852789B CN 201410091200 A CN201410091200 A CN 201410091200A CN 103852789 B CN103852789 B CN 103852789B
Authority
CN
China
Prior art keywords
data
seismic data
initial
parameters
processing
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
Application number
CN201410091200.2A
Other languages
English (en)
Other versions
CN103852789A (zh
Inventor
张恩嘉
金德刚
唐虎
吴战培
陈爱萍
陈三平
曹中林
陈波
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Petroleum Corp
BGP Inc
Original Assignee
Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd filed Critical Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd
Priority to CN201410091200.2A priority Critical patent/CN103852789B/zh
Publication of CN103852789A publication Critical patent/CN103852789A/zh
Application granted granted Critical
Publication of CN103852789B publication Critical patent/CN103852789B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供了一种用于地震数据的非线性层析方法和装置,所述方法包括:获取初始地震数据参数,其中,初始地震数据参数包括地表观测数据和初始速度模型;利用叠前深度偏移以及运动学反偏移对初始地震数据参数进行处理,得到处理后的第一数据,其中,第一数据包括目标函数值;通过第一数据得到第二数据,其中,第二数据包括迭代下降方向和迭代步长;当目标函数值下降时,输出与迭代下降方向相对应的速度更新模型。通过采用非线性层析的方法对速度模型进行更新,很大程度上提高了计算效率,且得到的速度更新模型更加稳定可靠。

Description

用于地震数据的非线性层析方法及其装置
技术领域
本发明涉及地震数据,更具体地讲,涉及一种用于地震数据的非线性层析方法和装置。
背景技术
随着地球物理勘探技术的广泛应用,出现了很对针对用于地震数据层析方法。现有的层析方法主要有走时层析和线性层析。走时层析利用走时的残差反投影到射线路径上进行模型更新,其对初始模型的要求相对较低,但反演速度的分辨率较低。线性层析在每次更新速度模型后,都需要重新进行偏移和拾取,才能得到更新后的模型的道集是否拉平以及剩余动校正量,理论上只有在局部线性化的假设下才能一次收敛到真解,但实际上需要反复多次更新才能反演出真实的速度,而在此过程中重复的偏移和拾取会大大降低层析的效率。
因此,需要一种提高更新效率且降低计算成本的层析方法。
发明内容
本发明的目的在于解决上述技术问题,因此,本发明的一方面提供了一种用于地震数据的非线性层析方法和装置。
为了实现上述目的,提供了一种用于地震数据的非线性层析方法,可包括以下步骤:获取初始地震数据参数,其中,初始地震数据参数包括地表观测数据和初始速度模型;利用叠前深度偏移以及运动学反偏移对初始地震数据参数进行处理,得到处理后的第一数据,其中,第一数据包括目标函数值;通过第一数据得到第二数据,其中,第二数据包括迭代下降方向和迭代步长;当目标函数值下降时,输出与迭代下降方向相对应的速度更新模型。
优选地,第一数据可还包括:剩余时差和Frechet微商,通过第一数据得到第二数据的步骤可包括:通过剩余时差和Frechet微商得到梯度参数,通过梯度参数得到第二数据。
优选地,利用叠前深度偏移以及运动学反偏移对初始地震数据参数进行处理得到处理后的第一数据的步骤可包括:初始地震数据参数经过道集拾取处理后得到地震数据运动学信息,地震数据运动学信息经过零射线偏移后得到反射点位置,再经过射线正演得到第一数据。
优选地,初始地震数据参数可经过偏移得到角度道集数据,角度道集数据通过道集拾取处理和反偏移处理后得到地震数据运动学信息。
为了实现上述目的,提出了一种用于地震数据的非线性层析装置,可包括:获取模块,获取初始地震数据参数,其中,初始地震数据参数包括地表观测数据和初始速度模型;第一处理模块,利用叠前深度偏移以及运动学反偏移对初始地震数据参数进行处理,得到处理后的第一数据,其中,第一数据包括目标函数值;第二处理模块,通过第一数据得到第二数据,其中,第二数据包括迭代下降方向和迭代步长;输出模块,当目标函数值下降时,输出与迭代下降方向相对应的速度更新模型。
优选地,第一数据可还包括:剩余时差和Frechet微商,第二处理模块可通过剩余时差和Frechet微商得到梯度参数,通过梯度参数得到第二数据。
优选地,第一处理模块可将初始地震数据参数经过道集拾取处理后得到地震数据运动学信息,地震数据运动学信息经过零射线偏移后得到反射点位置,再经过射线正演得到第一数据。
优选地,第一处理模块可将初始地震数据参数经过偏移得到角度道集数据,角度道集数据通过道集拾取处理和反偏移处理后得到地震数据运动学信息。
上述用于地震数据的非线性层析方法和装置,通过采用非线性层析的方法对速度模型进行更新,很大程度上提高了计算效率,且得到的速度更新模型更加稳定可靠。
附图说明
通过下面结合附图对本发明的示例性实施例进行的描述,本发明的上述和其他目的和特点将会变得更加清楚,其中:
图1是根据本发明第一示例性实施例的用于地震数据的非线性层析方法的流程图;
图2是根据本发明第二示例性实施例的用于地震数据的非线性层析方法的流程图;
图3是根据本发明第三示例性实施例的用于地震数据的非线性层析装置的框图。
具体实施方式
以下,参照附图来详细说明本发明的实施例。
总体来讲,根据本发明的用于地震数据的非线性层析方法包括以下步骤:获取初始地震数据参数,其中,初始地震数据参数包括地表观测数据和初始速度模型;利用叠前深度偏移以及运动学反偏移对初始地震数据参数进行处理,得到处理后的第一数据,其中,第一数据包括目标函数值;通过第一数据得到第二数据,其中,第二数据包括迭代下降方向和迭代步长;当目标函数值下降时,输出与迭代下降方向相对应的速度更新模型。
下面将结合本发明的具体实施例来对本发明的用于地震数据的非线性层析方法进行详细描述。
图1是根据本发明第一示例性实施例的用于地震数据的非线性层析方法的流程图。
在步骤S101,获取初始地震数据参数,其中,初始地震数据参数包括地表观测数据和初始速度模型。
在步骤S102,利用叠前深度偏移以及运动学反偏移对初始地震数据参数进行处理,得到处理后的第一数据,其中,第一数据包括目标函数值。
例如,初始地震数据参数可经过道集拾取处理后得到地震数据运动学信息,地震数据运动学信息经过零射线偏移后得到反射点位置,再经过射线正演得到第一数据。
例如,初始地震数据参数可经过偏移得到角度道集数据,角度道集数据通过道集拾取处理和反偏移处理后得到地震数据运动学信息。
在步骤S103,通过第一数据得到第二数据,其中,第二数据包括迭代下降方向和迭代步长。
此外,第一数据可还包括:剩余时差和Frechet微商(即射线路径)。可通过剩余时差和Frechet微商得到梯度参数,通过梯度参数得到第二数据,从而可通过第一数据来得到第二数据。
在步骤S104,当目标函数值下降时,输出与迭代下降方向相对应的速度更新模型。
图2是根据本发明第二示例性实施例的用于地震数据的非线性层析方法的流程图。
在步骤S201,获取初始地震数据参数中的地表观测数据dobs和初始速度模型m0。此时迭代次数k为1。
在步骤S202,利用叠前深度偏移以及运动学反偏移计算剩余时差Frechet微商以及目标函数值(Ck),设定Q=0,Q为目标函数值连续上升的累计次数。
具体地讲,地表观测数据dobs经过偏移后得到角道集参数,角道集参数通过道集拾取后得到拾取后的角道集,拾取后的角道集再经过反偏移得到地震数据运动学信息,地震数据运动学信息经过零偏射线偏移后确定反射点的位置,通过射线正演得到地表正演数据d和地震数据叠前走时tobs,通过地表正演数据d得到正演的射线走时t。
对于网格表达的模型,Frechet微商矩阵的元素可由下面的公式(1):
∂ t i ∂ m j = l ij - - - ( 1 )
其中,lij为第i根射线、第j个网格内的射线路径的片段长度,本发明采用样条表达的模型,通过下面公式(2)表示:
∂ t i ∂ m j = ∫ l β x ( x ) β y ( y ) β z ( z ) dl - - - ( 2 )
其中,βi(i=x,y,z)为三个方向上的样条基函数(样条基函数为已知量)。
由公式(1)和(2)即可得到第k次迭代中模型m的Frechet微商的值。通过剩余时差和Frechet微商得到第k次迭代的目标函数值的公式(如公式(3)所示)。
C ( m ) = 1 2 ( | | t ( m ) - t obs | | 2 + λ | | Dm | | 2 ) - - - ( 3 )
其中,C(m)为第k次迭代中模型m的目标函数值,即Ck,t(m)为模型m的正演的射线走时,tobs为地震数据叠前走时(tobs为不变量),λ为正则化因子(其为已知量),Dm为模型m的差分矩阵(其为已知量)。
在步骤S203,判断目标函数值是否下降,即比较Ck是否小于Ck-1,如果下降,则执行步骤S204,计算梯度参数g(k),并通过梯度参数计算得到迭代的下降方向p(k)和迭代步长α(k)。随后,在步骤S205,输出更新的速度模型m。具体公式如下面公式(4)所示。
▿ C ( m ) = ( ∂ t ∂ m ) T ( t - t obs ) | m = m 0 + λD T Dm | m = m 0 - - - ( 4 )
其中,为第k次迭代中模型m的梯度参数的值,即g(k)
通过采用共轭梯度法进行非线性优化,得到第k次迭代的迭代下降方向p(k)和迭代步长α(k),如公式(5)和(6)所示。
p(k)=g(k)(k)p(k-1) (5)
其中, β ( k ) = | | g ( k ) | | 2 2 | | g ( k - 1 ) | | 2 2 .
α ( k ) = C k | | g ( k ) | | 2 2 - - - ( 6 )
如果没有下降,则执行步骤S206,Q值增加1,然后执行步骤S207,判断Q值是否大于预定值Q0(一般情况下,预定值大于或等于三次),如果大于预定值,则执行步骤S208,停止迭代,即结束此次更新。
如果小于等于预定值,则执行步骤S202,重新计算剩余时差、Frechet微商以及目标函数值。具体地讲,利用已经得到的地震数据运动学信息,在更新后速度模型中进行零偏射线偏移确定反射点的位置,然后通过射线正演得到更新的地表正演数据d,将更新的地表正演数据与地震数据叠前走时tobs进行匹配,通过更新的地表正演数据d得到更新的正演的射线走时t。由于采用了更新后的速度拾取对数据进行一系列处理,使得通过公式(1)、(2)和(3)得到的剩余时差、Frechet微商以及目标函数值也是更新后的值。
图3是根据本发明第三示例性实施例的用于地震数据的非线性层析装置的框图。
如图3所示,本发明示例性实施例的用于地震数据的非线性层析装置包括:获取模块100、第一处理模块200、第二处理模块300和输出模块400。
获取模块100用于获取初始地震数据参数,其中,初始地震数据参数包括地表观测数据和初始速度模型。
第一处理模块200利用叠前深度偏移以及运动学反偏移对初始地震数据参数进行处理,得到处理后的第一数据,其中,第一数据包括目标函数值。
第二处理模块300通过第一数据得到第二数据,其中,第二数据包括迭代下降方向和迭代步长。
当目标函数值下降时,输出模块400输出与迭代下降方向相对应的速度更新模型。
进一步地,第一数据还包括:剩余时差和Frechet微商;第二处理模块通过剩余时差和Frechet微商得到梯度参数,通过梯度参数得到第二数据。
可选地,第一处理模块将初始地震数据参数经过道集拾取处理后得到地震数据运动学信息,地震数据运动学信息经过零射线偏移后得到反射点位置,再经过射线正演得到第一数据。
可选地,第一处理模块将初始地震数据参数经过偏移得到角度道集数据,角度道集数据通过道集拾取处理和反偏移处理后得到地震数据运动学信息。
应该理解,根据本发明示例性实施例的用于地震数据的非线性层析装置可执行以上参照图1至图2描述的用于地震数据的非线性层析方法,为了避免重复,在此不再赘述。
上述用于地震数据的非线性层析方法和装置,通过采用非线性层析的方法对速度模型进行更新,很大程度上提高了计算效率,且得到的速度更新模型更加稳定可靠。
上面已经结合具体实施例描述了本发明,但是本发明的实施不限于此。在本发明的精神和范围内,本领域技术人员可以进行各种修改和变型,这些修改和变型将落入权利要求限定的保护范围之内。

Claims (4)

1.一种用于地震数据的非线性层析方法,包括以下步骤:
获取初始地震数据参数,其中,初始地震数据参数包括地表观测数据和初始速度模型;
利用叠前深度偏移以及运动学反偏移对初始地震数据参数进行处理,得到处理后的第一数据,其中,第一数据包括目标函数值;
通过第一数据得到第二数据,其中,第二数据包括迭代下降方向和迭代步长;
当目标函数值下降时,输出与迭代下降方向相对应的速度更新模型,
其中,利用叠前深度偏移以及运动学反偏移对初始地震数据参数进行处理得到处理后的第一数据的步骤包括:初始地震数据参数经过道集拾取处理后得到地震数据运动学信息,地震数据运动学信息经过零射线偏移后得到反射点位置,再经过射线正演得到第一数据,
其中,初始地震数据参数经过偏移得到角度道集数据,角度道集数据通过道集拾取处理和反偏移处理后得到地震数据运动学信息。
2.如权利要求1所述的方法,其特征在于,第一数据还包括:剩余时差和Frechet微商,
通过第一数据得到第二数据的步骤包括:通过剩余时差和Frechet微商得到梯度参数,通过梯度参数得到第二数据。
3.一种用于地震数据的非线性层析装置,包括:
获取模块,获取初始地震数据参数,其中,初始地震数据参数包括地表观测数据和初始速度模型;
第一处理模块,利用叠前深度偏移以及运动学反偏移对初始地震数据参数进行处理,得到处理后的第一数据,其中,第一数据包括目标函数值;
第二处理模块,通过第一数据得到第二数据,其中,第二数据包括迭代下降方向和迭代步长;
输出模块,当目标函数值下降时,输出与迭代下降方向相对应的速度更新模型,
其中,第一处理模块将初始地震数据参数经过道集拾取处理后得到地震数据运动学信息,地震数据运动学信息经过零射线偏移后得到反射点位置,再经过射线正演得到第一数据,
其中,第一处理模块将初始地震数据参数经过偏移得到角度道集数据,角度道集数据通过道集拾取处理和反偏移处理后得到地震数据运动学信息。
4.如权利要求3所述的装置,其特征在于,第一数据还包括:剩余时差和Frechet微商,第二处理模块通过剩余时差和Frechet微商得到梯度参数,通过梯度参数得到第二数据。
CN201410091200.2A 2014-03-12 2014-03-12 用于地震数据的非线性层析方法及其装置 Active CN103852789B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410091200.2A CN103852789B (zh) 2014-03-12 2014-03-12 用于地震数据的非线性层析方法及其装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410091200.2A CN103852789B (zh) 2014-03-12 2014-03-12 用于地震数据的非线性层析方法及其装置

Publications (2)

Publication Number Publication Date
CN103852789A CN103852789A (zh) 2014-06-11
CN103852789B true CN103852789B (zh) 2017-02-01

Family

ID=50860668

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410091200.2A Active CN103852789B (zh) 2014-03-12 2014-03-12 用于地震数据的非线性层析方法及其装置

Country Status (1)

Country Link
CN (1) CN103852789B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105319598A (zh) * 2014-07-31 2016-02-10 中国石油化工股份有限公司 一种不需要拾取的叠前地震反射层析方法
CN104391323B (zh) * 2014-11-21 2015-11-18 中国石油大学(华东) 一种利用反射波信息反演速度场中低波数成分的方法
CN104360396B (zh) * 2014-12-04 2016-10-19 中国海洋石油总公司 一种海上井间tti介质三种初至波走时层析成像方法
CN107203002B (zh) * 2017-06-12 2019-05-24 中国科学院地质与地球物理研究所 反演速度模型的建立方法和地下结构的像的获得方法
CN111239822B (zh) * 2020-02-27 2022-08-23 广州海洋地质调查局 一种消除井震闭合差的断层下速度建模方法及处理终端

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102096099A (zh) * 2009-12-14 2011-06-15 武侯区巅峰机电科技研发中心 折射波射线层析成像方法
CN103217715B (zh) * 2012-01-18 2015-12-02 中国石油天然气集团公司 多尺度规则网格层析反演静校正方法
CN102841375A (zh) * 2012-09-06 2012-12-26 中国石油大学(华东) 一种复杂条件下基于角度域共成像点道集的层析速度反演方法
CN102841376A (zh) * 2012-09-06 2012-12-26 中国石油大学(华东) 一种基于起伏地表的层析速度反演方法
CN102937721B (zh) * 2012-11-07 2015-07-08 中国石油集团川庆钻探工程有限公司地球物理勘探公司 利用初至波走时的有限频层析成像方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Herve Chauris et al..Migration velocity analysis from locally coherent events in 2-D laterally heterogeneous media, Part I: Theoretical aspects.《GEOPHYSICS》.2002,第67卷(第4期), *
Nonlinear Diffraction Tomography Applied to Cross-Borehole Data;Delaine R. Thompson et al.;《1992 SEG Annual Meeting》;19921231;第775页第3段、第7段、第9段、第11段 *
倪瑶等.立体层析反演方法理论分析与应用测试.《石油物探》.2013,第52卷(第2期), *
段鹏飞等.TI介质局部角度域射线追踪与叠前深度偏移成像.《地球物理学报》.2013,第56卷(第1期), *
特征波立体层析速度反演方法及应用研究;孙小东等;《中国地球物理2013》;20131231;第895页 *

Also Published As

Publication number Publication date
CN103852789A (zh) 2014-06-11

Similar Documents

Publication Publication Date Title
CN103852789B (zh) 用于地震数据的非线性层析方法及其装置
CN104536043B (zh) 一种深度域整体速度模型融合方法及装置
CN106483559B (zh) 一种地下速度模型的构建方法
CA2961572C (en) Velocity tomography using property scans
CN105319589B (zh) 一种利用局部同相轴斜率的全自动立体层析反演方法
CN109917454B (zh) 基于双基准面的真地表叠前深度偏移成像方法及装置
CN107783186B (zh) 一种近地表速度建模方法及装置
US20160291177A1 (en) System and method for saltsurface updating via wavefield redatuming
CN109633749B (zh) 基于散射积分法的非线性菲涅尔体地震走时层析成像方法
CN106959467B (zh) 地震波速度反演方法和装置
CN111580163B (zh) 一种基于非单调搜索技术的全波形反演方法及系统
WO2014209779A1 (en) Processing survey data containing ghost data
CN104181593B (zh) 一种三维无射线追踪回折波层析成像方法及装置
CN112363211A (zh) 一种改进的sirt法射线走时层析成像方法
CN104749625B (zh) 一种基于正则化技术的地震数据倾角估计方法及装置
CN109655890A (zh) 一种深度域浅中深层联合层析反演速度建模方法及系统
CN105717538B (zh) 起伏地表地震数据偏移基准面转换方法及装置
CN113970789B (zh) 全波形反演方法、装置、存储介质及电子设备
CN115421195A (zh) 地震探测中速度场的生成方法、装置、设备及存储介质
CN113791447B (zh) 一种反射结构导引的反射波层析反演方法
CN112308206B (zh) 基于深度学习的大地电磁-地震波初至时联合反演方法
CN104570100B (zh) 多子波克希霍夫地震数据偏移方法
Kotsi et al. A statistical comparison of three 4d full-waveform inversion schemes
CN112465983A (zh) 一种高程地表模型点线状障碍物修正方法和装置
Przebindowska Acoustic full waveform inversion of marine reflection seismic data

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20180207

Address after: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189

Patentee after: BGP INC., CHINA NATIONAL PETROLEUM Corp.

Address before: Shuangliu County Huayang Huayang Road in Chengdu city of Sichuan Province in 610213 section of No. 216, Igawa geophysical exploration company of the Ministry of science and technology

Patentee before: GEOPHYSICAL EXPLORATION COMPANY OF CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20200921

Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No.

Patentee after: CHINA NATIONAL PETROLEUM Corp.

Patentee after: BGP Inc., China National Petroleum Corp.

Address before: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189

Patentee before: BGP Inc., China National Petroleum Corp.