CN110673212B - 一种模型约束的薄层多波ava联合反演方法 - Google Patents
一种模型约束的薄层多波ava联合反演方法 Download PDFInfo
- Publication number
- CN110673212B CN110673212B CN201911021825.0A CN201911021825A CN110673212B CN 110673212 B CN110673212 B CN 110673212B CN 201911021825 A CN201911021825 A CN 201911021825A CN 110673212 B CN110673212 B CN 110673212B
- Authority
- CN
- China
- Prior art keywords
- thin
- layer
- ava
- model
- thin layer
- 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 27
- 230000001419 dependent effect Effects 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000013016 damping Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract 1
- 239000003245 coal Substances 0.000 description 7
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 208000035126 Facies Diseases 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- AYEKOFBPNLCAJY-UHFFFAOYSA-O thiamine pyrophosphate Chemical compound CC1=C(CCOP(O)(=O)OP(O)(O)=O)SC=[N+]1CC1=CN=C(C)N=C1N AYEKOFBPNLCAJY-UHFFFAOYSA-O 0.000 description 1
- 238000002834 transmittance 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
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
- G01V2210/6242—Elastic parameters, e.g. Young, Lamé or Poisson
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
本发明涉及一种模型约束的薄层多波AVA联合反演方法,包括:根据薄层模型m,计算薄层频变反射系数RPP(f)、RPS(f);求取频率域地震子波WPP(f)、WPS(f);计算频率域的薄层AVA道集SPP(f)、SPS(f);计算时间域的薄层AVA道集sPP(t)、sPS(t),并带入建立的薄层AVA联合反演的目标函数Q(m);判断Q(m)是否小于误差阀值或者迭代次数是否大于最大迭代次数;若否,利用LM算法,计算模型更新量Δm,并重新计算薄层频变反射系数RPP(f)、RPS(f),进行迭代计算,直至精度或者迭代次数满足要求。若是,输出反演结果m,薄层模型参数m包含薄层三层介质的纵波速度、横波速度、密度及薄层厚度。本发明充分地利用了薄层多波信息,相比常规地震反演方法,可以获得更为准确的薄层弹性参数及厚度反演结果。
Description
技术领域
本发明涉及地震勘探技术领域,尤其涉及一种模型约束的薄层多波AVA联合反演方法。
背景技术
随着我国油气、煤炭等资源开发的不断深入,岩性储层、碳酸盐岩缝洞型储层、非常规油气等复杂油气藏及煤储层成为勘探的重点。在我国东部地区或一些西部盆地,绝大多数的中、新生代陆相含油盆地大都以薄层砂、泥岩沉积为主,并夹有少量的碳酸盐岩、页岩及膏盐层,这些地层岩性和厚度横向变化较大,而且厚度远远低于目前地震技术的垂向分辨率,以薄层的形式存在。我国煤系地层煤层厚度多在2-10m之间,且常与砂泥岩互层以薄互层的形式存在,反射响应并非单一界面产生,而是顶底反射、层间多次波和转换波叠加而成的复合波。建立在单阻抗差界面、厚层反射基础之上的地震技术受到越来越多的挑战,存在不可接受的误差。因此,开展基于薄层反射理论的反演方法研究,对油气、煤炭等地震勘探开发来说有着重要的现实意义。
发明内容
本发明的目的在于,解决现有技术中存在的上述不足之处。
为实现上述目的,本发明实施例提供了一种模型约束的薄层多波AVA联合反演方法,包括以下步骤:S1,输入薄层模型m,利用薄层精准反透射系数方程,计算薄层频变反射系数RPP(f)、RPS(f)。S2,根据实际地震记录,可以获得时间域的地震子波wPP(t)、wPS(t)。S3,根据wPP(t)、wPS(t),利用FFT变换,求取频率域地震子波WPP(f)、WPS(f)。S4,计算频率域的薄层AVA道集SPP(f)、SPS(f),SPP(f)=RPP(f)·WPP(f)、SPS(f)=RPS(f)·WPS(f)。S5,利用FFT逆变换,计算时间域的薄层AVA道集sPP(t)、sPS(t)。S6,建立薄层AVA联合反演的目标函数Q(m),并将sPP(t)、sPS(t)代入目标函数Q(m)。S7,判断Q(m)是否小于误差阀值或者迭代次数是否大于最大迭代次数。S8,若Q(m)大于或者等于误差阀值,且迭代次数小于或者等于最大迭代次数,利用LM算法,计算模型更新量Δm。S9,更新模型参数m=m+Δm,并将更新后的薄层模型m返回S1,进行迭代计算,直至精度满足要求或者迭代次数大于最大迭代次数为止。S10,若Q(m)小于误差阀值,或者迭代次数大于最大迭代次数,输出反演结果m,薄层模型参数m包含薄层三层介质的纵波速度、横波速度、密度及薄层厚度。
在一种可能的实施方式中,薄层模型参数m包含薄层三层介质的纵波速度、横波速度、密度及薄层厚度。
在一种可能的实施方式中,在步骤S6中,目标函数Q(m)=η||dPP-sPP||2+(1-η)||dPS-sPS||2,其中,η为PP波的权重,1-η为PS波的权重。
在一种可能的实施方式中,在步骤S8中,利用LM算法,计算模型更新量Δm,满足如下公式:
Δm=[H(m)+λI]-1JT(m)[s(m0+Δm)-d(m0)];
其中,H(m)为海森矩阵,且有H(m)=JT(m)·J(m),λ为阻尼系数,I为单位矩阵,s为拟合地震AVA道集,代表sPP或sPS,d为实际的AVA道集,代表dPP或dPS,J为雅各比矩阵,且有:
在一种可能的实施方式中,在步骤S3中,求取频率域地震子波WPP(f)、WPS(f)满足如下公式包括:
WPP(f)=FFT[wPP(t)];
WPS(f)=FFT[wPS(t)]。
在一种可能的实施方式中,在步骤S5中,计算时间域的薄层AVA道集sPP(t)、sPS(t)满足如下公式:
sPP(t)=FFT-1[SPP(f)];
sPS(t)=FFT-1[SPS(f)]。
本发明充分地利用了薄层多波信息,相比常规地震反演方法,可以获得更为准确的薄层弹性参数及厚度反演结果。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1为本发明实施例提供一种薄层多波AVA联合反演方法的流程示意图;
图2为本发明实施例提供一种薄层多波AVA联合反演流程图;
图3为本发明实施例提供的模型1的AVA道集;
图4为本发明实施例提供的模型1的反演结果;
图5为本发明实施例提供的模型2的AVA道集;
图6为本发明实施例提供的模型2的反演结果;
图7为本发明实施例提供的模型3的AVA道集;
图8为本发明实施例提供的模型3的反演结果;
图9为本发明实施例提供的模型4的AVA道集;
图10为本发明实施例提供的模型4的反演结果。
具体实施方式
本发明的说明书实施例和权利要求书及附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元。方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
下面结合附图和实施例,对本发明的技术方案做进一步的详细描述。
如图1-2,本发明实施例提供一种模型约束的薄层多波AVA联合反演方法,包括以下步骤:
S1,输入薄层模型m,利用薄层精准反透射系数方程,计算薄层频变反射系数RPP(f)、RPS(f)。
S2,根据实际地震记录,可以获得时间域的地震子波wPP(t)、wPS(t)。
S3,根据wPP(t)、wPS(t),利用FFT变换,求取频率域地震子波WPP(f)、WPS(f)。
S4,计算频率域的薄层AVA道集SPP(f)、SPS(f),SPP(f)=RPP(f)·WPP(f)、SPS(f)=RPS(f)·WPS(f)。
S5,利用FFT逆变换,计算时间域的薄层AVA道集sPP(t)、sPS(t)。
S6,建立薄层AVA联合反演的目标函数Q(m),并将sPP(t)、sPS(t)代入目标函数Q(m)。
S7,判断Q(m)是否小于误差阀值或者迭代次数是否大于最大迭代次数。
S8,若Q(m)大于或者等于误差阀值,且迭代次数小于或者等于最大迭代次数,利用LM算法,计算模型更新量Δm。
S9,更新模型参数m=m+Δm,并将更新后的薄层模型m返回S1,进行迭代计算,直至精度满足要求或者迭代次数大于最大迭代次数为止。
S10,若Q(m)小于误差阀值,或者迭代次数大于最大迭代次数,输出反演结果m,薄层模型参数m包含薄层三层介质的纵波速度、横波速度、密度及薄层厚度。
下面结合具体的示例进一步解释。
S1,输入薄层模型m,在已知薄层模型m的情况下,可以根据薄层精准反透射系数方程(Yang et al.,2016)计算薄层PP波及PS波频变反射系数RPP(f)、RPS(f):
M[RPP RPS TPP TPS]T=N, (1)
其中M=[mi1 mi2 mi3 mi4],N=[ni],i=1,2,3,4,且有:
[mi2]=A[-cosβ1 -sinβ1 -jρ1ωvS1sin2β1 -jρ1ωvS1 cos2β1]T, (2b)
[mi4]=[cosβ3 -sinβ3 jρ3ωvS3 sin2β3 -jρ3ωvS3cos2β3]T, (2d)
a11=a44=2sin2β2cosP+cos2β2cosQ,
a12=a34=-j(tanα2cos2β2sinP-sin2β2sinQ),
a22=a33=cos2β2cosP+2sin2β2cosQ,
a31=a42=2jρ2ωvS2sinβ2cos2β2(cosQ-cosP),
S2,根据实际地震记录,获得时间域的地震子波wPP(t)、wPS(t)。
S3,利用FFT变换,求取频率域地震子波WPP(f)、WPS(f):
WPP(f)=FFT[wPP(t)], (4a)
WPS(f)=FFT[wPS(t)]。 (4b)
S4,在频率域利用反射系数与地震子波乘积获得薄层AVA道集SPP(f)、SPS(f),即:
SPP(f)=RPP(f)·WPP(f), (5a)
SPS(f)=RPS(f)·WPS(f)。 (5b)
S5,通过FFT逆变换,可以获得薄层模型m的时间域AVA道集sPP(t)、sPS(t);
sPP(t)=FFT-1[SPP(f)], (6a)
sPS(t)=FFT-1[SPS(f)], (6b)
S6,建立薄层AVA联合反演的目标函数,并在已知实际薄层PP波及PS波AVA道集dPP(t)、dPS(t)的情况下,计算目标函数Q(m)=η||dPP-sPP||2+(1-η)||dPS-sPS||2。
其中,η为PP波的权重,1-η为PS波的权重。
S7,设定误差阀值ε及最大迭代次数,判断Q(m)是否小于ε或者迭代次数大于最大迭代次数。
S8,若Q(m)大于或者等于误差阀值,且迭代次数小于或者等于最大迭代次数,
则利用LM算法(Levenberg,1944;Marquardt,1963)计算模型更新量Δm:
Δm=[H(m)+λI]-1JT(m)[s(m0+Δm)-d(m0)]。 (7)
其中:H(m)为海森矩阵,且有H(m)=JT(m)·J(m);λ为阻尼系数;I为单位矩阵;s为拟合地震AVA道集,代表sPP或sPS;d为实际的AVA道集,代表dPP或dPS;J为雅各比矩阵,且有:
地震记录关于模型参数的一阶导数可以利用差分近似求解,如PP波地震记录对第一层纵波速度的偏导数为:
S9,更新模型参数m=m+Δm,并将更新后的薄层模型m返回S1,进行迭代计算。
S10,若Q(m)小于误差阀值,或者迭代次数大于最大迭代次数,则输出反演结果m,薄层模型参数m包含薄层三层介质的纵波速度、横波速度、密度及薄层厚度。
反演测试
为了测试该发明给出的薄层多波AVA联合反演方法的有效性,建立如表1所示的四个典型的薄层模型,包括:高阻抗薄层、低阻抗薄层、阻抗递增薄层和阻抗递减薄层。设定薄层厚度为10m,30Hz主频的雷克子波作为地震子波,入射角由5~25度,AVA道集时间采样间隔为1ms,采样长度为200ms。模型1~4的合成AVA道集如图3、5、7、9所示,其中PS波剖面已压缩到PP波剖面时刻上。利用本发明给出基于薄层反射理论的多波AVA联合反演方法对以上拟合的AVA道集进行反演,反演结果如图4、6、8、10所示。同时,为了展示本发明反演方法的优势,我们给出基于Zoeppritz方程的多波AVA反演结果。反演过程中采取了相同的模型硬约束条件,反演结果显示:相比基于Zoepprtiz方程的多波AVA反演方法,本发明给出的薄层多波AVA联合反演方法可以更精准的刻画目标薄层,获得更为准确的薄层弹性参数及厚度反演结果。
表1薄层模型的弹性参数,其中速度单位为km/s,密度单位为g/cm3。
以上的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种模型约束的薄层多波AVA联合反演方法,其特征在于,包括以下步骤:
S1,输入薄层模型m,利用薄层精准反透射系数方程,计算薄层频变反射系数RPP(f)、RPS(f);
S2,根据实际地震记录,可以获得时间域的地震子波wPP(t)、wPS(t);
S3,根据wPP(t)、wPS(t),利用FFT变换,求取频率域地震子波WPP(f)、WPS(f);
S4,计算频率域的薄层AVA道集SPP(f)、SPS(f),SPP(f)=RPP(f)·WPP(f)、SPS(f)=RPS(f)·WPS(f);
S5,利用FFT逆变换,计算时间域的薄层AVA道集sPP(t)、sPS(t);
S6,建立薄层AVA联合反演的目标函数Q(m),并将sPP(t)、sPS(t)代入目标函数Q(m);
S7,判断Q(m)是否小于误差阀值或者迭代次数是否大于最大迭代次数;
S8,若Q(m)大于或者等于误差阀值,且迭代次数小于或者等于最大迭代次数,利用LM算法,计算模型更新量Δm;
S9,更新薄层模型m=m+Δm,并将更新后的薄层模型m返回S1,进行迭代计算,直至精度满足要求或者迭代次数大于最大迭代次数为止;
S10,若Q(m)小于误差阀值,或者迭代次数大于最大迭代次数,输出反演结果m,薄层模型m包含薄层三层介质的纵波速度、横波速度、密度及薄层厚度。
2.根据权利要求1所述的模型约束的薄层多波AVA联合反演方法,其特征在于,在步骤S6中,所述目标函数Q(m)=η||dPP-sPP||2+(1-η)||dPS-sPS||2
其中,η为PP波的权重,1-η为PS波的权重。
4.根据权利要求1所述的模型约束的薄层多波AVA联合反演方法,其特征在于,在步骤S3中,所述求取频率域地震子波WPP(f)、WPS(f)满足如下公式包括:
WPP(f)=FFT[wPP(t)];
WPS(f)=FFT[wPS(t)]。
5.根据权利要求1所述的模型约束的薄层多波AVA联合反演方法,其特征在于,在步骤S5中,所述计算时间域的薄层AVA道集sPP(t)、sPS(t)满足如下公式:
sPP(t)=FFT-1[SPP(f)];
sPS(t)=FFT-1[SPS(f)]。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911021825.0A CN110673212B (zh) | 2019-10-25 | 2019-10-25 | 一种模型约束的薄层多波ava联合反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911021825.0A CN110673212B (zh) | 2019-10-25 | 2019-10-25 | 一种模型约束的薄层多波ava联合反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110673212A CN110673212A (zh) | 2020-01-10 |
CN110673212B true CN110673212B (zh) | 2021-06-18 |
Family
ID=69084400
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911021825.0A Active CN110673212B (zh) | 2019-10-25 | 2019-10-25 | 一种模型约束的薄层多波ava联合反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110673212B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102508293A (zh) * | 2011-11-28 | 2012-06-20 | 中国石油大学(北京) | 一种叠前反演的薄层含油气性识别方法 |
CN103984018A (zh) * | 2014-06-06 | 2014-08-13 | 中国地质大学(北京) | 多波联合振幅随入射角变化反演方法 |
WO2014209728A1 (en) * | 2013-06-25 | 2014-12-31 | Westerngeco Llc | Seismic wavefield deghosting and noise attenuation |
CN104614763A (zh) * | 2015-01-19 | 2015-05-13 | 中国石油大学(北京) | 基于反射率法的多波avo储层弹性参数反演方法及系统 |
CN107688197A (zh) * | 2017-07-12 | 2018-02-13 | 中国石油天然气集团公司 | 薄层预测方法和装置 |
CN107861153A (zh) * | 2017-10-25 | 2018-03-30 | 中国地质大学(北京) | 薄层pp波反射系数的快速求解方法 |
CN109188511A (zh) * | 2018-08-27 | 2019-01-11 | 中国地质大学(北京) | 一种砂泥岩薄互层介质多波avo联合反演方法 |
CN110187389A (zh) * | 2019-06-24 | 2019-08-30 | 中国地质大学(北京) | 一种基于薄层反射理论的ava反演方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8073624B2 (en) * | 2009-03-12 | 2011-12-06 | Schlumberger Technology Corporation | Technique and system to process a seismic signal to determine an attribute of a geological layer |
WO2018183346A1 (en) * | 2017-03-27 | 2018-10-04 | Ion Geophysical Corporation | Amplitude compensation of reverse time migration (rtm) gathers for avo/ava analysis |
-
2019
- 2019-10-25 CN CN201911021825.0A patent/CN110673212B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102508293A (zh) * | 2011-11-28 | 2012-06-20 | 中国石油大学(北京) | 一种叠前反演的薄层含油气性识别方法 |
WO2014209728A1 (en) * | 2013-06-25 | 2014-12-31 | Westerngeco Llc | Seismic wavefield deghosting and noise attenuation |
CN103984018A (zh) * | 2014-06-06 | 2014-08-13 | 中国地质大学(北京) | 多波联合振幅随入射角变化反演方法 |
CN104614763A (zh) * | 2015-01-19 | 2015-05-13 | 中国石油大学(北京) | 基于反射率法的多波avo储层弹性参数反演方法及系统 |
CN107688197A (zh) * | 2017-07-12 | 2018-02-13 | 中国石油天然气集团公司 | 薄层预测方法和装置 |
CN107861153A (zh) * | 2017-10-25 | 2018-03-30 | 中国地质大学(北京) | 薄层pp波反射系数的快速求解方法 |
CN109188511A (zh) * | 2018-08-27 | 2019-01-11 | 中国地质大学(北京) | 一种砂泥岩薄互层介质多波avo联合反演方法 |
CN110187389A (zh) * | 2019-06-24 | 2019-08-30 | 中国地质大学(北京) | 一种基于薄层反射理论的ava反演方法 |
Non-Patent Citations (2)
Title |
---|
"两种薄层弹性参数多波联合反演的对比";成联正 等;《煤田地质与勘探》;20140228;第42卷(第1期);第68-71页 * |
"基于低频软约束的叠前AVA稀疏层反演";张丰麒 等;《石油地球物理勘探》;20170831;第52卷(第4期);第770-782页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110673212A (zh) | 2020-01-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yu et al. | Determining crustal structure beneath seismic stations overlying a low‐velocity sedimentary layer using receiver functions | |
KR101861060B1 (ko) | 오프셋 의존 탄성 fwi를 통한 다중 파라미터 역산 | |
US8990053B2 (en) | Method of wavelet estimation and multiple prediction in full wavefield inversion | |
CN107783187B (zh) | 一种将测井速度和地震速度结合建立三维速度场的方法 | |
AU2005200918A1 (en) | Horizon-based residual depth migration velocity analysis | |
EP3028071A2 (en) | Method and device for the generation and application of anisotropic elastic parameters in horizontal transverse isotropic (hti) media | |
CN101634718A (zh) | 相对波阻抗约束条件下多井测井参数反演的储层预测方法 | |
CN104570067A (zh) | 一种地球物理勘探中相控地震反演方法 | |
CN110780351B (zh) | 纵波和转换波叠前联合反演方法及系统 | |
CN111522063B (zh) | 基于分频联合反演的叠前高分辨率流体因子反演方法 | |
CN107861153A (zh) | 薄层pp波反射系数的快速求解方法 | |
CN113740901A (zh) | 基于复杂起伏地表的陆上地震数据全波形反演方法及装置 | |
Zhang et al. | Static corrections in mountainous areas using Fresnel-wavepath tomography | |
CN110187389B (zh) | 一种基于薄层反射理论的ava反演方法 | |
Coléou et al. | Petrophysical seismic inversion | |
CN110673212B (zh) | 一种模型约束的薄层多波ava联合反演方法 | |
CN111474576B (zh) | 一种保持地层结构的叠前地震道集反演初始模型构建方法 | |
EP3929630B1 (en) | Computer implemented method for correcting a reservoir model of a reservoir geological formation based on seismic images | |
CN106842291B (zh) | 一种基于叠前地震射线阻抗反演的不整合圈闭储层岩性预测方法 | |
CN109143335B (zh) | 一种合成地震记录的制作方法、系统、介质及设备 | |
CN111077578B (zh) | 岩层分布预测方法和装置 | |
CN115373022A (zh) | 一种基于振幅相位校正的弹性波场Helmholtz分解方法 | |
CN115906559A (zh) | 一种基于混合网格的大地电磁自适应有限元正演方法 | |
Kontakis et al. | Efficient 1.5 D full waveform inversion in the Laplace-Fourier domain | |
CN111308550B (zh) | 一种页岩vti储层的各向异性参数多波联合反演方法 |
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 |