地震纵波多方位属性椭圆拟合预测裂缝的方法及装置制造方法
【专利摘要】本发明实施例提供了一种地震纵波多方位属性椭圆拟合预测裂缝的方法及装置，其中，该方法包括：将至少包括三个方位的地震属性数据按照测线方位进行中心对称处理，生成具有完备方位的地震属性数据；对生成的具有完备方位的地震属性数据进行椭圆拟合；将椭圆扁率确定为裂缝发育强度。本发明实施例确保了椭圆的中心点位于极坐标的原点，将纵波预测裂缝所需要反演的参数从五个降低到了三个，从而降低了噪音对裂缝特征反演的干扰；同时，将椭圆扁率确定为裂缝发育强度，比现有技术中用长短轴之比描述裂缝发育强度更加合理。
【专利说明】地震纵波多方位属性椭圆拟合预测裂缝的方法及装置
【技术领域】
[0001]本发明涉及石油勘探【技术领域】，特别涉及一种地震纵波多方位属性椭圆拟合预测裂缝的方法及装置。
【背景技术】
[0002]裂缝是岩石受力发生破裂、沿破裂面两侧的岩石没有发生明显位移的一种断裂构造现象。裂缝既是油、气、水等地下流体的运移通道，也是流体的储集空间，裂缝能极大地提升储层的渗透率。有效裂缝是指开启并充填流体的裂缝，有效裂缝的探测对于识别油气运移通道，寻找有效储集空间具有重要的指导意义。
[0003]目前比较成熟的裂缝预测技术是利用相干和曲率等叠后不连续属性、岩心和成像测井来预测裂缝，这些方法存在较多不足。相干反演能够反映小断裂或断裂发育带，岩性突变带，但主要是用于断层边缘的检测，并且检测所有的断裂体系。当地下介质受多期次构造运动或成岩作用时，早期的裂缝可能被泥质或者方解石等所充填而成为无效裂缝，而这种无效裂缝在相干切片上仍能显示，相干反演技术很难预测出有效裂缝的发育特征。采集岩心和成像测井成本高昂，导致资料数量少，只能在井点特定井段处检测裂缝，无法在横向大面积和纵向全深度预测裂缝。而地震纵波宽方位属性能够检测有效裂缝引起的各向异性特征，可大规模 检测裂缝，其中，纵横比大于0.5为宽方位。
[0004]各向异性介质的物理性质(如地震波的速度，振幅等)随测线方位改变而变化，因此，利用方位角信息，采用地震纵波宽方位属性椭圆拟合预测裂缝方法比叠后方法更加精确。裂缝型储层是一种典型的各向异性介质，是具有水平对称轴的横向各向同性介质。国内外纵多学者的研究表明，地下介质含有一组垂直排列且相互平行的裂缝时，地震波速度、振幅和层间旅行时等地震属性随地震测线方位角的不同而变化。当地震纵波在含有裂缝的各向异性介质中传播时，沿着裂缝方向传播速度快，垂直裂缝方向传播速度慢；沿着裂缝方向能量衰减小，垂直裂缝方向能量衰减大；沿着裂缝方向地震波射线层间旅行时间短，垂直裂缝方向地震波射线层间旅行时间长。地震波在裂缝介质中传播时显示的方位差异特征，可用于探测裂缝的发育特征，裂缝的发育特征包括裂缝发育强度和裂缝走向。
[0005]裂缝介质中地震属性在直角坐标系中是二倍方位角的余弦函数:
[0006]F(<j5)=A+Bcos2(<j5-<j5syni)(I)
[0007]其中，F是采集方位角为Φ时的地震属性，A, B是与裂缝介质参数、观测几何有关的变量，为裂缝走向。由于地震属性余弦函数的梯度B远小于均值Α，地震属性在极坐标系中近似为椭圆，椭圆半长轴等于|Α| + |Β|，半短轴等于|Α|-|Β|，前人应用椭圆的长短轴之比来表示裂缝强度，长轴方向表示裂缝走向。应用以上算法进行裂缝反演的不足之处在于:该裂缝反演方法进行椭圆拟合时至少需要五个方位的数据，拟合的椭圆中心点有可能不在极坐标系的坐标原点，使得受噪音干扰大。
【发明内容】
[0008]本发明实施例提供了一种地震纵波多方位属性椭圆拟合预测裂缝的方法及装置，解决了现有技术中裂缝反演方法受噪音干扰大的技术问题。
[0009]本发明实施例提供了一种地震纵波多方位属性椭圆拟合预测裂缝的方法，该方法包括:将至少包括三个方位的地震属性数据按照测线方位进行中心对称处理，生成具有完备方位的地震属性数据；对生成的具有完备方位的地震属性数据进行椭圆拟合；将椭圆扁率确定为裂缝发育强度。
[0010]在一个实施例中，将多方位的地震属性数据按照测线方位进行中心对称处理，包括:将各个测线方位的地震属性数据投影到与该测线方位关于原点中心对称的方位上。
[0011]在一个实施例中，对生成的具有完备方位的地震属性数据进行椭圆拟合，包括:将生成的具有完备方位的地震属性数据代入初始椭圆方程式:a(x-xc)2+b(x-xc)(y-yc) +c (y-yc)2=l ;计算获得拟合椭圆的方程式:ax2+bxy+cy2=l,其中，a、b、c均是由椭圆的半长轴长度、椭圆的半短轴长度以及椭圆的长轴所在的方位值组成的函数，xc, Yc分别为椭圆中心点的横纵坐标值，X，y是椭圆上点的横纵坐标值。
[0012]在一个实施例中，所述裂缝发育强度的值域范围是[0，1)。
[0013]在一个实施例中，在将所述椭圆的扁率确定为裂缝发育强度之后，还包括:计算每个面元的每个时间采样点上的裂缝发育强度，输出裂缝发育强度数据体。
[0014]在一个实施例中，在对生成的具有完备方位的地震属性数据进行椭圆拟合之后，还包括:根据椭圆拟合结果计算椭圆的长轴和短轴所在的方位值；将椭圆的长轴或短轴所指示的方位确定为裂缝走向。
[0015]在一个实施例中，将椭圆的长轴或短轴所指示的方位确定为裂缝走向，包括:在所述地震属性与速度正相关的情况下，将所述椭圆的长轴指示的方位确定为裂缝走向；在所述地震属性与时间正相关的情况下，将所述椭圆的短轴指示的方位确定为裂缝走向。
[0016]在一个实施例中，将椭圆的长轴或短轴所指示的方位确定为裂缝走向之后，还包括:计算每个面元的每个时间采样点上的裂缝走向，输出裂缝走向数据体。
[0017]本发明实施例还提供了一种地震纵波多方位属性椭圆拟合预测裂缝的装置，该装置包括:数据修正模块，用于将至少包括三个方位的地震属性数据按照测线方位进行中心对称处理，生成具有完备方位的地震属性数据；椭圆拟合模块，用于对生成的具有完备方位的地震属性数据进行椭圆拟合；裂缝预测模块，用于将椭圆扁率确定为裂缝发育强度。
[0018]在一个实施例中，所述数据修正模块，用于将各个测线方位的地震属性数据投影到与该测线方位关于原点中心对称的方位上。
[0019]在一个实施例中，所述椭圆拟合模块包括:椭圆拟合单元，用于将生成的具有完备方位的地震属性数据代入初始椭圆方程式:a(x-xe)2+b(x-x。)(y-yc) +c (y-yc)2=1 ;计算单元，用于计算获得拟合椭圆的方程式:aX2+bXy+cy2=l，其中，a、b、c均是由椭圆的半长轴长度、椭圆的半短轴长度以及椭圆的长轴所在的方位值组成的函数，xc, Yc分别为椭圆中心点的横纵坐标值，X，y是椭圆上点的横纵坐标值。
[0020]在一个实施例中，所述裂缝发育强度的值域范围是[0，1)。
[0021]在一个实施例中，还包括:第一计算模块，用于将所述椭圆的扁率确定为裂缝发育强度之后，计算每个面元的每个时间采样点上的裂缝发育强度，输出裂缝发育强度数据体。
[0022]在一个实施例中，所述裂缝预测模块，包括:计算单元，用于根据椭圆拟合结果计算椭圆的长轴和短轴所在的方位值；预测单元，用于将椭圆的长轴或短轴所指示的方位确定为裂缝走向。
[0023]在一个实施例中，所述预测单元还用于在所述地震属性与速度正相关的情况下，将所述椭圆的长轴指示的方位确定为裂缝走向；在所述地震属性与时间正相关的情况下，将所述椭圆的短轴指示的方位确定为裂缝走向。
[0024]在一个实施例中，还包括:第二计算模块，用于在将椭圆的长轴或短轴所指示的方位确定为裂缝走向之后，计算每个面元的每个时间采样点上的裂缝走向，输出裂缝走向数据体。
[0025]在本发明实施例中，将至少包括三个方位的地震属性数据按照测线方位进行中心对称处理，生成具有完备方位的地震属性数据，再对生成的具有完备方位的地震属性数据进行椭圆拟合，以确保椭圆的中心点位于极坐标的原点，将地震纵波多方位属性椭圆拟合预测裂缝所需要反演的参数从五个降低到了三个，从而降低了噪音对裂缝特征反演的干扰；同时，将椭圆扁率确定为裂缝发育强度，比现有技术中用长短轴之比描述裂缝发育强度更加合理。
【专利附图】
【附图说明】
[0026]此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中:
[0027]图1是本发明实施例提供的一种地震纵波多方位属性椭圆拟合预测裂缝的方法的流程图；
[0028]图2是本发明实施例提供的一种地震纵波多方位属性椭圆拟合示意图；
[0029]图3是本发明实施例提供的一种地震纵波多方位属性椭圆拟合预测裂缝的装置的结构框图。
【具体实施方式】
[0030]为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施方式和附图，对本发明做进一步详细说明。在此，本发明的示意性实施方式及其说明用于解释本发明，但并不作为对本发明的限定。
[0031]在本发明实施例中，提供了一种地震纵波多方位属性椭圆拟合预测裂缝的方法，如图1所示，该方法包括:
[0032]步骤101:将至少包括三个方位的地震属性数据按照测线方位进行中心对称处理，生成具有完备方位的地震属性数据；
[0033]步骤102:对生成的具有完备方位的地震属性数据进行椭圆拟合；
[0034]步骤103:将椭圆扁率确定为裂缝发育强度。
[0035]由图1所示的流程可知，在本发明实施例中，将至少包括三个方位的地震属性数据按照测线方位进行中心对称处理，生成具有完备方位的地震属性数据，再对生成的具有完备方位的地震属性数据进行椭圆拟合，以确保椭圆的中心点位于极坐标的原点，将地震纵波多方位属性椭圆拟合预测裂缝所需要反演的参数从五个降低到了三个，从而降低了噪音对裂缝特征反演的干扰；同时，将椭圆扁率确定为裂缝发育强度，比现有技术中用长短轴之比描述裂缝发育强度更加合理。
[0036]在具体实施过程中，将至少包括三个方位的地震属性数据按照测线方位进行中心对称处理的过程可以这样实现，将各个测线方位的地震属性数据投影到与该测线方位关于原点中心对称的方位上，如图2所示，该图为极坐标系，纵轴箭头方向为正北方向，顺时针为方位角增大方向。图中黑色圆点是地震属性数据在极坐标的投影，黑色圆点和坐标原点的连线与正北方向夹角表示方位角，黑色圆点与坐标原点的距离表示地震属性值的大小，将各个测线方位的地震属性数据投影到与该测线方位关于原点中心对称的方位上，例如，将方位角为90度的测线方位上的震属性数据投影到与其关于原点中心对称的270度的方位上。多条测线的地震属性值在极坐标下表现为椭圆曲线，由此可拟合椭圆长轴、短轴及长轴所在方位。在将地震属性数据按照测线方位进行扇区分类后，扇形区域如图2中阴影部分所示，每个扇形区域包括一定角度范围的测线方位角(如图2中所示的方位)，例如，测线方位角为O到30度的分为一个扇形区域，测线方位角为30到60度的分为一个扇形区域，然后可以计算获取每个扇形区域内每个测线方位角对应的地震属性数据，每个扇形区域内的所有地震属性值可以称为方位角道集，将各个测线方位的地震属性数据投影到与该测线方位关于原点中心对称的方位上，可以得到具有完备方位的地震属性数据体，由于各向异性介质的地震属性随测线方位角在极坐标下呈椭圆形态，对获得的具有完备方位的地震属性数据进行椭圆拟合，以确保椭圆的中心点位于极坐标的原点上。
[0037]在具体实施过程中，可以这样来获得拟合椭圆的方程式，获取初始椭圆方程式，初始椭圆的方程含有五个独立参数:
[0038]a(x-xc)2+b(x~xc) (y-yc) +c(y-yc)2=1(2)
[0039]其中，a、b、c均是由椭圆的半长轴长度、椭圆的半短轴长度以及椭圆的长轴所在的方位值组成的函数，Xc;，y。分别为椭圆中心点的横纵坐标值；本发明实施例结合地震属性方位中心对称性原理，将生成的具有完备方位的地震属性数据代入初始椭圆方程式，进行椭圆拟合，计算得到拟合椭圆的方 程式:
[0040]ax2+bxy+cy2=l(3)
[0041]其中，a、b、c均是由椭圆的半长轴长度、椭圆的半短轴长度以及椭圆的长轴所在的方位值组成的函数，X，Y是椭圆上点的横纵坐标值，拟合椭圆的中心点位于极坐标的原点上，即方程式(2)中的Xc=Ojc=Ot5需要反演的参数由方程(2)中的五个降为方程(3)中的三个，本发明实施例很好地压制噪音的不利影响，有助于准确地反演出地下储层的裂缝分布情况。
[0042]在对生成的具有完备方位的地震属性数据进行椭圆拟合后，为了方便统计区域性裂缝发育方向，使得裂缝发育强度更加合理，反演结果数量级一致，在本实施例中，根据椭圆拟合结果计算椭圆的扁率；将所述椭圆的扁率确定为裂缝发育强度，椭圆的扁率即为椭圆的长轴长度减去椭圆的短轴长度的差与椭圆的长轴长度的比值，即椭圆的扁率为裂缝发育强度。例如，该裂缝发育强度的值域范围可以是[0，I)。即当各个测线方位对应的地震属性值没有差异时(即为各向同性)，对应裂缝发育强度为零；此种裂缝发育强度定义方式比前人给出的利用长轴与短轴之比表示裂缝发育强度更加合理，反演结果数量级一致，在数量级上稳定，有利于制作描 述裂缝发育特征的玫瑰图，方便统计区域性裂缝发育方向。
[0043]在实施过程中，为了输出裂缝发育强度数据体，在本实施例中，在将所述椭圆的扁率确定为裂缝发育强度之后，计算每个面元的每个时间采样点上的裂缝发育强度，输出裂缝发育强度数据体，该面元可以是⑶P (Common Depth Point,共深度采集点)面元。
[0044]在对生成的具有完备方位的地震属性数据进行椭圆拟合后，为了准确地反演裂缝走向，在本发明实施中，根据椭圆拟合结果计算椭圆的长轴和短轴所在的方位值；充分考虑地震属性数据的物理意义，将椭圆的长轴或短轴所指示的方位确定为裂缝走向。具体实施时，在地震属性与速度正相关的情况下，将椭圆的长轴指示的方位确定为裂缝走向；在所述地震属性与时间正相关的情况下，将椭圆的短轴指示的方位确定为裂缝走向。即沿着裂缝平面的动校正速度最大，层间旅行时最小，垂直裂缝平面的动校正速度最小，层间旅行时最大。若地震属性与速度正相关，则椭圆长轴指示裂缝走向；若地震属性与时间正相关，则椭圆短轴指示裂缝走向。
[0045]在预测出裂缝发育强度和裂缝走向后，可以结合裂缝发育强度和裂缝走向，制作玫瑰图，以分析区域性的裂缝发育情况。[0046]在实施过程中，为了获得裂缝走向数据体，在本实施例中，在将椭圆的长轴或短轴所指示的方位确定为裂缝走向之后，计算每个面元的每个时间采样点上的裂缝走向，输出裂缝走向数据体，该面元可以是⑶P (Common Depth Point,共深度采集点)面元。
[0047]在具体实施时，地震属性可以包括以下之一:反射振幅、动校正速度、层间旅行时、振幅随偏移距的变化梯度、品质因子。
[0048]各向异性裂缝介质模型中含有一系列垂直排列且相互平行的裂缝。地震射线的入射角为地震波射线与反射面法向的夹角；方位角定义为地震波的激发点和接收点的连线与X轴的夹角。当地震波沿着不同的射线传播时，振幅不仅随入射角变化，还随着方位角而变化。不仅振幅存在方位差异，传播速度和AV0(Amplitude Versus Offset,振幅随偏移距的变化)梯度等地震属性也随方位变化。
[0049]以下结合具体实施例来描述上述地震纵波多方位属性椭圆拟合预测裂缝的方法，该方法包括如下步骤:
[0050]步骤1:采用三维地震资料采集技术，获得至少包括三个方位的叠前地震资料(地震属性数据)，也就是获得宽方位(也称多方位)的地震属性数据，其中，纵横比大于0.5为宽方位。
[0051]步骤2:将地震属性数据按照测线方位进行扇区分类，进行中心对称处理，生成具有完备方位的地震属性数据体。
[0052]步骤3:对生成的具有完备方位的地震属性数据进行椭圆拟合，椭圆拟合算法为:椭圆的一般方程含有五个独立参数:
[0053]a(x-xc)2+b(x~xc) (y-yc) +c(y-yc)2=1(2)
[0054]其中,a、b、c均是由椭圆的半长轴长度、椭圆的半短轴长度以及椭圆的长轴所在的方位值组成的函数，xc, yc分别为椭圆中心点的横纵坐标值，X，y是椭圆上点的横纵坐标值；本发明实施例的改进之处在于考虑了地震属性的方位对称性，使得地震属性椭圆的中心点位于极坐标的原点，椭圆方程由一般方程(2)简化为:
[0055]ax2+bxy+cy2=l(3)
[0056]需要反演的参数由方程(2)中的五个降为方程(3)中的三个，然后计算椭圆半长轴长度、半短轴长度、长轴所在方位值，进而反演裂缝强度和裂缝走向，能有效压制噪音的不利影响。
[0057]步骤4:根据椭圆计算最大和最小方位属性值，即平行和垂直裂缝方位的方位属性值。
[0058]步骤5:定义椭圆扁率(方位属性的相对变化率)为裂缝发育强度，即将椭圆的长轴长度减去椭圆的短轴长度的差与所述椭圆的长轴长度的比值作为所述裂缝发育强度，其值域范围为[0，I);当各个方位的地震属性值没有差异时(即为各向同性)，对应的裂缝发育强度为零；此种裂缝发育强度定义方式比前人给出的利用长轴与短轴之比表示裂缝发育强度更加合理，反演结果数量级一致，方便统计区域性裂缝发育方向；算出每个面元的每个时间采样点上的裂缝强度，输出裂缝发育强度数据体。
[0059]步骤6:根据地震属性的物理意义，定义椭圆长轴或者短轴方向为裂缝走向，由于地震属性的物理意义不一致，沿着裂缝平面的动校正速度最大，层间旅行时最小，垂直裂缝平面的动校正速度最小，层间旅行时最大。若地震属性与速度正相关，则椭圆长轴指示裂缝走向；若地震属性与时间正相关，则椭圆短轴指示裂缝走向；算出每个面元的每个时间采样点上的裂缝走向，输出裂缝走向数据体。
[0060]步骤7:提取油气勘探目的层的反演结果，反演裂缝发育区域以及裂缝走向。
[0061]基于同一发明构思，本发明实施例中还提供了一种地震纵波多方位属性椭圆拟合预测裂缝的装置，如下面的实施例所述。由于地震纵波多方位属性椭圆拟合预测裂缝的装置解决问题的原理与地震纵波多方位属性椭圆拟合预测裂缝的方法相似，因此地震纵波多方位属性椭圆拟合预测裂缝的装置的实施可以参见地震纵波多方位属性椭圆拟合预测裂缝的方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
[0062]图3是本发明实施例的地震纵波多方位属性椭圆拟合预测裂缝的装置的一种结构框图，如图3所示，包括:数据修正模块301，椭圆拟合模块302，以及裂缝预测模块303，下面对该结构进行说明。
[0063]数据修正模块301，用于将至少包括三个方位的地震属性数据按照测线方位进行中心对称处理，生成具有完备方位的地震属性数据；
[0064]椭圆拟合模块302，与数据修正模块301连接，用于对生成的具有完备方位的地震属性数据进行椭圆拟合；
[0065]裂缝预测模块303，与椭圆拟合模块302连接，用于将椭圆扁率确定为裂缝发育强度。
[0066]在一个实施例中，所述数据修正模块，用于将各个测线方位的地震属性数据投影到与该测线方位关于原点中心对称的方位上
[0067]在一个实施例中，所述椭圆拟合模块302包括:椭圆拟合单元，用于将生成的具有完备方位的地震属性数据代入初始椭圆方程式:a(x-xe)2+b(x-x。)(y-yc) +c (y-yc)2=l ;计算单元，与椭圆拟合单元连接，用于计算获得拟合椭圆的方程式:ax2+bxy+cy2=l,其中，a、b、c均是由椭圆的半长轴长度、椭圆的半短轴长度以及椭圆的长轴所在的方位值组成的函数，xc, yc分别为椭圆中心点的横纵坐标值，X，y是椭圆上点的横纵坐标值。
[0068]在一个实施例中，所述裂缝发育强度的值域范围是[0，I)。[0069]在一个实施例中，还包括:第一计算模块，用于将所述椭圆的扁率确定为裂缝发育强度之后，计算每个面元的每个时间采样点上的裂缝发育强度，输出裂缝发育强度数据体。
[0070]在一个实施例中，所述裂缝预测模块303还包括:计算单元，用于根据椭圆拟合结果计算椭圆的长轴和短轴所在的方位值；预测单元，与计算单元连接，用于将椭圆的长轴或短轴所指示的方位确定为裂缝走向。
[0071]在一个实施例中，所述预测单元还用于在所述地震属性与速度正相关的情况下，将所述椭圆的长轴指示的方位确定为裂缝走向；在所述地震属性与时间正相关的情况下，将所述椭圆的短轴指示的方位确定为裂缝走向。
[0072]在一个实施例中，还包括:第二计算模块，用于在将椭圆的长轴或短轴所指示的方位确定为裂缝走向之后，计算每个面元的每个时间采样点上的裂缝走向，输出裂缝走向数据体。
[0073]在一个实施例中，所述地震属性可以包括以下之一:地震波反射振幅、动校正速度、层间旅行时、振幅随偏移距的变化梯度、品质因子。
[0074]在本发明实施例中，将多方位的地震属性数据按照测线方位进行中心对称处理，生成具有完备方位的地震属性数据，再对生成的具有完备方位的地震属性数据进行椭圆拟合，以确保椭圆的中心点位于极坐标的原点，将纵波预测裂缝所需要反演的参数从五个降低到了三个，从而降低了噪音对裂缝特征反演的干扰；同时，将椭圆扁率确定为裂缝发育强度，比现有技术中用长短轴之比描述裂缝发育强度更加合理。
[0075]显然，本领域的技术人员应该明白，上述的本发明实施例的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，并且在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本发明实施例不限制于任何特定的硬件和软件结合。
[0076]以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
【权利要求】
1.一种地震纵波多方位属性椭圆拟合预测裂缝的方法，其特征在于，包括: 将至少包括三个方位的地震属性数据按照测线方位进行中心对称处理，生成具有完备方位的地震属性数据； 对生成的具有完备方位的地震属性数据进行椭圆拟合； 将椭圆扁率确定为裂缝发育强度。
2.如权利要求1所述地震纵波多方位属性椭圆拟合预测裂缝的方法，其特征在于，将至少包括三个方位的地震属性数据按照测线方位进行中心对称处理，包括: 将各个测线方位的地震属性数据投影到与该测线方位关于原点中心对称的方位上。
3.如权利要求1所述地震纵波多方位属性椭圆拟合预测裂缝的方法，其特征在于，对生成的具有完备方位的地震属性数据进行椭圆拟合，包括: 将生成的具有完备方位的地震属性数据代入初始椭圆方程式:3(11。)2+13(11。)(y-yc) +c (y-yc)2=l ； 计算获得拟合椭圆的方程式:ax2+bxy+cy2=l, 其中，a、b、c均是由椭圆的半长轴长度、椭圆的半短轴长度以及椭圆的长轴所在的方位值组成的函数，xc, yc分别为椭圆中心点的横纵坐标值，X，y是椭圆上点的横纵坐标值。
4.如权利要求1所述地震纵波多方位属性椭圆拟合预测裂缝的方法，其特征在于，所述裂缝发育强度的值域范围是[O，I)。
5.如权利要求1所述地震纵波多方位属性椭圆拟合预测裂缝的方法，其特征在于，在将所述椭圆的扁率确定为裂缝发育强度`之后，还包括: 计算每个面元的每个时间采样点上的裂缝发育强度，输出裂缝发育强度数据体。
6.如权利要求1至5中任一项所述地震纵波多方位属性椭圆拟合预测裂缝的方法，其特征在于，在对生成的具有完备方位的地震属性数据进行椭圆拟合之后，还包括: 根据椭圆拟合结果计算椭圆的长轴和短轴所在的方位值； 将椭圆的长轴或短轴所指示的方位确定为裂缝走向。
7.如权利要求6所述地震纵波多方位属性椭圆拟合预测裂缝的方法，其特征在于，将椭圆的长轴或短轴所指示的方位确定为裂缝走向，包括: 在所述地震属性与速度正相关的情况下，将所述椭圆的长轴指示的方位确定为裂缝走向； 在所述地震属性与时间正相关的情况下，将所述椭圆的短轴指示的方位确定为裂缝走向。
8.如权利要求6所述地震纵波多方位属性椭圆拟合预测裂缝的方法，其特征在于，将椭圆的长轴或短轴所指示的方位确定为裂缝走向之后，还包括: 计算每个面元的每个时间采样点上的裂缝走向，输出裂缝走向数据体。
9.一种地震纵波多方位属性椭圆拟合预测裂缝的装置，其特征在于，包括: 数据修正模块，用于将至少包括三个方位的地震属性数据按照测线方位进行中心对称处理，生成具有完备方位的地震属性数据； 椭圆拟合模块，用于对生成的具有完备方位的地震属性数据进行椭圆拟合； 裂缝预测模块，用于将椭圆扁率确定为裂缝发育强度。
10.如权利要求9所述地震纵波多方位属性椭圆拟合预测裂缝的装置，其特征在于，所述数据修正模块，用于将各个测线方位的地震属性数据投影到与该测线方位关于原点中心对称的方位上。
11.如权利要求9所述地震纵波多方位属性椭圆拟合预测裂缝的装置，其特征在于，所述椭圆拟合模块包括: 椭圆拟合单元，用于将生成的具有完备方位的地震属性数据代入初始椭圆方程式:a (x-xc) 2+b (x-xc) (y-yc) +c (y-yc) 2=1 ； 计算单元，用于计算获得拟合椭圆的方程式:ax2+bxy+cy2=l, 其中，a、b、c均是由椭圆的半长轴长度、椭圆的半短轴长度以及椭圆的长轴所在的方位值组成的函数，xc, yc分别为椭圆中心点的横纵坐标值，X，y是椭圆上点的横纵坐标值。
12.如权利要求9所述地震纵波多方位属性椭圆拟合预测裂缝的装置，其特征在于，所述裂缝发育强度的值域范围是[O，I)。
13.如权利要求9所述地震纵波多方位属性椭圆拟合预测裂缝的装置，其特征在于，还包括: 第一计算模块，用于将所述椭圆的扁率确定为裂缝发育强度之后，计算每个面元的每个时间采样点上的裂缝发育强度，输出裂缝发育强度数据体。
14.如权利要求9至13中任一项所述地震纵波多方位属性椭圆拟合预测裂缝的装置，其特征在于，所述裂缝预测模块，包括: 计算单元，用于根据椭圆拟合结果计算椭圆的长轴和短轴所在的方位值； 预测单元，用于将椭圆的长 轴或短轴所指示的方位确定为裂缝走向。
15.如权利要求14所述地震纵波多方位属性椭圆拟合预测裂缝的装置，其特征在于，所述预测单元还用于在所述地震属性与速度正相关的情况下，将所述椭圆的长轴指示的方位确定为裂缝走向；在所述地震属性与时间正相关的情况下，将所述椭圆的短轴指示的方位确定为裂缝走向。
16.如权利要求14所述地震纵波多方位属性椭圆拟合预测裂缝的装置，其特征在于，还包括: 第二计算模块，用于在将椭圆的长轴或短轴所指示的方位确定为裂缝走向之后，计算每个面元的每个时间采样点上 的裂缝走向，输出裂缝走向数据体。
【文档编号】G01V1/36GK103869359SQ201410064797
【公开日】2014年6月18日 申请日期:2014年2月25日 优先权日:2014年2月25日
【发明者】谢春辉, 杨午阳, 周春雷, 张闻亭, 鄢高韩, 王洪求, 王恩利, 杨庆, 李琳, 魏新建 申请人:中国石油天然气股份有限公司