CN112836350B - 一种煤矿井下管网瓦斯抽采参数实时解算方法 - Google Patents

一种煤矿井下管网瓦斯抽采参数实时解算方法 Download PDF

Info

Publication number
CN112836350B
CN112836350B CN202110033133.9A CN202110033133A CN112836350B CN 112836350 B CN112836350 B CN 112836350B CN 202110033133 A CN202110033133 A CN 202110033133A CN 112836350 B CN112836350 B CN 112836350B
Authority
CN
China
Prior art keywords
gas
pipe network
extraction
node
working medium
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
Application number
CN202110033133.9A
Other languages
English (en)
Other versions
CN112836350A (zh
Inventor
夏同强
李子龙
高可
孙敦帅
崔昊杰
李小林
魏连江
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Mining and Technology CUMT
Original Assignee
China University of Mining and Technology CUMT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Mining and Technology CUMT filed Critical China University of Mining and Technology CUMT
Priority to CN202110033133.9A priority Critical patent/CN112836350B/zh
Publication of CN112836350A publication Critical patent/CN112836350A/zh
Application granted granted Critical
Publication of CN112836350B publication Critical patent/CN112836350B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/12—Simultaneous equations, e.g. systems of linear equations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Operations Research (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Pipeline Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

一种煤矿井下管网瓦斯抽采参数实时解算方法，首先根据瓦斯抽采系统的节点、管道和钻孔布置建立系统节点的点位分布图；然后基于气体流动质量守恒定律、能量守恒方程建立管网瓦斯抽采参数解算模型；最后将管网传感器实时监测数据、抽采区域煤体基础数据和管网基础参数代入管网瓦斯抽采参数解算模型中进行解算，实现其它节点处混合流量、流速、浓度、压力、密度的解算数据。本发明通过采集有限测点的监测数据，能够实时分析和计算整个管网的瓦斯抽采参数，不仅为瓦斯抽采动态评价和瓦斯灾害预警提供充足的数据基础，而且能够有效减少监测系统的建设成本，提高瓦斯抽采的经济性。

Description

一种煤矿井下管网瓦斯抽采参数实时解算方法
技术领域
本发明涉及一种煤矿井下瓦斯抽采管路抽采参数在线监测、实时解算和数据深度挖掘技术，具体是一种煤矿井下管网瓦斯抽采参数实时解算方法，属于自动化和智能化瓦斯抽采技术领域。
背景技术
我国是世界上最大的煤炭生产国和消费国，约1/3的煤炭年产量来自于高瓦斯、突出矿井，瓦斯灾害治理和资源化利用主要通过施工煤层钻孔进行抽采，而瓦斯抽采系统是建设在煤层瓦斯和抽采管网之间的人工通道。在瓦斯抽采过程中，瓦斯抽采的监测监控对瓦斯抽采动态评价和瓦斯灾害预警有着极其重要的作用。然而，出于经济性考虑，布置过多的传感器会耗费过多的建设成本，有限的传感器监测数据又不足以提供充足的数据来对瓦斯抽采进行动态评价和灾害的预警。那么，如何利用有限测点的监测数据实时地解算出整个管网的抽采数据，已经成为当前急需解决的问题。
发明内容
本发明的目的是提供一种煤矿井下管网瓦斯抽采参数实时解算方法，通过采集有限测点的监测数据，能够实时分析和解算整个管网的瓦斯抽采参数，不仅为瓦斯抽采动态评价和瓦斯灾害预警提供充足的数据基础，而且能够有效减少监测系统的建设成本，提高瓦斯抽采的经济性。
为了实现上述目的，本发明提供一种煤矿井下管网瓦斯抽采参数实时解算方法，包括以下步骤：
①根据瓦斯抽采系统的节点、管道和钻孔布置建立系统节点的点位分布图；
②基于气体流动质量守恒定律、能量守恒方程建立管网瓦斯抽采参数解算模型；
③将管网传感器实时监测数据、抽采区域煤体基础数据和管网基础参数代入管网瓦斯抽采参数解算模型中，实现其它节点处混合流量、流速、浓度、压力、密度的解算数据。
作为本发明的进一步改进，步骤①具体为：依据瓦斯抽采系统的节点、管道和钻孔布置对瓦斯抽采系统节点的点位进行布设，该节点的点位布置是按照管网瓦斯流向和分叉节点的原则形成的点位编号；将抽采泵所在地面作为基准面对每个点位进行标高标定以确定该节点点位的位置势能。
作为本发明的进一步改进，将煤矿瓦斯抽采中所有的钻场(孔)编为一个节点，瓦斯泵出口编为一个节点，所有的三通、四通变径编为其它节点，测点布置在节点位置。
作为本发明的进一步改进，步骤②具体为：基于气体流动质量守恒定律、能量守恒方程以及理想气体状态方程建立同一管网不同点位间的质量与能量的关联性，根据瓦斯抽采管网各传感器的实时监测数据、抽采区域系统基础参数和管网瓦斯抽采参数解算模型建立联立方程，从而实现由单测点监测数据解算管网流动方向任意点位的抽采参数；管网瓦斯抽采参数解算模型主要是利用气体流动质量守恒定律、瓦斯气体质量守恒、能量守恒方程、理想气体状态方程、摩尔质量计算方程建立的联立方程，其进一步包括利用理想气体状态方程计算工质密度，利用流速面积法计算流速，利用混合气体摩尔质量方程计算混合气体摩尔质量。
步骤②中解算模型的建立方法如下：
n-1、n、n+1是布置的三个相邻的节点，
基于气体流动质量守恒定律，建立气体质量守恒方程，如式(1)：
Qn-1ρn-1＝Qnρn+Qn+1ρn+1 (1)
式中：
Qn-1为节点n-1处工质混合流量；
Qn为节点n处工质混合流量；
Qn+1为节点n+1处工质混合流量；
ρn-1为节点n-1处工质密度；
ρn为节点n处工质密度；
ρn+1为节点n+1处工质密度；
其中混合流量Qi(i＝n-1,n,n+1)可通过流速面积法求取，如式(2)，
式中：
vi(i＝n-1,n,n+1)为测点处工质流速；
di(i＝n-1,n,n+1)为测点处管路直径；
工质密度ρi(i＝n-1,n,n+1)可通过理想气体状态方程求取，如式(3)：
式中：
Pi(i＝n-1,n,n+1)为测点处压力；
Mi(i＝n-1,n,n+1)为测点处工质摩尔质量；
R为通用气体常数；
T为气体热力学温度；
其中工质摩尔质量Mi可通过混合气体摩尔质量计算公式计算，如式(4)，
Mi＝16ci+29(1-ci) (4)
式中，
ci(i＝n-1,n,n+1)为测点处瓦斯浓度；
基于气体流动中瓦斯质量守恒，建立瓦斯质量守恒方程，如式(5)，
Qn-1ρn-1cn-1＝Qnρncn+Qn+1ρn+1cn+1 (5)
基于能量守恒原理，建立瓦斯抽采管路能量守恒方程，如式(6)，
式中：
hw1为节点n到节点n-1之间管路的单位质量工质能量损失；
hw2为节点n+1到节点n-1之间管路的单位质量工质能量损失；
α为管网阻力校正系数，取值范围0-2；
Zi(i＝n-1,n,n+1)分别为节点n-1，n，n+1处单位质量工质位置势能(铅垂高度)；其中，管路的单位质量工质能量损失hwj由式(7)计算，
式中：
λj(j＝1,2)分别为节点n到n-1和节点n+1到n-1之间管路的沿程阻力系数；
dj(j＝1,2)分别为抽采支管和抽采干管管径；
Lj(j＝1,2)为点位间距，即节点n到n-1和节点n+1到n-1之间管路长度；
其中，沿程阻力系数λj(j＝1,2)由式(8)确定，
式中：
Δ为管道壁面的粗糙度，一般查表可得；
Rej(j＝1,2)分别为节点n到n-1和节点n+1到n-1之间管路气体流动雷诺数，由式(9)计算，
式中：
vj(j＝1,2)分别为节点n和节点n+1处气体流速；
μj(j＝1,2)分别为节点n到n-1和节点n+1到n-1之间管路内混合气体动力粘度系数，一般取1×10-5；
由伯努利方程得：
Pn＝Pn-1+(Zn-1-Zn)ρng-hw1ρng (10)
Pn+1＝Pn-1+(Zn-1-Zn+1)ρn+1g-hw2ρn+1g (11)
Qk为节点n处混合工质中空气流量，由下式确定：
Qn(1-cn)＝Qk (12)
式中：
cn为节点n处瓦斯浓度；
pn为节点n处工质压力；
β为漏风阻力系数；
Pr为大气压力。
作为本发明的进一步改进，步骤③中管网传感器实时监测数据具体为：在抽采干管上靠近主管的某个点位安装瓦斯浓度传感器、流速传感器、密度传感器、压力传感器，将瓦斯浓度传感器检测到的瓦斯浓度数据、流速传感器检测到的工质流速数据、密度传感器检测到的工质密度数据、压力传感器检测到的工质压力数据传给步骤②中的解算模型，所述抽采区域煤体基础数据包括抽采管路管径、长度、管壁粗糙度。
作为本发明的进一步改进，本发明还包括误差分析过程，具体为：在抽采干管上的部分点位安装瓦斯浓度传感器、流速传感器、密度传感器、压力传感器，将瓦斯浓度传感器检测到的瓦斯浓度数据、流速传感器检测到的工质流速数据、密度传感器检测到的工质密度数据、压力传感器检测到的工质压力数据代入实时解算模型进行解算，基于实测数据和解算数据进行误差分析，用以验证实时解算模型的准确性和修正误差，解算误差包括：解算模型自动提取传感器监测的瓦斯抽采数据，并与同一点位的解算数据进行实时对比，进行误差分析，分析流程如下：若实测数据和解算数据误差在允许误差范围内(-5％～+5％)，则该解算结果符合要求，若实测数据和解算数据误差超过允许误差范围，则管网阻力系数检查校正解算模型，再次对该点位的数据进行解算。
与现有技术相比，本发明的抽采参数实时结算方法，通过建立系统点位分布图，并根据气体流动质量守恒定律、能量守恒方程等建立解算模型，在本发明所建立的一个节点上安装各个传感器，将传感器监测到的实时数据、抽采区域煤体基础数据和管网技术参数代入管网瓦斯抽采参数解算模型进行求解，可以计算出管网其他任一节点的抽采参数，同时通过其他部分节点安装的各个传感器的实测数据可以对该点位实时解算数据进行解算数据的误差校验分析；本发明通过采集有限测点的监测数据，能够实时分析和计算整个管网的瓦斯抽采参数，不仅为瓦斯抽采动态评价和瓦斯灾害预警提供充足的数据基础，而且能够有效减少监测系统的建设成本，提高瓦斯抽采的经济性。
附图说明
图1是本发明瓦斯抽采系统点位布置结构示意图；
图2是本发明瓦斯抽采参数解算流程示意图。
图中：1、抽采泵，2、抽采主管，3、抽采干管，4、抽采支管，5、瓦斯站。
具体实施方式
下面结合附图对本发明作进一步说明。
如图1所示，一种煤矿井下瓦斯抽采管网抽采参数实时解算方法，本发明的抽采管网包括抽采泵1、抽采主管2、抽采干管3、抽采支管4，抽采泵1的进气口连接抽采主管2，抽采泵1的出气口连接瓦斯站5，抽采泵1与瓦斯站5设置在地面上，设置在地面下的抽采主管2连接若干个抽采干管3，每一个抽采干管3连接若干个抽采支管4，由抽采支管4分散出若干个抽采管，抽采管的抽采口嵌入煤层中进行瓦斯抽采作业。
如图2所示，本发明包括以下步骤：
①根据瓦斯抽采系统的节点、管道和钻孔布置建立系统节点的点位分布图；依据瓦斯抽采系统的节点、管道和钻孔布置对瓦斯抽采系统节点的点位进行布设，该节点的点位布置是按照管网瓦斯流向和分叉节点的原则形成的点位编号；将抽采泵所在地面作为基准面对每个点位进行标高标定用以计算该节点处的位置势能；测点布置在节点位置；
②基于气体流动质量守恒定律、能量守恒方程建立管网瓦斯抽采参数解算模型；基于气体流动质量守恒定律、能量守恒方程以及理想气体状态方程建立同一管网不同点位间的质量与能量的关联性，根据瓦斯抽采管网各传感器的实时监测数据、抽采区域系统基础参数和管网瓦斯抽采参数解算模型建立联立方程，从而实现由单测点监测数据解算管网流动方向任一点位的抽采参数；管网瓦斯抽采参数解算模型主要是利用气体流动质量守恒定律、瓦斯气体质量守恒、能量守恒方程、理想气体状态方程、摩尔质量计算方程建立的联立方程，其进一步包括利用理想气体状态方程计算工质密度，利用流速面积法计算流速，利用混合气体摩尔质量方程计算混合气体摩尔质量；
③将管网传感器实时监测数据、抽采区域煤体基础数据和管网基础参数代入管网瓦斯抽采参数解算模型进行求解；在抽采干管3上的点位上安装瓦斯浓度传感器、流速传感器、密度传感器、压力传感器，将瓦斯浓度传感器检测到的瓦斯浓度数据、流速传感器检测到的工质流速数据、密度传感器检测到的工质密度数据、压力传感器检测到的工质压力数据传给步骤②中的解算模型，可以实现其它节点处混合流量、流速、浓度、压力、密度的解算数据，所述抽采区域煤体基础数据包括抽采管路管径、长度、管壁粗糙度。
本发明的解算模型还包括误差分析过程，具体为：在抽采干管3上的部分点位上安装瓦斯浓度传感器、流速传感器、密度传感器、压力传感器，将瓦斯浓度传感器检测到的瓦斯浓度数据、流速传感器检测到的工质流速数据、密度传感器检测到的工质密度数据、压力传感器检测到的压力数据传给误差分析过程，基于实测数据和解算数据进行误差分析，用以验证实时解算模型的准确性和修正误差，解算误差包括：解算模型自动提取传感器监测的瓦斯抽采数据，并与同一点位的解算数据进行实时对比，进行误差分析，分析流程如下：若实测数据和解算数据误差在允许误差范围内(-5％～+5％)，则该解算结果符合要求，若实测数据和解算数据误差超过允许误差范围，则管网阻力系数检查校正解算模型，再次对该点位的数据进行解算。
步骤②中解算模型的建立方法如下：
如图1所示，n-1、n、n+1是布置的三个相邻的节点，
基于气体流动质量守恒定律，建立气体质量守恒方程，如式(1)：
Qn-1ρn-1＝Qnρn+Qn+1ρn+1 (1)
式中：
Qn-1为节点n-1处工质混合流量；
Qn为节点n处工质混合流量；
Qn+1为节点n+1处工质混合流量；
ρn-1为节点n-1处工质密度；
ρn为节点n处工质密度；
ρn+1为节点n+1处工质密度；
其中混合流量Qi(i＝n-1,n,n+1)可通过流速面积法求取，如式(2)，
式中：
vi(i＝n-1,n,n+1)为测点处工质流速；
di(i＝n-1,n,n+1)为测点处管路直径；
工质密度ρi(i＝n-1,n,n+1)可通过理想气体状态方程PV＝mRgT求取，其中Rg＝R/M，如式(3)：
式中：
Pi(i＝n-1,n,n+1)为测点处压力；
Mi(i＝n-1,n,n+1)为测点处工质摩尔质量；
R为通用气体常数；
Rg为气体常数；
T为气体热力学温度；
其中工质摩尔质量Mi可通过混合气体摩尔质量计算公式计算，如式(4)，
Mi＝16ci+29(1-ci) (4)
式中，
ci(i＝n-1,n,n+1)为测点处瓦斯浓度；
基于气体流动中瓦斯质量守恒，建立瓦斯质量守恒方程，如式(5)，
Qn-1ρn-1cn-1＝Qnρncn+Qn+1ρn+1cn+1 (5)
基于能量守恒原理，建立瓦斯抽采管路能量守恒方程，如式(6)，
式中：
hw1为节点n到n-1之间管路的单位质量工质能量损失；
hw2为节点n+1到n-1之间管路的单位质量工质能量损失；
α为管网阻力校正系数，取值范围0-2；
Zi(i＝n-1,n,n+1)分别为节点n-1，n，n+1处单位质量工质位置势能(铅垂高度)；其中，管路的单位质量工质能量损失hwj由式(7)计算，
式中：
λj(j＝1,2)分别为节点n到n-1和节点n+1到n-1之间管路的沿程阻力系数；
dj(j＝1,2)分别为抽采支管和抽采干管管径；
Lj(j＝1,2)为点位间距，即节点n到n-1和节点n+1到n-1之间管路长度；
其中，沿程阻力系数λj(j＝1,2)由式(8)确定，
式中：
Δ为管道壁面的粗糙度，一般查表可得；
Rej(j＝1,2)分别为节点n到n-1和节点n+1到n-1之间管路气体流动雷诺数，由式(9)计算，
式中：
vj(j＝1,2)分别为节点n和节点n+1处气体流速；
μj(j＝1,2)分别为节点n到n-1和节点n+1到n-1之间管路内混合气体动力粘度系数，一般取1×10-5；
由伯努利方程得：
Pn＝Pn-1+(Zn-1-Zn)ρng-hw1ρng (10)
Pn+1＝Pn-1+(Zn-1-Zn+1)ρn+1g-hw2ρn+1g (11)
Qk为节点n处混合工质中空气流量，由下式确定：
Qn(1-cn)＝Qk (12)
式中：
cn为节点n处瓦斯浓度；
pn为节点n处工质压力；
β为漏风阻力系数；
Pr为大气压力。
给出本发明的一个实施例
步骤1，建立完整的矿井抽采系统点位分布图并依次编号，编号方法为：抽采干管3和抽采支管4编制为一类编号，按照风流流动方向，将靠近抽采主管2一端的抽采干管3上的首个点位标记为节点G1，抽采支管4上的点位标记为节点G2，其他点位按从小到大顺序依次进行编号(如图1)；
步骤2，首先，在抽采干管3上节点G1处需布置上述传感器，其他测点可根据抽采系统建设成本任选几个点位布置上述传感器，用于验证实时解算模型的准确性和修正解算误差；
步骤3，根据瓦斯抽采管网抽采系统点位图建立管网瓦斯抽采参数解算模型，系统自动提取上述传感器监测的抽采数据，并传输至解算模型；
步骤4，根据上述联立方程利用软件MathCAD通过有限点位的抽采数据计算其他点位的抽采数据，从而实时解算出整个抽采管网的抽采数据；
步骤5，通过瓦斯抽采管网抽采系统自动提取的监测数据，实时对比上述传感器实测数据和解算数据，计算误差，并且修正解算模型，提高解算结果准确性；
步骤6，在整个抽采管网抽采数据实时解算完成后，系统实时输出实时解算的抽采数据，在用户终端可视化多层次显示；
本实施例中，节点G1的实际抽采负压pG1为20kPa，混合流量QG1为47m3/min，浓度cG1为21％，工质密度ρG1为0.85kg/m3，抽采支管4的管径d1＝dG2为250mm，抽采干管3的管径d2＝dG1＝dG3为500mm，节点G1与节点G3之间管道长度L1为1100m，节点G1与节点G2之间管道长度L2为550m，校正系数α取0.2，为简化计算忽略各个节点的位置势差，即ZG1＝ZG2＝ZG3＝0。
根据上述数学模型建立过程联立方程，包括节点G1，节点G2和节点G3之间的总工质质量守恒、瓦斯质量守恒、能量守恒和漏汽量计算公式；
通过节点G1的监测数据求解节点G2和节点G3的相关参数，将节点G1的数据pG1,QG1,cG1,ρG1代入下列方程，
采用数学通用计算软件MathCad求解上述方程，并经过相关计算解得节点G2处混合流量QG2＝0.43m3/s，流速vG2＝8.68m/s，浓度cG2为22.4％，压力PG2＝80810Pa，密度ρG2＝0.845kg/m3；节点G3处混合流量QG3＝0.36m3/s，流速vG3＝1.85m/s，浓度cG3为22.7％，压力PG3＝80620pa，密度ρG3＝0.842kg/m3。

Claims (4)

1.一种煤矿井下管网瓦斯抽采参数实时解算方法，其特征在于，包括以下步骤：
①根据瓦斯抽采系统的节点、管道和钻孔布置建立系统节点的点位分布图；
②建立管网瓦斯抽采参数解算模型；
③将管网传感器实时监测数据、抽采区域煤体基础数据和管网基础参数代入管网瓦斯抽采参数解算模型中，实现其它节点处混合流量、流速、浓度、压力、密度的解算数据；
步骤①具体为：依据瓦斯抽采系统的节点、管道和钻孔布置对瓦斯抽采系统节点的点位进行布设，节点的点位布置是按照管网瓦斯流向和分叉节点的原则形成点位编号；将抽采泵所在地面作为基准面对每个点位进行标高标定以确定该节点的位置势能；
将煤矿瓦斯抽采中所有的钻孔编为一个节点，瓦斯泵出口编为一个节点，所有的三通、四通变径编为其它节点，测点布置在节点位置；
步骤②中解算模型的建立方法如下：
n-1、n、n+1是布置的三个相邻的节点，
基于气体流动质量守恒定律，建立气体质量守恒方程，如式(1)：
Qn-1ρn-1＝Qnρn+Qn+1ρn+1 (1)
式中：
Qn-1为节点n-1处工质混合流量；
Qn为节点n处工质混合流量；
Qn+1为节点n+1处工质混合流量；
ρn-1为节点n-1处工质密度；
ρn为节点n处工质密度；
ρn+1为节点n+1处工质密度；
其中混合流量Qi，i＝n-1,n,n+1可通过流速面积法求取，如式(2)，
式中：
vi为测点处工质流速；
di为测点处管路直径；
工质密度ρi可通过理想气体状态方程求取，如式(3)：
式中：
Pi为测点处压力；
Mi为测点处工质摩尔质量；
R为通用气体常数；
T为气体热力学温度；
其中工质摩尔质量Mi可通过混合气体摩尔质量计算公式计算，如式(4)，Mi＝16ci+29(1-ci) (4)
式中：
ci为测点处瓦斯浓度；
基于气体流动中瓦斯质量守恒，建立瓦斯质量守恒方程，如式(5)，
Qn-1ρn-1cn-1＝Qnρncn+Qn+1ρn+1cn+1 (5)
基于能量守恒原理，建立瓦斯抽采管路能量守恒方程，如式(6)，
式中：
hw1为节点n到节点n-1之间管路的单位质量工质能量损失；
hw2为节点n+1到节点n-1之间管路的单位质量工质能量损失；
α为管网阻力校正系数，取值范围0-2；
Zi，i＝n-1,n,n+1分别为节点n-1，n，n+1处单位质量工质位置势能；其中，管路的单位质量工质能量损失hwj由式(7)计算，
式中：
λj，j＝1,2分别为节点n到n-1和节点n+1到n-1之间管路的沿程阻力系数；
dj分别为抽采支管和抽采干管管径；
Lj为点位间距，即节点n-1到n和节点n+1到n-1之间管路长度；
其中，沿程阻力系数λj由式(8)确定，
式中：
Δ为管道壁面的粗糙度；
Rej分别为节点n到n-1和节点n+1到n-1之间管路气体流动雷诺数，由式(9)计算，
式中：
vj分别为节点n和节点n+1处气体流速；
μj分别为节点n到n-1和节点n+1到n-1之间管路内混合气体动力粘度系数；由伯努利方程得：
Pn＝Pn-1+(Zn-1-Zn)ρng-hw1ρng (10)
Pn+1＝Pn-1+(Zn-1-Zn+1)ρn+1g-hw2ρn+1g (11)
Qk为节点n处混合工质中空气流量，由下式确定：
Qn(1-cn)＝Qk (12)
式中：
cn为节点n处瓦斯浓度；
pn为节点n处工质压力；
β为漏风阻力系数；
Pr为大气压力。
2.根据权利要求1所述的一种煤矿井下管网瓦斯抽采参数实时解算方法，其特征在于，步骤②具体为：管网瓦斯抽采参数解算模型主要是利用气体流动质量守恒定律、瓦斯气体质量守恒、能量守恒方程、理想气体状态方程、摩尔质量计算方程建立的联立方程，其利用理想气体状态方程计算工质密度，利用流速面积法计算流速，利用混合气体摩尔质量方程计算混合气体摩尔质量。
3.根据权利要求1或2所述的一种煤矿井下管网瓦斯抽采参数实时解算方法，其特征在于，步骤③中管网传感器实时监测数据具体为：在抽采干管上靠近主管的某个点位安装瓦斯浓度传感器、流速传感器、密度传感器、压力传感器，将瓦斯浓度传感器检测到的瓦斯浓度数据、流速传感器检测到的工质流速数据、密度传感器检测到的工质密度数据、压力传感器检测到的工质压力数据传给步骤②中的管网瓦斯抽采参数解算模型，所述抽采区域煤体基础数据包括抽采管路管径、长度、管壁粗糙度。
4.根据权利要求3所述的一种煤矿井下管网瓦斯抽采参数实时解算方法，其特征在于，解算模型还包括误差分析过程，具体为：在抽采干管上的部分点位安装瓦斯浓度传感器、流速传感器、密度传感器、压力传感器，将瓦斯浓度传感器检测到的瓦斯浓度数据、流速传感器检测到的工质流速数据、密度传感器检测到的工质密度数据、压力传感器检测到的工质压力数据代入管网瓦斯抽采参数解算模型进行解算，基于实测数据和解算数据进行误差分析，用以验证实时解算模型的准确性和修正误差，解算误差分析包括：管网瓦斯抽采参数解算模型自动提取传感器监测的瓦斯抽采数据，并与同一点位的解算数据进行实时对比，进行误差分析，分析流程如下：若实测数据和解算数据误差在允许误差范围内，所述误差范围为-5％～+5％，则该解算结果符合要求，若实测数据和解算数据误差超过允许误差范围，则通过管网阻力系数的检查校正管网瓦斯抽采参数解算模型，再次对该点位的数据进行解算。
CN202110033133.9A 2021-01-11 2021-01-11 一种煤矿井下管网瓦斯抽采参数实时解算方法 Active CN112836350B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110033133.9A CN112836350B (zh) 2021-01-11 2021-01-11 一种煤矿井下管网瓦斯抽采参数实时解算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110033133.9A CN112836350B (zh) 2021-01-11 2021-01-11 一种煤矿井下管网瓦斯抽采参数实时解算方法

Publications (2)

Publication Number Publication Date
CN112836350A CN112836350A (zh) 2021-05-25
CN112836350B true CN112836350B (zh) 2024-03-08

Family

ID=75929608

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110033133.9A Active CN112836350B (zh) 2021-01-11 2021-01-11 一种煤矿井下管网瓦斯抽采参数实时解算方法

Country Status (1)

Country Link
CN (1) CN112836350B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115749689B (zh) * 2022-08-30 2024-01-30 中国矿业大学(北京) 一种瓦斯抽采管网智能决策调控方法

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104612635A (zh) * 2014-12-10 2015-05-13 中煤科工集团重庆研究院有限公司 煤层群瓦斯联合抽采达标预评判方法
CN105257335A (zh) * 2015-10-16 2016-01-20 中国矿业大学 地面钻井抽采卸压煤层与采空区瓦斯效果评价及控制方法
CN106223995A (zh) * 2016-08-25 2016-12-14 中煤科工集团重庆研究院有限公司 基于监测历史数据的井下煤层区域瓦斯抽采效果分析方法
CN106640182A (zh) * 2016-12-30 2017-05-10 中国矿业大学(北京) 一种测定钻孔瓦斯微流量方法
CN107967559A (zh) * 2017-11-21 2018-04-27 辽宁工程技术大学 一种瓦斯抽采达标评价可视化管理系统及方法
CN108896438A (zh) * 2018-07-25 2018-11-27 重庆大学 孔底瓦斯含量测量装置和测量方法
CN109598482A (zh) * 2018-11-30 2019-04-09 中国矿业大学 基于网络化模式的煤矿瓦斯抽采设计平台及其服务方法
CN109657310A (zh) * 2018-12-07 2019-04-19 北京理工大学 一种煤矿井下瓦斯爆炸实时动态虚拟现实仿真方法
CN111238583A (zh) * 2018-11-29 2020-06-05 河北金广科技有限公司 一种采用掺混空气检测低流速瓦斯抽采管路流量的方法
CN111577256A (zh) * 2020-05-18 2020-08-25 中国矿业大学 一种穿层钻孔水力冲孔增透效果定量评价方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
MX2017002606A (es) * 2014-08-28 2017-09-19 Total Sa Sistema y metodo para extraer gas de un pozo.

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104612635A (zh) * 2014-12-10 2015-05-13 中煤科工集团重庆研究院有限公司 煤层群瓦斯联合抽采达标预评判方法
CN105257335A (zh) * 2015-10-16 2016-01-20 中国矿业大学 地面钻井抽采卸压煤层与采空区瓦斯效果评价及控制方法
CN106223995A (zh) * 2016-08-25 2016-12-14 中煤科工集团重庆研究院有限公司 基于监测历史数据的井下煤层区域瓦斯抽采效果分析方法
CN106640182A (zh) * 2016-12-30 2017-05-10 中国矿业大学(北京) 一种测定钻孔瓦斯微流量方法
CN107967559A (zh) * 2017-11-21 2018-04-27 辽宁工程技术大学 一种瓦斯抽采达标评价可视化管理系统及方法
CN108896438A (zh) * 2018-07-25 2018-11-27 重庆大学 孔底瓦斯含量测量装置和测量方法
CN111238583A (zh) * 2018-11-29 2020-06-05 河北金广科技有限公司 一种采用掺混空气检测低流速瓦斯抽采管路流量的方法
CN109598482A (zh) * 2018-11-30 2019-04-09 中国矿业大学 基于网络化模式的煤矿瓦斯抽采设计平台及其服务方法
CN109657310A (zh) * 2018-12-07 2019-04-19 北京理工大学 一种煤矿井下瓦斯爆炸实时动态虚拟现实仿真方法
CN111577256A (zh) * 2020-05-18 2020-08-25 中国矿业大学 一种穿层钻孔水力冲孔增透效果定量评价方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
煤矿瓦斯抽采效果定量评价理论模型与可视化软件开发;高可;《中国优秀硕士学位论文全文数据库》(第2期);B021-888 *
瓦斯抽采过程中参数演化及钻孔优化数值模拟研究;陈月霞;《中国优秀博士学位论文全文数据库》(第4期);B021-27 *

Also Published As

Publication number Publication date
CN112836350A (zh) 2021-05-25

Similar Documents

Publication Publication Date Title
CN105354257B (zh) 一种管道本体数据的里程与地理空间坐标拟合方法
CN114741831B (zh) 一种供水管网水力压强预警分布的模拟方法及系统
CN114239194A (zh) 一种大水量输供水管网漏损分析和漏点定位方法
CN112836350B (zh) 一种煤矿井下管网瓦斯抽采参数实时解算方法
CN107859527A (zh) 矿井测压测风装置及其测压测风的方法
CN103983313A (zh) 一种管网阻力系数确定方法及管网风量测量方法
CN106052773A (zh) 一种基于多孔探针的烟道排气流量校准装置及方法
CN107083983A (zh) 一种集约化矿井通风阻力测算方法
CN107013812A (zh) 一种三场耦合管道泄漏监测方法
CN114048694A (zh) 一种基于物联网和水力模型的漏损定位方法
CN106709995A (zh) 核电厂闭合管线三维测量与建模方法
CN206095879U (zh) 一种二维岩石试样裂隙网络渗流试验中的流量采集系统
CN113585982B (zh) 一种带压力流量温度检测的油管
CN213656332U (zh) 一种供水管道泄漏检测装置
CN116164241A (zh) 一种瓦斯抽采管网泄漏故障智能检测方法
CN116008585A (zh) 一种巷道风量确定方法
CN211476790U (zh) 一种卧式高压加热器换热效率在线监测系统
CN113901628B (zh) 用于热油管道仿真的方法
CN110987544B (zh) 单管式多点位火电厂矩形烟道烟气采样装置及采样方法
CN102419161B (zh) 一种用于充满不可压缩流体管道的管道长度测量装置
CN114370984B (zh) 一种长距离供水管道水压试验渗漏位置的智能识别方法
CN207180780U (zh) 地下用高精度气体超声流量测量装置
CN111188999A (zh) 一种供水管道泄漏检测方法及系统
CN110985105A (zh) 一种煤层瓦斯抽采半径测量系统及测量方法
CN110332962B (zh) 一种供水渠道分水口流量测算方法

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