发明名称 煤与瓦斯突出的预测方法
摘要 本发明提供一种煤与瓦斯突出的预测方法,建立在煤矿安全监控系统基础之上;预测方法主要通过矿井安全监控系统提取设置在距离预警目标掘进面40~60m位置处瓦斯传感器传输的并以安全监控系统进行相应处理后的最近三次瓦斯浓度平均值、瓦斯浓度最大值、瓦斯浓度95%置信最大值和瓦斯浓度差值;以及通过矿井安全监控系统提取最近三次的瓦斯防突检测得到的钻屑量指标值和钻屑瓦斯解吸指标值,按照特定算法计算出当前状态下的钻屑量指标值和钻屑瓦斯解吸指标值,与相应的预测判断指标进行比对,实现对煤与瓦斯突出的连续预测。本发明集接触式预测法和非接触式连续预测法的优点于一体,提高了预测的可靠性和可信度,且具有较强的自适应能力。
申请公布号 CN102926810B 申请公布日期 2014.09.10
申请号 CN201210464368.4 申请日期 2012.11.16
申请人 天地(常州)自动化股份有限公司;中煤科工集团常州自动化研究院 发明人 屈世甲;刘水文;李继来;吕鹏飞
分类号 E21F17/18(2006.01)I 主分类号 E21F17/18(2006.01)I
代理机构 常州市江海阳光知识产权代理有限公司 32214 代理人 张兢
主权项 一种煤与瓦斯突出的预测方法,建立在煤矿安全监控系统基础之上;所述煤矿安全监控系统存储有通过实际检测得到的钻屑量指标值和钻屑瓦斯解吸指标值;且煤矿安全监控系统实时接收和存储矿井掘进面的瓦斯数据并进行相应处理,其特征在于:主要包括以下步骤:①煤矿安全监控系统提取由设置在距离预警目标掘进面40~60m位置处的瓦斯传感器传输的并经煤矿安全监控系统进行相应处理后的下列实时瓦斯数据:瓦斯浓度平均值C:一个预测区间或者作业循环内瓦斯浓度的平均值;其中,最近三次值C<sub>1</sub>、C<sub>2</sub>、C<sub>3</sub>,当前值C<sub>0</sub>;瓦斯浓度最大值M:一个预测区间或者作业循环内瓦斯浓度的最大值;其中,最近三次值M<sub>1</sub>、M<sub>2</sub>、M<sub>3</sub>,当前值M<sub>0</sub>;瓦斯浓度95%置信最大值F:一个预测区间或者作业循环内瓦斯浓度95%置信区间的最大值;其中,最近三次值F<sub>1</sub>、F<sub>2</sub>、F<sub>3</sub>,当前值F<sub>0</sub>;瓦斯浓度差值D:一定时间内或者两次巡检的瓦斯浓度差值;其中,最近三次值D<sub>1</sub>、D<sub>2</sub>、D<sub>3</sub>,当前值D<sub>0</sub>;②煤矿安全监控系统提取预警目标掘进面最近三次的瓦斯防突检测指标值:其中,第一次钻屑量指标值S<sub>max1</sub>、第二次值S<sub>max2</sub>、第三次值S<sub>max3</sub>;当前钻屑量指标值S<sub>max</sub>;第一次钻屑瓦斯解吸指标值K<sub>11</sub>、第二次值K<sub>12</sub>、第三次值K<sub>13</sub>;当前钻屑瓦斯解吸指标值K<sub>1</sub>;③将由步骤①和步骤②所得数据进行如下计算:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='{' close='}'><mtable><mtr><mtd><msub><mi>A</mi><mn>1</mn></msub><msub><mi>C</mi><mn>1</mn></msub><mo>+</mo><msub><mi>A</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>M</mi><mn>1</mn></msub><mo>+</mo><msub><mi>F</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>A</mi><mn>3</mn></msub><msub><mi>D</mi><mn>1</mn></msub><mo>=</mo><msub><mi>K</mi><mn>11</mn></msub></mtd></mtr><mtr><mtd><msub><mi>A</mi><mn>1</mn></msub><msub><mi>C</mi><mn>2</mn></msub><mo>+</mo><msub><mi>A</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>M</mi><mn>2</mn></msub><mo>+</mo><msub><mi>F</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>A</mi><mn>3</mn></msub><msub><mi>D</mi><mn>2</mn></msub><mo>=</mo><msub><mi>K</mi><mn>12</mn></msub></mtd></mtr><mtr><mtd><msub><mi>A</mi><mn>1</mn></msub><msub><mi>D</mi><mn>3</mn></msub><mo>+</mo><msub><mi>A</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>M</mi><mn>3</mn></msub><mo>+</mo><msub><mi>F</mi><mn>3</mn></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>A</mi><mn>3</mn></msub><msub><mi>D</mi><mn>3</mn></msub><mo>=</mo><msub><mi>K</mi><mn>13</mn></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000509083720000011.GIF" wi="978" he="256" /></maths>由方程组(1)解得参数A<sub>1</sub>、A<sub>2</sub>、A<sub>3</sub>值,将该三个参数值带入当前预测区间内数值进行如下计算:A<sub>1</sub>C<sub>0</sub>+A<sub>2</sub>(M<sub>0</sub>+F<sub>0</sub>)+A<sub>3</sub>D<sub>0</sub>=K<sub>1</sub> (2)由式(2)解得当前状态下K<sub>1</sub>值;④将由步骤①和步骤②所得数据进行如下计算:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open='{' close='}'><mtable><mtr><mtd><msub><mi>B</mi><mn>1</mn></msub><msub><mi>C</mi><mn>1</mn></msub><mo>+</mo><msub><mi>B</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>M</mi><mn>1</mn></msub><mo>+</mo><msub><mi>F</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>B</mi><mn>3</mn></msub><msub><mi>D</mi><mn>1</mn></msub><mo>=</mo><msub><mi>S</mi><mrow><mi>max</mi><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>B</mi><mn>1</mn></msub><msub><mi>C</mi><mn>2</mn></msub><mo>+</mo><msub><mi>B</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>M</mi><mn>2</mn></msub><mo>+</mo><msub><mi>F</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>B</mi><mn>3</mn></msub><msub><mi>D</mi><mn>2</mn></msub><mo>=</mo><msub><mi>S</mi><mrow><mi>max</mi><mn>2</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>B</mi><mn>1</mn></msub><msub><mi>C</mi><mn>3</mn></msub><mo>+</mo><msub><mi>B</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>M</mi><mn>3</mn></msub><mo>+</mo><msub><mi>F</mi><mn>3</mn></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>B</mi><mn>3</mn></msub><msub><mi>D</mi><mn>3</mn></msub><mo>=</mo><msub><mi>S</mi><mrow><mi>max</mi><mn>3</mn></mrow></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000509083720000021.GIF" wi="1022" he="257" /></maths>由方程组(3)解得参数B<sub>1</sub>、B<sub>2</sub>、B<sub>3</sub>值,将该三个参数值带入当前预测区间内数值进行如下计算:B<sub>1</sub>C<sub>0</sub>+B<sub>2</sub>(M<sub>0</sub>+F<sub>0</sub>)+B<sub>3</sub>D<sub>0</sub>=S<sub>max</sub> (4)由式(4)解得当前状态下S<sub>max</sub>值;⑤对由步骤③得到的当前状态下的K<sub>1</sub>值和由步骤④得到的当前状态下的S<sub>max</sub>进行如下分析判断:当S<sub>max</sub><4.80kg/m(4.32L/m)且K<sub>1</sub><0.40ml/(g.min<sup>1/2</sup>)时,预测结果为安全状态;当4.80kg/m(4.32L/m)≤S<sub>max</sub><6.00kg/m(5.40L/m)或者0.40ml/(g.min<sup>1/2</sup>)≤K<sub>1</sub><0.50ml/(g.min<sup>1/2</sup>)时,预测结果为威胁状态;当6.00kg/m(5.40L/m)≤S<sub>max</sub>或者0.50ml/(g.min<sup>1/2</sup>)≤K<sub>1</sub>时,预测结果为危险状态。
地址 213000 江苏省常州市新北区黄河西路219号