CN113075727B - 基于振幅方位导数的储层裂缝分布特征确定方法及系统 - Google Patents

Abstract

本发明提供了一种基于振幅方位导数的储层裂缝分布特征确定方法及系统，所述方法包括：根据三维地震资料得到储层裂缝目的层同相轴的振幅方位导数；根据振幅方位导数构建反演模型并得到拟各向异性参数，并将拟各向异性参数转换为储层裂缝参数；提取储层裂缝参数的沿层切片和剖面进行储层裂缝分布特征分析，本发明可提高储层裂缝分布特征的预测精度。

Description

基于振幅方位导数的储层裂缝分布特征确定方法及系统

技术领域

本发明涉及石油勘探技术领域，尤其涉及一种基于振幅方位导数的储层裂缝分布特征确定方法及系统。

背景技术

裂缝型油气藏的产量占全世界石油天然气总产量的一半以上，在国内的油气勘探中，裂缝型油气藏是一个重要的勘探领域，如塔里木盆地奥陶系碳酸盐岩裂缝型油气藏，四川致密砂岩裂缝型油气藏。裂缝既是油、气、水等地下流体的运移通道，也是流体的储集空间，裂缝能极大地提升储层的渗透率。裂缝型储层描述包括预测裂缝分布特征和检测裂隙充填物，也就是探测开启并充填流体的有效裂缝。

当没有CRP道集(共反射点道集)或者资料信噪比很低时，业界一般用AVD方法预测裂缝。主要原理是纵波属性在极坐标中表现为椭圆形状，常用的属性有反射振幅、动校正速度、AVO梯度、层间旅行时和品质因子。利用椭圆拟合计算地震属性的长短轴，用椭圆扁率表征裂缝发育程度，用椭圆长轴方向表示裂缝发育方向。肖鹏飞[21]研究了倾角对裂缝密度反演的影响，当裂缝倾角小于80°时，椭圆拟合预测的裂缝密度存在较大偏差。且方位反射振幅在极坐标中不全是椭圆，也可能是“花生”形状，此时，椭圆拟合预测裂缝误差较大。

发明内容

本发明的一个目的在于提供一种基于振幅方位导数的储层裂缝分布特征确定方法，提高储层裂缝分布特征的预测精度。本发明的另一个目的在于提供一种基于振幅方位导数的储层裂缝分布特征确定系统。本发明的再一个目的在于提供一种计算机设备。本发明的还一个目的在于提供一种可读介质。

为了达到以上目的，本发明一方面公开了一种基于振幅方位导数的储层裂缝分布特征确定方法，包括：

根据三维地震资料得到储层裂缝目的层同相轴的振幅方位导数；

根据振幅方位导数构建反演模型并得到拟各向异性参数，并将拟各向异性参数转换为储层裂缝参数；

提取储层裂缝参数的沿层切片和剖面进行储层裂缝分布特征分析。

优选的，所述根据三维地震资料得到储层裂缝目的层同相轴的振幅方位导数具体包括：

根据三维地震资料得到地层分方位的叠前分方位道集；

根据所述叠前分方位道集提取储层裂缝目的层的振幅和方位角；

根据所述振幅和方位角计算所述目的层振幅的方位导数。

优选的，所述根据振幅方位导数构建反演模型并得到拟各向异性参数具体包括：

根据三维地震资料得到储层裂缝目的层同相轴的方位角，并根据所述方位角得到系数矩阵；

根据所述振幅的方位导数和系数矩阵构建反演模型；

求解所述反演模型得到拟各向异性参数。

优选的，所述将拟各向异性参数转换为储层裂缝参数具体包括：

对所述拟各向异性参数进行道积分得到储层裂缝参数。

优选的，所述提取储层裂缝参数的沿层切片和剖面进行储层裂缝分布特征分析具体包括：

提取所述储层裂缝参数的沿层切片并进行裂缝的横向分布特征；

提取所述储层裂缝参数的剖面并进行裂缝的纵向分布特征。

本发明还公开了一种基于振幅方位导数的储层裂缝分布特征确定系统，包括：

方位导数确定单元，用于根据三维地震资料得到储层裂缝目的层同相轴的振幅方位导数；

参数确定单元，用于根据振幅方位导数构建反演模型并得到拟各向异性参数，并将拟各向异性参数转换为储层裂缝参数；

分布特征分析单元，用于提取储层裂缝参数的沿层切片和剖面进行储层裂缝分布特征分析。

优选的，所述方位导数确定单元具体包括：

原始数据处理模块，用于根据三维地震资料得到地层分方位的叠前分方位道集；

振幅参数确定模块，用于根据所述叠前分方位道集提取储层裂缝目的层的振幅和方位角；

方位导数确定模块，用于根据所述振幅和方位角计算所述目的层振幅的方位导数。

优选的，所述参数确定单元具体包括：

系数矩阵确定模块，用于根据三维地震资料得到储层裂缝目的层同相轴的方位角，并根据所述方位角得到系数矩阵；

反演模型构建模块，用于根据所述振幅的方位导数和系数矩阵构建反演模型；

反演模型求解模块，用于求解所述反演模型得到拟各向异性参数。

优选的，所述参数确定单元进一步包括：

参数确定模块，用于对所述拟各向异性参数进行道积分得到储层裂缝参数。

优选的，所述分布特征分析单元具体包括：

横向分析模块，用于提取所述储层裂缝参数的沿层切片并进行裂缝的横向分布特征；

纵向分析模块，用于提取所述储层裂缝参数的剖面并进行裂缝的纵向分布特征。

本发明还公开了一种计算机设备，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，

所述处理器执行所述程序时实现如上所述方法。

本发明还公开了一种计算机可读介质，其上存储有计算机程序，

该程序被处理器执行时实现如上所述方法。

本发明通过三维地震资料得到目的层同相轴的振幅方位导数，根据振幅方位导数构建反演模型，求解反演模型可得到拟各向异性参数，进一步根据拟各向异性参数得到储层裂缝参数，根据储层裂缝参数可进行储层裂缝分布特征分析。本发明通过研究裂缝分布特征的有关因素，通过振幅的方位导数对储层裂缝的分布特征进行分析，可有效提升裂缝分布特征预测的准确度，从而有利于确定位置更优的储层裂缝勘探位置，提升油气勘探效果。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1示出本发明基于振幅方位导数的储层裂缝分布特征确定方法一个具体实施例的流程图；

图2示出本发明基于振幅方位导数的储层裂缝分布特征确定方法一个具体实施例裂缝介质反射振幅曲线；

图3示出本发明基于振幅方位导数的储层裂缝分布特征确定方法一个具体实施例S100的流程图；

图4示出本发明基于振幅方位导数的储层裂缝分布特征确定方法一个具体实施例裂缝介质振幅的方位导数曲线；

图5示出本发明基于振幅方位导数的储层裂缝分布特征确定方法一个具体实施例振幅的方位导数曲线与裂缝密度关系图；

图6示出本发明基于振幅方位导数的储层裂缝分布特征确定方法一个具体实施例含噪音的目的层方位角道集；

图7示出本发明基于振幅方位导数的储层裂缝分布特征确定方法一个具体实施例振幅方位导数与反演结果拟合曲线；

图8示出本发明基于振幅方位导数的储层裂缝分布特征确定方法一个具体实施例S200得到拟各向异性参数的流程图；

图9示出本发明基于振幅方位导数的储层裂缝分布特征确定方法一个具体实施例S200得到储层裂缝参数的流程图；

图10示出本发明基于振幅方位导数的储层裂缝分布特征确定方法一个具体实施例S300的流程图；

图11示出本发明基于振幅方位导数的储层裂缝分布特征确定方法一个具体实施例储层裂缝参数沿层切片和剖面图；

图12示出本发明基于振幅方位导数的储层裂缝分布特征确定系统一个具体实施例的结构图；

图13示出本发明基于振幅方位导数的储层裂缝分布特征确定系统一个具体实施例方位导数确定单元的结构图；

图14示出本发明基于振幅方位导数的储层裂缝分布特征确定系统一个具体实施例参数确定单元的结构图；

图15示出本发明基于振幅方位导数的储层裂缝分布特征确定系统一个具体实施例参数确定单元包括参数确定模块的结构图；

图16示出本发明基于振幅方位导数的储层裂缝分布特征确定系统一个具体实施例分布特征分析单元的结构图；

图17示出适于用来实现本发明实施例的计算机设备的结构示意图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

根据本发明的一个方面，本实施例公开了一种基于振幅方位导数的储层裂缝分布特征确定方法。如图1所示，本实施例中，所述方法包括：

S100：根据三维地震资料得到储层裂缝目的层同相轴的振幅方位导数。

S200：构建反演模型并得到拟各向异性参数，并将拟各向异性参数转换为储层裂缝参数。

S300：根据储层裂缝参数的沿层切片和剖面得到储层裂缝分布。

可以理解的是，含竖向平行排列裂缝的介质是典型的HTI介质，其反射系数不仅随入射角变化(AVO特征)，还随方位角变化(各向异性特征)。在一个具体例子中，依据Ruger反射系数公式，储层裂缝目的层同相轴的振幅可表示为含有四阶项的AVO公式：

R(θ，φ)＝A+Bsin²θ+Csin²θtan²θ

＝A^iso+(B^iso+B^anicos²φ)sin²θ+(C^iso+C^ani1cos⁴φ+C^ani2sin²φcos²φ)sin²θtan²θ⁽¹⁾

其中，R是反射系数，θ为入射角。φ为测线与裂缝面法向的夹角。A是AVO截距，B和C分别是AVO梯度和曲率，包含各向异性修正量。

A＝A^iso

B＝B^iso+B^ani cos²φ

C＝C^iso+C^ani1 cos⁴φ+C^ani2 sin²φcos²φ

其中，上标^iso和^ani表示对相关参数的各向异性修正。α，β，Z，G分别是纵波速度、横波速度、纵波阻抗和剪切模量。ε^v，δ^v，γ为各向异性参数。ˉ表示裂缝上下两层介质间参数的平均，Δ表示上下两层介质间参数的差。

把Ruger公式改写成方位角余弦的和函数，可表示为：

R^iso＝A^iso+B^iso sin²θ+C^iso sin²θtan²θ

其中，d₀是对称轴方向的反射系数，余弦项的系数(拟各向异性参数)d₂和d₄随各向异性参数ε^v、δ^v、γ和入射角θ变化。R^iso是各向同性背景介质的反射系数，g是横纵波速度比的平方。

Ruger的模型中，上层为各向同性介质，下层为各向异性介质。在一个具体例子中，界面参数见表1。下层介质裂缝走向为90°，对称轴方向为0°。观测系统的方位角从0°到360°，入射角等于30°。模型的全方位反射系数在极坐标中呈“花生”状，如图2实线1所示。振幅在裂缝法向达到最大，在裂缝方向处最小。径向表示反射振幅，逆时针表示方位角增大方向。

表1：ISO/HTI两层介质模型参数

则振幅的方位导数(ADA：the azimuthal derivative of the amplitude)，即反射振幅的方位变化率，可表示反射振幅在每个方位的变化规律。

本发明通过三维地震资料得到目的层同相轴的振幅方位导数，根据振幅方位导数构建反演模型，求解反演模型可得到拟各向异性参数，进一步根据拟各向异性参数得到储层裂缝参数，根据储层裂缝参数可进行储层裂缝分布特征分析。本发明通过研究裂缝分布特征的有关因素，通过振幅的方位导数对储层裂缝的分布特征进行分析，可有效提升裂缝分布特征预测的准确度，从而有利于确定位置更优的储层裂缝勘探位置，提升油气勘探效果。

在优选的实施方式中，如图3所示，所述S100可包括：

S110：根据三维地震资料得到地层分方位的叠前分方位道集。

S120：根据所述叠前分方位道集提取储层裂缝目的层的振幅和方位角。

S130：根据所述振幅和方位角计算所述目的层振幅的方位导数。

可以理解的是，在公式(2)的基础上，方位导数R_φ的解析函数是二倍角正弦加四倍角正弦，即公式(3)。

该模型的振幅方位导数如图4所示，呈四叶花瓣形状。极径表示方位导数大小，导数在一三象限为负，如图中曲线1所示；导数在二四象限为正，如图中曲线2所示。方位导数绝对值关于x轴、y轴和原点对称，在裂缝方向和法向等于0。在第一象限内，方位导数绝对值在45°方位达到最大。

裂缝密度e从0增大至0.1时，对应的振幅方位导数如图5所示。可见，目标储层不含裂缝时，反射系数不随方位变化，花瓣归于零点；裂缝密度增大时，反射振幅的方位间差异加剧，四叶花瓣逐步变大。

固定中心方位φ时，方位间隔dφ对应的圆弧长度等于方位步长与半径的乘积，完整花瓣的周长L是弧长在全方位的积分。利用方位导数的对称性，全方位积分等于4倍的第一象限积分。

在第一象限内，方位导数正负同号，所以有：

将上式代入公式(4)，可得四叶花瓣的周长：

L＝4|d₂| (5)

固定中心方位φ时，方位间隔dφ对应的扇区面积等于方位步长与半径平方的乘积的一半，完整花瓣的面积S是扇区面积在全方位的积分，等于π(d² ₂+d² ₄)/2。利用方位导数的对称性，全方位积分等于4倍的第一象限积分。

方位导数在第一象限内正负同号，所以有：

分别求得各项函数的积分：

由于函数sin2φsin4φ关于方位角(φ＝π/4,y＝0)中心对称，正负值抵消，其在第一象限内的积分等于0。

将以上三式代入公式(6)，可得四叶花瓣的面积：

综上，花瓣周长与参数d₂成正比，面积与参数成正比，因此拟各向异性参数d₂和/>可用于描述储层裂缝的分布特征。

在一个具体例子中，可通过地质实验过程中得到的三维地震资料，得到储层裂缝目的层的叠前分方位道集。如图6所示，覆盖次数为n，即道集的道数为n，可提取目的层每一道同相轴的方位角和振幅(φi，Ri)，i＝1,2…n。

由于地震道集的方位离散性，可根据所述振幅和方位角采用一阶中心差分法求取每一道振幅的方位导数如图7所示。

其中，Ri为第i道，φi为第i道方位角，φ下标i-1，i，i+1是道序号。

在优选的实施方式中，如图8所示，所述S200中根据振幅方位导数构建反演模型并得到拟各向异性参数具体包括：

S210：根据三维地震资料得到储层裂缝目的层同相轴的方位角，并根据所述方位角得到系数矩阵。

S220：根据所述振幅的方位导数和系数矩阵构建反演模型。

S230：求解所述反演模型得到拟各向异性参数。

在一个具体例子中，可依据公式(3)建立ADA反演模型(9)，计算拟各向异性参数d₂和d₄。

M_(n×2)x_(2×1)＝R'_φ(n×1) (9)

其中，M为系数矩阵，可根据方位角得到；x为待求参数向量；R'_φ为振幅的方位导数构成的向量。

在求解所述反演模型得到拟各向异性参数时，优选的可设定阻尼因子，然后利用奇异值分解法求解所述反演模型的超定方程组，求解的结果为：

X＝M⁺R'

其中，M⁺为矩阵M的广义逆，可由阻尼奇异值分解法求得。

在优选的实施方式中，如图9所示，S200将拟各向异性参数转换为储层裂缝参数具体可包括：

S240：对所述拟各向异性参数进行道积分得到储层裂缝参数。

可以理解的是，拟各向异性参数d₂是裂缝密度变化量Δe的函数(公式10)。

U_g＝(1-2g)U₃₃+gU₁₁

依据d₂可以求得储层裂缝密度变化量Δe，下式f^-1是公式(10)中函数的反函数：

Δe＝f^-1(d₂)

对裂缝密度变化量Δe进行累计求和得到地下储层的裂缝密度e，则对所述拟各向异性参数进行道积分得到储层裂缝参数e的值。

Δe_j＝f^-1(d_2j)

其中，Δe_j＝e_j-e_j-1；j＝1,2,...i；i＝1,2,...i_max。

在优选的实施方式中，如图10所示，所述S300具体可包括：

S310：提取所述储层裂缝参数的沿层切片并进行裂缝的横向分布特征。

S320：提取所述储层裂缝参数的剖面并进行裂缝的纵向分布特征。

可以理解的是，如图11所示，通过拟各向异性参数可得到目的层的储层裂缝参数，从而可以了解目的层每个区域的地质情况，提取储层裂缝参数的沿层切片即得到同一水平面的储层裂缝参数的分布情况，从而根据储层裂缝参数在水平方面上的变化情况可分析得到储层裂缝的横向分布特征。同理的，提取储层裂缝参数在剖面即得到同一竖直面的储层裂缝参数的分布情况，从而根据储层裂缝参数在竖直方面上的变化情况可分析得到储层裂缝的纵向分布特征。根据提取得到的储层裂缝参数不同沿层切片和剖面，可分析得到整个储层裂缝的分布特征。

基于相同原理，本实施例还公开了一种基于振幅方位导数的储层裂缝分布特征确定系统。如图12所示，本实施例中，所述系统包括方位导数确定单元11、参数确定单元12和分布特征分析单元13。

其中，方位导数确定单元11用于根据三维地震资料得到储层裂缝目的层同相轴的振幅方位导数。

参数确定单元12用于根据振幅方位导数构建反演模型并得到拟各向异性参数，并将拟各向异性参数转换为储层裂缝参数。

分布特征分析单元13用于提取储层裂缝参数的沿层切片和剖面进行储层裂缝分布特征分析。

在优选的实施方式中，如图13所示，所述方位导数确定单元11具体包括原始数据处理模块111、振幅参数确定模块112和方位导数确定模块113。

其中，原始数据处理模块111用于根据三维地震资料得到地层分方位的叠前分方位道集。

振幅参数确定模块112用于根据所述叠前分方位道集提取储层裂缝目的层的振幅和方位角。

方位导数确定模块113用于根据所述振幅和方位角计算所述目的层振幅的方位导数。

在优选的实施方式中，如图14所示，所述参数确定单元12具体包括系数矩阵确定模块121、反演模型构建模块122和反演模型求解模块123。

其中，系数矩阵确定模块121用于根据三维地震资料得到储层裂缝目的层同相轴的方位角，并根据所述方位角得到系数矩阵。

反演模型构建模块122用于根据所述振幅的方位导数和系数矩阵构建反演模型。

反演模型求解模块123用于求解所述反演模型得到拟各向异性参数。

在优选的实施方式中，如图15所示，所述参数确定单元12进一步包括参数确定模块124。参数确定模块124用于对所述拟各向异性参数进行道积分得到储层裂缝参数。

在优选的实施方式中，如图16所示，所述分布特征分析单元13具体包括横向分析模块131和纵向分析模块132。

其中，横向分析模块131用于提取所述储层裂缝参数的沿层切片并进行裂缝的横向分布特征。

纵向分析模块132用于提取所述储层裂缝参数的剖面并进行裂缝的纵向分布特征。

由于该系统解决问题的原理与以上方法类似，因此本系统的实施可以参见方法的实施，在此不再赘述。

本实施例3提供一种电子设备，该电子设备可以是台式计算机、平板电脑及移动终端等，本实施例不限于此。在本实施例中，该电子设备可以参照实施例1的方法的实施及实施例2所述的装置，其内容被合并于此，重复之处不再赘述。

图17为本发明实施例的电子设备600的系统构成的示意框图。如图17所示，该电子设备600可以包括中央处理器100和存储器140；存储器140耦合到中央处理器100。值得注意的是，该图是示例性的；还可以使用其他类型的结构，来补充或代替该结构，以实现电信功能或其他功能。

一实施例中，基于振幅方位导数的储层裂缝分布特征确定功能可以被集成到中央处理器100中。其中，中央处理器100可以被配置为进行如下控制：根据三维地震资料得到储层裂缝目的层同相轴的振幅方位导数；根据振幅方位导数构建反演模型并得到拟各向异性参数，并将拟各向异性参数转换为储层裂缝参数；提取储层裂缝参数的沿层切片和剖面进行储层裂缝分布特征分析。

在另一个实施方式中，基于振幅方位导数的储层裂缝分布特征确定系统可以与中央处理器100分开配置，例如可以将基于振幅方位导数的储层裂缝分布特征确定系统设置为与中央处理器100连接的芯片，通过中央处理器的控制来实现基于振幅方位导数的储层裂缝分布特征确定功能。

如图17所示，该电子设备600还可以包括：通信模块110、输入单元120、音频处理单元130、显示器160、电源170。值得注意的是，电子设备600也并不是必须要包括图17中所示的所有部件；此外，电子设备600还可以包括图17中没有示出的部件，可以参考现有技术。

如图17所示，中央处理器100有时也称为控制器或操作控件，可以包括微处理器或其他处理器装置和/或逻辑装置，该中央处理器100接收输入并控制电子设备600的各个部件的操作。

其中，存储器140，例如可以是缓存器、闪存、硬驱、可移动介质、易失性存储器、非易失性存储器或其它合适装置中的一种或更多种。可储存上述与失败有关的信息，此外还可存储执行有关信息的程序。并且中央处理器100可执行该存储器140存储的该程序，以实现信息存储或处理等。

输入单元120向中央处理器100提供输入。该输入单元120例如为按键或触摸输入装置。电源170用于向电子设备600提供电力。显示器160用于进行图像和文字等显示对象的显示。该显示器例如可为LCD显示器，但并不限于此。

该存储器140可以是固态存储器，例如，只读存储器(ROM)、随机存取存储器(RAM)、SIM卡等。还可以是这样的存储器，其即使在断电时也保存信息，可被选择性地擦除且设有更多数据，该存储器的示例有时被称为EPROM等。存储器140还可以是某种其它类型的装置。存储器140包括缓冲存储器141(有时被称为缓冲器)。存储器140可以包括应用/功能存储部142，该应用/功能存储部142用于存储应用程序和功能程序或用于通过中央处理器100执行电子设备600的操作的流程。

存储器140还可以包括数据存储部143，该数据存储部143用于存储数据，例如联系人、数字数据、图片、声音和/或任何其他由电子设备使用的数据。存储器140的驱动程序存储部144可以包括电子设备的用于通信功能和/或用于执行电子设备的其他功能(如消息传送应用、通讯录应用等)的各种驱动程序。

通信模块110即为经由天线111发送和接收信号的发送机/接收机110。通信模块(发送机/接收机)110耦合到中央处理器100，以提供输入信号和接收输出信号，这可以和常规移动通信终端的情况相同。

基于不同的通信技术，在同一电子设备中，可以设置有多个通信模块110，如蜂窝网络模块、蓝牙模块和/或无线局域网模块等。通信模块(发送机/接收机)110还经由音频处理器130耦合到扬声器131和麦克风132，以经由扬声器131提供音频输出，并接收来自麦克风132的音频输入，从而实现通常的电信功能。音频处理器130可以包括任何合适的缓冲器、解码器、放大器等。另外，音频处理器130还耦合到中央处理器100，从而使得可以通过麦克风132能够在本机上录音，且使得可以通过扬声器131来播放本机上存储的声音。

以上参照附图描述了本发明的优选实施方式。这些实施方式的许多特征和优点根据该详细的说明书是清楚的，因此所附权利要求旨在覆盖这些实施方式的落入其真实精神和范围内的所有这些特征和优点。此外，由于本领域的技术人员容易想到很多修改和改变，因此不是要将本发明的实施方式限于所例示和描述的精确结构和操作，而是可以涵盖落入其范围内的所有合适修改和等同物。

本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。

Claims (12)

1.一种基于振幅方位导数的储层裂缝分布特征确定方法，其特征在于，包括：

根据三维地震资料得到储层裂缝目的层同相轴的振幅方位导数R _Φ ；

其中，

R^iso＝A^iso+B^iso sin²θ+C^iso sin²θtan²θ

其中，R是反射系数，θ为入射角，Φ为测线与裂缝面法向的夹角，A是AVO截距，B和C分别是AVO梯度和曲率，包含各向异性修正量，上标^iso表示对相关参数的各向异性修正，α，β分别是纵波速度、横波速度，d₀是对称轴方向的反射系数，拟各向异性参数d₂和d₄随各向异性参数ε^v、δ^v、γ和入射角θ变化，R^iso是各向同性背景介质的反射系数，g是横纵波速度比的平方，-表示裂缝上下两层介质间参数的平均，Δ表示上下两层介质间参数的差；

根据振幅方位导数构建反演模型并得到拟各向异性参数，并将拟各向异性参数转换为储层裂缝参数；

提取储层裂缝参数的沿层切片和剖面进行储层裂缝分布特征分析。

2.根据权利要求1所述的储层裂缝分布特征确定方法，其特征在于，所述根据三维地震资料得到储层裂缝目的层同相轴的振幅方位导数具体包括：

根据三维地震资料得到地层分方位的叠前分方位道集；

根据所述叠前分方位道集提取储层裂缝目的层的振幅和方位角；

根据所述振幅和方位角计算所述目的层的振幅方位导数。

3.根据权利要求1所述的储层裂缝分布特征确定方法，其特征在于，所述根据振幅方位导数构建反演模型并得到拟各向异性参数具体包括：

根据三维地震资料得到储层裂缝目的层同相轴的方位角，并根据所述方位角得到系数矩阵；

根据所述振幅方位导数和系数矩阵构建反演模型；

求解所述反演模型得到拟各向异性参数。

4.根据权利要求1所述的储层裂缝分布特征确定方法，其特征在于，所述将拟各向异性参数转换为储层裂缝参数具体包括：

对所述拟各向异性参数进行道积分得到储层裂缝参数。

5.根据权利要求1所述的储层裂缝分布特征确定方法，其特征在于，所述提取储层裂缝参数的沿层切片和剖面进行储层裂缝分布特征分析具体包括：

提取所述储层裂缝参数的沿层切片并进行裂缝的横向分布特征；

提取所述储层裂缝参数的剖面并进行裂缝的纵向分布特征。

6.一种基于振幅方位导数的储层裂缝分布特征确定系统，其特征在于，包括：

方位导数确定单元，用于根据三维地震资料得到储层裂缝目的层同相轴的振幅方位导数R _Φ ；

其中，

R^iso＝A^iso+B^iso sin²θ+C^iso sin²θtan²θ

其中，R是反射系数，θ为入射角，Φ为测线与裂缝面法向的夹角，A是AVO截距，B和C分别是AVO梯度和曲率，包含各向异性修正量，上标^iso表示对相关参数的各向异性修正，α，β分别是纵波速度、横波速度，d₀是对称轴方向的反射系数，拟各向异性参数d₂和d₄随各向异性参数ε^v、δ^v、γ和入射角θ变化，R^iso是各向同性背景介质的反射系数，g是横纵波速度比的平方，-表示裂缝上下两层介质间参数的平均，Δ表示上下两层介质间参数的差；

参数确定单元，用于根据振幅方位导数构建反演模型并得到拟各向异性参数，并将拟各向异性参数转换为储层裂缝参数；

分布特征分析单元，用于提取储层裂缝参数的沿层切片和剖面进行储层裂缝分布特征分析。

7.根据权利要求6所述的储层裂缝分布特征确定系统，其特征在于，所述方位导数确定单元具体包括：

原始数据处理模块，用于根据三维地震资料得到地层分方位的叠前分方位道集；

振幅参数确定模块，用于根据所述叠前分方位道集提取储层裂缝目的层的振幅和方位角；

方位导数确定模块，用于根据所述振幅和方位角计算所述目的层的振幅方位导数。

8.根据权利要求6所述的储层裂缝分布特征确定系统，其特征在于，所述参数确定单元具体包括：

系数矩阵确定模块，用于根据三维地震资料得到储层裂缝目的层同相轴的方位角，并根据所述方位角得到系数矩阵；

反演模型构建模块，用于根据所述振幅方位导数和系数矩阵构建反演模型；

反演模型求解模块，用于求解所述反演模型得到拟各向异性参数。

9.根据权利要求6所述的储层裂缝分布特征确定系统，其特征在于，所述参数确定单元进一步包括：

参数确定模块，用于对所述拟各向异性参数进行道积分得到储层裂缝参数。

10.根据权利要求6所述的储层裂缝分布特征确定系统，其特征在于，所述分布特征分析单元具体包括：

横向分析模块，用于提取所述储层裂缝参数的沿层切片并进行裂缝的横向分布特征；

纵向分析模块，用于提取所述储层裂缝参数的剖面并进行裂缝的纵向分布特征。

11.一种计算机设备，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，

所述处理器执行所述程序时实现如权利要求1-5任一项所述方法。

12.一种计算机可读介质，其上存储有计算机程序，其特征在于，

该程序被处理器执行时实现如权利要求1-5任一项所述方法。