CN114167515A - 一种岩性圈闭有效性识别的方法 - Google Patents
一种岩性圈闭有效性识别的方法 Download PDFInfo
- Publication number
- CN114167515A CN114167515A CN202010953383.XA CN202010953383A CN114167515A CN 114167515 A CN114167515 A CN 114167515A CN 202010953383 A CN202010953383 A CN 202010953383A CN 114167515 A CN114167515 A CN 114167515A
- Authority
- CN
- China
- Prior art keywords
- target stratum
- stratum
- longitudinal wave
- wave impedance
- target
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种岩性圈闭有效性识别的方法,涉及油气开发技术领域,包括:S1:获取目标地层中盖层的封堵性表征关系式,S2:获取目标地层的盖层厚度分布图,S3:获得目标地层的封堵层厚度分布图,S4:获得目标地层的曲率分布图,S5:建立所述目标地层的岩性圈闭有效性参数计算式,并根据所述目标地层的盖层厚度有效界限值、封堵层厚度有效界限值、曲率有效界限值得到所述目标地层的岩性圈闭有效性阈值,得到所述目标地层的有效圈闭发育区分布图。本发明公开的一种岩性圈闭有效性识别的方法弥补了现有技术中岩性圈闭有效性识别时未考虑地层裂缝和断裂影响的问题,提高了岩性圈闭有效性识别的准确性,为油气勘探提供了更可靠的参考资料。
Description
技术领域
本发明涉及石油与天然气勘探开发研究领域,特别涉及一种岩性圈闭有效性识别的方法。
背景技术
在油气勘探中,为了提高钻井钻遇率、提升油气勘探开发经济效益,在钻井部署之前需要对已发现的岩性圈闭有效性进行识别。盖层和封堵层是岩性圈闭发育的必要条件,而它们的识别往往依赖于其突破压力和排驱压力的测定,对于钻井和取心较少的区,如何获取相关数据成为了一个难题;同时,获取了单井的测试数据以后,如何与地震资料结合,进而明确封堵层和盖层在平面上的展布特征,是岩性圈闭有效性判别的关键;而后期断裂和裂缝的发育情况对岩性圈闭的有效性也有重要影响。
但目前岩性圈闭有效性识别方法研究却少有开展裂缝对圈闭有效性影响的分析,不能很好的将地质-地震结合、没有考虑到裂缝对岩性圈闭的影响,影响岩性圈闭有效性识别的准确性。
发明内容
本申请的目的在于克服现有技术中所存在的岩性圈闭有效性识别中没有考虑裂缝对岩性圈闭的影响的问题,提供一种岩性圈闭有效性识别的方法。
为了实现上述发明目的,本申请提供了以下技术方案:一种岩性圈闭有效性识别的方法,包括以下步骤:
获取目标地层内不同地层深度对应的突破压力和排驱压力;
根据所述目标地层内已有钻井的纵波声波时差测井曲线和密度测井曲线,获得不同地层深度的纵波阻抗数值IMP;建立纵波阻抗数值与排驱压力的关系式、纵波阻抗数值与突破压力的关系式,并获得所述目标地层中盖层的封堵性表征关系式;
根据所述目标地层中盖层的岩性分布,获得不同岩性的盖层的纵波阻抗数值分布范围;
根据所述目标地层的三维地震数据体得到所述目标地层的纵波阻抗三维数据体,并根据所述目标地层中盖层的所述纵波阻抗数值的分布范围筛选不同岩性的盖层,得到所述目标地层中的盖层厚度分布图;
根据所述封堵性表征关系式和所述三维地震数据体获得封堵性参数三维数据体,并根据所述封堵性参数三维数据体,筛选出合格的封堵层,得到所述目标地层的封堵层厚度有效界限值和所述目标地层的封堵层厚度分布图;
根据所述目标地层的三维地震数据体,计算所述目标地层的曲率有效界限值和所述目标地层的曲率分布图;
建立所述目标地层的岩性圈闭有效性参数计算式,并根据所述目标地层的盖层厚度有效界限值、封堵层厚度有效界限值、曲率有效界限值得到所述目标地层的岩性圈闭有效性阈值,得到所述目标地层的有效圈闭发育区分布图;所述岩性圈闭有效性参数计算式为:
其中,所述f(C)为盖层厚度;所述f(S)为封堵层厚度,f(F)为曲率。
在上述技术方案中,在建立岩性圈闭有效性参数计算式过程中引入用于表征所述目标地层裂缝和断裂分布特征的曲率,弥补了现有技术中岩性圈闭有效性识别时未考虑地层裂缝和断裂影响的问题,提高了岩性圈闭有效性识别的准确性,为油气勘探提供了更可靠的参考资料;此外,上述技术方案中,还通过所述目标地层的三维地震数据体结合封堵性表征关系式来计算所述目标地层的封堵层分布、裂缝和断裂分布,减小了由于所述目标地层内的已有钻井数量较少使获得的突破压力和排驱压力数据量较少而导致的应用范围局限,解决了已有钻井数量较少、测井数据较缺乏的地层区域内难以开展岩性圈闭有效性评价工作的问题。
进一步地,所述封堵层厚度有效界限值为5m;所述盖层分为泥页岩层、膏岩层和致密灰岩层;当所述盖层为泥页岩层时,所述盖层厚度有效界限值为5m;当所述盖层为膏岩层时,所述盖层有效界限值为5m,当所述盖层为致密灰岩层时,所述盖层有效界限值为100m;所述曲率有效界限值为所述目标地层内断层处曲率的最小值。
进一步地,所述目标地层的有效圈闭发育区分布图获得步骤如下:
根据所述岩性圈闭有效性参数计算式,并带入所述目标地层的所述盖层厚度f(C)、所述封堵层厚度f(S)、所述曲率f(F),获得所述目标地层内的岩性圈闭有效性参数值分布图;
再通过所述目标地层的盖层厚度有效界限值、封堵层厚度有效界限值、曲率有效界限值计算得到所述目标地层的岩性圈闭有效性阈值;
判断所述岩性圈闭有效性参数值与所述岩性圈闭有效性阈值的大小,得到所述目标地层的有效圈闭发育区分布图。
进一步地,所述纵波阻抗数值IMP通过以下公式计算得到:
IMP=1000000/AC×DEN
其中:IMP为纵波阻抗数值,g/cm3*m/s;AC为纵波声波时差,μs/m;DEN为补偿密度,g/cm3。
进一步地,建立纵波阻抗数值与排驱压力的关系式、纵波阻抗数值与突破压力的关系式,筛选出相关性系数最高的关系式作为封堵性表征关系式,包括以下步骤:
获取不同所述地层深度的所述纵波阻抗数值和对应的排驱压力、获取不同所述地层深度的所述纵波阻抗数值和对应的突破压力;
分别采用不同的函数拟合方法拟合所述纵波阻抗数值和所述排驱压力、所述纵波阻抗数值和所述突破压力,根据拟合得到的函数关系式的相关性系数筛选出相关性系数最高的关系式,作为封堵性表征关系式;
所述函数拟合方法包括:线性拟合、多项式拟合、指数拟合、对数拟合、幂函数拟合。
进一步地,所述盖层分为泥页岩层、膏岩层和致密灰岩层;所述总盖层厚度为所述泥页岩层、膏岩层、致密灰岩层的厚度之和;所述目标地层内不同岩性的盖层分布根据所述目标地层的录井岩屑、钻井岩芯资料获取。
进一步地,所述目标地层的纵波阻抗三维数据体通过所述三维地震数据体开展纵波阻抗反演得到;所述三维地震数据体根据所述已有钻井的单井合成地震记录、地震采集数据得到。
进一步地,目标地层内不同地层深度对应的突破压力和排驱压力通过以下方法获取:沿所述目标地层纵向采集若干岩心样品,并测定每一岩心样品的突破压力和排驱压力,确定目标地层内不同地层深度对应的突破压力和排驱压力。
进一步地,所述岩心样品的取样间隔为≤0.2m,保证每一岩性的地层至少有一个岩心样品。
进一步地,测定每一岩心样品的突破压力和排驱压力之后,筛选出突破压力大于等于5MPa、排驱压力大于等于5MPa的岩心样品作为合格的岩心样品,并将合格的岩心样品对应的突破压力和排驱压力用于建立纵波阻抗数值与排驱压力的关系式和纵波阻抗数值与突破压力的关系式。通过排驱压力和突破压力的数据筛选,去掉不合格的排驱压力和突破压力数据值和不合格的岩心样品,筛选出有效数据,提高岩性圈闭有效性识别的准确性。
与现有技术相比,本发明的具有以下有益效果:本发明公开的岩性圈闭有效性识别的方法,建立了裂缝表征参数、盖层厚度和封堵层厚度与岩性圈闭识别的关系式,通过这一关系式来进行识别岩性圈闭发育区情况,实现了在建立岩性圈闭有效性参数计算式过程中引入用于表征所述目标地层裂缝和断裂分布特征的曲率,弥补了现有技术中岩性圈闭有效性识别时未考虑地层裂缝和断裂影响的问题,提高了岩性圈闭有效性识别的准确性,为油气勘探提供了更可靠的参考资料;此外,上述技术方案中,还通过所述目标地层的三维地震数据体结合封堵性表征关系式来计算所述目标地层的封堵层分布、裂缝和断裂分布,减小了由于所述目标地层内的已有钻井数量较少使获得的突破压力和排驱压力数据量较少而导致的应用范围局限,解决了已有钻井数量较少、测井数据较缺乏的地层区域内难以开展岩性圈闭有效性评价工作的问题。
附图说明
图1本发明公开的岩性圈闭有效性识别的方法的流程示意图;
图2本发明一些实施例中纵波阻抗数值IMP和所述突破压力Pb关系式拟合示意图;
图3本发明一些实施例中目标地层中的盖层厚度分布图;
图4本发明一些实施例中目标地层的封堵层厚度分布图;
图5本发明一些实施例中目标地层的曲率分布图;
图6本发明一些实施例中目标地层的有效圈闭发育区分布图。
具体实施方式
下面结合试验例及具体实施方式对本发明作进一步的详细描述。但不应将此理解为本发明上述主题的范围仅限于以下的实施例,凡基于本发明内容所实现的技术均属于本发明的范围。
传统的岩性圈闭有效性识别方法通常是采用目标地层内已有钻井的突破压力和排驱压力来进行评价,而忽略了地层裂缝和断层对岩性圈闭有效性的影响,使得到的有效岩性圈闭与实际情况有较大的出入,影响油气勘探。同时,由于有些区域内的已有钻井数量较少或者分布不均匀,缺少测井数据,导致测定的突破压力和排驱压力难以表征这一区域内的有效岩性圈闭情况,导致采用突破压力和排驱压力进行岩性圈闭有效性识别这一方法的适用性及较差。
基于上述技术问题的考虑,本申请公开了一种岩性圈闭有效性识别的方法,该方法采用曲率来表征地层裂缝和断层,在岩性圈闭有效性参数计算式中引入曲率,并通过三维地震数据体与纵波阻抗三维数据体来获得封堵层厚度、盖层厚度,获得岩性圈闭有效性参数计算式,进而判断目标地层的有效圈闭发育区分布图。相较于传统的岩性圈闭有效性识别方法,本申请公开的岩性圈闭有效性识别方法引入了表征目标地层裂缝和断裂的曲率,在进行岩性圈闭有效性识别过程中考虑了地层裂缝和断裂的影响因素,提高了岩性圈闭识别的准确性;同时,在计算盖层厚度和封堵层厚度的过程中加入了所述目标地层的地震数据,弥补了所述目标地层内可能存在的已有钻井数量较少、分布不均匀导致的岩性圈闭识别准确性差、无法开展的问题。
本申请公开的一种岩性圈闭有效性识别的方法,如图1所示,包括以下步骤:
S1:获取目标地层中盖层的封堵性表征关系式
获取目标地层内不同地层深度对应的突破压力Pb和排驱压力Pd;
根据所述目标地层内已有钻井的纵波声波时差测井曲线和密度测井曲线,获得不同地层深度的纵波阻抗数值IMP;建立纵波阻抗数值与排驱压力的关系式、纵波阻抗数值与突破压力的关系式,并获得所述目标地层中盖层的封堵性表征关系式;
S2:获取目标地层的盖层厚度分布图
根据所述目标地层中盖层的岩性分布,获得不同岩性的盖层的纵波阻抗数值分布范围;
根据所述目标地层的三维地震数据体得到所述目标地层的纵波阻抗三维数据体,并根据所述目标地层中盖层的所述纵波阻抗数值的分布范围筛选不同岩性的盖层,得到所述目标地层中的盖层厚度分布图;
S3:获得目标地层的封堵层厚度分布图
根据所述封堵性表征关系式和所述三维地震数据体获得封堵性参数三维数据体,并根据所述封堵性参数三维数据体,筛选出合格的封堵层,得到所述目标地层的封堵层厚度分布图;
S4:获得目标地层的曲率分布图
根据所述目标地层的三维地震数据体,计算所述目标地层的曲率有效界限值和所述目标地层的曲率分布图;
S5:建立所述目标地层的岩性圈闭有效性参数计算式,并根据所述目标地层的盖层厚度有效界限值、封堵层厚度有效界限值、曲率有效界限值得到所述目标地层的岩性圈闭有效性阈值,得到所述目标地层的有效圈闭发育区分布图;所述岩性圈闭有效性参数计算式为:
其中,所述f(C)为盖层厚度;所述f(S)为封堵层厚度,f(F)为曲率。
需要说明的是,目标地层内不同地层深度对应的突破压力和排驱压力通过以下方法获取:沿所述目标地层纵向采集若干岩心样品,并测定每一岩心样品的突破压力和排驱压力,确定目标地层内不同地层深度对应的突破压力和排驱压力。所述突破压力和排驱压力的测定通过岩心驱替模拟实验获得,其具体包括以下步骤:将采集到的岩心样品抽真空,再用流体进行充分饱和,并根据被测的岩心样品采集时所在的地层深度设置围压和温度,然后用氮气排驱,定时加压直至气体突破,即可得到这一岩心样品的突破压力和排驱压力;更换不同的岩心样品,即可获得所有采集的岩心样品的突破压力和排驱压力;并记载每一岩心样品被采集的已有钻井的位置、采集深度。所述流体可以为地层水、标准盐水或煤油中的任一种。
所述岩心样品的取样间隔为≤0.2m,保证每一岩性的地层至少有一个岩心样品。
在一些实施例中,测定每一岩心样品的突破压力和排驱压力之后,筛选出突破压力大于等于5MPa、排驱压力大于等于5MPa的岩心样品作为合格的岩心样品,并将合格的岩心样品对应的突破压力和排驱压力用于建立纵波阻抗数值与排驱压力的关系式和纵波阻抗数值与突破压力的关系式。
需要说明的是,所述纵波阻抗数值IMP通过以下公式计算得到:
IMP=1000000/AC×DEN (1)
其中:IMP为纵波阻抗数值,g/cm3*m/s;AC为纵波声波时差,us/m;DEN为补偿密度,g/cm3。
需要说明的是,所述已有钻井的纵波声波时差测井曲线通过补偿声波时差测井仪测得;所述已有钻井的密度测井曲线通过补偿密度测径仪测得。根据上述筛选获得的岩心样品所在地层深度,即可获得不同地层深度的纵波声波时差和补偿密度,并带入式(1)即可获得不同地层深度下的纵波阻抗数值。
在一些实施例中,建立纵波阻抗数值与排驱压力的关系式、纵波阻抗数值与突破压力的关系式,筛选出相关性系数最高的关系式作为封堵性表征关系式,包括以下步骤:
获取不同所述地层深度的所述纵波阻抗数值和对应的排驱压力、获取不同所述地层深度的所述纵波阻抗数值和对应的突破压力;
分别采用不同的函数拟合方法拟合所述纵波阻抗数值和所述排驱压力的关系式、所述纵波阻抗数值和所述突破压力的关系式,根据拟合得到的函数关系式的相关性系数筛选出相关性系数最高的关系式,作为封堵性表征关系式;
所述函数拟合方法包括:线性拟合、多项式拟合、指数拟合、对数拟合、幂函数拟合。
需要说明的是,所述步骤S2中,根据所述目标地层中已有钻井的录井岩屑、钻井岩芯资料即可获得所述目标地层中已有钻井在不同地层深度的岩性情况,并根据王少昌等人在1987年提出的盖层判别标准(王少昌,1987,《陕甘宁盆地上古生界泥质岩封盖能力研究》):即可判断已有钻井中泥页岩层、膏岩层、致密灰岩层所在的地层深度;并根据步骤S1中获得的在不同地层深度的纵波阻抗数值,即可获得不同岩性的盖层的纵波阻抗数值分布范围。
需要说明的是,所述三维地震数据体通过所述目标地层的已有钻井的单井合成地震记录、地震采集记录经过处理后得到。所述三维地震数据体(如《中华人民共和国石油天然气行业标准》中的《地震勘探数据处理成果验收规程》所示)包括若干三维地震数据点,每一三维地震数据点均包括一个三维坐标和这一三维坐标所在的地震合成记录的振幅;所述地震合成记录通过单井合成地震记录和地震采集记录获得。
在一些实施例中,所述目标地层的纵波阻抗三维数据体通过所述三维地震数据体开展纵波阻抗反演得到;所述三维地震数据体根据所述已有钻井的单井合成地震记录、地震采集数据得到。
需要说明的是,由于地层中,每一岩性的盖层可能并不是连续的,因此在计算不同的盖层厚度时,需要先根据不同岩性的盖层对应的纵波阻抗数值分布范围对纵波阻抗三维数据体进行筛选分类,然后将纵波阻抗三维数据体中同一岩性的盖层的厚度进行统计累加,得到这一岩性的盖层厚度;进而得到目标地层的盖层厚度分布特征,进而绘制成目标地层的盖层厚度分布图。
在步骤S3中,由于封堵性表征关系式是从建立纵波阻抗数值与排驱压力的关系式和纵波阻抗数值与突破压力的关系式中筛选出的相关性更好的关系式,即本发明中的目标地层的封堵性表征是通过排驱压力或突破压力中相关性更强的一个来进行表征。通过步骤S1确定的封堵性表征关系式将纵波阻抗三维数据体转化为封堵性参数三维数据体,然后根据突破压力或排驱压力的≥5MPa的筛选出合格的封堵层,然后累加封堵层厚度,进而得到所述目标地层的封堵层厚度分布图。
在一些实施例中,所述步骤S4中,所述曲率有效界限值通过以下方法获得:
根据所述目标地层的断层构造解释得到所述目标地层的断层平面展布特征,然后根据《中华人民共和国石油天然气行业标准》中的《地震勘探数据特殊处理技术规程》获得所述目标地层的地震曲率三维数据体,具体的计算公式如式(4):
式中,k为曲率,x,y,z为对应的三维坐标。
将获得的所述目标地层的地震曲率三维数据体与所述目标地层的断层平面展布特征的进行比对,筛选出所述目标地层的平面上所有断层处的曲率值,并所述目标地层的平面上的所有断层处的曲率值中的最小值作为所述曲率有效界限值。
需要说明的是,在步骤S5中,所述封堵层厚度有效界限值为5m;所述盖层分为泥页岩层、膏岩层和致密灰岩层;当所述盖层为泥页岩层时,所述盖层厚度有效界限值为5m;当所述盖层为膏岩层时,所述盖层有效界限值为5m,当所述盖层为致密灰岩层时,所述盖层有效界限值为100m;所述曲率有效界限值为所述目标地层内断层处曲率的最小值。
以下结合某地区的#112井所在地层为例对申请进行说明,其中,X1~X12均为该区域内的已有钻井,其具体实施过程如下:
S1:在该地区的已有钻井进行岩心样品取样,岩心样品的取样间隔为0.2m,并测定每一岩心样品的突破压力和排驱压力,得到不同地层深度对应的突破压力和排驱压力测试数据;
并根据封堵性评价标准,挑选排驱压力≥5MPa或突破压力≥5Mpa的岩心样品和对应的突破压力和排驱压力测试数据;
根据每一岩心样品对应的地层深度寻找对应的纵波声波时差测井曲线和密度测井曲线,根据纵波阻抗计算公式(式(1))计算出该深度的纵波阻抗数值,如表1所示;
表1 某地区#112井实验结果统计表
获取不同所述地层深度的所述纵波阻抗数值IMP和对应的排驱压力Pd、获取不同所述地层深度的所述纵波阻抗数值IMP和对应的突破压力Pb;
如图2所示,分别采用不同的函数拟合方法拟合所述纵波阻抗数值IMP和所述排驱压力Pd、所述纵波阻抗数值IMP和所述突破压力Pb,筛选出相关性系数最高的关系式;突破压力Pb与所述纵波阻抗数值IMP的关系式为作为封堵性表征关系式,突破压力Pb作为表征封堵性的参数。拟合得到的纵波阻抗数值IMP与突破压力Pb的关系式为:
Pb=8×10-08IMP2-0.0055IMP+86.902 (3)
S2:根据该区域内已有钻井的录井岩屑、钻井岩心资料,可获得所述目标地层中已有钻井在不同地层深度的岩性情况,并根据王少昌等人在1987年提出的盖层判别标准:即可判断已有钻井中泥页岩层、膏岩层、致密灰岩层所在的地层深度;并根据步骤S1中获得的在不同地层深度的纵波阻抗数值,即可获得不同岩性的盖层的纵波阻抗数值分布范围:致密灰岩盖层纵波阻抗数值范围为15500~17500g/cm3×m/s,泥岩盖层纵波阻抗数值范围为9690~13030g/cm3×m/s;
根据该区域内已有钻井的单井合成地震记录、地震采集记录获得这一区域内的三维地震数据体,并根据三维地震数据体和单井合成地震记录开展纵波阻抗反演,得到该区域内的纵波阻抗三维数据体,并根据所述目标地层中盖层的所述纵波阻抗数值的分布范围在纵波阻抗三维数据体中画出盖层厚度分布图,如图3所示,厚度大于100m的致密灰岩、泥岩盖层在全区均有分布。
S3:根据纵波阻抗数值IMP与突破压力Pb的关系式式(3),将纵波阻抗三维数据体转换为封堵性参数三维数据体,然后根据突破压力或排驱压力的≥5MPa的筛选出合格的封堵层,然后累加封堵层厚度,进而得到所述目标地层的封堵层厚度分布图,如图4所示;
S4:利用三维地震数据进行地震曲率属性的计算,得到工区内断裂及裂缝分布特征,如图5所示;
S6:建立该区域内的岩性圈闭有效性参数计算式,并根据所述目标地层的盖层厚度有效界限值、封堵层厚度有效界限值、曲率有效界限值得到该区域内的岩性圈闭有效性阈值,得到该区域内的有效圈闭发育区分布图;所述岩性圈闭有效性参数计算式为:
其中,所述f(C)为盖层厚度;所述f(S)为封堵层厚度,f(F)为曲率。
根据该区域的已有钻井实际情况,盖层厚度有效界限值为100m、封堵层厚度有效界限值为5m、曲率有效界限值为0.04,则根据式(2)计算得到该区域内的岩性圈闭有效性阈值为12500;利用岩性圈闭有效性阈值,对岩性圈闭的有效性进行综合分析与评价,明确有效圈闭发育区,如图6所示。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种岩性圈闭有效性识别的方法,其特征在于,包括以下步骤:
获取目标地层内不同地层深度对应的突破压力和排驱压力;
根据所述目标地层内已有钻井的纵波声波时差测井曲线和密度测井曲线,获得不同地层深度的纵波阻抗数值;建立纵波阻抗数值与排驱压力的关系式、纵波阻抗数值与突破压力的关系式,并获得所述目标地层中盖层的封堵性表征关系式;
根据所述目标地层中盖层的岩性分布,获得不同岩性的盖层的纵波阻抗数值分布范围;
根据所述目标地层的三维地震数据体得到所述目标地层的纵波阻抗三维数据体,并根据所述目标地层中盖层的所述纵波阻抗数值的分布范围筛选不同岩性的盖层,得到所述目标地层中的盖层厚度分布图;
根据所述封堵性表征关系式和所述三维地震数据体获得封堵性参数三维数据体,并根据所述封堵性参数三维数据体,筛选出合格的封堵层,得到所述目标地层的封堵层厚度有效界限值和所述目标地层的封堵层厚度分布图;
根据所述目标地层的三维地震数据体,计算所述目标地层的曲率有效界限值和所述目标地层的曲率分布图;
建立所述目标地层的岩性圈闭有效性参数计算式,并根据所述目标地层的盖层厚度有效界限值、封堵层厚度有效界限值、曲率有效界限值得到所述目标地层的岩性圈闭有效性阈值,得到所述目标地层的有效圈闭发育区分布图;
所述岩性圈闭有效性参数计算式为:
其中,f(C)为盖层厚度;f(S)为封堵层厚度,f(F)为曲率。
2.根据权利要求1所述的岩性圈闭有效性识别的方法,其特征在于,所述封堵层厚度有效界限值为5m;
所述盖层分为泥页岩层、膏岩层和致密灰岩层;当所述盖层为泥页岩层时,所述盖层厚度有效界限值为5m;当所述盖层为膏岩层时,所述盖层有效界限值为5m,当所述盖层为致密灰岩层时,所述盖层有效界限值为100m;
所述曲率有效界限值为所述目标地层内断层处曲率的最小值。
3.根据权利要求1所述的岩性圈闭有效性识别的方法,其特征在于,所述目标地层的有效圈闭发育区分布图获得步骤如下:
根据所述岩性圈闭有效性参数计算式,并带入所述目标地层的所述盖层厚度、所述封堵层厚度、所述曲率,获得所述目标地层内的岩性圈闭有效性参数值分布图;
再通过所述目标地层的盖层厚度有效界限值、封堵层厚度有效界限值、曲率有效界限值计算得到所述目标地层的岩性圈闭有效性阈值;
判断所述岩性圈闭有效性参数值与所述岩性圈闭有效性阈值的大小,得到所述目标地层的有效圈闭发育区分布图。
4.根据权利要求1所述的岩性圈闭有效性识别的方法,其特征在于,所述纵波阻抗数值通过以下公式计算得到:
IMP=1000000/AC×DEN
其中:IMP为纵波阻抗数值,g/cm3*m/s;AC为纵波声波时差,μs/m;DEN为补偿密度,g/cm3。
5.根据权利要求1或2所述的岩性圈闭有效性识别的方法,其特征在于,建立纵波阻抗数值与排驱压力的关系式、纵波阻抗数值与突破压力的关系式,筛选出相关性系数最高的关系式作为封堵性表征关系式,包括以下步骤:
获取不同所述地层深度的所述纵波阻抗数值和对应的排驱压力、获取不同所述地层深度的所述纵波阻抗数值IMP和对应的突破压力;
分别采用不同的函数拟合方法拟合所述纵波阻抗数值和所述排驱压力、所述纵波阻抗数值和所述突破压力,根据拟合得到的函数关系式的相关性系数筛选出相关性系数最高的关系式,作为封堵性表征关系式;
所述函数拟合方法包括:线性拟合、多项式拟合、指数拟合、对数拟合、幂函数拟合。
6.根据权利要求2所述的岩性圈闭有效性识别的方法,其特征在于,所述盖层厚度为所述泥页岩层、膏岩层、致密灰岩层的厚度之和;所述目标地层内不同岩性的盖层分布根据所述目标地层的录井岩屑、钻井岩芯资料获取。
7.根据权利要求1所述的岩性圈闭有效性识别的方法,其特征在于,所述目标地层的纵波阻抗三维数据体通过所述三维地震数据体开展纵波阻抗反演得到;所述三维地震数据体根据所述已有钻井的单井合成地震记录、地震采集数据得到。
8.根据权利要求1所述的岩性圈闭有效性识别的方法,其特征在于,目标地层内不同地层深度对应的突破压力和排驱压力通过以下方法获取:沿所述目标地层纵向采集若干岩心样品,并测定每一岩心样品的突破压力和排驱压力,确定目标地层内不同地层深度对应的突破压力和排驱压力。
9.根据权利要求8所述的岩性圈闭有效性识别的方法,其特征在于,所述岩心样品的取样间隔为≤0.2m,保证每一岩性的地层至少有一个岩心样品。
10.根据权利要求9所述的岩性圈闭有效性识别的方法,其特征在于,测定每一岩心样品的突破压力和排驱压力之后,筛选出突破压力大于或等于5MPa且排驱压力大于或等于5MPa的岩心样品作为合格的岩心样品,并将合格的岩心样品对应的突破压力和排驱压力用于建立纵波阻抗数值与排驱压力的关系式、纵波阻抗数值与突破压力的关系式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010953383.XA CN114167515B (zh) | 2020-09-11 | 2020-09-11 | 一种岩性圈闭有效性识别的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010953383.XA CN114167515B (zh) | 2020-09-11 | 2020-09-11 | 一种岩性圈闭有效性识别的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114167515A true CN114167515A (zh) | 2022-03-11 |
CN114167515B CN114167515B (zh) | 2023-03-28 |
Family
ID=80475944
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010953383.XA Active CN114167515B (zh) | 2020-09-11 | 2020-09-11 | 一种岩性圈闭有效性识别的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114167515B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115793053A (zh) * | 2022-12-22 | 2023-03-14 | 中国地质调查局油气资源调查中心 | 一种基于地震、钻井与录井的寒武系膏盐岩综合预测方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110320414A1 (en) * | 2010-06-28 | 2011-12-29 | Nhn Corporation | Method, system and computer-readable storage medium for detecting trap of web-based perpetual calendar and building retrieval database using the same |
CN103116186A (zh) * | 2011-11-16 | 2013-05-22 | 中国石油天然气集团公司 | 一种小尺度非均质储集体的容积确定方法 |
CN104637075A (zh) * | 2013-11-13 | 2015-05-20 | 中国石油化工股份有限公司 | 砂体自动快速成图方法 |
RU2672766C1 (ru) * | 2018-02-08 | 2018-11-19 | Общество с ограниченной ответственностью "Газпромнефть Научно-Технический Центр" (ООО "Газпромнефть НТЦ") | СПОСОБ ПРОГНОЗИРОВАНИЯ МОРФОМЕТРИЧЕСКИХ ПАРАМЕТРОВ РУСЛОВЫХ ТЕЛ (палеоканалов) |
CN109061765A (zh) * | 2018-09-26 | 2018-12-21 | 西南石油大学 | 非均质薄砂岩互层油藏的圈闭评价方法 |
US20190011584A1 (en) * | 2017-07-05 | 2019-01-10 | Petrochina Company Limited | Method and apparatus for identifying low permeable conglomerate diagenetic trap |
US20190080273A1 (en) * | 2017-09-08 | 2019-03-14 | Petrochina Company Limited | Method and apparatus for evaluating exploitation value of a geological resource |
US20190100997A1 (en) * | 2017-09-30 | 2019-04-04 | Petrochina Company Limited | Oil and gas zone effectiveness evaluation method and apparatus |
CN110414783A (zh) * | 2016-08-16 | 2019-11-05 | 中国石油化工股份有限公司 | 一种圈闭评价方法 |
CN110954944A (zh) * | 2019-10-31 | 2020-04-03 | 中国石油大港油田勘探开发研究院 | 一种断层圈闭含油高度地震预测方法 |
-
2020
- 2020-09-11 CN CN202010953383.XA patent/CN114167515B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110320414A1 (en) * | 2010-06-28 | 2011-12-29 | Nhn Corporation | Method, system and computer-readable storage medium for detecting trap of web-based perpetual calendar and building retrieval database using the same |
CN103116186A (zh) * | 2011-11-16 | 2013-05-22 | 中国石油天然气集团公司 | 一种小尺度非均质储集体的容积确定方法 |
CN104637075A (zh) * | 2013-11-13 | 2015-05-20 | 中国石油化工股份有限公司 | 砂体自动快速成图方法 |
CN110414783A (zh) * | 2016-08-16 | 2019-11-05 | 中国石油化工股份有限公司 | 一种圈闭评价方法 |
US20190011584A1 (en) * | 2017-07-05 | 2019-01-10 | Petrochina Company Limited | Method and apparatus for identifying low permeable conglomerate diagenetic trap |
US20190080273A1 (en) * | 2017-09-08 | 2019-03-14 | Petrochina Company Limited | Method and apparatus for evaluating exploitation value of a geological resource |
US20190100997A1 (en) * | 2017-09-30 | 2019-04-04 | Petrochina Company Limited | Oil and gas zone effectiveness evaluation method and apparatus |
RU2672766C1 (ru) * | 2018-02-08 | 2018-11-19 | Общество с ограниченной ответственностью "Газпромнефть Научно-Технический Центр" (ООО "Газпромнефть НТЦ") | СПОСОБ ПРОГНОЗИРОВАНИЯ МОРФОМЕТРИЧЕСКИХ ПАРАМЕТРОВ РУСЛОВЫХ ТЕЛ (палеоканалов) |
CN109061765A (zh) * | 2018-09-26 | 2018-12-21 | 西南石油大学 | 非均质薄砂岩互层油藏的圈闭评价方法 |
CN110954944A (zh) * | 2019-10-31 | 2020-04-03 | 中国石油大港油田勘探开发研究院 | 一种断层圈闭含油高度地震预测方法 |
Non-Patent Citations (1)
Title |
---|
杨占龙等: "利用地震信息评价陆相盆地岩性圈闭的关键点分析", 《岩性油气藏》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115793053A (zh) * | 2022-12-22 | 2023-03-14 | 中国地质调查局油气资源调查中心 | 一种基于地震、钻井与录井的寒武系膏盐岩综合预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114167515B (zh) | 2023-03-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3414566B1 (en) | Thermal maturity determination of rock formations using mud gas isotope logging | |
CN108374657B (zh) | 井断点自动识别方法 | |
KR101853334B1 (ko) | 순차층서 분석을 통한 비전통자원 개발 구간 예측 방법 | |
CN110954944A (zh) | 一种断层圈闭含油高度地震预测方法 | |
CN110320569B (zh) | 一种致密砂岩储层单井裂缝发育强度定量评价方法 | |
CN110927794A (zh) | 一种致密储层裂缝识别与孔隙度定量计算的方法 | |
CN112145165B (zh) | 一种微裂缝-孔隙型储层动静态渗透率转换方法 | |
CN114215513A (zh) | 潜山油藏模式的定量判别方法、装置、介质及设备 | |
CN112835124B (zh) | 基于成像测井和阵列声波测井资料的裂缝有效性评价方法 | |
CN115130934A (zh) | 基于zel模型和多源数据的区域井漏风险预评价方法 | |
CN114167515B (zh) | 一种岩性圈闭有效性识别的方法 | |
CN107861162A (zh) | 基于微电极测井资料的天然裂缝识别方法及系统 | |
CN110231407B (zh) | 一种判断碳酸盐岩盖层有效性的方法 | |
CN112746838A (zh) | 发育天然裂缝储层的水平井不同井段岩石可压性判定方法 | |
CN107831540A (zh) | 储层物性参数直接提取新方法 | |
CN112946782B (zh) | 一种致密油气储渗体地震精细刻画方法 | |
CN111381292A (zh) | 一种预测砂岩含烃储层的测井解释方法与装置 | |
CN109826623B (zh) | 一种致密砂岩储层层理缝的地球物理测井判识方法 | |
CN112198551A (zh) | 一种通用的定量化评价储层裂缝强度的新方法 | |
CN113050168B (zh) | 基于阵列声波测井和声波远探测测井资料的裂缝有效性评价方法 | |
CN113075748B (zh) | 基于成像测井和声波远探测测井资料的裂缝有效性评价方法 | |
CN115391739A (zh) | 一种裂缝渗透率的定量计算方法和系统 | |
CN113325468B (zh) | 储层缝洞分布范围预测方法及装置 | |
CN113792932A (zh) | 一种利用微震-损伤-渗流关系的页岩气产量预测方法 | |
CN106353813A (zh) | 基于阵列声波测井的流体性质识别方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |