★ FEATURED ARTICLE 搞瓦斯治理的同事普遍有个路径依赖离瓦斯源越近抽采效果一定越好。去年做某矿采空区抽采方案比选时工区一位老资格技术员极力推荐在垮落带预埋管路理由是那里就是破碎煤体堆出来的瓦斯库插管进去等于直接插进瓦斯源。我们当时没有直接拍板而是用COMSOL把采空区裂隙带高位钻孔抽采和垮落带埋管抽采放进了同一个数值实验台连续跑了30天模拟。结论让会议室安静了很久垮落带方案前期确实猛但衰减速度比裂隙带方案快得多最终无论是累计抽采纯量还是浓度稳定性都落了下风。这篇文章就把这段模拟研究的建模思路、参数设定和对比结论完整复盘一遍给正在做采空区瓦斯治理方案比选或者想用COMSOL做同类研究的朋友一个参考。1. 采空区里的两类瓦斯通道裂隙带与垮落带的流动个性1.1 三带不是上下三层抽屉而是随采动演化的动态系统采空区上覆岩层在回采后会形成垮落带、裂隙带和弯曲下沉带这个说法大家都很熟。但我发现很多方案讨论里这三带被理解成了三张叠好的蛋糕层垮落带在下面、裂隙带在中间、弯曲下沉带在上面各干各的。真实情况远不是这样。垮落带是直接顶破碎后堆积形成的岩石块度大小不一块体之间的大孔隙非常多但它的孔隙结构很不稳定。随着工作面向前推进采空区后方岩层不断被上覆载荷压实垮落带深处会逐渐重新压密孔隙率和渗透率快速下降。裂隙带则是在老顶弯曲下沉和岩层离层过程中形成的裂隙网络分布范围比垮落带更大裂隙以纵向和层间裂隙为主连通性更好而且不容易被完全压实。在COMSOL建模时最忌讳的就是把这几带设置为静态均质区块。我在第一版模型里就吃过这个亏当时把垮落带设成一个均匀高渗体、裂隙带设成一个均匀中渗体结果模拟出来的抽采流量几乎是线性上升和现场钻孔的衰减曲线完全对不上。后来改成渗透率随采空区深度连续衰减的模型才开始像样。1.2 垮落带并非所有位置都“好抽”裂隙带也并非处处“没量”很多人在直觉上把垮落带等同于高渗透区这是对的但只对了一半。垮落带新鲜垮落段确实高渗压实段却未必比裂隙带强。裂隙带的优势在于“面状连通”它覆盖的卸压范围和瓦斯解吸范围更广抽采负压可以通过裂隙网络向远处延伸这比垮落带局部的大孔隙通道更有利于长期抽采。我把典型参数整理过一张表建模初期直接按这个量级去赋值问题不大区域孔隙率参考区间渗透率参考区间流态特征垮落带新鲜垮落区0.25—0.4010^-9—10^-10 m²强透可能偏离达西线性垮落带压实区0.10—0.1810^-11—10^-12 m²达西流为主裂隙带裂隙率0.01—0.0510^-12—10^-13 m²达西流局部裂隙优势通道弯曲下沉带接近原岩孔隙率10^-13 m²以下近似不透气这张表只能当起点不同顶板岩性差别非常大。但我想强调一点垮落带的高渗是“块状高渗”范围小、位置浅裂隙带是“网络状中渗”范围宽、连通性强。拿城市交通来类比裂隙带相当于高架路网垮落带就像地面街道。高架路网虽然每条车道不是最宽的但整体通行效率往往更高。2. 从某矿资料到COMSOL几何简化、参数分区与边界条件2.1 几何怎么取材既省算力又不丢物理我们当时拿到的某矿条件大致是走向长度约400m采高3.2m埋深350m左右顶板以中硬砂岩为主。做第一轮方案比选时我没有一上来就建全三维模型而是先用走向剖面二维模型。为什么偷懒因为我们要回答的核心问题是“裂隙带抽采和垮落带抽采的流场差异到底在哪”二维剖面已经能清楚看到负压扩散路径、钻孔位置和两类瓦斯通道的关系。模型几何取走向长度350m、竖直方向60m煤层在底部上方依次是垮落带、裂隙带和弯曲下沉带。钻孔在二维剖面里用细长矩形表示宽度按孔径换算。有人会问钻孔那么细网格是不是要卡死实际上二维模型里钻孔宽度设成0.1m左右就可以重点不是孔壁细节而是孔周围压力梯度的合理性。到了第二阶段优化钻孔参数时再上三维模型把走向钻孔、倾向钻孔和群孔叠加效果拉出来。几何简化的原则是现场条件非常复杂模型不可能面面俱到。只要把煤层、采空区、顶板三带、抽采孔这一条主链拉住其他细节都可以先省略。省略的东西靠边界条件补而不是堆几何细节。2.2 渗透率分区不能拍脑袋要体现压实演化渗透率是这类模型里最敏感的变量。COMSOL里可以直接给常数但采空区抽采模拟如果这么做结果基本没有参考价值。我把采空区沿走向分成若干段每一段渗透率随距工作面的距离变化采用指数衰减式描述形式很简单κ(d) κ_残 (κ_初始 - κ_残) × exp(-d / L0)d是到采空区起点的距离L0是压实长度常数按顶板岩性和推进速度标定。我们初始取L035m垮落带新鲜区κ_初始取10^-9 m²压实后κ_残取10^-12 m²。裂隙带的渗透率也用类似方式处理但衰减幅度比垮落带小得多这样模型里就能体现出“垮落带快速压死、裂隙带长期导通”的物理现象。瓦斯源项也不能简单放一个常数。采空区遗煤的瓦斯解吸与压力有关煤层破碎后解吸速度会随压力下降而加快但解吸总量是有限的。我建议把源项按实测钻孔衰减曲线拟合成分段函数或者用一个随时间平滑衰减的源项而不是所有网格都给定同一个产气速率。否则模拟出来的抽采曲线会过于“饱满”掩盖两种方案的差异。2.3 边界条件取舍压力边界优先于流量边界边界条件是整类模拟最容易造假的地方。很多初学者习惯把进风、回风、采空区表面全给固定流量然后发现不管怎么调模型结果都好看因为边界已经决定了答案。我们模型的设定是进风侧和回风侧给大气压边界也就是表压0允许漏风从两端补充顶底板和远离采空区的侧面给零通量边界抽采孔口给负压边界比如-20kPa表压。这里有个关键选择钻孔不能用固定流量边界。为什么因为流量是模拟的输出结果你先把流量定死了还比什么方案顶板方面弯曲下沉带虽然渗透率极低但不等于完全没气源。我把它当成一层低渗透率薄层处理而不是直接删掉。这样裂隙带上部的高位钻孔也能间接感应到弯曲下沉带里的瓦斯补给比完全封闭更接近实际。提示边界越少人为给定流量模拟结果越可信。凡是能定压力的位置尽量用压力边界。3. 两套抽采方案的COMSOL实验设计不仅比最终数值3.1 方案A裂隙带高位钻孔抽采怎么布局方案A模拟的是顶板高位钻孔这是采空区瓦斯治理里非常成熟的手段。钻孔从回风巷钻场开孔仰角穿过垮落带上部终孔落在裂隙带的O形圈裂隙发育范围。模型里简化为一个水平走向的钻孔段高度距煤层顶板18m孔径94mm孔口负压20kPa。为什么要落在18m而不是更靠近垮落带因为太低了容易受到垮落带垮落扰动钻孔可能被切断而且抽采出来的气体很多直接从垮落带漏风里来甲烷浓度起不来太高了又进入弯曲下沉带裂隙不发育负压传不出去。裂隙带中下部是抽采浓度和纯量的平衡点。3.2 方案B垮落带埋管抽采怎么布置方案B模拟的是采空区埋管抽采即把大直径抽采管路引到采空区内管口埋在垮落带里。现场通常用DN300左右的埋管沿回风巷一侧向采空区延伸。模型里抽采口位于垮落带内距顶板1.0m处孔口负压同样给20kPa。从直觉上看管口直接放在破碎煤体里离遗煤瓦斯源最近初期应该效果非常好。这确实是方案B有吸引力的地方。但模型必须把两方面因素同时考虑进去一是垮落带破碎煤体的高渗段会不断被压实二是抽采后期管口周围容易形成负压死区漏风路径可能把浓度拉低。这两点不放到模型里方案B容易被人为美化。3.3 用哪些指标评价方案才能支撑决策比选不能只看一个指标。我们同时跟踪了四项抽采纯量也就是单位时间抽出的甲烷体积单位m³/min孔口甲烷体积分数反映抽采质量30天累计抽采纯量反映长期作战能力有效影响半径即负压扰动或者甲烷浓度变化能覆盖到多远。有效影响半径我要多解释一句。有人直接把压力降大于0的定义为影响范围这样算出来往往太大没有工程意义。我建议取“压力降低值超过1kPa”的范围或者依据现场实际能测到的浓度降低范围来定义。否则一群钻孔的负压场叠加起来整个采空区都算“影响区”钻机布置方案就没法优化了。4. 30天模拟结果对照裂隙带胜在稳定垮落带胜在开局4.1 压力场形态一个是宽漏斗一个是深漏斗模拟结果出来后第一个直观差异在压力场形态。方案A的负压场从高位钻孔沿着裂隙带向外扩散形成的漏斗很宽、坡度缓。裂隙带本身连通性好负压能传到20m外的垮落带上部相当于把整个裂隙网络变成了抽采管路。方案B的负压场则明显集中在管口附近漏斗很深但半径小。管口周围破碎煤体的渗透率高负压很容易到达但出了这个区域压实区的渗透率拦住了延伸负压传不出去。用工程语言说方案A是“有一张大网在兜瓦斯”方案B是“用一根粗吸管在喝眼前的汤”。网兜得慢但面广吸管局部猛但范围有限。4.2 抽采曲线对比前期反超后期被甩开我摘一组具有代表性的模拟结果两组方案都是30天连续抽采指标方案A裂隙带高位钻孔方案B垮落带埋管孔口甲烷体积分数第1天42%65%前5天平均抽采纯量1.62 m³/min2.08 m³/min30天平均抽采纯量1.87 m³/min1.26 m³/min第30天抽采纯量1.51 m³/min0.74 m³/min30天累计抽采纯量约8.08万m³约5.44万m³有效影响半径约20m约10m方案B第一天浓度冲到65%前5天平均纯量也高于方案A完全符合它“离源近”的直觉。但第10天之后压实效应和漏风作用逐渐压过初期优势纯量一路下滑第30天只剩0.74m³/min。方案A初期不算惊艳但衰减很慢第30天仍然有1.51m³/min30天累计纯量比方案B高出约48%。模拟里有个细节很有意思方案B的孔口浓度下降速度比纯量下降更猛。垮落带压实后管口周围形成一个局部低压区压力越低漏风越容易从上下端头倒灌进来把甲烷浓度一稀释纯量自然就塌了。4.3 这两种方案的差异本质瓦斯补给源面积不同把数据拆开看两套方案的核心差距不在钻孔设备而在补给源面积。方案A通过裂隙带网络衔接的是采空区上方更大范围的卸压瓦斯和遗煤解吸瓦斯补给路径多、面积大所以抽采负压能以较低流速抽很长一段时间。方案B的补给源基本局限在管口周边的塌落煤岩体初始浓度很高但这一坨煤的瓦斯量有限抽完就抽完后续补给难以跟上。这解释了很多现场现象为什么有些采空区埋管抽采一两个月后基本成了“干抽”负压一直开流量也不小但浓度掉到百分之十几甚至更低。不是设备不行而是补给通道没有建立起来。5. 建模和调试中容易翻车的几个细节5.1 别把孔口负压设得脱离抽采泵实际能力这个坑特别常见。COMSOL里给钻孔设一个-50kPa的边界条件模型会跑得很漂亮抽采量大得惊人但现场抽采泵的绝对真空都可能达不到这个孔口负压。孔口负压不是泵的额定负压它要扣除管路阻力、钻孔局部损失和沿途漏气实际到孔口可能只有泵出口负压的三分之一到一半。我们标定时取-20kPa这是某矿钻孔孔口实测负压附近的值。如果你要去作对照模拟建议先查清楚现场孔口负压计的真实读数再决定模型条件。不要拿“泵额定负压”直接当孔口边界。5.2 渗透率跳变会导致收敛困难需要平滑过渡COMSOL遇到渗透率分区突变时经常出现两种情况要么压力场在交界处出现锯齿要么非线性求解器卡住。垮落带和裂隙带之间的渗透率可能相差两个数量级如果不做过渡交界处数值上的“墙”会非常明显。我的做法是渗透率和孔隙率都不用阶跃函数而是用解析函数加光滑过渡带。过渡带取2—5m参数从一个区域平滑过渡到另一个区域。这么做收敛快得多结果也更接近岩体内部连续损伤演化的真实状态。5.3 压力场和浓度场要分开看待收敛达西压力场通常收敛很快真正麻烦的是浓度场。甲烷浓度在多孔介质里传输时如果网格不够细或者求解器没有开流线稳定化很容易出现局部浓度超过100%甚至负浓度这种明显失真。遇到这种问题不要着急调全局网格先看钻孔附近和压实区边界的网格尺寸。我在模型里对钻孔周围做了局部加密最小网格尺寸控制在0.02m左右其他区域允许0.5—1m。压力场和浓度场的求解器相对容差也都分开设置压力场取10^-5浓度场取10^-4这样既能保证精度又不至于让整个模型慢到跑不动。5.4 至少要用一组现场实测数据做锚点这是我最想强调的一条。没有现场数据锚定的采空区模拟只能算是数学游戏。我们当时拿一组真实高位钻孔的“孔口负压—混合流量—甲烷浓度”对应数据来反演渗透率把模型系数调到能复现现场测量值之后才敢拿它去做方案比选。反演不需要很复杂。先固定一个钻孔位置改变渗透率系数对比模拟出的孔口流量和浓度找到与实测最接近的一组参数把这组参数固定下来再跑不同方案。这一步没做后面所有的“优化结论”都可能建立在错误的渗透率地基上。提示模拟值不要求跟实测百分百吻合但趋势和量级必须对得上。如果模拟出的抽采纯量比现场高两倍以上先回去查渗透率不要急着调边界。6. 先单孔再群孔是这类模拟项目最值得推荐的推进路径6.1 单孔模型先把渗透率标定干净做群孔之前的第一个模型我只放了一个钻孔。一个孔、一个采空区、一套已知现场数据重点就是标定渗透率。参数少自由度低反演结果不容易漂移。有些项目组上来就建十几个高位钻孔加埋管的完整系统网格一两百万算一次要十几个小时结果参数错了所有孔位都在错的基础上优化既浪费时间又误导工程决策。单孔模型跑顺了后面群孔模型只是增加钻孔数量而不是重新找物理规律。6.2 群孔模型要注意负压叠加和孔口流量重分配群孔模拟不能简单把多个钻孔都设成同一个恒定负压边界。钻孔并联之后每个孔口的实际压力和流量会重新分配群孔总流量不等于单孔流量乘孔数。COMSOL里可以把钻孔设成线边界并给定一个共用的抽采总负压但真正要和现场抽采系统一致最好给每个钻孔口加上流量调节条件让模型自动分配各孔的流量。当然如果你的目的只是方案比选不是精确设计抽采系统也可以保守一点所有钻孔统一给一个中等负压看流场形态差异。但要清楚这只回答“哪种方案更优”回答不了“需要打几个孔”。6.3 数值模拟回答的是趋势问题不是预言问题我在实际项目里越来越觉得COMSOL对于采空区抽采比选的价值更多是帮我们把物理过程想清楚。它告诉我们裂隙带为什么稳定、垮落带为什么衰减快、负压该往哪个方向优先延伸这些趋势判断的价值远大于某个钻孔的精确流量预测。如果你想在某个矿套用这套思路建议按这个顺序来先收集一组真实钻孔数据再建立单孔模型标定渗透率然后跑压力场再耦合浓度场最后扩大到群孔方案比选。中间每一步都尽量用现场数据做锚点模型就会越跑越可信。我自己做这个项目时最深的感触就是COMSOL确实能算但真正让模拟结果可靠的还是建模前对采空区瓦斯流动规律的那份“笨功夫”。
阅读完成 · 觉得有帮助？
