CN116663210B - 三维巷道空间关系到通风网络解算模型的转换方法和系统 - Google Patents
三维巷道空间关系到通风网络解算模型的转换方法和系统 Download PDFInfo
- Publication number
- CN116663210B CN116663210B CN202310544884.6A CN202310544884A CN116663210B CN 116663210 B CN116663210 B CN 116663210B CN 202310544884 A CN202310544884 A CN 202310544884A CN 116663210 B CN116663210 B CN 116663210B
- Authority
- CN
- China
- Prior art keywords
- roadway
- node
- dimensional
- real
- loop
- 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.)
- Active
Links
- 238000009423 ventilation Methods 0.000 title claims abstract description 62
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000004364 calculation method Methods 0.000 title claims abstract description 28
- 238000006243 chemical reaction Methods 0.000 title claims abstract description 17
- 238000010586 diagram Methods 0.000 claims abstract description 31
- 230000008569 process Effects 0.000 claims abstract description 10
- 239000011159 matrix material Substances 0.000 claims description 33
- 238000004134 energy conservation Methods 0.000 claims description 20
- 238000012800 visualization Methods 0.000 claims description 9
- 238000004040 coloring Methods 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 2
- 238000013178 mathematical model Methods 0.000 abstract description 6
- 230000008901 benefit Effects 0.000 description 5
- 238000013523 data management Methods 0.000 description 5
- 238000007726 management method Methods 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 239000003245 coal Substances 0.000 description 2
- 238000013499 data model Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Mathematical Analysis (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Software Systems (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Databases & Information Systems (AREA)
- Computer Networks & Wireless Communication (AREA)
- Remote Sensing (AREA)
- Computer Graphics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种真实三维巷道空间关系到通风网络解算模型的转换方法和应用系统。通风网络解算需将真实的三维巷道转换成符合图论中图数学特征的通风网络解算数学模型进行方程建立和求解计算。本发明利用地理信息系统(GIS)中存有的巷道空间关系和属性数据，以及拓扑关系和属性数据的强大管理功能，在真实三维巷道上增加若干属性，使得程序能够自动完成这一转换，增加的属性包括巷道位置属性、巷道方向属性。算法通过属性自动计算出进回风井口，然后添加连接进回风井口的大气分支，使得构成闭环的网络图，进而生成通风网络解算模型。该方法使得可以在GIS的三维真实巷道上直接进行通风网络解算，避免了绘制网络拓扑图，添加大气分支的繁琐过程。
Description
技术领域
本发明涉及煤矿井下通风领域，尤其涉及一种真实三维巷道空间关系到通风网络解算模型的转换方法和应用系统。
背景技术
智能通风系统建设是煤矿智能化建设的重要组成部分，智能通风系统建设要求对通风系统实现三维动态可视化分析，那么如何在真实三维巷道模型上建立通风网络解算模型是关键问题。通风网络解算数学模型，是需要在一个闭环的基于图论的网络图中，建立节点流量守恒方程组和回路能量守恒方程组，传统的做法是在表达通风系统的二维网络图的进回风井口之间通过交互式手动增加大气分支(虚拟分支)来使通风系统形成一个闭环的网络图，而在基于真实三维坐标巷道模型的三维通风仿真应用中，如果添加虚拟分支，会破坏真实巷道拓扑和度量空间关系，因此并不适用。
目前矿山的一些图形成果是基于CAD技术的，没有拓扑关系的描述功能，属性数据管理功能也较弱，在转换时十分困难。与此相比，地理信息系统(GIS)具有更加强大的空间关系和属性数据的管理功能，能够存储、管理、处理和分析各种地理空间数据，可以将地理空间数据与属性数据进行关联，支持空间数据模型和算法，数据模型更加丰富，可以更加准确地描述地理空间关系。
发明内容
鉴于上述问题，本发明提出了一种真实三维巷道空间关系到通风网络解算模型的转换方法和应用系统。
本发明实施例提供了一种真实三维巷道空间关系到通风网络解算模型的转换方法，所述转换方法包括：
在所述真实三维巷道上添加位置属性和方向属性；
根据所述真实三维巷道的真实三维坐标计算巷道的连接关系，建立巷道拓扑关系，形成开环网络拓扑图；
在开环网络拓扑图中计算图的连通性，划分成多个子网；
针对每个子网，依据所述位置属性和方向属性查找进回风井口，添加虚拟分支连接回风井井口到进风井井口，形成闭环网络图；
在所述闭环网络图上，查找最小树，圈划回路，建立节点流量守恒方程组和回路能量守恒方程组，得到通风网络解算模型。
可选地，所述位置属性包括：连通大气、内部巷道、独头巷道、废弃巷道；
所述方向属性包括：正向或者反向；
其中，所述真实三维巷道的默认方向是正向，即按照构建巷道的坐标点顺序确定的方向，实际风流方向和巷道的默认方向一致的为正向，不一致的为反向。
可选地，根据所述真实三维巷道的真实三维坐标计算巷道的连接关系，建立巷道拓扑关系，形成开环网络拓扑图，包括：
依据所述真实三维巷道的拐点坐标数据，建立直接拓扑图；
进行简化处理，合并串联分支，并在分支上记录合并的拐点坐标数据，形成所述开环网络拓扑图。
可选地，依据所述真实三维巷道的拐点坐标数据，建立直接拓扑图，包括：
记网络拓扑图的节点集合为V，边集合为E，形成的拓扑图记作G＝(V,E)；所述每条巷道数据是一组按照先后顺序排列的三维点坐标，一个含有n个拐点的巷道记为P＝{p0,…,pn}，其中p0为巷道的始端点，pn为巷道的末端点；所有巷道组成的集合记为M＝{P1,P2,P3…}；
依据所有巷道组成的集合M，建立直接拓扑图G1＝(V1,E1)的具体方法：
遍历所述集合M，对于其中的每一条巷道P，依次遍历其拐点，对于拐点pi，在集合V中查找是否有位置重合的节点，如果没有，则向V1中添加一个新节点vi，如果有则记录已有的节点；对于拐点pi+1，在集合V1中查找是否有位置重合的节点，如果没有，则向V1中添加一个新节点vi+1，如果有则记录已有的节点；向边集合E1中添加边ei＝(vi,vi+1)。
可选地，进行简化处理，合并串联分支，并在分支上记录合并的拐点坐标数据，形成所述开环网络拓扑图，包括：
遍历集合V1，对于V1中的每一个点vi，其出度和入度之和deg(vi)＝2，则代表该节点是巷道的中间节点生成得到的，沿着该中间节点向两端搜索节点，直到搜索到deg(vj)≠2，deg(vk)≠2，则把节点vj、vk加入到集合V2，向集合E2中添加边e＝(vj,vk)，同时在边e上关联记录搜索中经过的节点坐标信息，用于后续分支长度及自然风压的计算；完成G1的所有节点遍历后即形成开环网络拓扑图G2，G2＝(V2,E2)。
可选地，在开环网络拓扑图中计算图的连通性，划分成多个子网，包括：
步骤S1：从所述开环网络拓扑图的一个节点出发，对所述开环网络拓扑图进行宽度优先搜索，搜索到的节点和边组成第一个子网；
步骤S2：选一个不属于所述第一个子网的节点重复步骤S1，形成第二个子网，直到所有节点都包含在子网中，搜索结束，得到多个所述子网。
可选地，从所述开环网络拓扑图的一个节点出发，对所述开环网络拓扑图进行宽度优先搜索，搜索到的节点和边组成第一个子网，包括：
从所述集合V2的任一节点vi出发，对所述开环网络拓扑图进行宽度优先搜索，将搜索到的节点加入到集合Vn，将搜索到的边加入到集合En，结束对所述集合V2的搜索后，即形成第一子网Gn＝(Vn,En)。
可选地，选一个不属于所述第一个子网的节点重复步骤S1，形成第二个子网，直到所有节点都包含在子网中，搜索结束，得到多个所述子网，包括：
从所述集合V2中挑选任一节点vi+1，并且继续从所述集合V2搜索所述开环网络拓扑图得到所述第二个子网Gn+1，直至所述集合V2中不存在不属于任何一个子网的节点，搜索结束，得到n个子网，并且所述开环网络拓扑图G2＝Gn∪Gn+1∪Gn+2…。
可选地，针对每个子网，依据所述位置属性和方向属性查找进回风井口，添加虚拟分支连接回风井井口到进风井井口，形成闭环网络图，包括：
在所述连通大气的分支里，依据方向属性，将一个入度为1的节点与其他节点按照特定方向用虚拟分支连接，形成所述闭环网络图；
设位置属性是连通大气的分支集合Vt、内部巷道的分支集合Vi、独头巷道的分支集合Ve、废弃巷道的分支集合Va，则在集合Vt中，依据方向属性找到一个入度indeg(v)＝1的节点，记作vin，以此为始点，创建虚拟分支连接Vt中出度outdeg(v)＝1的节点，以此为末点，创建虚拟分支连接Vt中其他入度indeg(v)＝1的节点，这些虚拟分支构成集合Ev，则一个闭环网络图Gb＝(Vn,En+Ev)；依次构建所有子网形成所述闭环网络图。
可选地，在所述闭环网络图上，查找最小树，圈划回路，建立节点流量守恒方程组和回路能量守恒方程组，得到通风网络解算模型，包括：
假设闭环网络中有m个节点，n条分支，基本回路矩阵记作C，基本关联矩阵记作B，分支流量列矩阵记作Q，则有质量守恒和能量守恒，具体包括：
对于所述质量守恒有：
闭环网络中流入和流出节点的质量守恒，网络m个节点，可建立m-1个质量守恒方程：
上式中：bij表示基本关联矩阵的元素；qi表示第i条分支的流量；
用矩阵表达为式：
BQ＝0
上式中，B表示基本关联矩阵，Q表示流量列矩阵。
对于所述能量守恒有：
对任一个稳态回路，回路中考虑风流方向，巷道消耗的能量和增加的能量相等，记分支沿程阻力为hi，自然风压为hz和机械风压为hf，分支阻力列矩阵H，则有能量守恒方程：
∑hi+hz+hf＝0
对所有的回路，表示成矩阵形式为：
CH＝0
上式中，C表示基本回路矩阵，H表示分支阻力列矩阵。
即共有m-1个节点质量守恒方程，n-m+1个回路能量守恒方程，共计n条分支，n个未知变量，所述通风网络解算模型构建完成，对所述质量守恒方程、所述能量守恒方程组求解得到各个巷道的风量。
本发明实施例提供了一种真实三维巷道空间关系到通风网络解算模型的应用系统，所述应用系统包括：
三维巷道可视化模块，用于对三维真实巷道进行可视化，包括对录入属性数据的着色可视化，通风网络解算结果数据的着色可视化；
属性录入模块，用于对三维巷道属性的录入和查询；
拓扑构建模块，用于根据权利要求1-10任一所述的转换方法得到通风网络解算模型；
通风网络解算模块，用于所述通风网络解算模型对应的拓扑网络图上构建并求解方程组，得到通风网络各个巷道的风量。
本发明提供的真实三维巷道空间关系到通风网络解算模型的转换方法，首先在真实三维巷道上添加位置属性和方向属性；再根据真实三维巷道的真实三维坐标计算巷道的连接关系，建立巷道拓扑关系，形成开环网络拓扑图；之后在开环网络拓扑图中计算图的连通性，划分成多个子网；针对每个子网，依据位置属性和方向属性查找进回风井口，添加虚拟分支连接回风井井口到进风井井口，形成闭环网络图；最后在闭环网络图上，查找最小树，圈划回路，建立节点流量守恒方程组和回路能量守恒方程组，得到通风网络解算模型。
本发明借助地理信息系统对拓扑关系和属性数据的管理功能，很好的解决了这个问题，在基于真实三维坐标巷道模型的三维通风仿真应用中，添加虚拟分支，不会破坏真实巷道拓扑和度量空间关系，用户只需要根据巷道的实际情况进行简单设置，程序即可自动建立用于网络解算的闭环图和数学模型。
相比传统的通风网络解算方法，本发明省略了绘制网络拓扑图和添加大气分支的繁琐过程，使得可以在GIS的三维真实巷道上直接进行通风网络解算。本发明的实现依赖于GIS强大的拓扑关系和属性数据管理功能，充分发挥了GIS在数据管理方面的优势，为通风网络解算提供了高效、准确的解决方案。
附图说明
通过阅读下文优选实施方式的详细描述，各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的，而并不认为是对本发明的限制。而且在整个附图中，用相同的参考符号表示相同的部件。在附图中：
图1是本发明实施例的一种真实三维巷道空间关系到通风网络解算模型的转换方法的流程图；
图2是本发明实施例中真实三维巷道空间关系到通风网络解算模型的转换方法的概述图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。应当理解，此处所描述的具体实施例仅用以解释本发明，仅仅是本发明一部分实施例，而不是全部的实施例，并不用于限定本发明。
参照图1，示出了本发明实施例一种真实三维巷道空间关系到通风网络解算模型的转换方法的流程图，该转换方法包括：
步骤101：在真实三维巷道上添加位置属性和方向属性。
本发明基于地理信息系统，以及地理信息系统对拓扑关系和属性数据的管理功能的基础上实现。首先在真实三维巷道上添加位置属性和方向属性，位置属性和方向属性即为属性数据。
在一种可能的实施例中，位置属性包括：连通大气、内部巷道、独头巷道、废弃巷道；方向属性包括：正向或者反向；其中，真实三维巷道的默认方向是正向，即按照构建巷道的坐标点顺序确定的方向一致，实际风流方向和巷道的默认方向一致的为正向，不一致的为反向。
步骤102：根据真实三维巷道的真实三维坐标计算巷道的连接关系，建立巷道拓扑关系，形成开环网络拓扑图。
添加属性后，再根据真实三维巷道的真实三维坐标计算巷道的连接关系，建立巷道拓扑关系，形成开环网络拓扑图。具体可以包括：
首先依据真实三维巷道的拐点坐标数据，建立直接拓扑图；在进行简化处理，合并串联分支，并在分支上记录合并的拐地坐标数据，即可形成开环网络拓扑图。
其中，建立直接拓扑图的一种较优方式包括：
记网络拓扑图的节点集合为V，边集合为E，形成的拓扑图记作G＝(V,E)；所述每条巷道数据是一组按照先后顺序排列的三维点坐标，一个含有n个拐点的巷道记为P＝{p0,…,pn}，其中p0为巷道的始端点，pn为巷道的末端点；所有巷道组成的集合记为M＝{P1,P2,P3…}。
依据所有巷道组成的集合M，建立直接拓扑图G1＝(V1,E1)的具体方法：
遍历集合M，对于其中的每一条巷道P，依次遍历其拐点，对于拐点pi，在集合V中查找是否有位置重合的节点，如果没有，则向V1中添加一个新节点vi，如果有则记录已有的节点；对于拐点pi+1，在集合V1中查找是否有位置重合的节点，如果没有，则向V1中添加一个新节点vi+1，如果有则记录已有的节点；向边集合E1中添加边ei＝(vi,vi+1)。
进行简化处理，合并串联分支，并在分支上记录合并的拐地坐标数据，形成开环网络拓扑图的一种较优的方式包括：
遍历集合V1，对于V1中的每一个点vi，其出度和入度之和deg(vi)＝2，则代表该节点是巷道的中间节点生成得到的，沿着该中间节点向两端搜索节点，直到搜索到deg(vj)≠2，deg(vk)≠2，则把节点vj、vk加入到集合V2，向集合E2中添加边e＝(vj,vk)，同时在边e上关联记录搜索中经过的节点坐标信息，用于后续分支长度及自然风压的计算；完成G1的所有节点遍历后即形成开环网络拓扑图G2，G2＝(V2,E2)。
步骤103：在开环网络拓扑图中计算图的连通性，划分成多个子网。
在得到开环网络拓扑图后，可以在开环网络拓扑图中计算图的连通性，划分成多个子网。具体包括如下两个步骤：
步骤S1：从开环网络拓扑图的一个节点出发，对开环网络拓扑图进行宽度优先搜索，搜索到的节点和边组成第一个子网；
步骤S2：选一个不属于第一个子网的节点重复步骤S1，形成第二个子网，直到所有节点都包含在子网中，搜索结束，得到多个子网。
其中，步骤S1的一种较优实现方式包括：
从集合V2的任一节点vi出发，对开环网络拓扑图进行宽度优先搜索，将搜索到的节点加入到集合Vn，将搜索到的边加入到集合En，结束对集合V2的搜索后，即形成第一子网Gn＝(Vn,En)。
步骤S2的一种较优实现方式包括：
从集合V2中挑选任一节点vi+1，并且即该任一节点vi+1不属于第一个子网的节点，继续从集合V2搜索开环网络拓扑图得到第二个子网Gn+1，直至集合V2中不存在不属于任何一个子网的节点，搜索结束，得到n个子网，并且开环网络拓扑图G2＝Gn∪Gn+1∪Gn+2…。
步骤104：针对每个子网，依据位置属性和方向属性查找进回风井口，添加虚拟分支连接回风井井口到进风井井口，形成闭环网络图。
得到多个子网后，针对每一个子网均有：依据位置属性和方向属性查找进回风井口，添加虚拟分支连接回风井井口到进风井井口，形成闭环网络图。一种较优的实现方式包括：
在连通大气的分支里，依据方向属性，将一个入度为1的节点与其他节点按照特定方向用虚拟分支连接，形成闭环网络图。
设位置属性是连通大气的分支集合Vt、内部巷道的分支集合Vi、独头巷道的分支集合Ve、废弃巷道的分支集合Va，则在集合Vt中，依据方向属性找到一个入度indeg(v)＝1的节点，记作vin，以此为始点，创建虚拟分支连接Vt中出度outdeg(v)＝1的节点，以此为末点，创建虚拟分支连接Vt中其他入度indeg(v)＝1的节点，这些虚拟分支构成集合Ev，则一个闭环网络图Gb＝(Vn,En+Ev)；依次构建所有子网形成闭环网络图。
步骤105：在闭环网络图上，查找最小树，圈划回路，建立节点流量守恒方程组和回路能量守恒方程组，得到通风网络解算模型。
构建所有子网形成闭环网络图后，在闭环网络图上，查找最小树，圈划回路，建立节点流量守恒方程组和回路能量守恒方程组，得到通风网络解算模型。一种较优的实现方式包括：
假设闭环网络中有m个节点，n条分支，基本回路矩阵记作C，基本关联矩阵记作B，分支流量列矩阵记作Q，则有质量守恒和能量守恒，具体包括：
对于质量守恒有：
闭环网络中流入和流出节点的质量守恒，网络m个节点，可建立m-1个质量守恒方程：
上式中：bij表示基本关联矩阵的元素；qi表示第i条分支的流量；
用矩阵表达为式：
BQ＝0
上式中，B表示基本关联矩阵，Q表示流量列矩阵。
对于能量守恒有：
对任一个稳态回路，回路中考虑风流方向，巷道消耗的能量和增加的能量相等，记分支沿程阻力为hi，自然风压为hz和机械风压为hf，分支阻力列矩阵H，则有能量守恒方程：
∑hi+hz+hf＝0
对所有的回路，表示成矩阵形式为：
CH＝0
上式中，C表示基本回路矩阵，H表示分支阻力列矩阵。
即共有m-1个节点质量守恒方程，n-m+1个回路能量守恒方程，共计n条分支，n个未知变量，通风网络解算模型构建完成，对质量守恒方程、能量守恒方程组求解得到各个巷道的风量。
上述整个转换方法可以用图2所示的概述图概括为：
获取到真实巷道数据，添加巷道位置属性、添加巷道方向属性，再计算巷道拓扑关系，形成开环网络拓扑图。紧接着计算图的连通性，划分子网络。对于每个子网络，即遍历每个子网的过程中，均需要：依据位置属性和方向属性查找进回风井口，在进回风井口间添加大气分支，形成闭环网络图，查找回路，建立节点流量守恒方程组，回路能量守恒方程组，形成网络解算数学模型。
本发明实施例中，基于上述真实三维巷道空间关系到通风网络解算模型的转换方法，还提出一种真实三维巷道空间关系到通风网络解算模型的应用系统，应用系统包括：
三维巷道可视化模块，用于对三维真实巷道进行可视化，包括对录入属性数据的着色可视化，通风网络解算结果数据的着色可视化；
属性录入模块，用于对三维巷道属性的录入和查询；
拓扑构建模块，用于根据上述步骤101～步骤105所述的转换方法得到通风网络解算模型；
通风网络解算模块，用于所述通风网络解算模型对应的拓扑网络图上构建并求解方程组，得到通风网络各个巷道的风量。
综上所述，本发明的提供的真实三维巷道空间关系到通风网络解算模型的转换方法，首先在真实三维巷道上添加位置属性和方向属性；再根据真实三维巷道的真实三维坐标计算巷道的连接关系，建立巷道拓扑关系，形成开环网络拓扑图；之后在开环网络拓扑图中计算图的连通性，划分成多个子网；针对每个子网，依据位置属性和方向属性查找进回风井口，添加虚拟分支连接回风井井口到进风井井口，形成闭环网络图；最后在闭环网络图上，查找最小树，圈划回路，建立节点流量守恒方程组和回路能量守恒方程组，得到通风网络解算模型。
本发明借助地理信息系统对拓扑关系和属性数据的管理功能，很好的解决了这个问题，在基于真实三维坐标巷道模型的三维通风仿真应用中，添加虚拟分支，不会破坏真实巷道拓扑和度量空间关系，用户只需要根据巷道的实际情况进行简单设置，程序即可自动建立用于网络解算的闭环图和数学模型。
相比传统的通风网络解算方法，本发明省略了绘制网络拓扑图和添加大气分支的繁琐过程，使得可以在GIS的三维真实巷道上直接进行通风网络解算。本发明的实现依赖于GIS强大的拓扑关系和属性数据管理功能，充分发挥了GIS在数据管理方面的优势，为通风网络解算提供了高效、准确的解决方案。
尽管已描述了本发明实施例的优选实施例，但本领域内的技术人员一旦得知了基本创造性概念，则可对这些实施例做出另外的变更和修改。所以，所附权利要求意欲解释为包括优选实施例以及落入本发明实施例范围的所有变更和修改。
最后，还需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者终端设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者终端设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者终端设备中还存在另外的相同要素。
上面结合附图对本发明的实施例进行了描述，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，而不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可做出很多形式，这些均属于本发明的保护之内。
Claims (6)
1.一种真实三维巷道空间关系到通风网络解算模型的转换方法，其特征在于，所述转换方法包括：
在所述真实三维巷道上添加位置属性和方向属性；
根据所述真实三维巷道的真实三维坐标计算巷道的连接关系，建立巷道拓扑关系，形成开环网络拓扑图；
从所述开环网络拓扑图中集合V2的任一节点vi出发，对所述开环网络拓扑图进行宽度优先搜索，将搜索到的节点加入到集合Vn，将搜索到的边加入到集合En，结束对所述集合V2的搜索后，即形成第一子网Gn＝(Vn,En)；
从所述集合V2中挑选任一节点vi+1，并且继续从所述集合V2搜索所述开环网络拓扑图得到第二个子网Gn+1，直至所述集合V2中不存在不属于任何一个子网的节点，搜索结束，得到n个子网，并且所述开环网络拓扑图G2＝Gn∪Gn+1∪Gn+2…；
在连通大气的分支里，依据方向属性，将一个入度为1的节点与其他节点按照特定方向用虚拟分支连接，形成闭环网络图；
设位置属性是连通大气的分支集合Vt、内部巷道的分支集合Vi、独头巷道的分支集合Ve、废弃巷道的分支集合Va，则在集合Vt中，依据方向属性找到一个入度indeg(v)＝1的节点，记作vin，以此为始点，创建虚拟分支连接Vt中出度outdeg(v)＝1的节点，以此为末点，创建虚拟分支连接Vt中其他入度indeg(v)＝1的节点，这些虚拟分支构成集合Ev，则一个闭环网络图Gb＝(Vn,En+Ev)；依次构建所有子网形成所述闭环网络图；
假设闭环网络中有m个节点，n条分支，基本回路矩阵记作C，基本关联矩阵记作B，分支流量列矩阵记作Q，则有质量守恒和能量守恒，具体包括：
对于所述质量守恒有：
闭环网络中流入和流出节点的质量守恒，网络m个节点，可建立m-1个质量守恒方程：
上式中：bij表示基本关联矩阵的元素；qi表示第i条分支的流量；
用矩阵表达为式：
BQ＝0
上式中，B表示基本关联矩阵，Q表示流量列矩阵；
对于所述能量守恒有：
对任一个稳态回路，回路中考虑风流方向，巷道消耗的能量和增加的能量相等，记分支沿程阻力为hi，自然风压为hz和机械风压为hf，分支阻力列矩阵H，则有能量守恒方程：
∑hi+hz+hf＝0
对所有的回路，表示成矩阵形式为：
CH＝0
上式中，C表示基本回路矩阵，H表示分支阻力列矩阵；
即共有m-1个节点质量守恒方程，n-m+1个回路能量守恒方程，共计n条分支，n个未知变量，所述通风网络解算模型构建完成，对所述质量守恒方程、所述能量守恒方程组求解得到各个巷道的风量。
2.根据权利要求1所述的转换方法，其特征在于，所述位置属性包括：连通大气、内部巷道、独头巷道、废弃巷道；
所述方向属性包括：正向或者反向；
其中，所述真实三维巷道的默认方向是正向，即按照构建巷道的坐标点顺序确定的方向，实际风流方向和巷道的默认方向一致的为正向，不一致的为反向。
3.根据权利要求1所述的转换方法，其特征在于，根据所述真实三维巷道的真实三维坐标计算巷道的连接关系，建立巷道拓扑关系，形成开环网络拓扑图，包括：
依据所述真实三维巷道的拐点坐标数据，建立直接拓扑图；
进行简化处理，合并串联分支，并在分支上记录合并的拐点坐标数据，形成所述开环网络拓扑图。
4.根据权利要求3所述的转换方法，其特征在于，依据所述真实三维巷道的拐点坐标数据，建立直接拓扑图，包括：
记网络拓扑图的节点集合为V，边集合为E，形成的拓扑图记作G＝(V,E)；每条巷道数据是一组按照先后顺序排列的三维点坐标，一个含有n个拐点的巷道记为P＝{p0,…,pn}，其中p0为巷道的始端点，pn为巷道的末端点；所有巷道组成的集合记为M＝{P1,P2,P3…}；
依据所有巷道组成的集合M，建立直接拓扑图G1＝(V1,E1)的具体方法：
遍历所述集合M，对于其中的每一条巷道P，依次遍历其拐点，对于拐点pi，在集合V中查找是否有位置重合的节点，如果没有，则向V1中添加一个新节点vi，如果有则记录已有的节点；对于拐点pi+1，在集合V1中查找是否有位置重合的节点，如果没有，则向V1中添加一个新节点vi+1，如果有则记录已有的节点；向边集合E1中添加边ei＝(vi,vi+1)。
5.根据权利要求4所述的转换方法，其特征在于，进行简化处理，合并串联分支，并在分支上记录合并的拐点坐标数据，形成所述开环网络拓扑图，包括：
遍历集合V1，对于V1中的每一个点vi，其出度和入度之和deg(vi)＝2，则代表该节点是巷道的中间节点生成得到的，沿着该中间节点向两端搜索节点，直到搜索到deg(vj)≠2，deg(vk)≠2，则把节点vj、vk加入到集合V2，向集合E2中添加边e＝(vj,vk)，同时在边e上关联记录搜索中经过的节点坐标信息，用于后续分支长度及自然风压的计算；完成G1的所有节点遍历后即形成开环网络拓扑图G2，G2＝(V2,E2)。
6.一种真实三维巷道空间关系到通风网络解算模型的应用系统，其特征在于，所述应用系统包括：
三维巷道可视化模块，用于对三维真实巷道进行可视化，包括对录入属性数据的着色可视化，通风网络解算结果数据的着色可视化；
属性录入模块，用于对三维巷道属性的录入和查询；
拓扑构建模块，用于根据权利要求1-5任一所述的转换方法得到通风网络解算模型；
通风网络解算模块，用于所述通风网络解算模型对应的拓扑网络图上构建并求解方程组，得到通风网络各个巷道的风量。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310544884.6A CN116663210B (zh) | 2023-05-15 | 2023-05-15 | 三维巷道空间关系到通风网络解算模型的转换方法和系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310544884.6A CN116663210B (zh) | 2023-05-15 | 2023-05-15 | 三维巷道空间关系到通风网络解算模型的转换方法和系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116663210A CN116663210A (zh) | 2023-08-29 |
| CN116663210B true CN116663210B (zh) | 2024-02-02 |
Family
ID=87723376
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310544884.6A Active CN116663210B (zh) | 2023-05-15 | 2023-05-15 | 三维巷道空间关系到通风网络解算模型的转换方法和系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116663210B (zh) |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103266906A (zh) * | 2013-05-03 | 2013-08-28 | 中国矿业大学 | 一种巷道风阻参数自调整的验证矿井通风系统状态的方法 |
| CN104265349A (zh) * | 2014-08-07 | 2015-01-07 | 玉溪矿业有限公司 | 一种基于dimine软件的矿井通风系统构建方法 |
| CN104564120A (zh) * | 2014-11-11 | 2015-04-29 | 中国矿业大学 | 一种矿井通风系统运行状态控制决策方法 |
| CN104636560A (zh) * | 2015-02-12 | 2015-05-20 | 西安科技大学 | 一种基于子图同构的高风险区域自动识别方法 |
| CN104732065A (zh) * | 2015-02-13 | 2015-06-24 | 中南大学 | 矿井通风网络解算算法及矿井通风网络预测系统 |
| CN104747224A (zh) * | 2015-01-28 | 2015-07-01 | 西安科技大学 | 一种基于矿井通风等效面积的通风难易程度分析方法 |
| CN107277827A (zh) * | 2017-06-30 | 2017-10-20 | 中国矿业大学 | 一种矿井无线传感器网络的网络结构及网络节点部署方法 |
| CN109615697A (zh) * | 2018-11-28 | 2019-04-12 | 中煤科工集团重庆研究院有限公司 | 实现基于WebGL三维巷道地理信息系统的方法 |
| CN113609630A (zh) * | 2021-08-10 | 2021-11-05 | 北京大学 | 一种全矿井自适应的一三维耦合通风网络解算方法 |
| CN114562334A (zh) * | 2022-03-31 | 2022-05-31 | 北京科技大学 | 一种矿井巷道通风状态实时可视化监测预警系统与方法 |
| CN114676490A (zh) * | 2022-05-27 | 2022-06-28 | 太原理工大学 | 一种矿井智能通风网络解算方法 |
| CN115130262A (zh) * | 2022-07-13 | 2022-09-30 | 太原理工大学 | 一种矿井通风系统最大阻力路径的计算方法 |
| CN115796081A (zh) * | 2023-01-04 | 2023-03-14 | 内蒙古科技大学 | 一种复杂通风网络解算方法及解算装置 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9053519B2 (en) * | 2012-02-13 | 2015-06-09 | TaKaDu Ltd. | System and method for analyzing GIS data to improve operation and monitoring of water distribution networks |
-
2023
- 2023-05-15 CN CN202310544884.6A patent/CN116663210B/zh active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103266906A (zh) * | 2013-05-03 | 2013-08-28 | 中国矿业大学 | 一种巷道风阻参数自调整的验证矿井通风系统状态的方法 |
| CN104265349A (zh) * | 2014-08-07 | 2015-01-07 | 玉溪矿业有限公司 | 一种基于dimine软件的矿井通风系统构建方法 |
| CN104564120A (zh) * | 2014-11-11 | 2015-04-29 | 中国矿业大学 | 一种矿井通风系统运行状态控制决策方法 |
| CN104747224A (zh) * | 2015-01-28 | 2015-07-01 | 西安科技大学 | 一种基于矿井通风等效面积的通风难易程度分析方法 |
| CN104636560A (zh) * | 2015-02-12 | 2015-05-20 | 西安科技大学 | 一种基于子图同构的高风险区域自动识别方法 |
| CN104732065A (zh) * | 2015-02-13 | 2015-06-24 | 中南大学 | 矿井通风网络解算算法及矿井通风网络预测系统 |
| CN107277827A (zh) * | 2017-06-30 | 2017-10-20 | 中国矿业大学 | 一种矿井无线传感器网络的网络结构及网络节点部署方法 |
| CN109615697A (zh) * | 2018-11-28 | 2019-04-12 | 中煤科工集团重庆研究院有限公司 | 实现基于WebGL三维巷道地理信息系统的方法 |
| CN113609630A (zh) * | 2021-08-10 | 2021-11-05 | 北京大学 | 一种全矿井自适应的一三维耦合通风网络解算方法 |
| CN114562334A (zh) * | 2022-03-31 | 2022-05-31 | 北京科技大学 | 一种矿井巷道通风状态实时可视化监测预警系统与方法 |
| CN114676490A (zh) * | 2022-05-27 | 2022-06-28 | 太原理工大学 | 一种矿井智能通风网络解算方法 |
| CN115130262A (zh) * | 2022-07-13 | 2022-09-30 | 太原理工大学 | 一种矿井通风系统最大阻力路径的计算方法 |
| CN115796081A (zh) * | 2023-01-04 | 2023-03-14 | 内蒙古科技大学 | 一种复杂通风网络解算方法及解算装置 |
Non-Patent Citations (1)
| Title |
|---|
| 矿井通风系统三维模型的构建与应用;张珂等;工矿自动化;第46卷(第2期);59-64 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116663210A (zh) | 2023-08-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112052547B (zh) | 一种基于人工智能的城市道路网络自动生成方法 | |
| CN104462190B (zh) | 一种基于海量空间轨迹挖掘的在线的位置预测方法 | |
| CN110992473B (zh) | 一种基于车载激光扫描点云的树木枝干建模方法及系统 | |
| CN110057362B (zh) | 有限单元地图的移动机器人路径规划方法 | |
| CN109993064B (zh) | 一种提取图片中道路网络节点间的连接路径的方法 | |
| CN103699654B (zh) | 一种跨比例尺矢量地图水网数据同名目标匹配方法 | |
| CN109145171A (zh) | 一种多尺度地图数据更新方法 | |
| CN111858810B (zh) | 一种面向道路dem构建的建模高程点筛选方法 | |
| CN111667101A (zh) | 融合高分辨率遥感影像与地形的个性化电力野外作业路径规划方法及系统 | |
| Kong et al. | Enhanced facade parsing for street-level images using convolutional neural networks | |
| CN110766792B (zh) | 基于ArcGIS水文分析工具网格演算次序编码方法 | |
| CN107203619A (zh) | 一种复杂网络下的核心子图提取算法 | |
| CN106981194B (zh) | 一种高速路网关键路段的识别方法 | |
| CN107492130B (zh) | 一种基于层面渲染的非连续等值云图生成方法 | |
| CN117875533B (zh) | 矿用安全逃生路径规划方法及系统 | |
| CN116663210B (zh) | 三维巷道空间关系到通风网络解算模型的转换方法和系统 | |
| CN104700453A (zh) | 一种基于gis生成电缆管网三维模型的方法 | |
| Borisov et al. | An automated process of creating 3D city model for monitoring urban infrastructures | |
| CN117473756A (zh) | 一种基于复杂地形电缆铺设优化设计方法 | |
| Uitermark et al. | Semantic and geometric aspects of integrating road networks | |
| CN104766367A (zh) | 一种计算三维模型处理中三维网格拓扑结构图构造方法 | |
| CN114969944A (zh) | 一种高精度道路dem构建方法 | |
| CN113327327A (zh) | 一种面状水系骨架线自动提取方法 | |
| Wu et al. | Extended maptree: a representation of fine-grained topology and spatial hierarchy of bim | |
| Dakowicz et al. | A unified spatial model for GIS |
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 |