基于等效毛细管的低渗透气藏液相侵入模型

谢晓永， 黄 敏

（中国石化西南油气分公司石油工程技术研究院，四川德阳 618001）

致密砂岩气藏具有孔渗特性差、水润湿性强、微裂缝发育和毛细管效应显著等特征，在钻井完井及储层改造过程中易遭受水锁损害，大幅度降低油气井产能[1-2]。针对水锁损害，众多学者开展了一系列实验与理论研究，建立了水锁预测模型[3-5]，探讨了水锁损害机理及防治措施[6-10]。赖南君等人[11]利用可对岩心加温加压的水锁损害评价仪研究了毛细管力与正压差对水锁损害的影响；刘建坤等人[12]将低磁场核磁共振T2谱技术与常规流动实验相结合，提出了水锁伤害核磁共振试验评价方法；唐洪明等人[13]利用核磁共振和T2谱技术开展了毛细管自吸与水相返排可视化实验；丁绍卿等人[14]将核磁共振技术应用于压裂液伤害机理研究，分析了黏土吸水效应及水锁效应对岩心渗透率的伤害程度。

以上方法多采用岩心驱替实验进行宏观规律研究，基于微观可视化的研究较少。为此，笔者建立了致密砂岩气藏孔隙网络模型，以期通过室内微观可视化实验分析液相侵入过程中孔隙网络内水相前缘的动态分布，揭示液相侵入过程中的微观流动机理；并在此基础上，建立了基于等效毛细管的低渗透气藏液相侵入微观流动模型，结合实验验证了模型的可行性，以期为低渗透气藏储层保护提供理论支撑。

1 孔隙网络微观可视化实验

1.1 微观流动实验装置

该实验中，将HC区块须家河组致密砂岩制成铸体薄片，利用图形扫描软件刻画孔隙网格，然后采用激光刻蚀致密砂岩孔隙网格。微观流动实验装置主要由计量管、微观流动实验装置、体视显微镜和量筒组成（见图1）。所选岩心孔隙度为9.23%，渗透率为0.27 mD。刻蚀材料为光学石英玻璃，尺寸30 mm×30 mm，刻蚀模型尺寸11 mm×8 mm。石英玻璃无涂层，孔隙流道亲水，且石英玻璃透光性好，便于采用光学显微镜观察水相的侵入与返排。实验流体为蒸馏水，由于孔隙网络微观模型尺寸小，实验开始时在刻蚀模型注入端预先注入10 cm高的水柱，由于水相侵入开始后不再补充水，随着实验的进行水柱高度下降。

图1 微观流动实验装置示意Fig. 1 Microscopic flow test device

1.2 微观流动实验结果

不同侵入时间下孔隙水相分布实验结果如图2所示。

从图2可以看出，侵入初期（＜8 s），在毛细管力和水柱压力作用下，水相侵入速度较快，迅速占据孔隙空间及其喉道；侵入中期，水相侵入速度减小，侵入深度缓慢增加，直至水侵前缘抵达模型右端出口；侵入后期，由于水柱压力降低且黏滞阻力增加，孔隙内水相流动能力下降，直至水侵呈稳定状态，水相侵入基本停止，侵入水相大多以残余水状态分布于孔隙网络模型中。由30 s时的水相分布可知，水相大部分位于孔喉处。根据实验结果，孔隙内水相侵入主要发生在流道中，孔道连接处的影响几乎可以忽略；同时，水相前缘推进与毛细管流动规律类似，可以为建立水相侵入模型提供实验依据。

水侵实验结束后，采用注射器返排侵入水相，当返排压力大于水相黏性阻力时，水相开始返排。返排后孔隙中的水相分布如图3所示。从图3可以看出，大部分孔隙中的水相得以返排，但仍有部分较小喉道中的水相未能返排。未能返排的水相以残余水状态存在，集中于孔喉处，阻碍气相流动。对于致密砂岩气藏，工作液与储层接触后，在井筒压差和毛细管力作用下沿孔隙侵入基质，由于致密砂岩通常具有亚束缚水饱和度特征，侵入水相可能在部分孔道壁面形成滞留水，即使后期采取负压差返排，孔喉处的滞留水仍不容易排出。

图3 水相返排后孔隙空间水相分布Fig. 3 Aqueous phase distribution in pore space after the flow back of aqueous phase

2 不同阶段水相侵入评价模型

由于水相侵入主要发生在孔隙流道中，孔喉处可以忽略，因此将致密砂岩孔隙网络简化为等效毛细管，忽略孔道壁面的摩擦系数。水相受到毛细管力、孔道壁面黏滞阻力、重力及惯性力作用[15]，根据主要作用力将水相侵入划分为初始阶段和平稳阶段，其中初始阶段包括惯性力侵入阶段、惯性力-黏滞阻力作用阶段及黏滞阻力作用阶段[16]。

2.1 惯性力侵入阶段

水相在惯性力侵入阶段开始进入毛细管，侵入量和侵入距离极小，主要作用力为毛细管力和惯性力，该阶段的瞬时侵入深度为[17]：

式中：l1为惯性力作用阶段水相侵入深度，m；t为时间，s；σ为表面张力，N/m；θ为接触角，（°）；ρ为流体密度，kg/m3；r为孔隙喉道半径，m；Fp为压差作用力，N。

等效毛细管半径采用Kozeny-Carman公式计算：

式中：K为渗透率，mD；ϕ为孔隙度。

惯性力作用阶段持续的时间为[17]：

式中：t1为惯性力作用阶段持续时间，s。

2.2 惯性力–黏滞阻力作用阶段

惯性力-黏滞阻力作用阶段的黏滞阻力不能忽略，此时水相侵入深度为[18]：

式中：l2为惯性力-黏滞阻力作用阶段的水相侵入深度，m。

该阶段的作用时间为[17]：

式中：t2为惯性力-黏滞阻力作用持续时间，s。

2.3 黏滞阻力作用阶段

水相侵入进入黏滞阻力作用阶段时，惯性效应可以忽略，动力平衡条件为[15]：

式中：α为自吸方向与水平方向的夹角，（°）；l3为平稳阶段水相侵入深度，m；为平稳阶段水相侵入速度，m/s；g 为重力加速度，m/s2。

由式（8）可得水相侵入速度的计算公式：

一维线性流动时，水相侵入过程中的压力分布为：

式中：L为水侵压力波及长度，m；pin为入口处的流体压力，Pa；pou为出口处的流体压力，Pa。

则压差作用力为：

将式（10）和式（11）代入式（9）求解。由于通常无法用解析方法求解，于是采用数值求解方法来求取近似解，笔者采用四阶Runge-Kutta差分格式求解：

式（12）即为建立的平稳阶段侵入模型，l3的初始值取t2时刻的l2。负压差条件下，水相侵入达到动力学平衡时的侵入深度可以认为是最大侵入深度：

式中：l3max为负压差下水相最大侵入深度，m。

3 模型验证

通过微观流动实验装置测量刻蚀模型毛细管束平均半径约为149.6 μm，表面张力取0.072 N/m，接触角取30°，黏度为1.0 mPa∙s，计算毛细管束不同水相侵入阶段水相自吸侵入深度，结果如图4所示。由式（2）和式（4）可知：惯性力和惯性力-黏滞阻力作用阶段持续时间均与毛细管半径成正比。由图4可知，惯性力作用阶段持续时间t1为0.52 ms时，水相自吸侵入深度为0.43 mm；惯性力-黏滞阻力作用阶段持续时间t2为47.0 ms时，水相自吸侵入深度为14.41 mm。之后水相侵入进入黏滞阻力作用阶段，水平方向自吸侵入不受重力的影响，侵入深度随时间增长不断增大；对于垂向自吸侵入，在重力作用下存在最大水相自吸侵入高度。算例条件下垂向水相最大自吸侵入高度为84.96 mm。致密砂岩孔喉半径通常介于0.03～2.00 μm[18]，可知惯性力和惯性力-黏滞阻力作用阶段持续时间短，水相自吸侵入深度主要取决于黏滞阻力作用阶段。

为了验证侵入模型（式（12））的可行性，采用西南石油大学研制的岩心自吸水测量仪，测试了致密砂岩垂向自吸侵入高度，实验装置及测试流程见文献[19]。岩心取自HC地区须家河组致密砂岩，渗透率0.24 mD，孔隙度9.3%，长度50.0 mm。实验结果表明，初始阶段水相侵入速度较快，之后侵入高度趋于平缓，实验数据与模型计算结果吻合度较高（见图5）。

图4 不同阶段毛细管水相自吸侵入深度Fig. 4 Self-absorption invasion depth of capillary aqueous phase at different stages

图5 模型预测结果与实验数据对比Fig. 5 Comparison of model prediction results with experimental data

通过式（13）可以计算负压差条件下最大水相侵入深度，此处采用负压差水相自吸实验数据进行验证。实验岩样取自川西蓬莱镇组和川中须家河组，岩样L54渗透率为0.17 mD，孔隙度为6.2%；岩样HE6渗透率为1.88 mD，孔隙度为11.8%[19]，实验欠压值1 MPa，实验结果与计算结果如图6所示。

图6 负压差条件下水相最大侵入深度Fig. 6 The maximum invasion depth of aqueous phase under negative pressure difference

从图6可以看出，负压差条件下，实验初期仍有水相侵入岩心，随着实验时间增长，水相侵入深度未明显增加。岩样L54和HE6的最大侵入深度的实验测试结果分别为3.1 和1.8 mm，计算结果分别为2.9和2.3 mm。可见，在相同实验条件下，岩石越致密，渗透率越低，水相最大侵入深度越大。

4 结 论

1）采用激光刻蚀技术雕刻了致密砂岩孔隙网络，开展了液相侵入微观可视化实验，分析了液相侵入过程中及返排后孔隙网络内的水相分布。实验结果表明，孔隙介质内的水相侵入主要发生在流道中，水相推进与毛细管流动规律类似，部分较小喉道中的水相不能返排。

2）孔隙网络内水相侵入过程包括惯性力作用、惯性力-黏滞阻力作用和黏滞阻力作用3个阶段。致密砂岩水相侵入的惯性力和惯性力-黏滞阻力作用阶段持续时间短，黏滞阻力对水相自吸侵入起主导作用。

3）将孔隙网络视为等效毛细管束，建立了致密砂岩液相侵入微观动力学模型，通过致密砂岩垂向自吸和负压差条件水平向侵入实验验证了模型的可行性。分析表明，负压差条件下水相仍能侵入岩心，且岩石越致密最大侵入深度越大。