基于粗糙离散裂隙网络模型的煤层气流动特性研究.pdf

29fracture networkmodelJSafettinCoalMines54(10):29-35.移动扫码阅读LI Wenbin,ZHANG Shujin,HENG Xianwei.Study on flow characteristics of CBM based on rough discrete全，2 0 2 3，10):29-35.李文彬，张离散裂隙网络模型的煤.煤矿安特件研饼SafetyinCoalMinesOct.20232023年10 月煤防发全No.10Vol1.54第54卷第10 期DOl:10.13347/ki.mkaq.2023.10.005基于粗糙离散裂隙网络模型的煤层气流动特性研究李文彬1.2，张书金1.2，衡献伟51,2,3（1.贵州省煤矿设计研究院有限公司，贵州贵阳550 0 2 5;2.贵州省矿山安全科学研究院有限公司，贵州贵阳550 0 2 5；3.贵州大学矿业学院，贵州贵阳550 0 2 5）摘要：煤体内随机分布的裂隙网络是煤层气的主要流动通道，而研究煤层气流动特性的关键在于准确模拟煤体中错综复杂的裂隙网络。基于像元概率分解理论，介绍了一种新的离散裂隙网络建模方法，相比传统方法多假设裂隙为平滑裂隙这一不足，该方法可充分考虑裂隙的粗糙特征，并结合数值模拟手段分析了不同裂隙特征参数对煤层渗透性的影响规律。结果表明：裂隙率（0.1 0.2 5）、分形维数（1.55 1.8 5）、最大裂隙长度（0.0 11 0.0 2 6 m）以及比例系数（0.0 0 6 0.0 8）与裂隙网络出口流速之间呈现较好的幂律增长关系，而迁曲度分形维数（1.0 1.0 2）、裂隙倾角（10 8 5)与裂隙网络出口流速之间呈现幂律降低的关系；其中，比例系数、裂隙倾角和最大裂隙长度对裂隙网络出口流速的影响最为显著。关键词：煤层气；裂隙网络；分形几何；流速；粗糙特征中图分类号：TD713文献标志码：A文章编号：10 0 3-496 X(2023)10-0029-07Study on flow characteristics of CBM based on rough discrete fracture network modelLI Wenbin-,ZHANG Shujin2,HENG Xianweil.2.3(1.Guizhou Coal Mine Design Research Institute Co.,Ltd.,Guiyang 550025,China;2.Guizhou Mine Safety Scientific ResearchInstitute Co.,Ltd.,Guiyang 550025,China;3.College of Mining Engineering,Guizhou University,Guiyang 550025,China)Abstract:The randomly distributed fracture network in coal body is the main flow channel of coalbed methane,and the key to studythe flow characteristics of coalbed methane is to accurately simulate the complex fracture network in coal body.Based on the pixelprobability decomposition theory,a new discrete fracture network modeling method is introduced.Compared with the traditionalmethod,which assumes that the fracture is a smooth fracture,this method can fully consider the rough characteristics of the fracture.Combined with numerical simulation,the influence of different fracture characteristic parameters on the permeability of the coalseam is analyzed.The results show that the porosity(0.1-0.25),fractal dimension(1.55-1.85),maximum fracture length(0.011-0.026 m)and proportion coefficient(0.006-0.08)have a good power-law growth relationship with the outlet velocity of the fracture network,while the tortuosity fractal dimension(1.0-1.02),fracture dip angle(10-85)and the outlet velocity of the fracture network have apower-law decrease relationship;among them,the proportion coefficient,fracture dip angle and maximum fracture length have themost significant effects on the outlet velocity of fracture network.收稿日期：2 0 2 2-0 6-2 4责任编辑：兰莹基金项目：贵州省科技支撑计划资助项目（黔科合支撑 2 0 2 1 一般348；黔科合支撑 2 0 2 2 一般0 0 5；黔科合支撑 2 0 2 0 2 Y030号）作者简介：李文彬（1990 一），男，河南焦作人，工程师，硕士，从事煤矿瓦斯灾害预测及防治、矿山压力与岩层控制工作。E-mail：通讯作者：张书金（198 5），男，河南南阳人，高级工程师，硕士，从事煤矿瓦斯灾害预测及防治工作。E-mail：32 8 8 3337 2 q.c o m30SafetyinCoal Mines2023年10 月Oct.2023Vol.54No.10煤码发全第10 期第54卷Key words:coalbed methane;fracture network;fractal geometry;velocity;rough feature煤层气是一种以甲烷为主要成分的清洁高效能源，对其开采不仅可以减少温室气体排放、接替常规天然气能源，同时也可减少煤矿灾害的发生，保障矿井安全生产等-2 。目前的研究认为 3-5,煤体是一种多孔介质材料，其内部广泛分布的裂隙结构是气体运移的重要通道，对产气量有显著影响。煤体内的裂隙分布呈现明显的随机性和无序性，现有模拟储层裂隙的方法主要分为等效连续介质模型 6 和离散裂隙网络模型 7 ，2 种方法各有优缺点，在实际应用时，等效连续介质模型对裂隙的假设过于理想化，无法描述真实裂隙的细节，而离散裂隙网络模型恰好弥补了这一不足，可以充分刻画裂隙的真实细节和分布特征，因而被广泛应用。其中，采用MonteCarlo方法生成离散裂隙网络模型是目前最多的方法，许多学者利用该方法分析了煤层气在裂隙内的流动规律。然而传统离散裂隙网络模型仍然是假设裂隙为平滑裂隙，对真实裂隙细节的描述仍然存在不足。而真实煤体内的裂隙绝大多数是凹凸不平的，气体在其内的流动也不同于平滑裂隙网络模型。因此，采用一种新的粗糙离散裂隙网络建模方法建立裂隙网络几何模型，该方法相比MonteCarlo方法可充分考虑裂隙粗糙特征，因而更加符合实际情况。为此，基于建立的裂隙网络模型，结合数值模拟分析煤层气在裂隙内的流动规律，探讨不同裂隙结构特征参数对煤层气流动特性的影响。1粗糙离散裂隙网络模型实际煤体内的裂隙往往是弯曲的，具有一定的粗糙特性，且裂隙相互交叉形成裂隙网络，组合形式多为网状和不规则网状，渗透率相对较好。为了更真实地建立符合实际工程情况的粗糙离散裂隙网络模型，WU等 8 和XIA等 9 等提出了基于像元概率分解算法建立粗糙离散裂隙网络模型的方法，根据其理论，在像元空间中，任意1条裂隙都可以由点集合的形式组成，即：Gp=(Di,Mi,ti),i=0,1,2,.,n(1)式中：G为裂隙点集合系统；D,为裂隙点的成长方向；M,为成长概率；t为成长步长。根据像元概率分解算法的基本理论，借助MATLAB软件编制二维粗糙离散裂隙网络生成程序，以实现像元空间范围内多条裂隙的随机生成。由式(1)可知，如果确定了裂隙的成长方向、成长概率、成长步长，裂隙点集合系统也就确定了。为了更好地说明其具体实现过程，假定裂隙的成长概率M1=0.67,0.33,0,0,0,0,0,0)，为二项随机，成长步长t=15，则1条在像元空间范围 10 0 10 0 中裂隙点集合，粗糙离散裂隙网络生成过程如图1。100100100909090808080707070606060505050404040303030202020101010010203040506070809010001020304050607080901000102030405060708090100XXX（a）单裂隙集合（b）150 条裂隙集合（c）30 0 条裂隙集合图1米粗糙离散裂隙网络生成过程Fig.1Generation process of rough discrete fracture network当需要生成具有同样特征的150 条裂隙点集合时，原理与生成1条裂隙一致，只需要在程序中设置裂隙总数，裂隙在像元空间中随机分布，结果如图1(b）。由于真实煤岩体中的裂隙相互交错，为生成更为符合真实情况的裂隙网络，这里裂隙的分布特征参考文献 10-11 获得的统计规律，即裂隙长度幂律符合分布特征，倾角分布符合Fisher分布特征。在图1(b)中增加150 条裂隙，成长概率M2=(0,0,0,0,0,0,0.33,0.67)，成长步长t=15,因此，具有不同成长概率，数量为30 0 条的二维像元空间中的裂隙网络集合系统如图1(c）。2粗糙离散裂隙网络渗透率模型煤体内的裂隙相互交叉形成裂隙网络，对于31SafetyinCoalMinesOct.20232023年10 月No.10Vol.54煤矿发全第10 期第54卷任意复杂的裂隙网络均是由一条条单裂隙所组成目前描述流体在单裂隙内的流动多采用经典的平行板立方定律 12 ，即：Q2(1-cos2 0;sin 2)Apq(l)=(2)12uLo式中：q为流过裂隙的总流量，m/s；为裂隙开度，m；1为裂隙长度，m；0 i、2 分别为裂隙方位角和倾角，（）；u为流体动力黏度，Pas;Ap为压差，Pa；Lo 为特征单元体长度，m。MIAO等!借鉴多孔介质模型中孔隙网络渗透率模型的推导方法，在上述平行板立方定律的基础上，结合分形结合推导得到了适用于随机分布的裂隙网络渗透率模型，即：B2(2-D)(1-cos*0 in2)Kf(3)2124-DI1-fmax式中：K为裂隙网络渗透率，m；为比例系数，取值范围一般为10-310 ；D为裂隙网络的分形维数；f为裂隙率；Imax为最大裂隙长度，m。然而，上述模型是假设煤岩体内的裂隙为平滑裂隙，没有考虑裂隙的粗糙或迁曲度特征。针对上述平行板裂隙网络模型存在的不足，LUO 等 12 改进了上述模型，在上述模型中增加了裂隙的粗糙特征，即：Kf=l+Dr(2-D)(1-cos2f sin0.)fLI-DT/l+DTr123+DTf-Di1-df0max(4)式中：Drr为裂隙迁曲度分形维数。式（4)考虑了裂隙的粗糙特征，因而也更加符合实际情况。因此，在后续章节将主要结合该渗透率模型和粗糙离散裂隙网络建模方法分析气体在裂隙内的流动特性。3煤层气流动数值模拟由于实际煤层气的流动是一个极其复杂的过程，涉及的影响因素众多，研究仅分析裂隙结构特征参数的变化对煤层渗透性的影响。因此，研究重点在于关注不同裂隙结构特征对渗流的影响，故不考虑其它因素的影响，即假设气体为单相流，遵循达西定律，忽略应力、温度等因素的变化对渗流的影响。3.1边界条件和网格划分及气体流动控制方程为采用像元概率分解算法建立的粗糙离散裂隙网络几何模型，然后将建立的裂隙网络模型导到COMSOL多物理场耦合软件中进行煤层气流动数值计算，煤层气流动数值模型如图2。模型尺寸为10 0 mm100mm，在模型的左侧设置压力人口，人口压力为10 0 0 Pa，模型右侧为出口压力，压力为0，将模型下边界固定，上边界为不渗流边界。然后对模型进行网格划分，这里采用自由三角形网格对其进行划分，对特别关注的裂隙区域进行网格加密处理。三角形单元数量为57 6 17 个，边单元数量为37 3个，最小单元质量0.3931mm，平均单元质量0.8 2 46 mm。压力压力入口出口（a）粗糙裂隙网络几何模型（b）边界条件及网格划分图2煤层气流动数值模型Fig.2Numerical modeling of coalbed methane flow数值计算过程中，煤层气的流动控制方程采用COMSOL内置的Navier-Stokes方程 13，即：uV2V=Vp(5)VV=0-=Kf Vp1V=(6)式中：V为气体流速，m/s；为哈密顿算子。3.2模拟结果设定边界条件及相关参数后，通过达西稳态流求解控制方程（式（6)），可得到气体压力分布云图以及气体在裂隙内的流速分布云图，气体压力分布云图如图3，气体在裂隙网络中的流速分布云图如图4。压力/kPa1.00.80.60.40.20图3气体压力分布云图Fig.3Cloud diagram of gas pressure distribution32Safetyin Coal MinesOct.20232023年10 月No.10Vol.54煤砺发全第54卷第10 期流速/10-(m：s-)1.41.21.00.80.60.40.200.20.40.60.8出口流速/10-3（ms-1)图4气体在裂隙网络中的流速分布云图Fig.4Cloud image of gas velocity distribution in fracture network由图3可知，由于煤体内裂隙的存在，使得压力分布呈现非均匀分布特征。由图4可知，气体的流动主要发生在煤体内相互连通的裂隙网络内，裂隙的分布特征影响着气体的流动特性。4参数分析及讨论由式(4)可知，裂隙网络渗透率主要与裂隙率、分形维数、迁曲度分形维数、裂隙倾角、最大裂隙长度以及比例系数等参数有关。因此，为了进一步分析裂隙特征参数对煤层渗透性的影响规律，采用像元概率分解算法建立具有不同结构特征的粗糙离散裂隙网络模型，结合数模模拟方法分析不同裂隙特征参数对煤层渗透性的影响，裂隙特征参数见表1。表1裂隙特征参数Table1Parameter values of fracture characteristics特征长迁曲度裂隙比例最大裂隙裂隙率分形维数度/mm分形维数倾角/()系数长度/mm1000.100.251.551.851.001.0210850.0060.08011264.1裂隙率和分形维数裂隙率和分形维数是表征裂隙复杂程度的重要参数，其中裂隙率为裂隙的总面积与测量煤岩样品的比值，其值越大表明裂隙数量越多，越容易形成更多的流动通道。同样，分形维数也是表征裂隙复杂程度常用的定量表征参数，是裂隙占有度量空间的多少，裂隙分形维数越大，裂隙分布越密集，越容易形成更多的渗流通道。为具体分析裂隙率和分形维数的变化对煤层渗透性的影响，固定其它裂隙参数不变(其中每个案例中裂隙的倾角为45、迁曲度分形维数1.0 15、比例系数为0.0 14以及最大裂隙长度为17 mm），采用像元概率分解算法生成具有不同裂隙结构参数的粗糙离散裂隙网络几何模型，裂隙率和分形维数的变化范围分别为0.1 0.2 5、1.55 1.8 5。裂隙率与裂隙网络模型出口流速的关系如图5，分形维数与裂隙网络模型出口流速的关系如图6。1.010-39.010-41S8.010-47.010-46.010-4y=0.002 87x:0.7675.010-4R2=0.9714.010-40.080.10 0.120.140.16 0.180.200.220.24裂隙率图5裂隙率与裂隙网络模型出口流速的关系Fig.5Relationship between fracturerate andoutlet velocity of fracture network model从图5可以看出：随着裂隙率的增大，模型出口流速也随之增大，通过对模拟数据拟合发现，二者基本呈现较好的幂律增长关系。由图6 可知：裂隙网络分形维数的增加使得模型出口速度也随之增大。对其数据进行拟合发现，分形维数与出口流速之间同样展现了较好的33SafetyinCoalMinesOct.20232023年10 月No.10Vol.54煤砺发全第54卷第10 期1.0510-39.0010-4IS7.5010-4S.U)/率6.0010-4y=3.453.x5.581R2-0.9924.5010-43.0010-41.551.601.651.701.751.801.85分形维数图6分形维数与裂隙网络模型出口流速的关系Fig.6Relationship between fractal dimension andoutlet velocity of fracture network model幂律增长关系4.2迁曲度分形维数和最大裂隙长度分析裂隙迁曲度分形维数的变化对煤层渗透性的影响所需参数见表1。固定其它裂隙参数不变，迁曲度分形维数和最大裂隙长度的变化范围分别为1.0 1.0 2 和11 2 6 mm。迁曲度分形维数与裂隙网络模型出口流速的关系如图7，最大裂隙长度与裂隙网络模型出口流速的关系如图8。1.6510-41.5010-4y=1.69410-4x-27.835R2=0.9321.3510-4S.U)/率41.2010-41.0510-49.001051.0001.0051.0101.0151.020迁曲度分形维数图7迁曲度分形维数与裂隙网络模型出口流速的关系Fig.7Relationship between tortuosity fractal dimensionand outlet velocity of fracture network model由图7 可知：迁曲度分形维数与流速之间呈现了负相关关系，即随着迁曲度分形维数的增大，裂隙网络模型出口流速逐渐减小，这是因为迁曲度分形维数表征的是裂隙的弯曲程度，迁曲度分形维数值越大，说明裂隙越弯曲（或粗糙），则气3.010-32.5x10-3(f-s.u2.010-3y=17.066.x4.1471.510-3R2=0.9751.010-35.010-400.0100.0140.0180.0220.026最大裂隙长度/m图8最大裂隙长度与裂隙网络模型出口流速的关系Fig.8Relationship between maximum fracture length andoutlet velocity of fracture network model体在裂隙内的流动阻力也就越大，因此，出口流速逐渐减小。由图8 可知：随着裂隙长度的逐渐增大，模型出口流速也随之增大。这是因为裂隙最大长度的增加意味着裂隙之间更易形成网络，使得裂隙的连通性整体性增强，进而影响着煤层渗透性。4.3比例系数和裂隙倾角同样，为分析比例系数和裂隙倾角对煤层渗透性的影响，固定裂隙长度、迁曲度分形维数等参数不变，比例系数变化范围为0.0 0 6 0.0 8，倾角变化范围为10 8 5。比例系数与裂隙网络模型出口流速的关系如图9，裂隙倾角与裂隙网络模型出口流速的关系如图10。0.300.25y=424.972x2.9620.20R2=0.988s.U)/率0.150.100.0500.020.040.060.080.10比例系数图9比例系数与裂隙网络模型出口流速的关系Fig.9Relationship between proportional coefficient andoutlet velocity of fracture networkmodel由图9可知：随着比例系数的增加，裂隙网络模型出口流速逐渐增大。其原因在于，当裂隙长度不变时，比例系数与裂隙的开度成正相关关34SafetyinCoal Mines2023年10 月Oct.2023No.10Vol.54煤防全第54卷第10 期0.0120.0100.008s.)/率y=0.195x-1.2250.006R2=0.9760.0040.00201020304050 60708090裂倾角/（）图10裂隙倾角与裂隙网络模型出口流速的关系Fig.10Relationship between fracture inclination angleand outlet velocity of fracture network model系，比例系数的增加意味着裂隙开度的增大，因此气体在裂隙内的流动更为容易，导致裂隙网络模型出口流速增大。此外，由图中数据的变化可知，比例系数的增大对裂隙网络流速的影响最为显著。由图10 可知：随着倾角的逐渐增大，裂隙网络模型出口流速逐渐减小，这是由于较小裂隙倾角将更有利于气体流动，而倾角越大，气体流动阻力越大，可见煤层中裂隙倾角分布方向决定裂隙网络的渗流方向。综上，为了更好地观察裂隙特征参数与流速之间的关系，对模拟结果进行非线性拟合，拟合结果表明，裂隙率、分形维数、最大裂隙长度以及比例系数与裂隙网络出口流速之间呈现较好的幂律增长关系，而迁曲度分形维数、裂隙倾角与裂隙网络出口流速之间呈现幂律降低的关系5结语1)详细介绍了粗糙裂隙离散裂隙网络的建模理论及其实现过程，并基于该方法建立了更符合工程实际的粗糙离散裂隙网络模型。2)由于煤体内裂隙网络的存在，气体压力分布呈现非均匀分布特征，气体的流动主要发生在相互连通的裂隙网络中，煤层渗透率主要为裂隙网络所控制。3)裂隙率、分形维数、最大裂隙长度以及比例系数与裂隙网络出口流速之间呈现较好的幂律增长关系，而迁曲度分形维数、裂隙倾角与裂隙网络出口流速之间呈现幂律降低的关系；在这些参数中，比例系数、裂隙倾角和最大裂隙长度对煤岩渗透率的影响最为显著。参考文献(References):1闫霞，徐凤银，张雷，等.微构造对煤层气的控藏机理与控产模式 J.煤炭学报，2 0 2 2,47(2)：8 93-90 5.YAN Xia,XU Fengyin,ZHANG Lei,et al.Reservoir-controlling mechanism and production-controlling pat-terns of microstructure to coalbed methaneJ.Journalof China Coal Society,2022,47(2):893-905.2叶建平，侯淞译，张守仁.“十三五”期间我国煤层气勘探开发进展及下一步勘探方向 J.煤田地质与勘探,2 0 2 2,50(3):15-2 2.YE Jianping,HOU Songyi,ZHANG Shouren.Progressof coalbed methane exploration and development inChina during the 13th Five-Year Plan period and thenext exploration directionJ.Coal Geology&Explora-tion,2022,50(3):15-22.3杨琪，于岩斌，程卫民，等.基于Micro-CT的受载煤岩裂隙孔隙时空演化规律 J.岩石力学与工程报，2 0 2 2,41(S1):2626-2638.YANG Qi,YU Yanbin,CHENG Weimin,et al.Micro-CT-based temporal and spatial evolution of fracturesand pores in loaded coalJ.Chinese Journal of RockMechanics and Engineering,2022,41(S1):2626-2638.4王磊，刘怀谦，谢广祥，等.含瓦斯煤孔裂隙结构精细表征及强度劣化机制 J.岩土力学，2 0 2 1，42(12)：3203-3216.WANG Lei,LIU Huaiqian,XIE Guangxiang,et al.Finecharacterization of the pore and fracture structure andstrength degradation mechanism of gas bearing coalJ.Rock and Soil Mechanics,2021,42(12):3203-3216.5王刚，王世斌，李怀兴，等.考虑煤体孔-裂隙介尺度特征的水相渗流演化模型研究 J.岩石力学与工程学报,2 0 2 1,40(8):1547-1558.WANG Gang,WANG Shibin,LI Huaixing,et al.Studyon water phase seepage evolution model consideringmesoscale characteristics of pore and fissure in coal J.Chinese Journal of Rock Mechanics and Engineering,2021,40(8):1547-1558.6徐轩，杨正明，祖立凯，等.多重介质储层渗流的等效连续介质模型及数值模拟 J.断块油气田，2 0 10,17(6):733-737.XU Xuan,YANG Zhengming,ZU Likai,et al.Equival-ent continuous medium model and numerical simula-tion of multimedia reservoirJ.Fault-block Oil&GasField,2010,17(6):733-737.7陈必光，宋二祥，程晓辉.二维裂隙岩体渗流传热的离35SafetyinCoalMinesOct.20232023年10 月煤砺发全No.10Vol.54第10 期第54卷散裂隙网络模型数值计算方法 J.岩石力学与工程学报,2 0 14,33(1):43-51.CHEN Biguang,SONG Erxiang,CHENG Xiaohui.Anumerical method for discrete fracture network modelfor flow and heat transfer in two-dimensional fracturedrocksJ.Chinese Journal of Rock Mechanics and En-gineering,2014,33(1):43-51.8WU Mingyang,WANG Wensong,ZHANG Dongming,et al.The pixel crack reconstruction method:From frac-ture image to fracture geological model for fractureevolution simulation J.Construction and Building Ma-terials,2021,273(1):121733.9XIA Binwei,LUO Yafei,PAN Chao,et al.Coalbedmethane flow characteristics based on fractal geometryand stochastic rough fracture networkJ.EnergySources,Part A:Recovery,Utilization,and Environ-mental Effects,2021(4):1-19.10 熊祥斌，张楚汉，王恩志.岩石单裂隙稳态渗流研究进展 J.岩石力学与工程学报，2 0 0 9,2 8（9)：18 39-18 47.XIONG Xiangbin,ZHANG Chuhan,WANG Enzhi.Areview of steady state seepage in a single fracture ofrockJ.Chinese Journal of Rock Mechanics and En-gineering,2009,28(9):1839-1847.11 MIAO Tongjun,YU Boming,YE Duan,et al.A fractalanalysis of permeability for fractured rocks J.Interna-tional Journal of Heat and Mass Transfer,2015,81:75-80.12 LUO Yafei,XIA Binwei,LI Honglian,et al.Fractal per-meability model for dual-porosity media embeddedwithnatural tortuous fractures J.Fuel,2021,295:120610.13张弛，吴剑锋，陈干，等.裂隙网络生成的随机模拟研究 J.水文地质工程地质，2 0 15,42(4)：12-17.ZHANG Chi,WU Jianfeng,CHEN Gan,et al.A com-parative study of stochastic modeling for generatingfracture networksJ.Hydrogeology&EngineeringGeology,2015,42(4):12-17.