一种基于透射方程的各向异性参数反演方法
一种基于透射方程的各向异性参数反演方法
【技术领域】
[0001] 本发明属于地球物理反演领域，具体涉及一种基于透射方程的各向异性参数反演 方法。
【背景技术】
[0002] 各向异性参数在裂缝定量预测中起着举足轻重的作用，而由于各向异性参数数量 级远远小于地震数据数量级，为其计算带来了极大的困难。目前计算各向异性参数主要是 基于叠前方位道集，主要分为W下Η类渠道；1.根据随地震方位变化的速度信息，速度分 析精度难W满足方法需要；2.根据随地震方位变化的剩余时差信息，拾取剩余时差时存在 较多误差；3. W简化Ruger(1998)方程为基础（未经简化或者近似的Ruger方程的解是计 算不出来的，必须要通过简化方程或者一些近似关系将Ruger方程简化，才能得到Ruger方 程的解。并且现有方法中使用的Ruger方程都是Ruger反射系数方程，而不是Ruger透射 方程），基于各种假设条件的弹性参数反演，虽然研究表明简化后的Ruger方程对计算反射 系数影响不大，但是由此产生的传递误差对于计算数量级极小的各向异性参数影响还是十 分巨大的。
【发明内容】
[0003] 本发明的目的在于解决上述现有技术中存在的难题，提供一种基于透射方程的各 向异性参数反演方法，基于完整的Ruger透射系数方程，W叠后数据体方位之间的差异为 基础，得到各向异性参数A δ W，Λ ε W的解。
[0004] 本发明是通过W下技术方案实现的：
[0005] -种基于透射方程的各向异性参数反演方法，包括：
[0006] (1)选取任意Η个不同方位角带1，竿2，牺对应的叠后地震数据体，通过稀疏脉 冲的反演方法计算得到对应该Η个不同方位角的反射系数Rpp(化），Rpp(CP2)，Rpp((P3);
[0007] 似分别计算Η个不同方位角的反射系数Rpp如i)，： Rpp(cp2)，Κρρ(φ〇对应的透 射系数而如1)，Τρ(恥)，TV如3);
[000引做分析叠前方位道集获得入射角信息，所述入射角信息包括最大入射角i2、最小 入射角ii和间隔；
[0009] (4)利用方位角信息斬，範，恥，遽射系数信息Τρ如1)，Tp(q>2)，Τρ(φ:;)计算 得到各向异性参数Δ ε W和Λ δ W。
[0010] 所述步骤（2)是送样实现的：
[0011] 分别对S个不同方位角的反射系数Rpp(化），Κρρ(ψ2)，Κρρ(φ3)采用公式（11)得 到对应的透射系数Τρ如1)，Tp((pj，Τρ(φ3):
[001 引 Rpp+Tp=l (11)。
[0013] 所述步骤（4)是送样实现的：
[0014] 利用公式（6)和（7)计算得到各向异性参数Δ ε W和Λ δ W :
[0015]
[0018] 与现有技术相比，本发明的有益效果是：
[0019] (1)采用的是完整的Ruger透射系数方程，未对方程进行简化，理论上不存在误 差；
[0020] 似具有叠后地震反演稳定性高、信噪比高的优点；
[0021] (3)输入数据范围可控，实际应用中易于实现；
[0022] (4)计算效率高，可直接获得各向异性参数Δ ε W和Λ δ W。
【附图说明】
[0023] 图1本发明方法的步骤框图
[0024] 图2a模型数据中的化图 [00巧]图化模型数据中的Vs图
[0026] 图2c模型数据中的Den图
[0027] 图2d模型数据中的δ图
[0028] 图2e模型数据中的ε图
[0029] 图2f模型数据中的Υ图
[0030] 图3a正演得到的0°方位的叠前数据
[0031] 图3b演得到的30°方位的叠前数据
[0032] 图3c正演得到的60°方位的叠前数据
[0033] 图4a反演得到的0°方位对应的反射系数
[0034] 图4b反演得到的30°方位对应的反射系数
[0035] 图4c反演得到的60。方位对应的反射系数
[0036] 图5a计算得到的0°方位对应的透射系数
[0037] 图化计算得到的30°方位对应的透射系数 [003引 图5c计算得到的60°方位对应的透射系数
[0039] 图6a最终结果各向异性参数Δ ε W
[0040] 图化最终结果各向异性参数Δ δ W。
【具体实施方式】
[0041] 下面结合附图对本发明作进一步详细描述：
[0042] 目前绝大多数的各向异性参数反演是基于简化的Ruger反射系数方程，其近似或 者简化的过程引入了误差，而Ruger透射方程与反射系数方程具有相同的各向异性特征。 因此本发明基于完整的Ruger透射系数方程得到各向异性参数Δ δ W，Λ ε W的解，从理 论上不引入误差，计算结果精度更高，并且具有叠后反演的优点，反演结果可用于裂缝预测 和页岩气预测等热点领域。
[0043] 本发明W完整的Ruger透射方程为基础，计算出各向异性参数ε W，S W，方法适 用于ΗΤΙ (横向各向异性介质）。采用的Ruger透射方程与反射系数方程具有相同的各向异 性特征，并且未对方程进行简化，理论上不存在误差，计算结果精度更高，并且具有叠后反 演的优点，反演结果可用于裂缝预测和页岩气预测等热点领域
[0044] Ruger(2002)借助弱各向异性的概念，推导出HTI介质的纵波近似透射系数关系 式：
[00451
[0046] 式中：面，β巧P分别为HTI介质上下两层的纵波、横波速度平均值W及密度平均 值；Δ δ W，Λ ε W分别为上下两层各向异性参数差值；i和φ分别代表入射角和方位角。
[0047] 对于方程1按照入射角i方向进行叠加，方程变为：
[0048]
别代表入射角和方位角。
[0052]
[0054] 解由方程（3) (4)妨组成的方程组，得
[00 巧]
i和^分别代表入射角和方位角， 如1)，Γρ(仍），馬知％)分别代表爭=单1，斬，恥时的透射系数，上述为HTI介质中近似 各向异性参数计算公式，在此基础上进行各向异性参数Δ ε W和Λ δ W的反演研究。
[005引反射系数Rpp((Pi.)，RPP(CP2)，Rpp如3)通过稀疏脉冲反演计算，常规的稀疏脉冲 流程可分为Η步：
[00则 （1)反射系数反演
[0060] 采用最大似然反權积进行反射系数的反演，最大似然反權积对地层的假设认为： 地层的反射系数是由较大的反射界面的反射和具有高斯背景的小反射叠加组合而成，导出 一个最小目标函数：
[0061]
[0062] 式中，R2和Ν2分别为反射系数和噪音的均方值，r(K)和η (Κ)表示第Κ个采样点 的反射系数和噪音，Μ表示反射层数， L表示采样总数，λ表示给定反射系数的似然值。通 过多次迭代，求取反射系数。
[0063] (2)根据反射系数的反演结果结合阻抗趋势计算一个初始的波阻抗
[0064] 根据最大似然反權积计算得到的反射系数，结合初始阻抗模型，采用递推算法，反 演得到初始的波阻抗模型：
[0065] (9)
[006引式中，Z(i)为第i层的波阻抗值，R(i)为第i层的反射系数。
[0067] (3)结合井的约束条件进行波阻抗反演
[0068] 约束稀疏脉冲反演对每一道依据目标函数对计算出的初始波阻抗进行调整，包括 对反射系数的调整（即根据初始波阻抗（送时有一个对应初始波阻抗的反射系数）对实际 波阻抗进行逼近，使用过迭代寻优的方法，最终得到一个近似于实际波阻抗的结果，送时对 应的反射系数就是所需要的反射系数）。目标优化函数为：
[0069] F = Lp(r) + ALq(s-d) + a iLiAZ (10)
[0070] 式中，r为反射系数序列，Δ z为与阻抗趋势的差序列，d为地震道序列，s为合成 地震道序列，λ为残差权重因子，α为趋势权重因子，p、q为L模因子。具体的，右式第一 项反映了反射系数的绝对值和，第二项反映了合成声波记录与原始地震数据的差值，第Η 项为趋势约束项。
[0071] 本方法需要的就是（3)中调整后的反射系数Rpp(cpJ，Rpp(cp2)，Rpp((P:.;)。
[0072] 根据zeoppritz方程，在垂直入射情况下，反射系数与透射系数存在如下关系： [007引 Rpp+Tp=l (11)
[0074] 而本方法反演是W叠后数据为基础得到的反射系数Rpp，因此是近似地符合垂直入 射条件，可W通过方程（11)计算得到透射系数Τρ。
[0075] 本发明基于透射方程的各向异性参数反演方法的实现步骤为：
[007引 （1)选取任意Η个不同方位角化，Φ2，φ:3.对应的叠后地震数据体，通过稀疏脉 冲的反演方法计算得到对应该Η个不同方位角的反射系数Rpp(cPi)，Rpp(cp2)，Rpp(q):J:
[0077] 似通过方程（11)分别计算对应的透射系数Τρ(φ.ι)，Τρ(φ;;)，Τρ如3);
[0078] (3)分析叠前方位道集获得入射角信息，所述入射角信息包括最大入射角i2、 最小入射角ii和间隔，然后通过公苗
Η十算出B、C、D、E的值；
[007引（4)利用方位角信息化，化，爭3,透射系数信息Τρ如1)，Τρ如2)，了,，如3)和参数 B、C、D、E带入方程化）（7)中，计算得到各向异性参数Δ ε W和Λ δ W。
[0080] 图1示出本发明一种基于透射方程的各向异性参数反演方法的具体流程。下面参 照附图并结合实施例对本发明做进一步说明。图2a-图2f为使用的各向异性模型数据。图 3a-图3c为正演得到的0°、30°、60°方位的叠前数据。
[0081] 首先根据步骤（1)对0°，30°，60°方位的叠后数据体分别进行稀疏脉冲反演， 计算得到对应的反射系数体Rpp(〇° )，Rpp(30° )，Rpp化0° )，如图4a-图4c所示，然后根 据方程（11)分别计算对应的透射系数Τρ(φι)，Τρ恤2)，Τρ(φ,)，如图5a-图5c所示。
[0082] 然后根据步骤（3)分析叠前方位道集的入射角信息，本次叠前地震数据沿入射 角叠加的范围选取为ii = 〇°，i2 = 30°，入射角间隔为Γ，因此根据
.分别计算出B = 1.3607， C = 0. 2636,D = 2. 7214,E = 2. 985。对于实际资料应用过程，由于目前国内绝大多数观测 系统为树状观测系统，导致从根本上讲叠前道集会存在近、远偏移距（近远入射角）能量分 布较弱，目前尚未有有效的方法可w解决叠前道集能量问题，比较常见的方法是对叠前道 集做能量均衡，但是由于缺乏理论的依据，往往容易使原本存在的AVO特征扭曲，为之后的 叠前反演带来巨大的影响。本发明在叠前道集叠加至叠后数据处对于入射角范围ii、i2可 根据叠前道集质量进行分选，虽然在一定程度上造成了能量的损失，但是可W获得更稳定 的叠后数据，并且在理论上获得验证，反演效果不受影响。
[008引接下来根据步骤（4)将已知的方位角信息化=0。，q>2 = 3估。，q>3 = 6估。， 透射系数Τρ(φ.ι)，Τρ(φ,), Τ;，(φ:,)和b、C、D、E的值，代入公式化）、（7)中，计算得到各向 异性参数Δ ε W和Λ δ W，如图6曰、图化所示，反演结果可W用与裂缝定量预测。
[0084] 上述技术方案只是本发明的一种实施方式，对于本领域内的技术人员而言，在本 发明公开了应用方法和原理的基础上，很容易做出各种类型的改进或变形，而不仅限于本 发明上述【具体实施方式】所描述的方法，因此前面描述的方式只是优选的，而并不具有限制 性的意义。
【主权项】
1. 一种基于透射方程的各向异性参数反演方法，其特征在于：所述方法包括： (1) 选取任意三个不同方位角φι，Φυ fPs对应的叠后地震数据体，通过稀疏脉冲的 反演方法计算得到对应该三个不同方位角的反射系数(2) 分别计算三个不同方位角的反射系数对应的透射系 数(3) 分析叠前方位道集获得入射角信息，所述入射角信息包括最大入射角i2、最小入射 角I1和间隔； ⑷利用方位角信息物:，φ2，φ3:，透射系数信肩计算得到 各向异性参数Λ ε (ν)和Λ S w。2. 根据权利要求1所述的基于透射方程的各向异性参数反演方法，其特征在于：所述 步骤（2)是这样实现的： 分别对三个不同方位角的反射系数I采用公式（11)得到对 应的透射系数RPP+TP = 1 (11)〇3. 根据权利要求2所述的基于透射方程的各向异性参数反演方法，其特征在于：所述 步骤（4)是这样实现的： 利用公式（6)和（7)计算得到各向异性参数Λ ε (v)和Λ δ (v):
【专利摘要】本发明提供了一种基于透射方程的各向异性参数反演方法，属于地球物理反演领域。所述方法包括：(1)选取任意三个不同方位角对应的叠后地震数据体，通过稀疏脉冲的反演方法计算得到对应该三个不同方位角的反射系数(2)分别计算三个不同方位角的反射系数对应的透射系数(3)分析叠前方位道集获得入射角信息，所述入射角信息包括最大入射角i2、最小入射角i1和间隔；(4)利用方位角信息透射系数信息计算得到各向异性参数Δε(V)和Δδ(V)。
【IPC分类】G01V1/28
【公开号】CN105487110
【申请号】CN201410471907
【发明人】张克非, 肖鹏飞, 李呈呈, 周单, 须振华, 许凯
【申请人】中国石油化工股份有限公司, 中国石油化工股份有限公司石油物探技术研究院
【公开日】2016年4月13日
【申请日】2014年9月16日
【技术领域】
[0001] 本发明属于地球物理反演领域，具体涉及一种基于透射方程的各向异性参数反演 方法。
【背景技术】
[0002] 各向异性参数在裂缝定量预测中起着举足轻重的作用，而由于各向异性参数数量 级远远小于地震数据数量级，为其计算带来了极大的困难。目前计算各向异性参数主要是 基于叠前方位道集，主要分为W下Η类渠道；1.根据随地震方位变化的速度信息，速度分 析精度难W满足方法需要；2.根据随地震方位变化的剩余时差信息，拾取剩余时差时存在 较多误差；3. W简化Ruger(1998)方程为基础（未经简化或者近似的Ruger方程的解是计 算不出来的，必须要通过简化方程或者一些近似关系将Ruger方程简化，才能得到Ruger方 程的解。并且现有方法中使用的Ruger方程都是Ruger反射系数方程，而不是Ruger透射 方程），基于各种假设条件的弹性参数反演，虽然研究表明简化后的Ruger方程对计算反射 系数影响不大，但是由此产生的传递误差对于计算数量级极小的各向异性参数影响还是十 分巨大的。
【发明内容】
[0003] 本发明的目的在于解决上述现有技术中存在的难题，提供一种基于透射方程的各 向异性参数反演方法，基于完整的Ruger透射系数方程，W叠后数据体方位之间的差异为 基础，得到各向异性参数A δ W，Λ ε W的解。
[0004] 本发明是通过W下技术方案实现的：
[0005] -种基于透射方程的各向异性参数反演方法，包括：
[0006] (1)选取任意Η个不同方位角带1，竿2，牺对应的叠后地震数据体，通过稀疏脉 冲的反演方法计算得到对应该Η个不同方位角的反射系数Rpp(化），Rpp(CP2)，Rpp((P3);
[0007] 似分别计算Η个不同方位角的反射系数Rpp如i)，： Rpp(cp2)，Κρρ(φ〇对应的透 射系数而如1)，Τρ(恥)，TV如3);
[000引做分析叠前方位道集获得入射角信息，所述入射角信息包括最大入射角i2、最小 入射角ii和间隔；
[0009] (4)利用方位角信息斬，範，恥，遽射系数信息Τρ如1)，Tp(q>2)，Τρ(φ:;)计算 得到各向异性参数Δ ε W和Λ δ W。
[0010] 所述步骤（2)是送样实现的：
[0011] 分别对S个不同方位角的反射系数Rpp(化），Κρρ(ψ2)，Κρρ(φ3)采用公式（11)得 到对应的透射系数Τρ如1)，Tp((pj，Τρ(φ3):
[001 引 Rpp+Tp=l (11)。
[0013] 所述步骤（4)是送样实现的：
[0014] 利用公式（6)和（7)计算得到各向异性参数Δ ε W和Λ δ W :
[0015]
[0018] 与现有技术相比，本发明的有益效果是：
[0019] (1)采用的是完整的Ruger透射系数方程，未对方程进行简化，理论上不存在误 差；
[0020] 似具有叠后地震反演稳定性高、信噪比高的优点；
[0021] (3)输入数据范围可控，实际应用中易于实现；
[0022] (4)计算效率高，可直接获得各向异性参数Δ ε W和Λ δ W。
【附图说明】
[0023] 图1本发明方法的步骤框图
[0024] 图2a模型数据中的化图 [00巧]图化模型数据中的Vs图
[0026] 图2c模型数据中的Den图
[0027] 图2d模型数据中的δ图
[0028] 图2e模型数据中的ε图
[0029] 图2f模型数据中的Υ图
[0030] 图3a正演得到的0°方位的叠前数据
[0031] 图3b演得到的30°方位的叠前数据
[0032] 图3c正演得到的60°方位的叠前数据
[0033] 图4a反演得到的0°方位对应的反射系数
[0034] 图4b反演得到的30°方位对应的反射系数
[0035] 图4c反演得到的60。方位对应的反射系数
[0036] 图5a计算得到的0°方位对应的透射系数
[0037] 图化计算得到的30°方位对应的透射系数 [003引 图5c计算得到的60°方位对应的透射系数
[0039] 图6a最终结果各向异性参数Δ ε W
[0040] 图化最终结果各向异性参数Δ δ W。
【具体实施方式】
[0041] 下面结合附图对本发明作进一步详细描述：
[0042] 目前绝大多数的各向异性参数反演是基于简化的Ruger反射系数方程，其近似或 者简化的过程引入了误差，而Ruger透射方程与反射系数方程具有相同的各向异性特征。 因此本发明基于完整的Ruger透射系数方程得到各向异性参数Δ δ W，Λ ε W的解，从理 论上不引入误差，计算结果精度更高，并且具有叠后反演的优点，反演结果可用于裂缝预测 和页岩气预测等热点领域。
[0043] 本发明W完整的Ruger透射方程为基础，计算出各向异性参数ε W，S W，方法适 用于ΗΤΙ (横向各向异性介质）。采用的Ruger透射方程与反射系数方程具有相同的各向异 性特征，并且未对方程进行简化，理论上不存在误差，计算结果精度更高，并且具有叠后反 演的优点，反演结果可用于裂缝预测和页岩气预测等热点领域
[0044] Ruger(2002)借助弱各向异性的概念，推导出HTI介质的纵波近似透射系数关系 式：
[00451
[0046] 式中：面，β巧P分别为HTI介质上下两层的纵波、横波速度平均值W及密度平均 值；Δ δ W，Λ ε W分别为上下两层各向异性参数差值；i和φ分别代表入射角和方位角。
[0047] 对于方程1按照入射角i方向进行叠加，方程变为：
[0048]
别代表入射角和方位角。
[0052]
[0054] 解由方程（3) (4)妨组成的方程组，得
[00 巧]
i和^分别代表入射角和方位角， 如1)，Γρ(仍），馬知％)分别代表爭=单1，斬，恥时的透射系数，上述为HTI介质中近似 各向异性参数计算公式，在此基础上进行各向异性参数Δ ε W和Λ δ W的反演研究。
[005引反射系数Rpp((Pi.)，RPP(CP2)，Rpp如3)通过稀疏脉冲反演计算，常规的稀疏脉冲 流程可分为Η步：
[00则 （1)反射系数反演
[0060] 采用最大似然反權积进行反射系数的反演，最大似然反權积对地层的假设认为： 地层的反射系数是由较大的反射界面的反射和具有高斯背景的小反射叠加组合而成，导出 一个最小目标函数：
[0061]
[0062] 式中，R2和Ν2分别为反射系数和噪音的均方值，r(K)和η (Κ)表示第Κ个采样点 的反射系数和噪音，Μ表示反射层数， L表示采样总数，λ表示给定反射系数的似然值。通 过多次迭代，求取反射系数。
[0063] (2)根据反射系数的反演结果结合阻抗趋势计算一个初始的波阻抗
[0064] 根据最大似然反權积计算得到的反射系数，结合初始阻抗模型，采用递推算法，反 演得到初始的波阻抗模型：
[0065] (9)
[006引式中，Z(i)为第i层的波阻抗值，R(i)为第i层的反射系数。
[0067] (3)结合井的约束条件进行波阻抗反演
[0068] 约束稀疏脉冲反演对每一道依据目标函数对计算出的初始波阻抗进行调整，包括 对反射系数的调整（即根据初始波阻抗（送时有一个对应初始波阻抗的反射系数）对实际 波阻抗进行逼近，使用过迭代寻优的方法，最终得到一个近似于实际波阻抗的结果，送时对 应的反射系数就是所需要的反射系数）。目标优化函数为：
[0069] F = Lp(r) + ALq(s-d) + a iLiAZ (10)
[0070] 式中，r为反射系数序列，Δ z为与阻抗趋势的差序列，d为地震道序列，s为合成 地震道序列，λ为残差权重因子，α为趋势权重因子，p、q为L模因子。具体的，右式第一 项反映了反射系数的绝对值和，第二项反映了合成声波记录与原始地震数据的差值，第Η 项为趋势约束项。
[0071] 本方法需要的就是（3)中调整后的反射系数Rpp(cpJ，Rpp(cp2)，Rpp((P:.;)。
[0072] 根据zeoppritz方程，在垂直入射情况下，反射系数与透射系数存在如下关系： [007引 Rpp+Tp=l (11)
[0074] 而本方法反演是W叠后数据为基础得到的反射系数Rpp，因此是近似地符合垂直入 射条件，可W通过方程（11)计算得到透射系数Τρ。
[0075] 本发明基于透射方程的各向异性参数反演方法的实现步骤为：
[007引 （1)选取任意Η个不同方位角化，Φ2，φ:3.对应的叠后地震数据体，通过稀疏脉 冲的反演方法计算得到对应该Η个不同方位角的反射系数Rpp(cPi)，Rpp(cp2)，Rpp(q):J:
[0077] 似通过方程（11)分别计算对应的透射系数Τρ(φ.ι)，Τρ(φ;;)，Τρ如3);
[0078] (3)分析叠前方位道集获得入射角信息，所述入射角信息包括最大入射角i2、 最小入射角ii和间隔，然后通过公苗
Η十算出B、C、D、E的值；
[007引（4)利用方位角信息化，化，爭3,透射系数信息Τρ如1)，Τρ如2)，了,，如3)和参数 B、C、D、E带入方程化）（7)中，计算得到各向异性参数Δ ε W和Λ δ W。
[0080] 图1示出本发明一种基于透射方程的各向异性参数反演方法的具体流程。下面参 照附图并结合实施例对本发明做进一步说明。图2a-图2f为使用的各向异性模型数据。图 3a-图3c为正演得到的0°、30°、60°方位的叠前数据。
[0081] 首先根据步骤（1)对0°，30°，60°方位的叠后数据体分别进行稀疏脉冲反演， 计算得到对应的反射系数体Rpp(〇° )，Rpp(30° )，Rpp化0° )，如图4a-图4c所示，然后根 据方程（11)分别计算对应的透射系数Τρ(φι)，Τρ恤2)，Τρ(φ,)，如图5a-图5c所示。
[0082] 然后根据步骤（3)分析叠前方位道集的入射角信息，本次叠前地震数据沿入射 角叠加的范围选取为ii = 〇°，i2 = 30°，入射角间隔为Γ，因此根据
.分别计算出B = 1.3607， C = 0. 2636,D = 2. 7214,E = 2. 985。对于实际资料应用过程，由于目前国内绝大多数观测 系统为树状观测系统，导致从根本上讲叠前道集会存在近、远偏移距（近远入射角）能量分 布较弱，目前尚未有有效的方法可w解决叠前道集能量问题，比较常见的方法是对叠前道 集做能量均衡，但是由于缺乏理论的依据，往往容易使原本存在的AVO特征扭曲，为之后的 叠前反演带来巨大的影响。本发明在叠前道集叠加至叠后数据处对于入射角范围ii、i2可 根据叠前道集质量进行分选，虽然在一定程度上造成了能量的损失，但是可W获得更稳定 的叠后数据，并且在理论上获得验证，反演效果不受影响。
[008引接下来根据步骤（4)将已知的方位角信息化=0。，q>2 = 3估。，q>3 = 6估。， 透射系数Τρ(φ.ι)，Τρ(φ,), Τ;，(φ:,)和b、C、D、E的值，代入公式化）、（7)中，计算得到各向 异性参数Δ ε W和Λ δ W，如图6曰、图化所示，反演结果可W用与裂缝定量预测。
[0084] 上述技术方案只是本发明的一种实施方式，对于本领域内的技术人员而言，在本 发明公开了应用方法和原理的基础上，很容易做出各种类型的改进或变形，而不仅限于本 发明上述【具体实施方式】所描述的方法，因此前面描述的方式只是优选的，而并不具有限制 性的意义。
【主权项】
1. 一种基于透射方程的各向异性参数反演方法，其特征在于：所述方法包括： (1) 选取任意三个不同方位角φι，Φυ fPs对应的叠后地震数据体，通过稀疏脉冲的 反演方法计算得到对应该三个不同方位角的反射系数(2) 分别计算三个不同方位角的反射系数对应的透射系 数(3) 分析叠前方位道集获得入射角信息，所述入射角信息包括最大入射角i2、最小入射 角I1和间隔； ⑷利用方位角信息物:，φ2，φ3:，透射系数信肩计算得到 各向异性参数Λ ε (ν)和Λ S w。2. 根据权利要求1所述的基于透射方程的各向异性参数反演方法，其特征在于：所述 步骤（2)是这样实现的： 分别对三个不同方位角的反射系数I采用公式（11)得到对 应的透射系数RPP+TP = 1 (11)〇3. 根据权利要求2所述的基于透射方程的各向异性参数反演方法，其特征在于：所述 步骤（4)是这样实现的： 利用公式（6)和（7)计算得到各向异性参数Λ ε (v)和Λ δ (v):
【专利摘要】本发明提供了一种基于透射方程的各向异性参数反演方法，属于地球物理反演领域。所述方法包括：(1)选取任意三个不同方位角对应的叠后地震数据体，通过稀疏脉冲的反演方法计算得到对应该三个不同方位角的反射系数(2)分别计算三个不同方位角的反射系数对应的透射系数(3)分析叠前方位道集获得入射角信息，所述入射角信息包括最大入射角i2、最小入射角i1和间隔；(4)利用方位角信息透射系数信息计算得到各向异性参数Δε(V)和Δδ(V)。
【IPC分类】G01V1/28
【公开号】CN105487110
【申请号】CN201410471907
【发明人】张克非, 肖鹏飞, 李呈呈, 周单, 须振华, 许凯
【申请人】中国石油化工股份有限公司, 中国石油化工股份有限公司石油物探技术研究院
【公开日】2016年4月13日
【申请日】2014年9月16日
最新回复(0)