CN104155701B - 一种利用叠前地震资料及井信息的多尺度裂缝预测方法 - Google Patents
一种利用叠前地震资料及井信息的多尺度裂缝预测方法 Download PDFInfo
- Publication number
- CN104155701B CN104155701B CN201410315202.5A CN201410315202A CN104155701B CN 104155701 B CN104155701 B CN 104155701B CN 201410315202 A CN201410315202 A CN 201410315202A CN 104155701 B CN104155701 B CN 104155701B
- Authority
- CN
- China
- Prior art keywords
- fracture
- crack
- scale
- orientation
- spacing
- 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.)
- Active
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种利用叠前地震资料及井信息的多尺度裂缝预测方法,包括:将叠前地震数据分选为分方位道集;利用每一方位上计算的衰减品质因子形成新的道集;利用每一新的道集拟合椭圆,根据椭圆的长轴方向、长轴与短轴之比确定裂缝的预测裂缝方位和预测裂缝密度;对勘探目标区域进行成像测井解释,得到测井井点处各裂缝的实际裂缝密度;将实际裂缝密度投影到裂缝尺度‑‑理论裂缝密度‑‑平均裂缝孔隙度对应关系图中,根据投影的结果确定裂缝密度的有效分布范围,进而对预测裂缝密度进行优化,得到勘探目标区域的有效裂缝密度。本发明能够实现对不同尺度裂缝的密度和方位信息的分辨,在裂缝预测中利用井信息实现裂缝密度、方位和尺度的综合预测。
Description
技术领域
本发明涉及石油勘探技术领域,具体地,涉及一种利用叠前地震资料及井信息的多尺度裂缝预测方法。
背景技术
地下岩石中广泛存在裂缝,裂缝型油藏是重要的油气勘探目标。在油气勘探中的主要技术是地震勘探,Ruger等研究了众多方法可以通过反射地震预测地下裂缝分布,这些方法被实际应用证明能够对地下裂缝分布进行预测。但是,实际上地下裂缝分布是呈现多组、多尺度分布特征,而现有技术只能体现所有裂缝的综合信息,不能有效的区分不同尺度裂缝。因此,现有技术预测裂缝实际上并不是地下裂缝的真实分布特征,而是一种对地下裂缝特征的等效显示。
目前,常用的裂缝预测技术分为横波分裂、微地震和纵波裂缝预测技术,由于横波采集成本高,微地震监测范围有限等问题,目前最为常用的技术是纵波裂缝预测技术,该技术通常根据地震资料类型分为叠前和叠后两类,叠后主要通过相干、曲率等几何属性研究地震道横向不连续性,通过相干等不连续属性研究由于断层、小断裂等大裂缝特征。另外我们可以通过叠前地震资料方位特征进行裂缝预测,叠前纵波方位各向异性裂缝预测方法主要通过计算不同方位属性,然后对不同方位属性进行椭圆拟合,椭圆长轴对应裂缝方位,椭圆短长轴比对应裂缝密度的方式进行裂缝预测,根据方位属性的不同,可以将叠前裂缝预测方法分为速度各向异性技术(VVAZ),振幅各向异性技术(AVAZ)以及衰减各向异性技术(QVAZ),这三类技术分别通过研究不同方位的速度、振幅和衰减属性(品质因子Q)的差异来预测裂缝方位和密度。但是目前所有这些技术的实现都只是将井信息(如成像测井得到裂缝方位等)作为最终验证结果正确性的参考,而没有将井信息加入到实际裂缝预测中进行有效利用,另外没有区分出所预测的裂缝究竟尺度为多大(延伸长度),因此我们有必要在现有技术基础上做进一步的发展。
发明内容
本发明实施例的主要目的在于提供一种利用叠前地震资料及井信息的多尺度裂缝预测方法,以解决现有的叠前裂缝预测技术所预测裂缝信息不能分辨尺度的问题,以及井信息无法参与裂缝预测的问题。
为了实现上述目的,本发明实施例提供一种利用叠前地震资料及井信息的多尺度裂缝预测方法,包括:
获取勘探目标区域的叠前地震数据,并将所述叠前地震数据分选为分方位道集;
在每一方位上计算共中心点道集的衰减品质因子,所有方位执行之后,将同一中心反射点在不同方位的衰减品质因子形成新的道集;
利用所述每一新的道集拟合椭圆,将拟合得到的椭圆的长轴方向确定为该新的道集所对应的中心反射点处裂缝的预测裂缝方位,将椭圆的长轴与短轴之比确定为该新的道集所对应的中心反射点处裂缝的预测裂缝密度,利用所有新的道集拟合椭圆之后,得到勘探目标区域中各个中心反射点处裂缝的预测裂缝方位和预测裂缝密度;
获取勘探目标区域所有井的成像测井资料,针对每一测井井点,对其成像测井资料进行成像测井解释,得到该测井井点处各裂缝的实际裂缝密度;将每一测井井点处裂缝的实际裂缝密度投影到该勘探目标区域的裂缝尺度--理论裂缝密度--平均裂缝孔隙度对应关系图中,根据投影结果,确定该勘探目标区域中裂缝密度的有效分布范围,进而,利用该有效分布范围对该勘探目标区域中各个中心反射点处裂缝的预测裂缝密度进行优化,得到勘探目标区域的有效裂缝密度;
其中,所述勘探目标区域的裂缝尺度--理论裂缝密度--平均裂缝孔隙度对应关系图是按照如下步骤确定:
获取勘探目标区域所有井的成像测井资料;针对每一测井井点,对其成像测井资料进行成像测井解释,得到该测井井点处各裂缝的孔隙度,求取该测井井点所有裂缝的孔隙度的平均值,得到平均裂缝孔隙度;
针对每一测井井点,利用其成像测井资料,采用多尺度裂缝岩石物理模型,计算每一尺度的裂缝在平行裂缝方位的衰减系数和垂直裂缝方位的衰减系数;然后将每一尺度的裂缝在平行裂缝方位的衰减系数和垂直裂缝方位的衰减系数的比值,确定为该尺度裂缝的理论裂缝密度;
将勘探目标区域所有测井井点中裂缝的各种尺度、每种尺度裂缝的理论裂缝密度、平均裂缝孔隙度记录至以裂缝尺度、裂缝密度和裂缝孔隙度为坐标的图系中,得到该勘探目标区域的裂缝尺度--理论裂缝密度--平均裂缝孔隙度对应关系图。
借助于上述技术方案,本发明能够实现对不同尺度裂缝的密度和方位信息的分辨,可以在裂缝预测中综合利用井信息,实现裂缝密度、方位和尺度信息的综合预测,提高了裂缝预测的精度和可靠度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本发明提供的多尺度裂缝预测方法流程图;
图2是本发明实施例提供的分方位道集;
图3是本发明实施例提供的椭圆拟合用新道集;
图4(a)、4(b)是本发明实施例提供的勘探目标区域中各个中心反射点处裂缝的预测裂缝密度和预测裂缝方位图;
图5是本发明实施例提供的某油田H区块成像测井解释成果统计的裂缝孔隙度值;
图6是本发明实施例提供的采用实际地震数据H区块Ha601-4井平均裂缝孔隙度计算得到裂缝尺度从0.01m-1.4m时的理论裂缝密度;
图7是本发明实施例提供的裂缝尺度--理论裂缝密度--平均裂缝孔隙度对应关系图;
图8是本发明实施例提供的各个成像测井井点位置处的实际裂缝密度;
图9是本发明实施例提供的优化后的预测裂缝密度;
图10是本发明实施例提供的预测裂缝密度与裂缝尺度的交会图;
图11是本发明实施例提供的根据有效裂缝密度与裂缝尺度的关系函数转换得到的裂缝尺度分布图;
图12(a)、12(b)、12(c)是本发明提供的纵波振幅、速度及衰减随方位角变化的示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在研究本发明的过程中,为了利用叠前地震资料进行多尺度裂缝预测,需要选择合适的纵波属性以区分不同尺度裂缝的响应特征。如图1所示,研究人员设计了一组模型实验进行对比分析,图12中固定入射角30°,定向排列的方位为50°裂缝,分别得到如图12(a)、(b)、(c)所示的纵波振幅、速度及衰减随方位角变化的情况。根据图1可知,振幅、速度极大值与衰减极小值都可指示裂缝方位,这说明利用振幅、速度和衰减的各向异性来分析裂缝方位是可靠的;此外,图1还显示出衰减随方位变化的特征随着裂缝尺度增大而逐渐增强,而且衰减的各向异性强度随尺度增加也逐渐增强,相对而言裂缝尺度越大,衰减的各向异性越强,裂缝尺度越小,衰减的各向异性特征越不明显,但是对于不同尺度的裂缝,振幅、速度的各向异性差异很小。这说明衰减各向异性对不同尺度裂缝所体现出的差异更加明显,也就是说,如果在叠前计算衰减属性,然后通过椭圆拟合得到裂缝的各向异性强度(裂缝密度),不同尺度裂缝所表征出的各向异性强度值差异是比较明显的,如果通过确定不同各向异性强度所对应的裂缝尺度,就可以进一步的得到裂缝尺度的分布范围。
研究表明,不同尺度的裂缝会表现出不同的波场特征,大于波长尺度的裂缝更多的表现为反射特征,尺度小于或者接近波长的裂缝更多的表现为散射特征。因此,通过偏移和叠加可能会消除尺度小于或者接近波长裂缝引起的波场变化,也就是说叠后地震资料预测更多的是大尺度裂缝特征,而叠前地震资料包含了不同尺度裂缝的信息,采用叠前地震资料进行裂缝预测会得到更丰富的预测结果。另外,现有裂缝预测技术更多的只是采用了地震资料本身所提取的属性进行预测,并没有井信息的参与,而井的分辨率和精确度要明显高于地震资料,但是井信息覆盖面积小,地震信息虽然精度不如井资料,但是覆盖面积大。因此,通过井信息与地震资料联合可以综合两者的优势,提高地震预测的精度和可靠性。岩石物理技术是连接岩石物理属性和地震波属性的桥梁,井信息能够直接对岩石的物理属性进行测量,可以通过岩石物理分析,将井信息转化为地震预测可用的信息,对地震预测结果进行精细化标定,明确地震预测信息更深层次的含义。
基于以上研究,本发明提供一种利用叠前地震资料和井信息,综合预测裂缝尺度、裂缝密度以及裂缝方位的方法。以下结合附图与实例,对本发明提供的利用叠前地震资料及井信息的多尺度裂缝预测方法进行说明。
本发明提供的利用叠前地震资料及井信息的多尺度裂缝预测方法,如图1所示,包括如下步骤:
步骤S1,获取勘探目标区域的叠前地震数据,并将这些叠前地震数据分选为分方位道集。如附图2所示为实际地震数据Inline990原始振幅方位角道集。
步骤S2,在每一方位上计算共中心点道集的衰减品质因子,所有方位执行之后,将同一中心反射点在不同方位的衰减品质因子形成新的道集,即,每一新的道集为对应中心反射点在不同方位的衰减品质因子。
衰减品质因子是地震勘探中用于表征地震波衰减属性的一个指标,衰减品质因子与衰减系数呈反比,衰减品质因子越大,地震波衰减越小,衰减品质因子越小,地震波衰减越大。最初衰减品质因子计算主要通过VSP和叠后地震记录实现,Clark在1999年提出在叠前计算衰减品质因子的方法,称为QVO方法,该方法考虑到品质因子随偏移距的变化,通过分析不同偏移距处品质因子与偏移距平方呈线性关系去计算每一层某一固定方位的衰减品质因子,这样得到每一方位的衰减品质因子就可以用于分析各向异性特征。
步骤S2即采用QVO方法计算衰减品质因子,具体过程如下:
针对每一方位,在每一共中心点道集中,选择目标层及其对应的参照道,较佳的,选择目标层顶部反射波能量强、波形好的道(通常为第一道)作为参照道。
步骤S21,计算目标层顶界面以上介质的衰减品质因子平均值,包括如下三步:
(1)将目标层顶界面反射波的频谱与参照道的频谱进行比值运算,频谱比与衰减品质因子具有如下关系式:
其中,A1(f)是目标层顶界面反射波的振幅值,Ar(f)是参照道的振幅值,f是频率,R是反射系数,G示几何发散因子,t1是目标位置旅行时目标层顶界面旅行时,tr是参照道的旅行时,Q是衰减品质因子。
(2)采用如下公式,将该共中心点道集中每道的频谱比值进行线性拟合,求出各自对数谱比值曲线的斜率:
k=-π(tr-t1)/Q (公式2)
其中,k是对数谱比值曲线的斜率。
为了将公式2中的tr和t1变为计算更方便准确的零偏移距旅行时,根据如下公式分析斜率:
其中,t0,r、t0,1分别是目标层顶界面和参照道的零炮检距旅行时,即自激自收时间,Vr、V1分别是目标层和参照道均方根速度,x表示炮检距。
(3)由公式3可知,通过对随炮检距平方变化的对数谱比值曲线的斜率k进行线性拟合求出截距l,利用以下公式就可以计算得到目标层顶界面以上的衰减品质因子平均值:
l=-π(t0,r-t0,1)/Q (公式4)
步骤S22,对目标层底界面重复步骤S21中的(1)、(2)、(3)三步,可以得到目标层底界面以上的衰减品质因子平均值。
步骤S23,利用已经计算出的目标层顶界面/底界面以上的衰减品质因子平均值,采用以下公式可求出目标层的衰减品质因子:
其中,Qn是目标层顶界面以上的衰减品质因子平均值,Qn-1是目标层底界面以上的衰减品质因子平均值,Qi是目标层的衰减品质因子,tn是目标层顶界面的零炮检距旅行时,tn-1是目标层底界面的零炮检距旅行时。
通过以上步骤可以得到一个方位上一个共中心点道集的目标层衰减品质因子,重复此过程到所有方位的所有共中心点道集进行计算,所有方位执行之后,再将同一中心反射点在不同方位的衰减品质因子形成新的道集,如附图3所示,图3中每一道集为同一中心反射点在不同方位的衰减品质因子。
步骤S3,利用步骤S2中形成的每一新的道集拟合椭圆,将拟合得到的椭圆的长轴方向确定为该新的道集所对应的中心反射点处裂缝的预测裂缝方位,将椭圆的长轴与短轴之比确定为该新的道集所对应的中心反射点处裂缝的预测裂缝密度,利用所有新的道集拟合椭圆之后,得到勘探目标区域中各个中心反射点处裂缝的预测裂缝方位和预测裂缝密度。
如附图4所示,附图4(a)、4(b)分别为勘探目标区域中各个中心反射点处裂缝的预测裂缝密度和预测裂缝方位。
目前常规叠前裂缝预测方法大多就认为此结果就是最终的裂缝预测结果,但是由于图形显示不直观,预测结果裂缝分布范围过于广泛,而且无法确认哪一区域裂缝是我们所需要的储层裂缝,因此给后期的解释造成了困难。本专利通过成像测井解释成果与地震预测结果综合应用来解决该问题。
步骤S4,获取勘探目标区域所有井的成像测井资料;针对每一测井井点,对其成像测井资料进行成像测井解释,得到该测井井点处各裂缝的实际裂缝密度;将每一测井井点处裂缝的实际裂缝密度投影到该勘探目标区域的裂缝尺度--理论裂缝密度--平均裂缝孔隙度对应关系图中,根据投影结果,确定该勘探目标区域中裂缝密度的有效分布范围,进而,利用该有效分布范围对步骤S3所得的勘探目标区域中各个中心反射点处裂缝的预测裂缝密度进行优化,得到勘探目标区域的有效裂缝密度,相比于步骤S3所得的预测裂缝密度,有效裂缝密度更加精确。
将图8所示的各个成像测井井点处的实际裂缝密度投影到各向异性强度-裂缝尺度-裂缝孔隙度量板中(如图7所示的散点),由于井点位置处可以得知地震计算的各向异性强度和裂缝孔隙度,这样就可以从量板中确定出裂缝各向异性强度的有效分布范围(如果各向异性强度不能投影到该量板中视为无效),可以根据该分布范围确定门槛值将原来预测的各向异性强度结果进行优化,所得到的就是有效裂缝各向异性强度分布图。在本例中,通过测井裂缝值投影确定门槛值大约为1.2-2(图7中黑色点确定的彩色范围最大值和最小值),然后对步骤S4中的叠前裂缝各向异性强度预测结果确定色标有效范围为1.2-2,其他结果属于冗余信息,统一将色标变为1(因为短轴与长轴比为1表示不具有各向异性),所得结果如图9所示,与图4相比,裂缝的分布范围得到了优化,去掉了冗余信息,使得结果更加精确,便于后期解释。
步骤S4中所用的勘探目标区域的裂缝尺度--理论裂缝密度--平均裂缝孔隙度对应关系图是按照如下步骤确定:
步骤S41,获取勘探目标区域所有井的成像测井资料;针对每一测井井点,对其成像测井资料进行成像测井解释,得到该测井井点处各裂缝的孔隙度,求取该测井井点所有裂缝的孔隙度的平均值,得到平均裂缝孔隙度。
选取勘探目标区域内的成像测井资料(如微电阻率扫描成像FWI),成像测井资料能够对井点位置岩石特征进行精确成像,并且能够从中提取中更为精确的裂缝信息(裂缝方位和裂缝密度),通过成像测井解释可以得到裂缝孔隙度,考虑到成像测井单个裂缝孔隙度值覆盖深度较小,地震资料难以达到,因此考虑将所测储层所有深度裂缝孔隙度进行平均,以得到平均裂缝孔隙度,该裂缝孔隙度范围可以覆盖储层段地震预测结果,如附图5为某油田H区块(横坐标为井点名称、纵坐标为裂缝孔隙度)成像测井解释成果统计的裂缝孔隙度值。
步骤S42,针对每一测井井点,利用其成像测井资料,采用多尺度裂缝岩石物理模型,计算每一尺度的裂缝在平行裂缝方位的衰减系数和垂直裂缝方位的衰减系数;然后将每一尺度的裂缝在平行裂缝方位的衰减系数和垂直裂缝方位的衰减系数的比值,确定为该尺度裂缝的理论裂缝密度。
该步骤采用多尺度裂缝岩石物理模型(Chapman模型)计算裂缝在每一尺度条件下平行裂缝方位和垂直裂缝方位的衰减系数,Chapman提出的多尺度定向排列裂缝模型(Chapman模型)考虑了两种尺度裂缝流体的相互作用,该模型弹性刚度张量可以写为如下形式:
(公式6)
其中,C是弹性刚度张量,是裂缝孔隙度,ec是微裂隙密度,ef是裂缝密度,C0是弹性刚度矩阵,C1、C2、C3分别表征孔隙、裂缝、小裂纹的贡献。利用该模型就可以通过改变裂缝弹性张量C3和裂缝密度ef得到不同尺度裂缝的弹性模量,从而计算出裂缝介质的衰减系数、速度等参数。
平行裂缝方位属性最大,垂直裂缝方位属性最小,两者之比能够等效于椭圆拟合后的长轴与短轴之比,也就是可以表征裂缝的各向异性强度,因此采用岩石物理模型计算得到的垂直裂缝方位衰减系数与平行裂缝方位衰减系数之比表征理论计算的各向异性强度值(即理论裂缝密度)。如图6所示,采用前述实际地震数据H区块Ha601-4井平均裂缝孔隙度计算得到裂缝尺度从0.01m-1.4m时的平行和垂直裂缝方向衰减之比,该比值视为固定裂缝孔隙度条件下不同尺度裂缝的裂缝密度值。
步骤S43,将勘探目标区域所有测井井点中裂缝的各种尺度、每种尺度裂缝的理论裂缝密度、平均裂缝孔隙度记录至以裂缝尺度、裂缝密度和裂缝孔隙度为坐标的图系中,得到该勘探目标区域的裂缝尺度--理论裂缝密度--平均裂缝孔隙度对应关系图,也就是说,该对应关系图记录了该勘探目标区域所有测井井点中裂缝的各种尺度、每种尺度裂缝的理论裂缝密度、平均裂缝孔隙度的对应关系。
如图7所示为裂缝尺度--理论裂缝密度--平均裂缝孔隙度对应关系图,其中颜色代表理论裂缝密度,横坐标代表裂缝尺度,纵坐标代表裂缝孔隙度。
此外,本发明提供的利用叠前地震资料及井信息的多尺度裂缝预测方法还包括如下步骤:
步骤5,根据裂缝尺度--理论裂缝密度--平均裂缝孔隙度对应关系图,确定勘探目标区域的有效裂缝密度与裂缝尺度的对应关系,对该对应关系进行曲线拟合得到有效裂缝密度与裂缝尺度的关系函数;根据该关系函数将步骤S3得到的预测裂缝密度转换为预测裂缝尺度,进而确定该勘探目标区域的裂缝尺度分布范围。
首先在附图7中确定出测井井点位置裂缝的预测裂缝密度与裂缝尺度,交会于图10所示的图中并进行曲线拟合,可以得到地震预测裂缝密度与裂缝尺度的函数关系以及裂缝尺度有效预测范围,这样就可以根据该函数关系将裂缝密度分布图(附图9)转换为裂缝尺度分布图。如附图11所示,根据附图10分析结果,裂缝尺度大致范围为1.1m-4m,裂缝密度与裂缝尺度对应关系为:
y=-0.9174x2+6.533x-7.818
其中,x表示裂缝密度,y表示裂缝尺度。
至此,综合利用测井资料、地震资料综合分析得到了裂缝密度、裂缝尺度和裂缝方位信息。通过综合对比三种预测结果可以指导裂缝型储层的勘探工作,本专利的实施也为提高裂缝预测精度,增加裂缝预测信息,综合利用岩石物理和井信息进行裂缝预测提供了技术方案。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (2)
1.一种利用叠前地震资料及井信息的多尺度裂缝预测方法,其特征在于,包括:
获取勘探目标区域的叠前地震数据,并将所述叠前地震数据分选为分方位道集;
在每一方位上计算共中心点道集的衰减品质因子,所有方位执行之后,将同一中心反射点在不同方位的衰减品质因子形成新的道集;
利用所述每一新的道集拟合椭圆,将拟合得到的椭圆的长轴方向确定为该新的道集所对应的中心反射点处裂缝的预测裂缝方位,将椭圆的长轴与短轴之比确定为该新的道集所对应的中心反射点处裂缝的预测裂缝密度,利用所有新的道集拟合椭圆之后,得到勘探目标区域中各个中心反射点处裂缝的预测裂缝方位和预测裂缝密度;
获取勘探目标区域所有井的成像测井资料,针对每一测井井点,对其成像测井资料进行成像测井解释,得到该测井井点处各裂缝的实际裂缝密度;将每一测井井点处裂缝的实际裂缝密度投影到该勘探目标区域的裂缝尺度--理论裂缝密度--平均裂缝孔隙度对应关系图中,根据投影结果,确定该勘探目标区域中裂缝密度的有效分布范围,进而,利用该有效分布范围对该勘探目标区域中各个中心反射点处裂缝的预测裂缝密度进行优化,得到勘探目标区域的有效裂缝密度;
其中,所述勘探目标区域的裂缝尺度--理论裂缝密度--平均裂缝孔隙度对应关系图是按照如下步骤确定:
获取勘探目标区域所有井的成像测井资料;针对每一测井井点,对其成像测井资料进行成像测井解释,得到该测井井点处各裂缝的孔隙度,求取该测井井点所有裂缝的孔隙度的平均值,得到平均裂缝孔隙度;
针对每一测井井点,利用其成像测井资料,采用多尺度裂缝岩石物理模型,计算每一尺度的裂缝在平行裂缝方位的衰减系数和垂直裂缝方位的衰减系数;然后将每一尺度的裂缝在平行裂缝方位的衰减系数和垂直裂缝方位的衰减系数的比值,确定为该尺度裂缝的理论裂缝密度;
将勘探目标区域所有测井井点中裂缝的各种尺度、每种尺度裂缝的理论裂缝密度、平均裂缝孔隙度记录至以裂缝尺度、理论裂缝密度和平均裂缝孔隙度为坐标的图系中,得到该勘探目标区域的裂缝尺度--理论裂缝密度--平均裂缝孔隙度对应关系图。
2.根据权利要求1所述的方法,其特征在于,所述方法还包括:
根据裂缝尺度--理论裂缝密度--平均裂缝孔隙度对应关系图,确定勘探目标区域的有效裂缝密度与裂缝尺度的对应关系,对该对应关系进行曲线拟合得到有效裂缝密度与裂缝尺度的关系函数;根据该关系函数将所述预测裂缝密度转换为预测裂缝尺度,进而确定该勘探目标区域的裂缝尺度分布范围。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410315202.5A CN104155701B (zh) | 2014-07-03 | 2014-07-03 | 一种利用叠前地震资料及井信息的多尺度裂缝预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410315202.5A CN104155701B (zh) | 2014-07-03 | 2014-07-03 | 一种利用叠前地震资料及井信息的多尺度裂缝预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104155701A CN104155701A (zh) | 2014-11-19 |
CN104155701B true CN104155701B (zh) | 2016-08-31 |
Family
ID=51881251
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410315202.5A Active CN104155701B (zh) | 2014-07-03 | 2014-07-03 | 一种利用叠前地震资料及井信息的多尺度裂缝预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104155701B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105445814B (zh) * | 2015-11-09 | 2018-07-03 | 恒泰艾普集团股份有限公司 | 海上窄方位角资料的方位各向异性属性差裂缝预测方法 |
CN105425292B (zh) * | 2015-11-20 | 2019-04-12 | 中国石油天然气集团公司 | 一种油气预测方法及装置 |
CN106772680B (zh) * | 2016-12-30 | 2019-03-15 | 中国石油天然气股份有限公司 | 基于二维叠后地震资料微观裂缝表征确定方法和装置 |
CN109655908A (zh) * | 2017-10-11 | 2019-04-19 | 中国石油化工股份有限公司 | 基于HTI介质的方位AVAz的正演模拟方法及系统 |
CN109165737B (zh) * | 2018-08-29 | 2020-09-01 | 电子科技大学 | 基于条件随机场和bp神经网络的孔隙度预测方法 |
CN111679318B (zh) * | 2019-03-11 | 2022-12-02 | 中国石油天然气股份有限公司 | 储层中断裂系统的识别方法及系统 |
CN110656933B (zh) * | 2019-10-08 | 2022-12-13 | 核工业北京地质研究院 | 一种破碎带裂隙密度确定方法及系统 |
CN112731556B (zh) * | 2019-10-28 | 2023-06-30 | 中国石油化工股份有限公司 | 一种裂缝发育区的预测方法及用于预测裂缝发育区的计算机存储介质 |
CN112068197B (zh) * | 2020-08-24 | 2022-03-18 | 西南科技大学 | 基于地震波动力学参数的裂缝破裂尺度的描述方法 |
CN117951857A (zh) * | 2022-10-21 | 2024-04-30 | 中国石油天然气股份有限公司 | 预测原位构造裂缝延伸长度的方法及装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102654586A (zh) * | 2012-05-18 | 2012-09-05 | 北京师范大学 | 一种多尺度裂缝建模方法 |
CN103869359A (zh) * | 2014-02-25 | 2014-06-18 | 中国石油天然气股份有限公司 | 地震纵波多方位属性椭圆拟合预测裂缝的方法及装置 |
-
2014
- 2014-07-03 CN CN201410315202.5A patent/CN104155701B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102654586A (zh) * | 2012-05-18 | 2012-09-05 | 北京师范大学 | 一种多尺度裂缝建模方法 |
CN103869359A (zh) * | 2014-02-25 | 2014-06-18 | 中国石油天然气股份有限公司 | 地震纵波多方位属性椭圆拟合预测裂缝的方法及装置 |
Non-Patent Citations (5)
Title |
---|
P波资料反演裂缝方法及实例;马中高;《石油地球物理勘探》;20031031;第38卷(第5期);517-521 * |
The analysis of frequency-dependent characteristics for fl uid detection: a physical model experiment;Chen Shuang-Quan,et al.;《APPLIED GEOPHYSICS》;20120630;第9卷(第2期);195-206 * |
利用电成像评价塔里木山前致密砂岩裂缝有效性;陈胜;《中国优秀硕士学位论文全文数据库 基础科学辑》;20130115(第01期);正文第20-22页 * |
基于叠前方位振幅的大港—埕海地区奥陶系风化壳裂缝储层的叠前预测;刘军迎,等;《地球物理学进展》;20120831;第27卷(第4期);第1588-1597页 * |
松辽盆地南部营城组火山岩裂缝的叠前地震识别;李瑞磊,等;《成都理工大学学报(自然科学版)》;20121231;第39卷(第6期);第611-616页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104155701A (zh) | 2014-11-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104155701B (zh) | 一种利用叠前地震资料及井信息的多尺度裂缝预测方法 | |
US10386513B2 (en) | Seismic survey analysis | |
CN104635269B (zh) | 一种基于叠前方位角道集预测火成岩裂缝型储层的方法 | |
CN104237940B (zh) | 一种基于动力学特征的绕射波成像方法及装置 | |
CN105182424B (zh) | 一种基于斑块饱和模型定量预测储层孔隙度的方法和装置 | |
CN103733089B (zh) | 用于包括不确定性估计的地下表征的系统和方法 | |
CN107045145B (zh) | 地震层序控制下的叠前振幅随偏移距变化检测缝洞方法 | |
CN105425292A (zh) | 一种油气预测方法及装置 | |
CN103487831B (zh) | Avo地震正演计算方法 | |
CN104142516B (zh) | 一种薄单砂层厚度预测方法 | |
CN104678434A (zh) | 一种预测储层裂缝发育参数的方法 | |
CN104636980A (zh) | 针对河道砂油藏类型油气汇集条件的地球物理表征方法 | |
CN105445814B (zh) | 海上窄方位角资料的方位各向异性属性差裂缝预测方法 | |
CN105005080A (zh) | 一种利用振幅比率属性识别地层圈闭尖灭线的方法 | |
CN101276001B (zh) | 地下非均匀介质地震探测复杂性定量评估方法 | |
CN107728204A (zh) | 基于叠前纵波各向异性的裂缝预测方法及系统 | |
CN107807410A (zh) | 储层地球物理敏感参数优选方法及系统 | |
CN104834008B (zh) | 计算碳酸盐岩缝洞型储层充填程度的方法 | |
CN104678432A (zh) | 砂砾岩裂缝识别的方法 | |
CN102096098A (zh) | 地震属性提取方法 | |
CN106033127B (zh) | 基于横波速度变化率的地应力方位地震预测方法 | |
CN108957527A (zh) | 岩层网状裂缝的地震预测方法 | |
CN112946751B (zh) | 一种预测不同储层平面分布的方法及系统 | |
Mahani et al. | Determination of local magnitude for induced earthquakes in the Western Canada Sedimentary Basin: an update | |
CN106353807A (zh) | 裂缝识别方法和装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |