CN104749625A - 一种基于正则化技术的地震数据倾角估计方法及装置 - Google Patents
一种基于正则化技术的地震数据倾角估计方法及装置 Download PDFInfo
- Publication number
- CN104749625A CN104749625A CN201510106959.8A CN201510106959A CN104749625A CN 104749625 A CN104749625 A CN 104749625A CN 201510106959 A CN201510106959 A CN 201510106959A CN 104749625 A CN104749625 A CN 104749625A
- Authority
- CN
- China
- Prior art keywords
- geological data
- sigma
- inclination angle
- regularization
- alpha
- 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.)
- Granted
Links
Abstract
本发明公开了一种基于正则化技术的地震数据倾角估计方法及装置,根据地震数据,选取其中任意一点作为地震数据倾角估计位置;依据全通数字滤波器计算出邻近地震道滤波器系数,进而计算出当前点地震数据预测值,其中包含待求地震数据倾角;循环执行上述未选取的地震数据位置,计算出所有位置地震数据预测值;根据所获取的地震数据、预测的地震数据及地震数据倾角构建正则化模型;利用非线性迭代反演算法获得地震数据倾角。解决现有技术中无法精确估计地震数据倾角难题。
Description
技术领域
本发明属于地震数据处理技术领域,涉及一种基于正则化技术的地震数据倾角估计方法,本发明还涉及一种基于正则化技术的地震数据倾角估计装置。
背景技术
估计地震数据或图像的倾角是地震数据处理重要研究内容,并广泛应用于地震波场分离与去噪(Harlan et al.,1984;Fomel et al.,2007)、反假频插值方法(Bardan,1987)、共反射面元叠加(CRS),法入射点(NIP)层析建模等。地震数据倾角属性不仅包含了重要的地震地层学信息,而且还可以直接用于工区构造解释。地震倾角属于运动学属性,其估计方法一般需要两个步骤。首先是检测步骤,主要是根据局部相干分析;其次是提取步骤,该过程主要依据相干水平与同相轴连续性。在地震数据倾角估计研究上,Ottolini提出了局部倾斜叠加方法;Barnes提出了复数道地震数据分析方法;Fehmers给出了一种利用局部构造张量估计倾角的方法;Marfurt等通过截取地震数据构造分析窗口计算地震波形相干值,并将最大相干值对应的倾角作为局部地层倾角;Fomel等人利用平面波破坏方法研究地层倾角,该方法通过局部平面波近似局部波场。但是现有的技术均无法精确估计地震数据倾角的问题。
考虑了地震数据通常被噪声污染及地层倾角平滑特性,本发明通过构建正则化模型,提出了一种快速精确的倾角估计方法及装置。
发明内容
本发明的目的是提供一种基于正则化技术的地震数据倾角估计方法,解决现有技术中无法精确估计地震数据倾角难题。
本发明的另一目的是提供一种基于正则化技术的地震数据倾角估计装置。
本发明所采用的技术方案是,一种基于正则化技术的地震数据倾角估计方法,包括以下步骤:
步骤101):根据地震数据,选取其中任意一点作为地震数据倾角估计位置;
步骤102):依据全通数字滤波器计算出邻近地震道滤波器系数,进而计算出当前点地震数据预测值,其中包含待求地震数据倾角;
步骤103):循环执行上述未选取的地震数据位置,计算出所有位置地震数据预测值;
步骤104):根据所获取的地震数据、预测的地震数据及地震数据倾角构建正则化模型;
步骤105):利用非线性迭代反演算法获得地震数据倾角。
本发明所采用的另一技术方案是,一种基于正则化技术的地震数据倾角估计装置,包括:
地震数据获取器件,用于从储存介质获取地震数据;
地震数据预测器件,用于依据全通数字滤波器计算出邻近地震道滤波器系数,进而计算出当前点地震数据预测值;
循环执行器件,用于循环执行所述未选取的地震数据位置,计算出所有位置地震数据预测值;
正则化模型构建器件,用于根据所述获取的地震数据、预测的地震数据及地震数据倾角构建正则化模型;
非线性求解器件,用于根据非线性迭代反演算法获得地震数据倾角。
本发明的有益效果是,本发明通过最小二乘化获取地震数据与预测地震数据残差,并约束地震数据倾角,构建了正则化模型,该模型充分考虑了地震倾角平滑特性,而且非线性迭代算法能够快速精确逼近真实值,适用于大规模地震数据计算。通过正则化模型,提供了一种精确的地震数据倾角估计方法,该技术能够更好的解决地震数据处理问题,如偏移成像孔径优选,速度建模,反射预测、噪声去除等,可以广泛应用于地震资料分析的各个环节。
附图说明
图1为本发明提出的地震数据倾角估计方法流程图。
图2为本发明的地震数据倾角估计装置结构框图。
图3为三维地震成像数据主测线方向示意图。
图4为三维地震成像数据地震数据倾角示意图。
图中,201.地震数据获取器件,202.地震数据预测器件,203.循环执行器件,204.正则化模型构建器件,205.非线性求解器件。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
一种基于正则化技术的地震数据倾角估计方法,流程如图1所示,包括以下步骤:
步骤101):根据地震数据,选取其中任意一点作为地震数据倾角估计位置;
步骤102):依据全通数字滤波器计算出邻近地震道滤波器系数,进而计算出当前点地震数据预测值,其中包含待求地震数据倾角;
步骤103):循环执行上述未选取的地震数据位置,计算出所有位置地震数据预测值;
步骤104):根据所获取的地震数据、预测的地震数据及地震数据倾角构建正则化模型;
步骤105):利用非线性迭代反演算法获得地震数据倾角。
其中,在上述步骤102中计算出当前点地震数据预测值包括以下步骤:
由全通数字滤波器设计出不同阶数的滤波器系数,并根据滤波器系数、地震数据倾角和获取的地震数据,计算出当前点地震数据预测值。
在上述步骤104中,构建正则化模型包括以下步骤:
根据所获取的地震数据和模拟数据预测值求最小二乘残差,作为模型构建第一项;约束的地震数据倾角作为第二项,构建正则化模型,如下所示:
其中,Jα(σ)为目标函数,d表示获取的地震数据,σ为待估计的地震数据倾角,P为全通数字滤波器破坏算子,α为正则化因子,L为拉普拉斯算子,min表示最小化。
在上述步骤105中,非线性迭代反演算法,通过计算梯度、方向参数实现迭代,包括以下步骤:
对上式Jα(σ)求梯度,得出:
其中,为梯度,且
且当仅且得出倾角σ临界点σ*:
其中,步长参数ξ>0,PΓ表示投影,由上式,给出以下固定点迭代公式:
迭代过程如下所示:
步骤1:输入初始倾角σ=0,设定迭代序号k:=0;
步骤2:k次迭代,并判断收敛条件;
步骤3:计算负梯度
步骤4:线性搜索得出,
步骤5:更新倾角:σk+1=PΓ(σk+ξksk);
步骤6:设置k:=k+1,返回步骤2.
其中,步长参数ξ,可通过Powell-Wolfe线性搜索准则获得。
一种基于正则化技术的地震数据倾角估计装置,结构如图2所示,
包括:
地震数据获取器件201,用于从储存介质获取地震数据;
地震数据预测器件202,用于依据全通数字滤波器计算出邻近地震道滤波器系数,进而计算出当前点地震数据预测值;
循环执行器件203,用于循环执行所述未选取的地震数据位置,计算出所有位置地震数据预测值;
正则化模型构建器件204,用于根据所述获取的地震数据、预测的地震数据及地震数据倾角构建正则化模型;
非线性求解器件205,用于根据非线性迭代反演算法获得地震数据倾角。
其中,地震数据预测器件202件包括:
滤波器系数计算单元,由全通数字滤波器设计出不用阶数的滤波器系数;
数值计算单元,根据滤波器系数、地震数据倾角和获取的地震数据,计算出当前点地震数据预测值。
正则化模型构建器件204包括:
最小二乘单元,根据所获取的地震数据和模拟预测值求最小二乘残差,作为模型构建第一项;
倾角单元,根据约束的地震数据倾角作为第二项,构建正则化估计模型,如下所示:
其中,Jα(σ)为目标函数,d表示获取的地震数据,σ为待估计的地震数据倾角,P为全通数字滤波器破坏算子,α为正则化因子,L为拉普拉斯算子,min表示最小化。
非线性求解器件205,通过计算梯度、方向等参数实现,迭代过程如下:
对上式Jα(σ)求梯度,得出:
其中,为梯度,且
且当仅且得出倾角σ临界点σ*:
其中,步长参数ξ>0,PΓ表示投影,由上式,给出以下固定点迭代公式:
迭代过程如下所示:
步骤1:输入初始倾角σ=0,设定迭代序号k:=0;
步骤2:k次迭代,并判断收敛条件;
步骤3:计算负梯度
步骤4:线性搜索得出,
步骤5:更新倾角:σk+1=PΓ(σk+ξksk);
步骤6:设置k:=k+1,返回步骤2.
其中,步长参数ξ,可通过Powell-Wolfe线性搜索准则获得。
实施例
通过实际采集的三维地震数据实例分析,说明一种基于正则化技术的地震数据倾角估计方法在倾角估计上的应用。图3给出的是,三维叠前时间偏移成像结果,图中2.5秒位置为风化壳界面,3秒为“串珠状”溶洞发育位置。对三维成像地震数据,进行正则化地震倾角估计后,得到的结果如图4所示,图中带状分布地震数据倾角场表示平滑反射层,而“串珠状”溶洞、断层等不连续地质体位置,倾角变化剧烈。可见,估计出的地震数据倾角能够在一定程度上,反映地层连续性情况。
Claims (8)
1.一种基于正则化技术的地震数据倾角估计方法,其特征在于,包括以下步骤:
步骤101):根据地震数据,选取其中任意一点作为地震数据倾角估计位置;
步骤102):依据全通数字滤波器计算出邻近地震道滤波器系数,进而计算出当前点地震数据预测值,其中包含待求地震数据倾角;
步骤103):循环执行上述未选取的地震数据位置,计算出所有位置地震数据预测值;
步骤104):根据所获取的地震数据、预测的地震数据及地震数据倾角构建正则化模型;
步骤105):利用非线性迭代反演算法获得地震数据倾角。
2.根据权利要求1所述的一种基于正则化技术的地震数据倾角估计方法,其特征在于,所述步骤102中计算出当前点地震数据预测值包括以下步骤:
由全通数字滤波器设计出不同阶数的滤波器系数,并根据滤波器系数、地震数据倾角和获取的地震数据,计算出当前点地震数据预测值。
3.根据权利要求1所述的一种基于正则化技术的地震数据倾角估计方法,其特征在于,所述步骤104中,构建正则化模型包括以下步骤:
根据所获取的地震数据和模拟数据预测值求最小二乘残差,作为模型构建第一项;约束的地震数据倾角作为第二项,构建正则化模型,如下所示:
其中,Jα(σ)为目标函数,d表示获取的地震数据,σ为待估计的地震数据倾角,P为全通数字滤波器破坏算子,α为正则化因子,L为拉普拉斯算子,min表示最小化。
4.根据权利要求1所述的一种基于正则化技术的地震数据倾角估计方法,其特征在于,所述步骤105中,非线性迭代反演算法,通过计算梯度、方向参数实现迭代,包括以下步骤:
对上式Jα(σ)求梯度,得出:
其中,为梯度,且
且当仅且得出倾角σ临界点σ*:
其中,步长参数ξ>0,PΓ表示投影,由上式,给出以下固定点迭代公式:
迭代过程如下所示:
步骤1:输入初始倾角σ=0,设定迭代序号k:=0;
步骤2:k次迭代,并判断收敛条件;
步骤3:计算负梯度
步骤4:线性搜索得出,
步骤5:更新倾角:σk+1=PΓ(σk+ξksk);
步骤6:设置k:=k+1,返回步骤2.
其中,步长参数ξ,可通过Powell-Wolfe线性搜索准则获得。
5.一种基于正则化技术的地震数据倾角估计装置,其特征在于,包括:
地震数据获取器件(201),用于从储存介质获取地震数据;
地震数据预测器件(202),用于依据全通数字滤波器计算出邻近地震道滤波器系数,进而计算出当前点地震数据预测值;
循环执行器件(203),用于循环执行所述未选取的地震数据位置,计算出所有位置地震数据预测值;
正则化模型构建器件(204),用于根据所述获取的地震数据、预测的地震数据及地震数据倾角构建正则化模型;
非线性求解器件(205),用于根据非线性迭代反演算法获得地震数据倾角。
6.根据权利要求5所述的一种基于正则化技术的地震数据倾角估计装置,其特征在于,所述地震数据预测器件(202)件包括:
滤波器系数计算单元,由全通数字滤波器设计出不用阶数的滤波器系数;
数值计算单元,根据滤波器系数、地震数据倾角和获取的地震数据,计算出当前点地震数据预测值。
7.根据权利要求5所述的一种基于正则化技术的地震数据倾角估计装置,其特征在于,所述正则化模型构建器件(204)包括:
最小二乘单元,根据所获取的地震数据和模拟预测值求最小二乘残差,作为模型构建第一项;
倾角单元,根据约束的地震数据倾角作为第二项,构建正则化估计模型,如下所示:
其中,Jα(σ)为目标函数,d表示获取的地震数据,σ为待估计的地震数据倾角,P为全通数字滤波器破坏算子,α为正则化因子,L为拉普拉斯算子,min表示最小化。
8.根据权利要求5所述的一种基于正则化技术的地震数据倾角估计装置,其特征在于,所述非线性求解器件(205),通过计算梯度、方向的参数实现,迭代过程如下:
对上式Jα(σ)求梯度,得出:
其中,为梯度,且
且当仅且得出倾角σ临界点σ*:
其中,步长参数ξ>0,PΓ表示投影,由上式,给出以下固定点迭代公式:
迭代过程如下所示:
步骤1:输入初始倾角σ=0,设定迭代序号k:=0;
步骤2:k次迭代,并判断收敛条件;
步骤3:计算负梯度
步骤4:线性搜索得出,
步骤5:更新倾角:σk+1=PΓ(σk+ξksk);
步骤6:设置k:=k+1,返回步骤2.
其中,步长参数ξ,可通过Powell-Wolfe线性搜索准则获得。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510106959.8A CN104749625B (zh) | 2015-03-11 | 2015-03-11 | 一种基于正则化技术的地震数据倾角估计方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510106959.8A CN104749625B (zh) | 2015-03-11 | 2015-03-11 | 一种基于正则化技术的地震数据倾角估计方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104749625A true CN104749625A (zh) | 2015-07-01 |
CN104749625B CN104749625B (zh) | 2016-10-05 |
Family
ID=53589595
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510106959.8A Active CN104749625B (zh) | 2015-03-11 | 2015-03-11 | 一种基于正则化技术的地震数据倾角估计方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104749625B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109387835A (zh) * | 2018-08-28 | 2019-02-26 | 国家能源投资集团有限责任公司 | 探地雷达绕射波弱信号提取方法及电子设备 |
CN109581497A (zh) * | 2018-10-29 | 2019-04-05 | 中国石油天然气集团有限公司 | 层析反演中的双向构造倾角约束方法及装置 |
CN109581496A (zh) * | 2018-10-29 | 2019-04-05 | 中国石油天然气集团有限公司 | 层析反演中的双向构造倾角约束方法及装置 |
CN113009564A (zh) * | 2021-03-24 | 2021-06-22 | 中国石油大学(北京) | 地震数据处理方法和装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060122780A1 (en) * | 2002-11-09 | 2006-06-08 | Geoenergy, Inc | Method and apparatus for seismic feature extraction |
US7548487B2 (en) * | 2007-01-17 | 2009-06-16 | Pgs Geophysical As | Diagonal gather trace interpolation |
CN103364833A (zh) * | 2013-07-01 | 2013-10-23 | 西安交通大学 | 一种高精度地层倾角估计方法 |
CN104216011A (zh) * | 2013-06-05 | 2014-12-17 | 上海青凤致远地球物理地质勘探科技有限公司 | 一种TTI介质稳定的qP波逆时偏移方法 |
-
2015
- 2015-03-11 CN CN201510106959.8A patent/CN104749625B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060122780A1 (en) * | 2002-11-09 | 2006-06-08 | Geoenergy, Inc | Method and apparatus for seismic feature extraction |
US7548487B2 (en) * | 2007-01-17 | 2009-06-16 | Pgs Geophysical As | Diagonal gather trace interpolation |
US7672195B2 (en) * | 2007-01-17 | 2010-03-02 | Pgs Geophysical As | Diagonal gather trace interpolation |
CN104216011A (zh) * | 2013-06-05 | 2014-12-17 | 上海青凤致远地球物理地质勘探科技有限公司 | 一种TTI介质稳定的qP波逆时偏移方法 |
CN103364833A (zh) * | 2013-07-01 | 2013-10-23 | 西安交通大学 | 一种高精度地层倾角估计方法 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109387835A (zh) * | 2018-08-28 | 2019-02-26 | 国家能源投资集团有限责任公司 | 探地雷达绕射波弱信号提取方法及电子设备 |
CN109581497A (zh) * | 2018-10-29 | 2019-04-05 | 中国石油天然气集团有限公司 | 层析反演中的双向构造倾角约束方法及装置 |
CN109581496A (zh) * | 2018-10-29 | 2019-04-05 | 中国石油天然气集团有限公司 | 层析反演中的双向构造倾角约束方法及装置 |
CN113009564A (zh) * | 2021-03-24 | 2021-06-22 | 中国石油大学(北京) | 地震数据处理方法和装置 |
CN113009564B (zh) * | 2021-03-24 | 2022-05-17 | 中国石油大学(北京) | 地震数据处理方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN104749625B (zh) | 2016-10-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Nangia et al. | A DLM immersed boundary method based wave-structure interaction solver for high density ratio multiphase flows | |
Fuhrman et al. | Tsunami generation, propagation, and run-up with a high-order Boussinesq model | |
Zhuang et al. | Accurate fracture modelling using meshless methods, the visibility criterion and level sets: formulation and 2D modelling | |
Umeyama | Eulerian–Lagrangian analysis for particle velocities and trajectories in a pure wave motion using particle image velocimetry | |
US20170338802A1 (en) | Actually-measured marine environment data assimilation method based on sequence recursive filtering three-dimensional variation | |
CN109633745B (zh) | 一种三维构造图的制图方法及装置 | |
CN107407736B (zh) | 生成无多次波的数据集的多阶段全波场反演处理 | |
JP7142968B2 (ja) | フルウェーブフォームインバージョン方法、装置及び電子設備 | |
CN107016649A (zh) | 一种基于局部低秩张量估计的视觉数据补全方法 | |
CN111859748B (zh) | 一种基于垂向混合坐标的海洋内波模拟方法 | |
Cui et al. | Optimal dispersion with minimized Poisson equations for non-hydrostatic free surface flows | |
CN104749625B (zh) | 一种基于正则化技术的地震数据倾角估计方法及装置 | |
Gao et al. | ChannelSeg3D: Channel simulation and deep learning for channel interpretation in 3D seismic images | |
CN105844637A (zh) | 基于非局部cv模型的sar图像变化检测方法 | |
CN115659848B (zh) | 一种基于深度学习网络快速预测二维盆地基底界面的方法 | |
Chen et al. | Estimation of spatially varying open boundary conditions for a numerical internal tidal model with adjoint method | |
Galley et al. | Geophysical inversion for 3D contact surface geometry | |
Asai et al. | Coupled tsunami simulations based on a 2d shallow-water equation-based finite difference method and 3d incompressible smoothed particle hydrodynamics | |
Yang et al. | A robust solver for incompressible high-Reynolds-number two-fluid flows with high density contrast | |
Chen et al. | Theoretical and experimental study of particle trajectories for nonlinear water waves propagating on a sloping bottom | |
EP3217354B1 (en) | Seismic modeling system providing seismic survey data frequency domain inpainting and related methods | |
Pudjaprasetya et al. | A nonhydrostatic two-layer staggered scheme for transient waves due to anti-symmetric seabed thrust | |
Liao et al. | Data assimilation for strongly nonlinear problems by transformed ensemble Kalman filter | |
CN107479091A (zh) | 一种提取逆时偏移角道集的方法 | |
Tarwidi et al. | A reduced two-layer non-hydrostatic model for submarine landslide-generated tsunamis |
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 |