裂缝型多孔介质气液两相渗流的分形分叉网络模型.pdf

Advances in Porous Flow 渗流力学进展渗流力学进展,2022,12(2),11-20 Published Online June 2022 in Hans.http:/www.hanspub.org/journal/apf https:/doi.org/10.12677/apf.2022.122002 文章引用文章引用:范子晨,王春玲,高妍,孟谨.裂缝型多孔介质气液两相渗流的分形分叉网络模型J.渗流力学进展,2022,12(2):11-20.DOI:10.12677/apf.2022.122002 裂缝型多孔介质气液两相渗流的分形分叉网络裂缝型多孔介质气液两相渗流的分形分叉网络模型模型 范子晨范子晨，王春玲王春玲，高高 妍妍，孟，孟 谨谨 中国计量大学理学院，浙江 杭州 收稿日期：2023年2月21日；录用日期：2023年3月21日；发布日期：2023年3月30日 摘摘 要要 裂缝型多孔介质的气液两相渗流对于页岩气水平井的产量预测具有重要的理论意义，因此，为了深入理裂缝型多孔介质的气液两相渗流对于页岩气水平井的产量预测具有重要的理论意义，因此，为了深入理解裂缝型多孔介质的气液两相渗流机理，本文采用分形分叉网络定量表征裂缝系统拓扑结构，发展了气解裂缝型多孔介质的气液两相渗流机理，本文采用分形分叉网络定量表征裂缝系统拓扑结构，发展了气液两相渗流的双重介质模型，推导了裂缝型多孔介质的有效渗透率和相对渗透率，并基于有限元方法数液两相渗流的双重介质模型，推导了裂缝型多孔介质的有效渗透率和相对渗透率，并基于有限元方法数值模拟了裂缝型多孔介质的气液两相渗流。结果表明，分形双重介质模型的预测结果与数值模拟以及实值模拟了裂缝型多孔介质的气液两相渗流。结果表明，分形双重介质模型的预测结果与数值模拟以及实验数据吻合较好，气相和液相相对渗透率随着液相饱和度的增加分别降低和提高，孔隙率和分形维数对验数据吻合较好，气相和液相相对渗透率随着液相饱和度的增加分别降低和提高，孔隙率和分形维数对相对渗透率具有显著影响。相对渗透率具有显著影响。关键词关键词 裂缝型多孔介质，页岩气，双裂缝型多孔介质，页岩气，双重介质，分形分叉网络，相对渗透率重介质，分形分叉网络，相对渗透率 A Fractal Bifurcating Network Model for Gas-Liquid Flow in Fractured Porous Media Zichen Fan,Chunling Wang,Yan Gao,Jin Meng College of Science,China Jiliang University,Hangzhou Zhejiang Received:Feb.21st,2023;accepted:Mar.21st,2023;published:Mar.30th,2023 Abstract The gas-liquid two-phase permeability in fractured porous media has important theoretical signi-ficance for the production prediction of shale gas horizontal wells.Therefore,in order to further understand the gas-liquid two-phase flow mechanism through fractured porous media,the fractal bifurcation network is used to quantitatively characterize the fracture topology.A fractal dual-medium model for gas-liquid two-phase flow was therefore developed to derive the analyti-范子晨 等 DOI:10.12677/apf.2022.122002 12 渗流力学进展 cal expressions of effective permeability and relative permeability of fractured porous media.And the gas-liquid two-phase flow in fractured porous media is also simulated based on the finite ele-ment method.The predicted results by the present fractal dual-medium model are in good agree-ment with the numerical simulation and experimental data.The results show that the relative permeability of gas phase and liquid phase decreases and increases with the increase of liquid phase saturation,respectively.The porosity and fractal dimension indicate significant effect on the relative permeability of gas phase.Keywords Fractured Porous Media,Shale Gas,Dual Medium,Fractal Bifurcating Network,Relative Permeability Copyright 2022 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/1.引言引言 页岩气作为重要的能源资源，对于国家工业、经济社会等的发展具有重要影响和战略意义。目前，水平井分段压裂是页岩储层成功改造的关键技术，国内外学者针对压裂后的页岩储层建立了大量的数学模型，描述页岩气的多重输运物理机制，分析页岩储层压裂后的产气规律等1 2。然而，多数页岩气井都具有较长的产水期，但大部分数学模型忽略了返排阶段压裂液对产气规律的影响，造成产气量预测结果偏大。不仅如此，压裂后的页岩储层包含天然裂缝和人工裂缝以及基质孔隙，介质类型复杂，孔缝尺度跨越几个量级。因此，裂缝型多孔介质的多相渗流模型是一项非常具有挑战性的课题，对于理解渗流机理和压裂后产气量预测具有重要的理论和实践价值。实验方法可以直接获得裂缝型多孔介质多相渗流的相对渗透率，但一般需要较长的时间和较高的成本，且对于一些低渗多孔介质往往很难直接测得相对渗透率，测量结果缺乏普适性3。随着计算机技术的发展，数值模拟方法也被应用到多孔介质气液多相渗流的研究中，但得出的结论通常含有经验参数，物理机理不清晰4。基于欧氏几何的传统分析方法一般仅适用于简单结构，很难构建实用的、符合客观规律的多孔介质多相渗流解析模型5。大量研究表明，多孔介质的裂缝网络、固相质量、孔隙面积、弯曲流道等均表现出统计自相似分形标度行为，分形几何为表征多孔介质的复杂微细结构提供了一种有效的数学工具和方法6 7 8。虽然多个学者提出不同的分形模型用以研究多孔介质气液两相渗流特性，但大多分形模型忽略了裂缝之间的交叉连通特征，而页岩储层中的压裂裂缝显示出分形分叉网络拓扑结构，扮演着页岩气的主要运移路径9。因此，本文采用分形分叉网络来定量表征裂缝型页岩储层的裂缝网络结构，其中分叉网络拓扑结构利用最小化流阻原理10进行结构优化，并将分形分叉网络嵌入各向同性均质多孔基质发展一种分形双重介质模型，推导裂缝型多孔介质气液两相渗流的相对渗透率，并采用有限元方法构建裂缝型多孔介质气液两相渗流的数学模型，研究裂缝网络中的局域流动特性和规律。2.气液两相渗流的分形双重介质模型气液两相渗流的分形双重介质模型 2.1.裂缝裂缝网络渗流模型网络渗流模型 考虑如图 1 所示的立方体代表性单元(REV)，采用分形分叉网络模型表征裂缝拓扑结构。假设分形分叉网络的数目具有统计自相似分形标度律11：Open AccessOpen Access范子晨 等 DOI:10.12677/apf.2022.122002 13 渗流力学进展 ()()00,max0fDN rrrr=(1)式中，0,maxr表示 0 级母管的最大半径，N 是 0 级母管半径大于0r的数目，fD是裂缝母管半径分形维数，其范围是02fD。假设方程(1)是连续的和可微的，0 级母管孔径从0r到00drr+的数目可以表示为：()10,max 00ddffDDfND rrr+=(2)因此，到第 0 级裂缝母管半径0r的概率密度函数为：()()100,minffDDff rD rr+=(3)裂缝母管半径分形维数fD与孔隙率之间存在如下标度关系11：()0,min0,maxEfDDrr=(4)式中下标 min 和 max 分别表示最小和最大孔径，ED为欧氏维数，此处2ED=。由于分形概率密度函数需满足归一化条件，因此裂缝网络 0 级母管的尺度满足()0,min0,max0fDrr=。实践中，通常将 20,min0,max10rr=作为满足分形标度律的条件。由此可得裂缝母管横截面积为：()0,max0,min220 max1d2rfprfDArNrD=，(5)则垂直于流动方向的横截面积可以表示为pAA=。Figure 1.Diagram of a fractal dual-medium model for fractured porous media 图图 1.裂缝型多孔介质分形双重介质模型示意图 典型的分形分叉网络可以通过迭代算法生成。通常单根分叉网络是由一母管按照一定的分叉角度分叉出两个或多个子管，然后每个子管按照相同的规则分叉出更多的分叉管道，不断重复这一分叉过程形成的，假设分叉网络的母管长度0l相同，且每一级分叉管道均是圆形直管且管壁厚度可以忽略不计。图2所示的是一个典型的双分叉分形网络。对于其中任意一级分叉而言，分别用kd和kl表示第k(0,1,2,m)级分叉管道的直径和长度，1kd+和1kl+表示第1k+级分叉管道的直径和长度，分叉角为，最大分叉级数为 m，分叉数为kn(图 2 中2kn=)。引入两个描述分叉网络几何结构的尺度因子来表示相邻两级分叉的长度比和直径比：范子晨 等 DOI:10.12677/apf.2022.122002 14 渗流力学进展 1kkkll+=(6)1kkkdd+=(7)由于网络具有分形特征，要求分叉数和分叉比保持常数(与分叉级数无关)10，即knn=、k=、k=。Figure 2.The kth single bifurcation of a fractal bifurcating network 图图 2.分形分叉网络的第 k 级单分叉结构 许多研究表明，当两种不相混溶流体同时流过多孔介质时，它们往往会有不同的流动路径，气相流体会占据大孔隙，而液相流体会占据小孔径3 12。为了区别气相和液相流体，假设存在一个临界半径0,cr，母管半径小于该临界半径的所有分叉网络都被液体占据，而母管半径大于该临界半径的所有分叉网络都被气体占据13。考虑一个单元横截面，则液体的饱和度可以表示为：()()0,0,min0,max0,min22200,0,min2220,max0,min0ddcffffrDDrcwrDDrrNrrSrrrN=(8)由此可以得到临界半径0,cr的表达式：()()1 20,0,max1fDcwrSr=+(9)根据 Hagen-Poiseuille 方程，单根分形分叉网络中气液两相流的流量可以表示为：()440,140181ggmgrnqPln+=(10a)()440,140181wwmwrnqPln+=(10b)对方程(10)进行积分，结合方程(9)得到裂缝网络气液两相总流量为：()()()0,max4440,max21401d11841ffcDrfDggwmrgfD rPnQqNSlDn+=+(11a)()()()0,min4440,max21401d1841fcfDrfDwwwmrwfD rPnQqNSlDn+=+(11b)在上述推导中，使用了0,min0,max1rr和()40,min0,max0fDrr=。运用广义达西定律：范子晨 等 DOI:10.12677/apf.2022.122002 15 渗流力学进展 LQKPA=(12)其中特征长度 L 定义如下：101cos1mLl+=+(13)联立方程(5)和(11)(13)，得到裂缝分叉网络气液两相流的渗透率：()()()()()()24410,max2,142111cos1118 141ffDmfwDn gwmfrDSnKSDn+=+(14a)()()()()()24410,max2,14211cos118 141ffDmfwDn wwmfrDSnKSDn+=+(14b)对于饱和多孔介质，每根分叉管道里只充满一种流体，因此，饱和裂缝网络的绝对渗透率可以表示为：()()()()2410,max14211cos18 141mfnmfrDnKDn+=+(15)联立方程(14)和(15)可以得到裂缝网络的相对渗透率：()()4,2,111ffDn gDn rgwwnKKSSK=+(16a)()4,2,1ffDn wDn rwwwnKKSSK=+(16b)2.2.双重介质双重介质渗流模型渗流模型 将各向同性多孔材料作为基质，利用分形分叉网络嵌入到基质材料中形成双重介质模型，用以研究裂缝型多孔介质的气液两相渗流特性。通常分叉裂缝起主要的渗流作用，而大量的流体储存在低渗透的母体基质。Journel 等人14提出有效渗透率eK可以表示成局域渗透率的幂平均：()111deVKkk xVV=(17)其中 V 是平均体积，幂指数符合11，其值取决于材料的性质，一般可以由实验数据拟合得到。实际上，多孔介质中的渗流管流化越明显，低渗流区域越少，有效渗透率越接近算术平均值(1=)。因此，为了简化计算，假定有效渗透率可以近似为网络和母体基质渗透率的算术平均值。对于双重多孔介质气液两相有效渗透率，方程(17)变为：1,e gmm gnn gKf Kf K=+(18a)1,e wmm wnn wKf Kf K=+(18b)其中,m gK、,m wK分别表示基质多孔材料的气相、液相渗透率；,n gK、,n wK分别表示裂缝网络的气相、液相渗透率，其值可通过方程(14)计算得出；mf和nf是母体基质和分叉网络的体积分数(1mnff+=)。对于总分叉级数为 m 的单根分叉网络，网络体积为：范子晨 等 DOI:10.12677/apf.2022.122002 16 渗流力学进展 ()12221002011mmkkkknVnr lr ln +=(19)对母管半径在0,min0,max,rr范围内积分可得分叉网络的总体积nV，与渗流区域的总体积 V 相比即可得到裂缝网络的体积分数：()()12220 0,max0,min20,max1112fmDfnfnD l rrfrnVD +=(20)相应地，母体基质材料的体积分数1mnff=。3.有限有限元模拟元模拟 由于裂缝在渗流等流体输运过程中通常扮演流体流动的首选和主要通道，对于流体渗流过程起主导作用，母体基质多孔材料的渗透率通常较小，故母体基质的渗透率远远小于裂缝网络的渗透率9。因此，本文使用 COMSOL Multiphysics 建立树状裂缝的气液两相流动模型。通过最小化流阻原理得到最优化结构参数10，即1 32=和37.45，设置第 0 级管道的长度和半径分别为05 mml=和00.5 mmr=。采用以下边界条件和初始条件：(1)选择甲烷作为气相，水作为液相。网络在初始条件下充满气态甲烷，左侧被液体置换，流体从右侧流出；(2)模型左侧为压力入口，边界符合狄利克雷边界条件的要求，压强为 10 Pa，流体禁止回流，无切向速度；(3)模型右侧为出口，并设定出口压强为 0 Pa，允许流体为气液混合，允许切向速度；(4)模型管壁为封闭边界，与外围无流体交换。通过 Navier-Stokes 方程结合水平集方法构建气液两相渗流数学模型，其动量守恒方程如下：()()()Tp=+uuIuuF (21)0=u (22)其中为流体密度，p 为速度矢量，u 为 t 时刻的速度分量，I 为单位对角矩阵，为粘度，F 为单位体积流体受到的外力。水平集方法是界面跟踪和形状建模中常用的一种数值模拟方法，其基本思想是将流体流动界面抽象为一个高维空间，通常用水平集函数的零水平集表示。通过界面水平集方程的不断迭代演化，直至稳定，得到流动界面的形状。在裂缝网络气液两相流模拟过程中，水平集函数可以描述气液两相界面形状的动态变化特征，其界面变化的方程可写成：()1lst+=u (23)流体的性质控制方程如下：()121=+(24)()121=+(25)其中，为水平集变量，一般取值为01，0=表示气相，1=表示液相，通常将取值为 0.5 的等范子晨 等 DOI:10.12677/apf.2022.122002 17 渗流力学进展 值面作为相界面；ls为界面厚度控制参数；是水平集函数初始化参数，一般为流体流动的最大速度；310.657 kg m=和321000 kg m=分别为气液两相的密度；511.1067 10Pa s=和322.98 10Pa s=分别为气液两相的粘度。4.结果与讨论结果与讨论 4.1.结结果验证果验证 在验证的过程中，将最小最大孔径比设定为0.001=，孔隙率为0.2=。根据方程(4)可以得到分形维数1.77fD=。数值计算结果表明液体从左侧入口进入分叉管道，沿着管道对气相进行驱替，其体积分数如表 1 所示。Table 1.The volume fraction of fluid at different times in a fracture network 表表 1.裂缝网络不同时刻下流体体积分数 t(s)0.01 0.10 0.30 0.50 液相 体积 分数 根据简化后的达西方程计算渗透率：oDLQLKvPAP=(26)其中oQ为出口流体的流速，Dv为达西速度。将数值模拟计算得到的气液两相相对渗透率与方程(16)的分形模型进行对比分析，如图 3 所示。进一步地，将分形模型预测的相对渗透率与实验数据(2%=,5%=)15 16对比，如图 4 所示。在很低的液体饱和度下，液相相对渗透率接近于零，而当饱和度达到 1.0 时，液相相对渗透率则接近于 1.0；气相的相对渗透率则与液相相对渗透率相反。由图 3 和图 4 可以发现，分形模型的预测结果与有限元模拟和实验数据较为吻合，验证了分形模型的合理性和正确性。0.00.20.40.60.81.00.00.20.40.60.81.00.00.20.40.60.81.0Df=1.77SwKn,rwKn,rg Kn,rg Kn,rw Fractal Model Numerical Simulation Figure 3.A comparison of relative permeability by fractal model and finite element simulation 图图 3.分形模型预测值与有限元模拟对比图 范子晨 等 DOI:10.12677/apf.2022.122002 18 渗流力学进展 0.00.20.40.60.81.00.00.20.40.60.81.00.00.20.40.60.81.0Kn,rwKn,rgSw(a)Df=1.43 Kn,rg Kn,rw Fractal model Experiment(=2%)150.00.20.40.60.81.00.00.20.40.60.81.00.00.20.40.60.81.0SwKn,rwKn,rg(b)Df=1.57 Kn,rg Kn,rw Fractal model Experiment(=5%)16 Figure 4.A comparison between predicted relative permeability by fractal model and experimental data with porosity of(a)2%=and(b)5%=图图 4.分形模型预测结果与实验数据对比图：(a)2%=；(b)5%=4.2.参数化分析参数化分析 0.00.10.20.30.40.50.60.70.80.91.01.01.11.21.31.41.51.61.71.81.92.0Df=0.01=0.005=0.001 Figure 5.The relationship between fractal dimension and porosity as well as the ratio of maximum to minimum pore size 图图 5.分形维数与孔隙率以及最大最小孔径比之间的关系 0.00.20.40.60.81.00.00.20.40.60.81.00.00.20.40.60.81.0SwKn,rgKn,rw(a)=0.001 Kn,rg Kn,rw Df=1.1,=0.002 Df=1.3,=0.008 Df=1.5,=0.032 Df=1.7,=0.126 Df=1.9,=0.5010.00.20.40.60.81.00.00.20.40.60.81.00.00.20.40.60.81.0Kn,rgKn,rwSw Kn,rg Kn,rw =0.001,Df=1.67 =0.002,Df=1.63 =0.005,Df=1.57 =0.010,Df=1.50 (b)=0.10 Figure 6.The effect of porosity and fractal dimension on relative permeability 图图 6.孔隙率和分形维数对相对渗透率的影响规律 范子晨 等 DOI:10.12677/apf.2022.122002 19 渗流力学进展 根据方程(4)可知孔隙分形维数fD、孔隙率以及孔径比之间存在内在关联，如图 5 所示，在固定的孔径比，分形维数随孔隙率的增大而快速变大，当孔隙率趋向于 1.0 时，此时孔隙几乎填满整个平面，因此孔隙分形维数趋向于 2.0；在相同孔隙率条件下，分形维数随孔径比的减小(孔隙范围变大)而增大，孔隙结构越复杂。如图 6(a)所示，在特定的液体饱和度下，分形维数的增加意味着孔隙率的提高和孔隙尺度的增加，而气相占据较大尺寸孔隙，因此气相相对渗透率随着分形维数的增加而提高，而液相相对渗透率随着孔隙分形维数的增加而显著降低。而图 6(b)显示，在相同孔隙率和液相饱和度条件下，分形维数对于气相相对渗透率的影响较小，而液相相对渗透率随着孔隙分形维数的增加而降低，但幅度也较小。由此可见，裂缝网络的结构对于气液两相渗流具有显著的影响。5.结结论论 本文基于分形几何理论建立了裂缝型多孔介质气液两相渗流模型，得到了气液两相相对渗透率的解析表达式，并利用有限元方法建立了裂缝型多孔介质的气液两相渗流的数学模型。结果表明，裂缝型多孔介质气液两相渗流的分形模型预测结果与有限元模拟以及实验数据吻合较好，相对渗透率主要依赖于裂缝结构特征和气液两相分布特征。气相相对渗透率和液相相对渗透率随着液相饱和度的增加而降低和提高，且在特定液体饱和度的情况下，增大裂缝孔隙分形维数可以提高气相相对渗透率和降低液相相对渗透率。相比较而言，孔隙率对于相对渗透率的影响更为显著，但在相同孔隙率条件下，孔隙微细结构的差别(分形维数不同)仍然对于气液渗流产生较为明显的影响。本文提出的分形双重介质模型为裂缝型多孔介质渗流机理研究以及页岩气高效开采等领域提供了一种有效的数学方法和理论依据。参考文献参考文献 1 Li,Y.M.,Jiang,Y.S.,Zhao,J.Z.,et al.(2015)Extended Finite Element Method for Analysis of Multi-Scale Flow in Fractured Shale Gas Reservoirs.Environmental Earth Sciences,73,6035-6045.https:/doi.org/10.1007/s12665-015-4367-x 2 赵金洲,李志强,胡永全.考虑页岩储层微观渗流的压裂产能数值模拟J.天然气工业,2015,35(6):53-58.3 Abaci,S.and Edwards,J.S.(1992)Relative Permeability Measurements for Two Phase Flow in Unconsolidated Sands.Mine Water and the Environment,11,11-26.https:/doi.org/10.1007/BF02919583 4 Andres,R.V.and Bernardo,M.R.(2020)Uncertainty Quantification and Sensitivity Analysis for Relative Permeability Models of Two-Phase Flow in Porous Media.Journal of Petroleum Science and Engineering,192,Article ID:107297.https:/doi.org/10.1016/j.petrol.2020.107297 5 Snow,D.T.(1969)Anisotropic Permeability of Fractured Media.Water Resources Research,5,1273-1289.https:/doi.org/10.1029/WR005i006p01273 6 郁伯铭.多孔介质输运性质的分形分析研究进展J.力学进展,2003,33(3):333-346.7 Xu,P.(2015)A Discussion on Fractal Models for Transport Physics of Porous Media.Fractals,23,Article ID:1530001.https:/doi.org/10.1142/S0218348X15300019 8 Qiu,S.X.,Yang,M.and Xu,P.(2020)A New Fractal Model for Porous Media Based on Low-Field Nuclear Magnetic Resonance.Journal of Hydrology,586,Article ID:124890.https:/doi.org/10.1016/j.jhydrol.2020.124890 9 Luo,Y.F.,Xia,B.W.and Li,H.L.(2021)Fractal Permeability Model for Dual-Porosity Media Embedded with Natural Tortuous Fractures.Fuel,295,Article ID:120610.https:/doi.org/10.1016/j.fuel.2021.120610 10 徐鹏.树状分形分叉网络的输运特性D:博士学位论文.武汉:华中科技大学,2008.11 Yu,B.M.and Cheng,P.(2002)A Fractal Permeability Model for Bi-Dispersed Porous Media.International Journal of Heat and Mass Transfer,45,2983-2993.https:/doi.org/10.1016/S0017-9310(02)00014-5 12 Burdine,N.T.,Magnolia,P.and Dallas(1953)Relative Permeability Calculations from Pore Size Distribution Data.Petroleum Transactions,198,71-78.https:/doi.org/10.2118/225-G 13 Xu,P.,Qiu,S.X.and Yu,B.M.(2013)Prediction of Relative Permeability in Unsaturated Porous Media with a Fractal Approach.International Journal of Heat and Mass Transfer,64,829-837.https:/doi.org/10.1016/j.ijheatmasstransfer.2013.05.003 范子晨 等 DOI:10.12677/apf.2022.122002 20 渗流力学进展 14 Journel,A.G.,Deutsch,C.V.and Desbarats,A.J.(1986)Power Averaging for Block Effective Permeability.Society of Petroleum Engineers,Los Angeles,329-334.https:/doi.org/10.2118/15128-MS 15 Kenton,A.R.,Wooyong,U.and Sean,M.C.(2019)Relative Permeability for Water and Gas through Fractures in Ce-ment.PLOS ONE,14,e0210741.https:/doi.org/10.1371/journal.pone.0210741 16 郭晓茜.煤储层相对渗透率动态变化规律研究D:博士学位论文.北京:中国地质大学,2016.