低煤化欠饱和煤心含气量预测数值模拟

韦 波 ，张 冀 ，李 鑫 ，胡 永 ，苏红梅

（1.新疆维吾尔自治区煤田地质局 一五六煤田地质勘探队，新疆 乌鲁木齐 830009；2.新疆维吾尔自治区地质学会，新疆 乌鲁木齐 830009；3.新疆大学 新疆中亚造山带大陆动力学与成矿预测自治区重点实验室，新疆 乌鲁木齐 830017）

煤储层含气量是煤层气资源勘探测试的关键参数[1-4]，也是表征煤储层开发潜力和确保矿井瓦斯安全的关键参数之一[4-11]；煤储层含气量测试可为煤层气资源量、储量估算和煤层气开发设计提供重要依据[4]。煤储层含气量测值的不准确性是导致煤层气资源计算出现差异的重要原因[12]。现阶段，我国煤炭资源勘探和煤层气资源勘探开发煤层含气量测试方法主要采用GB/T19559—2021《煤层气含量测定方法》提供的解吸法及其矫正计算方法[13-17]，其中损失气量计算采用最初10 个地面实测解吸气量数据，由损失气时间与解吸时间和的平方根与累计解吸气量之间线性关系倒推零时间解吸气量得出。

应用上述方法，多年来我国在多个煤层气勘探开发区块获取了大量含气量数据。然而，目前煤层含气量计算仍有以下不足：吸附气欠饱和储层现场含气量测试的可靠性缺乏数值模拟验证，特别是损失气量估算的准确程度尚未有数值模型验证；另外，上述方法对含气饱和储层（吸附气、游离气均饱和，下同）含气量测试的准确性尚不清楚，饱和储层煤心含气量测试过程模拟尚未开展。基于此，选取新疆低煤化煤（本次为长焰煤）煤心为研究对象，构建了储层煤心含气量解吸-扩散数值模型，并通过模型计算分析了吸附气欠饱和与含气饱和储层煤心解吸动态，对比分析了数值模型与现场测试损失气量、解吸气量、残余气量构成的差异性，以期为我国低煤化储层煤层气勘探提供含气性分析新思路。

1 测试方法

1.1 现场含气量相关测试

按照GB/T 19559—2021《煤层气含量测定方法》国家标准提供的方法，开展自然煤心采样，记录采样及装样时间，开展自然解吸气量连续测试和残余气测试；根据GB/T 19560—2008《煤的高压等温吸附试验方法》、GB/T 212—2008《煤的工业分析方法》等开展煤等温吸附实验和工业分析。现场煤心采样及上述相关工作开展于新疆某长焰煤储层煤层气开发先导试验区，煤心样品为长焰煤（属于低煤化度煤），且煤心对应的实测含气饱和度分别为54.77%、77.51%、99.79%。工业分析及煤岩组分分析成果见表1。

表1 煤心工业分析及煤岩组分分析成果Table 1 Coal core industrial analysis and coal rock composition analysis results

1.2 储层煤心含气量测试数值模拟流程

1.2.1 假设条件

煤是一种复杂的多孔介质，为了方便求解，通常对煤心样品做出以下假设：①煤屑由球形颗粒组成；②煤颗粒为均质、各向同性体；③CH4解吸-扩散遵从连续性原理；④扩散系数与浓度、时间和坐标无关；⑤煤屑瓦斯解吸为等温条件下的解吸过程；⑥煤心含气量测试过程中孔隙度不变。

1.2.2 控制方程

单位体积吸附气欠饱和储层煤心基质CH4质量可表述为（CH4全部吸附在煤基质表面，在储层压力条件下吸附解吸平衡，孔隙表面外不含游离CH4）：

式中：mmu为单位体积煤基质中赋存的瓦斯质量，g；φun为欠饱和储层含气饱和度，%；ϕm为基质孔隙度；Mad为空气干燥基水分含量，%；Aad为空气干燥基灰分产率，%；VL为朗缪尔体积，单分子层最大的吸附量，cm3/g；pm为基质孔隙中的CH4压力，MPa；pL为朗缪尔压力，吸附量为最大吸附量1/2 时间的吸附平衡压力，MPa；ρcoal为煤视密度，kg/m3；ρsc为标准状态下的CH4密度，g/cm3。

煤心解吸至大气环境控制方程：

式中：t为时间，s；p为煤心孔隙气体压力，MPa；Mc为甲烷分子摩尔质量，kg/mol；VM为CH4分子摩尔体积，m3/mol；R为理想气体常数，J/（mol·K）；T为煤层温度，K；D为扩散系数，10-9m2/s。

单位体积饱和储层煤心基质CH4质量可表述为（CH4在煤基质表面吸附饱和，且在储层压力条件下游离CH4充满孔隙）：

式中：mmo为初始条件下单位体积饱和储层煤心基质CH4质量，g。

饱和储层煤心含气饱和度：

1.2.3 几何模型和边界条件

几何模型为实际圆柱体煤心物理模型，几何尺寸与现场测量一致；初始煤心各点孔隙压力为根据煤层气井试井储层压力推算得出，由于解吸时煤心已与大气接触，认为边界条件煤心柱面和断面表面压力为大气压力0.1 MPa。煤心几何模型如图1。

图1 煤心几何模型 （单位：m）Fig.1 Geometric model of coal core

2 结果与讨论

2.1 欠饱和储层煤心自然解吸测试与数值模拟

采用与煤心参数近似数据开展模拟并拟合解吸体积曲线，还原现场解吸参数和煤样特性参数。储层煤心参数实测值与拟合值对比见表2，煤心xy中心截面单位体积煤心含气性变化如图2。BF-1煤心累计解吸体积与解吸时间关系如图3。

图2 煤心xy 中心截面单位体积煤心含气性变化Fig.2 Change of gas content per unit volume of coal core at each position of xy central section of coal core

图3 煤心累计解吸体积与解吸时间关系Fig.3 Relationship between cumulative desorption volume and desorption time

表2 储层煤心参数实测值与拟合值对比Table 2 Comparison between measured values and fitting values of coal core parameters of reservoir

由图2 可知：解吸开始24 h，含气量快速下降，各时刻煤心中心位置含气量最高而边缘最低。

由图3 可知：BF-1 扣除损失气时间拟合累计解吸气量（6 427.66 cm3）与现场煤心累计解吸气量（6 346.73 cm3）接近，相差2.00%；未扣损失气时间模拟累计解吸量（7 361.4 cm3）大于上述二者，显示损失气的存在；同理，BF-2 和BF-3 也有类似结果。

吸附气饱和储层煤心现场解吸与数值模拟含气量构成对比见表3，BF-1、BF-2、BF-3 煤心损失气时间（T）与解吸时间（t）和的平方根和解吸初期累计解吸气量回归关系如图4。

图4 煤心损失气时间与解吸时间和的平方根和解吸初期累计解吸气量回归关系Fig.4 Linear regression relationship between the square root of the sum of gas loss time and desorption time and cumulative amount of desorbed gas at initial stage of desorption

表3 吸附气饱和储层煤心现场解吸与数值模拟含气量构成对比Table 3 Comparison of adsorbed gas saturated reservoirs between on-site desorption and numerical simulation of cores

对于BF-1 煤心，现场测试BF-1 煤心损失气量结果为979.43 cm3，与数值模拟结果前15 min解吸气量（880.26 cm3）（图4（a））接近。现场BF-1煤心解吸气量和残余气量与数值模型计算的解吸气量和残余气量也较接近。

BF-2 和BF-3 煤心现场测试煤心损失气量结果分 别 为823.04、793.71 cm3，略 低 于 与 数 值 模拟结果前15 min 解吸气量（952.40、1 015.80 cm3）（图4（b）、图4（c）），说明吸附气饱和度较高煤心解吸初期损失气量时间平方根法计算结果可能略偏小。这可能是吸附气饱和度较高，煤心内CH4浓度高，扩散作用更显著导致。用损失气时间与解吸时间和的平方根与解吸初期累计解吸气量的线性回归关系计算BF-2 和BF-3 煤心损失气量，发现基于数值模拟计算的损失气量结果（755.37、841.77 cm3）与基于现场实测的损失气量计算结果（832.04、793.71 cm3）较为接近。说明基于数值模拟数据利用时间平方根法计算的损失气量与实测时间平方根法估算的损失气量较为接近。

现场实测BF-1、BF-2 和BF-3 煤心解吸气量和残余气量与数值模型计算出的解吸气量和残余气量也较为接近，误差不超过10.12%（表3）。说明数值模型可以近似反映吸附气欠饱和储层煤心真实解吸过程和损失气量、解吸气量、残余气量特征。

2.2 饱和储层煤心含气量测试与数值模拟

以BF-1 煤心为例，采用现场实测参数数据（表2），假设煤心处于饱和状态（即吸附气饱和、游离气在储层压力条件下充满孔隙也达到饱和状态），对现阶段现场含气量测试难度较大的含气饱和煤心解吸-扩散过程进行模拟。含气饱和煤心xy中心截面各位置单位体积煤心含气性变化如图5。

图5 含气饱和煤心xy 中心截面各位置单位体积煤心含气性变化Fig.5 Change of gas content per unit volume of coal core at each position of xy center section of gas bearing saturated coal core

由图5 可知：解吸前24 h 煤心含气量快速下降，各时刻煤心中心位置含气量最高而煤心边缘含气量最低。

损失气时间与解吸时间和的平方根和累计解吸气量的线性回归关系如图6。饱和储层煤心解吸-扩散数值模拟含气量构成预测见表4。

图6 损失气时间与解吸时间和的平方根和累计解吸气量的线性回归关系Fig.6 Linear regression relationship between the square root of the sum of gas loss time and desorption time and cumulative amount of desorbed gas

表4 饱和储层煤心解吸-扩散数值模拟含气量构成预测Table 4 Prediction of gas content composition of coal core desorption-diffusion numerical simulation in saturated core

数值模拟显示，损失气时间内（15 min）假设的饱和煤心解吸-扩散气量达到4 441.3 cm3，相同损失气时间的情况下，约为相应吸附气欠饱和储层煤心此时间段解吸气量（即本例损失气量）的4 倍。

损失气时间与解吸时间和的平方根与解吸初期累计解吸气量的线性回归关系截距为正数（R2=0.998 8），暗示采用解吸时间平方根与解吸初期累计解吸气量线性回归方法并不适用于饱和煤心损失气量的估算；采用解吸时间平方根与解吸初期累计解吸气量多项式回归方法，计算得煤心损失气量仅为433.11 cm3（R2=0.998 8），甚至小于吸附气欠饱和储层煤心损失气量估算值，说明基于时间平方根与累计解吸气量的相关性分析方法可能无法准确估算损失气量。同时测试初期时间平方根与累计解吸气量的回归关系（图6 蓝色标记）和损失气时间与累计气量的回归关系（图6 红色标记，实测测不出部分）也存在较大差异，说明饱和煤心损失气时间解吸-扩散规律较为复杂，常规回归分析法可能难以实现对损失气量的准确估算。

由表4 可知：饱和煤心损失气、解吸气、残余气占比分别为18.64%、80.90%、0.46%，其损失气占比高于吸附气欠饱和煤心损失气占比、解吸气占比低于吸附气欠饱和煤心解吸气占比。

数值模型可以近似反映吸附气欠饱和储层煤心真实解吸过程和损失气量、解吸气量、残余气量特征。目前，含气量现场测试手段很难达到对饱和煤心含气量的准确测试和估算，希望通过数值模拟为含气饱和煤心含气量分析提供思路。

3 结 语

1）构建的吸附气欠饱和煤心解吸数值模型计算的损失气量、解吸气量和残余气量与现场测试相应结果接近（误差<10.12%），可近似反映吸附气欠饱和储层煤心含气量构成。

2）对于构建的含气饱和储层煤心数值模型，损失气时间与解吸时间和的平方根与解吸初期累计解吸气量的回归分析法预测损失气量误差较大。饱和煤心损失气时间解吸-扩散规律较为复杂，常规回归分析可能难以实现对损失气的准确估算。

3）相同损失气时间条件下，本例饱和煤心损失气总含气量占比（18.64%）高于吸附气欠饱和煤心损失气占比（11.95%），饱和煤心解吸气占比（80.90%）低于吸附气欠饱和煤心解吸气占比（87.32%），饱和煤心残余气占比（0.46%）低于吸附气欠饱和煤心残余气占比（0.63%）。