CN103869359A - 地震纵波多方位属性椭圆拟合预测裂缝的方法及装置 - Google Patents
地震纵波多方位属性椭圆拟合预测裂缝的方法及装置 Download PDFInfo
- Publication number
- CN103869359A CN103869359A CN201410064797.1A CN201410064797A CN103869359A CN 103869359 A CN103869359 A CN 103869359A CN 201410064797 A CN201410064797 A CN 201410064797A CN 103869359 A CN103869359 A CN 103869359A
- Authority
- CN
- China
- Prior art keywords
- orientation
- ellipse
- fracture
- crack
- ellipse fitting
- 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.)
- Pending
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明实施例提供了一种地震纵波多方位属性椭圆拟合预测裂缝的方法及装置,其中,该方法包括:将至少包括三个方位的地震属性数据按照测线方位进行中心对称处理,生成具有完备方位的地震属性数据;对生成的具有完备方位的地震属性数据进行椭圆拟合;将椭圆扁率确定为裂缝发育强度。本发明实施例确保了椭圆的中心点位于极坐标的原点,将纵波预测裂缝所需要反演的参数从五个降低到了三个,从而降低了噪音对裂缝特征反演的干扰;同时,将椭圆扁率确定为裂缝发育强度,比现有技术中用长短轴之比描述裂缝发育强度更加合理。
Description
技术领域
本发明涉及石油勘探技术领域,特别涉及一种地震纵波多方位属性椭圆拟合预测裂缝的方法及装置。
背景技术
裂缝是岩石受力发生破裂、沿破裂面两侧的岩石没有发生明显位移的一种断裂构造现象。裂缝既是油、气、水等地下流体的运移通道,也是流体的储集空间,裂缝能极大地提升储层的渗透率。有效裂缝是指开启并充填流体的裂缝,有效裂缝的探测对于识别油气运移通道,寻找有效储集空间具有重要的指导意义。
目前比较成熟的裂缝预测技术是利用相干和曲率等叠后不连续属性、岩心和成像测井来预测裂缝,这些方法存在较多不足。相干反演能够反映小断裂或断裂发育带,岩性突变带,但主要是用于断层边缘的检测,并且检测所有的断裂体系。当地下介质受多期次构造运动或成岩作用时,早期的裂缝可能被泥质或者方解石等所充填而成为无效裂缝,而这种无效裂缝在相干切片上仍能显示,相干反演技术很难预测出有效裂缝的发育特征。采集岩心和成像测井成本高昂,导致资料数量少,只能在井点特定井段处检测裂缝,无法在横向大面积和纵向全深度预测裂缝。而地震纵波宽方位属性能够检测有效裂缝引起的各向异性特征,可大规模检测裂缝,其中,纵横比大于0.5为宽方位。
各向异性介质的物理性质(如地震波的速度,振幅等)随测线方位改变而变化,因此,利用方位角信息,采用地震纵波宽方位属性椭圆拟合预测裂缝方法比叠后方法更加精确。裂缝型储层是一种典型的各向异性介质,是具有水平对称轴的横向各向同性介质。国内外纵多学者的研究表明,地下介质含有一组垂直排列且相互平行的裂缝时,地震波速度、振幅和层间旅行时等地震属性随地震测线方位角的不同而变化。当地震纵波在含有裂缝的各向异性介质中传播时,沿着裂缝方向传播速度快,垂直裂缝方向传播速度慢;沿着裂缝方向能量衰减小,垂直裂缝方向能量衰减大;沿着裂缝方向地震波射线层间旅行时间短,垂直裂缝方向地震波射线层间旅行时间长。地震波在裂缝介质中传播时显示的方位差异特征,可用于探测裂缝的发育特征,裂缝的发育特征包括裂缝发育强度和裂缝走向。
裂缝介质中地震属性在直角坐标系中是二倍方位角的余弦函数:
F(φ)=A+Bcos2(φ-φsym) (1)
其中,F是采集方位角为φ时的地震属性,A,B是与裂缝介质参数、观测几何有关的变量,φsym为裂缝走向。由于地震属性余弦函数的梯度B远小于均值A,地震属性在极坐标系中近似为椭圆,椭圆半长轴等于|A|+|B|,半短轴等于|A|-|B|,前人应用椭圆的长短轴之比来表示裂缝强度,长轴方向表示裂缝走向。应用以上算法进行裂缝反演的不足之处在于:该裂缝反演方法进行椭圆拟合时至少需要五个方位的数据,拟合的椭圆中心点有可能不在极坐标系的坐标原点,使得受噪音干扰大。
发明内容
本发明实施例提供了一种地震纵波多方位属性椭圆拟合预测裂缝的方法及装置,解决了现有技术中裂缝反演方法受噪音干扰大的技术问题。
本发明实施例提供了一种地震纵波多方位属性椭圆拟合预测裂缝的方法,该方法包括:将至少包括三个方位的地震属性数据按照测线方位进行中心对称处理,生成具有完备方位的地震属性数据;对生成的具有完备方位的地震属性数据进行椭圆拟合;将椭圆扁率确定为裂缝发育强度。
在一个实施例中,将多方位的地震属性数据按照测线方位进行中心对称处理,包括:将各个测线方位的地震属性数据投影到与该测线方位关于原点中心对称的方位上。
在一个实施例中,对生成的具有完备方位的地震属性数据进行椭圆拟合,包括:将生成的具有完备方位的地震属性数据代入初始椭圆方程式:a(x-xc)2+b(x-xc)(y-yc)+c(y-yc)2=1;计算获得拟合椭圆的方程式:ax2+bxy+cy2=1,其中,a、b、c均是由椭圆的半长轴长度、椭圆的半短轴长度以及椭圆的长轴所在的方位值组成的函数,xc,yc分别为椭圆中心点的横纵坐标值,x,y是椭圆上点的横纵坐标值。
在一个实施例中,所述裂缝发育强度的值域范围是[0,1)。
在一个实施例中,在将所述椭圆的扁率确定为裂缝发育强度之后,还包括:计算每个面元的每个时间采样点上的裂缝发育强度,输出裂缝发育强度数据体。
在一个实施例中,在对生成的具有完备方位的地震属性数据进行椭圆拟合之后,还包括:根据椭圆拟合结果计算椭圆的长轴和短轴所在的方位值;将椭圆的长轴或短轴所指示的方位确定为裂缝走向。
在一个实施例中,将椭圆的长轴或短轴所指示的方位确定为裂缝走向,包括:在所述地震属性与速度正相关的情况下,将所述椭圆的长轴指示的方位确定为裂缝走向;在所述地震属性与时间正相关的情况下,将所述椭圆的短轴指示的方位确定为裂缝走向。
在一个实施例中,将椭圆的长轴或短轴所指示的方位确定为裂缝走向之后,还包括:计算每个面元的每个时间采样点上的裂缝走向,输出裂缝走向数据体。
本发明实施例还提供了一种地震纵波多方位属性椭圆拟合预测裂缝的装置,该装置包括:数据修正模块,用于将至少包括三个方位的地震属性数据按照测线方位进行中心对称处理,生成具有完备方位的地震属性数据;椭圆拟合模块,用于对生成的具有完备方位的地震属性数据进行椭圆拟合;裂缝预测模块,用于将椭圆扁率确定为裂缝发育强度。
在一个实施例中,所述数据修正模块,用于将各个测线方位的地震属性数据投影到与该测线方位关于原点中心对称的方位上。
在一个实施例中,所述椭圆拟合模块包括:椭圆拟合单元,用于将生成的具有完备方位的地震属性数据代入初始椭圆方程式:a(x-xc)2+b(x-xc)(y-yc)+c(y-yc)2=1;计算单元,用于计算获得拟合椭圆的方程式:ax2+bxy+cy2=1,其中,a、b、c均是由椭圆的半长轴长度、椭圆的半短轴长度以及椭圆的长轴所在的方位值组成的函数,xc,yc分别为椭圆中心点的横纵坐标值,x,y是椭圆上点的横纵坐标值。
在一个实施例中,所述裂缝发育强度的值域范围是[0,1)。
在一个实施例中,还包括:第一计算模块,用于将所述椭圆的扁率确定为裂缝发育强度之后,计算每个面元的每个时间采样点上的裂缝发育强度,输出裂缝发育强度数据体。
在一个实施例中,所述裂缝预测模块,包括:计算单元,用于根据椭圆拟合结果计算椭圆的长轴和短轴所在的方位值;预测单元,用于将椭圆的长轴或短轴所指示的方位确定为裂缝走向。
在一个实施例中,所述预测单元还用于在所述地震属性与速度正相关的情况下,将所述椭圆的长轴指示的方位确定为裂缝走向;在所述地震属性与时间正相关的情况下,将所述椭圆的短轴指示的方位确定为裂缝走向。
在一个实施例中,还包括:第二计算模块,用于在将椭圆的长轴或短轴所指示的方位确定为裂缝走向之后,计算每个面元的每个时间采样点上的裂缝走向,输出裂缝走向数据体。
在本发明实施例中,将至少包括三个方位的地震属性数据按照测线方位进行中心对称处理,生成具有完备方位的地震属性数据,再对生成的具有完备方位的地震属性数据进行椭圆拟合,以确保椭圆的中心点位于极坐标的原点,将地震纵波多方位属性椭圆拟合预测裂缝所需要反演的参数从五个降低到了三个,从而降低了噪音对裂缝特征反演的干扰;同时,将椭圆扁率确定为裂缝发育强度,比现有技术中用长短轴之比描述裂缝发育强度更加合理。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,并不构成对本发明的限定。在附图中:
图1是本发明实施例提供的一种地震纵波多方位属性椭圆拟合预测裂缝的方法的流程图;
图2是本发明实施例提供的一种地震纵波多方位属性椭圆拟合示意图;
图3是本发明实施例提供的一种地震纵波多方位属性椭圆拟合预测裂缝的装置的结构框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施方式和附图,对本发明做进一步详细说明。在此,本发明的示意性实施方式及其说明用于解释本发明,但并不作为对本发明的限定。
在本发明实施例中,提供了一种地震纵波多方位属性椭圆拟合预测裂缝的方法,如图1所示,该方法包括:
步骤101:将至少包括三个方位的地震属性数据按照测线方位进行中心对称处理,生成具有完备方位的地震属性数据;
步骤102:对生成的具有完备方位的地震属性数据进行椭圆拟合;
步骤103:将椭圆扁率确定为裂缝发育强度。
由图1所示的流程可知,在本发明实施例中,将至少包括三个方位的地震属性数据按照测线方位进行中心对称处理,生成具有完备方位的地震属性数据,再对生成的具有完备方位的地震属性数据进行椭圆拟合,以确保椭圆的中心点位于极坐标的原点,将地震纵波多方位属性椭圆拟合预测裂缝所需要反演的参数从五个降低到了三个,从而降低了噪音对裂缝特征反演的干扰;同时,将椭圆扁率确定为裂缝发育强度,比现有技术中用长短轴之比描述裂缝发育强度更加合理。
在具体实施过程中,将至少包括三个方位的地震属性数据按照测线方位进行中心对称处理的过程可以这样实现,将各个测线方位的地震属性数据投影到与该测线方位关于原点中心对称的方位上,如图2所示,该图为极坐标系,纵轴箭头方向为正北方向,顺时针为方位角增大方向。图中黑色圆点是地震属性数据在极坐标的投影,黑色圆点和坐标原点的连线与正北方向夹角表示方位角,黑色圆点与坐标原点的距离表示地震属性值的大小,将各个测线方位的地震属性数据投影到与该测线方位关于原点中心对称的方位上,例如,将方位角为90度的测线方位上的震属性数据投影到与其关于原点中心对称的270度的方位上。多条测线的地震属性值在极坐标下表现为椭圆曲线,由此可拟合椭圆长轴、短轴及长轴所在方位。在将地震属性数据按照测线方位进行扇区分类后,扇形区域如图2中阴影部分所示,每个扇形区域包括一定角度范围的测线方位角(如图2中所示的方位),例如,测线方位角为0到30度的分为一个扇形区域,测线方位角为30到60度的分为一个扇形区域,然后可以计算获取每个扇形区域内每个测线方位角对应的地震属性数据,每个扇形区域内的所有地震属性值可以称为方位角道集,将各个测线方位的地震属性数据投影到与该测线方位关于原点中心对称的方位上,可以得到具有完备方位的地震属性数据体,由于各向异性介质的地震属性随测线方位角在极坐标下呈椭圆形态,对获得的具有完备方位的地震属性数据进行椭圆拟合,以确保椭圆的中心点位于极坐标的原点上。
在具体实施过程中,可以这样来获得拟合椭圆的方程式,获取初始椭圆方程式,初始椭圆的方程含有五个独立参数:
a(x-xc)2+b(x-xc)(y-yc)+c(y-yc)2=1 (2)
其中,a、b、c均是由椭圆的半长轴长度、椭圆的半短轴长度以及椭圆的长轴所在的方位值组成的函数,xc,yc分别为椭圆中心点的横纵坐标值;本发明实施例结合地震属性方位中心对称性原理,将生成的具有完备方位的地震属性数据代入初始椭圆方程式,进行椭圆拟合,计算得到拟合椭圆的方程式:
ax2+bxy+cy2=1 (3)
其中,a、b、c均是由椭圆的半长轴长度、椭圆的半短轴长度以及椭圆的长轴所在的方位值组成的函数,x,y是椭圆上点的横纵坐标值,拟合椭圆的中心点位于极坐标的原点上,即方程式(2)中的xc=0,yc=0。需要反演的参数由方程(2)中的五个降为方程(3)中的三个,本发明实施例很好地压制噪音的不利影响,有助于准确地反演出地下储层的裂缝分布情况。
在对生成的具有完备方位的地震属性数据进行椭圆拟合后,为了方便统计区域性裂缝发育方向,使得裂缝发育强度更加合理,反演结果数量级一致,在本实施例中,根据椭圆拟合结果计算椭圆的扁率;将所述椭圆的扁率确定为裂缝发育强度,椭圆的扁率即为椭圆的长轴长度减去椭圆的短轴长度的差与椭圆的长轴长度的比值,即椭圆的扁率为裂缝发育强度。例如,该裂缝发育强度的值域范围可以是[0,1)。即当各个测线方位对应的地震属性值没有差异时(即为各向同性),对应裂缝发育强度为零;此种裂缝发育强度定义方式比前人给出的利用长轴与短轴之比表示裂缝发育强度更加合理,反演结果数量级一致,在数量级上稳定,有利于制作描述裂缝发育特征的玫瑰图,方便统计区域性裂缝发育方向。
在实施过程中,为了输出裂缝发育强度数据体,在本实施例中,在将所述椭圆的扁率确定为裂缝发育强度之后,计算每个面元的每个时间采样点上的裂缝发育强度,输出裂缝发育强度数据体,该面元可以是CDP(Common Depth Point,共深度采集点)面元。
在对生成的具有完备方位的地震属性数据进行椭圆拟合后,为了准确地反演裂缝走向,在本发明实施中,根据椭圆拟合结果计算椭圆的长轴和短轴所在的方位值;充分考虑地震属性数据的物理意义,将椭圆的长轴或短轴所指示的方位确定为裂缝走向。具体实施时,在地震属性与速度正相关的情况下,将椭圆的长轴指示的方位确定为裂缝走向;在所述地震属性与时间正相关的情况下,将椭圆的短轴指示的方位确定为裂缝走向。即沿着裂缝平面的动校正速度最大,层间旅行时最小,垂直裂缝平面的动校正速度最小,层间旅行时最大。若地震属性与速度正相关,则椭圆长轴指示裂缝走向;若地震属性与时间正相关,则椭圆短轴指示裂缝走向。
在预测出裂缝发育强度和裂缝走向后,可以结合裂缝发育强度和裂缝走向,制作玫瑰图,以分析区域性的裂缝发育情况。
在实施过程中,为了获得裂缝走向数据体,在本实施例中,在将椭圆的长轴或短轴所指示的方位确定为裂缝走向之后,计算每个面元的每个时间采样点上的裂缝走向,输出裂缝走向数据体,该面元可以是CDP(Common Depth Point,共深度采集点)面元。
在具体实施时,地震属性可以包括以下之一:反射振幅、动校正速度、层间旅行时、振幅随偏移距的变化梯度、品质因子。
各向异性裂缝介质模型中含有一系列垂直排列且相互平行的裂缝。地震射线的入射角为地震波射线与反射面法向的夹角;方位角定义为地震波的激发点和接收点的连线与x轴的夹角。当地震波沿着不同的射线传播时,振幅不仅随入射角变化,还随着方位角而变化。不仅振幅存在方位差异,传播速度和AVO(Amplitude Versus Offset,振幅随偏移距的变化)梯度等地震属性也随方位变化。
以下结合具体实施例来描述上述地震纵波多方位属性椭圆拟合预测裂缝的方法,该方法包括如下步骤:
步骤1:采用三维地震资料采集技术,获得至少包括三个方位的叠前地震资料(地震属性数据),也就是获得宽方位(也称多方位)的地震属性数据,其中,纵横比大于0.5为宽方位。
步骤2:将地震属性数据按照测线方位进行扇区分类,进行中心对称处理,生成具有完备方位的地震属性数据体。
步骤3:对生成的具有完备方位的地震属性数据进行椭圆拟合,椭圆拟合算法为:椭圆的一般方程含有五个独立参数:
a(x-xc)2+b(x-xc)(y-yc)+c(y-yc)2=1 (2)
其中,a、b、c均是由椭圆的半长轴长度、椭圆的半短轴长度以及椭圆的长轴所在的方位值组成的函数,xc,yc分别为椭圆中心点的横纵坐标值,x,y是椭圆上点的横纵坐标值;本发明实施例的改进之处在于考虑了地震属性的方位对称性,使得地震属性椭圆的中心点位于极坐标的原点,椭圆方程由一般方程(2)简化为:
ax2+bxy+cy2=1 (3)
需要反演的参数由方程(2)中的五个降为方程(3)中的三个,然后计算椭圆半长轴长度、半短轴长度、长轴所在方位值,进而反演裂缝强度和裂缝走向,能有效压制噪音的不利影响。
步骤4:根据椭圆计算最大和最小方位属性值,即平行和垂直裂缝方位的方位属性值。
步骤5:定义椭圆扁率(方位属性的相对变化率)为裂缝发育强度,即将椭圆的长轴长度减去椭圆的短轴长度的差与所述椭圆的长轴长度的比值作为所述裂缝发育强度,其值域范围为[0,1);当各个方位的地震属性值没有差异时(即为各向同性),对应的裂缝发育强度为零;此种裂缝发育强度定义方式比前人给出的利用长轴与短轴之比表示裂缝发育强度更加合理,反演结果数量级一致,方便统计区域性裂缝发育方向;算出每个面元的每个时间采样点上的裂缝强度,输出裂缝发育强度数据体。
步骤6:根据地震属性的物理意义,定义椭圆长轴或者短轴方向为裂缝走向,由于地震属性的物理意义不一致,沿着裂缝平面的动校正速度最大,层间旅行时最小,垂直裂缝平面的动校正速度最小,层间旅行时最大。若地震属性与速度正相关,则椭圆长轴指示裂缝走向;若地震属性与时间正相关,则椭圆短轴指示裂缝走向;算出每个面元的每个时间采样点上的裂缝走向,输出裂缝走向数据体。
步骤7:提取油气勘探目的层的反演结果,反演裂缝发育区域以及裂缝走向。
基于同一发明构思,本发明实施例中还提供了一种地震纵波多方位属性椭圆拟合预测裂缝的装置,如下面的实施例所述。由于地震纵波多方位属性椭圆拟合预测裂缝的装置解决问题的原理与地震纵波多方位属性椭圆拟合预测裂缝的方法相似,因此地震纵波多方位属性椭圆拟合预测裂缝的装置的实施可以参见地震纵波多方位属性椭圆拟合预测裂缝的方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
图3是本发明实施例的地震纵波多方位属性椭圆拟合预测裂缝的装置的一种结构框图,如图3所示,包括:数据修正模块301,椭圆拟合模块302,以及裂缝预测模块303,下面对该结构进行说明。
数据修正模块301,用于将至少包括三个方位的地震属性数据按照测线方位进行中心对称处理,生成具有完备方位的地震属性数据;
椭圆拟合模块302,与数据修正模块301连接,用于对生成的具有完备方位的地震属性数据进行椭圆拟合;
裂缝预测模块303,与椭圆拟合模块302连接,用于将椭圆扁率确定为裂缝发育强度。
在一个实施例中,所述数据修正模块,用于将各个测线方位的地震属性数据投影到与该测线方位关于原点中心对称的方位上
在一个实施例中,所述椭圆拟合模块302包括:椭圆拟合单元,用于将生成的具有完备方位的地震属性数据代入初始椭圆方程式:a(x-xc)2+b(x-xc)(y-yc)+c(y-yc)2=1;计算单元,与椭圆拟合单元连接,用于计算获得拟合椭圆的方程式:ax2+bxy+cy2=1,其中,a、b、c均是由椭圆的半长轴长度、椭圆的半短轴长度以及椭圆的长轴所在的方位值组成的函数,xc,yc分别为椭圆中心点的横纵坐标值,x,y是椭圆上点的横纵坐标值。
在一个实施例中,所述裂缝发育强度的值域范围是[0,1)。
在一个实施例中,还包括:第一计算模块,用于将所述椭圆的扁率确定为裂缝发育强度之后,计算每个面元的每个时间采样点上的裂缝发育强度,输出裂缝发育强度数据体。
在一个实施例中,所述裂缝预测模块303还包括:计算单元,用于根据椭圆拟合结果计算椭圆的长轴和短轴所在的方位值;预测单元,与计算单元连接,用于将椭圆的长轴或短轴所指示的方位确定为裂缝走向。
在一个实施例中,所述预测单元还用于在所述地震属性与速度正相关的情况下,将所述椭圆的长轴指示的方位确定为裂缝走向;在所述地震属性与时间正相关的情况下,将所述椭圆的短轴指示的方位确定为裂缝走向。
在一个实施例中,还包括:第二计算模块,用于在将椭圆的长轴或短轴所指示的方位确定为裂缝走向之后,计算每个面元的每个时间采样点上的裂缝走向,输出裂缝走向数据体。
在一个实施例中,所述地震属性可以包括以下之一:地震波反射振幅、动校正速度、层间旅行时、振幅随偏移距的变化梯度、品质因子。
在本发明实施例中,将多方位的地震属性数据按照测线方位进行中心对称处理,生成具有完备方位的地震属性数据,再对生成的具有完备方位的地震属性数据进行椭圆拟合,以确保椭圆的中心点位于极坐标的原点,将纵波预测裂缝所需要反演的参数从五个降低到了三个,从而降低了噪音对裂缝特征反演的干扰;同时,将椭圆扁率确定为裂缝发育强度,比现有技术中用长短轴之比描述裂缝发育强度更加合理。
显然,本领域的技术人员应该明白,上述的本发明实施例的各模块或各步骤可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,并且在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本发明实施例不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (16)
1.一种地震纵波多方位属性椭圆拟合预测裂缝的方法,其特征在于,包括:
将至少包括三个方位的地震属性数据按照测线方位进行中心对称处理,生成具有完备方位的地震属性数据;
对生成的具有完备方位的地震属性数据进行椭圆拟合;
将椭圆扁率确定为裂缝发育强度。
2.如权利要求1所述地震纵波多方位属性椭圆拟合预测裂缝的方法,其特征在于,将至少包括三个方位的地震属性数据按照测线方位进行中心对称处理,包括:
将各个测线方位的地震属性数据投影到与该测线方位关于原点中心对称的方位上。
3.如权利要求1所述地震纵波多方位属性椭圆拟合预测裂缝的方法,其特征在于,对生成的具有完备方位的地震属性数据进行椭圆拟合,包括:
将生成的具有完备方位的地震属性数据代入初始椭圆方程式:a(x-xc)2+b(x-xc)(y-yc)+c(y-yc)2=1;
计算获得拟合椭圆的方程式:ax2+bxy+cy2=1,
其中,a、b、c均是由椭圆的半长轴长度、椭圆的半短轴长度以及椭圆的长轴所在的方位值组成的函数,xc,yc分别为椭圆中心点的横纵坐标值,x,y是椭圆上点的横纵坐标值。
4.如权利要求1所述地震纵波多方位属性椭圆拟合预测裂缝的方法,其特征在于,所述裂缝发育强度的值域范围是[0,1)。
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=1,
其中,a、b、c均是由椭圆的半长轴长度、椭圆的半短轴长度以及椭圆的长轴所在的方位值组成的函数,xc,yc分别为椭圆中心点的横纵坐标值,x,y是椭圆上点的横纵坐标值。
12.如权利要求9所述地震纵波多方位属性椭圆拟合预测裂缝的装置,其特征在于,所述裂缝发育强度的值域范围是[0,1)。
13.如权利要求9所述地震纵波多方位属性椭圆拟合预测裂缝的装置,其特征在于,还包括:
第一计算模块,用于将所述椭圆的扁率确定为裂缝发育强度之后,计算每个面元的每个时间采样点上的裂缝发育强度,输出裂缝发育强度数据体。
14.如权利要求9至13中任一项所述地震纵波多方位属性椭圆拟合预测裂缝的装置,其特征在于,所述裂缝预测模块,包括:
计算单元,用于根据椭圆拟合结果计算椭圆的长轴和短轴所在的方位值;
预测单元,用于将椭圆的长轴或短轴所指示的方位确定为裂缝走向。
15.如权利要求14所述地震纵波多方位属性椭圆拟合预测裂缝的装置,其特征在于,所述预测单元还用于在所述地震属性与速度正相关的情况下,将所述椭圆的长轴指示的方位确定为裂缝走向;在所述地震属性与时间正相关的情况下,将所述椭圆的短轴指示的方位确定为裂缝走向。
16.如权利要求14所述地震纵波多方位属性椭圆拟合预测裂缝的装置,其特征在于,还包括:
第二计算模块,用于在将椭圆的长轴或短轴所指示的方位确定为裂缝走向之后,计算每个面元的每个时间采样点上的裂缝走向,输出裂缝走向数据体。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410064797.1A CN103869359A (zh) | 2014-02-25 | 2014-02-25 | 地震纵波多方位属性椭圆拟合预测裂缝的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410064797.1A CN103869359A (zh) | 2014-02-25 | 2014-02-25 | 地震纵波多方位属性椭圆拟合预测裂缝的方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103869359A true CN103869359A (zh) | 2014-06-18 |
Family
ID=50908092
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410064797.1A Pending CN103869359A (zh) | 2014-02-25 | 2014-02-25 | 地震纵波多方位属性椭圆拟合预测裂缝的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103869359A (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104155701A (zh) * | 2014-07-03 | 2014-11-19 | 中国石油天然气集团公司 | 一种利用叠前地震资料及井信息的多尺度裂缝预测方法 |
CN104818982A (zh) * | 2015-04-20 | 2015-08-05 | 中国石油天然气股份有限公司 | 一种缝洞型储层连通性定量分析方法及装置 |
CN105676289A (zh) * | 2016-03-31 | 2016-06-15 | 中国地质大学(武汉) | 一种通过前积角度玫瑰花图判断前积体主物源方向的方法 |
CN106501851A (zh) * | 2016-09-30 | 2017-03-15 | 中国石油天然气集团公司 | 一种地震属性优化方法及装置 |
CN106842299A (zh) * | 2016-12-19 | 2017-06-13 | 中国石油天然气股份有限公司 | 一种基于地震属性的裂缝定量化预测的方法 |
CN107203003A (zh) * | 2017-07-27 | 2017-09-26 | 河北煤炭科学研究院 | 一种矿井水害微震监测时空簇分析方法 |
CN108957527A (zh) * | 2017-05-27 | 2018-12-07 | 中国石油化工股份有限公司 | 岩层网状裂缝的地震预测方法 |
CN109188519A (zh) * | 2018-09-17 | 2019-01-11 | 中国石油大学(华东) | 一种极坐标下的弹性波纵横波速度反演系统及方法 |
CN109655908A (zh) * | 2017-10-11 | 2019-04-19 | 中国石油化工股份有限公司 | 基于HTI介质的方位AVAz的正演模拟方法及系统 |
CN109725350A (zh) * | 2018-12-05 | 2019-05-07 | 中国石油天然气集团有限公司 | 一种确定纵波方位各向异性参数的方法、装置及系统 |
CN111158053A (zh) * | 2019-12-20 | 2020-05-15 | 中石化石油工程技术服务有限公司 | 裂缝预测方法及装置 |
CN112649853A (zh) * | 2019-10-10 | 2021-04-13 | 中国石油化工股份有限公司 | 基于叠后数据的裂缝参数预测方法及系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101907725A (zh) * | 2010-06-12 | 2010-12-08 | 中国石油天然气股份有限公司 | 裂缝预测方法和装置 |
-
2014
- 2014-02-25 CN CN201410064797.1A patent/CN103869359A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101907725A (zh) * | 2010-06-12 | 2010-12-08 | 中国石油天然气股份有限公司 | 裂缝预测方法和装置 |
Non-Patent Citations (6)
Title |
---|
LI X Y ET AL: "Fracture detection using land 3D seismic data from the Yellow River Delta, China", 《THE LEADING EDGE》 * |
MARIA A.PEREZ ET AL: "Fracture detection in a carbonate reservoir using a variety of seismic methods", 《GEOPHYSICS》 * |
于海生等: "方位纵波裂缝检测技术在深层火山岩中的应用", 《新疆石油地质》 * |
刘军迎等: "基于叠前方位振幅的大港-埕海地区奥陶系风化壳裂缝储层的叠前预测", 《地球物理学进展》 * |
刘卫华等: "纵波地震叠前方位数据裂缝预测方法及应用", 《中国地球物理学会第二十七届年会论文集》 * |
钱忠平等: "HTI介质中地震各向异性模拟与流体预测", 《中国地球物理学会第二十八届年会论文集》 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104155701B (zh) * | 2014-07-03 | 2016-08-31 | 中国石油天然气集团公司 | 一种利用叠前地震资料及井信息的多尺度裂缝预测方法 |
CN104155701A (zh) * | 2014-07-03 | 2014-11-19 | 中国石油天然气集团公司 | 一种利用叠前地震资料及井信息的多尺度裂缝预测方法 |
CN104818982B (zh) * | 2015-04-20 | 2017-12-05 | 中国石油天然气股份有限公司 | 一种缝洞型储层连通性定量分析方法及装置 |
CN104818982A (zh) * | 2015-04-20 | 2015-08-05 | 中国石油天然气股份有限公司 | 一种缝洞型储层连通性定量分析方法及装置 |
CN105676289A (zh) * | 2016-03-31 | 2016-06-15 | 中国地质大学(武汉) | 一种通过前积角度玫瑰花图判断前积体主物源方向的方法 |
CN106501851A (zh) * | 2016-09-30 | 2017-03-15 | 中国石油天然气集团公司 | 一种地震属性优化方法及装置 |
CN106842299A (zh) * | 2016-12-19 | 2017-06-13 | 中国石油天然气股份有限公司 | 一种基于地震属性的裂缝定量化预测的方法 |
CN108957527A (zh) * | 2017-05-27 | 2018-12-07 | 中国石油化工股份有限公司 | 岩层网状裂缝的地震预测方法 |
CN107203003A (zh) * | 2017-07-27 | 2017-09-26 | 河北煤炭科学研究院 | 一种矿井水害微震监测时空簇分析方法 |
CN109655908A (zh) * | 2017-10-11 | 2019-04-19 | 中国石油化工股份有限公司 | 基于HTI介质的方位AVAz的正演模拟方法及系统 |
CN109188519A (zh) * | 2018-09-17 | 2019-01-11 | 中国石油大学(华东) | 一种极坐标下的弹性波纵横波速度反演系统及方法 |
CN109725350A (zh) * | 2018-12-05 | 2019-05-07 | 中国石油天然气集团有限公司 | 一种确定纵波方位各向异性参数的方法、装置及系统 |
CN109725350B (zh) * | 2018-12-05 | 2020-09-08 | 中国石油天然气集团有限公司 | 一种确定纵波方位各向异性参数的方法、装置及系统 |
CN112649853A (zh) * | 2019-10-10 | 2021-04-13 | 中国石油化工股份有限公司 | 基于叠后数据的裂缝参数预测方法及系统 |
CN112649853B (zh) * | 2019-10-10 | 2023-08-11 | 中国石油化工股份有限公司 | 基于叠后数据的裂缝参数预测方法及系统 |
CN111158053A (zh) * | 2019-12-20 | 2020-05-15 | 中石化石油工程技术服务有限公司 | 裂缝预测方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103869359A (zh) | 地震纵波多方位属性椭圆拟合预测裂缝的方法及装置 | |
Cvetič et al. | Photon spheres and sonic horizons in black holes from supergravity and other theories | |
CN101551466B (zh) | 一种利用与偏移距有关的地震属性提高油气储层预测精度的方法 | |
KR102021276B1 (ko) | 진폭 보존을 갖는 fwi 모델 도메인 각도 스택들 | |
CN104375182B (zh) | 一种裂缝型储层流体的识别方法及装置 | |
CN105629303A (zh) | 基于岩石物理的叠前裂缝定量预测方法及系统 | |
CN105425292A (zh) | 一种油气预测方法及装置 | |
CN104155693A (zh) | 储层流体流度的角道集地震响应数值计算方法 | |
CN102414581A (zh) | 包括利用多次波来确定地下区域的地下特性的基于波场的数据处理方法 | |
Wang et al. | Seismic prediction method of multiscale fractured reservoir | |
US10948617B2 (en) | Generating a velocity model for a subsurface structure using refraction travel time tomography | |
Zong et al. | Young’s modulus variation with azimuth for fracture-orientation estimation | |
WO2021126814A1 (en) | Mapping near-surface heterogeneities in a subterranean formation | |
CN103558637B (zh) | 基于三分量传感器的远探测方法 | |
CN106353807B (zh) | 裂缝识别方法和装置 | |
CN102053262A (zh) | 获得地震转换波的方位速度的方法及处理地震数据的方法 | |
Gao et al. | An Efficient Multiscale Finite‐Element Method for Frequency‐Domain Seismic Wave Propagation | |
CN104122584A (zh) | 根据地震数据确定方向性的方法及装置 | |
Protasov et al. | Characterization of fractured zones via topological analysis of 3D seismic diffraction images | |
CN112684498A (zh) | 一种基于宽方位地震数据的储层裂缝预测方法及系统 | |
US11353609B2 (en) | Identifying geologic features in a subterranean formation using angle domain gathers sampled in a spiral coordinate space | |
US11460595B2 (en) | Unified continuous seismic reservoir monitoring | |
Ursin et al. | Fresnel volume and interface Fresnel zone for reflected and transmitted waves from a curved interface in anisotropic media | |
Grechka | On the nonuniqueness of traveltime inversion in elliptically anisotropic media | |
Al-Gawas et al. | Fracture detection via correlating P-wave amplitude variation with offset and azimuth analysis and well data in eastern central Saudi Arabia |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20140618 |