煤层瓦斯含量的直接拟合测定方法
【专利摘要】本发明公开了一种煤层瓦斯含量的拟合测定方法，首先建立拟合计算模型，根据计算推导得出煤样瓦斯解吸规律的一般表达式；然后从井下煤芯取样，依次测定井下测定地点的大气压力p1和温度T1、井下自然解吸瓦斯量；进行井下自然解吸瓦斯量数据校正，经温度、压力校正到标准状态后即可求出单位质量煤芯的瓦斯损失量。本发明利用直接拟合测定方法得到的瓦斯含量与间接法测定的瓦斯含量相对误差δ小于10％，符合现场实际要求，该方法具有可行性；该方法可进行大量测点的测量，可更好地用于煤层突出危险性工作面及区域预测、预抽瓦斯效果评价以及矿井煤层瓦斯涌出量预测等。
【专利说明】
煤层瓦斯含量的直接拟合测定方法
技术领域
[0001] 本发明涉及煤矿瓦斯含量的测定技术领域，具体涉及一种煤层瓦斯含量的直接拟 合测定方法。
【背景技术】
[0002] 煤层瓦斯含量决定矿井瓦斯涌出量的大小和矿井瓦斯等级，是煤与瓦斯突出危险 性预测和评价的重要指标。目前，国内外采用的煤层瓦斯含量测定方法包括间接法和直接 法：间接法以煤吸附瓦斯理论为基础，需测定煤层原始瓦斯压力、孔隙率、瓦斯吸附常数等 参数，耗时较长且受地质条件的影响较大;作为在现场实际测定煤层瓦斯含量时的首选方 法，直接法具有简单、省时、直观、适用性强、受地质条件影响小等优点，其缺点是需根据煤 样瓦斯解吸规律对采集煤样过程中逸散的瓦斯含量进行推算，易产生较大误差，影响煤层 瓦斯含量测定结果的准确性。
【发明内容】
[0003] 针对上述现有技术中存在的技术问题，本发明提出了一种煤层瓦斯含量的直接拟 合测定方法，该方法首先通过建立拟合计算模型;然后从井下煤芯取样，依次测量井下测定 地点的大气压力和温度、井下自然解吸瓦斯量;然后将井下自然解吸瓦斯量数据经温度、压 力校正到标准状态下的数据，计算出单位质量煤芯的自然解吸瓦斯量;最后依据拟合计算 模型，利用IstOpt软件对校正后的井下自然解吸瓦斯数据进行非线性拟合计算，得到煤层 瓦斯含量。
[0004] 本发明通过下述技术方案来实现：
[0005] -种煤层瓦斯含量的直接拟合测定方法，依次包括以下步骤：
[0006] 步骤1:拟合计算模型的建立
[0007] 作出如下假设：
[0008] ①煤肩为球形颗粒，瓦斯在煤肩中的扩散过程看作是瓦斯在煤肩球形模型的扩 散;②煤肩为均质、各向同性体;③煤肩球形模型的表面浓度为常数;④扩散系数和坐标无 关，忽略瓦斯浓度和时间对扩散系数的影响;⑤瓦斯在煤肩中的扩散过程遵从质量守恒定 律和连续性原理；
[0009] 由于瓦斯在煤肩中是沿径向方向扩散，故瓦斯浓度C只是关于半径r和时间t的函 数，则可得到球坐标系下的Fi ck第二扩散定律，令u = cr，即：
[0011] 式中:r为极坐标半径;D为扩散系数；
[0012] 则煤肩瓦斯扩散方程变为一维线性流动方程，设表面浓度为常数C1，这时初始条 件和边界条件为：
[0014] 式中：c〇为球体内部初始时刻的浓度；
[0015] 根据初始条件和边界条件，采用分离变量法求解，则：
[0019]用Qt表示扩散时间为t时累计的瓦斯扩散量，如表示扩散时间时极限瓦斯扩 散量(即煤肩吸附瓦斯含量），根据式(2)，可得：
[0024] 对式(3)经过拉普拉斯变换，可得：
[0026]忽略上式中的误差函数方程可变为：
[0028] 式(4)即煤样瓦斯解吸规律的一般表达式，是煤样瓦斯解吸量Qt关于时间4的一 元二次方程，即：这=抒而，式中〇〇〇、D、r a均为常数；
[0029] 步骤2:现场测定，依次包括以下子步骤：
[0030] 步骤2.1、井下煤芯取样；
[0031]步骤2.2、测定井下测定地点的大气压力?1和温度T1;
[0032]步骤2.3、测定井下自然解吸瓦斯量；
[0033]步骤3:校正井下自然解吸瓦斯量数据:经温度、压力校正到标准状态后即可计算 出单位质量煤芯的自然解吸瓦斯量；
[0036]式中，V1Q为换算为标准状态下瓦斯体积，ml;Pl为井下测试地点大气压力，KPa;p〇 为标准大气压力，KPa;!^*井下测定温度/CAiS单位质量煤芯的自然解吸瓦斯量，ml/g;G 为采集分析煤芯质量，g;
[0037] 步骤4:拟合计算煤层瓦斯含量，以式(4)为基本拟合方程，利用IstOpt软件对校正 后的井下自然解吸瓦斯数据进行非线性拟合计算，得到煤层瓦斯含量。
[0038] 作为本发明的一个优选方案，步骤2.1中，在未经过瓦斯抽采的石门、岩石巷道或 新暴露的采掘工作面向煤层打钻，记录起钻时间、钻进结束时间；快速取出取芯管，将所取 煤样进行适当分选，装进煤样筒，记录煤样筒密封时间；
[0039] 作为本发明的另一个优选方案，步骤2.2中，用空盒气压计和温度计测试并记录井 下测试地点大气压力P1和温度!\。
[0040] 优选的步骤2.3中，连接煤样筒与解吸仪进行井下自然解吸瓦斯量测定，用排水集 气法将瓦斯收集在量管内，当实测解吸瓦斯体积达到单根测量管最大量程85%时，打开转 换开关用第二根测量管测量;在打开煤样筒的气阀的同时记录读数和时间，读完第一个数 后立刻启动秒表，每隔30s记录解吸量管内的瓦斯量，解吸30min后终止测试。
[0041] 相比现有技术，本发明带来了以下有益技术效果：
[0042]煤样瓦斯解吸量Qt与时间t具有较好的拟合关系(相关系数R2均在0.9以上），式(4) 能够较合理的反应出煤样的瓦斯解吸过程，即：式(4)作为煤样瓦斯解吸规律的一般表达式 是合理的；测定工作量小，测定成本低，可在60min内快速测定煤层内的瓦斯含量；可在 60min内快速测定煤层内的瓦斯含量，工作量小，利用直接拟合测定方法计算煤层瓦斯含量 时无需再按照《AQ1066-2008煤层瓦斯含量井下直接测定方法》中所述测定粉碎后解吸瓦斯 量、不可解吸瓦斯量等参数或者间接法推算瓦斯含量时测定多个参数数值。
[0043] 本发明方法具有方便快捷、节约成本等特点；利用直接拟合测定方法得到的瓦斯 含量与间接法测定的瓦斯含量相对误差S(即：直接拟合测定方法得到的瓦斯含量与间接法 测定的瓦斯含量的差值占间接法测定的瓦斯含量的比例，％ )小于10%，符合现场实际要 求，该方法具有可行性;该方法可进行大量测点的测量，可更好地用于煤层突出危险性工作 面及区域预测、预抽瓦斯效果评价以及矿井煤层瓦斯涌出量预测等。
【附图说明】
[0044] 图1为本发明拟合测定方法的流程图。
【具体实施方式】
[0045] 本发明指的"煤层瓦斯含量"是指单位质量的煤在20°C和一个大气压下条件下所 直接测定和计算出的煤层瓦斯解吸量，不包括瓦斯残存量，单位为m 3/t，其表达基准为原煤 基。采用保持或代表煤层实际情况的煤样，不考虑原始煤层中的水分和非煤物质的存在。
[0046] 本发明指的"井下自然解吸瓦斯量"是指煤单位质量的煤芯从装入煤样筒开始到 解析量测定完成之前所解吸出的瓦斯量，单位为m 3/t，其表达基准为原煤基。采用保持或代 表煤层实际情况的煤样，不考虑原始煤层中的水分和非煤物质的存在。
[0047] 本发明，如图1所示，一种煤层瓦斯含量的直接拟合测定方法，依次包括以下步骤： [0048]步骤1:拟合计算模型的建立
[0049]各种采掘工艺条件下采落煤的瓦斯涌出、突出发展过程中已破碎煤的瓦斯涌出、 测定煤层瓦斯含量时的瓦斯解吸等问题，皆可归结为煤肩瓦斯的扩散问题。结合目前研究 成果，作出如下假设:①煤肩为球形颗粒，瓦斯在煤肩中的扩散过程看作是瓦斯在煤肩球形 模型的扩散;②煤肩为均质、各向同性体;③煤肩球形模型的表面浓度为常数;④扩散系数 和坐标无关，忽略瓦斯浓度和时间对扩散系数的影响;⑤瓦斯在煤肩中的扩散过程遵从质 量守恒定律和连续性原理；
[0050]由于瓦斯在煤肩中是沿径向方向扩散，故瓦斯浓度c只是关于半径r和时间t的函 数，则可得到球坐标系下的Fi ck第二扩散定律，令u = cr，即：
[0052]式中:r为极坐标半径;D为扩散系数；
[0053]则煤肩瓦斯扩散方程变为一维线性流动方程，设表面浓度为常数C1，这时初始条 件和边界条件为：
[0055] 式中：co为球体内部初始时刻的浓度；
[0056] 根据初始条件和边界条件，采用分离变量法求解，则：
[0060]用Qt表示扩散时间为t时累计的瓦斯扩散量，如表示扩散时间时极限瓦斯扩 散量(即煤肩吸附瓦斯含量），根据式(2)，可得：
[0065] 对式(3)经过拉普拉斯变换，可得：
[0067]忽略上式中的误差函数方程可变为：
[0069]式(4)即煤样瓦斯解吸规律的一般表达式，是煤样瓦斯解吸量Qt关于时间4的一 元二次方程，即：式中Q〇〇、D、ra均为常数；
[0070] 根据煤样瓦斯解吸过程的实际意义，可知式（4)的定义域
(t)取得的最大值应该是如;根据一元二次方程的特点
，Qt = F (t)取得最大 值(Qt)max = 3ftV3i;此处出现了误差，但因（Qt)max与Q〇〇有明确的关系（（Q t)max = 3Q〇〇/3i)，故仍 可用式(4)来表述煤样瓦斯解吸规律，并作为拟合计算模型；
[0071] 步骤2:现场测定
[0072] 步骤2.1-一井下煤芯取样:在未经过瓦斯抽采的石门、岩石巷道或新暴露的采掘 工作面向煤层打钻，记录起钻时间、钻进结束时间；快速取出取芯管，将所取煤样进行适当 分选，装进煤样筒，记录煤样筒密封时间；
[0073] 步骤2.2-一测定井下测定地点的大气压力和温度：用空盒气压计和温度计测试 并记录井下测试地点大气压力？1和温度?\;
[0074] 步骤2.3-一测定井下自然解吸瓦斯量:连接煤样筒与解吸仪进行井下自然解吸 瓦斯量测定，用排水集气法将瓦斯收集在量管内，当实测解吸瓦斯体积达到单根测量管最 大量程85%时，打开转换开关用第二根测量管测量;在打开煤样筒的气阀的同时记录读数 和时间，读完第一个数后立刻启动秒表，每隔30s记录解吸量管内的瓦斯量，解吸30min后终 止测试。
[0075]步骤3:校正井下自然解吸瓦斯量数据
[0076]经温度、压力校正到标准状态后即可计算出单位质量煤芯的自然解吸瓦斯量；
[0079] 式中，V1Q为换算为标准状态下瓦斯体积，ml;Pl为井下测试地点大气压力，KPa;p〇 为标准大气压力，KPa;!^*井下测定温度/CAiS单位质量煤芯的自然解吸瓦斯量，ml/g;G 为采集分析煤芯质量，g;
[0080] 步骤4:拟合计算煤层瓦斯含量
[0081]
基本拟合方程，利用IstOpt软件对校正后的井下自然解吸瓦 斯数据进行非线性拟合计算，得到煤层瓦斯含量如，算法设置:选用优化算法"1:麦夸特法： Levenberg-Marquardt(LM)"，收敛判断指标1.00E-10，最大迭代数1000，实时输出控制数 20,选用"标准(LM)+通用全局优化法"模式，重复数400,控制迭代数20,收敛判断迭代数15。 [0082] 实施例1:
[0083] 按照《GB 474-2008煤样的制备方法》、《GB/T 212-2008煤的工业分析方法》、《GB/T 19560-2008煤的高压等温吸附试验方法》等方法对陕西宋家沟煤矿（SX)、内蒙古巴彦高勒 煤矿(NMG)、新疆金塔大黄山煤矿(XJ)、肥城曹庄煤矿(CZ)、菏泽新巨龙煤矿(XJL)、菏泽郓 城煤矿(YC)、山西玉泉煤业(YQ)、淮南潘一东矿(PYD)、贵州桐梓煤矿(TZ)的煤样进行处理 及分析，得到煤样基本参数(表1)。
[0084]表1煤样基本参数
[0085]
[0086] 为验证煤样瓦斯解吸规律一般表达式的合理性和煤层瓦斯含量的直接拟合测定 方法的可行性，对9种不同变质程度煤样的前30min瓦斯解吸数据进行非线性拟合，得到的 瓦斯含量与间接法测定的瓦斯含量进行比较(表2,表中（Q t)max表示Qt = F(t)取得的最大值， 0)^表示直接拟合测定方法得到的瓦斯含量，表示间接法测定的瓦斯含量）。
[0087] 表2瓦斯含量测定结果对照表
[0088]
[0089] 从表1、2可知：
[0090] 1)煤样瓦斯解吸量Qt与时间t具有较好的拟合关系（相关系数R2均在0.9以上），式 (4)能够较合理的反应出煤样的瓦斯解吸过程，即：式(4)作为煤样瓦斯解吸规律的一般表 达式是合理的；
[0091] 2)利用直接拟合测定方法得到的瓦斯含量0κ合与间接法测定的瓦斯含量Qi?相对 误差S(即:直接拟合测定方法得到的瓦斯含量Of哈与间接法测定的瓦斯含量Qfs?的差值占间 接法测定的瓦斯含量的比例，％)小于10%，符合现场实际要求，该方法具有可行性； [0092] 3)利用直接拟合测定方法计算煤层瓦斯含量时无需再按照《AQ1066-2008煤层瓦 斯含量井下直接测定方法》中所述测定粉碎后解吸瓦斯量、不可解吸瓦斯量等参数，该方法 具有方便快捷、节约成本等特点。
[0093] 代码举例：
[0094] 附l:lst0pt软件程序代码样本
【主权项】
1. 一种煤层瓦斯含量的直接拟合测定方法，其特征在于，依次包括以下步骤： 步骤1:拟合计算模型的建立， 作出如下假设： ①煤肩为球形颗粒，瓦斯在煤肩中的扩散过程看作是瓦斯在煤肩球形模型的扩散;② 煤肩为均质、各向同性体;③煤肩球形模型的表面浓度为常数;④扩散系数和坐标无关，忽 略瓦斯浓度和时间对扩散系数的影响;⑤瓦斯在煤肩中的扩散过程遵从质量守恒定律和连 续性原理； 由于瓦斯在煤肩中是沿径向方向扩散，故瓦斯浓度C只是关于半径r和时间t的函数，则 可得到球坐标系下的Fi ck第二扩散定律，令u = cr，即：式中:r为极坐标半径;D为扩散系数； 则煤肩瓦斯扩散方程变为一维线性流动方程，设表面浓度为常数C1，这时初始条件和边 界条件为：式中：c〇为球体内部初始时刻的浓度； 根据初始条件和边界条件，采用分离变量法求解，则： ^0 ? 炉1· a 当r-Ο时，上式极限为：用Qt表示扩散时间为t时累计的瓦斯扩散量，如表示扩散时间时极限瓦斯扩散量 (即煤肩吸附瓦斯含量），根据式(2)，可得：对式(3)经过拉普拉斯变换，可得：忽略上式中的误差函数方程可变为： (4) 式(4)即煤样瓦斯解吸规律的一般表达式，是煤样瓦斯解吸量Qt关于时间i的一元二次 方程，即：式中Q〇〇、D、ra均为常数； 步骤2:现场测定，依次包括以下子步骤： 步骤2.1、井下煤芯取样； 步骤2.2、测定井下测定地点的大气压力?1和温度T1; 步骤2.3、测定井下自然解吸瓦斯量； 步骤3:校正井下自然解吸瓦斯量数据， 经温度、压力校正到标准状态后即可计算出单位质量煤芯的自然解吸瓦斯量；式中，V1Q为换算为标准状态下瓦斯体积，ml;Pl为井下测试地点大气压力，KPa;p Q为标准 大气压力，KPa; 1\为井下测定温度，°C ; 〇:为单位质量煤芯的自然解吸瓦斯量，ml/g ;G为采集 分析煤芯质量，g; 步骤4:拟合计算煤层瓦斯含量， 以式(4)为基本拟合方程，利用IstOpt软件对校正后的井下自然解吸瓦斯数据进行非 线性拟合计算，得到煤层瓦斯含量。2. 根据权利要求1所述的煤层瓦斯含量的直接拟合测定方法，其特征在于:步骤2.1中， 在未经过瓦斯抽采的石门、岩石巷道或新暴露的采掘工作面向煤层打钻，记录起钻时间、钻 进结束时间，快速取出取芯管，将所取煤样进行适当分选，装进煤样筒，记录煤样筒密封时 间。3. 根据权利要求1所述的煤层瓦斯含量的直接拟合测定方法，其特征在于:步骤2.2中， 用空盒气压计和温度计测试并记录井下测试地点大气压力 ？1和温度!\。4. 根据权利要求1所述的煤层瓦斯含量的直接拟合测定方法，其特征在于，步骤2.3井 下自然解吸瓦斯量的测定包括:连接煤样筒与解吸仪进行井下自然解吸瓦斯量测定，用排 水集气法将瓦斯收集在量管内，当实测解吸瓦斯体积达到单根测量管最大量程85%时，打 开转换开关用第二根测量管测量;在打开煤样筒的气阀的同时记录读数和时间，读完第一 个数后立刻启动秒表，每隔30s记录解吸量管内的瓦斯量，解吸30min后终止测试。
【文档编号】G01N7/14GK105865970SQ201610183993
【公开日】2016年8月17日
【申请日】2016年3月28日
【发明人】颜国强, 杜文州, 黄启铭, 张孝强, 孙路路
【申请人】山东科技大学