孙振介¹，郑 军¹，刘鸿博²，左银辉¹

¹成都理工大学，四川 成都

²四川水利职业技术学院，四川 成都

收稿日期：2019年1月29日；录用日期：2019年2月13日；发布日期：2019年2月21日

岩石渗透率是一个典型的矢量，渗透率各向异性的测试和评价方法越来越受到人们重视。本文针对传统渗透率各向异性评价方法不能充分体现极值点影响、物理意义不明确、不便于岩心间的对比和归类等缺点，为了能够更加明确地表示储层岩石渗透率各向异性的特点，通过对前人评价岩石渗透率方法的研究，结合对不同岩性储层岩石的渗透率测试，提出了新的评价方法。包括用三向渗透率的几何平均代表岩心的视渗透率，用三角图表示各向异性渗透率的相对大小关系，用三向渗透率几何平均与算术平均比值来定量描述各向异性程度。通过对实际不同岩性岩心渗透率各向异性的具体应用，结果表明该新方法能科学、全面评价岩石渗透率各向异性。

关键词 :渗透率各向异性，视渗透率，三角图，各向异性程度

Copyright © 2019 by authors and beplay安卓登录

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

油气储层岩石受自身沉积、成岩作用，以及受地下非均质应力环境等影响，其渗透率往往表现出明显的各向异性，渗透率是一个矢量。储层在地下环境中，渗透率的各向异性是油气田开发的重要岩石物性参数，它将影响油藏布井方式、井距大小、水平井水平段延展方向、人工压裂方案等。目前有关油气藏数值模拟的软件基本都能处理渗透率各向异性模型。采用渗透率张量对油藏进行描述，能更准确反映油气等流体在地下储层中流动规律，有效提高油藏数值模拟研究精度，为油气藏开发方案及提高采收率措施等提供理论基础，对于指导油气田开发具有重要的意义。

国外早在上世纪30时代就已经开展了对于渗透率方向性的研究 [ 1 ] ，上世纪60年代国外已有专门测试渗透率方向性的具体实验方法 [ 2 ] [ 3 ] [ 4 ] [ 5 ] 。我国对渗透率方向性研究起步较晚，国内学者们一直都在探索新型的渗透率方向性描述与评价方法 [ 6 ] - [ 16 ] 。一般认为岩石在一个点不同方向渗透率可以构成如图1所示的椭球体，目前有多种实验测试手段可以得到岩心分别在x，y，z方向的渗透率K_x、K_y、K_z，渗透率各向异性描述方法也是针对这三个方向渗透率值展开的，主要涉及三方面的内容：其一为渗透率各向异性的图示方法，采用恰当的作图方式对各向异性渗透率进行直观描述(图1)。其二为矢量渗透率的平均方法，由于岩心在每一个方向上渗透率都不同，因此需要提取一个特征渗透率值(或视渗透率)，代表整个岩心渗透率的大小。最后是渗透率各向异性程度的表示方法，表征不同方向渗透率值的差异程度。

图1. 矢量渗透率所构成的椭球体示意图

现有技术一般用如图2所示的柱状图，把三向渗透率放在一起，直观描述渗透率的大小和差异程度 [ 17 ] [ 18 ] 。

图2. 渗透率各向异性图示方法示意图

特征渗透率采用三向渗透率的算术平均表示，即：

K ¯ = K x + K y + K z 3 (1)

定义 δ 1 为两向渗透率各向异性系数

δ 1 = K i / K j (2)

式中K_i，K_j分别为K_x、K_y和K_z中任意两个，用其来描述渗透率各向异性的程度。

定义 δ 2 为三向渗透率各向异性系数

δ 2 = ( ( K x − K ¯ K ¯ ) 2 + ( K y − K ¯ K ¯ ) 2 + ( K z − K ¯ K ¯ ) 2 ) / 3 (3)

用其来描述三向渗透率各向异性的程度 [ 19 ] [ 20 ] 。

传统渗透率各向异性图示方法由于参数较多，渗透率变化范围大，不利于不同岩心间的比较。 δ 1 只适用于评价两向渗透率极差。基于算术平均的特征渗透率及各向异性描述方法适用于无相互影响、独立存在的数据体评价，该方式不能体现极值数据点的影响，不适用于存在几何关系的三向渗透率的评价。

针对矢量渗透率的特点，本文提出了新的各向异性渗透率的描述方法，以克服传统方法不能充分体现三向渗透率极值点的影响，物理意义不够明确，且不便于不同岩心间的对比和归类等缺点。

储层岩石的渗透率是评价储层优劣的主要指标，渗透率的大小则能最直观的反应储层性质。由于岩石的渗透率方向性在地层中受应力影响较为明显，渗透率与应力同为矢量，且一般情况下，K_x、K_y、K_z为不同值，如何用一个特征值来表征渗透率的大小呢？

以往提出的渗透率方向性特征值大都是运用三向渗透率的算术平均值进行计算。我们知道，算术平均值只适合于描述相互独立的多个数据点之间的平均。运用算术平均值推算渗透率方向性特征值时，受极端值影响较大，得出的结果不能够准确表示渗透率的大小。而三向渗透率K_x、K_y、K_z之间是存在相互关系且代表着椭球体在三个轴向上的值。把岩石样品三向渗透率数据K_x、K_y、K_z的几何平均值作为渗透率方向性特征值，受极端值影响较小，能够更加准确的表示渗透率的大小。三向渗透率的几何平均可用下式表示：

K a = K x K y K z 3 (4)

几何平均数考虑了数据间的几何关系，且能更有效反应数据极值的影响，当已知数据呈倍数关系或不对称分布时，几何平均数的效果更好。除此之外，三向渗透率的几何平均还具有明确的物理意义，即为以几何平均数K_a为半径的球体与以三向渗透率K_x、K_y、K_z为半径的椭球体具有相同的体积，如图3所示，而矢量渗透率所包罗的体积能反应该岩石的综合渗透能力，因此用几何平均K_a代表各向异性渗透率的特征值更具有科学性。

图3. 特征渗透率Ka为半径所构成的球体示意图

三角形统计图是采用数字坐标形式来表现三项要素的数字信息图像。三角形坐标图常用百分数(%)来表示某项要素与整体的结构比例。三条边分别表示三个不同的要素，三个顶点可以看作是三个原点。三角形统计图的适用范围比较广，可用来表示三个项目所占整体的百分比。

已知岩石样品的三向渗透率分别为K_x、K_y、K_z，对其进行归一化可得：

P x = K x K x + K y + K z × 100 %

P y = K y K x + K y + K z × 100 % (5)

P z = K z K x + K y + K z × 100 %

根据归一化后的三向渗透率值P_x、P_y、P_z，可在三角图中进行标识，把三个要素体现在同一点上，便于不同岩心之间的比较，表征渗透率各向异性的三角图如图4所示。根据三向渗透率大小之间的关系，渗透率各向异性的三角图可分为3大类区域，每个大类区域又可再划分为两个区域(图4)，每个区域点内的三向渗透率大小特征如表1所示。利用这样的原则相比于以往的方法，可以更加有效的描述储层岩石渗透率各向异性情况。

表1. 渗透率各向异性三角图分区及特征

图4. 描述储层渗透率方向性的三角网状图

在判断渗透率各向异性的同时，我们还需要知道渗透率各向异性程度的大小。可以利用三向渗透率几何平均与算术平均比值d来代表各向异性程度，即：

d = K x K y K z 3 1 3 ( K x + K y + K z ) (5)

在三角图内分别作出d = 1，d = 0.9，d = 0.7以及d = 0的等值线图，如图5所示。从图中可以看出，在三角形中心点d = 1，此时归一化后的K_x= K_y= K_z，各向同性，在三角形边线上的点d = 0，此时K_x、K_y、K_z中至少有一个为0，三向渗透率变为两向或单向渗透率。在过中心点同一直线上，离中心点距离越远，d值越小，非均质程度越高。可按表2标准对各向异性程度进行评价。

图5. 渗透率各项异性三角等值线图

表2. 渗透率各向异性程度三角图分区及特征图

本文从大小和方向两个方面评价了方向渗透率，该方法克服了无法将三向渗透率统一描述的困难，建立了全新的渗透率方向性评价与描述方法，具有以下优势：

1) 综合岩石的三向渗透率大小提取出描述渗透率大小的特征值，可对不同地层条件的储层或相同地层不同层位的储层渗透率大小进行比较判断；

2) 用三角图对岩石渗透率各向异性进行表征，通过对标点位置可以分析出样品渗透性特征；

3) 利用各向异性程度d和渗透率各项异性三角等值线图可以判断各样品的各向异性程度。

这种表征方法不仅可以判断各向渗透率相对大小，推断优势渗流的方向，还能对比各岩石间渗透率各向异性强弱。相比于利用渗透率的算数平均值描述渗透率各向异性的方法，本文提出的评价方法更加清晰直观，可操作性强，是岩石渗透率各向异性表征方法中较为科学的方法之一，值得使用和推广。

孙振介,郑 军,刘鸿博,左银辉. 岩石渗透率各向异性评价新方法 A New Method to Evaluate the Anisotropy of Rock Permeability[J]. 地球科学前沿, 2019, 09(02): 44-53. https://doi.org/10.12677/AG.2019.92006