CN101329405B - 一种简单的多参数地震反演方法 - Google Patents
一种简单的多参数地震反演方法 Download PDFInfo
- Publication number
- CN101329405B CN101329405B CN200710119093XA CN200710119093A CN101329405B CN 101329405 B CN101329405 B CN 101329405B CN 200710119093X A CN200710119093X A CN 200710119093XA CN 200710119093 A CN200710119093 A CN 200710119093A CN 101329405 B CN101329405 B CN 101329405B
- Authority
- CN
- China
- Prior art keywords
- impedance
- seismic
- wave
- data
- elastic
- 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
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明是一种借助于地震波弹性阻抗来提高常规叠前地震反演的精度和可靠性,进而同时提取多种地层参数的简单的多参数地震反演方法。步骤是:利用全波测井资料;确定声波阻抗AI和弹性波阻抗EI;将CMP道集变换成角道集并进行部分角度叠加;利用常规叠前弹性阻抗反演流程分别对小角度叠加数据及中等入射角的部分叠加数据进行过井测线的AI和EI反演;对反演得到的EI和AI数据进行道运算,得到弹性密度比和层岩性变化和孔隙流体成分密切相关的地震属性;确定储层岩性及孔隙流体性质。本发明仅需要近、中角度的地震数据,其它参数是由弹性密度比与声波阻抗之间的关系通过道运算获得,提高了计算效率,可用于绝大多数采集条件获得的地震数据的叠前反演。
Description
技术领域
本发明涉及石油地球物理勘探技术,是一种借助于地震波弹性阻抗来提高常规叠前地震反演的精度和可靠性,进而同时提取多种地层参数的简单的多参数地震反演方法。
背景技术
地震勘探是通过人工激发地震波,在地表以一定的方式记录地震波的传播行为,并研究它们在地层中的传播规律,以查明地下的地质构造和岩性特征,进而寻找油气田的地球物理勘探方法。
地震勘探的最终目的是综合利用地震波的运动学和动力学特征来尽可能的获取地层全面的信息以及与油气藏的构造部位有关的岩性信息。早先,由于油气勘探和开发还处于大规模的勘探阶段,利用地震波在地球内部传播的旅行时信息就可以基本圈定构造的形态和位置。但是,随着勘探精度和难度的不断提高,常规的地震成像理论往往不能满足各种不同的地层条件,水平叠加技术也掩盖了反射振幅随偏移距变化所反映的重要信息,仅利用地震波的运动学特征很难准确地实现油气储层预测。因此,从地震资料中提取介质的弹性参数,并将这些参数与岩性和流体成分联系起来,在油藏描述和油气监测中发挥着重要作用。目前,利用地震信息提取岩性参数主要是通过地震资料反演来实现的,这主要包括常规叠后波阻抗反演、叠前振幅随入射角变化(AVO)反演、弹性波阻抗反演等。
1.波阻抗反演
波阻抗反演是从反射振幅求取波阻抗(速度和密度的乘积),进而得到波阻抗剖面或速度剖面。假定地震剖面上的地震道是法向入射道,即入射射线与岩层分界面垂直,其法向入射的反射系数为:
其中,ri为地层的反射系数,ρ为地层密度,v为地层的纵波速度,已知反射系数后,可由(1)式逐层递推出地下各层的波阻抗:
当地层密度变化不大时,可近似看成常数,这样便可直接得到速度剖面,当密度变化较大时,也可以利用Gardner经验关系式ρ=0.31v0.25进一步求得速度和密度剖面。
由波阻抗反演的原理可知,叠后波阻抗反演方法是以地震波垂直入射为假设前提的,只有当地下有波阻抗差异、且共中心点道集的炮检距变化范围较小或者反演的目的层较深时才能够取得比较好的地质效果。然而,实际的叠后地震资料并非是自激自收的地震记录,而是共反射点道集反射振幅叠加平均的结果,尤其在炮检距变化范围较大时,垂直入射假设的反射系数就与实际叠加道集的反射系数有很大不同,不能反映真实的叠加振幅信息。
2.AVO反演
AVO技术是利用叠前振幅随入射角变化关系来估算岩石弹性参量的方法。当入射角不太大且地层界面两侧岩性变化差异不大时,描述平面波反射和透射的Zoeppritz方程可以近似为:
R=Af1(θ)+Bf2(θ)+Cf3(θ)+L (3)
其中,R为实际地震记录的反射振幅,A,B,C为要反演的地层参数或AVO属性,fi(θ)为与角度有关的权系数。理论上,对不同入射角的角度道集进行曲线拟合就可以得到多种AVO属性。
由于AVO分析假设地震记录是一次反射波振幅与地震子波的褶积。但是,实际地震记录同时记录了模式转换波、多次反射波及各种噪音,随着偏移距(入射角大于25°)的增大(Mallick,2006),用(3)描述实际地震记录的振幅将会引起很大的误差。另外,由于地震子波是时变和空变的,而且在处理过程中由于正常时差较正(NMO)拉伸畸变也会使得近、远道的频带、子波形状不一致,因此,利用(3)式同时估计3个属性是十分困难的,特别是当资料的信噪比较低时,反演出高阶项属性将会完全扭曲真正的岩性信息所反映的事实。目前,常规的AVO分析方法是提取截距和梯度属性,然而,由于截距和梯度两个属性使得纵波速度、横波速度和地层密度三个弹性参数捆绑在一起而无法分离,只能通过各种经验关系得到一些地震属性剖面。
3.弹性波阻抗反演
设函数EI具有和波阻抗相似的性质,当平面波以任一角度入射时,反射系数可以近似表示成与垂直入射时相似的波阻抗形式:
R(θ)=(EI2-EI1)/(EI2+EI1) (4)
结合Zoeppritz方程的Aki&Richards(1980)近似公式,Connolly(1999)给出了第一个PP波弹性波阻抗公式:
上式中,R为实际地震记录的反射振幅,EI为弹性波阻抗,α、β、ρ分别为地层的纵、横波速度和密度, θ为P波入射角。
从理论上来看,利用声波阻抗(AI)和大角度的弹性阻抗(EI)不仅可以定性区分孔隙流体和岩性,而且还可以由一些特殊角度的弹性阻抗来构建孔隙度、渗透率、泥质含量等储层物性参数;从反演的实现过程来看,弹性阻抗反演是通过对近、中、远等多个部分角度叠加剖面来联合反演出多个地层参数,一般需要3—5个部分角度叠加剖面来稳定反演的结果,而且同时考虑了子波的空变特征,这比常规AVO反演有着更好的应用前景;从表达形式上看,非零偏移距的阻抗表达形式同时包含了纵、横波速度、密度等信息,而且对应于不同表述形式的反射系数近似公式,弹性波阻抗表达式也具有不同的形式。因此,弹性波阻抗反演能较常规地震道反演获得更多、更可靠的流体、孔隙度和泥岩含量等信息,有助于解释常规地震道反演和道积分剖面中的假象,降低反演的多解性,提高储层预测的精度。
但是,由于弹性波阻抗是对反射系数近似公式形式的改造,与AVO反演相类似,为了可靠的同时获得多个岩性参数,必须要参与反演的角道集尽可能的包含大角度(远道)信息。首先,就目前的采集标准,一般很难得到中、深层大角度的地震资料,而浅层的大角度信息受NMO拉伸影响很严重,实际上可用于反演的资料往往包含极为有限的远道信息;其次,对较为复杂的地层,入射角很可能比较大甚至会超过临界角,这时反射系数近似公式已经不再成立,再用弹性阻抗来描述反射振幅的变化就会带来更大的误差;再次,与AVO技术相类似,弹性波阻抗技术也是用一次反射波振幅与地震子波的褶积来描述实际地震振幅的变化特征,但当入射角比较大时,地震记录上远偏移距的模式转换波、多次反射波与有效波互相干涉,使得地震记录上的实际振幅与弹性波阻抗描述的振幅具有很大差异,另外,随着偏移距的增大,入射角估计引起的误差也越来越大(远道误差可达7—10°,Barens,2006),加上由于资料处理引起的振幅也越来越不保真,而且远道地震资料的信噪比较低,反演中过多的应用远道信息反而会降低反演结果的可靠性(Mallick,2001,2006)。
发明内容
本发明目的是提供一种仅需要近、中角度的地震数据,反演结果可靠稳定的简单的多参数地震反演方法。
本发明提供以下技术方案:
本发明具体步骤包括:
(1)利用全波列测井进行纵、横波、密度测井记录测井资料;
所述的测井资料利用已有的常规声波、密度测井得到声波和密度测井资料。
所述的测井资料中的横波记录,可以通过全波列测井得到或由已知岩石物理学关系确定。
(2)利用密度、纵波和横波资料确定声波阻抗AI和弹性波阻抗EI,用于步骤(4)中地震数据反演的低频模型约束;
所述的声波阻抗AI和弹性波阻抗EI通过下面公式得到:
AI=ρα (6)
上式中,α为纵波速度,β为横波速度,ρ为密度,θ为纵波入射角,同步骤(3)中中等角度叠加剖面的中心角。
(3)激发地震波并记录地震波,按常规地震资料处理流程对记录的地震资料进行振幅保持处理,形成供振幅随偏移距变化分析的正常时差校正后(NMO)的共中心点(CMP)道集;然后按常规叠前弹性反演数据处理方式将CMP道集变换成角道集并进行部分角度叠加,首先,将角道集中入射角较小的部分道集叠加形成用于步骤(4)中AI反演的小角度叠加数据,其次,选择以入射角θ为中心角,将中等入射角的角道集叠加得到用于步骤(4)中EI反演的中等角度叠加数据;
(4)以步骤(2)中AI和EI为低频约束,利用常规叠前弹性阻抗反演流程分别对小角度叠加数据及中等入射角的部分叠加数据进行过井测线的AI和EI反演;
(5)采用下式对步骤(4)中反演得到的EI和AI数据进行道运算,得到弹性密度比Rd=EI/AI (8)
(6)将Rd归一化到纵、横波速度比的尺度上或利用下式得到速度比参数,
上式中,α为纵波速度,β为横波速度,Rd为弹性密度比,θ为中等入射角度叠加剖面的中心角;
(7)利用步骤(5)中得到的横、纵波速度比β/α和AI,用以下公式得到与储层岩性变化和孔隙流体成分密切相关的地震属性;
SI=βρ=(β/α)αρ (10.1)
λρ=(αρ)2-2(βρ)2=AI2-2SI2 (10.3)
μρ=(βρ)2=SI2 (10.4)
上式中SI为横波阻抗,σ泊松比,ρ为密度,λ和μ分别为拉梅常数和剪切模量;
(8)确定储层岩性及孔隙流体性质。
所述的确定储层岩性及孔隙流体性质是:当弹性密度比Rd较低时,地层岩石为孔隙中含气或地层为孔隙度较大的疏松砂岩;
所述的确定储层岩性及孔隙流体性质是:当弹性密度比Rd较高时,地层为孔隙度较小的泥岩或孔隙中含水;
所述的确定储层岩性及孔隙流体性质是:利用横波阻抗SI、泊松比σ、弹性参数λρ、μρ、λ/μ地层参数进行拉梅常数油气检测。
本发明有以下效果:
仅需要近、中角度的地震数据,可用于绝大多数采集条件获得的地震数据的叠前反演。
仅对近、中角度的地震数据进行叠后反演,其它参数是由弹性密度比与声波阻抗之间的关系通过道运算获得,提高了计算效率。
可以同时估计多种地层参数,利用估计的弹性密度比不仅可以得到地下岩性和流体的变化特征,而且还提供了多分量地震联合反演、解释的连接桥梁。
不需要大角度(远道)的地震叠加数据,可以:
减小由于射角较大时反射系数近似公式不成立导致的反演结果的不确定性;
降低由于入射角估计的偏差、NMO拉伸、处理中振幅不保真、远道信噪比低等引起的反演结果的不可靠性;
改善由于远道模式转换波、多次反射波及有效波相互干涉造成反演结果的不稳定性。
附图说明
图1为实际纵、横波及密度测井曲线;
图2.1为声波阻抗与不同角度弹性阻抗的交会图;
图2.2为声波阻抗与不同角度弹性阻抗构建的弹性密度比的交会图;
图3.1为不同角度弹性阻抗构建的弹性密度比与实际泊松比的对比;
图3.2为不同角度弹性阻抗构建的弹性密度比与实际纵横波速度比的对比;
图4为用10°的弹性阻抗构建的弹性密度比与声波阻抗提取的多种地层参数与理论值的对比。
具体实施方式
本发明实施例具体实施方式为:
(1)利用测井所得的纵、横波、密度测井资料计算声波阻抗和中等入射角的弹性阻抗;
(2)同常规叠前反演类似,对角道集数据进行部分角度叠加;
(3)分别对近、中角度叠加剖面实现过井声阻抗和弹性阻抗反演;
(4)利用(8)式计算弹性密度比,利用公式(9)计算速度比,或直接将弹性密度比归一化到纵横波速度比的尺度上;
(5)利用(4)中所得的速度比参数与(3)中得到的声波阻抗通过道运算得到多种地层参数;
(6)利用(5)中所得的参数进行油气检测。
上述步骤还可得出地震属性进行综合储层预测,首先,弹性密度比Rd和纵、横波速度比成线性关系。由岩石物理学可知,速度比与地层岩性变化及孔隙流体成分密切相关,当孔隙中含气或地层为孔隙度较大的疏松砂岩时速度比比较低,而地层为孔隙度较小的泥岩或空隙中含水时则具有较高的速度比,因此弹性密度比Rd就直接反映了地层岩性变化和孔隙流体成分。
其次,λ主要与孔隙流体性质有关,而μ与孔隙流体的性质无关仅与岩石骨架有关,多种地层参数可以进行拉梅常数油气检测(λρ-μρ-λ/μ,简称LMR,Goodway(1997))。
从实施例图1可见,声波阻抗与弹性阻抗的交会图很难区分岩性和孔隙流体的变化,而弹性密度比与声波阻抗的交会图则很容易反映这一变化,即使在小角度时也能完全区分。
图3中实线为实际地层参数,虚线为估计值,可以看出,弹性密度比与泊松比及速度比吻合的非常好,相关系数达到0.99。
图4中实线为实际地层参数,虚线为估计值,很明显,用本发明方法估计的弹性参数与对应的实际地层参数匹配的很好。
Claims (6)
1.一种简单的多参数地震反演方法,其特征在于步骤包括:
(1)利用全波列测井进行纵、横波、密度测井记录测井资料;
(2)利用所述测井资料确定声波阻抗AI和弹性波阻抗EI,用于步骤(4)中地震数据反演的低频模型约束;
(3)激发地震波并记录地震波,按常规地震资料处理流程对记录的地震资料进行振幅保持处理,形成供振幅随偏移距变化分析的正常时差校正后的共中心点道集;然后按常规叠前弹性反演数据处理方式将共中心点道集变换成角道集并进行部分角度叠加,首先,将角道集中入射角较小的部分道集叠加形成用于步骤(4)中声波阻抗AI反演的小角度叠加数据,其次,选择位于中心的入射角θ,将中等入射角的角道集叠加得到用于步骤(4)中弹性波阻抗EI反演的中等入射角叠加数据;
(4)以步骤(2)中声波阻抗AI和弹性波阻抗EI为低频约束,利用常规叠前弹性阻抗反演流程分别对小角度叠加数据及中等入射角的部分叠加数据进行过井测线的声波阻抗AI和弹性波阻抗EI反演;
(5)采用下式对步骤(4)中反演得到的声波阻抗AI和弹性波阻抗EI数据进行道运算,得到弹性密度比:
Rd=EI/AI; (8)
(6)将弹性密度比Rd归一化到纵、横波速度比的尺度上或利用下式得到速度比参数,
上式中,α为纵波速度,β为横波速度,Rd为弹性密度比,θ为入射角;
(7)利用步骤(6)中得到的横、纵波速度比β/α和弹性波阻抗EI,用以下公式得到与储层岩性变化和孔隙流体成分密切相关的地震属性;
SI=βρ=(β/α)αρ (10.1)
λρ=(αρ)2-2(βρ)2=AI2-2SI2 (10.3)
μρ=(βρ)2=SI2 (10.4)
上式中SI为横波阻抗,σ泊松比,ρ为密度,λ和μ分别为拉梅常数和剪切模量;
(8)确定储层岩性及孔隙流体性质。
2.根据权利要求1所述的简单的多参数地震反演方法,其特征在于:步骤(1)所述的测井资料利用已有的常规声波、密度测井得到声波和密度测井资料。
3.根据权利要求1所述的简单的多参数地震反演方法,其特征在于:步骤(1)所述的测井资料中的横波记录,通过全波列测井得到或由已知岩石物理学关系确定。
4.根据权利要求1所述的简单的多参数地震反演方法,其特征在于:步骤(8)所述的确定储层岩性及孔隙流体性质是:当弹性密度比Rd较低时,地层岩石为孔隙中含气或地层为孔隙度较大的疏松砂岩。
5.根据权利要求1所述的简单的多参数地震反演方法,其特征在于:步骤(8)所述的确定储层岩性及孔隙流体性质是:当弹性密度比Rd较高时,地层为孔隙度较小的泥岩或孔隙中含水。
6.根据权利要求1所述的简单的多参数地震反演方法,其特征在于:步骤(8)所述的确定储层岩性及孔隙流体性质是:利用横波阻抗SI、泊松比σ、弹性参数λρ、μρ、λ/μ地层参数进行拉梅常数油气检测。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200710119093XA CN101329405B (zh) | 2007-06-20 | 2007-06-20 | 一种简单的多参数地震反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200710119093XA CN101329405B (zh) | 2007-06-20 | 2007-06-20 | 一种简单的多参数地震反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101329405A CN101329405A (zh) | 2008-12-24 |
CN101329405B true CN101329405B (zh) | 2011-02-09 |
Family
ID=40205306
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200710119093XA Active CN101329405B (zh) | 2007-06-20 | 2007-06-20 | 一种简单的多参数地震反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101329405B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2615591C1 (ru) * | 2013-05-24 | 2017-04-05 | Эксонмобил Апстрим Рисерч Компани | Многопараметрическая инверсия через зависящую от сдвига упругую полноволновую инверсию (fwi) |
Families Citing this family (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101644781B (zh) * | 2009-07-28 | 2011-08-17 | 刘学伟 | 一种利用纵横波波阻抗增量比识别天然气水合物的方法 |
CN102109616B (zh) * | 2009-12-23 | 2012-08-15 | 中国石油天然气集团公司 | 一种沉积旋回约束的波阻抗反演方法 |
CN101907729B (zh) * | 2010-06-11 | 2012-09-05 | 北京诺克斯达石油科技有限公司 | 射线弹性参数反演方法 |
CN102455437B (zh) * | 2010-10-14 | 2014-06-18 | 中国石油天然气集团公司 | 一种确定储层和流体的方法 |
CN102455272A (zh) * | 2010-10-15 | 2012-05-16 | 中国石油化工股份有限公司 | 确定流体密度及确定固体骨架密度的方法 |
CN102454399B (zh) * | 2010-10-26 | 2015-01-14 | 中国石油化工股份有限公司 | 测井声波时差信号校正方法 |
CN102129086A (zh) * | 2010-12-10 | 2011-07-20 | 中国石油天然气集团公司 | 用pp、ps波叠前数据联合反演纵、横波阻抗的方法 |
CN102147478B (zh) * | 2010-12-29 | 2013-06-26 | 中国海洋大学 | 复杂油藏的叠前低频信号识别方法 |
CN102169190B (zh) * | 2011-01-06 | 2012-10-24 | 中国科学院地质与地球物理研究所 | 一种调制补充子空间的井约束叠前弹性参数反演方法 |
CN102288993B (zh) * | 2011-05-16 | 2013-03-20 | 中国石油大学(北京) | 基于砂岩油藏叠前地震数据的流体替换方法 |
CN102928870B (zh) * | 2012-09-21 | 2015-12-02 | 中国石油天然气股份有限公司勘探开发研究院廊坊分院 | 基于正则化的非线性地震叠前弹性参数反演方法 |
CN103076629B (zh) * | 2012-09-21 | 2014-05-14 | 中国石油天然气集团公司 | 基于纵横地震波速度比进行油气勘探的方法及设备 |
CN103176211B (zh) * | 2013-03-07 | 2016-04-20 | 中国石油集团东方地球物理勘探有限责任公司 | 基于多敏感弹性参数的含气储层预测方法及装置 |
CN103257361B (zh) * | 2013-05-24 | 2015-07-29 | 中国石油天然气集团公司 | 基于Zoeppritz方程近似式的油气预测方法及*** |
CN104280767A (zh) * | 2013-07-12 | 2015-01-14 | 中国石油天然气集团公司 | 一种基于柯西分布的稀疏脉冲反演方法 |
CN103485758A (zh) * | 2013-08-12 | 2014-01-01 | 中国石油天然气股份有限公司 | 一种利用多种资料拟合反演复杂岩性解释方法及装置 |
CN103412330B (zh) * | 2013-08-28 | 2016-05-11 | 中国科学院地质与地球物理研究所 | 一种地震属性数据处理方法及装置 |
CN104090302B (zh) * | 2014-07-02 | 2017-08-25 | 中国海洋石油总公司 | 工区地下介质频率域异常分析的方法 |
CN106932819B (zh) * | 2017-02-23 | 2019-01-01 | 河海大学 | 基于各向异性马尔科夫随机域的叠前地震参数反演方法 |
CN107193044B (zh) * | 2017-05-31 | 2019-01-22 | 河海大学 | 一种混合全局优化算法的叠前地震多参数反演方法 |
CN109521221B (zh) * | 2018-12-10 | 2020-09-29 | 东北大学 | 一种钻爆法施工硬岩隧道微震波波速实时获取方法 |
CN112817042B (zh) * | 2019-11-15 | 2024-01-30 | 中国石油天然气股份有限公司 | 沉积储层中油气识别方法及装置 |
CN113311482B (zh) * | 2020-06-09 | 2022-12-23 | 中国海洋石油集团有限公司 | 一种基于叠前谱反演优化的高分辨率中深层储层预测方法 |
CN112051609B (zh) * | 2020-08-21 | 2023-01-06 | 成都理工大学 | 地震波成像方法、***、存储介质、计算机程序、终端 |
CN112987088B (zh) * | 2021-02-22 | 2023-04-18 | 成都理工大学 | 一种渗流介质地震横波数值模拟和成像方法 |
CN113917538B (zh) * | 2021-09-18 | 2023-06-09 | 中国石油大学(华东) | 一种地震波形约束的叠前地震反演方法、设备及存储介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2616920A1 (fr) * | 1987-06-19 | 1988-12-23 | Schlumberger Prospection | Inversion d'un profil sismique vertical en minimisant une fonction du type entropie |
CN1164651A (zh) * | 1996-10-29 | 1997-11-12 | 中国石油天然气总公司西北地质研究所 | 地震-测井多信息储层参数反演方法 |
US6058073A (en) * | 1999-03-30 | 2000-05-02 | Atlantic Richfield Company | Elastic impedance estimation for inversion of far offset seismic sections |
US6219621B1 (en) * | 1998-06-30 | 2001-04-17 | Exxonmobil Upstream Research Co. | Sparse hyperbolic inversion of seismic data |
CN1504764A (zh) * | 2002-12-05 | 2004-06-16 | 大庆石油管理局 | 数据驱动岩性参数反演方法 |
CN1710446A (zh) * | 2005-06-21 | 2005-12-21 | 中国石油大学(北京) | 一种利用叠前地震波形反演构建虚拟井数据的方法 |
-
2007
- 2007-06-20 CN CN200710119093XA patent/CN101329405B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2616920A1 (fr) * | 1987-06-19 | 1988-12-23 | Schlumberger Prospection | Inversion d'un profil sismique vertical en minimisant une fonction du type entropie |
CN1164651A (zh) * | 1996-10-29 | 1997-11-12 | 中国石油天然气总公司西北地质研究所 | 地震-测井多信息储层参数反演方法 |
US6219621B1 (en) * | 1998-06-30 | 2001-04-17 | Exxonmobil Upstream Research Co. | Sparse hyperbolic inversion of seismic data |
US6058073A (en) * | 1999-03-30 | 2000-05-02 | Atlantic Richfield Company | Elastic impedance estimation for inversion of far offset seismic sections |
CN1504764A (zh) * | 2002-12-05 | 2004-06-16 | 大庆石油管理局 | 数据驱动岩性参数反演方法 |
CN1710446A (zh) * | 2005-06-21 | 2005-12-21 | 中国石油大学(北京) | 一种利用叠前地震波形反演构建虚拟井数据的方法 |
Non-Patent Citations (8)
Title |
---|
印兴耀.一种叠前地震记录的全波场正反演方法.石油物探.2004,43(3),217-223. |
尚作萍 |
张宏兵 |
张宏兵;尚作萍;杨长春;段秋梁.波阻抗反演正则参数估计.地球物理学报.2005,48(1),181-188. * |
张繁昌 |
张繁昌;印兴耀.一种叠前地震记录的全波场正反演方法.石油物探.2004,43(3),217-223. * |
杨长春 |
段秋梁.波阻抗反演正则参数估计.地球物理学报.2005,48(1),181-188. |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2615591C1 (ru) * | 2013-05-24 | 2017-04-05 | Эксонмобил Апстрим Рисерч Компани | Многопараметрическая инверсия через зависящую от сдвига упругую полноволновую инверсию (fwi) |
Also Published As
Publication number | Publication date |
---|---|
CN101329405A (zh) | 2008-12-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101329405B (zh) | 一种简单的多参数地震反演方法 | |
CN101551466B (zh) | 一种利用与偏移距有关的地震属性提高油气储层预测精度的方法 | |
CN101281253B (zh) | 一种利用振幅随偏移距变化特征提高油气检测精度的方法 | |
CN101598805B (zh) | 一种多分量地震资料层位对比和标定方法 | |
Hornby | Imaging of near-borehole structure using full-waveform sonic data | |
CN101446645B (zh) | 一种利用地震流体阻抗进行流体确定的方法 | |
US5596548A (en) | Seismic imaging using wave equation extrapolation | |
US8379482B1 (en) | Using seismic attributes for data alignment and seismic inversion in joint PP/PS seismic analysis | |
US6055482A (en) | Method of seismic signal processing | |
EP0060029B1 (en) | A method of determining the ratio of the velocities of compressional and shear waves in subterranean formations | |
CN103116184B (zh) | 一种表层调查中计算横波静校正量的方法 | |
CA2386568A1 (en) | System for estimating azimuthal variations in seismic data | |
CN101354444A (zh) | 一种确定地层岩性和孔隙流体的方法 | |
CN106556861A (zh) | 一种基于全方位地震资料的方位avo反演方法 | |
CN103675911A (zh) | 一种基于纵波和转换横波联合反演截距和梯度的方法 | |
Lynn et al. | Relationship of P-wave seismic attributes, azimuthal anisotropy, and commercial gas pay in 3-DP-wave multiazimuth data, Rulison Field, Piceance Basin, Colorado | |
CA2387760A1 (en) | Transfer function method of seismic signal processing and exploration | |
CN102305942B (zh) | 基于三参量的非线性avo流体判别方法 | |
US20120269035A1 (en) | Evaluating Prospects from P-Wave Seismic Data Using S-Wave Vertical Shear Profile Data | |
Zhang et al. | Seismic azimuthal impedance anisotropy in the Barnett Shale | |
Fortin et al. | Reflection imaging around a well with the EVA full-waveform tool | |
CN110244383A (zh) | 基于近地表数据的地质岩性综合模型创建方法 | |
CN109239777B (zh) | 一种利用联合反演方法检测构造煤发育的方法 | |
CN104965226A (zh) | 一种岩体中油气信息的分析方法 | |
Yang et al. | Pure S-waves in land P-wave source VSP 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 |