CN105487110B - 一种基于透射方程的各向异性参数反演方法 - Google Patents
一种基于透射方程的各向异性参数反演方法 Download PDFInfo
- Publication number
- CN105487110B CN105487110B CN201410471907.6A CN201410471907A CN105487110B CN 105487110 B CN105487110 B CN 105487110B CN 201410471907 A CN201410471907 A CN 201410471907A CN 105487110 B CN105487110 B CN 105487110B
- Authority
- CN
- China
- Prior art keywords
- anisotropic parameters
- calculated
- reflectance factor
- different orientations
- equation
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种基于透射方程的各向异性参数反演方法,属于地球物理反演领域。所述方法包括:(1)选取任意三个不同方位角对应的叠后地震数据体,通过稀疏脉冲的反演方法计算得到对应该三个不同方位角的反射系数(2)分别计算三个不同方位角的反射系数 对应的透射系数(3)分析叠前方位道集获得入射角信息,所述入射角信息包括最大入射角i2、最小入射角i1和间隔;(4)利用方位角信息透射系数信息计算得到各向异性参数Δε(V)和Δδ(V)。
Description
技术领域
本发明属于地球物理反演领域,具体涉及一种基于透射方程的各向异性参数反演方法。
背景技术
各向异性参数在裂缝定量预测中起着举足轻重的作用,而由于各向异性参数数量级远远小于地震数据数量级,为其计算带来了极大的困难。目前计算各向异性参数主要是基于叠前方位道集,主要分为以下三类渠道:1.根据随地震方位变化的速度信息,速度分析精度难以满足方法需要;2.根据随地震方位变化的剩余时差信息,拾取剩余时差时存在较多误差;3.以简化Ruger(1998)方程为基础(未经简化或者近似的Ruger方程的解是计算不出来的,必须要通过简化方程或者一些近似关系将Ruger方程简化,才能得到Ruger方程的解。并且现有方法中使用的Ruger方程都是Ruger反射系数方程,而不是Ruger透射方程),基于各种假设条件的弹性参数反演,虽然研究表明简化后的Ruger方程对计算反射系数影响不大,但是由此产生的传递误差对于计算数量级极小的各向异性参数影响还是十分巨大的。
发明内容
本发明的目的在于解决上述现有技术中存在的难题,提供一种基于透射方程的各向异性参数反演方法,基于完整的Ruger透射系数方程,以叠后数据体方位之间的差异为基础,得到各向异性参数Δδ(V),Δε(V)的解。
本发明是通过以下技术方案实现的:
一种基于透射方程的各向异性参数反演方法,包括:
(1)选取任意三个不同方位角对应的叠后地震数据体,通过稀疏脉冲的反演方法计算得到对应该三个不同方位角的反射系数
(2)分别计算三个不同方位角的反射系数对应的透射系数
(3)分析叠前方位道集获得入射角信息,所述入射角信息包括最大入射角i2、最小入射角i1和间隔;
(4)利用方位角信息透射系数信息计算得到各向异性参数Δε(V)和Δδ(V)。
所述步骤(2)是这样实现的:
分别对三个不同方位角的反射系数采用公式(11)得到对应的透射系数
RPP+Tp=1 (11)。
所述步骤(4)是这样实现的:
利用公式(6)和(7)计算得到各向异性参数Δε(V)和Δδ(V):
其中,
与现有技术相比,本发明的有益效果是:
(1)采用的是完整的Ruger透射系数方程,未对方程进行简化,理论上不存在误差;
(2)具有叠后地震反演稳定性高、信噪比高的优点;
(3)输入数据范围可控,实际应用中易于实现;
(4)计算效率高,可直接获得各向异性参数Δε(V)和Δδ(V)。
附图说明
图1本发明方法的步骤框图
图2a模型数据中的Vp图
图2b模型数据中的Vs图
图2c模型数据中的Den图
图2d模型数据中的δ图
图2e模型数据中的ε图
图2f模型数据中的γ图
图3a正演得到的0°方位的叠前数据
图3b演得到的30°方位的叠前数据
图3c正演得到的60°方位的叠前数据
图4a反演得到的0°方位对应的反射系数
图4b反演得到的30°方位对应的反射系数
图4c反演得到的60°方位对应的反射系数
图5a计算得到的0°方位对应的透射系数
图5b计算得到的30°方位对应的透射系数
图5c计算得到的60°方位对应的透射系数
图6a最终结果各向异性参数Δε(V)
图6b最终结果各向异性参数Δδ(V)。
具体实施方式
下面结合附图对本发明作进一步详细描述:
目前绝大多数的各向异性参数反演是基于简化的Ruger反射系数方程,其近似或者简化的过程引入了误差,而Ruger透射方程与反射系数方程具有相同的各向异性特征。因此本发明基于完整的Ruger透射系数方程得到各向异性参数Δδ(V),Δε(V)的解,从理论上不引入误差,计算结果精度更高,并且具有叠后反演的优点,反演结果可用于裂缝预测和页岩气预测等热点领域。
本发明以完整的Ruger透射方程为基础,计算出各向异性参数ε(V),δ(V),方法适用于HTI(横向各向异性介质)。采用的Ruger透射方程与反射系数方程具有相同的各向异性特征,并且未对方程进行简化,理论上不存在误差,计算结果精度更高,并且具有叠后反演的优点,反演结果可用于裂缝预测和页岩气预测等热点领域
Ruger(2002)借助弱各向异性的概念,推导出HTI介质的纵波近似透射系数关系式:
式中:和分别为HTI介质上下两层的纵波、横波速度平均值以及密度平均值;Δδ(V),Δε(V)分别为上下两层各向异性参数差值;i和分别代表入射角和方位角。
对于方程1按照入射角i方向进行叠加,方程变为:
其中: i和分别代表入射角和方位角。
分别令
解由方程(3)(4)(5)组成的方程组,得
其中 i和分别代表入射角和方位角, 分别代表时的透射系数,上述为HTI介质中近似各向异性参数计算公式,在此基础上进行各向异性参数Δε(V)和Δδ(V)的反演研究。
反射系数通过稀疏脉冲反演计算,常规的稀疏脉冲流程可分为三步:
(1)反射系数反演
采用最大似然反褶积进行反射系数的反演,最大似然反褶积对地层的假设认为:地层的反射系数是由较大的反射界面的反射和具有高斯背景的小反射叠加组合而成,导出一个最小目标函数:
式中,R2和N2分别为反射系数和噪音的均方值,r(K)和n(K)表示第K个采样点的反射系数和噪音,M表示反射层数,L表示采样总数,λ表示给定反射系数的似然值。通过多次迭代,求取反射系数。
(2)根据反射系数的反演结果结合阻抗趋势计算一个初始的波阻抗
根据最大似然反褶积计算得到的反射系数,结合初始阻抗模型,采用递推算法,反演得到初始的波阻抗模型:
式中,Z(i)为第i层的波阻抗值,R(i)为第i层的反射系数。
(3)结合井的约束条件进行波阻抗反演
约束稀疏脉冲反演对每一道依据目标函数对计算出的初始波阻抗进行调整,包括对反射系数的调整(即根据初始波阻抗(这时有一个对应初始波阻抗的反射系数)对实际波阻抗进行逼近,使用过迭代寻优的方法,最终得到一个近似于实际波阻抗的结果,这时对应的反射系数就是所需要的反射系数)。目标优化函数为:
F=Lp(r)+λLq(s-d)+α-1L1ΔZ (10)
式中,r为反射系数序列,Δz为与阻抗趋势的差序列,d为地震道序列,s为合成地震道序列,λ为残差权重因子,α为趋势权重因子,p、q为L模因子。具体的,右式第一项反映了反射系数的绝对值和,第二项反映了合成声波记录与原始地震数据的差值,第三项为趋势约束项。
本方法需要的就是(3)中调整后的反射系数
根据zeoppritz方程,在垂直入射情况下,反射系数与透射系数存在如下关系:
RPP+Tp=1 (11)
而本方法反演是以叠后数据为基础得到的反射系数RPP,因此是近似地符合垂直入射条件,可以通过方程(11)计算得到透射系数Tp。
本发明基于透射方程的各向异性参数反演方法的实现步骤为:
(1)选取任意三个不同方位角对应的叠后地震数据体,通过稀疏脉冲的反演方法计算得到对应该三个不同方位角的反射系数
(2)通过方程(11)分别计算对应的透射系数
(3)分析叠前方位道集获得入射角信息,所述入射角信息包括最大入射角i2、最小入射角i1和间隔,然后通过公式 计算出B、C、D、E的值;
(4)利用方位角信息透射系数信息和参数B、C、D、E带入方程(6)(7)中,计算得到各向异性参数Δε(V)和Δδ(V)。
图1示出本发明一种基于透射方程的各向异性参数反演方法的具体流程。下面参照附图并结合实施例对本发明做进一步说明。图2a-图2f为使用的各向异性模型数据。图3a-图3c为正演得到的0°、30°、60°方位的叠前数据。
首先根据步骤(1)对0°,30°,60°方位的叠后数据体分别进行稀疏脉冲反演,计算得到对应的反射系数体RPP(0°),RPP(30°),RPP(60°),如图4a-图4c所示,然后根据方程(11)分别计算对应的透射系数如图5a-图5c所示。
然后根据步骤(3)分析叠前方位道集的入射角信息,本次叠前地震数据沿入射角叠加的范围选取为i1=0°,i2=30°,入射角间隔为1°,因此根据 分别计算出B=1.3607,C=0.2636,D=2.7214,E=2.985。对于实际资料应用过程,由于目前国内绝大多数观测系统为树状观测系统,导致从根本上讲叠前道集会存在近、远偏移距(近远入射角)能量分布较弱,目前尚未有有效的方法可以解决叠前道集能量问题,比较常见的方法是对叠前道集做能量均衡,但是由于缺乏理论的依据,往往容易使原本存在的AVO特征扭曲,为之后的叠前反演带来巨大的影响。本发明在叠前道集叠加至叠后数据处对于入射角范围i1、i2可根据叠前道集质量进行分选,虽然在一定程度上造成了能量的损失,但是可以获得更稳定的叠后数据,并且在理论上获得验证,反演效果不受影响。
接下来根据步骤(4)将已知的方位角信息透射系数和B、C、D、E的值,代入公式(6)、(7)中,计算得到各向异性参数Δε(V)和Δδ(V),如图6a、图6b所示,反演结果可以用与裂缝定量预测。
上述技术方案只是本发明的一种实施方式,对于本领域内的技术人员而言,在本发明公开了应用方法和原理的基础上,很容易做出各种类型的改进或变形,而不仅限于本发明上述具体实施方式所描述的方法,因此前面描述的方式只是优选的,而并不具有限制性的意义。
Claims (2)
1.一种基于透射方程的各向异性参数反演方法,其特征在于:所述方法包括:
(1)选取任意三个不同方位角对应的叠后地震数据体,通过稀疏脉冲的反演方法计算得到对应该三个不同方位角的反射系数
(2)分别计算三个不同方位角的反射系数对应的透射系数
(3)分析叠前方位道集获得入射角信息,所述入射角信息包括最大入射角i2、最小入射角i1和间隔;
(4)利用方位角信息透射系数信息计算得到各向异性参数Δε(V)和Δδ(V),其中利用公式(6)和(7)计算得到各向异性参数Δε(V)和Δδ(v):
其中,
2.根据权利要求1所述的基于透射方程的各向异性参数反演方法,其特征在于:所述步骤(2)是这样实现的:
分别对三个不同方位角的反射系数采用公式(11)得到对应的透射系数
Rpp+Tp=1 (11)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410471907.6A CN105487110B (zh) | 2014-09-16 | 2014-09-16 | 一种基于透射方程的各向异性参数反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410471907.6A CN105487110B (zh) | 2014-09-16 | 2014-09-16 | 一种基于透射方程的各向异性参数反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105487110A CN105487110A (zh) | 2016-04-13 |
CN105487110B true CN105487110B (zh) | 2018-05-04 |
Family
ID=55674216
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410471907.6A Active CN105487110B (zh) | 2014-09-16 | 2014-09-16 | 一种基于透射方程的各向异性参数反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105487110B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106353807B (zh) * | 2016-08-08 | 2018-08-14 | 中国石油天然气集团公司 | 裂缝识别方法和装置 |
CN110516358A (zh) * | 2019-08-28 | 2019-11-29 | 电子科技大学 | 一种基于稀疏表示的地震各向异性参数反演方法 |
CN113009571B (zh) * | 2021-02-18 | 2022-03-08 | 中国矿业大学(北京) | 双相介质中的水平裂缝的反射系数和透射系数确定方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102053267A (zh) * | 2010-10-22 | 2011-05-11 | 中国石油化工股份有限公司 | 一种地震剖面资料处理中基于参数反演的vsp波场分离方法 |
CN102455436A (zh) * | 2010-11-02 | 2012-05-16 | 中国石油大学(北京) | 有限方位角条件下压噪叠前纵波各向异性裂缝检测方法 |
CN102540250A (zh) * | 2010-12-08 | 2012-07-04 | 同济大学 | 基于方位保真角度域成像的裂缝型油气储层地震探测方法 |
CN102540251A (zh) * | 2010-12-16 | 2012-07-04 | 中国石油天然气集团公司 | 二维的横向各向异性介质hti叠前深度偏移建模方法及装置 |
CN102749645A (zh) * | 2012-03-14 | 2012-10-24 | 中国石油天然气股份有限公司 | 一种利用角度阻抗梯度进行储层烃类检测的方法及装置 |
CN102854527A (zh) * | 2012-07-13 | 2013-01-02 | 孙赞东 | 基于纵波方位avo的裂缝流体识别方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6904368B2 (en) * | 2002-11-12 | 2005-06-07 | Landmark Graphics Corporation | Seismic analysis using post-imaging seismic anisotropy corrections |
-
2014
- 2014-09-16 CN CN201410471907.6A patent/CN105487110B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102053267A (zh) * | 2010-10-22 | 2011-05-11 | 中国石油化工股份有限公司 | 一种地震剖面资料处理中基于参数反演的vsp波场分离方法 |
CN102455436A (zh) * | 2010-11-02 | 2012-05-16 | 中国石油大学(北京) | 有限方位角条件下压噪叠前纵波各向异性裂缝检测方法 |
CN102540250A (zh) * | 2010-12-08 | 2012-07-04 | 同济大学 | 基于方位保真角度域成像的裂缝型油气储层地震探测方法 |
CN102540251A (zh) * | 2010-12-16 | 2012-07-04 | 中国石油天然气集团公司 | 二维的横向各向异性介质hti叠前深度偏移建模方法及装置 |
CN102749645A (zh) * | 2012-03-14 | 2012-10-24 | 中国石油天然气股份有限公司 | 一种利用角度阻抗梯度进行储层烃类检测的方法及装置 |
CN102854527A (zh) * | 2012-07-13 | 2013-01-02 | 孙赞东 | 基于纵波方位avo的裂缝流体识别方法 |
Non-Patent Citations (4)
Title |
---|
"Plane-wave reflection and transmission coefficients for a transversely isotropic solid";Mark Graebner;《Geophysics》;19921130;第57卷(第11期);第1512-1519页 * |
"Weak-contrast reflection-transmission coefficients in a generally anisotropic background";Ludek Klimes;《Geophysics》;20031231;第68卷(第6期);第2063-2073页 * |
"基于各向异性AVO的裂缝弹性参数叠前反演方法";张广智 等;《吉林大学学报(地球科学版)》;20120531;第42卷(第3期);第845-851、871页 * |
"裂缝介质的纵波方位AVO反演研究";朱兆林 等;《石油勘探》;20050930;第44卷(第5期);第499-503页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105487110A (zh) | 2016-04-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103293552B (zh) | 一种叠前地震资料的反演方法及系统 | |
CN103105624B (zh) | 基于数据库技术的纵横波时差定位方法 | |
CN102759746B (zh) | 一种变偏移距垂直地震剖面数据反演各向异性参数方法 | |
CN103399346B (zh) | 一种井震联合初始波阻抗建模方法 | |
CN105319590B (zh) | 一种基于hti介质的各向异性单参数反演方法 | |
CN104570079A (zh) | 一种纵波、转换横波地震资料的时间匹配方法 | |
CN106443775A (zh) | 高分辨率转换波裂缝预测方法 | |
CN104502997A (zh) | 一种利用裂缝密度曲线预测裂缝密度体的方法 | |
CN105089652A (zh) | 一种拟声波曲线重构与稀疏脉冲联合反演方法 | |
CN110780351B (zh) | 纵波和转换波叠前联合反演方法及系统 | |
CN106646601A (zh) | 多信息联合约束的浅中深层三维q体建立方法 | |
CN102866426B (zh) | 一种利用avo大角度道集分析岩体油气信息的方法 | |
CN103412324B (zh) | 一种估计介质品质因子的epifvo方法 | |
CN104570076A (zh) | 一种基于二分法的地震波初至自动拾取方法 | |
CN105487110B (zh) | 一种基于透射方程的各向异性参数反演方法 | |
CN104155687A (zh) | 一种相控叠后声波阻抗反演方法 | |
CN104237937A (zh) | 叠前地震反演方法及其系统 | |
CN105093278A (zh) | 基于激发主能量优化算法的全波形反演梯度算子的提取方法 | |
CN108398720A (zh) | 一种基于杨氏模量、泊松比的两项式地震叠前反演方法 | |
CN104122581A (zh) | 一种叠后声波阻抗反演方法 | |
Du et al. | Pre-stack seismic inversion using SeisInv-ResNet | |
CN104965222A (zh) | 三维纵波阻抗全波形反演方法及装置 | |
CN106501872B (zh) | 一种裂缝储层地应力特征的计算方法及装置 | |
CN104297800A (zh) | 一种自相控叠前反演方法 | |
CN103984018A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |