密度梯度驱动的煤粒瓦斯解吸扩散模型及试验研究.pdf

第 卷第 期煤 炭 科 学 技 术 年 月 移动扫码阅读 秦跃平，徐 浩，毋 凡，等密度梯度驱动的煤粒瓦斯解吸扩散模型及试验研究煤炭科学技术，， （） ， ， ， ，，（） 密度梯度驱动的煤粒瓦斯解吸扩散模型及试验研究 秦跃平，徐 浩，毋 凡，张凤杰 （中国矿业大学（北京） 应急管理与安全工程学院，北京 ） 摘 要为解决菲克扩散模型在煤粒瓦斯解吸扩散过程中存在较大偏差的问题，更好地解释瓦斯在煤 粒中的解吸扩散机制，首先开展了封闭空间煤粒瓦斯解吸试验，得到不同初始解吸压力下的煤样瓦斯 解吸数据，并评估了 种解吸经验公式的精确性；随后回顾了浓度梯度驱动的菲克扩散模型，提出了 游离瓦斯密度梯度驱动的新扩散模型；在此基础上，建立了相关的解吸扩散数学方程，并采用有限差 分数值方法进行数值解算；最后将新模型与菲克模型的解算结果以及试验数据进行比较，讨论了新瓦 斯扩散模型的优势。 结果表明不同的解吸经验公式拟合试验数据的精确度差异较大，优选出了一种 能够在全试验时间尺度下拟合情况良好的经验公式； 种不同初始解吸压力条件下的新扩散模型的 数值预测结果在整个试验时间尺度内与试验数据基本一致；相对来说，基于菲克扩散模型数值解算的 结果准确性不高；在保证理论预测精度的前提下，新扩散模型中的微孔道扩散系数不随时间和压力而 改变；而建立在常扩散系数的假设之上的经典菲克扩散模型的解析解和数值解会产生较大误差。 新 模型可以更合理和准确地描述瓦斯解吸扩散行为。 研究结果旨在为瓦斯解吸扩散机理的认识、瓦斯 突出解吸指标及煤层瓦斯含量的标定等工作提供一些新的思路和参考。 关键词煤粒；瓦斯解吸；菲克扩散；游离瓦斯密度梯度；扩散系数 中图分类号 文献标志码 文章编号（） ， ， ， （ ， ， ，） 收稿日期；责任编辑郭 鑫 基金项目国家自然科学基金资助项目（，）；中央高校基本科研业务费专项资金资助项目（） 作者简介秦跃平（），男，山西夏县人，教授，博士。 ， ， ， ， ， ； ， ， ， ， ， ， ， 中国煤炭行业知识服务平台 w w w . c h in a c a j . n et 年第 期煤 炭 科 学 技 术第 卷 ， ； ； ； ； 引 言 随着我国煤炭资源开采深度的逐步增加，瓦斯 涌出量及强度增大，会诱发更多的瓦斯超限、瓦斯突 出等瓦斯灾害问题。 在深部煤层高强度开采阶 段，煤体的瓦斯解吸速率会更快，瓦斯释放特性会更 加复杂。 因此，深入认识并理解煤体中的瓦斯解吸 扩散规律对于精确计算预测井下瓦斯含量、瓦斯涌 出量以及瓦斯解吸指标等工作至关重要。 为了从理论上解释煤粒中的瓦斯解吸扩散机 理，菲克扩散理论首先被提出来，并且产生了广 泛而深刻的影响。 在此基础上，出现了单孔扩散模 型、双孔模型，用来表述瓦斯扩散特性及规 律。 经典菲克单孔扩散模型由于其计算简单、物理 意义明确的特点受到了人们的青睐，也是目前沿袭 最为普及的扩散模型。 国外学者 、 等根据菲克扩散模型建立了相应的数学 方程，并且提出了有关解吸扩散量的解析解公式。 后续国内学者杨其銮等、聂百胜等基于此对 经典菲克扩散模型进一步简化推导，得到了模型的 近似解析解，由此可以便利地计算常扩散系数。 这 也是目前获取煤粒中的瓦斯扩散系数的常用办法。 但是这些由经典扩散模型简化导出的近似解析解公 式与解吸扩散试验数据存在较大偏差，并不能很好 地描述全时段瓦斯解吸扩散行为。 如此一来，像 扩散系数这类工程参数的测定准确性也就无法保 证。 不少学者试图将扩散系数随时间变化的表达式 引入经典菲克扩散模型，用以保证预测结果与试验 数据相一致。 但是这种方式违背了菲克模型 最初的常扩散系数的基本假设，而且文献也指 出目前还没有足够的物理基础来证明空间内均匀的 扩散系数确实是随时间下降的。 文献主张 用达西定律或者反常扩散理论来描述煤粒瓦斯扩散 过程。 虽然也获得了不错的试验匹配结果，但是在 适用性、物理意义及推广应用等方面也存在一些争 议。 另外，根据实验室条件下的煤粒瓦斯解吸数 据表现出的扩散规律，涌现出了许多形式简单的解 吸经验公式。 各解吸经验公式由于适用条件 不同，所造成的误差也各异。 就目前来看，没有一个 经验公式能够准确表征各类条件下的全时段瓦斯解 吸过程。 上述研究已经从一定程度上表明菲克理论 模型似乎不再适合用来描述煤粒瓦斯扩散行为。 因 此，亟需建立一个更加科学的而且能在整个时间尺 度下准确描述煤粒瓦斯扩散过程的理论模型。 实质上，菲克扩散模型认为瓦斯在煤粒中的流 动过程是由浓度梯度驱动的，其所谓的浓度就是煤 粒瓦斯含量，并且没有区分开吸附态和游离态瓦斯 含量。 而在煤粒瓦斯解吸扩散试验中，主要参与流 动的实际上是游离瓦斯。 在此基础上，笔者首 先在 种初始解吸压力下开展瓦斯解吸试验，并提 出了由游离瓦斯密度梯度驱动的煤粒瓦斯扩散新模 型。 将新扩散模型的数值解算结果与试验数据对 比，并探讨了新扩散模型相对于经典菲克模型的优 势，目的是进一步验证新扩散模型的精确性、普适性 以及先进性。 解吸扩散试验 煤样选取与制备 试验采用的是云南宣威县大菁煤矿的煤样。 尽 可能地将现场采集到的新鲜大块煤样密封保存良 好，并运送到实验室以开展试验。 将煤块在 的温度下烘干后做初步干燥处理，之后冷却至室温。 将煤样用粉碎机粉碎成煤粒后，使用标准样品筛筛 选出 粒径的煤粒来进行试验。 每 组试验均统一称取 煤样，并放入真空干燥箱中 去除水分，以备试验使用。 试验过程 使用 型高温高压吸附仪，具体试验系 统及结构如图 所示。 该装置主要由样品预处理区 和煤样等温吸附试验区组成。 装置通过高精度软性 控温装置实现恒温条件下的瓦斯解吸试验。 另外该 装置的数据收集系统经由 编程，具有较高的灵 敏度，能够精确记录并计算瓦斯压力、吸附量及解吸 量等数据。 具体的试验步骤如下气密性检验对 试验装置进行气密性检测，直至气密性达到要求之 后开展试验；自由空间体积计算采用真空充氦气 的措施计算出样品罐中包含煤粒孔隙在内的自由空 间体积；瓦斯吸附过程对整个试验系统进行抽真 空处理，接着向参考罐中充入瓦斯气体，打开参考罐 和样品罐的控制阀门，向样品罐通入瓦斯。 当样品 中国煤炭行业知识服务平台 w w w . c h in a c a j . n et 秦跃平等密度梯度驱动的煤粒瓦斯解吸扩散模型及试验研究 年第 期 图 试验系统结构 罐的压力维持稳定时，煤粒中的瓦斯达到吸附解吸 平衡状态，记录此时的平衡压力值（即初始解吸压 力）。 瓦斯解吸过程瓦斯吸附达到平衡状态后， 对样品罐抽气至一个大气压，关闭样品罐的通气阀 门，使煤样在封闭空间自由解吸。 观察并且记录试 验过程中的压力变化，待到样品罐的瓦斯压力保持 平稳时，解吸试验完成。 根据相邻 个时刻的瓦斯 压力数据可以计算单位时间内单位质量煤样瓦斯解 吸量，将其相加就能得到累计瓦斯解吸量。 试验结果与分析 解吸量随时间的变化规律 煤样的累计解吸量随时间的变化趋势如图 所 示。 在试验的时间尺度范围内，煤样中的瓦斯解吸 量总体随初始吸附瓦斯压力的升高而增大。 初始吸 附瓦斯压力越大，煤粒中瓦斯解吸速率也越快。 图 不同初始压力下各煤样解吸量变化 解吸经验公式拟合效果 众多学者对煤粒中的瓦斯解吸量随时间的变化 规律进行了试验研究，并且依托各自所获得的试验 数据提出了许多解吸经验公式。 某些经验公式存在 一些局限性，难以全面并精确地描述解吸量随时间 的变化趋势。 基于试验数据， 选取了巴雷尔 式、王佑安式、孙重旭式以及秦跃平式 等部分经验公式来拟合瓦斯解吸数据，所用经验公 式见表 。 表 瓦斯解吸经验公式 公式表达式备注 巴式 为煤样暴露 内的瓦斯 解吸量， （） 王式 （ ） 为极限瓦斯解吸量， ； 为表征煤质的变化的常数 项， 孙式 和 为与煤的瓦斯含量和 结构有关的常数， 秦式 （ ） 为表征极限瓦斯解吸量的 常数， ； 为与瓦斯解 吸速率有关的常数， ； 为 的常数 为避免冗余，仅以煤样在 压力下的试 验数据为例，来展示各经验公式的拟合情况，如图 所示。 通过各自的拟合结果发现，秦跃平式的拟合 相关性系数 最接近于 ，说明其能够较准确表征 煤粒瓦斯解吸过程。 相对来说其他经验公式的拟合 效果要差一些。 图 解吸试验数据的拟合过程 为进一步直观观察秦跃平式与瓦斯解吸试验数 据的拟合效果以及精确度，将表 中的秦跃平式变 形为 （） 依照式（）对解吸试验数据进行处理，即以解 吸量的倒数为 轴，以时间函数倒数（ ）为 轴 作图，具体如图 所示。 式中 的取值为 ，就试验煤样来看，图 中 的取为 。 通过 图 更直观地看出，不同压力尺度下的秦跃平式 的拟合相关性系数 均大于 ，可以说其是能 够较为精确描述瓦斯解吸量随时间变化的经验公 式之一。 中国煤炭行业知识服务平台 w w w . c h in a c a j . n et 年第 期煤 炭 科 学 技 术第 卷 图 解吸量倒数与时间函数倒数关系曲线 煤粒瓦斯扩散模型及解算 浓度梯度驱动的菲克扩散模型 煤粒中瓦斯的解吸扩散行为一直是一个广泛而 深远的研究课题。 最初瓦斯在煤粒这种多孔介质中 的扩散被认为是符合菲克扩散理论，表述为单位时 间单位面积的分子扩散流通量与浓度梯度成 正比 （） 式中 为扩散质量流量，即单位时间单位面积上流 过的瓦斯质量， （）； 为扩散系数，为常数 值， ； 为单位体积煤体的游离和吸附瓦斯含 量， 。 许多学者对以菲克扩散理论为基础的模型进行 了深入探索，并推导出了球形颗粒经典菲克扩散模 型的理论解及近似解析解，依此能便利地求解扩散 系数， 。 但是这些解算结果只能在特定的时间 段与试验数据保持一致，在对整个时间段的瓦斯解 吸扩散特性表征方面却无能为力。 有文献指出 为得到解析解，菲克模型放宽了许多假设条件，这也 有可能导致预测结果不准确，数值解相对来说会比 解析解更合理。 笔者基于有限差分的数值方法 对菲克解吸扩散数学模型进行了解算，封闭空间内 的瓦斯解吸扩散微分方程可表示为 （） 式中 为极坐标半径，。 初始条件和边界条件为 ， ， ， ， ， ， （） 文献已经采用有限差分法并基于质量守 恒定律以及式（）和式（）建立了详细的有限差分 方程，也编写了相应的数值解算程序来计算瓦斯解 吸扩散量，这里不再赘述。 以此为参照，在本文的试 验和模拟条件下来计算菲克扩散模型，并得到煤粒 中的瓦斯解吸扩散的数值模拟数据。 游离瓦斯密度梯度驱动的新扩散模型 数学模型 菲克扩散模型主要考虑的是，煤粒瓦斯流动是 由瓦斯含量的差异引起的，并没有区分开吸附态瓦 斯和游离态瓦斯的赋存状态， 。 像一些微孔中 的吸附态瓦斯大部分情况下是处于相对静止的状 态，实际上瓦斯流动过程的主要参与者是游离瓦斯 （包含原始游离瓦斯以及表面解吸形成的游离瓦 斯）。 因此，参考菲克定律的假设，可以认为煤粒瓦 斯扩散是由游离瓦斯密度梯度驱动的。 基于传质的 思想，考虑瓦斯质量流通量与游离瓦斯密度梯度成 正比，具体数学方程为 （） 式中为煤粒中单位时间内通过单位面积的瓦斯 质量， （）；为游离瓦斯的微孔道扩散系 数， ；为游离态瓦斯密度， 。 煤粒中的瓦斯含量用 方程表示 （） 式中， 为朗格缪尔吸附常数，单位分别为 ， ；为煤的视密度， ； 为瓦斯压力，； 为孔隙率；为标准情况下的温度， ；为标 准情况下的压力， ； 为试验温度，。 将煤粒中的游离瓦斯看作是理想气体，则有 （） 式中 为瓦斯的摩尔质量， ；为通用气 体常数， （）。 假设煤粒为各向同性球形多孔介质，孔隙结构 不受气体压力的影响，即不考虑吸附引起的煤体膨 胀或解吸引起的煤体收缩。 瓦斯流动过程遵循 质量守恒定律 （ ） （） 式中为瓦斯标准密度， ； 为沿 球形煤粒半径方向的单元壳的厚度，。 联立式（）式（）可以推导出球形煤粒瓦斯 解吸扩散流动的连续性方程 中国煤炭行业知识服务平台 w w w . c h in a c a j . n et 秦跃平等密度梯度驱动的煤粒瓦斯解吸扩散模型及试验研究 年第 期 （） 瓦斯解吸流动方程的初始和边界条件为 ， ， ， ， ， ， （） 式中为解吸开始时的煤粒内部的初始瓦斯压力， ；为煤粒外表面的压力，。 由于试验是在封闭空间内进行的，瓦斯随着解 吸进程不断从煤粒扩散到外部自由空间，导致煤粒 外表面的瓦斯压力逐渐增大。 因此煤粒外表面的压 力边界处于动态变化之中，表述为 （） 式中为煤粒外表面的初始压力，试验设置的大 气压，； 为煤粒的质量，；为多样品管中自 由空间的体积，；为 时刻单位质量煤粒的累 计瓦斯解吸质量， 。 的计算公式如下 （） 有限差分数值解算 在此使用有限差分的数值方法来对上述数学方 程进行解算。 首先将球形煤粒从球心至球面划分成 个节点，各节点间的距离是逐渐减小的，符合等比 数列，如图 所示。 以 个相邻节点（用实线表示） 间的中心作同心球面（用虚线表示）。 如此相邻 个虚线球面之间便形成球壳，在中心处形成 个实 心球体。 最终得到以 点为中心的实心球和包含各 节点的 个球壳。 通过质量守恒定律可以推导出节点 到 的瓦斯流动差分方程 （ ） （ ） （） （ ，，， ； ，，，）（） 图 球形煤粒节点划分 （ ） （） 式中， ， 为划分节点编号； 为时间节点编 号；为第 个时间步长。 节点 处的差分方程为 （） （） （ ） （） 节点 的边界方程为 （）（） 其中， （）为第 时刻的单位质量煤粒累 计解吸瓦斯质量，可由 节点和 节点的压力值 的差计算 （） （） 由式子（）式（）构成了第 时刻以 个 节点瓦斯压力为未知量的完整差分方程组。 文献 ， 已经由高斯迭代的思想开发了多个非 线性方程组解算的程序代码，也介绍了相关的解算 流程。 在此基础上，对建立的差分方程组进行解算， 具体解算步骤不再赘述。 试验获得的是单位质量煤粒解吸出来的瓦斯体 积，为方便将模拟结果与试验数据进行比较，这里将 模拟得到的单位质量煤粒解吸出来的瓦斯质量统一 转化为瓦斯体积含量，则 时刻单位质量煤粒累计 解吸的瓦斯体积 为 中国煤炭行业知识服务平台 w w w . c h in a c a j . n et 年第 期煤 炭 科 学 技 术第 卷 （） 数值解算中的参数选取见表。 其中 、、、、 等均由试验手段测得，而 是在解算程序中不断 调试反算而确定的。 表 数值解算参数取值 参数 取值 单位 （） 来源实测实测实测实测实测通用通用通用通用通用 结果与讨论 试验与解算结果的对比 将 种压力尺度下的瓦斯解吸试验数据与对应 的密度新模型和菲克模型解算的模拟结果进行对 比。 以初始解吸压力为 为例进行说明，为更 加直观地观察模拟结果与试验数据的匹配程度，这 里将横坐标取为对数坐标进行展示，如图 所示。 图 种模型的模拟曲线与试验数据的匹配（ ） （ ） 由图 可以看出，在试验的整个时间范围内，基 于游离瓦斯密度梯度驱动的新扩散模型计算的解吸 量随时间变化曲线与试验数据数据基本相符。 有些 试验点与模拟曲线存在些许偏差，但是从整个时间 尺度来看，这些小误差是可以接受的。 而菲克模型 的计算结果只能在初始或者后期小时间段与试验保 持一致，在整个时间尺度下的预测误差较大。 总体 来说，游离瓦斯密度梯度驱动的新扩散模型更精确， 是可以用来描述煤粒中的解吸扩散流动过程的。 另 外从图 还能得出，菲克模型中的常扩散系数不能 使试验数据与解算结果相吻合；而单独的常微孔道 扩散系数值可使模拟结果与试验数据保持一致，并 且不受时间的影响。 图 展现了不同初始解吸压力情况下的密度模 型解算结果与试验数据的对比情况。 由此可知， 种初始解吸压力下的模拟结果均与试验数据吻合良 好，这进一步验证了密度模型的准确性。 其中的关键 参数微孔道扩散系数 始终为 ，即 不同初始解吸压力尺度下的微孔道扩散系数是相同 的。 虽然取了相同的 可能会导致有些试验点 与模拟曲线有轻微偏差，但从整体看这种误差是可 忽略的，也表明了初始解吸压力对煤样的微孔道扩 散系数影响不大。 概括来说，密度模型的煤粒的微 孔道扩散系数与时间、压力等状态参数无关。 图 密度模型解算与试验数据的对比（ ） （ ） 密度模型与菲克模型的讨论 近几十年有各种条件下的瓦斯扩散试验都被用 来检验上述基于菲克扩散模型推导出的解析解、近 似解以及数值解的准确性。 但是通过大量的试验数 据发现，经典菲克扩散模型的理论曲线与试验点之 间存在较大的偏差，菲克模型的常扩散系数不能够 精确表征整个试验时间尺度的瓦斯扩散过程。 在特 定的试验条件下，扩散系数必须随时间发生变化才 能确保菲克模型的解算结果与试验数据相吻合。 但 是扩散系数随时间变化的原理及物理意义并没有得 到很好的阐释。 也就是说，目前没有强有力的证据 能够表明煤粒的扩散系数是要随着某种特定的方程 表达式而改变的。 再一点就是，实际煤体中煤粒外 部瓦斯压力随时间千变万化，一旦外界条件变化，扩 散系数随时间变化的函数也将会发生变化，这也不 中国煤炭行业知识服务平台 w w w . c h in a c a j . n et 秦跃平等密度梯度驱动的煤粒瓦斯解吸扩散模型及试验研究 年第 期 方便实际应用。 因此，一些动态扩散系数模型虽然 能够保证与特定试验条件下的试验数据的吻合精 度，也不能简单移植到复杂外界条件下的煤粒瓦斯 解吸扩散的建模工作之中。 尽管菲克定律模型简 单，能够得到解析解，但无数的试验数据已证实了菲 克定律的确不适用于煤粒中的瓦斯流动。 新密度扩散模型经由有限差分数值方法解算的 结果，在整个试验解吸时间尺度能够很好地描述瓦 斯的解吸扩散行为。 并且其关键比例系数微孔道 扩散系数是一个常数值，与时间或者压力这种状态 参数无关，只与煤的孔隙结构有关。 除此之外，游离 瓦斯密度梯度驱动的新扩散模型也具有较明确的物 理意义，在表征瓦斯扩散特性方面比经典菲克模型 更加精确合理。 结论与展望 ）封闭空间条件下煤粒瓦斯解吸扩散量及解 吸速率随压力的增大而增大。 不同的解吸经验公式 与试验数据的匹配程度差异较大，相对来说，秦跃平 式解吸经验公式与试验数据的拟合度更高，更准确 一些。 ）基于游离瓦斯密度梯度驱动的新扩散模型 数值解算的结果与解吸试验数据相吻合，该模型可 用来表征煤粒瓦斯解吸扩散特性。 据此反算出的相 应的微孔道扩散系数不受时间和压力的影响。 ）菲克扩散模型的数值和解析解算结果均与 试验数据偏差较大，其常扩散系数只适用于描述解 吸过程的某些特定时间段。 而新扩散模型中的微孔 道扩散系数也是常数，其解算结果却能准确描述整 个时间尺度下的瓦斯解吸扩散特性，由此较好地解 决了菲克扩散模型误差较大的缺陷。 ）游离瓦斯密度梯度驱动的新扩散模型是对 瓦斯解吸扩散机理的丰富及充实。 本文试验选取的 煤样及试验条件是比较特殊的，今后尚需开展多种 条件下的大量重复试验工作来验证所提出的扩散 模型。 参考文献（） 齐庆新，潘一山，李海涛，等 煤矿深部开采煤岩动力灾害防控 理论基础与关键技术 煤炭学报，，（） ， ， ， ， ，（） 杨 军，闵铁军，刘斌慧，等 深部开采灾害及防治研究进展 科学技术与工程，，（） ， ， ， ，，（） ，，（） ， （） ， ， ， ，， ， ，， 赵 伟，王 亮，陈向军，等 受限空间煤粒瓦斯吸附扩散特性 及模型适配差异分析 煤炭科学技术，，（） ， ， ， ，，（） 李青松，李国红，王恩元，等 基于经典扩散模型不同粒径粒煤 瓦斯扩散特征实验研究 中国安全生产科学技术，， （） ， ， ， ，，（） 杨其銮，王佑安 煤屑瓦斯扩散理论及其应用 煤炭学报， ，（） ， ，，（） 聂百胜，郭勇义，吴世跃，等 煤粒瓦斯扩散的理论模型及其 解析解 中国矿业大学学报，，（） ， ， ， ，， （） 李志强，刘 勇，许彦鹏，等 煤粒多尺度孔隙中瓦斯扩散机 理及动扩散系数新模型 煤炭学报，，（） ， ， ， ，，（） ， ，，（） 张路路，魏建平，温志辉，等 基于动态扩散系数的煤粒瓦斯 扩散模型 中国矿业大学学报，，（） 中国煤炭行业知识服务平台 w w w . c h in a c a j . n et 年第 期煤 炭 科 学 技 术第 卷 ， ， ， ， ， （ ） ， ， ， ，（）， 秦跃平，王健，郑赟，等 煤粒瓦斯变压吸附数学模型及 数值解算 煤炭学报，，（） ， ， ， ，，（） ， ， ， ， ， ， ， ， ，，（） 张慧杰，张浪，汪东，等 构造煤的瓦斯放散特征及孔隙 结构微观解释 煤炭学报，，（） ， ， ， ，， （） 李修磊 不同解吸公式对混合粒度颗粒煤瓦斯解吸的适用性 研究 煤矿安全，，（） ， ， （） 秦跃平，郝永江，刘鹏，等 封闭空间内煤粒瓦斯解吸实验 与数值模拟 煤炭学报，，（） ， ， ， ，，（） ， ， ， ，， （） 王佑安，杨思敬 煤和瓦斯突出危险煤层的某些特征 煤 炭学报，，（） ， ，， （） 孙重旭 煤样解吸瓦斯泄出的研究及其突出煤层煤样解吸的 特点 重庆煤炭科学研究总院重庆分院， ， ， ， ， ，， 贾宏福，安丰华，彭信山 巴雷尔法计算煤瓦斯扩散系数的局 限性分析 中国安全科学学报，，（） ， ， ，，（） 徐 浩，秦跃平，毋 凡，等煤粒瓦斯定压吸附数学模型及数 值解算 矿业科学学报，，（） ， ， ， ，， （） 中国煤炭行业知识服务平台 w w w . c h in a c a j . n et