CN116357374A - 一种多风机联合运转干扰下采空区漏风的智能调控方法 - Google Patents

一种多风机联合运转干扰下采空区漏风的智能调控方法 Download PDF

Info

Publication number
CN116357374A
CN116357374A CN202310340916.0A CN202310340916A CN116357374A CN 116357374 A CN116357374 A CN 116357374A CN 202310340916 A CN202310340916 A CN 202310340916A CN 116357374 A CN116357374 A CN 116357374A
Authority
CN
China
Prior art keywords
air
air leakage
difference
goaf
upper computer
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.)
Pending
Application number
CN202310340916.0A
Other languages
English (en)
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.)
Taiyuan University of Technology
Original Assignee
Taiyuan University of Technology
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 Taiyuan University of Technology filed Critical Taiyuan University of Technology
Priority to CN202310340916.0A priority Critical patent/CN116357374A/zh
Publication of CN116357374A publication Critical patent/CN116357374A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21FSAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
    • E21F1/00Ventilation of mines or tunnels; Distribution of ventilating currents
    • E21F1/10Air doors
    • E21F1/12Devices for automatically opening air doors
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21FSAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
    • E21F1/00Ventilation of mines or tunnels; Distribution of ventilating currents
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02BCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
    • Y02B30/00Energy efficient heating, ventilation or air conditioning [HVAC]
    • Y02B30/70Efficient control or regulation technologies, e.g. for control of refrigerant flow, motor or heating

Landscapes

  • Engineering & Computer Science (AREA)
  • Mining & Mineral Resources (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geology (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明公开关于一种多风机联合运转干扰下采空区漏风的智能调控方法。该方法基于控制箱接收到传感器实时监测到的风量、瓦斯浓度数据并反馈到上位机中,上位机实时计算工作面向采空区的实际漏风量,由控制箱中设有变论域模糊PID系统的PLC做出相应的伺服电机信号,伺服电机控制风窗面积进行调控风量,上位机根据调控后的漏风量做出进一步指令,直到达到设定的漏风量要求。本方法可以根据需求智能调控工作面向采空区的漏风量,最终有效治理在多风机联合运转干扰下造成的采空区漏风问题,消除漏风带来的隐患,保证矿井安全生产。

Description

一种多风机联合运转干扰下采空区漏风的智能调控方法
技术领域
本发明涉及煤炭开采技术领域,具体涉及一种多风机联合运转干扰下采空区漏风的智能调控方法。
背景技术
由于原来矿井开采时存在私挖乱采和不规范性等现象,造成资源整合后的矿井采空区相互连通。当相邻矿井开采同一煤层的邻近工作面时,由于其上覆近距离煤层开采已形成大面积的采空区,随着该层工作面的开采和后方顶板的垮落,导致现采空区与上覆煤层采空区连通;同时受多风机负压联合干扰影响,导致相邻矿井采空区之间存在压能差。由于此压能差的存在极易造成工作面风量流向采空区,形成采空区漏风,导致通风系统不稳定、工作面风量不足、采空区遗煤自燃、回风隅角瓦斯容易积聚等现象,给矿井安全生产带来隐患。
目前治理采空区漏风常采用的方法有:①注浆堵漏:浆液注入或者直接喷洒在采空区内,浆液能够渗透到煤和矸石的缝隙中,缩小煤与氧气接触的反应面,降低采空区内的孔隙率,增强碎煤和矸石的胶结性,提供区域内的气密性,从而减少采空区的漏风。②高分子材料填充:向回风隅角采空区漏风通道或工作面液压支架后新形成采空区漏风通道注入高分子材料,使其与产生裂隙的煤岩体紧密结合,达到减少漏风效果。③挡风墙技术:在进、回风隅角处垒砌墙,中间填砂,有的采用填河沙、黄土等不燃性材料装入编织袋堆砌构筑堵漏风墙,尽量使两个隅角处封严,以减小向采空区漏风。
上述几种方法对减少工作面向采空区漏风有一定效果,但需随着工作面的推进及时进行施工,耗费大量的人力物力,智能化程度也低,并且对存在大面积的采空区漏风治理效果不明显,再加上矿井负压受气候条件变化影响较大,在不同的时间段会造成井下采空区漏风量不一样,在多风机联合运转干扰下采空区漏风情况就更加复杂,增加了采空区漏风治理难度。因此提供一种多风机联合运转干扰下采空区漏风的智能调控方法成为必要。
发明内容
当相邻矿井开采同一煤层的邻近工作面时,由于其上覆近距离煤层开采已形成大面积的采空区,随着该层工作面的开采和后方顶板的垮落,导致现采空区与上覆煤层采空区连通;同时受多风机负压联合干扰影响,导致相邻矿井采空区之间存在漏风现象,为解决此问题本发明提出一种多风机联合运转干扰下采空区漏风的智能调控方法,减少采空区漏风的方法。
为实现上述目的,本发明采用如下技术方案:
一种多风机联合运转干扰下采空区漏风的智能调控方法,利用传感器实时监测风量、瓦斯浓度,控制箱接收传感器的监测数据并反馈至上位机,上位机实时计算工作面向采空区的实际漏风量,控制箱中的伺服电机控制风窗面积进行调控风量,上位机根据调控后的漏风量做出进一步指令,直到达到设定的漏风量要求,具体步骤如下:
步骤1.在采空区漏风区域工作面两侧的进风巷与回风巷安装超声波风速传感器,进风巷安装智能风门-风窗,在智能风门-风窗附近安装配套控制箱;
步骤2.上位机根据工作面两侧进风巷与回风巷风量实时计算实际漏风量Q,根据工况要求在上位机中设定工作面漏风量Q,上位机实时计算实际漏风量Q与设定工作面漏风量Q的差值e,以及差值变化率ec,并将差值e和差值变化率ec传输到控制箱中;
步骤3.将上位机计算所得的实际漏风量Q与设定工作面漏风量Q的差值e和差值变化率ec作为输入变量传输给控制箱中的变论域模糊PID调控系统;经过变论域模块得出伸缩因子,伸缩因子对初始比例因子和初始量化因子进行调控;经过模糊调节模块得出PID的三个参数KI、KP、KD;将漏风量差值e、差值变化率ec、KI、KP、KD作为PID控制器模块的输入变量,PID控制器模块输出伺服电机的控制信号;
步骤4.伺服电机根据接收的控制信号调节风门-风窗过风面积,改变进风量,平衡所监测的采空区与周边采空区之间的压差,直到实际漏风量Q低于或者相等设定漏风量Q;若上位机计算得出漏风量未达到所设定的数值,将重新进入步骤3中开始计算,直到实际漏风量低于或等于设定漏风量。
进一步,控制箱为KXJ127矿用隔爆兼本安型PLC控制箱,在PLC控制箱中设置变论域模糊PID调控系统,控制箱的信号接收端连接用于检测风量的超声波风速传感器、采煤工作面回风隅角甲烷传感器T0和两个采煤工作面回风巷甲烷传感器;控制箱的信号输出端连接伺服电机的控制端;控制箱与上位机、伺服电机信号双向传输。
进一步,控制箱的变论域模糊PID调控系统包括变论域模块、模糊调节模块和PID控制器模块。
进一步,超声波风速传感器设置在采空区漏风工作面的进风巷和回风巷,距离工作面切眼60m开始向采空区方向等间距安装一组用于监测所漏风工作面进风巷和回风巷风量的超声波风速传感器,且进风巷和回风巷两侧的超声波风速传感器的数量均不少于3组。
进一步,步骤2中,上位机实时计算实际漏风量Q与设定工作面漏风量Q的差值e与差值变化率ec,计算方法如下:
步骤2.1设进风巷超声波风速传感器监测的风量为Q进i,表示第i个传感器监测的风量,同理回风巷超声波传感器监测的风量为Q回i,进风巷风量为
Figure BDA0004158046780000031
其中n表示传感器的数量,同理计算/>
Figure BDA0004158046780000032
工作面实际漏风量Q=Q-Q
步骤2.2实际漏风量与设定漏风量差值e=Q-Q,差值变化率
Figure BDA0004158046780000033
进一步,上位机中设定工作面最低需风量,将回风巷风量视为工作面风量,若工作面风量低于工作面最低需风量需重新设定漏风量Q
进一步,变论域模糊PID调控系统具体包括以下内容:
S1,设输入变量差值e与差值变化率ec初始论域分别为[-e0,e0]、[-ec0,ec0],输出变量ΔKP、ΔKI、ΔKD的初始论域分别为[-Kp0,Kp0]、[-Ki0,Ki0]、[-Kd0,Kd0],并设输入变量差值e与差值变化率ec以及输出变量ΔKP、ΔKI、ΔKD的模糊论域均为[-M,M];
S2,将实际漏风量Q与设定工作面漏风量Q的差值e与差值变化率ec作为变论域模块的输入变量,采用伸缩因子函数模型得出输入变量差值e的伸缩因子αe、差值变化率ec的伸缩因子αec和三个输出变量ΔKP、ΔKI、ΔKD的伸缩因子βP、βI、βD;其中伸缩因子函数模型如下:
输入伸缩因子:
Figure BDA0004158046780000034
输出伸缩因子:
Figure BDA0004158046780000041
式中λ∈(0,1),k>0;KI为比例因子;Pi为输入权重系数;ei(τ)为输入变量误差;β(0)为输出论域伸缩因子初值;
S3,根据S1中所设定的初始论域和模糊论域得出初始量化因子Ke0、Kec0以及初始比例因子LP0、LI0、LD0,即
Figure BDA0004158046780000042
Figure BDA0004158046780000043
将S2中所得的伸缩因子αe、αec对初始量化因子Ke0、Kec0进行调整得出量化因子Ke与Kec,即Ke=αe·Ke0、Kec=αec·Kec0;差值e与差值变化率ec作为模糊调节模块输入变量且经过量化因子作用模糊化处理得到模糊量e*与ec *,即e*=Ke·e、ec *=Kec·ec;参照专家领域知识设置模糊集和模糊规则以及采用重心法对模糊量e*与ec *进行模糊处理得到输出值ΔKP *、ΔKI *、ΔKD *;利用所输出变量的伸缩因子βP、βI、βD对初始比例因子进行调节得到比例因子LP、LI、LD,即LP=βP·LP0、LI=βI·LI0、LD=βD·LD0;将输出值ΔKP *、ΔKI*、ΔKD*对三个比例因子对应进行比例调节得到参数调整量ΔKP、ΔKI、ΔKD,即ΔKP=ΔKP *·LP、ΔKI=ΔKI *·LI、ΔKD=ΔKD *·LD
S4,通过PID控制器参数整定方法得出初始值KP0、KI0、KD0,根据上述得出的调整量ΔKP、ΔKI、ΔKD对初始值进行调整,获得PID控制的3个参数KP=KP0+ΔKP、KI=KI0+ΔKI、KD=KD0+ΔKD;将差值e与差值变化率ec以及3个参数KP、KI、KD作为PID控制器模块的输入变量,PID控制器模块输出控制信号传输到伺服电机;其中关于PID控制算法的离散化公式为:
Figure BDA0004158046780000044
式中un为调节n次后输出的结果,en为第n次实际漏风量Q与设定漏风量Q差值,Ts为调控时间,KP为比例常数、TI积分时间常数、TD微分时间常数。
进一步,上位机中设定工作面最低需风量,若调节后工作面风量低于最低需风量,重新设定漏风量Q;上位机同时实时监测回风隅角T0与回风巷的瓦斯浓度数据,当瓦斯浓度数据超限时上位机开启气体浓度警告,发出警告指令到控制箱进行语音警告,并控制箱输出控制信号到伺服电机,控制智能风门-风窗开度恢复到最大状态,待瓦斯浓度正常后,根据实时监测的风量,再重新设定漏风量Q并保证瓦斯浓度不再超限,继续进行调节。
综上,发明具有以下有益效果:
本发明主要能够有效平衡本工作面采空区与周边工作面采空区存在的压能差,降低采空区漏风对通风系统的影响,同时能随时监测工作面两侧进风巷和回风巷风量情况以及采空区漏风情况,并且能够自动调节进风巷风量,有效控制工作面漏风量,利用变论域模糊PID控制方法实现对风窗的优化控制,实现了调整时间短、稳态精度高、稳态误差小、超调量小等优点,实现了采空区漏风量实时智能化调控。
附图说明
图1为本方法实例说明图
图2为本方法的结构示意图
图3变论域模糊PID关系图
图中:1-智能风门-风窗;2-控制箱;3-超声波风速传感器;4-A煤业进风巷;5-A煤业回风巷;6-A煤业采空区;7-上覆采空区;8-B煤业采空区;9-伺服电机。
具体实施案例
下面结合附图和实施例对本发明进行清楚、完整地说明,显然,所描述的实施案例仅是本发明一部分实施案例,而不是全部的实施例。基于本发明的实施案例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施案例,都属于本发明的保护范围。
如图1~3所示,本发明公开了一种多风机联合运转干扰下采空区漏风的智能调控方法,其特征在于:利用传感器实时监测风量、瓦斯浓度,控制箱接收传感器的监测数据并反馈至上位机,上位机实时计算工作面向采空区的实际漏风量,控制箱中的伺服电机控制风窗面积进行调控风量,上位机根据调控后的漏风量做出进一步指令,直到达到设定的漏风量要求,具体步骤如下:
步骤1,在采空区漏风所在区域工作面两侧的进风巷与回风巷布置超声波风速传感器。进风巷安装智能风门-风窗,在智能风门-风窗附近安装配套控制箱。上述传感器实时监测工作面两侧进风巷和回风巷的风量。所述控制箱分别与上位机和超声波风速传感器连接,且将采煤工作面回风隅角甲烷传感器T0和采煤工作面回风巷甲烷传感器T1、T2接入控制箱。
步骤2,上位机根据工作面两侧进风巷与回风巷风量实时计算实际漏风量Q,根据工况要求在上位机中设定工作面漏风量Q,上位机实时计算实际漏风量Q与设定工作面漏风量Q的差值e,以及差值变化率ec,并将数据传输到控制箱中。
步骤3,将上位机计算所得的实际漏风量Q与设定工作面漏风量Q的差值e和差值变化率ec作为输入变量传输给控制箱中的变论域模糊PID调控系统;经过变论域模块得出伸缩因子,伸缩因子对初始比例因子和初始量化因子进行调控;经过模糊调节模块得出PID的三个参数KI、KP、KD;将漏风量差值e、差值变化率ec、KI、KP、KD作为PID控制器模块的输入变量,PID控制器模块输出伺服电机的控制信号。
步骤4,伺服电机根据信号智能调节风门-风窗过风面积,改变进风量,平衡本采空区与周边采空区之间的压差,直到实际漏风量Q低于或者相等设定漏风量Q。上位机中设定工作面最低需风量,将回风巷风量视为工作面风量,若工作面风量低于工作面最低需风量需重新设定漏风量Q
步骤1中,智能风门-风窗安设于进风巷中,前提条件为本工作面采空区向周边采空区漏风。
步骤1中,智能风门-风窗设有伺服电机控制风窗开度大小,其中伺服电机与控制箱连接并接收控制箱的调节信号;控制箱与上位机、伺服电机均可双向传递信号。控制箱接收传感器监测的数据且反馈到上位机以及接收上位机传输的指令和数据。其中超声波风速传感器布置于所在采空区漏风工作面的进风巷和回风巷,距离工作面切眼60m开始向采空区方向等间距安装一组超声波传感器且两侧的超声波风速传感器的数量均不少于3组,在本实施例中,传感器数量为3组,用于监测所漏风工作面进风巷和回风巷的风量。
所述步骤2中,上位机实时监测进风巷和回风巷的超声波风速传感器的数据且实时计算进风巷和回风巷风量取其平均值作为实际风量Q,根据进风巷和回风巷风量推算实际漏风量Q,计算方法如下:
步骤2.1设进风巷3组超声波风速传感器监测的风量为Q进i,表示第i个传感器监测的风量,同理回风巷超声波传感器监测的风量为Q回i,在本实施例中,分别为Q进1、Q进2、Q进3,同理回风巷3组超声波传感器监测的风量为Q回1、Q回2、Q回3,即:
进风巷风量为
Figure BDA0004158046780000071
同理计算/>
Figure BDA0004158046780000072
工作面实际漏风量Q=Q-Q
步骤2.2实际漏风量与设定漏风量差值e=Q-Q,差值变化率
Figure BDA0004158046780000073
控制箱为KXJ127矿用隔爆兼本安型PLC控制箱,在PLC控制箱设有变论域模糊PID调控系统的PLC。变论域模糊PID调控系统包括变论域模块、模糊调节模块和PID控制器模块。变论域模糊PID调控系统具体包括以下内容:
S1,设输入变量差值e与差值变化率ec初始论域分别为[-e0,e0]、[-ec0,ec0],输出变量ΔKP、ΔKI、ΔKD的初始论域分别为[-Kp0,Kp0]、[-Ki0,Ki0]、[-Kd0,Kd0],并设输入变量差值e与差值变化率ec以及输出变量ΔKP、ΔKI、ΔKD的模糊论域均为[-M,M]。
S2,将实际漏风量Q与设定工作面漏风量Q的差值e与差值变化率ec作为变论域模块的输入变量,在变论域模块中经过伸缩因子函数模型调整得出输入变量差值e与差值变化率ec的伸缩因子αe、αec和输出变量ΔKP、ΔKI、ΔKD的伸缩因子βP、βI、βD。其中伸缩因子函数模型如下:
输入伸缩因子:α(x)=1-λe-kx2
输出伸缩因子:
Figure BDA0004158046780000074
式中λ∈(0,1),k>0;KI为比例因子;Pi为输入权重系数;ei(τ)为输入变量误差;β(0)为输出论域伸缩因子初值。
S3,根据S1中所设定的初始论域和模糊论域得出初始量化因子Ke0、Kec0以及初始比例因子LP0、LI0、LD0,即
Figure BDA0004158046780000075
Figure BDA0004158046780000081
将所得的输入变量差值e与差值变化率ec的伸缩因子αe、αec对初始量化因子Ke0、Kec0进行调整得到量化因子Ke与Kec,即Ke=αe·Ke0、Kec=αec·Kec0;差值e与差值变化率ec作为模糊调节模块输入变量且经过量化因子作用模糊化处理得到模糊量e*与ec *,即e*=Ke·e、ec *=Kec·ec;参照专家领域知识设置模糊集和模糊规则以及采用重心法对模糊量e*与ec *进行模糊处理得到输出值ΔKP *、ΔKI *、ΔKD *;利用所输出变量的伸缩因子βP、βI、βD对初始比例因子进行调节得到比例因子LP、LI、LD,即LP=βP·LP0、LI=βI·LI0、LD=βD·LD0;将输出值ΔKP *、ΔKI*、ΔKD*对其进行比例调节得到参数调整量ΔKP、ΔKI、ΔKD,即ΔKP=ΔKP *·LP、ΔKI=ΔKI *·LI、ΔKD=ΔKD *·LD
S4,通过PID控制器参数整定方法得出初始值KP0、KI0、KD0,根据上述得出的调整量ΔKP、ΔKI、ΔKD对初始值进行调整,获得PID的3个参数KP=KP0+ΔKP、KI=KI0+ΔKI、KD=KD0+ΔKD。将差值e与差值变化率ec以及3个参数KP、KI、KD作为PID控制器模块的输入变量,PID控制器模块输出控制信号。其中关于PID控制算法的离散化公式为:
Figure BDA0004158046780000082
式中un为调节n次后输出的结果,en为第n次实际漏风量Q与设定漏风量Q差值,Ts为调控时间,KP为比例常数、TI积分时间常数、TD微分时间常数。
所述步骤4中,伺服电机接收到控制箱的调节信号对智能风门-风窗开度进行调节,平衡本采空区与周边采空区之间存在的压差,以此达到控制漏风量的目的。若上位机计算得出漏风量未达到所设定的数值,将重新进入步骤二中开始计算,直到实际漏风量低于或等于设定漏风量。上位机中设定工作面最低需风量,若调节后工作面风量低于最低需风量,将需重新设定漏风量Q。在本实施例中,以检测甲烷浓度为例,上位机同时实时监测回风隅角T0与回风巷T1、T2的甲烷浓度数据,当上传甲烷浓度数据超限时上位机会进行自动警告,并且发出指令到控制箱,将智能风门-风窗的开度恢复到原状态,待瓦斯浓度正常后,再重新设定漏风量Q,继续进行调节。
实施例1:布置相关装置
在采空区漏风所在区域工作面两侧的进风巷与回风巷布置超声波风速传感器,在本实施例中,距离工作面切眼60m开始向采空区方向每隔10m布置一组超声波传感器且每侧各布置3组,用于监测所漏风工作面进风巷和回风巷风量。进风巷安装智能风门-风窗,智能风门-风窗设有伺服电机控制风窗开度大小。在智能风门-风窗附近安装配套控制箱,其中伺服电机与控制箱连接并接收控制箱的调节信号。所述控制箱分别与上位机和超声波风速传感器连接,且将采煤工作面回风隅角甲烷传感器T0和采煤工作面回风巷甲烷传感器T1、T2接入控制箱,控制箱接收传感器监测的数据且反馈到上位机以及接收上位机传输的指令和数据。
步骤二:设定工作面漏风量Q,计算变论域模糊PID调控系统的输入变量
上位机实时监测进风巷和回风巷的超声波风速传感器的数据且实时计算进风巷和回风巷风量取其平均值作为实际风量,根据进风巷和回风巷风量推算实际漏风量Q;根据要求在上位机中设定工作面漏风量Q,上位机实时计算实际漏风量Q与设定工作面漏风量Q的差值e与差值变化率ec,并将数据传输到控制箱中,计算方法如下:
①设进风巷超声波风速传感器监测的风量分别为Q进1、Q进2、Q进3,同理回风巷超声波传感器监测的风量为Q回1、Q回2、Q回3
②进风巷风量为
Figure BDA0004158046780000091
同理计算/>
Figure BDA0004158046780000092
工作面实际漏风量Q=Q-Q
③实际漏风量与设定漏风量差值e=Q-Q,差值变化率
Figure BDA0004158046780000093
步骤三:变论域模糊PID系统输出控制信号
将实际漏风量Q与设定工作面漏风量Q的差值e与差值变化率ec作为输入变量传输给控制箱中的变论域模糊PID调控系统。经过变论域模块得出伸缩因子,伸缩因子对初始比例因子和初始量化因子进行调控;经过模糊调节模块得出PID的三个参数KI、KP、KD。将漏风量差值e、差值变化率ec、KI、KP、KD作为PID控制器模块的输入变量,PID控制器模块输出伺服电机的控制信号。
步骤四:调节风窗面积使实际漏风量Q达到设定漏风量Q
控制箱接收到控制信号并传输给伺服电机,伺服电机根据信号调节智能风门-风窗开度的大小,通过调节风窗面积大小在一定程度改变进风量,平衡本采空区与周边采空区之间的压差。上位机实时监测实际漏风量Q,若上位机得出实际漏风量Q未达到所设定漏风量Q,将重新进入步骤二中开始计算,直到实际漏风量达到设定漏风量。上位机中设定工作面最低需风量,若调节后工作面风量低于最低需风量,将需重新设定漏风量Q。上位机同时实时监测回风隅角T0与回风巷T1、T2的甲烷浓度数据,当上传甲烷浓度数据超限时上位机会进行自动警告,并且发出指令到控制箱,将智能风门-风窗的开度恢复到原状态,待瓦斯浓度正常后,再重新设定漏风量Q,继续进行调节。
实施例2:A、B煤业均为资源整合相邻矿井,A、B煤业所处区域的2号煤层已开采完毕并形成大面积的连通采空区,现均开采3号煤层,2号煤层距离3号煤层平均间距为7m。其中A煤业3304工作面与B煤业30303工作面在同时开采3号煤层时受其顶板垮落影响导致采空区与上覆近距离老采空区导通,形成漏风通道。经过现场调研和测试,A煤业3304工作面向B煤业方向(采空区)漏风,造成A煤业3304工作面风量不足,并容易引发采空区自然发火,具有严重安全隐患。
为了控制A煤业3304工作面漏风情况,提出一种多风机联合运转干扰下采空区漏风的智能调控方法具体步骤如下:
步骤一:布置相关装置
在A煤业3304工作面两侧的3304进风巷4与3304回风巷5向工作面切眼60m方向开始每隔10m布置一组超声波风速传感器3,实时监测工作面两侧进风巷和回风巷的风量。将A煤业3304工作面回风隅角甲烷传感器T0和3304回风巷5甲烷传感器T1、T2接入控制箱2。经过现场勘测A煤业3304工作面向采空区漏风,为了控制漏风则在3304进风巷4布置智能风门-风窗1。在智能风门-风窗1附近安装配套控制箱2,控制箱2连接控制风窗开度大小的伺服电机9、超声波风速传感器3和上位机。控制箱2将接收传感器所监测的风量、瓦斯的数据传输上位机。
步骤二:设定工作面漏风量Q,计算变论域模糊PID调控系统的输入变量。
上位机根据3304进风巷4与3304回风巷5布置的超声波风速传感器3的数据实时计算实际漏风量Q。监测到3304进风巷4风量为1408.37m3/min,3304回风巷5风量为890.7m3/min,即A煤业3304工作面实际漏风量Q为517.67m3/min。上位机监测瓦斯浓度T0、T1、T2,以防止瓦斯浓度过高造成安全事故。根据A煤业实际需求将漏风量设定到200m3/min时符合要求,即将其在上位机中设定为工作面漏风量Q,上位机计算计算实际漏风量Q与设定工作面漏风量Q的差值e与差值变化率ec作为变论域模糊PID调控系统输入变量。
步骤三:变论域模糊PID系统输出控制信号
将步骤二中上位机计算的变论域模糊PID调控系统输入变量传输到控制箱2。经过变论域模块得出伸缩因子,伸缩因子进行调控初始比例因子和初始量化因子,经过模糊调节模块得出PID的三个参数KI、KP、KD,将漏风量差值e、差值变化率ec、KI、KP、KD作为PID控制器模块的输入,PID控制器模块输出伺服电机的控制信号,控制箱2接收到控制信号并传输给伺服电机9。
步骤四:调节风窗面积使实际漏风量Q达到设定漏风量Q
伺服电机9根据控制箱2传递的信号调节风窗开度的大小,减小风窗面积,平衡B煤业30303工作面采空区8与A煤业3304工作面采空区6压差,在一定程度减少3304进风巷4进风量。经过多次计算和调节最终将A煤业3304进风巷4风量为890.21m3/min、3304回风巷5风量为767.32m3/min,则A煤业3304工作面漏风量降低到122.9m3/min,且工作面需风量高于最低需风量以及瓦斯浓度T0、T1、T2并无超限,达到本次控制漏风量的目的。
以上所述仅是本发明的优选实施方式,本发明的保护范围并不仅局限于上述实施例,凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理前提下的若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (8)

1.一种多风机联合运转干扰下采空区漏风的智能调控方法,其特征在于:利用传感器实时监测风量、瓦斯浓度,控制箱接收传感器的监测数据并反馈至上位机,上位机实时计算工作面向采空区的实际漏风量,控制箱中的伺服电机控制风窗面积进行调控风量,上位机根据调控后的漏风量做出进一步指令,直到达到设定的漏风量要求,具体步骤如下:
步骤1.在采空区漏风区域工作面两侧的进风巷与回风巷安装超声波风速传感器,进风巷安装智能风门-风窗,在智能风门-风窗附近安装配套控制箱;
步骤2.上位机根据工作面两侧进风巷与回风巷风量实时计算实际漏风量Q,根据工况要求在上位机中设定工作面漏风量Q,上位机实时计算实际漏风量Q与设定工作面漏风量Q的差值e,以及差值变化率ec,并将差值e和差值变化率ec传输到控制箱中;
步骤3.将上位机计算所得的实际漏风量Q与设定工作面漏风量Q的差值e和差值变化率ec作为输入变量传输给控制箱中的变论域模糊PID调控系统;经过变论域模块得出伸缩因子,伸缩因子对初始比例因子和初始量化因子进行调控;经过模糊调节模块得出PID的三个参数KI、KP、KD;将漏风量差值e、差值变化率ec、KI、KP、KD作为PID控制器模块的输入变量,PID控制器模块输出伺服电机的控制信号;
步骤4.伺服电机根据接收的控制信号调节风门-风窗过风面积,改变进风量,平衡所监测的采空区与周边采空区之间的压差,直到实际漏风量Q低于或者相等设定漏风量Q;若上位机计算得出漏风量未达到所设定的数值,将重新进入步骤3中开始计算,直到实际漏风量低于或等于设定漏风量。
2.根据权利要求1所述的多风机联合运转干扰下采空区漏风的智能调控方法,其特征在于:所述控制箱为KXJ127矿用隔爆兼本安型PLC控制箱,在PLC控制箱中设置变论域模糊PID调控系统,控制箱的信号接收端连接用于检测风量的超声波风速传感器、采煤工作面回风隅角甲烷传感器T0和两个采煤工作面回风巷甲烷传感器;控制箱的信号输出端连接伺服电机的控制端;控制箱与上位机、伺服电机信号双向传输。
3.根据权利要求2所述的多风机联合运转干扰下采空区漏风的智能调控方法,其特征在于:所述控制箱的变论域模糊PID调控系统包括变论域模块、模糊调节模块和PID控制器模块。
4.根据权利要求1或2所述的多风机联合运转干扰下采空区漏风的智能调控方法,其特征在于:所述超声波风速传感器设置在采空区漏风工作面的进风巷和回风巷,距离工作面切眼60m开始向采空区方向等间距安装一组用于监测所漏风工作面进风巷和回风巷风量的超声波风速传感器,且进风巷和回风巷两侧的超声波风速传感器的数量均不少于3组。
5.根据权利要求1所述的多风机联合运转干扰下采空区漏风的智能调控方法,其特征在于:所述步骤2中,上位机实时计算实际漏风量Q与设定工作面漏风量Q的差值e与差值变化率ec,计算方法如下:
步骤2.1设进风巷超声波风速传感器监测的风量为Q进i,表示第i个传感器监测的风量,同理回风巷超声波传感器监测的风量为Q回i,进风巷风量为
Figure FDA0004158046770000021
其中n表示传感器的数量,同理计算/>
Figure FDA0004158046770000022
工作面实际漏风量Q=Q-Q
步骤2.2实际漏风量与设定漏风量差值e=Q-Q,差值变化率
Figure FDA0004158046770000023
6.根据权利要求1所述的多风机联合运转干扰下采空区漏风的智能调控方法,其特征在于:所述上位机中设定工作面最低需风量,将回风巷风量视为工作面风量,若工作面风量低于工作面最低需风量需重新设定漏风量Q
7.根据权利要求1所述的多风机联合运转干扰下采空区漏风的智能调控方法,其特征在于:变论域模糊PID调控系统具体包括以下内容:
S1,设输入变量差值e与差值变化率ec初始论域分别为[-e0,e0]、[-ec0,ec0],输出变量ΔKP、ΔKI、ΔKD的初始论域分别为[-Kp0,Kp0]、[-Ki0,Ki0]、[-Kd0,Kd0],并设输入变量差值e与差值变化率ec以及输出变量ΔKP、ΔKI、ΔKD的模糊论域均为[-M,M];
S2,将实际漏风量Q与设定工作面漏风量Q的差值e与差值变化率ec作为变论域模块的输入变量,采用伸缩因子函数模型得出输入变量差值e的伸缩因子αe、差值变化率ec的伸缩因子αec和三个输出变量ΔKP、ΔKI、ΔKD的伸缩因子βP、βI、βD;其中伸缩因子函数模型如下:
输入伸缩因子:
Figure FDA0004158046770000024
输出伸缩因子:
Figure FDA0004158046770000031
式中λ∈(0,1),k>0;KI为比例因子;Pi为输入权重系数;ei(τ)为输入变量误差;β(0)为输出论域伸缩因子初值;
S3,根据S1中所设定的初始论域和模糊论域得出初始量化因子Ke0、Kec0以及初始比例因子LP0、LI0、LD0,即
Figure FDA0004158046770000032
Figure FDA0004158046770000033
将S2中所得的伸缩因子αe、αec对初始量化因子Ke0、Kec0进行调整得出量化因子Ke与Kec,即Ke=αe·Ke0、Kec=αec·Kec0;差值e与差值变化率ec作为模糊调节模块输入变量且经过量化因子作用模糊化处理得到模糊量e*与ec *,即e*=Ke·e、ec *=Kec·ec;参照专家领域知识设置模糊集和模糊规则以及采用重心法对模糊量e*与ec *进行模糊处理得到输出值ΔKP *、ΔKI *、ΔKD *;利用所输出变量的伸缩因子βP、βI、βD对初始比例因子进行调节得到比例因子LP、LI、LD,即LP=βP·LP0、LI=βI·LI0、LD=βD·LD0;将输出值ΔKP *、ΔKI*、ΔKD*对三个比例因子对应进行比例调节得到参数调整量ΔKP、ΔKI、ΔKD,即ΔKP=ΔKP *·LP、ΔKI=ΔKI *·LI、ΔKD=ΔKD *·LD
S4,通过PID控制器参数整定方法得出初始值KP0、KI0、KD0,根据上述得出的调整量ΔKP、ΔKI、ΔKD对初始值进行调整,获得PID控制的3个参数KP=KP0+ΔKP、KI=KI0+ΔKI、KD=KD0+ΔKD;将差值e与差值变化率ec以及3个参数KP、KI、KD作为PID控制器模块的输入变量,PID控制器模块输出控制信号传输到伺服电机;其中关于PID控制算法的离散化公式为:
Figure FDA0004158046770000034
式中un为调节n次后输出的结果,en为第n次实际漏风量Q与设定漏风量Q差值,Ts为调控时间,KP为比例常数、TI积分时间常数、TD微分时间常数。
8.根据权利要求1所述的多风机联合运转干扰下采空区漏风的智能调控方法,其特征在于:所述上位机中设定工作面最低需风量,若调节后工作面风量低于最低需风量,重新设定漏风量Q;上位机同时实时监测回风隅角T0与回风巷的瓦斯浓度数据,当瓦斯浓度数据超限时上位机开启气体浓度警告,发出警告指令到控制箱进行语音警告,并控制箱输出控制信号到伺服电机,控制智能风门-风窗开度恢复到最大状态,待瓦斯浓度正常后,根据实时监测的风量,再重新设定漏风量Q并保证瓦斯浓度不再超限,继续进行调节。
CN202310340916.0A 2023-04-03 2023-04-03 一种多风机联合运转干扰下采空区漏风的智能调控方法 Pending CN116357374A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310340916.0A CN116357374A (zh) 2023-04-03 2023-04-03 一种多风机联合运转干扰下采空区漏风的智能调控方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310340916.0A CN116357374A (zh) 2023-04-03 2023-04-03 一种多风机联合运转干扰下采空区漏风的智能调控方法

Publications (1)

Publication Number Publication Date
CN116357374A true CN116357374A (zh) 2023-06-30

Family

ID=86906881

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310340916.0A Pending CN116357374A (zh) 2023-04-03 2023-04-03 一种多风机联合运转干扰下采空区漏风的智能调控方法

Country Status (1)

Country Link
CN (1) CN116357374A (zh)

Similar Documents

Publication Publication Date Title
US11466568B2 (en) Optimized mine ventilation system
CN109252884B (zh) 一种沿空留巷工作面y型通风调压防灾方法
CN102797451B (zh) 一种深水表层钻井井下环空压力自动控制系统及控制方法
CN104329108B (zh) 一种矿井多工作面风压整体动态平衡控制系统的控制方法
CN114738031B (zh) 自然和机械通风耦合的施工隧道数字孪生智能通风方法
CN202850970U (zh) 一种深水表层钻井井下环空压力自动控制系统
CN101210495A (zh) 全断面大型掘进装备的动态协调控制系统
CN100507207C (zh) 通过井生产油气混合物的方法和系统
CN103306701B (zh) 大区域自动控制均压防灭火系统
CN112228136A (zh) 一种快速排出密封巷道有毒有害气体系统及使用方法
CN115013026A (zh) 一种用于隧道火灾烟气控制的智能型重点排烟方法和系统
CN112177657A (zh) 长距离隧道用分段式通风系统及其实施方法
CN116357374A (zh) 一种多风机联合运转干扰下采空区漏风的智能调控方法
CN109578048B (zh) 矿井智能通风系统与控制方法
CN111594252A (zh) 一种极易自燃煤层开采期间采空区煤炭自燃的防治方法
CN213360154U (zh) 长距离隧道用分布式通风系统
CN111119968B (zh) 基于耗能降速的煤矿顶板大面积来压暴风危害消除方法
CN110374655B (zh) 一种煤矿掘进工作面局部通风智能控制方法及控制系统
CN112685926A (zh) 一种浅埋煤层防治漏风工程的施工周期判别方法
CN113202486A (zh) 一种富水岩层samson系统辅助保压推进方法
CN113513361B (zh) 一种高瓦斯煤矿回采工作面瓦斯浓度调控方法及系统
CN112539065B (zh) 一种顶管开挖面稳定性的监控系统及方法
CN113719295B (zh) 一种泥水盾构法隧道掘进全过程稳定性智能控制系统
CN114837724A (zh) 金属矿山深井开采均压通风调控装置及调控方法
Wei et al. Control Mode of Mine Ventilation System and Its Implementation Based on Monitoring Data and Safety Information

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