CN108491663A - 基于质量守恒的煤矿井下空间瓦斯分布规律计算方法 - Google Patents
基于质量守恒的煤矿井下空间瓦斯分布规律计算方法 Download PDFInfo
- Publication number
- CN108491663A CN108491663A CN201810296297.9A CN201810296297A CN108491663A CN 108491663 A CN108491663 A CN 108491663A CN 201810296297 A CN201810296297 A CN 201810296297A CN 108491663 A CN108491663 A CN 108491663A
- Authority
- CN
- China
- Prior art keywords
- gas
- network
- node
- network node
- mine
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Abstract
本发明公开了一种基于质量守恒的煤矿井下空间瓦斯分布规律计算方法,该计算方法基于通风网络,充分利用煤矿现有监控系统,获得矿井部分巷道瓦斯浓度,能实时计算出矿井巷道瓦斯分布;若发生瓦斯事故,能准确地判断其产生的灾害气体可能波及的范围,为矿井迅速采取有效的应急处理措施提供决策依据,从而避免发生连锁瓦斯事故,将已发生的瓦斯灾害带来的损失控制在最低程度;为煤矿掌握井下巷道瓦斯分布提供技术支持,为矿井通风瓦斯安全管理及安全生产提供科学管理工具和先进技术手段。
Description
技术领域
本发明涉及一种基于质量守恒的煤矿井下空间瓦斯分布规律计算方法。
背景技术
随着我国煤矿开采向深部空间发展和矿井开采规模的不断扩大,煤矿瓦斯涌出量大幅增长,是制约矿井安全生产的重要瓶颈。由于日常瓦斯巡检工作量大,监控传感器布置覆盖范围小、数据挖掘分析欠缺等原因,目前依然对全矿井空间瓦斯分布情况无法全面直观把控,不能给系统布置、采掘工作和瓦斯治理提供科学的技术支撑。因此,应采取有效手段以明晰煤矿井下空间瓦斯分布规律,为矿井安全生产提供技术保障。
在正常生产时期,需要及时准确掌握井下通风系统运行状况,并对各条巷道瓦斯浓度即矿井巷道瓦斯分布准确把控,为矿井风量调控、瓦斯抽采、采煤掘进等工作提供指导;若矿井发生瓦斯事故,需要准确判断其产生的灾害气体可能波及的范围及波及区域瓦斯浓度,为矿井迅速采取有效的应急处理措施提供决策依据,从而避免发生连锁瓦斯事故,将已发生的瓦斯灾害带来的损失控制在最低程度。
发明内容
本发明的目的是提供一种基于质量守恒的煤矿井下空间瓦斯分布规律计算方法,该方法可获得矿井部分巷道瓦斯浓度,能实时计算出矿井巷道瓦斯分布。
为解决上述技术问题,本发明提供一种基于质量守恒的煤矿井下空间瓦斯分布规律计算方法,包括以下步骤:
S1:获取矿井通风网络,根据所述矿井通风网络构建几何网络,所述几何网络包括邻接矩阵A,网络节点集合V,有向边集合E;
S2:创建未赋值有向边集合E′,E′为可变集合;
S3:选取有向边集合E′中一条有向边e,获取有向边e的始节点v1的前驱网络节点集合V′,并根据矿井通风理论与瓦斯流动质量质量守恒原理得出未计算的网络节点流出的瓦斯浓度C与其前驱网络节点流入的瓦斯浓度Ci′之间的关系式,同时计算得出网络节点v单位时间段流出的瓦斯浓度C。
S4:网络节点v单位时间段流出的瓦斯浓度C赋值给有向边e;
S5:从未赋值的有向边集合E′中移除已赋值有向边e;
S6:判断集合E′中是否存在未赋值的有向边,若存在,则返回步骤S3直至所有有向边赋值完成。
进一步地,所述步骤S3具体包括如下步骤:
S31:选取有向边集合E′中一条有向边e;
S32:判断该有向边e是否为已知瓦斯浓度C′的有向边,若判断结果为是,则将有向边e瓦斯浓度赋值为C′;若判断结果为否,则搜索邻接矩阵A,得到有向边e的始节点v1的前驱网络节点集合V′;
S33:判断步骤S32中所述网络节点集合V′中是否有未计算的网络节点,若判断没有未计算的网络节点,则从向边集合E′中重新获取一条有向边e,返回上述步骤S31;若判断有未计算的网络节点,则根据矿井通风理论与瓦斯流动质量质量守恒原理得出未计算的网络节点单位时间流出的瓦斯浓度C与其前驱网络节点流入的瓦斯浓度Ci′之间的关系式,进一步计算得出网络节点单位时间段流出的瓦斯浓度C;
其中,网络节点单位时间流出的瓦斯浓度C为:
Qi为单位时间内节点前驱网络节点流入风量;Ci′为节点前驱网络节点瓦斯浓度;T为单位时间;n为该节点前驱节点数。
进一步地,所述有向边记为:e<v1,v2>e,v1为始节点,v2为末节点,且令有向边方向为矿井通风网络中巷道风向。
进一步地,所述邻接矩阵A包括一个顺序存储节点信息的节点表和一个存储节点间相互关系的关系矩阵组成,用于反映几何网络中节点间的连通关系。
进一步地,所述矿井通风理论具体为:网络节点流入风量与网路节点流出风量相等;所述瓦斯流动质量质量守恒原理具体为:同时间段内节点流入瓦斯质量与流出瓦斯质量相等。
本发明的有益效果为:该计算方法基于通风网络,充分利用煤矿现有监控系统,获得矿井部分巷道瓦斯浓度,能实时计算出矿井巷道瓦斯分布;若发生瓦斯事故,能准确地判断其产生的灾害气体可能波及的范围,为矿井迅速采取有效的应急处理措施提供决策依据,从而避免发生连锁瓦斯事故,将已发生的瓦斯灾害带来的损失控制在最低程度;为煤矿掌握井下巷道瓦斯分布提供技术支持,为矿井通风瓦斯安全管理及安全生产提供科学管理工具和先进技术手段。
附图说明
此处所说明的附图用来提供对本申请的进一步理解,构成本申请的一部分,在这些附图中使用相同的参考标号来表示相同或相似的部分,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。在附图中:
图1是本发明算法流程图;
图2是本发明算法数据框架图;
图3为是巷道瓦斯分布计算范围示意图。
具体实施方式
如图1所示,本申请提供了一种基于质量守恒的煤矿井下空间瓦斯分布规律计算方法,包括以下步骤:
S1:获取矿井通风网络,根据所述矿井通风网络构建几何网络,所述几何网络包括邻接矩阵A,网络节点集合V,有向边集合E;将节点集合V和有向边集合E记为:有向图G=(V,E);
S2:创建未赋值有向边集合E′,E′为可变集合;
S3:选取有向边集合E′中一条有向边e,其中,所述有向边记为:e<v1,v2>e,v1为始节点,v2为末节点,且令有向边方向为矿井通风网络中巷道风向;
S4:判断该有向边e是否为已知瓦斯浓度C′的有向边,若判断结果为是,则将有向边e瓦斯浓度赋值为C′;若判断结果为否,则执行步骤S5;
S5:搜索邻接矩阵A,得到有向边e的始节点v1的前驱网络节点集合V′;其中,所述邻接矩阵A包括一个顺序存储节点信息的节点表和一个存储节点间相互关系的关系矩阵组成,用于反映几何网络中节点间的连通关系。
S6:判断步骤S5中所述网络节点集合V′中是否有未计算的网络节点,若判断没有未计算的网络节点,则从向边集合E′中重新获取一条有向边e,返回上述步骤S3;若判断有未计算的网络节点,则执行步骤S7;
S7:根据矿井通风理论与瓦斯流动质量质量守恒原理得出未计算的网络节点单位时间流出的瓦斯浓度C与其前驱网络节点流入的瓦斯浓度Ci′之间的关系式,进一步计算得出网络节点单位时间段流出的瓦斯浓度C;
其中,网络节点单位时间流出的瓦斯浓度C为:
Qi为单位时间内节点前驱网络节点流入风量,Ci′为节点前驱网络节点瓦斯浓度;T为单位时间;n为该节点前驱节点数。
S8:网络节点v单位时间段流出的瓦斯浓度C赋值给有向边e;
S9:从未赋值的有向边集合E′中移除已赋值有向边e;
S10:判断集合E′中是否存在未赋值的有向边,若存在,则返回步骤S3直至所有向边赋值完成。
如图3所示,本发明中的初始数据来源于通风网络解算所得的矿井风量分布情况、矿井监控系统对关键巷道的瓦斯浓度监控数据、收集的矿井日常巡检瓦斯浓度数据;将对矿井通风网络解算得到的各个巷道的风量值与获得的瓦斯浓度值相乘得到对应巷道的瓦斯量;根据瓦斯在通风网络中流动的质量守恒原理和算法所遵循的质量守恒原理:①节点流入风量与节点流出风量相等,②同时间段内节点流入瓦斯质量与流出瓦斯质量相等,得到未知瓦斯浓度巷道的瓦斯浓度数值,计算出全矿井巷道空间瓦斯浓度分布,计算结果分布如图3所示。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (5)
1.一种基于质量守恒的煤矿井下空间瓦斯分布规律计算方法,其特征在于,包括以下步骤:
S1:获取矿井通风网络,根据所述矿井通风网络构建几何网络,所述几何网络包括邻接矩阵A,网络节点集合V,有向边集合E;
S2:创建未赋值有向边集合E′,E′为可变集合;
S3:选取有向边集合E′中一条有向边e,获取有向边e的始节点v1的前驱网络节点集合V′,并根据矿井通风理论与瓦斯流动质量质量守恒原理得出未计算的网络节点流出的瓦斯浓度C与其前驱网络节点流入的瓦斯浓度Ci′之间的关系式,同时计算得出网络节点v单位时间段流出的瓦斯浓度C。
S4:网络节点v单位时间段流出的瓦斯浓度C赋值给有向边e;
S5:从未赋值的有向边集合E′中移除已赋值有向边e;
S6:判断集合E′中是否存在未赋值的有向边,若存在,则返回步骤S3直至所有有向边赋值完成。
2.根据权利要求1所述的基于质量守恒的煤矿井下空间瓦斯分布规律计算方法,其特征在于,所述步骤S3具体包括如下步骤:
S31:选取有向边集合E′中一条有向边e;
S32:判断该有向边e是否为已知瓦斯浓度C′的有向边,若判断结果为是,则将有向边e瓦斯浓度赋值为C′;若判断结果为否,则搜索邻接矩阵A,得到有向边e的始节点v1的前驱网络节点集合V′;
S33:判断步骤S32中所述网络节点集合V′中是否有未计算的网络节点,若判断没有未计算的网络节点,则从向边集合E′中重新获取一条有向边e,返回上述步骤S31;若判断有未计算的网络节点,则根据矿井通风理论与瓦斯流动质量质量守恒原理得出未计算的网络节点单位时间流出的瓦斯浓度C与其前驱网络节点流入的瓦斯浓度Ci′之间的关系式,进一步计算得出网络节点单位时间段流出的瓦斯浓度C;
其中,网络节点单位时间流出的瓦斯浓度C为:
Qi为单位时间内节点前驱网络节点流入风量,Ci′为节点前驱网络节点瓦斯浓度;T为单位时间;n为该节点前驱节点数。
3.根据权利要求1所述的基于质量守恒的煤矿井下空间瓦斯分布规律计算方法,其特征在于,所述有向边记为:e<v1,v2>e,v1为始节点,v2为末节点,且令有向边方向为矿井通风网络中巷道风向。
4.根据权利要求1所述的基于质量守恒的煤矿井下空间瓦斯分布规律计算方法,其特征在于,所述邻接矩阵A包括一个顺序存储节点信息的节点表和一个存储节点间相互关系的关系矩阵组成,用于反映几何网络中节点间的连通关系。
5.根据权利要求1-4任一所述的基于质量守恒的煤矿井下空间瓦斯分布规律计算方法,其特征在于,所述矿井通风理论具体为:网络节点流入风量与网路节点流出风量相等;所述瓦斯流动质量质量守恒原理具体为:同时间段内节点流入瓦斯质量与流出瓦斯质量相等。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810296297.9A CN108491663B (zh) | 2018-03-30 | 2018-03-30 | 基于质量守恒的煤矿井下空间瓦斯分布规律计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810296297.9A CN108491663B (zh) | 2018-03-30 | 2018-03-30 | 基于质量守恒的煤矿井下空间瓦斯分布规律计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108491663A true CN108491663A (zh) | 2018-09-04 |
CN108491663B CN108491663B (zh) | 2022-03-22 |
Family
ID=63318105
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810296297.9A Active CN108491663B (zh) | 2018-03-30 | 2018-03-30 | 基于质量守恒的煤矿井下空间瓦斯分布规律计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108491663B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112989533A (zh) * | 2021-02-22 | 2021-06-18 | 中煤科工集团重庆研究院有限公司 | 一种基于分元质子法的矿井串联通风判识方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104775841A (zh) * | 2015-03-02 | 2015-07-15 | 陕西陕煤黄陵矿业有限公司 | 面向矿井用风点的通风网络安全分区划分方法 |
CN106503363A (zh) * | 2016-10-27 | 2017-03-15 | 中国矿业大学(北京) | 一种含高浓度瓦斯井巷瞬态风流稳定性计算方法 |
AU2017100115A4 (en) * | 2017-01-31 | 2017-03-23 | Thomas, Graham Gordon MR | Carbon Dioxide CO2 as renewable source of energy through bio-sequestration ex coal or gas fired electricity generators. |
CN107609312A (zh) * | 2017-10-17 | 2018-01-19 | 上海燃气市北销售有限公司 | 城市输气管网供气能力模拟方法及模拟系统 |
-
2018
- 2018-03-30 CN CN201810296297.9A patent/CN108491663B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104775841A (zh) * | 2015-03-02 | 2015-07-15 | 陕西陕煤黄陵矿业有限公司 | 面向矿井用风点的通风网络安全分区划分方法 |
CN106503363A (zh) * | 2016-10-27 | 2017-03-15 | 中国矿业大学(北京) | 一种含高浓度瓦斯井巷瞬态风流稳定性计算方法 |
AU2017100115A4 (en) * | 2017-01-31 | 2017-03-23 | Thomas, Graham Gordon MR | Carbon Dioxide CO2 as renewable source of energy through bio-sequestration ex coal or gas fired electricity generators. |
CN107609312A (zh) * | 2017-10-17 | 2018-01-19 | 上海燃气市北销售有限公司 | 城市输气管网供气能力模拟方法及模拟系统 |
Non-Patent Citations (1)
Title |
---|
魏引尚等: ""瓦斯在通风巷道中流动分布情况研究"", 《西安科技大学学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112989533A (zh) * | 2021-02-22 | 2021-06-18 | 中煤科工集团重庆研究院有限公司 | 一种基于分元质子法的矿井串联通风判识方法 |
CN112989533B (zh) * | 2021-02-22 | 2022-09-23 | 中煤科工集团重庆研究院有限公司 | 一种基于分元质子法的矿井串联通风判识方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108491663B (zh) | 2022-03-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Archibald et al. | An aggregate stochastic dynamic programming model of multireservoir systems | |
CN103603691B (zh) | 一种矿井通风动态解算和分析预警方法及系统 | |
CN103984994B (zh) | 一种城市轨道交通客流高峰持续时间预测方法 | |
CN108197747B (zh) | 城市轨道交通应急救援站选址方法及装置 | |
CN107633340B (zh) | 一种基于压力监测的供水管网漏失区域快速识别方法 | |
CN106777528B (zh) | 矿井需风量的全息预测方法 | |
CN107786369A (zh) | 基于irt层次分析和lstm的电力通信网络安全态势感知和预测方法 | |
CN103987054A (zh) | 一种无线网络传感器网络覆盖的分布式方法 | |
CN104809877A (zh) | 基于特征参数加权gefcm算法的高速公路地点交通状态估计方法 | |
CN109828578A (zh) | 一种基于YOLOv3的仪表巡检机器人最优路线规划方法 | |
CN108764551A (zh) | 基于系统级寿命信息的成组维修决策方法 | |
CN109724606A (zh) | 基于改进的a*算法的矿井突水避灾路径规划方法和装置 | |
CN105844384A (zh) | 一种道路安全评价方法及装置 | |
CN104407688A (zh) | 基于树回归的虚拟化云平台能耗测量方法及系统 | |
CN104636560A (zh) | 一种基于子图同构的高风险区域自动识别方法 | |
CN104775841A (zh) | 面向矿井用风点的通风网络安全分区划分方法 | |
CN108491663A (zh) | 基于质量守恒的煤矿井下空间瓦斯分布规律计算方法 | |
CN104539601A (zh) | 动态网络攻击过程可靠性分析方法及系统 | |
CN104732304A (zh) | 基于灰色人工神经网络组合模型导水裂隙带高度预测方法 | |
Osorio et al. | An efficient simulation-based optimization algorithm for large-scale transportation problems | |
Cheng et al. | A reliability allocation model and application in designing a mine ventilation system | |
CN109992579A (zh) | 一种公路基础设施多源异质数据的数据修复方法及系统 | |
CN104515591A (zh) | 一种高精度的磨煤机振动检测工艺方法 | |
CN102968883B (zh) | 自适应泥石流预警方法 | |
Long et al. | Urban road traffic congestion prediction based on knowledge graph. |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |