煤层气解吸曲线的拟合模型马飞英1，王永清1，王林2，张强3，蒋睿1（1.“油气藏地质及开发工程”国家重点实验室·西南石油大学 四川 成都 610500；2.中石油浙江油田分公司西南采气厂 四川 宜宾 645250；3.中石化西南油气分公司广西采油厂 广西 田东 531599）摘要：煤层气解吸曲线反映了降压解吸的关系针对目前煤层气解吸理论存在的问题，开展了对煤层气解吸理论修正的必要性分析结果表明：对解吸曲线进行拟合，建立一个精确的解吸模型，不但可以解决煤层气解吸理论存在的问题，而且能指导煤层气产量预测及开发于是，在现有模型的基础上，提出了一个新的煤层气解吸曲线的拟合模型，并借用马东民教授及其学生对采集的20个煤样进行吸附/解吸实验所测得的实验数据进行了验证对比分析发现，新模型更加符合解吸曲线的趋势，20个煤样的解吸数据拟合结果的相关系数的平方更高关键词：煤层气，解吸曲线，拟合模型中图分类号：TE31 文献标识码：A引 言与常规天然气藏相比，煤层气藏不同的是，大部分煤层气吸附在基质孔隙内表面上，只有少量的甲烷处于游离状态存在于裂缝或溶解于地层水中煤层气的产出，首先需要通过降压将吸附态的甲烷气解吸下来变成游离态。
解吸是煤层气产出的基础，降压解吸规律对煤层气开采非常重要[1]起初，人们都认为解吸是吸附的逆过程，解吸等温线和吸附等温线应该是吻合的，并利用吸附等温方程来描述解吸过程2005年，张遂安等[2]学者进行了煤层气吸附/解吸可逆性实验研究，发现煤层气的解吸曲线与吸附曲线虽然很相似，但存在着差异2008年，马东民教授[3]利用自行设计的AST-1000型煤层气吸附/解吸测定仪做了大量的煤层气吸附/解吸实验，结果表明：解吸曲线要滞后于吸附曲线煤层气的解吸过程既然不同于吸附过程，那么如何对解吸曲线进行描述，目前还没有一种被大家认可的解吸曲线拟合模型于是，本文针对如何拟合解吸曲线这个问题，提出对已有模型进行修正，并借用马东民教授及其学生做的吸附/解吸实验所测得的实验数据进行了验证1煤层气解吸理论1.1存在的问题目前，国内外煤层气井产气机理主要是：解吸-扩散-渗流，可见解吸作用是煤层气井产气的关键[4]相对于吸附理论来说，煤层气解吸理论方面的研究较少[5-11]，认识不足，有待完善用吸附理论来描述解吸过程主要存在的问题有：第一，对煤层气井临界解吸压力的预测存在误差[12]第二，对煤层气井产量的计算以及采收率的预测存在影响。
第三，对最大解吸量的计算存在错误1.2修正的必要性煤层气的解吸过程既然不同于吸附过程，具有“可逆和滞后”的双重特性，那么对之前人们用吸附理论来描述解吸过程进行修正是必要的因为建立一个精确的解吸模型，不但可以解决煤层气解吸理论存在的问题，而且能指导煤层气产量预测及开发；利用解吸曲线能够较准确地判断临界解吸压力、计算开采过程中的煤层气解吸量和预测理论采收率[13]，以便指导地面煤层气井排采方面的工作3煤层气解吸曲线的拟合模型最初，人们认为解吸是吸附的完全可逆过程，并用Langmuir等温吸附方程来描述煤层气的解吸过程，即Langmuir方程： （1）式中，为储层压力，；为煤层气解吸到压力时的煤层气残余吸附量，；为最大吸附量，也就是Langmuir体积，；，为Langmuir压力，；后来，马东民等学者[14]通过实验研究得出：煤层气的吸附/解吸过程具有可逆性和解吸过程的滞后性，说明煤层气的等温吸附曲线与解吸曲线并不重合通过将解吸曲线与标准的Langmuir曲线进行对比，马东民教授认为其差值相当于枯竭压力下的残余吸附量据此，马东民教授在Langmuir等温吸附方程后加入一常数项来表征煤层气藏枯竭压力下的残余吸附量，即Langmuir解吸式： （2）然而实际的实验数据表明：吸附曲线与解吸曲线并不是马东民教授所说的那样，只相差一个常数，两曲线平行；而是随着压力的增加，吸附曲线和解吸曲线之间的差异越来越小，最终趋于零，两曲线慢慢相交成一条曲线，如图1所示。
于是，在Langmuir解吸式的基础上，对常数进行修正，用一指数函数取而代之，命名为Langmuir指数解吸式： （3）式中，为煤层气藏的残余吸附量，，当取零时，即为Langmuir等温吸附方程；为煤层气藏枯竭压力，；为曲线差异递减速率常数，即越大吸附与解吸曲线的差异随着压力增加就减小的越快[15]，4实验数据验证马东民教授的学生刘永彬从我国中高阶烟煤矿区采集了20个煤样，利用高压容量法对煤样进行了吸附/解吸实验，获得了吸附过程与解吸过程中过剩吸附量的实验数据，并按照质量守恒原理对过剩吸附量进行了校正，得到了校正吸附量（绝对吸附量的近似值）[16]吸附实验数据用Langmuir公式拟合，对于解吸实验数据，利用Langmuir公式、Langmuir解吸式和Langmuir指数解吸式对其进行拟合，拟合参数值见表1，表2从表1和表2可以看出，Langmuir公式对吸附数据的拟合，效果很好，误差极小，拟合结果的相关系数的平方基本上都在99%以上，拟合图见图1所示；而对于解吸数据的拟合，Langmuir公式效果最差，Langmuir解吸式次之，修正模型Langmuir指数解吸式效果最好（拟合结果的相关系数的平方都在99%以上），拟合图见图1所示。
通过对比分析发现，修正模型Langmuir指数解吸式能够更好地描述解吸曲线，揭示煤层气的压降解吸规律图1 不同岩芯的等温吸附/解吸实验数据拟合曲线Figure 1 The fitting curves of the different coal cores adsorption/desorption experiment data5解吸滞后原因早在1992年，A T 艾鲁尼[ ]认为甲烷在煤储层中的赋存方式除了游离态、吸附态、溶解态之外，还有可能存在固溶态2004年，Alexeev等[ ]学者通过对煤样进行核磁共振测量，发现煤中甲烷的存在形态不仅仅呈吸附态和游离态，还存在固溶态，且固溶态甲烷是以晶体形式赋存于煤基质中的由于在降压解吸条件下，固溶态是解吸不出来的，从而导致了煤层气解吸滞后6结论（1）降压条件下，降压解吸曲线更加符合，而不是或2）煤储层中的固溶态甲烷导致了解吸滞后3）修正模型Langmuir指数解吸式能够更好地描述解吸曲线，揭示煤层气的压降解吸规律，对煤层气的开采具有一定的指导意义表1 吸附/解吸实验数据拟合参数值Table 1 The fitting parameters value of adsorption/desorption experiment data煤样名称平衡水分/%吸附过程：Langmuir公式拟合参数值解吸过程：Langmuir解吸式拟合参数值伊16.526.76680.09650.98203.58810.19350.98480.9899杭15.298.02250.08480.99336.81730.08020.77190.9973乌14.8624.50870.18600.996722.87550.11804.16280.9984杭24.739.84130.07100.99066.23400.11860.80410.9993贵94.4528.07530.32070.995625.46950.14246.80490.9967贵81.0021.37460.32660.997918.35040.35332.27710.9994伊23.446.53650.15640.98964.21310.25921.23520.9957平10.7027.12820.23810.998721.87050.15616.15410.9979平20.8231.51490.24760.998928.40360.12897.14220.9987贵21.3616.91520.21490.995714.39740.20272.05500.9980和12.4632.92310.42810.996029.60230.34274.82170.9990和34.3637.06420.33060.996833.94000.25554.60960.9990和22.0630.50650.45660.997026.75270.37354.73090.9990贵64.0535.42820.32040.996231.49320.30714.24500.9990贵44.1436.66540.34430.998433.18440.30333.75560.9996贵54.4736.69630.32250.998435.41950.26762.65940.9996贵72.2326.87400.38760.994222.63340.35994.60120.9987贵11.5242.76110.27990.994937.33760.20477.02040.9981贵32.3432.90360.46930.994928.03520.32826.26520.9969贵103.1048.48810.31810.996444.45810.20637.18370.9986表2 吸附/解吸实验数据拟合参数值Table 2 The fitting parameters value of adsorption/desorption experiment data煤样名称平衡水分/%解吸过程：Langmuir公式拟合参数值解吸过程：Langmuir指数解吸式拟合参数值伊16.523.71400.62740.96364.37060.32620.90100.04300.63040.9900杭15.294.38060.38060.95186.16540.15680.82670.00590.67190.9979乌14.8617.83590.53740.947823.41920.23213.50100.31120.72050.9987杭24.735.19800.32900.97946.65420.16820.75480.11450.44810.9993贵94.4522.88050.78740.924029.66670.27285.58940.30490.62320.9967贵81.0018.68200.62580.992120.86580.39142.09030.00030.52430.9996伊23.444.53640.86380.95995.07650.50491.02950.25031.36650.9967平10.7019.83610.84040.931524.75210.34934.93620.30120.89520.9981平20.8223.93630.71430.936230.68240.29215.27740.45140.79670.9992贵21.3613.15960.52760.979316.71980.23211.94250.04440.37650.9984和12.4631.43360.63040.991434.11050.43603.31370.34870.73080.9991和34.3634.38540.47070.991439.43230.27554.27590.00020.34690.9992和22.0628.79590.70950.990431.44880.46343.70670.15490.66930.9992贵64.0533.39390.48700.994634.79960.26134.66030.14600.05130.9992贵44.1433.31150.52000.993637.36450.33403.44580.00130.46580.9997贵54.4734.76340.40160.996438.37290.28732.31660.15190.43790.9997贵72.2325.32880.66320.990526.13440.28095.13070.00280.05360.9990贵11.5236.01770.57940.971245.22850.24205.65830.45680.36560.9986贵32.3430.67810.73770.983935.42500.36124.05070.87440.39370.9977贵103.1043.16650.50930.979052.66110.24366.01680.42580.37770.9989参考文献[1]马东民,李卫波,蔺亚兵.降压解吸关系式在中高阶煤煤层气排采中的应用[J].西安科技大学学报, 2010, 30(6):697-701.[2]张遂安,叶建平,唐书恒等.煤对甲烷气体吸附-解吸机理的可逆性实验研究[J].天然气工业,2005,25(1):44-46.[3]马东民.煤层气吸附/解吸机理研究[D].西安科技大学博士学位论文.2008.[4]Gamson P, Beamish B, Johnson D. Coal microstructure and mineralization：their effect on methane recovery. In: Coalbed methane and coal geology. (eds. by Gayer R. and Harrie I.),Geology Society Specil Publication 1996；(109): 165-179.[5]宋岩,张新民,柳少波. 中国煤层气基础研究和勘探开发技术新进展[J]. 天然气工业, 2005, 25(1):1-7.[6]马东民,韦波,蔡忠勇. 煤层气解吸特征的实验研究[J]. 地质学报, 2008, 82(10):1432-1436.[7]杨新乐,张永利,李成全,等. 考虑温度影响下煤层气解吸渗流规律试验研究[J]. 岩土工程学报, 2008, 30(12):1811-1814.[8]曾社教,马东民,王鹏刚. 温度变化对煤层气解吸效果的影响[J]. 西安科技大学学报, 2009, 29(4):449-453.[9]马东民,张遂安,王鹏刚,等. 煤层气解吸的温度效应机制[J]. 煤田地质与勘探, 2010, 38(6):39-42.[10]马东民,张遂安,王鹏刚,等. 煤层气解吸的温度效应[J]. 煤田地质与勘探, 2011, 39(1):20-23.[11]马东民,蔺亚兵,张遂安. 煤层气升温解吸特征分析与应用[J]. 中国煤层气, 2011, 8(3):11-15.[12]马东民,马薇,蔺亚兵.煤层气解吸滞后特征分析[J].煤炭学报,2012,37(11):1885-1889.[13]李小彦,司胜利.我国煤储层煤层气解吸特征[J].煤田地质与勘探,2004,32(3):27-29.[14]马东民,张遂安,蔺亚兵.煤的等温吸附-解吸实验及其精确拟合[J].煤炭学报,2011,36(3):477-480.[15]张强,章双龙,杨玲智.应力敏感煤层气藏水平井压力动态分析[J].特种油气藏,2012,19(6):108-111.[16]刘永彬.中高阶烟煤对甲烷的吸附/解吸特征研究[D].西安科技大学,2010.The fitting model of coalbed methane desorption curveMa Feiying1，Wang Yongqing1，Wang Lin2，Zhang Qiang3，Jiang Rui1(1.State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Southwest Petroleum University, Chengdu, Sichuan 610500, China; 2.Southwest gas production plant, Zhejiang Oil Field Branch Company, CNPC, Yibin Sichuan 645250, China; 3.Guangxi Oil Production Plant, Southwest Oil and Gas Company, SINOPEC, guangxi tiandong 531599, China)Abstract：Coalbed methane desorption curve reflects the relationship between drop pressure and desorption. In view of the coalbed methane desorption theory problems, we developed the analysis of the necessity of theory revision of coalbed methane desorption. The results show: to fit the desorption curve, set up an accurate desorption model, which not only can solve the problems of coalbed methane desorption theory, but also guide the prediction of coalbed methane production and development. Therefore, on the basis of the existing models, we put forward a new fitting model for coalbed methane desorption curve, and it was verified by the measured experimental data of Professor Ma Dongmin and his students on the acquisition of 20 coal samples adsorption/desorption experiments. Through the comparative analysis, we found that the new model is more consistent with the desorption curve, the square of the correlation coefficient of the desorption data fitting results of 20 coal samples is higher.Key words: coalbed methane, desorption curve, fitting model。