CN112796738A - 一种阵列声波测井和常规测井结合的地层渗透率计算方法 - Google Patents
一种阵列声波测井和常规测井结合的地层渗透率计算方法 Download PDFInfo
- Publication number
- CN112796738A CN112796738A CN202110154244.5A CN202110154244A CN112796738A CN 112796738 A CN112796738 A CN 112796738A CN 202110154244 A CN202110154244 A CN 202110154244A CN 112796738 A CN112796738 A CN 112796738A
- Authority
- CN
- China
- Prior art keywords
- permeability
- core
- stratum
- model
- time difference
- 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
- 230000035699 permeability Effects 0.000 title claims abstract description 283
- 238000004364 calculation method Methods 0.000 title claims abstract description 108
- 238000000034 method Methods 0.000 claims abstract description 17
- 239000011435 rock Substances 0.000 claims abstract description 17
- 239000011148 porous material Substances 0.000 claims abstract description 15
- 238000010586 diagram Methods 0.000 claims description 36
- 230000015572 biosynthetic process Effects 0.000 claims description 21
- 238000005755 formation reaction Methods 0.000 claims description 21
- 238000011160 research Methods 0.000 claims description 21
- 238000004458 analytical method Methods 0.000 claims description 12
- 230000000704 physical effect Effects 0.000 claims description 12
- 238000000611 regression analysis Methods 0.000 claims description 5
- 238000007796 conventional method Methods 0.000 abstract description 6
- 238000011161 development Methods 0.000 abstract description 5
- 238000013461 design Methods 0.000 abstract description 3
- 238000005481 NMR spectroscopy Methods 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/12—Means for transmitting measuring-signals or control signals from the well to the surface, or from the surface to the well, e.g. for logging while drilling
- E21B47/14—Means for transmitting measuring-signals or control signals from the well to the surface, or from the surface to the well, e.g. for logging while drilling using acoustic waves
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Theoretical Computer Science (AREA)
- Geochemistry & Mineralogy (AREA)
- General Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Geophysics (AREA)
- Business, Economics & Management (AREA)
- Environmental & Geological Engineering (AREA)
- Acoustics & Sound (AREA)
- Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Computer Hardware Design (AREA)
- Remote Sensing (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种阵列声波测井和常规测井结合的地层渗透率计算方法,包括(1)数据准备;(2)岩心孔渗关系分类;(3)渗透率敏感参数优选;(4)渗透率计算模型建立;(5)渗透率模型优选;(6)渗透率计算与误差计算。本发明不同于常规的利用孔隙度计算渗透率的方法,本方法主要结合阵列声波测井和常规测井,考虑多方面影响因素,并针对不同孔渗情况的地层设计相应的渗透率计算模型,提高地层渗透率计算精度;利用其计算结果可以在油气勘探阶段定量分析地层的渗透性,帮助后续油气藏的勘探与开发。
Description
技术领域
本发明属于油气勘探开发领域,具体涉及一种阵列声波测井和常规测井结合的地层渗透率计算方法。
背景技术
现有的地层渗透率计算模型主要是基于微观参数和常规测井曲线建立的,包括与孔隙度建立相关关系,也大量使用了孔隙或者颗粒微观参数(孔隙半径、颗粒半径、比表面积、弯曲度、束缚水含量、流动单元等),广泛运用了核磁共振测井(T2、核磁共振孔隙度)和电阻率测井(孔隙度、电阻率),其次是自然伽马测井。
阵列声波测井是具有多探头声系、测量多波列的声波测井方法,通过记录多条曲线进行相关和叠加处理,可以有效地压制干扰,准确提取纵波、横波和斯通利波的各种信息。目前少有基于阵列声波测井的渗透率计算方法,仅有利用横波时差或者纵波时差计算孔隙度,从而计算地层渗透率,其精度远不足以对地层渗透率进行准确定量描述。以前只能利用阵列声波测井定性评价地层渗透性,通过进一步的研究目前已经开始利用阵列声波测井定量确定地层渗透性。
通过调研大量的文献、分析相关测井及实验资料,发现地层渗透率与阵列声波之间具有较强的相关性,本发明在总结前人研究得出的常规渗透率计算公式及结论的基础上,提出了不同于以往的渗透率计算方法——结合了阵列声波测井和常规测井资料的渗透率计算方法,即通过已有的常规测井、阵列声波测井和岩心物性分析资料建立渗透率与各影响因素之间的关系模型,应用到不同孔渗关系的地层,极大地提高了地层渗透率计算精度。利用本发明方法可以在油气勘探阶段定量分析地层的渗透性,进而帮助后续油气藏的勘探与开发。
发明内容
本发明主要是克服现有技术中的缺点,本发明提供一种阵列声波测井和常规测井结合的地层渗透率计算方法,本发明通过综合分析常规测井、阵列声波测井和岩心物性分析资料,研究渗透率与常规测井、阵列声波测井数据之间的相关性,针对不同孔隙度的地层建立不同的高精度渗透率计算模型,从而提高地层渗透率的计算精确率。
本发明解决上述技术问题,所提供的技术方案是:一种阵列声波测井和常规测井结合的地层渗透率计算方法,根据研究区的岩心物性分析资料、常规测井资料和阵列声波测井资料获取岩心孔隙度、岩心渗透率、自然伽马曲线、泥值含量曲线、测井计算孔隙度、横波时差、声波时差和斯通利波时差;
根据研究区的岩心物性分析资料,并利用岩心渗透率和测井计算孔隙度绘制孔隙度-岩心渗透率散点图;
根据孔隙度-岩心渗透率散点图中数据点分布的不同特征,将研究区的岩心地层划分成不同孔渗类别的地层;
根据每类地层的岩心渗透率,每类地层分别建立波时差、孔隙度、自然伽马、泥质含量、横波时差和斯通利波时差与岩心渗透率的拟合关系以及回归系数;并根据回归系数在每类地层中优选出三个渗透率模型敏感参数;
根据每类地层中优选出的三个渗透率模型敏感参数分别建立指数乘积形式渗透率计算模型、幂形式渗透率计算模型;
根据每类地层建立的指数乘积形式渗透率计算模型、幂形式渗透率计算模型分别计算出结果,将每类地层的两个计算结果分别与每类地层的岩心渗透率进行回归分析,比较回归系数的大小,选择回归系数大的模型作为每类地层的最优渗透率计算模型;
根据每类地层的最优渗透率计算模型确定研究区结合阵列声波和常规测井的地层渗透率计算模型;
根据研究区结合阵列声波和常规测井的地层渗透率计算模型以及研究区其他任意没有岩心物性分析资料的地层计算其孔隙度,应用相应的渗透率计算模型得到地层渗透率。
进一步的技术方案是,根据孔隙度-岩心渗透率散点图中数据点分布的不同特征,将研究区的岩心地层划分成不同孔渗类别的地层包括:
根据孔隙度-岩心渗透率散点图图中分布的不同特征确定岩心孔渗分类结果和分类标准;
根据岩心孔渗分类结果和分类标准将研究区的岩心地层划分成不同孔渗类别的地层。
进一步的技术方案是,根据每类地层的岩心渗透率,每类地层分别建立波时差、孔隙度、自然伽马、泥质含量、横波时差和斯通利波时差与岩心渗透率的拟合关系以及回归系数;并根据回归系数每类地层优选出三个渗透率模型敏感参数包括:
根据每类地层的岩心渗透率,每类地层分别绘制岩心渗透率-声波时差散点图、绘制岩心渗透率-孔隙度散点图、岩心渗透率-自然伽马散点图、绘制岩心渗透率-泥质含量散点图、绘制岩心渗透率-横波时差散点图、绘制岩心渗透率-斯通利波时差散点图;并分别拟合计算得到岩心渗透率-声波时差回归系数、岩心渗透率-孔隙度回归系、岩心渗透率-自然伽马回归系数、岩心渗透率-泥质含量回归系数、岩心渗透率-横波时差回归系数、岩心渗透率-斯通利波时差回归系数;
根据岩心渗透率-声波时差回归系数、岩心渗透率-孔隙度回归系,比较这两个回归系数大小,选择回归系数更大的因素作为渗透率模型敏感参数之一;
根据岩心渗透率-自然伽马回归系数、岩心渗透率-泥质含量回归系数,比较这两个回归系数大小,选择回归系数更大的因素作为渗透率模型敏感参数之一;
根据岩心渗透率-横波时差回归系数、岩心渗透率-斯通利波时差回归系数,比较这两个回归系数大小,选择回归系数更大的因素作为渗透率模型敏感参数之一。
进一步的技术方案是,所述指数乘积形式渗透率计算模型为:
式中:K计算为计算渗透率;X1、X2、X3分别为三个渗透率模型敏感参数的值;a1、b1、c1、d1为拟合系数。
进一步的技术方案是,所述幂形式渗透率计算模型为:
式中:K计算为计算渗透率;X1、X2、X3分别为三个渗透率模型敏感参数的值;a2、b2、c2、d2为拟合系数。
进一步的技术方案是,所述指数乘积形式渗透率计算模型的建立过程包括:
首先计算岩心渗透率的对数值和三个渗透率模型敏感参数的对数值,然后建立这四者的多元回归方程;
lgK岩心=a1+b1×lgX1+c1×lgX2+d1×lgX3
式中:K岩心为岩心渗透率;X1、X2、X3分别为三个渗透率模型敏感参数的值;a1、b1、c1、d1为拟合系数;
将上述多元回归方程的等号左右两边同时取作10的指数,得到指数乘积形式的渗透率模型;
式中:K计算为计算渗透率;X1、X2、X3分别为三个渗透率模型敏感参数的值;a1、b1、c1、d1为拟合系数。
进一步的技术方案是,所述幂形式渗透率计算模型的建立过程为:
首先计算岩心渗透率对数值,然后建立岩心渗透率对数值与三个渗透率模型敏感参数的多元回归方程;
lgK岩心=a2+b2×X1+c2×X2+d2×X3
式中:K岩心为岩心渗透率;X1、X2、X3分别为三个渗透率模型敏感参数的值;a2、b2、c2、d2为拟合系数;
将上述多元回归方程的等号左右两边同时取作10的指数,得到幂形式的渗透率模型如下;
式中:K计算为计算渗透率;X1、X2、X3分别为三个渗透率模型敏感参数的值;a2、b2、c2、d2为拟合系数。
本发明具有以下有益效果:本发明不同于常规的利用孔隙度计算渗透率的方法,本方法主要结合阵列声波测井和常规测井,考虑多方面影响因素,并针对不同孔渗情况的地层设计相应的渗透率计算模型,提高地层渗透率计算精度;利用其计算结果可以在油气勘探阶段定量分析地层的渗透性,帮助后续油气藏的勘探与开发。
附图说明
图1为本发明的流程框图;
图2为岩心孔渗关系分类结果图;
图3为实施例1的岩心孔渗关系分类结果图;
图4为实施例1的岩心渗透率与声波时差、孔隙度的拟合关系对比图;
图5为实施例1的岩心渗透率与自然伽马、泥质含量的拟合关系对比图;
图6为实施例1的岩心渗透率与横波时差、斯通利波的拟合关系对比图;
图7为高孔渗地层下不同模型的渗透率计算结果与岩心渗透率回归结果对比图
图8为低孔渗地层下不同模型的渗透率计算结果与岩心渗透率回归分析对比图;
图9为低孔渗地层下本发明方法(结合阵列声波测井和常规测井)与常规方法(利用孔隙度)计算渗透率结果对比图;
图10为高孔渗地层下本发明方法(结合阵列声波测井和常规测井)与常规方法(利用孔隙度)计算渗透率结果对比图。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
此外,下面所描述的本发明不同实施方式中所涉及的技术特征只要彼此之间未构成冲突就可以相互结合。
如图1所示,本发明的一种阵列声波测井和常规测井结合的地层渗透率计算方法,包括以下步骤:
(1)数据准备:准备研究区的岩心物性分析资料、常规测井资料和阵列声波测井资料,根据上述资料获得岩心孔隙度、岩心渗透率、自然伽马曲线、泥值含量曲线、测井计算孔隙度、横波时差和斯通利波时差等参数;
(2)岩心孔渗关系分类:选取步骤(1)中有岩心物性分析资料的地层,利用岩心渗透率和测井计算的孔隙度,绘制孔隙度-岩心渗透率散点图,根据数据点分布的不同特征(如两部分不同孔隙度区域的孔渗拟合曲线有明显的差别),将其划分成不同的孔渗类别,如图2所示,得到岩心孔渗分类结果并确定分类标准:低孔渗区域孔隙度≤A%,A%<中孔渗区域孔隙度≤B%,B%<高孔渗区域孔隙度(其中A<B);其中图2中仅划分了三个类别孔渗关系,实际应用中还可以按照具体情况,划分更多或更少的孔渗关系类别;
(3)渗透率模型敏感参数优选:基于步骤(2)的孔渗关系分类结果,选择其中具有某一类孔渗关系的地层,分别建立声波时差、孔隙度、自然伽马、泥质含量、横波时差和斯通利波时差与岩心渗透率的拟合关系和回归系数,并优选出三个渗透率模型敏感参数:
具体的是:绘制岩心渗透率-声波时差散点图,拟合计算得到岩心渗透率-声波时差回归系数,绘制岩心渗透率-孔隙度散点图,拟合计算得到岩心渗透率-孔隙度回归系数,比较这两个回归系数大小,选择回归系数更大的因素作为渗透率模型敏感参数之一;
绘制岩心渗透率-自然伽马散点图,拟合计算得到岩心渗透率-自然伽马回归系数,绘制岩心渗透率-泥质含量散点图,拟合计算得到岩心渗透率-泥质含量回归系数,比较这两个回归系数大小,选择回归系数更大的因素作为渗透率模型敏感参数之一;
绘制岩心渗透率-横波时差散点图,拟合计算得到岩心渗透率-横波时差回归系数,绘制岩心渗透率-斯通利波时差散点图,拟合计算得到岩心渗透率-斯通利波时差回归系数,比较这两个回归系数大小,选择回归系数更大的因素作为渗透率模型敏感参数之一;
(4)渗透率模型建立:针对步骤(3)中选取的某一类孔渗关系的地层,利用步骤(3)所选择的三个渗透率模型敏感参数建立该类地层的渗透率计算模型;
本发明采用了两种形式的渗透率计算模型:指数乘积形式渗透率计算模型和幂形式渗透率计算模型,以下分别说明建立指数乘积形式渗透率计算模型和幂形式渗透率计算模型的操作步骤:
①指数乘积形式渗透率计算模型
建立指数乘积形式渗透率计算模型的操作步骤如下:
首先计算岩心渗透率的对数值和步骤(3)中所确定的三个渗透率模型敏感参数的对数值,然后建立这四者的多元回归方程(1):
lgK岩心=a1+b1×lgX1+c1×lgX2+d1×lgX3 (1)
式中:K岩心为岩心渗透率;X1、X2、X3分别为三个渗透率模型敏感参数的值;a1、b1、c1、d1为拟合系数;
将公式(1)等号左右两边同时取作10的指数,得到指数乘积形式渗透率计算模型(公式2):
式中:K计算为计算渗透率;X1、X2、X3分别为三个渗透率模型敏感参数的值;a1、b1、c1、d1为拟合系数;
②幂形式渗透率计算模型
建立幂形式渗透率模型的操作步骤如下:
首先计算岩心渗透率对数值,然后建立岩心渗透率对数值与步骤(3)中所确定的三个渗透率模型敏感参数的多元回归方程(3):
lgK岩心=a2+b2×X1+c2×X2+d2×X3 (3)
式中:K岩心为岩心渗透率;X1、X2、X3分别为三个渗透率模型敏感参数的值;a2、b2、c2、d2为拟合系数;
将公式(3)等号左右两边同时取作10的指数,得到幂形式渗透率计算模型(4)如下:
式中:K计算为计算渗透率;X1、X2、X3分别为三个渗透率模型敏感参数的值;a2、b2、c2、d2为拟合系数;
(5)渗透率模型优选:利用步骤(4)建立的两个渗透率计算模型(公式2和公式4),分别计算出某一类孔渗关系地层指数乘积形式和幂形式渗透率模型的计算结果,将计算结果分别与岩心渗透率进行回归分析,比较回归系数的大小,选择回归系数更大的模型作为该类孔渗关系地层的最优渗透率计算模型;
至此,就完成了步骤(2)中某一类孔渗关系地层的最优渗透率模型建立,只需重复步骤(3)到步骤(5)的操作,直至所有类别的孔渗关系均建立出最优的渗透率模型,这样就得到了结合阵列声波和常规测井的地层渗透率计算模型;
(6)渗透率与误差计算:选择研究区其他任意没有岩心物性分析资料的地层计算其孔隙度,根据步骤(5)得到的不同孔渗关系的渗透率计算模型,应用相应的渗透率模型计算得到地层渗透率;
误差计算:利用绝对误差平均值计算公式(5),比较本发明模型的计算精度与常规模型的计算精度;
实施例1
本发明的一种阵列声波测井和常规测井结合的地层渗透率计算方法,具体包括以下步骤:
(1)数据准备:准备研究区的岩心物性分析资料、常规测井资料和阵列声波测井资料,根据资料获取岩心孔隙度、岩心渗透率、自然伽马曲线、泥值含量曲线、测井计算孔隙度、横波时差和斯通利波时差;
(2)岩心孔渗关系分类:选取步骤(1)中有岩心物性分析资料的所有地层,利用岩心渗透率和测井计算的孔隙度,绘制孔隙度-岩心渗透率散点图,根据其分布的不同特征(如两部分不同孔隙度区域的孔渗拟合曲线有明显的差别),将其划分成不同的孔渗类别,如图3所示,确定分类标准:低孔渗区域孔隙度≤9%,9%<高孔渗区域孔隙度;
(3)渗透率模型敏感参数优选:基于步骤(2)的孔渗关系分类结果,选择高孔渗区域孔隙度的地层;
(31)通过绘制岩心渗透率-声波时差散点图(如图4(a)所示),拟合计算得到岩心渗透率-声波时差回归系数RAC=0.433,绘制岩心渗透率-孔隙度散点图(如图4(b)所示),拟合计算得到岩心渗透率-孔隙度的回归系数Rφ=0.462,显然Rφ>RAC,即孔隙度与岩心渗透率的相关性比声波时差更好,故选择孔隙度作为本实施例的渗透率模型敏感参数之一;
(32)通过绘制岩心渗透率-自然伽马散点图(如图5(a)所示),拟合计算得到岩心渗透率-自然伽马回归系数RGR=0.178,绘制岩心渗透率-泥质含量散点图(如图5(b)所示),拟合计算得到岩心渗透率-泥质含量回归系数RSH=0.182,显然RSH>RGR,即泥质含量与岩心渗透率的相关性比自然伽马更好,故选择泥质含量作为本实施例的渗透率模型敏感参数之一;
(33)通过绘制岩心渗透率-横波时差散点图(如图6(a)所示),拟合计算得到岩心渗透率-横波时差回归系数RDTD=0.395,绘制岩心渗透率-斯通利波时差散点图(如图6(b)所示),拟合计算得到岩心渗透率-斯通利波时差回归系数RDTST=0.627,显然RDTST>RDTD,即斯通利波时差与岩心渗透率的相关性比横波时差更好,故选择斯通利波时差作为本实施例的渗透率模型敏感参数之一;
(34)综上确定出本实施例的高孔渗类型地层的渗透率计算模型的三个渗透率模型敏感参数包括孔隙度、泥质含量和斯通利波时差;
(4)渗透率模型建立:针对步骤(3)中选取的高孔渗关系的地层,利用步骤(3)所选择的三个渗透率模型敏感参数建立渗透率计算模型;指数乘积形式和幂形式,以下分别说明建立指数乘积形式模型和幂形式模型的操作步骤:
①指数形式
建立指数乘积形式的渗透率模型的操作步骤如下:
首先计算岩心渗透率对数值、孔隙度对数值、斯通利波时差对数值和泥质含量对数值,然后建立这四者的多元回归方程(6):
lgK岩心=-174.502+7.826×lgφ+73.085×lgDTST+2.268×lgSH (6)
式中:K岩心为岩心渗透率;φ为孔隙度;DTST为斯通利波时差;SH为泥质含量;
最后得到指数乘积形式的渗透率计算模型(7)如下:
K计算=10-174.502×φ7.826×DTST73.085×SH-2.268 (7)
式中:K计算为计算渗透率;φ为孔隙度;DTST为斯通利波时差;SH为泥质含量;
②幂形式
建立幂形式的渗透率模型的操作步骤如下:
首先计算岩心渗透率对数值、孔隙度、斯通利波时差和泥质含量,然后建立这四者的多元回归方程(8):
lgK岩心=-34.461+0.229φ+0.16DTST+0.054SH (8)
式中:K岩心为岩心渗透率;φ为孔隙度;DTST为斯通利波时差;SH为泥质含量;
最后得到幂形式的渗透率计算模型(9)如下:
K计算=10-34.461+0.229φ+0.16DTST+0.054SH (9)
式中:K计算为计算渗透率;φ为孔隙度;DTST为斯通利波时差;SH为泥质含量;
(5)渗透率模型优选:利用步骤(4)中建立的两个渗透率计算模型(公式7和公式9),计算得到所选的高孔渗关系地层的指数乘积形式和幂形式两类模型的渗透率计算结果,将两者分别与所选地层的岩心渗透率进行回归分析,如图7所示,比较得到的岩心渗透率与两个模型计算渗透率的回归系数,故选择线性回归系数更大的指数乘积形式的计算模型(公式7)作为本实施例的高孔渗关系的最优渗透率计算模型;
至此,就完成了步骤(2)所划分的多个孔渗关系的其中一个——高孔渗关系地层的最优渗透率计算模型的建立,只需重复步骤(3)到步骤(5)的操作得到低孔渗地层的两个渗透率模型,如图8所示,通过比较两个模型的回归系数选出低孔渗地层的最优渗透率模型——幂形式的渗透率计算模型,于是就得到了所有孔渗关系地层的渗透率模型如表1所示:
表1不同孔渗关系的渗透率计算模型
孔渗关系类别 | 渗透率计算模型 |
低孔渗(φ<9%) | K<sub>计算</sub>=10<sup>11.041+0.713φ-0.068DTST+0.0097SH</sup> |
高孔渗(9%≤φ) | K<sub>计算</sub>=10<sup>-174.502</sup>×φ<sup>7.826</sup>×DTST<sup>73.085</sup>×SH<sup>-2.268</sup> |
(6)渗透率与误差计算
选择研究区其他任意没有岩心物性分析资料的地层计算其孔隙度,根据步骤(5)得到的不同孔渗关系的渗透率计算模型,应用相应的渗透率模型计算得到地层渗透率,即只需根据所选地层的具体孔隙度情况应用所建立的渗透率计算模型,就可以得到该目标地层高精度的渗透率计算结果。
如图9、图10所示,为不同孔渗关系地层下本发明方法(结合阵列声波测井和常规测井)与常规方法(利用孔隙度)计算渗透率结果对比;
误差计算:利用绝对误差平均值计算公式(5),如表2所示,得到本发明的模型计算渗透率的相对误差与常规模型计算渗透率的相对误差。
表2在不同孔渗情况下不同模型计算渗透率的相对误差对比
本发明设计模型 | 常规渗透率计算模型 | |
高孔渗地层绝对误差 | 0.36 | 0.47 |
低孔渗地层绝对误差 | 0.23 | 0.58 |
结合表1所展示的在不同孔渗情况下不同模型计算渗透率的平均绝对误差(利用公式5计算平均绝对误差)对比结果,以及图9和图10所示的本发明方法(结合阵列声波测井和常规测井数据)与常规方法(利用孔隙度)计算渗透率的结果对比,可以发现利用本发明所建模型计算的地层渗透率结果在精度上有进一步的提高,利用其计算结果可以在油气勘探阶段定量分析地层的渗透性,帮助后续油气藏的勘探与开发。
以上所述,并非对本发明作任何形式上的限制,虽然本发明已通过上述实施例揭示,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,可利用上述揭示的技术内容作出些变动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (7)
1.一种阵列声波测井和常规测井结合的地层渗透率计算方法,其特征在于,包括以下步骤:
根据研究区的岩心物性分析资料、常规测井资料和阵列声波测井资料获取岩心孔隙度、岩心渗透率、自然伽马曲线、泥值含量曲线、测井计算孔隙度、横波时差、声波时差和斯通利波时差;
根据研究区的岩心渗透率和测井计算孔隙度绘制孔隙度-岩心渗透率散点图;
根据孔隙度-岩心渗透率散点图中数据点分布的不同特征,将研究区的地层划分成不同的孔渗类别;
根据每类地层的岩心渗透率,每类地层分别建立声波时差、孔隙度、自然伽马、泥质含量、横波时差和斯通利波时差与岩心渗透率的拟合关系以及回归系数;并根据回归系数为每类地层优选出三个渗透率模型敏感参数;
根据每类地层优选出的三个渗透率模型分别建立指数乘积形式渗透率计算模型、幂形式渗透率计算模型;
根据每类地层建立的指数乘积形式渗透率计算模型、幂形式渗透率计算模型分别计算出结果,将每类地层的两个计算结果分别与每类地层的岩心渗透率进行回归分析,比较回归系数的大小,选择回归系数大的模型作为每类地层的最优渗透率计算模型;
根据每类地层的最优渗透率计算模型确定研究区结合阵列声波和常规测井的地层渗透率计算模型;
根据研究区结合阵列声波和常规测井的地层渗透率计算模型以及研究区其他任意没有岩心物性分析资料的地层计算其孔隙度,应用相应的渗透率计算模型得到地层渗透率。
2.根据权利要求1所述的一种阵列声波测井和常规测井结合的地层渗透率计算方法,其特征在于,根据孔隙度-岩心渗透率散点图图中分布的不同特征,将研究区的岩心地层划分成不同孔渗类别的地层包括:
根据孔隙度-岩心渗透率散点图中数据点分布的不同特征确定岩心孔渗分类结果和分类标准;
根据岩心孔渗分类结果和分类标准将研究区的岩心地层划分成不同孔渗类别的地层。
3.根据权利要求1所述的一种阵列声波测井和常规测井结合的地层渗透率计算方法,其特征在于,根据每类地层的岩心渗透率,每类地层分别建立声波时差、孔隙度、自然伽马、泥质含量、横波时差和斯通利波时差与岩心渗透率的拟合关系以及回归系数;并根据回归系数每类地层优选出三个渗透率模型敏感参数包括:
根据每类地层的岩心渗透率,每类地层分别绘制岩心渗透率-声波时差散点图、绘制岩心渗透率-孔隙度散点图、岩心渗透率-自然伽马散点图、绘制岩心渗透率-泥质含量散点图、绘制岩心渗透率-横波时差散点图、绘制岩心渗透率-斯通利波时差散点图;并分别拟合计算得到岩心渗透率-声波时差回归系数、岩心渗透率-孔隙度回归系、岩心渗透率-自然伽马回归系数、岩心渗透率-泥质含量回归系数、岩心渗透率-横波时差回归系数、岩心渗透率-斯通利波时差回归系数;
根据岩心渗透率-声波时差回归系数、岩心渗透率-孔隙度回归系,比较这两个回归系数大小,选择回归系数更大的因素作为渗透率模型敏感参数之一;
根据岩心渗透率-自然伽马回归系数、岩心渗透率-泥质含量回归系数,比较这两个回归系数大小,选择回归系数更大的因素作为渗透率模型敏感参数之一;
根据岩心渗透率-横波时差回归系数、岩心渗透率-斯通利波时差回归系数,比较这两个回归系数大小,选择回归系数更大的因素作为渗透率模型敏感参数之一。
6.根据权利要求4所述的一种阵列声波测井和常规测井结合的地层渗透率计算方法,其特征在于,所述指数乘积形式渗透率计算模型的建立过程包括:
首先计算岩心渗透率的对数值和三个渗透率模型敏感参数的对数值,然后建立这四者的多元回归方程;
lgK岩心=a1+b1×lgX1+c1×lgX2+d1×lgX3
式中:K岩心为岩心渗透率;X1、X2、X3分别为三个渗透率模型敏感参数的值;a1、b1、c1、d1为拟合系数;
将上述多元回归方程的等号左右两边同时取作10的指数,得到指数乘积形式的渗透率模型;
式中:K计算为计算渗透率;X1、X2、X3分别为三个渗透率模型敏感参数的值;a1、b1、c1、d1为拟合系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110154244.5A CN112796738A (zh) | 2021-02-04 | 2021-02-04 | 一种阵列声波测井和常规测井结合的地层渗透率计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110154244.5A CN112796738A (zh) | 2021-02-04 | 2021-02-04 | 一种阵列声波测井和常规测井结合的地层渗透率计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112796738A true CN112796738A (zh) | 2021-05-14 |
Family
ID=75814162
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110154244.5A Pending CN112796738A (zh) | 2021-02-04 | 2021-02-04 | 一种阵列声波测井和常规测井结合的地层渗透率计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112796738A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113202459A (zh) * | 2021-06-11 | 2021-08-03 | 中国矿业大学 | 一种利用钻孔自然伽马测井预测断层渗透性的方法 |
CN113836747A (zh) * | 2021-10-18 | 2021-12-24 | 大连理工大学 | 一种水合物渗透率关联分析和模型预测的实时可视化方法 |
CN114458310A (zh) * | 2022-02-16 | 2022-05-10 | 西南石油大学 | 一种断层低速率运动条件下原油侧向分流定向评价的方法 |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2612515A1 (en) * | 2005-06-24 | 2007-01-04 | Exxonmobil Upstream Research Company | Method for determining reservoir permeability from borehole stoneley-wave attenuation using biot's poroelastic theory |
CN102400679A (zh) * | 2010-09-10 | 2012-04-04 | 北京锐致聪科技有限公司 | 一种用于多级阵列声波测井中的切向极化单极子换能器 |
US20120318500A1 (en) * | 2011-06-15 | 2012-12-20 | Esg Solutions Inc. | Methods and systems for monitoring and modeling hydraulic fracturing of a reservoir field |
CN104374685A (zh) * | 2014-11-26 | 2015-02-25 | 中国石油天然气集团公司 | 一种数字岩心三孔隙组分求取储层渗透率的方法 |
CN104500049A (zh) * | 2014-10-20 | 2015-04-08 | 成都创源油气技术开发有限公司 | 页岩气地球物理快速评价方法 |
CN104948176A (zh) * | 2015-05-08 | 2015-09-30 | 西南石油大学 | 一种基于渗透增大率识别碳酸盐岩储层裂缝的方法 |
CN106154351A (zh) * | 2016-08-09 | 2016-11-23 | 中国石油天然气集团公司 | 一种低孔渗储层渗透率的估算方法 |
CN107622139A (zh) * | 2016-07-15 | 2018-01-23 | 中国石油天然气股份有限公司 | 裂缝渗透率的计算方法 |
GB201816802D0 (en) * | 2017-08-15 | 2018-11-28 | Petrochina Co Ltd | Method and apparatus for determining permeability of reservoirs |
CN109753755A (zh) * | 2019-01-25 | 2019-05-14 | 中国石油天然气集团有限公司 | 一种确定储层含水饱和度的方法 |
CN110348135A (zh) * | 2019-07-15 | 2019-10-18 | 中国石油大学(华东) | 一种随钻声波测井评价地层渗透率的方法 |
WO2020048028A1 (zh) * | 2018-09-03 | 2020-03-12 | 中国石油大学(华东) | 一种基于压裂潜力的待压裂水平井压裂设计方法和装置 |
CN110955982A (zh) * | 2019-12-18 | 2020-04-03 | 长江大学 | 一种变质岩储层渗透率计算方法、装置及计算机存储介质 |
CN110984970A (zh) * | 2019-10-09 | 2020-04-10 | 中国海洋石油集团有限公司 | 一种利用地层测试确定启动压力梯度的方法 |
CN111827995A (zh) * | 2020-07-31 | 2020-10-27 | 中国石油天然气集团有限公司 | 一种基于核磁共振的超压砾岩储层渗透率的计算方法 |
-
2021
- 2021-02-04 CN CN202110154244.5A patent/CN112796738A/zh active Pending
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2612515A1 (en) * | 2005-06-24 | 2007-01-04 | Exxonmobil Upstream Research Company | Method for determining reservoir permeability from borehole stoneley-wave attenuation using biot's poroelastic theory |
CN102400679A (zh) * | 2010-09-10 | 2012-04-04 | 北京锐致聪科技有限公司 | 一种用于多级阵列声波测井中的切向极化单极子换能器 |
US20120318500A1 (en) * | 2011-06-15 | 2012-12-20 | Esg Solutions Inc. | Methods and systems for monitoring and modeling hydraulic fracturing of a reservoir field |
CN104500049A (zh) * | 2014-10-20 | 2015-04-08 | 成都创源油气技术开发有限公司 | 页岩气地球物理快速评价方法 |
CN104374685A (zh) * | 2014-11-26 | 2015-02-25 | 中国石油天然气集团公司 | 一种数字岩心三孔隙组分求取储层渗透率的方法 |
CN104948176A (zh) * | 2015-05-08 | 2015-09-30 | 西南石油大学 | 一种基于渗透增大率识别碳酸盐岩储层裂缝的方法 |
CN107622139A (zh) * | 2016-07-15 | 2018-01-23 | 中国石油天然气股份有限公司 | 裂缝渗透率的计算方法 |
CN106154351A (zh) * | 2016-08-09 | 2016-11-23 | 中国石油天然气集团公司 | 一种低孔渗储层渗透率的估算方法 |
GB201816802D0 (en) * | 2017-08-15 | 2018-11-28 | Petrochina Co Ltd | Method and apparatus for determining permeability of reservoirs |
WO2020048028A1 (zh) * | 2018-09-03 | 2020-03-12 | 中国石油大学(华东) | 一种基于压裂潜力的待压裂水平井压裂设计方法和装置 |
CN109753755A (zh) * | 2019-01-25 | 2019-05-14 | 中国石油天然气集团有限公司 | 一种确定储层含水饱和度的方法 |
CN110348135A (zh) * | 2019-07-15 | 2019-10-18 | 中国石油大学(华东) | 一种随钻声波测井评价地层渗透率的方法 |
CN110984970A (zh) * | 2019-10-09 | 2020-04-10 | 中国海洋石油集团有限公司 | 一种利用地层测试确定启动压力梯度的方法 |
CN110955982A (zh) * | 2019-12-18 | 2020-04-03 | 长江大学 | 一种变质岩储层渗透率计算方法、装置及计算机存储介质 |
CN111827995A (zh) * | 2020-07-31 | 2020-10-27 | 中国石油天然气集团有限公司 | 一种基于核磁共振的超压砾岩储层渗透率的计算方法 |
Non-Patent Citations (8)
Title |
---|
FENG WU: "Predicting gasewater relative permeability using Nuclear Magnetic Resonance and Mercury Injection Capillary Pressure measurements", 《JOURNAL OF NATURAL GAS SCIENCE AND ENGINEERING》 * |
何广军主编: "《现代测试技术原理与应用》", 31 December 2012, 国防工业出版社 * |
张建欣: "裂缝性储层关键参数测井计算方法", 《西安科技大学学报》 * |
张蕾等: "应用阵列声波测井资料定量评价碳酸盐岩储层渗透性", 《测井技术》 * |
张鹏浩: "低渗碳酸盐岩储层核磁共振渗透率计算新方法", 《当代化工》 * |
李玮: "利用等效孔隙度法估算火成岩储层渗透率", 《测井技术》 * |
程梦薇: "基于储层分类的渗透率测井计算方法", 《长江大学学报(自科版)》 * |
范宜仁等: "基于核磁共振刻度流动单元复杂砂岩储层渗透率建模方法", 《测井技术》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113202459A (zh) * | 2021-06-11 | 2021-08-03 | 中国矿业大学 | 一种利用钻孔自然伽马测井预测断层渗透性的方法 |
CN113836747A (zh) * | 2021-10-18 | 2021-12-24 | 大连理工大学 | 一种水合物渗透率关联分析和模型预测的实时可视化方法 |
CN113836747B (zh) * | 2021-10-18 | 2024-02-23 | 大连理工大学 | 一种水合物渗透率关联分析和模型预测的实时可视化方法 |
CN114458310A (zh) * | 2022-02-16 | 2022-05-10 | 西南石油大学 | 一种断层低速率运动条件下原油侧向分流定向评价的方法 |
CN114458310B (zh) * | 2022-02-16 | 2024-05-24 | 西南石油大学 | 一种断层低速率运动条件下原油侧向分流定向评价的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112796738A (zh) | 一种阵列声波测井和常规测井结合的地层渗透率计算方法 | |
CN111425193B (zh) | 一种基于聚类分析测井岩石物理相划分的储层可压性评价方法 | |
CN103775075A (zh) | 一种全井段岩性识别方法 | |
CN104047598A (zh) | 非均质古岩溶碳酸盐岩储层产能预测方法 | |
CN111665560A (zh) | 油气层识别方法、装置、计算机设备及可读存储介质 | |
CN112698399B (zh) | 一种砂砾岩井测震联动约束高效储层定量预测方法与*** | |
CN111598440A (zh) | 多角度驱动的复杂介质储层渗透率定量评价方法及*** | |
CN108982320A (zh) | 一种利用粒度参数进行复杂孔隙结构储层渗透率计算方法 | |
CN112253087A (zh) | 一种基于多源测井资料的生物扰动储集层物性计算方法 | |
CN107762483B (zh) | 一种基于测井曲线的相关系数及包络面积的流体识别方法 | |
CN114114414A (zh) | 一种页岩储层“甜点”信息人工智能预测方法 | |
CN113216945B (zh) | 一种致密砂岩储层渗透性定量评价方法 | |
CN113552624B (zh) | 孔隙度预测方法及装置 | |
CN109655394B (zh) | 一种孔喉结构参数约束下的核磁共振t2谱渗透率计算方法 | |
CN111984928A (zh) | 一种测井资料计算页岩油储层有机碳含量的方法 | |
CN113219531A (zh) | 致密砂岩气水分布的识别方法及装置 | |
CN111812716A (zh) | 页岩气储层总有机碳含量叠前定量预测方法、装置及设备 | |
CN115012903B (zh) | 判别泥页岩层理构造发育的测井评价方法 | |
CN113065286B (zh) | 一种基于常规测井的致密砂砾岩储层孔隙结构评价方法 | |
CN112505154B (zh) | 泥页岩储层矿物成分含量解析与岩相识别表征方法 | |
CN114565116A (zh) | 一种低渗储层物性参数的预测方法和装置 | |
CN112541523B (zh) | 一种泥质含量计算方法和装置 | |
CN116930023A (zh) | 致密砂岩相控分类孔隙度测井精细解释方法和装置 | |
CN112346147A (zh) | 一种基于中子密度孔隙度差的储层评价方法 | |
CN113050191A (zh) | 基于双参数的页岩油toc预测方法及装置 |
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 |