对于孔隙结构复杂的储层,现有的基于常规测井的渗透率模型计算精度低,无法满足生产需要。基于孔隙结构对储层进行分类,并根据测井资料的丰富程度提出两类渗透率计算模型:对于常规测井曲线齐全的井,利用测井曲线计算储层流动单元指数FZI,并进行储层分类识别,每一类储层单独建立渗透率模型;对于有阵列声波测井曲线的井,通过交会图筛选渗透率敏感参数,并建立基于阵列声波的渗透率模型。研究结果表明:声波时差(AC)、补偿密度(DEN)、补偿中子(CNL)测井曲线可准确计算流动单元指数FZI并识别I类、II类、III类储层,常规测井和阵列声波测井获得的密度曲线、泥质含量、斯通利波时差和泊松比等参数与渗透率相关性较高,基于这些参数建立的三因子指数模型对I类储层和II类储层适用性较好,双因子幂乘模型对III类储层适用性较好。基于流动单元和基于阵列声波的渗透率模型提高了复杂孔隙结构储层的渗透率计算精度,有效地支撑了气藏的高效勘探与开发。 For reservoir with complex pore structure, the calculation accuracy of existing permeability model based on conventional logging is low and cannot meet production needs. The reservoir is classified based on pore structure, and two permeability calculation models are proposed according to the abundance of logging data. For Wells with complete conventional logging curves, the reservoir flow unit index FZI is calculated by logging curves, and the reservoir is classified and identified. For Wells with array acoustic logging curves, permeability sensitive parameters are screened by cross plot and permeability model based on array acoustic logging is established. The results show that: Acoustic time difference (AC), compensated density (DEN), and compensated neutron (CNL) logs can accurately calculate flow unit index FZI and identify class I, II, and III reservoirs. Density curves, shale content, Stoneley wave time difference, poisson’s ratio and other parameters obtained from conventional and array acoustic logs are highly correlated with permeability. The three-factor exponential model based on these parameters is suitable for class I reservoir and class II reservoir, and the two-factor power model is suitable for class III reservoir. The permeability model based on flow unit and array acoustic wave improves the calculation accuracy of permeability of complex pore structure reservoir and supports efficient exploration and development of gas reservoir effectively.
对于孔隙结构复杂的储层,现有的基于常规测井的渗透率模型计算精度低,无法满足生产需要。基于孔隙结构对储层进行分类,并根据测井资料的丰富程度提出两类渗透率计算模型:对于常规测井曲线齐全的井,利用测井曲线计算储层流动单元指数FZI,并进行储层分类识别,每一类储层单独建立渗透率模型;对于有阵列声波测井曲线的井,通过交会图筛选渗透率敏感参数,并建立基于阵列声波的渗透率模型。研究结果表明:声波时差(AC)、补偿密度(DEN)、补偿中子(CNL)测井曲线可准确计算流动单元指数FZI并识别I类、II类、III类储层,常规测井和阵列声波测井获得的密度曲线、泥质含量、斯通利波时差和泊松比等参数与渗透率相关性较高,基于这些参数建立的三因子指数模型对I类储层和II类储层适用性较好,双因子幂乘模型对III类储层适用性较好。基于流动单元和基于阵列声波的渗透率模型提高了复杂孔隙结构储层的渗透率计算精度,有效地支撑了气藏的高效勘探与开发。
流动单元,阵列声波,渗透率,斯通利波
Ling Fan1, Yanping Xi2, Jin Dai3, Yingying Luo3, Menghan Wen4, Yun Liang3, Feng Wu3
1Middle Sichuan Oil and Gas Field, PetroChina Southwest Oil & Gas Field Company, Suining Sichuan
2Institute of Geological Exploration and Development, Petrochina Chuanqing Drilling Engineering Company, Chengdu Sichuan
3School of Geosciences and Technology, Southwest Petroleum University, Chengdu Sichuan
4Technical Consulting Center, Petrochina Southwest Oil and Gas Field Company, Chengdu Sichuan
Received: Nov. 1st, 2021; accepted: Nov. 19th, 2021; published: Nov. 29th, 2021
For reservoir with complex pore structure, the calculation accuracy of existing permeability model based on conventional logging is low and cannot meet production needs. The reservoir is classified based on pore structure, and two permeability calculation models are proposed according to the abundance of logging data. For Wells with complete conventional logging curves, the reservoir flow unit index FZI is calculated by logging curves, and the reservoir is classified and identified. For Wells with array acoustic logging curves, permeability sensitive parameters are screened by cross plot and permeability model based on array acoustic logging is established. The results show that: Acoustic time difference (AC), compensated density (DEN), and compensated neutron (CNL) logs can accurately calculate flow unit index FZI and identify class I, II, and III reservoirs. Density curves, shale content, Stoneley wave time difference, poisson’s ratio and other parameters obtained from conventional and array acoustic logs are highly correlated with permeability. The three-factor exponential model based on these parameters is suitable for class I reservoir and class II reservoir, and the two-factor power model is suitable for class III reservoir. The permeability model based on flow unit and array acoustic wave improves the calculation accuracy of permeability of complex pore structure reservoir and supports efficient exploration and development of gas reservoir effectively.
Keywords:Flow Unit, Array Acoustic, Permeability, Stoneley Wave
Copyright © 2021 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
渗透率是评价储层好坏的重要参数之一,现有的地层渗透率计算模型主要是基于微观参数和常规测井曲线建立的。为了提高储层渗透率评价的精度,大量学者尝试了新的方法或运用了大量新的测井数据。流动单元法是近年来用在渗透率评价上较多的一种方法,并且取的了一定的成果 [
另外,阵列声波测井资料可以获得井壁地层更多的信息,通过对接收到的波列处理,可获得纵波、横波、斯通利波时差及泊松比,并用于评价储层的岩石力学特性、含气性和渗透性。前人研究表明在岩心渗透率标定的基础上,可根据斯通利波来进行地层渗透的简单估算,但在储层孔隙结构复杂的地层中,仅依靠斯通利波进行渗透率计算精度较低 [
针对研究区地层孔隙结构复杂导致的孔渗关系精度低、现有流动单元法无法实现连续渗透率计算、阵列声波测井资料利用程度低这三大问题,本文拟先通过核磁资料对研究区储层进行分类,再利用测井曲线计算FZI (流动单元指数)进行储层的测井分类,并利用流动单元理论和阵列声波测井资料,分别建立渗透率评价模型,为研究区后续的作业决策、储量计算和开发方案提供依据。
研究区发育三角洲平原~三角洲前缘沉积,沉积物粒度相对较粗,岩性包括中砂岩、粗砂岩和含砾砂岩等。储层孔隙度主要在6%~20%之间,渗透率主要在0.01~1000 mD之间。由于储层孔隙度和渗透率的范围非常大,导致储层渗透率的准确计算非常困难。
图1. I类、II类、III类储层孔隙结构特征。(a)、(b)、(c)为岩石饱和水状态T2分布;(d)、(e)、(f)为岩石束缚水状态T2分布;(g)、(h)、(i)为岩石可动水状态T2分布;(j)、(k)、(l)为铸体薄片照片
采用全直径核磁共振岩心分析仪AniMR-Hole Core对岩心样品开展核磁共振测量,得到岩心饱和水状态、岩心束缚水状态的核磁T2分布。结合岩心核磁资料和铸体薄片资料,将研究区储层划分为三类储层:I类储层(中高渗储层)、II类储层(中低渗储层)、III类储层(低渗储层) (图1)。I类储层孔隙类型主要为粒间孔,孔径主要区间0.45~0.8 mm,T2谱主峰值约600 ms,孔隙度 > 13%,渗透率 > 10 mD;Ⅱ类储层孔隙类型主要为粒内溶孔和粒间孔,孔径主要区间0.25~0.45 mm,核磁共振的T2谱主峰值约70 ms,孔隙度为10%~13%,渗透率为1~10 mD;Ⅲ类储层孔隙类型主要为粒内溶孔,孔径主要区间0.01~0.25 mm,核磁共振的T2谱主峰值约20 ms,孔隙度 < 10%,渗透率 < 1 mD。
根据平均水单元动力概念,Amaefule等人通过将孔隙空间看成很多个毛细管,应用了Darcy定律以及Poisscuille定律,进而得到了可以表征不同流动单元间渗透率与孔隙度关系的方程 [
K = ϕ e 3 ( 1 − ϕ e ) 2 1 F s τ 2 S g v 2 (1)
式中,K为渗透率,mD; ϕ e 为有效孔隙度,小数; F s 为孔隙的几何形状指数,对于圆柱形,其值为2; S g v 是矿物颗粒比表面积,小数; τ 为流动路径的弯曲度,小数。
在Kozeny-Carman方程中,一般将 F s τ 2 称为Kozeny常数。在确定的流动单元内它是一个常数,但若在不同的流动单元之间,该值也会发生相应的改变。流动单元指数FZI的定义如下:
F Z I = 1 τ S g v F s = 0.0314 1 − ϕ e ϕ e K ϕ e (2)
式中,FZI为流动单元指数。
通过研究流动单元指数和通过实验分析的矿物颗粒表面、束缚水饱和度和矿物颗粒比重间的关系,对于比较粗的,分选性比较好的砂岩,它的流动单元指数相对较高,而对于具有高弯曲度的砂岩,它的分选性一般来说比较差,流动单元指数也比较低,基于流动单元指数能够更好地计算岩石的渗透率。
将公式2中的部分表达式提取出来,分别定义为储层品质因子RQI和孔隙度指标 ϕ z (公式3、公式4),则流动单元指数FZI可表示为储层品质因子RQI和孔隙度指标 ϕ z 的形式 [
R Q I = K ϕ e (3)
式中,RQI为储层品质因子。
ϕ z = ϕ e 1 − ϕ e (4)
式中, ϕ z 为孔隙度指标。
F Z I = 0.0314 R Q I ϕ z (5)
对公式5等号两边分别取对数,可得:
lg R Q I = a lg ϕ z + lg F Z I (6)
式中,a为系数。
以岩心物性分析资料进行标定,通过声波时差(AC)、补偿密度(DEN)、补偿中子(CNL)测井曲线建立流动单元指数FZI的测井计算模型(公式7)。基于FZI值划分I类、II类、III类储层的FZI界限:I类储层(中高渗储层),FZI ≥ 2.5;II类储层(中渗储层),2.5 > FZI ≥ 1;III类储层(低渗储层),FZI < 1 (图2)。根据公式6绘制以lg(RQI)为横坐标、 lg ( ϕ z ) 为纵坐标的交会图(图3(a)),数据点呈现明显的分带特征,研究区三类储层的数据点可拟合得到三条斜率相近但截距不同的平行线,由此实现了流动单元的测井分类与预测。
F Z I = 0.182 A C − 14.578 D E N − 0.389 C N L + 30.07 (7)
式中,AC为补偿声波测井值,DEN为补偿密度测井值,CNL为补偿中子测井值。
图2. 基于FZI的储层流动单元划分结果
图3. FZI分类与孔渗关系。(a) FZI分类;(b) 基于FZI分类的孔渗关系
在测井计算FZI值并进行I类、II类、III类储层划分的基础上,对这三类储层的数据点进行孔隙度~渗透关系拟合,分别得到I类、II类、III类储层的孔隙度~渗透率拟合关系(图3(b))。在相同孔隙度条件下,I类储层的渗透率最高,II类储层的渗透率其次,III类储层的渗透率最低。基于流动单元FZI分类的I类、II类、III类储层渗透率模型如表1所示。
储层类型 | FZI分类标准 | 渗透率模型 |
---|---|---|
I类储层 | FZI ≥ 3 | K = 0.0072 × Φ 3.7553 |
II类储层 | 3 > FZI ≥ 1 | K = 0.0009 × Φ 3.4991 |
III类储层 | FZI < 1 | K = 0.0003 × Φ 3.1661 |
表1. 基于流动单元FZI分类的渗透率模型
目前基于测井曲线建立的渗透率模型广泛运用了核磁共振、电阻率和自然伽马测井等 [
对比分析现有的经典渗透率模型计算公式,大多采用指数乘积形式,也有采用指数和、幂乘积 [
K = m T 1 t 1 T i t i ⋯ T n t n (7)
K = 10 t 1 T 1 + t i T i + ⋯ + t n T n + m (8)
K = t 1 T 1 + t i T i + ⋯ + t n T n + m (9)
式中, K 为渗透率; m 为系数; T n 为影响因素,影响因素包括孔隙度、泥质含量等; T i 为第 i 个影响因素, t n 为影响因素系数; t i 是第 i 个影响因素 T i 的系数; t i 为常数, i 为整数, n 为整数, 0 < i < n 。
测井 | 研究人员 | 渗透率模型 | 影响因素 |
---|---|---|---|
核磁共振 | Swanson (1981) | k = 4.6 ∅ 4 T 2 lm 2 | 孔隙度、T2对数平均值 |
Coates (1991) | k = [ ( ∅ / C ) 2 ( F F I / B V I ) ] | 孔隙度、自由流体指数 | |
Quintero et al. (1999) | k = C p f 4.6 ∅ 4 T 2 lm 2 | 孔隙度、T2对数平均值 | |
Sen et al. (1990) | k = 10 − 0.1 ( ∅ m T 2 ) 2.15 | 孔隙度、T2值 | |
电阻率 | Tixier (1949) | k = C [ 2.3 R o ( ρ w − ρ o ) Δ R Δ D ] 2 | 电阻率、地层水密度、 地层油密度 |
Coates and Dummannoir (1974) | k = C [ C ∅ 2 w w 4 ( R w / R t ) ] 2 | 电阻率、孔隙度、胶结和 弯曲度因素 | |
Yao and Holditch (1993) | k = ∅ e 1 ( 1 − I G R ) e 2 R i l d e 3 ( R i l d / R s f l ) e 4 | 自然伽马、深浅侧向电阻率、孔隙度 | |
自然伽马 | Saner et al. (1997) | log ( K ) = 7.04 − 4.19 ( F ) | 地层因素(与电阻率有关) |
Mohaghegh et al. (1997) | k = 126.5 + 0.0011 GR − 50.3 DEN + 0.625 R D | 自然伽马、密度、电阻率 | |
Xue et al. (1997) | log ( K ) = 0.151 Δ t − 0.019 ∅ d − 0.0392 G R + 0.0222 R D − 7.7 | 自然伽马、孔隙度、声波、 电阻率 |
表2. 已有基于测井曲线的渗透率模型对比
类型 | 影响因素 | 单相关关系 | 相关系数R2 |
---|---|---|---|
常规测井曲线 | 声波时差曲线/AC | K = 5 E − A C 40.311 | 0.6310 |
密度曲线/DEN | K = 2 E + 43 e − 40.16 ∗ D E N | 0.8163 | |
自然伽马曲线/GR | K = 5 E + 24 ∗ G R − 12.59 | 0.6629 | |
深侧向电阻率/RD | K = 4 E − 06 ∗ R D 4.6512 | 0.4722 | |
储层参数 | 孔隙度/POR | K = 0.0002 ∗ e 0.9669 ∗ P O R | 0.8359 |
泥质含量/VSH | K = 1 E + 07 ∗ V S H − 4.91 | 0.6361 | |
阵列声波曲线 | 斯通利波时差/DTST | K = 4 E − 99 ∗ e 1.0805 ∗ D T S T | 0.7233 |
横波时差/DTD | K = 10.656 ∗ D T D − 1013.4 | 0.0180 | |
泊松比/POSI | K = 53803 ∗ e − 55.16 ∗ P O S I | 0.7812 |
表3. 不同影响因素类型与渗透率的关系
形式 | 研究人员 | 渗透率模型 | 影响因素 |
---|---|---|---|
指数乘积 | Krumbein and Monk (1943) | k = 760 D g 2 e − 1.3 ∗ σ D | 晶粒直径 |
Berg (1970) | k = 80.8 ∅ 5.1 D 2 e − 1.385 p | 孔隙度、中值粒径 | |
Timur (1968) | k = 0.136 ∅ 4.4 / S w i r 2 | 孔隙度、束缚水饱和度 | |
Coates and Denoo (1981) | k = [ 100 ∅ e 2 ( 1 − S w i ) 2 / S w i ] 2 | 孔隙度、束缚水饱和度 | |
Sen et al. (1990) | k = 10 6.59 ( ∅ m V p / S ) 2.08 | 孔隙度、体积、表面比 | |
Sen et al. (1990) | k = 10 2.65 ( ∅ m / Q v ) 2.11 | 孔隙度、阳离子交换量 | |
Amaefule et al. (1993) | k = 1014 ( F Z I ) 2 [ ∅ e 3 / ( 1 − ∅ e ) 2 ] | 孔隙度、流动带指标 | |
Katz and Thomson (1986) | k = 226 − 1 l c 2 ( σ / σ o ) | 孔隙直径、电导率 | |
Martys et al. (1994) | k = ( 2 ∅ 2 ∗ / S 2 ) ( ∅ 1 − ∅ 1 c ) f | 孔隙度、比表面积 | |
Van Baaren (1979) | k = 10 D d 2 ∅ 3.64 + m C − 3.64 | 孔隙度、颗粒尺寸 | |
Hansen and Skieltorp (1988) | k = c ′ ∅ ( l 1 / l n ) 2 [ ( E s − D s ) + ( D v − E v ) ] | 孔隙度、长度、欧氏维数、分型维数 | |
Mavko and Nur (1997) | k ∞ c ( ∅ − ∅ c ) 3 d 2 | 孔隙度、颗粒大小 | |
指数和 | Pape et al. (1999) | k = A ∅ + B ∅ exp 1 + C ( 10 ∅ ) exp 2 | 孔隙度 |
幂乘积 | Garrison et al. (1993) | k = 10 − 1.7 − 0.6 D ′ s ( c o n t r o l ) + 113 S ′ a ( c o n t r o l ) | 视表面分形维数、面积形状因子 |
表4. 已有经典渗透率模型对比
综合渗透率影响因素筛选和渗透率模型形式选择结果,分两组参数建立渗透率模型:1) 选取孔隙度(POR)、斯通利波时差(DTST)两个参数,分别采用指数乘积形式、幂乘积形式和线性形式建立渗透率模型;2) 选取孔隙度(POR)、泥质含量(VSH)、斯通利波时差(DTST)三个参数,分别采用指数乘积形式、幂乘积形式和线性形式建立渗透率模型(表5)。
选取参数 | 公式形式 | 渗透率模型 | 相关系数R2 | 公式编号 |
---|---|---|---|---|
孔隙度(POR) | 幂 | K = 0.0009 ∗ e 0.84 ∗ Φ | 0.9163 | POR-K |
孔隙度(POR) 斯通利声波时差(DTST) | 指数 | K = 10 − 235.566 ∗ Φ 4.4364 ∗ D T S T 99.7999 | 0.9142 | K-1-1 |
幂 | K = 10 0.3239 ∗ Φ + 0.1173 ∗ D T S T − 27.4206 | 0.9264 | K-1-2 | |
线性 | K = − 39.8632 ∗ Φ + 114.4274 ∗ D T S T − 23594.13 | 0.5662 | K-1-3 | |
孔隙度(POR) 斯通利声波时差(DTST) 泥质含量(VSH) | 指数 | K = 10 − 285.594 ∗ Φ 1.99146 ∗ D T S T 123.506 ∗ V S H − 2.1313 | 0.9434 | K-2-1 |
幂 | K = 10 0.1773 ∗ Φ + 0.1950 ∗ D T S T − 0.0277 ∗ V S H − 41.6369 | 0.9428 | K-2-2 | |
线性 | K = − 58.4742 ∗ Φ + 121.2763 ∗ D T S T − 3.6354 ∗ V S H − 24757.2 | 0.6915 | K-2-3 |
表5. 基于阵列声波测井的渗透率模型
图4. 不同影响因素与渗透率的相关性
根据表5建立的渗透率模型,对研究区X井进行渗透率计算(图5)。从图5可以看出双因子及三因子的的渗透率模型较单因子孔隙度建立的渗透率模型有更高的预测精度。结合前述储层分类结果,对于I类储层和II类储层,三因子指数模型的效果较好;对于III类储层,双因子幂乘模型的效果较好。采用流动单元拟合法计算的渗透率比采用原始回归模型进行渗透率预测的效果更好。采用流动单元分类拟合法适用于常规测井曲线齐全的井,基于阵列声波的渗透率模型适用于测量了阵列声波测井曲线的井,这两类模型比不分类的孔隙度~渗透率拟合模型,计算渗透率的精度更高(图6)。
图5. 不同渗透率模型计算结果对比
图6. 不同渗透率模型计算精度对比
1) 研究区储层可分为三类:I类储层孔隙类型主要为粒间孔,T2谱主峰值约600 ms,孔隙度 > 13%,渗透率 > 10 mD;Ⅱ类储层孔隙类型主要为粒内溶孔和粒间孔,T2谱主峰值约70 ms,孔隙度为10%~13%,渗透率为1~10 mD;Ⅲ类储层孔隙类型主要为粒内溶孔,核磁共振的T2谱主峰值约20 ms,孔隙度 < 10%,渗透率 < 1 mD。
2) 基于流动单元的渗透率模型适用于常规测井曲线齐全的井,基于测井曲线建立流动单元指数FZI计算公式,可将在全井段实现I类储层、II类储层和III类储层的流动单元分类,最后通过流动单元分类的渗透率模型可提高渗透率计算精度。
3) 常规测井和阵列声波测井获得的密度曲线、泥质含量、斯通利波时差和泊松比等参数与渗透率相关性较高,基于这些参数建立的三因子指数模型更适合I类储层和II类储层、双因子幂乘模型更适合III类储层。
四川省重点研发计划(重大科技专项) (编号:2020YFSY0039);西南石油大学“油气藏地质及开发工程”国家重点实验室开放基金(编号:PLN201719)。
范 玲,习研平,代 槿,罗莹莹,温梦晗,梁 芸,吴 丰. 基于流动单元和阵列声波测井的储层渗透率计算Flow Unit and Array Acoustic Wave Logging-Based Reservoir Permeability Calculation[J]. 地球科学前沿, 2021, 11(11): 1530-1541. https://doi.org/10.12677/AG.2021.1111147