CN105549080A - 一种基于辅助坐标系的起伏地表波形反演方法 - Google Patents
一种基于辅助坐标系的起伏地表波形反演方法 Download PDFInfo
- Publication number
- CN105549080A CN105549080A CN201610037642.8A CN201610037642A CN105549080A CN 105549080 A CN105549080 A CN 105549080A CN 201610037642 A CN201610037642 A CN 201610037642A CN 105549080 A CN105549080 A CN 105549080A
- Authority
- CN
- China
- Prior art keywords
- eta
- delta
- big gun
- centerdot
- gun record
- 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
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000004088 simulation Methods 0.000 claims description 27
- 238000000354 decomposition reaction Methods 0.000 claims description 8
- 230000009466 transformation Effects 0.000 claims description 7
- 238000013507 mapping Methods 0.000 claims description 6
- 239000002245 particle Substances 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 239000004576 sand Substances 0.000 claims description 3
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 abstract description 3
- 239000003208 petroleum Substances 0.000 abstract description 2
- 230000008569 process Effects 0.000 description 7
- 238000013508 migration Methods 0.000 description 5
- 230000005012 migration Effects 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000012937 correction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000000844 transformation Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000007634 remodeling Methods 0.000 description 1
- 230000001131 transforming effect Effects 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
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)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
Abstract
本发明公开了一种基于辅助坐标系的起伏地表波形反演方法,属于石油地球物理勘探技术领域,本发明首先采用辅助坐标系下的早至波波形反演方法对曲网格区域进行反演,以得到更新的近地表速度,再采用辅助坐标系下的全波形反演全局速度场进行速度更新,以克服近地表速度不准确对深部速度场反演的影响;采用时间域多尺度反演方法从低频到高频进行速度反演,以克服波形反演方法对初始速度模型的依赖性。本发明能够很好地反演剧烈起伏地表的速度场,为高精度成像方法提供准确的偏移速度场。
Description
技术领域
本发明属于石油地球物理勘探领域,具体涉及一种基于辅助坐标系的起伏地表波形反演方法。
背景技术
剧烈起伏地表对地震资料采集、处理带来了严重的影响,地球物理研究人员发展了一系列方法来克服这个问题。目前,针对起伏地表的处理主要有两种策略:一是对表层波场进行校正,二是基于起伏地表进行深度偏移成像。但复杂条件下静校正量难以准确计算,且静校正无法彻底消除地表起伏对地震波场造成的扭曲,因此直接基于起伏地表的深度偏移成像方法逐渐成为研究的热点。
由于全波形反演高精度以及高分辨率的特点,使其成为速度建模的一种有力工具,逐渐成为研究的热点。全波形反演是一个非线性数据拟合的过程,通过减少观测数据与预测数据的之间的差值来更新参数模型,这个过程以迭代的方式重复下去,直到数据差值足够小为止。
发明内容
针对现有技术中存在的上述技术问题,本发明提出了一种基于辅助坐标系的起伏地表波形反演方法,设计合理,具有良好的推广价值。
为了实现上述目的,本发明采用如下技术方案:
一种基于辅助坐标系的起伏地表波形反演方法,按照如下步骤进行:
步骤1:输入初始全局速度场、常规叠前炮记录、起伏高程及震源子波,并建立观测系统;
步骤2:根据初始全局速度场及起伏高程进行网格剖分,将近地表附近的网格剖分成曲网格,深层的网格剖分成矩形网格;
步骤3:将初始全局速度场变换到辅助坐标系下的矩形网格,采用下式所示的变换格式:
其中,x和z表示笛卡尔坐标系下的横纵坐标;ξ和η表示辅助坐标系下的横纵坐标;zi-1(ξ)和zi(ξ)是笛卡尔坐标系下第i层顶界面、底界面的高程,定义最深层为零;ηi-1(ξ)和ηi(ξ)是辅助坐标系第i层顶界面、底界面的纵向采样点数,定义最深层为零;
步骤4:对常规叠前炮记录划分时窗,在辅助坐标系下对近地表附近的曲网格区域速度场应用早至波波形反演,更新近地表的速度场,早至波波形反演的梯度方向如下:
其中,g为梯度;xs为炮点坐标;v为速度场的速度;p为声压;p*是残差波场的反向传播;t为时间;Te表示早至波时窗;
步骤5:判断使用近地表速度场正演模拟的炮记录与步骤4中划分时窗后的常规叠前炮记录之差是否满足误差条件;
若:判断结果是使用近地表速度场正演模拟的炮记录与步骤4中划分时窗后的常规叠前炮记录之差满足误差条件,则近地表速度场更新完成,然后执行步骤6;
或判断结果是使用近地表速度场正演模拟的炮记录与步骤4中划分时窗后的常规叠前炮记录之差不满足误差条件,则执行步骤4;
步骤6:将步骤1中输入的常规叠前炮记录分解成不同主频的多尺度炮记录,采用如下式所示的分解公式:
其中,f是维纳滤波器;Wo表示初始炮记录子波;Wt是生成的炮记录子波;ω是角频率;ε为一个小数;*表示共轭转置;
步骤7:将更新过近地表速度的全局速度场作为初始速度场,在辅助坐标系下应用全波形反演从低频到高频更新全局速度场,全波形反演的梯度方向如下:
其中,f表示主频;f1和fmax为多尺度分解的最低频率和最高频率;Tmax表示常规叠前炮记录的最大记录时间;表示主频是f的残差波场反传;
步骤8:判断使用全局速度场正演模拟的炮记录与常规叠前炮记录之差是否满足误差条件;
若:判断结果是使用全局速度场正演模拟的炮记录与常规叠前炮记录之差满足误差条件,则全局速度场更新完成,然后执行步骤9;
或判断结果是使用全局速度场正演模拟的炮记录与常规叠前炮记录之差不满足误差条件,则执行步骤7;
步骤9:将更新完成的全局速度场反变换到笛卡尔坐标系下,反变换公式如下:
步骤10:输出反演的速度场。
优选地,在步骤4中,具体包括
步骤4.1:定义目标函数:
其中,u(t,xr,xs)代表模拟波场u=(u,w,p)T,其中T表示转置;R为限定算子;dobs(t,xr,xs)是常规叠前炮记录;xs和xr表示震源点和检波点的位置坐标;t表示时间;E为目标函数值;
步骤4.2:将目标函数变分,得到变分表达式:
步骤4.3:定义变换格式:
通过变换格式(3)与锁链法则得到下面的映射公式:
通过映射公式(4),得到辅助坐标系下的一阶方程:
其中,p是声压;u和w分别是水平方向和垂直方向的质点速度;v是介质速度;f表示震源项;
速度的扰动δv可以引起地震波场的扰动δu,δu=(δu,δw,δp)T,将v+δv→u+δu代入一阶方程(5)后与一阶方程(5)相减得到如下方程:
进一步可得:
其中,L代表正演算子;
步骤4.4:将方程(7)代入变分表达式(2),可得:
其中,L*R*(Ru-dobs)表示波场的逆时传播;
利用伴随状态法,L*R*(Ru-dobs)可由下式求得:
则L*R*(Ru-dobs)=p*;
步骤4.5:求取目标函数对速度模型的梯度,得到早至波波形反演的梯度方向:
其中,g为梯度;xs为炮点坐标;v为速度场的速度;p为声压;p*是残差波场的反向传播;t为时间;Te表示早至波时窗。
本发明所带来的有益技术效果:
本发明提出了一种基于辅助坐标系的起伏地表波形反演方法,与现有技术相比,一种基于辅助坐标系的起伏地表波形反演方法,采用一阶方程变密度方程克服传统常密度二阶方程在处理密度变化比较大地区的速度反演不准确的缺点;同时在起伏地表处采用曲网格剖分,在深部区域采用矩形网格剖分,并统一变换到辅助坐标系下的矩形网格进行计算,克服起伏地表对地震波传播的影响;首先采用辅助坐标系下的早至波波形反演方法对曲网格区域进行反演,以得到更新的近地表速度,再采用辅助坐标系下的全波形反演全局速度场进行速度更新,以克服近地表速度不准确对深部速度场反演的影响;采用时间域多尺度反演方法从低频到高频进行速度反演,以克服波形反演方法对初始速度模型的依赖性。本发明能够很好地反演剧烈起伏地表的速度场,为高精度成像方法提供准确的偏移速度场。
附图说明
图1为本发明的一种基于辅助坐标系的起伏地表波形反演方法的流程图。
图2为本发明使用的真实起伏地表速度模型。
图3为变换到辅助坐标系下的速度模型。
图4为采用统一曲网格变换到辅助坐标系下的速度模型。
图5为初始全局速度场。
图6为输入的常规叠前炮记录。
图7为采用统一曲网格正演模拟得到的炮记录。
图8为多尺度分解得到的5Hz叠前炮记录。
图9为多尺度分解得到的15Hz叠前炮记录。
图10为采用本发明方法得到的第一尺度反演速度场(主频为5Hz)。
图11为采用本发明方法得到的第二尺度反演速度场(主频为15Hz)。
图12为采用本发明方法得到的第三尺度反演速度场(主频为25Hz)。
图13为采用单尺度波形反演方法得到的速度场。
图14为采用统一曲网格得到的波形反演速度场。
具体实施方式
下面结合附图以及具体实施方式对本发明作进一步详细说明:
一种基于辅助坐标系的起伏地表波形反演方法的流程图(如图1所示),包括如下步骤:
输入初始全局速度场、常规叠前炮记录、起伏高程及震源子波,并建立观测系统;
根据初始全局速度场及起伏高程进行网格剖分,近地表附近的网格剖分成曲网格,深层的网格剖分成矩形网格;
将初始全局速度场变换到辅助坐标系下的矩形网格;
对常规叠前炮记录划分时窗,在辅助坐标系下对近地表附近的曲网格区域速度场应用早至波波形反演,更新近地表的速度场,判断使用近地表速度场正演模拟的炮记录与划分时窗的常规叠前炮记录的差是否满足误差条件,如果满足则近地表速度场更新完成,如果不满足再次应用早至波波形反演更新近地表速度场;
将输入的常规叠前炮记录分解成不同主频的多尺度炮记录,将更新过近地表速度的全局速度场作为初始速度场,在辅助坐标系下应用全波形反演从低频到高频更新全局速度场,判断使用全局速度场正演模拟的炮记录与常规叠前炮记录的差是否满足误差条件,如果满足则全局速度场更新完成,如果不满足再次应用全波形波形反演更新全局速度场;
将更新完成的全局速度场变换到笛卡尔坐标系下,并输出反演的速度场。
具体步骤为:
步骤1:输入初始全局速度场、常规叠前炮记录、起伏高程及震源子波,并建立观测系统;
步骤2:根据初始全局速度场及起伏高程进行网格剖分,将近地表附近的网格剖分成曲网格,深层的网格剖分成矩形网格;
步骤3:将初始全局速度场变换到辅助坐标系下的矩形网格,采用下式所示的变换格式:
其中,x和z表示笛卡尔坐标系下的横纵坐标;ξ和η表示辅助坐标系下的横纵坐标;zi-1(ξ)和zi(ξ)是笛卡尔坐标系下第i层顶界面、底界面的高程,定义最深层为零;ηi-1(ξ)和ηi(ξ)是辅助坐标系第i层顶界面、底界面的纵向采样点数,定义最深层为零。
通过变换格式(1)与锁链法则得到下面的映射公式:
步骤4:对常规叠前炮记录划分时窗,在辅助坐标系下对近地表附近的曲网格区域速度场应用早至波波形反演,更新近地表的速度场,早至波波形反演的梯度方向如下:
通过映射公式(2),得到辅助坐标系下的一阶方程:
其中,p是声压;u和w分别是水平方向和垂直方向的质点速度;v是介质速度;f表示震源项。
对于波形反演需要构建目标函数。
本发明采用常规叠前炮记录与模拟炮记录残差的L2模作为目标函数:
其中,u(t,xr,xs)代表模拟波场u=(u,w,p)T,其中T表示转置;R为限定算子;dobs(t,xr,xs)是常规叠前炮记录;xs和xr表示震源点和检波点的位置坐标;t表示时间;E为目标函数值;
将目标函数变分,得到变分表达式:
速度的扰动δv可以引起地震波场的扰动δu,δu=(δu,δw,δp)T,将v+δv→u+δu代入辅助坐标系下的一阶方程(3)中并与之相减得到如下方程:
从方程(6),我们可以得到其中L代表正演算子,方程(5)进一步可表示为:
其中,L*R*(Ru-dobs)表示波场的逆时传播;
利用伴随状态法,L*R*(Ru-dobs)可由下式求得:
则L*R*(Ru-dobs)=p*;
早至波波形反演的梯度方向如下:
其中,g为梯度;xs为炮点坐标;v为速度场的速度;p为声压;p*是残差波场的反向传播;t为时间;Te表示早至波时窗;
步骤5:判断使用近地表速度场正演模拟的炮记录与步骤4中划分时窗后的常规叠前炮记录之差是否满足误差条件;
若:判断结果是使用近地表速度场正演模拟的炮记录与步骤4中划分时窗后的常规叠前炮记录之差满足误差条件,则近地表速度场更新完成,然后执行步骤6;
或判断结果是使用近地表速度场正演模拟的炮记录与步骤4中划分时窗后的常规叠前炮记录之差不满足误差条件,则执行步骤4;
步骤6:将步骤1中输入的常规叠前炮记录分解成不同主频的多尺度炮记录,采用如下式所示的分解公式:
其中,f是维纳滤波器;Wo表示初始炮记录子波;Wt是生成的炮记录子波;ω是角频率;ε为一个小数;*表示共轭转置;
步骤7:将更新过近地表速度的全局速度场作为初始速度场,在辅助坐标系下应用全波形反演从低频到高频更新全局速度场,全波形反演的梯度方向如下:
其中,f表示主频;f1和fmax为多尺度分解的最低频率和最高频率;Tmax表示常规叠前炮记录的最大记录时间;表示主频是f的残差波场反传;
步骤8:判断使用全局速度场正演模拟的炮记录与常规叠前炮记录之差是否满足误差条件;
若:判断结果是使用全局速度场正演模拟的炮记录与常规叠前炮记录之差满足误差条件,则全局速度场更新完成,然后执行步骤9;
或判断结果是使用全局速度场正演模拟的炮记录与常规叠前炮记录之差不满足误差条件,则执行步骤7;
步骤9:将更新完成的全局速度场反变换到笛卡尔坐标系下,反变换公式如下:
步骤10:输出反演的速度场。
本发明的一种基于辅助坐标系的起伏地表波形反演方法,能够很好地反演剧烈起伏地表的速度场,为高精度成像方法提供准确的偏移速度场。
应用实验
本发明一种基于辅助坐标系的起伏地表波形反演方法,应用于复杂起伏地表模型数据。图2为本发明使用的真实起伏地表速度模型;图3为变换到辅助坐标系下速度模型;图4为采用统一曲网格变换到辅助坐标系下的速度模型;应用本发明的网格剖分策略对近地表附近的网格剖分成曲网格,深层的网格剖分成矩形网格,再变换到辅助坐标系下的速度模型(图3)与应用传统统一曲网格剖分,再变换到辅助坐标系下的速度模型(图4)相比,本发明的结果更好地保存了深部构造的原始形态。图5为输入初始全局速度场,图6为输入的常规叠前炮记录,也是采用本发明网格剖分策略正演模拟得到的炮记录,与采用统一曲网格剖分得到的叠前炮记录(图7)对比,本发明方法得到的炮记录同相轴更加清晰,而且虚假的绕射波和散射波信息更少,因为波形反演的基础是正演模拟,通过正演模拟的结果对比也能体现本发明的优势。将输入的常规叠前炮记录分解成不同主频的多尺度炮记录,图8为主频为5Hz的炮记录,图9为主频为15Hz的炮记录。在辅助坐标系下应用全波形反演从低频到高频更新全局速度场,图10为采用本发明方法得到的第一尺度反演速度场(主频为5Hz),图11采用本发明方法得到的第二尺度反演速度场(主频为15Hz),图12采用本发明方法得到的第三尺度反演速度场(主频为25Hz),从三图的对比可以看出,采用第一尺度得到的反演速度场,低频信息得到了很好地恢复,随着反演频率逐渐升高,高频信息也得到了很好地改善,反演结果非常接近于真实速度模型(图2),作为对比,这里给出采用单尺度波形反演方法得到的速度场(图13),因为初始速度模型(图5)极不准确,所以很难得到较好的反演结果,层位信息与速度信息都没有很好地反演出来。本发明的反演结果(图12)相比于采用统一曲网格得到的波形反演速度场(图14)也更加准确。
为此本发明提供了一种基于辅助坐标系的起伏地表波形反演方法,能够很好地反演剧烈起伏地表的速度场,为高精度成像方法提供准确的偏移速度场。
当然,上述说明并非是对本发明的限制,本发明也并不仅限于上述举例,本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换,也应属于本发明的保护范围。
Claims (2)
1.一种基于辅助坐标系的起伏地表波形反演方法,其特征在于:按照如下步骤进行:
步骤1:输入初始全局速度场、常规叠前炮记录、起伏高程及震源子波,并建立观测系统;
步骤2:根据初始全局速度场及起伏高程进行网格剖分,将近地表附近的网格剖分成曲网格,深层的网格剖分成矩形网格;
步骤3:将初始全局速度场变换到辅助坐标系下的矩形网格,采用下式所示的变换格式:
其中,x和z表示笛卡尔坐标系下的横纵坐标;ξ和η表示辅助坐标系下的横纵坐标;zi-1(ξ)和zi(ξ)是笛卡尔坐标系下第i层顶界面、底界面的高程,定义最深层为零;ηi-1(ξ)和ηi(ξ)是辅助坐标系第i层顶界面、底界面的纵向采样点数,定义最深层为零;
步骤4:对常规叠前炮记录划分时窗,在辅助坐标系下对近地表附近的曲网格区域速度场应用早至波波形反演,更新近地表的速度场,早至波波形反演的梯度方向如下:
其中,g为梯度;xs为炮点坐标;v为速度场的速度;p为声压;p*是残差波场的反向传播;t为时间;Te表示早至波时窗;
步骤5:判断使用近地表速度场正演模拟的炮记录与步骤4中划分时窗后的常规叠前炮记录之差是否满足误差条件;
若:判断结果是使用近地表速度场正演模拟的炮记录与步骤4中划分时窗后的常规叠前炮记录之差满足误差条件,则近地表速度场更新完成,然后执行步骤6;
或判断结果是使用近地表速度场正演模拟的炮记录与步骤4中划分时窗后的常规叠前炮记录之差不满足误差条件,则执行步骤4;
步骤6:将步骤1中输入的常规叠前炮记录分解成不同主频的多尺度炮记录,采用如下式所示的分解公式:
其中,f是维纳滤波器;Wo表示初始炮记录子波;Wt是生成的炮记录子波;ω是角频率;ε为一个小数;*表示共轭转置;
步骤7:将更新过近地表速度的全局速度场作为初始速度场,在辅助坐标系下应用全波形反演从低频到高频更新全局速度场,全波形反演的梯度方向如下:
其中,f表示主频;f1和fmax为多尺度分解的最低频率和最高频率;Tmax表示常规叠前炮记录的最大记录时间;表示主频是f的残差波场反传;
步骤8:判断使用全局速度场正演模拟的炮记录与常规叠前炮记录之差是否满足误差条件;
若:判断结果是使用全局速度场正演模拟的炮记录与常规叠前炮记录之差满足误差条件,则全局速度场更新完成,然后执行步骤9;
或判断结果是使用全局速度场正演模拟的炮记录与常规叠前炮记录之差不满足误差条件,则执行步骤7;
步骤9:将更新完成的全局速度场反变换到笛卡尔坐标系下,反变换公式如下:
步骤10:输出反演的速度场。
2.根据权利要求1所述的基于辅助坐标系的起伏地表波形反演方法,其特征在于:在步骤4中,具体包括
步骤4.1:定义目标函数:
其中,u(t,xr,xs)代表模拟波场u=(u,w,p)T,其中T表示转置;R为限定算子;dobs(t,xr,xs)是常规叠前炮记录;xs和xr表示震源点和检波点的位置坐标;t表示时间;E为目标函数值;
步骤4.2:将目标函数变分,得到变分表达式:
步骤4.3:定义变换格式:
通过变换格式(3)与锁链法则得到下面的映射公式:
通过映射公式(4),得到辅助坐标系下的一阶方程:
其中,p是声压;u和w分别是水平方向和垂直方向的质点速度;v是介质速度;f表示震源项;
将v+δv→u+δu代入一阶方程(5)后与一阶方程(5)相减得到如下方程:
进一步可得:
其中,L代表正演算子;
步骤4.4:将方程(7)代入变分表达式(2),可得:
其中,L*R*(Ru-dobs)表示波场的逆时传播;
利用伴随状态法,L*R*(Ru-dobs)可由下式求得:
则L*R*(Ru-dobs)=p*;
步骤4.5:求取目标函数对速度模型的梯度,得到早至波波形反演的梯度方向:
其中,g为梯度;xs为炮点坐标;v为速度场的速度;p为声压;p*是残差波场的反向传播;t为时间;Te表示早至波时窗。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610037642.8A CN105549080B (zh) | 2016-01-20 | 2016-01-20 | 一种基于辅助坐标系的起伏地表波形反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610037642.8A CN105549080B (zh) | 2016-01-20 | 2016-01-20 | 一种基于辅助坐标系的起伏地表波形反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105549080A true CN105549080A (zh) | 2016-05-04 |
CN105549080B CN105549080B (zh) | 2017-08-25 |
Family
ID=55828380
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610037642.8A Expired - Fee Related CN105549080B (zh) | 2016-01-20 | 2016-01-20 | 一种基于辅助坐标系的起伏地表波形反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105549080B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106054244A (zh) * | 2016-06-16 | 2016-10-26 | 吉林大学 | 截断时窗的低通滤波多尺度全波形反演方法 |
CN106898045A (zh) * | 2017-02-24 | 2017-06-27 | 郑州大学 | 一种基于sgog瓦块的大区域真三维地理场景自适应构建方法 |
CN106950596A (zh) * | 2017-04-11 | 2017-07-14 | 中国石油大学(华东) | 一种基于子波迭代估计的有限差分对比源全波形反演方法 |
CN106383361B (zh) * | 2016-08-31 | 2018-12-11 | 中国石油集团东方地球物理勘探有限责任公司 | 一种速度数据网格更新的方法 |
CN109031413A (zh) * | 2018-05-15 | 2018-12-18 | 中国石油大学(华东) | 一种基于起伏海底电缆数据的矢量波逆时偏移系统及方法 |
CN109085643A (zh) * | 2018-07-30 | 2018-12-25 | 中国石油化工股份有限公司 | 早至波的分步联合反演方法 |
CN109188519A (zh) * | 2018-09-17 | 2019-01-11 | 中国石油大学(华东) | 一种极坐标下的弹性波纵横波速度反演系统及方法 |
CN113552625A (zh) * | 2021-06-21 | 2021-10-26 | 中国地质科学院地球物理地球化学勘查研究所 | 一种用于常规陆域地震数据的多尺度全波形反演方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102183790A (zh) * | 2011-02-12 | 2011-09-14 | 中国石油大学(华东) | 基于时空双变网格的弹性波正演模拟技术 |
CN104391323A (zh) * | 2014-11-21 | 2015-03-04 | 中国石油大学(华东) | 一种利用反射波信息反演速度场中低波数成分的方法 |
CN104977607A (zh) * | 2014-04-09 | 2015-10-14 | 中国石油集团东方地球物理勘探有限责任公司 | 利用变步长网格声波波场模拟的时间域全波形反演方法 |
US20150362622A1 (en) * | 2014-06-17 | 2015-12-17 | Huseyin Denli | Fast Viscoacoustic and Viscoelastic Full Wavefield Inversion |
-
2016
- 2016-01-20 CN CN201610037642.8A patent/CN105549080B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102183790A (zh) * | 2011-02-12 | 2011-09-14 | 中国石油大学(华东) | 基于时空双变网格的弹性波正演模拟技术 |
CN104977607A (zh) * | 2014-04-09 | 2015-10-14 | 中国石油集团东方地球物理勘探有限责任公司 | 利用变步长网格声波波场模拟的时间域全波形反演方法 |
US20150362622A1 (en) * | 2014-06-17 | 2015-12-17 | Huseyin Denli | Fast Viscoacoustic and Viscoelastic Full Wavefield Inversion |
WO2015199800A1 (en) * | 2014-06-17 | 2015-12-30 | Exxonmobil Upstream Research Company | Fast viscoacoustic and viscoelastic full-wavefield inversion |
CN104391323A (zh) * | 2014-11-21 | 2015-03-04 | 中国石油大学(华东) | 一种利用反射波信息反演速度场中低波数成分的方法 |
Non-Patent Citations (5)
Title |
---|
YINGMING QU ET AL.: "Elastic wave modeling and wave field separation of irregular free-surface based on multi-block mapping method", 《NEAR-SURFACE ASIA PACIFIC CONFERENCE》 * |
YINGMING QU ET AL.: "Multisource elastic full waveform inversion method for irregular surface", 《2015 WORKSHOP:DEPTH MODEL BUILDING:FULL-WAVEFORM INVERSION》 * |
曲英铭等: "分层坐标变换法起伏自由地表弹性波叠前逆时偏移", 《地球物理学报》 * |
曲英铭等: "基于单元交错网格的变坐标系正演模拟方法在声-弹介质中的应用", 《石油物探》 * |
曲英铭等: "时间域多重双变网格精细全波形反演", 《中国地球科学联合学术年会 2014》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106054244A (zh) * | 2016-06-16 | 2016-10-26 | 吉林大学 | 截断时窗的低通滤波多尺度全波形反演方法 |
CN106054244B (zh) * | 2016-06-16 | 2017-11-28 | 吉林大学 | 截断时窗的低通滤波多尺度全波形反演方法 |
CN106383361B (zh) * | 2016-08-31 | 2018-12-11 | 中国石油集团东方地球物理勘探有限责任公司 | 一种速度数据网格更新的方法 |
CN106898045A (zh) * | 2017-02-24 | 2017-06-27 | 郑州大学 | 一种基于sgog瓦块的大区域真三维地理场景自适应构建方法 |
CN106898045B (zh) * | 2017-02-24 | 2020-02-07 | 郑州大学 | 一种基于sgog瓦块的大区域真三维地理场景自适应构建方法 |
CN106950596A (zh) * | 2017-04-11 | 2017-07-14 | 中国石油大学(华东) | 一种基于子波迭代估计的有限差分对比源全波形反演方法 |
CN109031413A (zh) * | 2018-05-15 | 2018-12-18 | 中国石油大学(华东) | 一种基于起伏海底电缆数据的矢量波逆时偏移系统及方法 |
CN109031413B (zh) * | 2018-05-15 | 2019-12-06 | 中国石油大学(华东) | 一种基于起伏海底电缆数据的矢量波逆时偏移系统及方法 |
CN109085643A (zh) * | 2018-07-30 | 2018-12-25 | 中国石油化工股份有限公司 | 早至波的分步联合反演方法 |
CN109188519A (zh) * | 2018-09-17 | 2019-01-11 | 中国石油大学(华东) | 一种极坐标下的弹性波纵横波速度反演系统及方法 |
CN113552625A (zh) * | 2021-06-21 | 2021-10-26 | 中国地质科学院地球物理地球化学勘查研究所 | 一种用于常规陆域地震数据的多尺度全波形反演方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105549080B (zh) | 2017-08-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105549080A (zh) | 一种基于辅助坐标系的起伏地表波形反演方法 | |
CN103713315B (zh) | 一种地震各向异性参数全波形反演方法及装置 | |
CN103703391B (zh) | 使用频谱整形的全波场反演的系统和计算机实施的方法 | |
CN108710153B (zh) | 一种磁全张量梯度反演地下三维磁性分布的波数域方法 | |
CN103238158B (zh) | 利用互相关目标函数进行的海洋拖缆数据同时源反演 | |
CN107843925B (zh) | 一种基于修正相位的反射波波形反演方法 | |
CN103499835A (zh) | 初至波波形反演近地表速度模型方法 | |
CN107203002A (zh) | 反演速度模型及其建立方法和地下结构的像的获得方法 | |
CN105676277B (zh) | 一种提高高陡构造速度反演效率的全波形联合反演方法 | |
CN104614763A (zh) | 基于反射率法的多波avo储层弹性参数反演方法及系统 | |
CN104122585A (zh) | 基于弹性波场矢量分解与低秩分解的地震正演模拟方法 | |
CN110058307B (zh) | 一种基于快速拟牛顿法的全波形反演方法 | |
CN105005076A (zh) | 基于最小二乘梯度更新速度模型的地震波全波形反演方法 | |
CN104570082A (zh) | 一种基于格林函数表征的全波形反演梯度算子的提取方法 | |
CN111239819B (zh) | 一种基于地震道属性分析的带极性直接包络反演方法 | |
CN101021568A (zh) | 基于最大能量旅行时计算的三维积分叠前深度偏移方法 | |
CN110579795A (zh) | 基于被动源地震波形及其逆时成像的联合速度反演方法 | |
CN104965222A (zh) | 三维纵波阻抗全波形反演方法及装置 | |
CN109490978B (zh) | 一种起伏地层的频率域快速高精度正演方法 | |
CN111505714B (zh) | 基于岩石物理约束的弹性波直接包络反演方法 | |
CN112630830B (zh) | 一种基于高斯加权的反射波全波形反演方法及系统 | |
CN111273346B (zh) | 去除沉积背景的方法、装置、计算机设备及可读存储介质 | |
CN111999764A (zh) | 基于时频域目标函数的盐下构造最小二乘逆时偏移方法 | |
CN111190224B (zh) | 一种基于三维地震波反向照明的动态采样全波形反演系统及方法 | |
CN105319594A (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 | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20170613 Address after: 266580 Qingdao Changjiang Road, Huangdao District, Shandong, No. 66 Applicant after: CHINA University OF PETROLEUM (EAST CHINA) Address before: 266580 Qingdao Changjiang Road, Huangdao District, Shandong, No. 66 Applicant before: CHINA University OF PETROLEUM (EAST CHINA) Applicant before: Qu Yingming |
|
TA01 | Transfer of patent application right | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170825 |
|
CF01 | Termination of patent right due to non-payment of annual fee |