CN110261896B - 稳定的各向异性ti介质正演模拟方法 - Google Patents
稳定的各向异性ti介质正演模拟方法 Download PDFInfo
- Publication number
- CN110261896B CN110261896B CN201910343187.8A CN201910343187A CN110261896B CN 110261896 B CN110261896 B CN 110261896B CN 201910343187 A CN201910343187 A CN 201910343187A CN 110261896 B CN110261896 B CN 110261896B
- Authority
- CN
- China
- Prior art keywords
- medium
- equation
- anisotropic
- wave
- phase velocity
- 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
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000004088 simulation Methods 0.000 claims abstract description 6
- 230000006870 function Effects 0.000 claims description 8
- 238000012360 testing method Methods 0.000 claims description 4
- 239000011435 rock Substances 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 2
- 238000001228 spectrum Methods 0.000 abstract 1
- 238000011161 development Methods 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 6
- 239000006185 dispersion Substances 0.000 description 5
- 238000007796 conventional method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000010835 comparative analysis Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
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
- G01V1/282—Application of seismic models, synthetic seismograms
-
- 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
- G01V1/30—Analysis
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)
Abstract
本发明提供一种稳定的各向异性TI介质正演模拟方法,该稳定的各向异性TI介质正演模拟方法包括:步骤1,建立各向异性地质模型;步骤2,建立VTI介质中的纯拟声波相速度方程;步骤3,建立VTI和TTI介质的纯拟声波控制方程;步骤4,分析拟声波相速度曲线的精度,选择近似阶数;步骤5,求解选定阶数的控制方程,得到TI介质拟声波正演模拟的结果。该稳定的各向异性TI介质正演模拟方法从准确的相速度公式开始,并通过使用二次范数中的最佳近似来开发TI介质的精确纯拟声波控制方程,并采用伪谱法来求解方程,非常接近耦合的伪声波方程,产生的波场不仅稳定,而且可以很好地适应复杂的TI介质。
Description
技术领域
本发明涉及石油地球物理勘探技术领域,特别是涉及到一种稳定的各向异性TI介质正演模拟方法。
背景技术
在石油地球物理勘探技术领域,地震正演模拟是逆时偏移和全波形反演的基础。准确的各向异性正演模拟能够更精确地模拟地下构造地震波的传播形态。
随着科学技术的迅猛发展、大型计算机的广泛应用以及地震勘探仪器动态范围的不断扩大,简单、浅中层常规油气资源的开发利用日渐殆尽,勘探开发的重点逐步向构造复杂的油气藏以及岩性、裂缝性等非常规、隐蔽油气藏转移。复杂油气藏的勘探开发促进了地球科学领域中一系列新概念、新理论的产生,在对地球内部介质的结构与属性的研究中,人们对地球介质各向异性的认识逐步深化。
基于传统耦合伪声波方程的正演模拟方法不仅存在SV波串扰和数值频散,而且当各向异性参数和倾斜横向各向同性(TTI)介质具有高度变化的倾角和方位角时存在不稳定性。各向异性包括TTI,VTI,HTI等,TTI是具有倾斜对称轴的横向各向同性,VTI是具有垂直对称轴的横向各向同性,HTI是具有水平对称轴的横向各向同性。
因此,考虑到上述各向异性正演模拟技术存在的缺陷,本发明提出了一种稳定的各向异性TI(横向各向同性)介质正演模拟方法。
发明内容
本发明的目的是提供一种用于解决目前的各向异性正演模拟技术存在横波SV波串扰和数值频散及当各向异性参数和倾斜横向各向同性(TTI)介质具有高度变化的倾角和方位角时存在不稳定性的问题的稳定的各向异性TI介质正演模拟方法。
本发明的目的可通过如下技术措施来实现:稳定的各向异性TI介质正演模拟方法,该稳定的各向异性TI介质正演模拟方法包括:步骤1,建立各向异性地质模型;步骤2,建立VTI介质中的纯拟声波相速度方程;步骤3,建立VTI和TTI介质的纯拟声波控制方程;步骤4,分析拟声波相速度曲线的精度,选择近似阶数;步骤5,求解选定阶数的控制方程,得到TI介质拟声波正演模拟的结果。
本发明的目的还可通过如下技术措施来实现:
在步骤1中,结合地震和测井资料及其解释结果包括层位、断层,建立模型框架,之后加上各层的速度参数,建立各向异性速度场;根据测井资料或者岩石物理测试建立各向异性参数场;根据各向异性速度场,通过倾角估计建立倾角场。
在步骤2,基于精确的VTI介质相速度公式:
其中,相位速度,θ是沿轴对称方向测量的相位角,vp0和vs0分别是拟声波和qSV波的速度;ε和δ是Thomsen参数;正号对应拟声波,负号对应qSV波公式;通过设置qSV波速度沿着对称轴方向等于0,也就是f≡1;则VTI介质拟声波相速度方程可写为:
这里,引入了一个新的符号f(θ)表示方程(2)的平方根部分,并使得用二次范数中的最佳逼近来逼近方程(2)的平方根部分,f(θ)由下式求得:
二次范数最佳逼近的思想是通过寻找λ1,λ2,…,λn以满足如下方程:
四阶函数系统的逼近格式为:{sin4θ,sin2θ,1};越高阶函数系统将产生越高的近似精度,但它会牺牲更多的计算量和存储,特别是具有倾斜对称轴的横向各向同性TTI介质的情况;
通过二次范数最佳逼近,VTI介质拟声波相速度方程(2)可以近似为
方程(5)、(6)、(7)分别是二阶、四阶、六阶近似的VTI介质纯拟声波相速度方程。
在步骤3,将VTI介质的对称轴从垂直角度旋转到倾斜角度,得到TTI介质中的纯拟声波相速度方程;对于TTI介质,将对称轴从垂直旋转到倾斜角,根据方程(5)、(6)、(7)得到不同近似阶数的TTI介质拟声波相速度方程:
其中φ表示TTI介质的对称轴旋转角度;
把方程(6)和(9)与傅里叶域中的波场UP(kx,kz,t)相乘,再进行傅里叶逆变换,建立四阶近似下VTI介质和TTI介质的新的纯拟声波控制方程如下:
同样可建立其它近似精度的VTI和TTI介质控制方程;
其中,UP表示拟声波波场,kx、kz分别表示x和z方向的波数,t表示时间,φ表示TTI介质的对称轴旋转角度,其余参数意义与前相同。
在步骤4中,根据各向异性介质的准确相速度方程以及不同阶数的近似相速度方程,绘制准确的以及不同近似阶数的拟声波相速度随相位角变化的曲线,并进行对比分析,选择合适的近似阶数。
在步骤5中,输入各向异性速度场、参数场和倾角场,求解选定阶数的控制方程,得到TI介质拟声波正演模拟的结果。
本发明中的稳定的各向异性TI介质正演模拟方法,用于解决目前的各向异性正演模拟技术存在SV波串扰和数值频散及当各向异性参数和倾斜横向各向同性(TTI)介质具有高度变化的倾角和方位角时存在不稳定性的问题。该稳定的各向异性TI介质正演模拟方法从准确的相速度公式开始,并通过使用二次范数中的最佳近似来开发TI介质的精确纯拟声波控制方程,并采用伪谱法来求解方程。新波动方程的运动学非常接近耦合的伪声波方程,完全没有SV波伪影。并且新波动方程产生的波场不仅稳定,而且可以很好地适应复杂的TI介质。
附图说明
图1是本发明的实施方式中不同精度近似及真实的拟声波相速度曲线图;
图2是本发明的实施方式中均匀TTI介质波场快照图;
图3是本发明的实施方式中HESS VTI模型的示意图;
图4是本发明的实施方式中HESS VTI模型得到的0.56s的波场快照图;
图5是本发明的实施方式中BP TTI模型的示意图;
图6是本发明的实施方式中BP TTI模型得到的0.56s的波场快照图;
图7为本发明的稳定的各向异性TI介质正演模拟方法的一具体实施例的流程图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合附图所示,作详细说明如下。
如图7所示,图7为本发明的稳定的各向异性TI介质正演模拟方法的结构图。
步骤101,建立各向异性地质模型。结合地震和测井资料及其解释结果包括层位、断层等,建立模型框架,之后加上各层的速度参数,建立各向异性速度场;根据测井资料或者岩石物理测试建立各向异性参数场(ε场和δ场);根据各向异性速度场,通过倾角估计建立倾角场(VTI不需要倾角场,TTI需要)。
步骤102,建立具有垂直对称轴的横向各向同性VTI介质中的纯拟声波相速度方程。
基于精确的VTI介质相速度公式:
这里是v相位速度,θ是沿轴对称方向测量的相位角,vp0和vs0分别是拟声波和qSV波的速度。ε和δ是Thomsen参数。正号对应拟声波,负号对应qSV波公式。在这里,只讨论拟声波公式。通过设置qSV波速度沿着对称轴方向等于0,也就是f≡1。则VTI介质拟声波相速度方程可写为:
这里,引入了一个新的符号f(θ)表示方程(2)的平方根部分,并使得用二次范数中的最佳逼近来逼近方程(2)的平方根部分,f(θ)由下式求的
二次范数最佳逼近的思想是通过寻找λ1,λ2,…,λn以满足如下方程:
n是近似阶数,特别地,四阶函数系统的逼近格式为:{sin4θ,sin2θ,1}。当然,越高阶函数系统将产生越高的近似精度,但它会牺牲更多的计算量和存储,特别是具有倾斜对称轴的横向各向同性TTI介质的情况。
通过二次范数最佳逼近,VTI介质拟声波相速度方程(2)可以近似为
方程(5)、(6)、(7)分别是二阶、四阶、六阶近似的VTI介质纯拟声波相速度方程。
步骤103,建立VTI和TTI介质的纯拟声波控制方程。将VTI介质的对称轴从垂直角度旋转到倾斜角度,得到TTI介质中的纯拟声波相速度方程。对于TTI介质,将对称轴从垂直旋转到倾斜角,根据方程(5)、(6)、(7)可以得到不同近似阶数的TTI介质拟声波相速度方程:
把方程(6)和(9)与傅里叶域中的波场UP(kx,kz,t)相乘,再进行傅里叶逆变换,可以建立四阶近似下VTI介质和TTI介质的新的纯拟声波控制方程如下:
同样可建立其它近似精度的VTI和TTI介质控制方程。
步骤104,分析拟声波相速度曲线的精度,选择近似阶数。根据各向异性介质的准确相速度方程以及不同阶数的近似相速度方程,绘制准确的以及不同近似阶数的拟声波相速度随相位角变化的曲线,并进行对比分析,选择合适的近似阶数;
步骤105,输入各向异性速度场、参数场和倾角场,求解选定阶数的控制方程,得到TI介质拟声波正演模拟的结果。
在应用本发明的一具体实施例中,应用于HESS VTI模型和BP TTI模型数据,取得了理想的计算效果。推导VTI介质中的纯拟声波方程,将对称轴从垂直角度旋转到倾斜角度,推导出TTI介质中的纯拟声波方程,并分析拟声波相速度曲线的精度,选择一组典型的模型参数ε=0.24,δ=0.18,vp0=2000m/s进行了测试,不同精度近似及真实的拟声波相速度曲线如图1所示。通过对比可以看出,二阶近似的结果与真实相速度具有较大的误差,而四阶近似的结果很接近于真实结果,因此不再需要进行六阶近似了,因为六阶近似需要大量的计算量和内存。基于上述讨论,选择了四阶近似进行各向异性正演模拟。
首先对均匀TTI介质模型进行测试。图2所示的为三种均匀模型的波场快照。图2a-2c所示的为常规方法的波场快照;图2d-2f所示的为本发明方法的波场快照。图2a和2d所示的为vp0=2500m/s,ε=0.24,δ=0.12,θ=0°的结果,图2b和2e为vp0=2500m/s,ε=0.24,δ=0.12,θ=30°的结果,图2c和2f所示的为vp0=2500m/s,ε=0.12,δ=0.24,θ=30°的结果。同图2a-2d的结果可以看出,传统方法存在明显的SV波的串扰噪音,而且当ε<δ时,波场快照出现了不稳定的现象。而本发明的方法压制了SV串扰噪音,而且在ε<δ时仍然稳定。
接下来对HESS VTI介质模型进行测试。输入速度场(图3a)和各向异性参数场(3b、3c),得到的0.56s的波场快照图如图4b所示,与常规方法得到的波场快照进行对比可以看出,箭头所示的区域,常规方法存在明显的不稳定现象,而本发明的结果是稳定的。最后对BP TTI介质模型进行测试,输入速度场(图5a)、各向异性参数场(5b、5c)和倾角场(5d),得到0.56s的波场快照图如图6所示,传统方法的结果图6a中存在明显的不稳定现象,而本发明的结果图6b是稳定的。证明了本发明提出的一种稳定的各向异性TI介质正演模拟方法,能够解决目前的各向异性正演模拟技术存在SV波串扰和数值频散及当各向异性参数和倾斜横向各向同性(TTI)介质具有高度变化的倾角和方位角时存在不稳定性的问题。
地震正演模拟是逆时偏移和全波形反演的基础。准确的各向异性正演模拟能够更精确地模拟地下构造地震波的传播形态。随着科学技术的迅猛发展、大型计算机的广泛应用以及地震勘探仪器动态范围的不断扩大,简单、浅中层常规油气资源的开发利用日渐殆尽,勘探开发的重点逐步向构造复杂的油气藏以及岩性、裂缝性等非常规、隐蔽油气藏转移。复杂油气藏的勘探开发促进了地球科学领域中一系列新概念、新理论的产生,在对地球内部介质的结构与属性的研究中,人们对地球介质各向异性的认识逐步深化。基于传统耦合伪声波方程的正演模拟方法不仅存在SV波串扰和数值频散,而且当各向异性参数和倾斜横向各向同性(TTI)介质具有高度变化的倾角和方位角时存在不稳定性。
因此,考虑到上述各向异性正演模拟技术存在的缺陷,本发明提出了一种稳定的各向异性TI介质正演模拟方法,为研究复杂地下介质的波场传播规律与特征及反演成像提供基础。
当然,上述说明并非是对本发明的限制,本发明也并不仅限于上述举例,本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换,也应属于本发明的保护范围。
Claims (6)
1.稳定的各向异性TI介质正演模拟方法,其特征在于,该稳定的各向异性TI介质正演模拟方法包括:
步骤1,建立各向异性地质模型;
步骤2,建立VTI介质中的纯拟声波相速度方程;
步骤3,建立VTI和TTI介质的纯拟声波控制方程;
步骤4,分析拟声波相速度曲线的精度,选择近似阶数;
步骤5,求解选定阶数的控制方程,得到TI介质拟声波正演模拟的结果。
2.根据权利要求1所述的稳定的各向异性TI介质正演模拟方法,其特征在于,在步骤1中,结合地震和测井资料及其解释结果包括层位、断层,建立模型框架,之后加上各层的速度参数,建立各向异性速度场;根据测井资料或者岩石物理测试建立各向异性参数场;根据各向异性速度场,通过倾角估计建立倾角场。
3.根据权利要求1所述的稳定的各向异性TI介质正演模拟方法,其特征在于,在步骤2,基于精确的VTI介质相速度公式:
其中,vp 是 相位速度,θ是沿轴对称方向测量的相位角,vp0和vs0分别是拟声波和qSV波的速度;ε和δ是Thomsen参数;正号对应拟声波,负号对应qSV波公式;通过设置qSV波速度沿着对称轴方向等于0,也就是f≡1;则VTI介质拟声波相速度方程写为:
这里,引入了一个新的符号f(θ)表示方程(2)的平方根部分,并使得用二次范数中的最佳逼近来逼近方程(2)的平方根部分,f(θ)由下式求得:
二次范数最佳逼近的思想是通过寻找λ1,λ2,…,λn以满足如下方程:
四阶函数系统的逼近格式为:{sin4θ,sin2θ,1};越高阶函数系统将产生越高的近似精度,但它会牺牲更多的计算量和存储,特别是具有倾斜对称轴的横向各向同性TTI介质的情况;
通过二次范数最佳逼近,VTI介质拟声波相速度方程(2)可以近似为
方程(5)、(6)、(7)分别是二阶、四阶、六阶近似的VTI介质纯拟声波相速度方程。
4.根据权利要求3所述的稳定的各向异性TI介质正演模拟方法,其特征在于,在步骤3,将VTI介质的对称轴从垂直角度旋转到倾斜角度,得到TTI介质中的纯拟声波相速度方程;对于TTI介质,将对称轴从垂直旋转到倾斜角,根据方程(5)、(6)、(7)得到不同近似阶数的TTI介质拟声波相速度方程:
其中φ表示TTI介质的对称轴旋转角度;
把方程(6)和(9)与傅里叶域中的波场UP(kx,kz,t)相乘,再进行傅里叶逆变换,建立四阶近似下VTI介质和TTI介质的新的纯拟声波控制方程如下:
同样可建立其它近似精度的VTI和TTI介质控制方程;
其中,UP表示拟声波波场,kx、kz分别表示x和z方向的波数,t表示时间,φ表示TTI介质的对称轴旋转角度。
5.根据权利要求1所述的稳定的各向异性TI介质正演模拟方法,其特征在于,在步骤4中,根据各向异性介质的准确相速度方程以及不同阶数的近似相速度方程,绘制准确的以及不同近似阶数的拟声波相速度随相位角变化的曲线,并进行对比分析,选择合适的近似阶数。
6.根据权利要求1所述的稳定的各向异性TI介质正演模拟方法,其特征在于,在步骤5中,输入各向异性速度场、参数场和倾角场,求解选定阶数的控制方程,得到TI介质拟声波正演模拟的结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910343187.8A CN110261896B (zh) | 2019-04-26 | 2019-04-26 | 稳定的各向异性ti介质正演模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910343187.8A CN110261896B (zh) | 2019-04-26 | 2019-04-26 | 稳定的各向异性ti介质正演模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110261896A CN110261896A (zh) | 2019-09-20 |
CN110261896B true CN110261896B (zh) | 2021-07-20 |
Family
ID=67913908
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910343187.8A Active CN110261896B (zh) | 2019-04-26 | 2019-04-26 | 稳定的各向异性ti介质正演模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110261896B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112764105B (zh) * | 2020-10-16 | 2022-07-12 | 中国石油大学(华东) | Hti介质准纵波正演模拟方法、装置、存储介质及处理器 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104597484A (zh) * | 2013-10-31 | 2015-05-06 | 中国石油天然气集团公司 | 一种三维tti地震各向异性介质逆时偏移成像方法及装置 |
CN105044771A (zh) * | 2015-08-05 | 2015-11-11 | 北京多分量地震技术研究院 | 基于有限差分法的三维tti双相介质地震波场数值模拟方法 |
CN106842306A (zh) * | 2017-04-18 | 2017-06-13 | 中国科学院地质与地球物理研究所 | 一种全局优化的交错网格有限差分正演模拟方法和装置 |
CN107092029A (zh) * | 2017-04-26 | 2017-08-25 | 中国石油大学(北京) | 一种地震反演方法和装置 |
CN108303736A (zh) * | 2017-12-07 | 2018-07-20 | 东华理工大学 | 各向异性ti介质最短路径射线追踪正演方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9759797B2 (en) * | 2012-03-20 | 2017-09-12 | Oxford University Innovation Limited | Motion sensitized and motion suppressed imaging using dante prepared pulse trains |
CN104948174A (zh) * | 2015-06-30 | 2015-09-30 | 中国石油天然气集团公司 | 一种用于模拟声波测井的ti介质缩尺模型井及其制作方法 |
-
2019
- 2019-04-26 CN CN201910343187.8A patent/CN110261896B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104597484A (zh) * | 2013-10-31 | 2015-05-06 | 中国石油天然气集团公司 | 一种三维tti地震各向异性介质逆时偏移成像方法及装置 |
CN105044771A (zh) * | 2015-08-05 | 2015-11-11 | 北京多分量地震技术研究院 | 基于有限差分法的三维tti双相介质地震波场数值模拟方法 |
CN106842306A (zh) * | 2017-04-18 | 2017-06-13 | 中国科学院地质与地球物理研究所 | 一种全局优化的交错网格有限差分正演模拟方法和装置 |
CN107092029A (zh) * | 2017-04-26 | 2017-08-25 | 中国石油大学(北京) | 一种地震反演方法和装置 |
CN108303736A (zh) * | 2017-12-07 | 2018-07-20 | 东华理工大学 | 各向异性ti介质最短路径射线追踪正演方法 |
Non-Patent Citations (1)
Title |
---|
"TI介质地震波传播特征与正演方法研究";梁锴;《中国博士学位论文全文数据库基础科学辑》;20091215;第五、六章 * |
Also Published As
Publication number | Publication date |
---|---|
CN110261896A (zh) | 2019-09-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Min et al. | Improved frequency-domain elastic wave modeling using weighted-averaging difference operators | |
Yin et al. | Bayesian inversion for effective pore-fluid bulk modulus based on fluid-matrix decoupled amplitude variation with offset approximation | |
EP3028071B1 (en) | Method and device for the generation and application of anisotropic elastic parameters in horizontal transverse isotropic (hti) media | |
US8830788B2 (en) | Sensitivity kernal-based migration velocity analysis in 3D anisotropic media | |
Zong et al. | Broadband seismic amplitude variation with offset inversion | |
CN109557582B (zh) | 一种二维多分量地震资料偏移成像方法及系统 | |
CN105652321A (zh) | 一种粘声各向异性最小二乘逆时偏移成像方法 | |
CN111399044A (zh) | 一种储层渗透率预测方法、装置及存储介质 | |
Chen et al. | A k-space operator-based least-squares staggered-grid finite-difference method for modeling scalar wave propagation | |
Zhang et al. | Viabilities of seismic ray impedance and elastic impedance for hydrocarbon-sand discrimination | |
Mu et al. | Modeling viscoacoustic wave propagation using a new spatial variable-order fractional Laplacian wave equation | |
CN109946742A (zh) | 一种TTI介质中纯qP波地震数据模拟方法 | |
CN110737018B (zh) | Vsp地震资料各向异性建模方法 | |
CN110261896B (zh) | 稳定的各向异性ti介质正演模拟方法 | |
CN112462427B (zh) | 多分量地震资料保幅角度域共成像点道集提取方法及系统 | |
BRPI0415559B1 (pt) | Method of pre-stacking migration in kirchhoff time for ps waves | |
Zhao et al. | Combined inversion of first‐arrival travel times and reflection travel times | |
CN112684498A (zh) | 一种基于宽方位地震数据的储层裂缝预测方法及系统 | |
US12000971B2 (en) | Method and system for seismic processing using virtual trace bins based on offset attributes and azimuthal attributes | |
Anthony et al. | An optimized staggered-grid finite-difference operator for seismic wave simulation in poroelastic media | |
Yang et al. | AVA inversion for novel fluid indicator considering consolidation parameter | |
Li et al. | Joint elastic and petrophysical inversion using prestack seismic and well log data | |
Gibson Jr et al. | Modeling and velocity analysis with a wavefront-construction algorithm for anisotropic media | |
Wang et al. | Three-parameter prestack nonlinear inversion constrained by gradient structure similarity | |
Jin et al. | Improved generalized moveout approximation with a novel parameterization |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |