CN104820239A - 一种方位叠前地震属性解耦提取方法 - Google Patents

Abstract

本发明提供一种方位叠前地震属性解耦提取方法，基于水平裂缝非均质地层假设，从页岩工区多方位叠前角度道集中解耦提取各向同性梯度属性及不同对称面内的各向异性梯度属性，具体步骤：首先充分分析页岩工区水平裂缝引起的方位各向异性特征，推导由水平裂缝引起的方位各向异性两项AVO近似方程，提取AVO截距属性及方位AVO梯度属性；然后，基于提取的方位AVO梯度属性，解耦提取工区各向同性梯度属性及不同对称面内的各向异性梯度属性；最后，运用提取的不同对称面内的各向异性梯度属性，识别地下水平裂缝位置及其各向异性强度，有效地指导后续地震反演工作。该技术可增强页岩水平裂缝发育位置识别的可靠性，提高页岩储层预测及流体识别的精度。

Description

一种方位叠前地震属性解耦提取方法

技术领域

本发明属于地震资料解释领域，具体的说，涉及一种基于水平裂缝非均质地层的方位叠前地震属性解耦提取方法。

背景技术

页岩油气储层较碳酸盐岩储层埋深较浅，地层的构造应力相对于地层的上覆压力较大，而生成的裂缝将沿着最大主应力、垂直于最小主应力的方向扩展，因此裂缝将平行于地层的层理面，此时垂向裂缝发育极少，而受粘土矿物颗粒定向排列影响的水平裂缝则成为地层各向异性的主要因素。同时，不同方位的水平裂缝造成地层的各向异性强弱程度存在差别，页岩地层将产生方位各向异性特征。常规AVO技术特征分析的基础是Zoeppritz方程的Shuey两项线性近似式，未考虑地层各向异性的影响，提取的AVO梯度属性是不存在方位信息的，各向同性梯度属性与不同对称面内的各向异性梯度属性是耦合的，无法突出地层各向异性的影响，为此，发展用于解耦提取实际地层的各向同性梯度及各向异性梯度属性，综合考虑方位各向异性的水平裂缝页岩地层特征属性分析的AVO技术便显得尤为重要。

发明内容

本发明的目的在于至少解决现有技术中存在的技术问题之一。

为此，本发明的目的在于，提供一种基于水平裂缝非均质地层的方位叠前地震属性解耦提取方法。

为解决上述技术问题，本发明采用以下技术方案予以实现：

为提高AVO技术的可靠性和稳定性，增强页岩地层水平裂缝强度的识别，并指导地下裂缝位置的预测和裂缝描述，本发明以方位各向异性理论为基础，综合成像测井资料和岩心解释结果，从水平裂缝非均质地层假设出发，提出一种直接从多方位叠前角度道集中提取不同对称面内的各向异性梯度属性的新方法。

本发明是一种以方位各向异性理论为基础，综合成像测井资料和岩心解释结果，基于水平裂缝非均质地层假设，直接从多方位叠前角度道集中提取不同对称面内的各向异性梯度属性而进行裂缝发育位置的识别方法。

本发明所采用的技术方案是：

一种方位叠前地震属性解耦提取方法，包括如下步骤：

首先，综合成像测井资料和岩心解释结果，基于水平裂缝非均质地层假设推导由水平裂缝引起的方位各向异性AVO两项近似方程如下：

R_PP(θ,φ)＝P+(G_iso+Γ_xcos²φ+Γ_ysin²φ)sin²θ

其中：R_PP(θ,φ)为由水平裂缝引起的方位各向异性纵波反射系数，P为截距属性，G_iso为各向同性梯度属性，Γ_x与Γ_y分别表征不同对称面内的各向异性梯度属性，

$P=\frac{1}{2}(\frac{\mathrm{\Δ\α}}{\stackrel{\‾}{\mathrm{\α}}}+\frac{\mathrm{\Δ\ρ}}{\stackrel{\‾}{\mathrm{\ρ}}}),$

$G_{\mathrm{iso}}=\frac{1}{2}=[\frac{\mathrm{\Δ\α}}{\stackrel{\‾}{\mathrm{\α}}}-{\left(\frac{2\stackrel{\‾}{\mathrm{\β}}}{\stackrel{\‾}{\mathrm{\α}}}\right)}^2\frac{\mathrm{\Δ}\left(\mathrm{\ρ}{\mathrm{\β}}^2\right)}{\left(\mathrm{\ρ}{\mathrm{\β}}^2\right)}],$

${\mathrm{\Γ}}_x=\frac{1}{2}(\mathrm{\Δ}{\mathrm{\δ}}_x-8{\left(\frac{\stackrel{\‾}{\mathrm{\β}}}{\stackrel{\‾}{\mathrm{\α}}}\right)}^2\mathrm{\Δ}{\mathrm{\γ}}_x),$

${\mathrm{\Γ}}_y=\frac{1}{2}(\mathrm{\Δ}{\mathrm{\δ}}_y-8{\left(\frac{\stackrel{\‾}{\mathrm{\β}}}{\stackrel{\‾}{\mathrm{\α}}}\right)}^2\mathrm{\Δ}{\mathrm{\γ}}_y).$

其中，θ为入射角，φ为方位角，α、β与ρ分别为各向同性背景部分的纵、横波速度及密度；为各向同性背景部分的横纵波速度比；与分别为各向同性背景部分的纵波反射系数、密度反射系数及剪切模量反射系数；Δδ_x与Δδ_y、Δγ_x与Δγ_y分别表征不同对称面内的上下两层介质Thomsen各向异性参数差值。

首先，基于推导的方位各向异性AVO两项近似方程，提取不同方位角的AVO截距和梯度属性；

然后，使用提取的方位AVO梯度属性直接从多方位叠前角度道集中解耦提取各向同性梯度属性及不同对称面内的各向异性梯度属性；

最后，运用提取的不同对称面内的各向异性梯度属性及其强弱差异，进行水平裂缝发育位置及各向异性强度的识别，指导地下裂缝位置的预测和裂缝强度的描述。

进一步地，小角度入射假设下，各向异性AVO近似公式写成由截距、各向同性梯度以及不同对称面内的各向异性梯度表征的形式。

进一步地，基于推导的水平裂缝引起的方位各向异性AVO近似公式，并假设已知裂缝发育方向为方位角0度，使用提取的方位AVO梯度属性，解耦提取各向同性梯度及不同对称面内的各向异性梯度，用于页岩储层水平裂缝位置的预测和各向异性强度的描述。

据此，可以准确判定地下水平裂缝发育较强的位置，为页岩储层水平裂缝位置的预测和裂缝强度的描述提供了一种较为可靠的研究思路。

该方法可增强页岩水平裂缝发育位置识别的可靠性，提高页岩储层预测及流体识别的精度。

具体分为下面两部分：

1)页岩地层不同方位AVO截距及梯度属性提取

基于推导的受水平裂缝影响的方位各向异性两项AVO近似公式，对已处理的页岩工区不同方位叠前角度道集提取不同方位的AVO截距和梯度属性。

2)页岩地层各向同性梯度及各向异性梯度属性解耦提取

假设已知裂缝发育方向为方位角0度，使用提取的方位AVO梯度属性解耦提取各向同性梯度及不同对称面内的各向异性梯度属性，用于识别地下水平裂缝的各向异性强度。

下面详细论述该方法的基本原理：

1)页岩地层不同方位AVO截距及梯度属性提取

根据弱各向异性近似理论，各向异性介质中反射透射系数矩阵可以描述为各向同性背景介质反射透射系数矩阵与各向异性扰动项之和的形式。Psencik和Martins(2001)和Bachrach等(2014)根据弱各向异性理论，详细研究了不同方位水平裂缝地层中地震波反射透射近似关系式，同样表示成各向同性背景部分反射系数和各向异性扰动项反射系数加和的形式，在此基础上，假设已知裂缝发育方向为方位角0度，对其形式进行了改写：

$\begin{array}{c}{R}_{\mathrm{PP}}(\mathrm{\θ},\mathrm{\φ})=R_{\mathrm{PP}}^{\mathrm{iso}}\left(\mathrm{\θ}\right)+R_{\mathrm{PP}}^{\mathrm{ani}}(\mathrm{\θ},\mathrm{\φ})\\ =R_{\mathrm{PP}}^{\mathrm{iso}}\left(\mathrm{\θ}\right)+\left\{\begin{array}{c}{[b}_1{\cos }^2\mathrm{\φ}+b_2{\sin }^2\mathrm{\φ}]{\sin }^2\mathrm{\θ}+\\ \left[\begin{array}{c}{b}_3{\cos }^4\mathrm{\φ}+b_4{\sin }^4\mathrm{\φ}\\ +b_5{\sin }^2\mathrm{\φ}{\cos }^2\mathrm{\φ}\end{array}\right]{\sin }^2\mathrm{\θ}{\tan }^2\mathrm{\θ}\end{array}\right\}\end{array}---\left(1\right)$

其中，表征各向同性背景部分反射系数，表征各向异性扰动项反射系数，θ为入射角，φ为方位角，定义为假设0度方位线与裂缝对称轴平行时观测测线的方位角，或者直接表示观测测线方位角与裂缝对称轴的夹角，Γ_x与Γ_y分别表征不同对称面内的各向异性梯度，Δε_x与Δε_y、Δδ_x、Δδ_y与Δδ_z、Δγ_x与Δγ_y分别表征不同对称面内的上下两层介质Thomsen各向异性参数差值，与表征各向同性背景部分的纵、横波速度均值。

$b_1={\mathrm{\Γ}}_x=\frac{1}{2}(\mathrm{\Δ}{\mathrm{\δ}}_x-8{\left(\frac{\stackrel{\‾}{\mathrm{\β}}}{\stackrel{\‾}{\mathrm{\α}}}\right)}^2\mathrm{\Δ}{\mathrm{\γ}}_x),$

$b_2={\mathrm{\Γ}}_y=\frac{1}{2}(\mathrm{\Δ}{\mathrm{\δ}}_y-8{\left(\frac{\stackrel{\‾}{\mathrm{\β}}}{\stackrel{\‾}{\mathrm{\α}}}\right)}^2\mathrm{\Δ}{\mathrm{\γ}}_y),$

$b_3=\frac{\mathrm{\Δ}{\mathrm{\ϵ}}_x}{2},$

$b_4=\frac{\mathrm{\Δ}{\mathrm{\ϵ}}_y}{2},$

$b_5=\frac{\mathrm{\Δ}{\mathrm{\δ}}_z}{2}.$

公式中各个弱各向异性参数的定义如下：

${\mathrm{\ϵ}}_x=\frac{A_{11}-{\mathrm{\α}}^2}{2{\mathrm{\α}}^2},{\mathrm{\ϵ}}_y=\frac{A_{22}-{\mathrm{\α}}^2}{2{\mathrm{\α}}^2},$

${\mathrm{\γ}}_x=\frac{A_{55}-{\mathrm{\β}}^2}{{\mathrm{\β}}^2},{\mathrm{\γ}}_y=\frac{A_{44}-{\mathrm{\β}}^2}{{\mathrm{\β}}^2},$

${\mathrm{\δ}}_x=\frac{A_{13}+2A_{55}-{\mathrm{\α}}^2}{{\mathrm{\α}}^2},{\mathrm{\δ}}_y=\frac{A_{23}+2A_{44}-{\mathrm{\α}}^2}{{\mathrm{\α}}^2},$

${\mathrm{\δ}}_z=\frac{A_{12}+2A_{66}-{\mathrm{\α}}^2}{2{\mathrm{\α}}^2}.$

其中，各个弱各向异性参数定义式中C_ij为OA介质的弹性系数。

当入射角小于30度时，舍掉sin²θtan²θ项。公式(1)化简为：

R_PP(θ,φ)＝P+(G_iso+Γ_xcos²φ+Γ_ysin²φ)sin²θ   (2)

其中：P为截距属性，G_iso为各向同性梯度属性，Γ_x与Γ_y分别表征不同对称面内的各向异性梯度属性，

$P=\frac{1}{2}(\frac{\mathrm{\Δ\α}}{\stackrel{\‾}{\mathrm{\α}}}+\frac{\mathrm{\Δ\ρ}}{\stackrel{\‾}{\mathrm{\ρ}}}),$

$G_{\mathrm{iso}}=\frac{1}{2}=[\frac{\mathrm{\Δ\α}}{\stackrel{\‾}{\mathrm{\α}}}-{\left(\frac{2\stackrel{\‾}{\mathrm{\β}}}{\stackrel{\‾}{\mathrm{\α}}}\right)}^2\frac{\mathrm{\Δ}\left(\mathrm{\ρ}{\mathrm{\β}}^2\right)}{\left(\mathrm{\ρ}{\mathrm{\β}}^2\right)}],$

${\mathrm{\Γ}}_x=\frac{1}{2}(\mathrm{\Δ}{\mathrm{\δ}}_x-8{\left(\frac{\stackrel{\‾}{\mathrm{\β}}}{\stackrel{\‾}{\mathrm{\α}}}\right)}^2\mathrm{\Δ}{\mathrm{\γ}}_x),$

${\mathrm{\Γ}}_y=\frac{1}{2}(\mathrm{\Δ}{\mathrm{\δ}}_y-8{\left(\frac{\stackrel{\‾}{\mathrm{\β}}}{\stackrel{\‾}{\mathrm{\α}}}\right)}^2\mathrm{\Δ}{\mathrm{\γ}}_y).$

其中，θ为入射角，φ为方位角，α、β与ρ分别为各向同性背景部分的纵、横波速度及密度；为各向同性背景部分的横纵波速度比；与分别为各向同性背景部分的纵波反射系数、密度反射系数及剪切模量反射系数。

公式(2)也可基于Shuey两项AVO近似公式，写成简单的两项AVO近似形式：

R_PP(θ,φ)＝P+G(φ)sin²θ   (3)

其中，P为截距属性，G(φ)为方位梯度属性，其表示形式如下：

G(φ)＝G_iso+Γ_xcos²φ+Γ_ysin²φ   (4)

基于公式(3)，便可提取截距属性及不同方位的梯度属性。

2)页岩地层各向同性梯度及各向异性梯度属性解耦提取

Γ_x和Γ_y为不同对称面内的各向异性梯度。假设已知裂缝发育方向为方位角0度，那么不同方向各向异性梯度Γ_x和Γ_y的计算方法如下：

G₁＝G_iso+Γ_xcos²(φ₁)+Γ_ysin²(φ₁)

G₂＝G_iso+Γ_xcos²(φ₂)+Γ_ysin²(φ₂)    (5)

G₃＝G_iso+Γ_xcos²(φ₃)+Γ_ysin²(φ₃)

先将上式两两相减，

B₁＝Γ_x[cos²(φ₂)-cos²(φ₁)]+Γ_y[sin²(φ₂)-sin²(φ₁)]

        (6)

B₂＝Γ_x[cos²(φ₃)-cos²(φ₂)]+Γ_y[sin²(φ₃)-sin²(φ₂)]

因此，求得不同方向的各向异性梯度为：

${\mathrm{\Γ}}_x=\frac{B_2[{\sin }^2\left({\mathrm{\φ}}_2\right)-{\sin }^2\left({\mathrm{\φ}}_1\right)]-B_1[{\sin }^2\left({\mathrm{\φ}}_3\right)-{\sin }^2\left({\mathrm{\φ}}_2\right)]}{[{\cos }^2\left({\mathrm{\φ}}_3\right)-{\cos }^2\left({\mathrm{\φ}}_2\right)][{\sin }^2\left({\mathrm{\φ}}_2\right)-{\sin }^2\left({\mathrm{\φ}}_1\right)]-[{\cos }^2\left({\mathrm{\φ}}_2\right)-{\cos }^2\left({\mathrm{\φ}}_1\right)][{\sin }^2\left({\mathrm{\φ}}_3\right)-{\sin }^2\left({\mathrm{\φ}}_2\right)]}---\left(7\right)$

${\mathrm{\Γ}}_y=\frac{B_2[{\cos }^2\left({\mathrm{\φ}}_2\right)-{\cos }^2\left({\mathrm{\φ}}_1\right)]-B_1[{\cos }^2\left({\mathrm{\φ}}_3\right)-{\cos }^2\left({\mathrm{\φ}}_2\right)]}{[{\sin }^2\left({\mathrm{\φ}}_3\right)-{\sin }^2\left({\mathrm{\φ}}_2\right)][{\cos }^2\left({\mathrm{\φ}}_2\right)-{\cos }^2\left({\mathrm{\φ}}_1\right)]-[{\sin }^2\left({\mathrm{\φ}}_2\right)-{\sin }^2\left({\mathrm{\φ}}_1\right)][{\cos }^2\left({\mathrm{\φ}}_3\right)-{\cos }^2\left({\mathrm{\φ}}_2\right)]}---\left(8\right)$

在此基础上，各向同性梯度便可简单的计算，形式如下：

G_iso＝G₁-Γ_xcos²(φ₁)-Γ_ysin²(φ₁)    (9)

与现有技术相比，本发明的优点和积极效果是：

本发明综合成像测井资料和岩心解释结果，从水平裂缝非均质地层假设出发，解耦提取各向同性梯度及各向异性梯度属性，指导地下水平裂缝位置的预测和水平裂缝强度的描述。

附图说明

图1为基于水平裂缝非均质地层的方位叠前地震属性解耦提取流程图；

图2a-f为提取的AVO截距P属性剖面及方位AVO梯度G(φ)属性剖面图，其中图2a-b、图2c-d和图2e-f分别提取的方位1、方位2和方位3的AVO截距和梯度属性图；

图3a-c为解耦提取的各向同性梯度属性及不同对称面内的各向异性梯度属性图，其中图3a为解耦提取的各向异性梯度Γ_x属性剖面，图3b解耦提取的各向异性梯度为Γ_y属性剖面，图3c为解耦提取的各向同性梯度G_iso属性剖面；

图4a-c为整个研究区域的各向同性梯度切片及不同对称面内的各向异性梯度切片，其中图4a为解耦提取的各向异性梯度Γ_x属性切片，图4b为解耦提取的各向异性梯度Γ_y属性切片，图4c为解耦提取的各向同性梯度G_iso属性切片；

图5a-b为不同对称面内的各向异性梯度差属性剖面及其切片，其中图5a为不同对称面内的各向异性梯度差Γ_y-Γ_x属性剖面，图5b为同对称面内的各向异性梯度差Γ_y-Γ_x属性切片；

图6为某个方向各向异性梯度与裂缝密度的测井解释结果。

具体实施方式

下面结合附图和具体实施方式对本发明作进一步详细的说明。

本发明方位叠前地震属性解耦提取方法，采用以下步骤，

首先，在充分分析水平裂缝方位各向异性的基础上，对实际水平裂缝型页岩工区不同方位叠前角度道集进行处理，基于推导的由水平裂缝引起的方位各向异性AVO两项近似方程，提取不同方位角的AVO截距和梯度属性；

然后，基于提取的方位AVO梯度属性，解耦提取工区各向同性梯度及不同对称面内的各向异性梯度；

最后，运用提取的不同对称面内的各向异性梯度属性及其差异属性，圈定较强的各向异性区域，识别地下水平裂缝位置及其各向异性强度，指导地下裂缝位置的预测和裂缝强度的描述，有效地指导后续地震反演工作。

参考图1所示，本发明的实施流程如下：

首先从不同方位水平裂缝引起各向异性的两项AVO近似公式推导出发，对经过处理的实际水平裂缝型页岩工区不同方位叠前角度道集提取AVO截距属性及方位AVO梯度属性；

然后，假设已知裂缝发育方向为方位角0度，使用提取的方位AVO梯度属性解耦提取各向同性梯度及不同对称面内的各向异性梯度属性；

最后，基于提取的不同对称面内的各向异性梯度属性，指导地下裂缝位置的预测和裂缝强度的描述。

下面结合某页岩工区实际地震数据的实例(附图2a-f至附图6)对本发明的具体实施方式做进一步说明。

1)页岩地层不同方位AVO截距及梯度属性提取

以某页岩工区实际地震数据为例，进行AVO截距属性及方位AVO梯度属性的提取试算。图2为提取的不同方位的AVO截距和梯度属性。

图2a-f分别提取的方位1、方位2和方位3的AVO截距和梯度属性。从图中可以看出，不同方位的AVO截距项基本相同。

2)页岩地层各向同性梯度及各向异性梯度属性解耦提取

根据前面研究的基于水平裂缝非均质地层的方位叠前地震属性解耦提取的流程，利用方位AVO梯度属性解耦提取的该条测线各向同性梯度及不同对称面内的各向异性梯度结果如图3所示。

图3a、图3b和图3c分别为G_iso剖面、Γ_x剖面和Γ_y剖面，从色标数值范围变化上可以看出，x对称面内的各向异性梯度值Γ_x略小于y对称面内的各向异性梯度值Γ_y，且y对称面内的各向异性梯度剖面显示的层位明显多于x对称面内。

图4为整个研究区域的不同对称面内的各向异性梯度切片，其中图4a为估测的各向异性梯度Γ_x切片，图4b为估测的各向异性梯度Γ_y切片。

从图中可以清晰地看出，如果x对称面内的各向异性梯度Γ_x呈现低值，那么y对称面内的各向异性梯度Γ_y同样呈现低值，两者的差别表现在Γ_x值略小于Γ_y值。

图5a-b为不同对称面内的各向异性梯度差属性剖面及其切片，其中图5a为不同对称面内的各向异性梯度差Γ_y-Γ_x属性剖面，图5b为同对称面内的各向异性梯度差Γ_y-Γ_x属性切片。

通过对比不同对称面内的各向异性梯度变化，可定量解释不同对称面内的属性差异，圈定较强的各向异性区域，能够有效地指导后续地震反演工作。

图6为某个方向各向异性梯度与裂缝密度的测井解释结果。

由图中可以看出，各向异性梯度的峰值和裂缝密度的峰值存在良好的对应关系，因而可以认为各向异性梯度能用来指示水平裂缝型储层，指导地下水平裂缝位置的预测和水平裂缝强度的描述。

以上所述，仅是本发明的较佳实施例而已，并非是对本发明作其它形式的限制，任何熟悉本专业的技术人员可能利用上述揭示的技术内容加以变更或改型为等同变化的等效实施例。但是凡是未脱离本发明技术方案内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与改型，仍属于本发明技术方案的保护范围。

Claims (3)

1.一种方位叠前地震属性解耦提取方法，其特征在于包括如下步骤：

综合成像测井资料和岩心解释结果，基于水平裂缝非均质地层假设推导由水平裂缝引起的方位各向异性AVO两项近似方程如下：

R_PP(θ,φ)＝P+(G_iso+Γ_xcos²φ+Γ_ysin²φ)sin²θ

其中：R_PP(θ,φ)为由水平裂缝引起的方位各向异性纵波反射系数，P为截距属性，G_iso为各向同性梯度属性，Γ_x与Γ_y分别表征不同对称面内的各向异性梯度属性，

$P=\frac{1}{2}(\frac{\mathrm{\Δ\α}}{\stackrel{\‾}{\mathrm{\α}}}+\frac{\mathrm{\Δ\ρ}}{\stackrel{\‾}{\mathrm{\ρ}}}),$

$G_{\mathrm{iso}}=\frac{1}{2}[\frac{\mathrm{\Δ\α}}{\stackrel{\‾}{\mathrm{\α}}}-{\left(\frac{2\stackrel{\‾}{\mathrm{\β}}}{\stackrel{\‾}{\mathrm{\α}}}\right)}^2\frac{\mathrm{\Δ}\left({\mathrm{\ρ\β}}^2\right)}{\stackrel{\‾}{\left({\mathrm{\ρ\β}}^2\right)}}],$

${\mathrm{\Γ}}_x=\frac{1}{2}({\mathrm{\Δ\δ}}_x-8{\left(\frac{\stackrel{\‾}{\mathrm{\β}}}{\stackrel{\‾}{\mathrm{\α}}}\right)}^2\mathrm{\Δ}{\mathrm{\γ}}_x),$

${\mathrm{\Γ}}_y=\frac{1}{2}({\mathrm{\Δ\δ}}_y-8{\left(\frac{\stackrel{\‾}{\mathrm{\β}}}{\stackrel{\‾}{\mathrm{\α}}}\right)}^2\mathrm{\Δ}{\mathrm{\γ}}_y),$

其中，θ为入射角，φ为方位角，α、β与ρ分别为各向同性背景部分的纵、横波速度及密度；为各向同性背景部分的横纵波速度比；与分别为各向同性背景部分的纵波反射系数、密度反射系数及剪切模量反射系数；△δ_x与△δ_y、△γ_x与△γ_y分别表征不同对称面内的上下两层介质Thomsen各向异性参数差值；

首先，基于推导的方位各向异性AVO两项近似方程，提取不同方位角的AVO截距和梯度属性；

然后，使用提取的方位AVO梯度属性直接从多方位叠前角度道集中解耦提取各向同性梯度属性及不同对称面内的各向异性梯度属性；

最后，运用提取的不同对称面内的各向异性梯度属性及其差异属性，进行水平裂缝发育位置及各向异性强度的识别，指导地下裂缝位置的预测和裂缝强度的描述。

2.根据权利要求1所述方位叠前地震属性解耦提取方法，其特征在于：小角度入射假设下，各向异性AVO近似公式写成由截距、各向同性梯度以及不同对称面内的各向异性梯度表征的形式。

3.根据权利要求1所述方位叠前地震属性解耦提取方法，其特征在于：基于推导的水平裂缝引起的方位各向异性AVO近似公式，并假设已知裂缝发育方向为方位角0度，使用提取的方位AVO梯度属性，解耦提取各向同性梯度及不同对称面内的各向异性梯度，用于页岩储层水平裂缝位置的预测和各向异性强度的描述。