基于页岩气井返排特征的闷井时间优化方法

翁定为，江 昀，易新斌，何春明，车明光，朱怡晖

（1.中国石油勘探开发研究院，北京 100083；2.中国石油油气和新能源分公司，北京 100028；3.中国石油西南油气田分公司页岩气研究院，四川成都 610051）

页岩气作为一种重要的非常规油气资源，在美国已实现商业开发，2021 年美国页岩气产量达到7 924.47×108m3[1]。2007 年至今，通过不断积累页岩气开发经验，我国基本实现了页岩气资源效益开发，产量呈现高速增长态势，2022 年页岩气产量达230×108m3[2-3]。页岩气资源高效开发必须依靠大规模水力压裂技术[4-6]，多位学者研究发现[7-12]，压裂后关井一段时间，增大滞留液量，有利于促进渗吸过程，提高压后裂缝网络复杂程度，最终提高油气井产量。E.Ghanbari 等人[8，13]研究发现，延长闷井时间有利于提高早期产量，但闷井时间过长不利于后期产量，这是由于基质渗透率太低或裂缝复杂程度不够，导致渗吸速率太慢，压裂液在近井地带滞留（即“水锁”），裂缝附近油/气相对渗透率降低，油/气排出阻力增大[14]。因此，合适的闷井时间对于充分发挥渗吸置换和降低“水锁”伤害至关重要。目前，确定闷井时间方法大致分为以下3 种：第1 种是经验公式法[15]。该方法是根据特定井或特定储层特征，构建计算闷井时间经验公式，但其不具有普遍适用性。第2 种是解析法[16-19]。该方法是通过求解毛细管内活塞式驱替模型，得到无因次时间。第3 种是数值法。该方法是通过采用数值法求解非活塞式驱替模型，得到产能与闷井时间的关系[20-23]，但未考虑闷井过程中动态相渗规律。上述3 种方法中的第2 种方法虽然存在一定不足，但应用最简便。

页岩气井返排资料显示，返排液的矿化度显著高于注入压裂液的矿化度[24-25]。基于这一特征，A.Zolfaghari 等人[26]研究了返排液盐度与累计产液量的关系，提出了基于离子扩散理论的缝宽分布反演方法。Yang Liu 等人[27]进一步提出了基于裂缝-基质三维渗吸模型的返排率模型，推导出无因次时间与返排率的函数关系。笔者选用泸州地区龙马溪组深层页岩气井和威远地区龙马溪组中深层页岩气井页岩岩心，开展自发渗吸试验，确定渗吸量达到稳定阶段的时间（即实验室尺度闷井时间）；统计返排液矿化度和返排率数据，利用A.Zolfaghari 等人[26-27]理论模型反演计算裂缝宽度和特征长度；结合实验室尺度闷井时间、无因次时间模型及矿场尺度特征长度等数据，形成了基于返排液相关数据的闷井时间优化方法。

1 自发渗吸试验

选取四川盆地泸州地区深层页岩泸A 井（垂深3 986.12～3 989.18 m）和威远地区中深层页岩威B 井（垂深2 705.18～2 710.13 m）岩样进行自发渗吸试验，以确定渗吸量达到稳定阶段的时间，即实验室尺度闷井时间。

1.1 试验准备

岩样基础物性参数测试结果表明，静态容量法（测试介质为氮气）测定的孔隙度为5.45%～8.41%，脉冲衰减法（测试介质为氮气）测定的渗透率为0.012～0.145 mD（见表1）。

表1 页岩样品物性参数Table 1 Physical parameters of shale samples

测试流体为蒸馏水，其室温下的密度为1.0 g/cm3，黏度为1.0 mPa·s。

1.2 试验方法

1）对岩样进行烘干（105 ℃密闭烘箱烘干48 h）处理，冷却后取出，测量干燥岩样的质量；

2）采用MiniMR-VTP 低场核磁共振分析仪（磁场强度0.5 T，磁体温度32 ℃，回波时间0.3 ms，间隔时间3 000 ms，回波数量8 000 个），测试干燥岩样的T2谱；

3）用蒸馏水清洗烧杯和电导率仪的电极，在清洗后的烧杯内加入蒸馏水200 mL，测试蒸馏水的电导率，若电导率大于2 μS/cm，则再次清洗烧杯，直至电导率测试值符合要求；

4）将干燥岩样放入装有蒸馏水的烧杯内，使用保鲜膜封闭烧杯杯口，以减少水分蒸发，降低试验误差；间隔一段时间后取出岩样，擦拭其表面液体，使用ME204E 梅特勒分析天平（精度为0.000 1 g，量程为220 g）称其质量，使用核磁共振分析仪测试岩样的T2谱，同时使用电导率仪测量烧杯内液体的电导率，记录试验数据，将岩样放回烧杯中，继续进行自发渗吸试验；

5）间隔一段时间后取出岩样，重复试验步骤4），直至测试数次后岩样的质量和T2谱不再出现显著变化时，结束试验。

1.3 试验结果与讨论

为对比深层与中深层页岩的渗吸-离子扩散规律，开展自发渗吸-离子扩散同步试验，渗吸过程中的核磁共振T2谱如图1 所示。由图1 可以看出，T2谱均为双峰形态，分别对应弛豫时间0.1～10 和10～100 ms。Jiang Yun 等人[28]认为，弛豫时间介于0.1～100 ms 时对应岩心中纳米孔（0.01～1 μm），弛豫时间大于100 ms 时对应微孔和大孔（大于1 μm）。因此，渗吸过程吸入的水主要赋存于纳米孔内。渗吸过程中，岩样A1 和A2 的T2谱的峰值右移，说明基质中的孔隙直径增大，这可能与黏土矿物吸水膨胀及盐类溶解有关。岩样B1 和B2 的T2谱形状基本不变，这说明岩心内部的孔隙结构未出现明显变化。此外，不同时间测定的T2谱结果显示，干燥岩样的信号幅值最小；随着时间增长，信号幅值上升，之后基本不变，说明岩样中的含水量基本不变，达到饱和状态。

图1 岩样自发渗吸试验T2 谱测试结果Fig.1 Results of T2 spectrum for rock samples spontaneous imbibition

对比渗吸试验前后不同弛豫时间的信号分布频率可以看出，干燥状态下的氢质子信号很弱，主要来自黏土水或含结晶水的矿物（见图2）。渗吸之后，4 块岩样的T2主要分布在0.1～100.0 ms，占比达98.5%以上。其中，威B 井2 块岩样的频率分布基本一致，T2分布在0.1～1.0 ms 的占比在65%左右，明显高于泸A 井2 块岩样，说明吸入的水主要分布在纳米微孔中。泸A 井2 块岩样渗吸试验后氢质子信号分布特征不同，岩样A1 吸入的水主要分布在纳米中孔（T2分布在1.0～10.0 ms 的占比约70%），岩样A2 吸入的水主要分布在纳米微孔和纳米中孔（T2分别为0.1～1.0 和1.0～10.0 ms，占比分别为48.3%和43.1%）。

试验过程中的吸水体积（根据称重法测定的质量换算得到）、核磁共振T2谱面积（不同时刻测定的T2谱面积与干燥岩样T2谱面积差值）及电导率随时间平方根的变化结果如图3 所示。

图3 不同岩样的渗吸体积、T2 谱面积及电导率随时间变化曲线Fig.3 Changes of imbibition volume,T2 spectrum area and conductivity with time of different rock samples

从图3 可以看出，渗吸体积与T2谱面积的变化趋势基本一致，说明2 种测试方法的差异较小。渗吸初期，岩样的吸水体积与时间平方根呈正比关系，与L.L.Handy 等人[19]的研究结果一致。之后，曲线出现拐点（岩样A1，A2，B1 和B2 自发渗吸试验中拐点处对应时间分别为22，29，3 和4 h），曲线与时间轴接近平行，岩样已经没有明显的吸水特征，这与T2谱测试结果一致；电导率与时间平方根同样呈正比关系，且整个渗吸过程中保持线性相关性。

2 矿场尺度特征长度求解

页岩气井特殊的“产液-产盐”返排动态曲线可为评价体积压裂缝网形态提供丰富的信息。利用返排数据进行压后裂缝参数反演，可以为确定矿场尺度特征长度提供依据。以Yang Liu 等人[27]的模型为基础，根据无因次时间与返排率关系，结合A.Zolfaghari 等人[26]的模型，反演得到裂缝宽度，求解矿场尺度特征长度。

2.1 理论模型

根据Yang Liu 等人[27]提出的基于裂缝-基质三维渗吸理论的返排率预测模型，将大规模体积压裂后“打碎”的储层简化为带有裂缝的立方体（见图4），并在岩石裂缝中铺设支撑剂。

图4 裂缝-基质示意图Fig.4 Illustration for fracture-matrix

根据质量守恒原理，注入地层压裂液的体积等于水力裂缝的体积（忽略流体损失），通过压裂作业形成m 个基质块：

式中：m为压裂形成的基质块数量；Vinj为注入地层压裂液的体积，m3；a为基质块的长度，m；d为支撑剂颗粒的直径，m；n为裂缝中支撑剂铺设层数。

吸入岩石基质块压裂液与注入地层压裂液的体积之比为：

式中：Vimb为吸入岩石基质块压裂液的体积；x为水相渗吸距离，m；LD为无量纲长度标度；Swf为最终含水饱和度；Swi为初始含水饱和度；ϕ为孔隙度。

Yang Liu 等人[27]经过推导，得到的返排率计算公式为：

式中：Rre为压裂液返排率；tD为无因次时间；LC为特征长度，与基质块尺寸和边界条件有关；Vb为岩心基质体积，cm3；Aj为j方向上渗吸接触面的面积，cm2；LAj为j方向上渗吸前缘沿开启面到封闭边界的距离，cm；K为渗透率，mD；μw为黏度，mPa·s。

根据式（6），计算得到基质块长度与特征长度的关系为：

由式（3）—（7）可知，可根据自发渗吸试验结果求取无因次时间，可根据现场返排液数据求得返排率，因此，只需确定裂缝宽度，即可确定特征长度。

根据A.Zolfaghari 等[26]提出的多级树状网络裂缝模型确定裂缝宽度。假定气藏为均质无限大，且具有各向同性；流体为牛顿流体，流动满足达西定律；流体中盐的运移是由基质到裂缝，储层中流体的矿化度远大于压裂液的矿化度。利用Fick 扩散定律来表示盐离子由基质到裂缝的流动方程：

式中：Ji为第i条裂缝的离子扩散通量，mol/s；Af,i为基质和第i条裂缝之间的界面面积，m2；Cm为基质内盐的质量浓度，mg/L；Cf,i为第i条裂缝内盐的质量浓度，mg/L；Lm为人工裂缝与基质中某一点的距离，反映盐质量浓度梯度变化特征，m；D为扩散系数，m2/s。

第i条裂缝内的平均盐质量浓度为：

式中：Δt为裂缝面与压裂液接触的时间，s；Wf,i为第i条裂缝宽度，m。

基于物质平衡方程，A.Zolfaghari 等人[26]推导得到裂缝宽度与返排液中盐质量浓度的函数关系式：

式中：f（Wf）为裂缝宽度分布，mm；NP,w为归一化的返排率；Vf,i为第i条裂缝体积，m3；Qw为返排液累计返排量，m3。

现场应用中，基质内的盐质量浓度根据现场返排液总矿化度数据得到，扩散系数根据返排液离子类型确定，Lm取10 倍Wf[26]，归一化返排率根据压裂液注入体积和返排后收集的压裂液体积确定。

2.2 理论模型计算实例

统计分析泸A 井和威B 井的产气量、产水量和套压等基础资料，结果如图5 所示。其中，泸A 井完成压裂22 段，压裂水平段长1 436 m，累计注入地层压裂液42 323 m3，见气前返排率9.95%；威B 完成压裂18 段，压裂水平段长1 655 m，累计注入地层压裂液45 440 m3，见气前返排率1.49%。

统计分析压裂后的返排液矿化度和返排液离子类型，结果如图6 所示。从图6 可以看出，随着时间的推移，返排液盐质量浓度逐渐升高，其中，钠离子与氯离子是主要的矿化度来源，其他离子含量较少。基质内盐质量浓度（Cm）应取总矿化度达到稳定阶段对应的值，威B 井返排时间（图6（b））达到100 d 左右，总矿化度接近不变，Cm取35 000 mg/L；泸A 井实测返排液矿化度时间为60 d，因此，为了与威B 井进行对比，拟合总矿化度曲线，并求解得到返排时间为100 d 对应的总矿化度，得到Cm为20 000 mg/L。

盐离子主要为氯离子和钠离子，扩散系数均为1.484×10-9m2/s，推到出泸A 井和威B 井第i条裂缝内的平均盐质量浓度（Cf,i）与裂缝宽度（Wf,i）的关系式：

根据返排数据和盐质量浓度变化，计算得到Wf对应的f（Wf）；再分别对各个盐质量浓度对应缝宽范围内的f（Wf）求和，得到该缝宽范围内的裂缝体积分数的分布情况。据此，可以得到泸A 和威B 井不同缝宽的裂缝体积分数分布结果（见图7）。

根据不同裂缝宽度的裂缝体积分数计算结果（见图7），可知威B 井的产盐量更高，计算得出的裂缝宽度主要分布在1.0～2.0 mm，而泸A 井裂缝宽度主要分布在2.0～2.5 mm。并且，泸A 井和威B 井加权平均缝宽分别为1.73 mm 和1.30 mm。将泸A 井和威B 井的岩样测定的实验室尺度闷井时间、现场返排率测试结果及缝宽分布计算结果等参数代入式（3），计算得到基质块长度a（即等效裂缝间距）的概率分布（见图8）。结果表明，泸A 井的基质块加权平均长度为0.20 m，根据式（6），得到对应的特征长度LC为0.10 m；威B 井的基质块加权平均长度为0.32 m，特征长度LC为0.16 m。

图8 基质块长度分布Fig.8 Distribution of matrix length

3 实例计算

以返排率预测模型（式（3））为基础，结合自发渗吸试验结果、返排液矿化度测试结果以及裂缝宽度反演结果，建立了矿场尺度闷井时间计算流程。详细步骤如下：

1）开展自发渗吸室内试验，确定岩心尺度闷井时间；

2）将孔隙度、渗透率、表面张力和岩心尺度闷井时间等参数代入式（3），计算无因次时间；

3）每隔一段时间收集并测定返排液总矿化度、离子类型、返排液总量和返排率；

4）根据式（10），计算裂缝宽度分布；

5）将步骤2）—4）计算得到的无因次时间、返排率、裂缝宽度等参数代入式（3），计算矿场尺度基质块长度；

6）根据式（6）计算矿场尺度特征长度，并代入式（5）计算矿场尺度闷井时间。

根据上述步骤，确定了泸A 井岩样A1 的实验室尺度闷井临界时间为22 h，特征长度为0.10 m，计算得到泸A 井矿场尺度闷井时间为12.5 d；威B 井岩样B1 的实验室尺度闷井时间为3.0 h，特征长度为0.16 m，计算得到威B 井矿场尺度闷井时间为16.7 d。

由式（4）可知，矿场尺度闷井时间与岩心尺度闷井时间和特征长度的平方均呈正比。从矿化度测试结果看，这样的矿场尺度闷井时间计算结果似乎是矛盾的，因为威B 井的水力裂缝更窄，裂缝网络更复杂。为了验证威B 井水力裂缝网络的复杂程度，采用A.Zolfaghari 等人[29]提出的压后裂缝表面积计算方法，并结合渗吸-离子扩散试验结果（见图3），计算得到泸A 井和威B 井的裂缝表面积分别为59 744.6 m2和125 960.1 m2。结合图7 裂缝宽度计算结果（即泸A 井和威B 井加权平均缝宽分别为1.73 mm 和1.30 mm），可知，威B 井产盐量更高，水力裂缝更窄，裂缝表面积和裂缝体积更大。因此，认为此时对应的裂缝网络更密集，基质单元间的间距更小，缝网形态更复杂。这样会造成威B 井特征长度短，并且岩心尺度闷井时间短，最终导致矿场尺度闷井时间也短。因此，矿场尺度闷井时间不是简单的线性关系，需要综合考虑返排液矿化度、返排率和渗吸速度等因素影响，进一步研究其内在机理。上述矿场尺度闷井时间计算方法适用于储层离子扩散特征明显及返排液矿化度较高的页岩气储层，并且需要区块内部分井已取得返排资料；而对于新区块而言，该方法并不适用，建立综合考虑“产液-产盐”剖面、闷井压降曲线及支撑剂回流等因素的闷井时间优化方法，是未来研究方向之一。

4 结论与建议

1）页岩气井压后闷井是提高采收率的主要技术手段，自发渗吸是微观渗流的主要方式之一，开展自发渗吸试验，有助于深入理解矿场尺度闷井过程中的压裂液侵入过程。泸州深层页岩气井和威远中深层页岩气井岩心自发渗吸试验结果表明，渗吸初期，岩心吸水量与时间平方根正相关，吸水量达到稳定阶段对应的时间即为岩心尺度闷井时间。

2）页岩气井特殊的“产液-产盐”返排动态曲线为评价体积压裂缝网形态提供了丰富的信息。根据现场返排液中矿化度和返排率变化对比了泸州深层页岩气井和威远中深层页岩气井压后裂缝宽度分布特征，结果表明，中深层页岩气井返排液的矿化度更高，水力裂缝更窄，裂缝表面积和裂缝体积更大，水力裂缝网络更密集，缝网形态更复杂。

3）基于无因次时间模型确定矿场尺度闷井时间是一个常规研究方向，但矿场尺度特征长度的选取仍是难点。矿场尺度闷井时间与岩心尺度闷井时间并不一定呈正相关关系，其结果受渗吸速率、返排液矿化度和返排率等因素影响，需持续跟踪闷井过程中的压裂液滤失规律和支撑剂回流规律，进一步完善矿场尺度闷井时间优化方法。