CN112145165B - 一种微裂缝-孔隙型储层动静态渗透率转换方法 - Google Patents
一种微裂缝-孔隙型储层动静态渗透率转换方法 Download PDFInfo
- Publication number
- CN112145165B CN112145165B CN201911363854.5A CN201911363854A CN112145165B CN 112145165 B CN112145165 B CN 112145165B CN 201911363854 A CN201911363854 A CN 201911363854A CN 112145165 B CN112145165 B CN 112145165B
- Authority
- CN
- China
- Prior art keywords
- reservoir
- microcrack
- permeability
- pore
- static
- 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
- 230000035699 permeability Effects 0.000 title claims abstract description 165
- 239000011148 porous material Substances 0.000 title claims abstract description 128
- 230000003068 static effect Effects 0.000 title claims abstract description 115
- 238000006243 chemical reaction Methods 0.000 title claims abstract description 29
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000012545 processing Methods 0.000 claims abstract description 11
- 238000003384 imaging method Methods 0.000 claims description 32
- 238000004364 calculation method Methods 0.000 claims description 20
- 208000013201 Stress fracture Diseases 0.000 claims description 11
- 238000013145 classification model Methods 0.000 claims description 10
- 239000004215 Carbon black (E152) Substances 0.000 claims description 8
- 229930195733 hydrocarbon Natural products 0.000 claims description 8
- 150000002430 hydrocarbons Chemical class 0.000 claims description 8
- 230000004044 response Effects 0.000 claims description 8
- 238000011161 development Methods 0.000 claims description 7
- 230000000704 physical effect Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 2
- 238000012360 testing method Methods 0.000 abstract description 13
- 230000008901 benefit Effects 0.000 abstract description 3
- 230000015572 biosynthetic process Effects 0.000 description 12
- 239000012530 fluid Substances 0.000 description 4
- 229910052925 anhydrite Inorganic materials 0.000 description 3
- OSGAYBCDTDRGGQ-UHFFFAOYSA-L calcium sulfate Chemical compound [Ca+2].[O-]S([O-])(=O)=O OSGAYBCDTDRGGQ-UHFFFAOYSA-L 0.000 description 3
- 229910000514 dolomite Inorganic materials 0.000 description 3
- 239000010459 dolomite Substances 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 238000007635 classification algorithm Methods 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 238000003066 decision tree Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 210000002615 epidermis Anatomy 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP 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
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Abstract
本发明提供了一种微裂缝‑孔隙型储层动静态渗透率转换方法，其包括以下内容：1)确定所发育微裂缝的类型和特征；2)通过测井资料处理和解释得到微裂缝孔隙度；3)建立每一类储层中有效孔隙度和静态渗透率的函数关系，得到全井段的静态渗透率；4)基于试井或测试作业得到每口井的动态渗透率；5)以微裂缝孔隙度的大小为分类标准，区分储层类型；6)分别建立静态渗透率与动态渗透率之间的转换关系；7)基于这一转换关系，在无法测得动态渗透率的井中，计算得到动态渗透率。本发明方法一种微裂缝‑孔隙型储层动静态渗透率转换方法，提供了一种通用性强、经济性好的有效方法，可以广泛应用于微裂缝‑孔隙型储层动静态渗透率转换中。
Description
技术领域
本发明涉及一种储层动静态渗透率转换方法，尤其涉及一种微裂缝-孔隙型储层动静态渗透率转换方法。
背景技术
渗透率是指在一定压差下，岩石允许流体通过的能力。静态渗透率主要是根据岩石物理实验的分析数据，或测井资料解释结果得到，主要反映流体在静止状态下时储层的渗流能力。一般情况下，静态渗透率的计算主要是基于孔隙度和渗透率之间高精度的函数关系，采用孔隙度去建立静态渗透率的计算公式。动态渗透率主要是根据试井、测试等资料计算得到，主要反映流体在流动状态下时岩石的渗流能力。一般情况下，动态渗透率的计算主要是依据达西定律计算得到。
在均质储层，静态渗透率与动态渗透率之间的差异并不大，也存在较好的对应关系。因此，在均质储层，静态渗透率可以准确表征油气藏的产能。而在微裂缝-孔隙型储层，由于储层的强非均质性，静态渗透率与动态渗透率存在较大差异，且对应关系较差，致使静态渗透率难以准确表征油气藏的产能。仅通过现有技术中的对于微裂缝-孔隙型储层的静态渗透率的计算无法得知微裂缝-孔隙型储层的动态渗透率。
同时，由于油气产能直接影响勘探、开发效果和经济效益，因此，针对微裂缝-孔隙型储层形成一种动静态渗透率转换方法，能有效地提高微裂缝-孔隙型油气藏产能评价的精度，并节省高昂的测试成本是十分必要的。
发明内容
针对上述问题，本发明的目的是提供一种微裂缝-孔隙型储层静态渗透率与动态渗透率之间的转换方法。
为实现上述目的，本发明拟采取以下技术方案：
一种微裂缝-孔隙型储层动、静态渗透率转换方法，包括以下步骤：
1)基于地层某一深度下微裂缝-孔隙型油气藏的岩心样品，对岩心柱子和岩心薄片中所发育的微裂缝进行综合分析，得出所发育的微裂缝的数据信息；
2)基于微裂缝-孔隙型油气藏的目的层段的电成像或阵列声波测井资料和步骤1中的微裂缝的数据信息建立不同类型微裂缝及其特征在电成像或阵列声波测井资料上的发育模式；
3)分析微裂缝在电成像或阵列声波测井资料上的响应特征，并通过测井资料处理和解释得到微裂缝孔隙度；
4)建立每一类储层中有效孔隙度和静态渗透率的函数关系；
5)建立基于不同测井曲线数值划分不同微裂缝-孔隙型储层类别的分类模型；
6)分别建立微裂缝-孔隙型储层、孔隙型储层静态渗透率与动态渗透率之间的转换关系；
7)基于步骤6)中的转换关系，采用静态渗透率便可以计算得到动态渗透率。
其中，所述步骤1)中的微裂缝的相关数据包括：微裂缝的类型、微裂缝的密度、微裂缝的宽度及微裂缝延伸长度；所述微裂缝的类型包括：高角度缝、斜交缝、诱导缝、充填-半充填缝和斜交网状缝。
其中，所述步骤2)中具体包括利用测井仪器在微裂缝-孔隙型油气藏的目的层段采集电成像或阵列声波测井资料；基于所采集的电成像或阵列声波测井资料，分析步骤1)中岩心柱子和岩心薄片上所发育微裂缝的类型、密度、宽度和延伸长度在电成像或阵列声波测井资料上的响应特征，并建立不同类型微裂缝及其特征在电成像或阵列声波测井资料上的发育模式。
其中，所述步骤3)中：通过对微裂缝-孔隙型油气藏目的层段所采集的电成像或阵列声波测井资料的处理与解释，得到微裂缝-孔隙型油气藏目的层段每一个测量深度点的微裂缝孔隙度；
基于岩心样品的物性分析结果，采用岩心分析的有效孔隙度、静态渗透率计算微裂缝-孔隙型油气藏目的层段不同深度储层的储集空间指标RSI，计算公式为：
式中，Kc为岩心分析的静态渗透率，mD；为岩心分析的有效孔隙度，小数。
其中，所述步骤4)中，每一类微裂缝-孔隙型储层有效孔隙度与静态渗透率之间的函数关系如下：
式中，Ks为静态渗透率；为有效孔隙度；g、h为常数。
其中，所述步骤4)中将步骤3)计算得到的不同深度储层的储集空间指标RSI细分为多个分布范围，并使得每一储集空间指标RSI的范围区间对应着一类微裂缝-孔隙型储层，分类方式如下：
Ι＝f(a＜RSI≤b),ΙΙ＝f(c＜RSI≤d),ΙΙΙ＝f(e＜RSI≤f)···
式中，Ι、ΙΙ、ΙΙΙ为微裂缝-孔隙型储层所属的类别号；RSI为储集空间指标；a、b、c、d、e、f为储集空间指标数值；
当每个分布范围内的有效孔隙度与静态渗透率之间的函数关系均满足所对应的一类微裂缝-孔隙型储层静态渗透率计算的精度时，则以该细分范围为准；如若不满足所对应的一类微裂缝-孔隙型储层静态渗透率计算的精度时，则需重新对不同深度储层的储集空间指标RSI细分，直到细分后每个范围的有效孔隙度与静态渗透率之间的函数关系均满足所对应的一类微裂缝-孔隙型储层静态渗透率计算的精度为止。其中，每一类微裂缝-孔隙型储层有效孔隙度与静态渗透率之间的函数关系如下：
式中，Ks为静态渗透率；为有效孔隙度；g、h为常数。
其中，步骤4)中，当每个分布范围内的有效孔隙度与静态渗透率之间的函数关系能够达到所对应的一类微裂缝-孔隙型储层渗透率精度的80％，即认定有效孔隙度与静态渗透率之间的函数关系均满足所对应的一类微裂缝-孔隙型储层静态渗透率计算的精度。
其中，所述步骤5)中分类模型可通过Geolog、GeoFrame、Forward、Lead等软件进行建立。
其中，还包括静态渗透率的获取：基于不同测井曲线数值的组合，采用微裂缝-孔隙型储层的分类模型去划分一口井整个井段的类别。针对每个井段微裂缝-孔隙型储层所属的类别，采用相应的有效孔隙度与静态渗透率之间的函数关系来计算该类别中微裂缝-孔隙型储层的静态渗透率。
其中，所述步骤6)中以微裂缝孔隙度的大小为分类标准，当微裂缝孔隙度大于某一数值时，井中所发育的储层为微裂缝-孔隙型储层；当微裂缝孔隙度小于某一数值时，井中所发育的储层则为孔隙型储层；然后分别建立微裂缝-孔隙型储层、孔隙型储层静态渗透率与动态渗透率之间的转换关系。
其中，所述步骤6)中的微裂缝-孔隙型储层、孔隙型储层静态渗透率与动态渗透率之间的转换关系为：
Kd＝mKs+n
式中，Kd为动态渗透率；Ks为静态渗透率；m、n为常数。
本发明由于采取以上技术方案，其具有以下优点：1、本发明的微裂缝-孔隙型储层动静态渗透率转换方法，能准确地匹配微裂缝-孔隙型储层静态渗透率与动态渗透率，大幅提高产能评价的精度2、本发明方法避免开展大量的试井、测试作业，能够有效地节约成本，具有较强的经济性。3、本发明为微裂缝-孔隙型储层动静态渗透率的转换，提供了一种通用性强、经济性好的有效方法，并可以广泛应用于微裂缝-孔隙型储层动静态渗透率的转换中。
附图说明
图1是微裂缝-孔隙型储层岩心柱子的微裂缝分布图；
图2是微裂缝-孔隙型储层岩心镜下的微裂缝分布图；
图3是某微裂缝-孔隙型油藏A井部分地层段微裂缝在电成像和阵列声波测井上的响应特征；
图4是某微裂缝-孔隙型油藏A井部分地层段基于电成像测井资料处理与解释的微裂缝孔隙度成果图；
图5是某微裂缝-孔隙型油藏B井部分地层段静态渗透率计算成果图；
图6是某微裂缝-孔隙型油藏36口井静态渗透率与动态渗透率的关系图；
图7是某微裂缝-孔隙型油藏30口井微裂缝-孔隙型储层静态渗透率与动态渗透率的关系图；
图8是某微裂缝-孔隙型油藏6口井孔隙型储层静态渗透率与动态渗透率的关系图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
本发明涉及一种微裂缝-孔隙型储层动静态渗透率转换方法，其包括以下步骤：
1)基于地层某一深度下微裂缝-孔隙型油气藏的岩心样品，对岩心柱子和岩心薄片中所发育的微裂缝进行综合分析，得出所发育微裂缝的类型(高角度缝、斜交缝、诱导缝、充填-半充填缝、斜交网状缝)、微裂缝的密度、微裂缝的宽度及微裂缝延伸长度等信息。
2)利用测井仪器在微裂缝-孔隙型油气藏的目的层段采集电成像或阵列声波测井资料；基于所采集的电成像或阵列声波测井资料，分析步骤1)中岩心柱子和岩心薄片上所发育微裂缝的类型、密度、宽度和延伸长度在电成像或阵列声波测井资料上的响应特征，并建立不同类型微裂缝及其特征在电成像或阵列声波测井资料上的发育模式。
3)通过对微裂缝-孔隙型油气藏目的层段所采集的电成像或阵列声波测井资料的处理与解释，可以得到微裂缝-孔隙型油气藏目的层段每一个测量深度点的微裂缝孔隙度。
4)基于岩心样品的物性分析结果，采用岩心分析的有效孔隙度、静态渗透率计算微裂缝-孔隙型油气藏目的层段不同深度储层的储集空间指标RSI，计算公式为：
式中，Kc为岩心分析的静态渗透率，mD；为岩心分析的有效孔隙度，小数。
5)将计算得到的不同深度储层的储集空间指标RSI细分为多个分布范围，并使得每一储集空间指标RSI的范围区间对应着一类微裂缝-孔隙型储层，分类方式如下：
Ι＝f(a＜RSI≤b),ΙΙ＝f(c＜RSI≤d),ΙΙΙ＝f(e＜RSI≤f)···
式中，Ι、ΙΙ、ΙΙΙ为微裂缝-孔隙型储层所属的类别号；RSI为储集空间指标；a、b、c、d、e、f为储集空间指标数值。
当每个分布范围内的有效孔隙度与静态渗透率之间的函数关系均满足所对应的一类微裂缝-孔隙型储层静态渗透率计算的精度时，则以该细分范围为准；如若不满足所对应的一类微裂缝-孔隙型储层静态渗透率计算的精度时，则需重新对不同深度储层的储集空间指标RSI细分，直到细分后每个范围的有效孔隙度与静态渗透率之间的函数关系均满足所对应的一类微裂缝-孔隙型储层静态渗透率计算的精度为止。其中，每一类微裂缝-孔隙型储层有效孔隙度与静态渗透率之间的函数关系如下：
式中，Ks为静态渗透率；为有效孔隙度；g、h为常数。
6)针对细分后每一类储集空间指标RSI的范围区间所对应的微裂缝-孔隙型储层的类别，采用分类算法，建立基于不同测井曲线数值划分不同微裂缝-孔隙型储层类别的分类模型。
通过反复调整微裂缝-孔隙型储层的类别和不同测井曲线数值的组合，在保证每一类微裂缝-孔隙型储层有效孔隙度与静态渗透率之间的函数关系满足静态渗透率计算精度的同时，最大程度地提高分类模型的精度。
7)基于不同测井曲线数值的组合，采用微裂缝-孔隙型储层的分类模型去划分一口井整个井段的类别。针对每个井段微裂缝-孔隙型储层所属的类别，采用相应的有效孔隙度与静态渗透率之间的函数关系来计算该类别中微裂缝-孔隙型储层的静态渗透率。利用不同类别微裂缝-孔隙型储层静态渗透率的组合，便完成了整个井段微裂缝-孔隙型储层静态渗透率的定量评价。
8)以微裂缝孔隙度的大小为分类标准，当微裂缝孔隙度大于某一数值时，井中所发育的储层为微裂缝-孔隙型储层；当微裂缝孔隙度小于某一数值时，井中所发育的储层则为孔隙型储层；然后分别建立微裂缝-孔隙型储层、孔隙型储层静态渗透率与动态渗透率之间的转换关系。
9)基于微裂缝-孔隙型储层、孔隙型储层静态渗透率与动态渗透率之间的转换关系，采用步骤7)所得到的静态渗透率，分别可以计算得到微裂缝-孔隙型储层、孔隙型储层的动态渗透率，从而便得到了该口井的动态渗透率。
上述步骤5)中，大多数情况下，当每个分布范围内的有效孔隙度与静态渗透率之间的函数关系能够达到所对应的一类微裂缝-孔隙型储层渗透率精度的80％，即满足要求。
上述步骤6)中，分类模型可通过Geolog、GeoFrame、Forward、Lead等软件进行建立。
下面通过具体实施例对本发明方法做进一步说明。
实施例1：以某微裂缝-孔隙型油藏M为例。
1)基于地层某一深度下微裂缝-孔隙型油藏M的岩心样品，对岩心柱子和岩心薄片中所发育的微裂缝进行综合分析，得出所发育微裂缝的类型主要是高角度缝，裂缝的密度主要分布在2—5条之间，微裂缝的宽度约0.02mm，微裂缝延伸长度约1—2m(图1和图2)。
2)利用测井仪器在微裂缝-孔隙型油藏M的目的层段采集电成像和阵列声波测井资料；基于A井所采集的电成像和阵列声波测井资料，分析步骤1)中岩心柱子和岩心薄片上所发育微裂缝的类型、密度、宽度和延伸长度在电成像和阵列声波测井资料上的响应特征，电成像和阵列声波测井资料均显示微裂缝的存在(图3)。
如图3所示，在A井部分层段测井响应图中，第1道是地层深度；第2道是自然伽马、自然电位、井径测井曲线，表示该地层的岩性特征；第3道是深、中、浅电阻率测井曲线，刻画该地层的电性特征；第4道是体积密度、中子孔隙度和纵波时差测井曲线，反映该地层的物性特征；第5道是电成像测井资料；第6道是阵列声波测井资料；第7道是硬石膏、白云岩和泥岩在地层中的百分含量。
3)结合岩心柱子和岩心薄片上微裂缝分布特征和电成像测井资料上微裂缝的发育特征，主要基于A井所采集的电成像测井资料，处理和解释可以得到A井采集了电成像测井资料地层段的微裂缝孔隙度(图4)。
如图4所示，在A井部分层段测井资料处理与解释成果图中，第1道是地层深度；第2道是自然伽马、自然电位、井径测井曲线，表示该地层的岩性特征；第3道是深、中、浅电阻率测井曲线，刻画该地层的电性特征；第4道是体积密度、中子孔隙度和纵波时差测井曲线，反映该地层的物性特征；第5道是电成像测井资料；第6道是测井解释得到的有效孔隙度；第7道基于电成像测井资料处理和解释得到的微裂缝孔隙度；第8道是硬石膏、白云岩和泥岩在地层中的百分含量。
4)在微裂缝-孔隙型油藏M，基于岩心分析的有效孔隙度和静态渗透率，可以计算得到储集空间指标RSI。对不同深度储层储集空间指标RSI的分布范围进行细分，分别建立每一类储层有效孔隙度和静态渗透率之间的函数关系(表1)。
表1不同类型储层有效孔隙度与静态渗透率之间的函数关系
基于不同测井曲线数值，采用决策树分类算法针对表1中的分类方案建立分类模型；采用分类模型对全井段进行储层类型的划分，将其分为五类。基于表1中每一类储层有效孔隙度和静态渗透率之间的函数关系，可以计算得到每一类储层的静态渗透率。将所计算五类储层的渗透率进行组合，便能得到全井段的静态渗透率(图5)。
如图5所示，在B井部分层段测井资料处理与解释成果图中第1道是地层深度；第2道是自然伽马测井曲线，表示该地层的岩性特征；第3道是深、中、浅电阻率测井曲线，刻画该地层的电性特征；第4道是体积密度、中子孔隙度和纵波时差测井曲线，反映该地层的物性特征；第5道是岩心分析的有效孔隙度、测井计算的有效孔隙度；第6道是岩心分析的静态渗透率、测井计算的静态渗透率；第7道是硬石膏、白云岩和泥岩在地层中的百分含量。
5)基于微裂缝-孔隙型油藏M中36口井的试井或测试作业，采用式(1)中的达西定律可以计算得到36口井的动态渗透率；
式中：Kd为动态渗透率，mD；Q为流量，m3/s；μ为流体粘度，Pa·s；Re为有效供给半径，m；Rw为井眼半径，m；S为表皮因子，无因次；He为有效厚度，m；Pe为边界压力，Pa；Pwf为井底流压，Pa。
另外，可以统计得到36口井射孔层段静态渗透率的平均值(表2)。
表2微裂缝-孔隙型油藏M中36口井的动态渗透率
6)基于表2中微裂缝-孔隙型油藏M中36口井静态渗透率与动态渗透率的关系图(图6)，以微裂缝孔隙度0.005为分类标准，将微裂缝孔隙度大于0.005的井中所发育的储层归为微裂缝-孔隙型储层，将微裂缝孔隙度小于或等于0.005的井中所发育的储层归为孔隙型储层，其中，此油藏中30口井属于微裂缝-孔隙型储层，其静态渗透率与动态渗透率的关系图如图7所示；此油藏中6口井属于孔隙型储层，其静态渗透率与动态渗透率的关系图如图8所示。
7)微裂缝-孔隙型油藏M中，微裂缝-孔隙型储层静态渗透率Ks与动态渗透率Kd之间的函数关系如下：
Kd＝26.273Ks+2.8554
孔隙型储层静态渗透率Ks与动态渗透率Kd之间的函数关系如下：
Kd＝7.9152Ks-64.022。
其中，上述公式中的系数m和n可以基于表2中静态渗透率与动态渗透率线性回归得到。
8)在微裂缝-孔隙型油藏M没有开展试井、测试作业的井中，基于步骤7)所得到微裂缝-孔隙型储层、孔隙型储层静态渗透率Ks与动态渗透率Kd之间的转换关系，采用静态渗透率计算得到微裂缝-孔隙型储层、孔隙型储层的动态渗透率，从而便可以在没有开展试井、测试作业的井中得到动态渗透率，从而为产能评价奠定基础。
上述实施例仅用于说明本发明，其中各部件的结构、连接方式等都是可以有所变化的，凡是在本发明技术方案的基础上进行的等同变换和改进，均不应排除在本发明的保护范围之外。
Claims (8)
1.一种微裂缝-孔隙型储层动静态渗透率转换方法，其特征在于，包括以下步骤：
1)基于地层某一深度下微裂缝-孔隙型油气藏的岩心样品，对岩心柱子和岩心薄片中所发育的微裂缝进行综合分析，得出所发育的微裂缝的数据信息；
2)基于微裂缝-孔隙型油气藏的目的层段的电成像或阵列声波测井资料和步骤1中的微裂缝的数据信息建立不同类型微裂缝及其特征在电成像或阵列声波测井资料上的发育模式；
3)分析微裂缝在电成像或阵列声波测井资料上的响应特征，并通过测井资料处理和解释得到微裂缝孔隙度；
4)建立每一类储层中有效孔隙度和静态渗透率的函数关系；
5)建立基于不同测井曲线数值划分不同微裂缝-孔隙型储层类别的分类模型；
6)分别建立微裂缝-孔隙型储层、孔隙型储层静态渗透率与动态渗透率之间的转换关系；
7)基于步骤6)中的转换关系，采用静态渗透率便可以计算得到动态渗透率；
其中，所述步骤3)中：通过对微裂缝-孔隙型油气藏目的层段所采集的电成像或阵列声波测井资料的处理与解释，得到微裂缝-孔隙型油气藏目的层段每一个测量深度点的微裂缝孔隙度；
基于岩心样品的物性分析结果，采用岩心分析的有效孔隙度、静态渗透率计算微裂缝-孔隙型油气藏目的层段不同深度储层的储集空间指标RSI，计算公式为：
式中，Kc为岩心分析的静态渗透率，mD；为岩心分析的有效孔隙度，小数；
所述步骤4)中，每一类微裂缝-孔隙型储层有效孔隙度与静态渗透率之间的函数关系如下：
式中，Ks为静态渗透率；为有效孔隙度；g、h为常数。
2.根据权利要求1所述的微裂缝-孔隙型储层动静态渗透率转换方法，其特征在于，所述步骤1)中的微裂缝的相关数据包括：微裂缝的类型、微裂缝的密度、微裂缝的宽度及微裂缝延伸长度；所述微裂缝的类型包括：高角度缝、斜交缝、诱导缝、充填-半充填缝和斜交网状缝。
3.根据权利要求2所述的微裂缝-孔隙型储层动静态渗透率转换方法，其特征在于，所述步骤2)中具体包括利用测井仪器在微裂缝-孔隙型油气藏的目的层段采集电成像或阵列声波测井资料；基于所采集的电成像或阵列声波测井资料，分析步骤1)中岩心柱子和岩心薄片上所发育微裂缝的类型、密度、宽度和延伸长度在电成像或阵列声波测井资料上的响应特征，并建立不同类型微裂缝及其特征在电成像或阵列声波测井资料上的发育模式。
4.根据权利要求1所述的微裂缝-孔隙型储层动静态渗透率转换方法，其特征在于，所述步骤4)中将步骤3)计算得到的不同深度储层的储集空间指标RSI细分为多个分布范围，并使得每一储集空间指标RSI的范围区间对应着一类微裂缝-孔隙型储层，分类方式如下：
Ι＝f(a＜RSI≤b),ΙΙ＝f(c＜RSI≤d),ΙΙΙ＝f(e＜RSI≤f)···
式中，Ι、ΙΙ、ΙΙΙ为微裂缝-孔隙型储层所属的类别号；RSI为储集空间指标；a、b、c、d、e、f为储集空间指标数值；
当每个分布范围内的有效孔隙度与静态渗透率之间的函数关系均满足所对应的一类微裂缝-孔隙型储层静态渗透率计算的精度时，则以该细分范围为准；如若不满足所对应的一类微裂缝-孔隙型储层静态渗透率计算的精度时，则需重新对不同深度储层的储集空间指标RSI细分，直到细分后每个范围的有效孔隙度与静态渗透率之间的函数关系均满足所对应的一类微裂缝-孔隙型储层静态渗透率计算的精度为止,其中，每一类微裂缝-孔隙型储层有效孔隙度与静态渗透率之间的函数关系如下：
式中，Ks为静态渗透率；为有效孔隙度；g、h为常数。
5.根据权利要求4所述的微裂缝-孔隙型储层动静态渗透率转换方法，其特征在于，步骤4)中，当每个分布范围内的有效孔隙度与静态渗透率之间的函数关系能够达到所对应的一类微裂缝-孔隙型储层渗透率精度的80％，即认定有效孔隙度与静态渗透率之间的函数关系均满足所对应的一类微裂缝-孔隙型储层静态渗透率计算的精度。
6.根据权利要求1所述的微裂缝-孔隙型储层动静态渗透率转换方法，其特征在于，还包括静态渗透率的获取：基于不同测井曲线数值的组合，采用微裂缝-孔隙型储层的分类模型去划分一口井整个井段的类别；针对每个井段微裂缝-孔隙型储层所属的类别，采用相应的有效孔隙度与静态渗透率之间的函数关系来计算该类别中微裂缝-孔隙型储层的静态渗透率。
7.根据权利要求1所述的微裂缝-孔隙型储层动静态渗透率转换方法，其特征在于，所述步骤6)中以微裂缝孔隙度的大小为分类标准，当微裂缝孔隙度大于某一数值时，井中所发育的储层为微裂缝-孔隙型储层；当微裂缝孔隙度小于某一数值时，井中所发育的储层则为孔隙型储层；然后分别建立微裂缝-孔隙型储层、孔隙型储层静态渗透率与动态渗透率之间的转换关系。
8.根据权利要求7所述的微裂缝-孔隙型储层动静态渗透率转换方法，其特征在于，所述步骤6)中的微裂缝-孔隙型储层、孔隙型储层静态渗透率与动态渗透率之间的转换关系为：
Kd＝mKs+n
式中，Kd为动态渗透率；Ks为静态渗透率；m、n为常数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911363854.5A CN112145165B (zh) | 2019-12-26 | 2019-12-26 | 一种微裂缝-孔隙型储层动静态渗透率转换方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911363854.5A CN112145165B (zh) | 2019-12-26 | 2019-12-26 | 一种微裂缝-孔隙型储层动静态渗透率转换方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112145165A CN112145165A (zh) | 2020-12-29 |
| CN112145165B true CN112145165B (zh) | 2024-04-23 |
Family
ID=73892085
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911363854.5A Active CN112145165B (zh) | 2019-12-26 | 2019-12-26 | 一种微裂缝-孔隙型储层动静态渗透率转换方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112145165B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112966383B (zh) * | 2021-03-11 | 2023-10-20 | 中海石油(中国)有限公司 | 一种裂缝性储层的产能预测方法 |
| CN113109890A (zh) * | 2021-04-30 | 2021-07-13 | 中海石油(中国)有限公司 | 一种裂缝有效性评价方法 |
| CN116661019B (zh) * | 2023-06-02 | 2024-03-08 | 中国石油天然气股份有限公司 | 微裂缝的评价方法及装置 |
Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2268979A1 (en) * | 1997-01-31 | 1998-08-06 | Phillips Petroleum Company | Method for determining distribution of reservoir permeability, porosity and pseudo relative permeability |
| CA2438134A1 (en) * | 2002-09-30 | 2004-03-30 | William Daniel Gunter | Process for predicting porosity and permeability of a coal bed |
| CA2615392A1 (en) * | 2003-08-07 | 2005-02-07 | Baker Hughes Incorporated | Method of compensating measurements made by a downhole survey instrument |
| WO2010148628A1 (zh) * | 2009-06-22 | 2010-12-29 | 中国石油天然气股份有限公司 | 裂缝储层含油气饱和度定量计算方法 |
| CN102852516A (zh) * | 2012-04-19 | 2013-01-02 | 北京大学 | 用于油气藏开发的全缝长三维压裂数值模拟的方法和装置 |
| CN104500050A (zh) * | 2014-10-31 | 2015-04-08 | 中国石油大学(华东) | 一种裂缝性储层渗透率张量及各向异性定量预测方法 |
| CN104948176A (zh) * | 2015-05-08 | 2015-09-30 | 西南石油大学 | 一种基于渗透增大率识别碳酸盐岩储层裂缝的方法 |
| WO2016011064A2 (en) * | 2014-07-15 | 2016-01-21 | Petroleum Fractured Reservoir Solutions, Llc | Discrete irregular cellular models for simulating the development of fractured reservoirs |
| CN105319603A (zh) * | 2015-11-06 | 2016-02-10 | 中国石油大学(华东) | 致密砂岩储层复杂网状裂缝的预测方法 |
| CN107622139A (zh) * | 2016-07-15 | 2018-01-23 | 中国石油天然气股份有限公司 | 裂缝渗透率的计算方法 |
| CN108694264A (zh) * | 2017-04-11 | 2018-10-23 | 中国石油化工股份有限公司 | 一种确定页岩气储层渗透率的方法 |
| CN110043254A (zh) * | 2019-05-11 | 2019-07-23 | 中国海洋石油集团有限公司 | 一种基于电缆地层测试资料地层有效渗透率的获取方法 |
| CN110471127A (zh) * | 2019-08-06 | 2019-11-19 | 中国石油集团川庆钻探工程有限公司 | 一种测井渗透率解释方法 |
| RU2707311C1 (ru) * | 2019-09-06 | 2019-11-26 | федеральное государственное автономное образовательное учреждение высшего образования "Российский государственный университет нефти и газа (национальный исследовательский университет) имени И.М. Губкина" | Способ оценки профиля фазовой проницаемости в нефтяных и газовых эксплуатационных скважинах |
| CN110530910A (zh) * | 2019-08-15 | 2019-12-03 | 中国石油大学(北京) | 一种模拟致密岩石微纳孔隙环境的油气赋存相态测定方法 |
-
2019
- 2019-12-26 CN CN201911363854.5A patent/CN112145165B/zh active Active
Patent Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2268979A1 (en) * | 1997-01-31 | 1998-08-06 | Phillips Petroleum Company | Method for determining distribution of reservoir permeability, porosity and pseudo relative permeability |
| CA2438134A1 (en) * | 2002-09-30 | 2004-03-30 | William Daniel Gunter | Process for predicting porosity and permeability of a coal bed |
| CA2615392A1 (en) * | 2003-08-07 | 2005-02-07 | Baker Hughes Incorporated | Method of compensating measurements made by a downhole survey instrument |
| WO2010148628A1 (zh) * | 2009-06-22 | 2010-12-29 | 中国石油天然气股份有限公司 | 裂缝储层含油气饱和度定量计算方法 |
| CN102852516A (zh) * | 2012-04-19 | 2013-01-02 | 北京大学 | 用于油气藏开发的全缝长三维压裂数值模拟的方法和装置 |
| WO2016011064A2 (en) * | 2014-07-15 | 2016-01-21 | Petroleum Fractured Reservoir Solutions, Llc | Discrete irregular cellular models for simulating the development of fractured reservoirs |
| CN104500050A (zh) * | 2014-10-31 | 2015-04-08 | 中国石油大学(华东) | 一种裂缝性储层渗透率张量及各向异性定量预测方法 |
| CN104948176A (zh) * | 2015-05-08 | 2015-09-30 | 西南石油大学 | 一种基于渗透增大率识别碳酸盐岩储层裂缝的方法 |
| CN105319603A (zh) * | 2015-11-06 | 2016-02-10 | 中国石油大学(华东) | 致密砂岩储层复杂网状裂缝的预测方法 |
| CN107622139A (zh) * | 2016-07-15 | 2018-01-23 | 中国石油天然气股份有限公司 | 裂缝渗透率的计算方法 |
| CN108694264A (zh) * | 2017-04-11 | 2018-10-23 | 中国石油化工股份有限公司 | 一种确定页岩气储层渗透率的方法 |
| CN110043254A (zh) * | 2019-05-11 | 2019-07-23 | 中国海洋石油集团有限公司 | 一种基于电缆地层测试资料地层有效渗透率的获取方法 |
| CN110471127A (zh) * | 2019-08-06 | 2019-11-19 | 中国石油集团川庆钻探工程有限公司 | 一种测井渗透率解释方法 |
| CN110530910A (zh) * | 2019-08-15 | 2019-12-03 | 中国石油大学(北京) | 一种模拟致密岩石微纳孔隙环境的油气赋存相态测定方法 |
| RU2707311C1 (ru) * | 2019-09-06 | 2019-11-26 | федеральное государственное автономное образовательное учреждение высшего образования "Российский государственный университет нефти и газа (национальный исследовательский университет) имени И.М. Губкина" | Способ оценки профиля фазовой проницаемости в нефтяных и газовых эксплуатационных скважинах |
Non-Patent Citations (3)
| Title |
|---|
| 基于CT的数字岩心三维建模;林承焰;王杨;杨山;任丽华;由春梅;吴松涛;吴玉其;张依旻;;吉林大学学报(地球科学版)(第01期);全文 * |
| 基于海上砂砾岩低渗透率成因分析及测井评价;张冲;;测井技术(第05期);全文 * |
| 电缆地层测试资料在低渗气藏产能预测中的应用;杨冬;张海荣;吴一雄;张恒荣;谭伟;;中国海上油气(第06期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112145165A (zh) | 2020-12-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101487390B (zh) | 一种确定油层原始含油饱和度的阿尔奇模式方法 | |
| CN102175832B (zh) | 一种确定典型储层最佳饱和度计算模型的方法 | |
| CN112145165B (zh) | 一种微裂缝-孔隙型储层动静态渗透率转换方法 | |
| CN110847901B (zh) | 一种变矿化度地层水下的致密砂岩储层流体识别方法 | |
| CN108694264B (zh) | 一种确定页岩气储层渗透率的方法 | |
| US20070112518A1 (en) | Method and Apparatus for Measuring the Wettability of Geological Formations | |
| CN104213899B (zh) | 一种地层岩石骨架的测井识别方法 | |
| CN109138975B (zh) | 一种基于时移测井数据的求解相渗特征曲线的新方法 | |
| CN104345346A (zh) | 一种获取裂缝宽度的方法 | |
| CN107829731B (zh) | 一种黏土蚀变的火山岩孔隙度校正方法 | |
| MX2014012042A (es) | Evaluacion volumetrica de una formacion usando datos diferenciales normalizados. | |
| CN107762483B (zh) | 一种基于测井曲线的相关系数及包络面积的流体识别方法 | |
| CN111381292B (zh) | 一种预测砂岩含烃储层的测井解释方法与装置 | |
| CN106777649B (zh) | 一种裂缝型储集层孔隙结构定量评价方法 | |
| Jasim et al. | Specifying quality of a tight oil reservoir through 3-d reservoir modeling | |
| CN112115592A (zh) | 一种基于岩石物理实验的大斜度/水平井电阻率校正方法 | |
| CN115857047B (zh) | 一种地震储层综合预测方法 | |
| CN112835124B (zh) | 基于成像测井和阵列声波测井资料的裂缝有效性评价方法 | |
| CN115099014A (zh) | 一种基于随钻录井的天然气井地质探明储量计算方法 | |
| CN103675945A (zh) | 一种测定孔洞型储层的饱和度的方法及设备 | |
| CN109994161B (zh) | 趋势基线法结合动态联动法计算地层有机碳含量的方法 | |
| CN116146170A (zh) | 一种三维感应各向异性泥页岩储层饱和度定量计算方法和装置 | |
| CN113720745A (zh) | 含碳屑碎屑岩储层地球物理测井计算孔隙度的方法 | |
| KR101818098B1 (ko) | 물리검층 자료를 이용한 지층 내 점토 함량 산출방법 | |
| CN115079261B (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 |