矿井通风系统测定巷道中风流速度时多采用的方法为走线法及点测法，走线法主要使用翼式风表在测风巷道沿既定路线测定单位时间内巷道断面的平均风速[1]；点测法则使用热敏式风表对巷道内多个特定位置风速进行测定，并以多个测点风速的平均值作为实际风速测量值，而巷道风流运动时各位置风速不完全一样，部分线路及特定点不能代表实际平均风速[2-3]。因此，走线法和点测法在实际测量时均会出现较大误差[4-5]。
为减少实际测量误差，此次研究拟对实际巷道断面风流分布进行研究，运用数学模型表示其分布规律，并在此基础上实现对巷道实际风量进行更精确的计算。
1 风速试验测试
1.1 测定系统
试验系统采用如图1所示的结构，装置截面为圆形[6-7]，装置动力由具有变频控制风速装置的风机提供，测试系统的风流在变频风机作用下，经装置右侧的整流格栅进入系统。为方便试验测试，装置中部设有高速和低速测风口，以便试验测定巷道断面各处风速值。
图1 测试系统结构Fig.1 The structure of test system
1.2 试验流程
测试利用装置低速测风口，测定该处断面各位置风速，并利用装置的位置标尺记录各测点与巷道中心轴线距离，从而记录巷道各不同位置处的实际风速值，风速测量装置选用宜于操作的热敏式风表。同时，为测试不同风速下巷道不同位置处的风速变化，试验测试时通过变频器设置风机档位为低、中、高三档，重复测试步骤测定低、中、高档位各点风速测定。
2 研究测试
2.1 测试分析
以测试断面中轴线为参考位置，测定距离测试断面中轴线不同位置处的风速值并记录，测试圆形断面直径为0.5 m，沿中轴线对称测定各点风速值，测点间隔小于5 cm。不同风速条件下(高速、中速、低速)重复试验测试，测试数据整理结果见表1，根据整理数据以风速值为纵坐标，距测定断面中轴线距离为横坐标，绘制表1对应的风速变化曲线图见图2。
表1 测试数据整理Table 1 Test data arrangement与中轴线距离/cm实测风速值/(m·s-1)低档中档高档..512.143..563.524..633.675..693.845..783.905..
图2 风速变化曲线Fig.2 Wind speed curves
2.2 测试结果分析
从上述试验结果并结合图2风速变化曲线，可知各风速条件下风速分布规律基本一致，即风流分布呈现距离巷道中轴线越近，风速越高的整体趋势[8]。结合图2所示这种分布规律区别于简单的线性关系，更类似于指数形式的数学关系，分析其原因是由于巷道壁并非绝对光滑，存在一定粗糙度[9-10]，靠近边壁的风流受巷道壁的黏滞，这种黏滞作用越靠近边壁越明显，实际测定时应考虑巷道壁对风速的黏滞作用，其作用范围与巷道实际的粗糙度有关[11]。
2.3 风流变化规律数学模型表示
根据2.1及2.2节的试验结果及分析，指出巷道风流分布规律大体趋势为类指数形式。因此，为将风流分布规律利用数学模型表征，研究拟以实际测定风速值为目标函数值V(x)，对应各实测点的位置(距巷道壁距离)为自变量x，利用以自然常数e为底的指数函数形式，对其进行非线性回归分析，得出相对应的回归方程，以建立巷道断面风流变化的数学模型。鉴于实际风流变化是以巷道中轴线对称分布，为简化回归数学模型，研究仅以巷道轴线对称的一侧测定数据(自变量取值范围0～25 cm)做回归分析。采用以e为底的指数函数形式，对不同风速(高、中、低)下巷道风流变化回归分析，其回归方程表示具体如下：
1)低风速回归方程
V(x)=-2.99e(-x/1.96)+2.89
2)中风速回归方程
V(x)=-4.23e(-x/2.01)+4.09
3)高风速回归方程
V(x)=-5.89e(-x/1.96)+5.73
绘制各回归方程曲线如图3所示。
图3 三种风速条件下回归方程曲线Fig.3 Regression equation curves under three wind speed conditions
各回归方程的相关度与均方根误差值如表2所示。
表2 回归方程可信度分析Table 2 Reliability analysis of regression equation可信度指标风速条件/(m·s-1)低档中档高档相关度均方根误差
回归相关度和均方根误差作为回归可信度分析的重要指标，从表2各指标数据分析，该回归方程可信度较高，可以作为风流变化规律的数学表示模型。
2.4 巷道风量计算
根据式(1)～(3)，将所测同一巷道断面不同位置风速的离散点，用数学模型拟合出风速与巷道壁距离间关系连续变化的曲线。由于研究的模拟巷道为圆形巷道，具有统一对称性，即上述研究所揭示的风流分布规律在圆形巷道各个方向上均一致。因此，在距离巷道壁同样距离处风速相同，整个巷道的风速分布变化应是形成距巷道壁距离相关的等速圈的形式，如图4所示。根据上述研究数学关系方程做轴线，某点风速V与其在轴向增量△x所示圆环面积的乘积即为该环内的风量，由此对整个圆形区域风量做积分，得出圆形巷道实际精确风量。根据图4有如下巷道断面内风速与风量间数学关系：
Q=2π(25-x)×V(x)dx
式中：Q为巷道风量值；V(x)为风速与距巷道壁距离x的关系函数。
将(1)～(3)代入(4)式得出如下关系式：
低速条件计算风量：
Q=2π(25-x)×[-2.99exp(-x/1.96)+2.89]dx
中速条件计算风量：
Q=2π(25-x)×[-4.23exp(-x/2.01)+4.09]dx
高速条件计算风量：
Q=2π(25-x)×[-5.89exp(-x/1.96)+5.73]dx
图4 巷道风流分布变化示意图Fig.4 Sketch map of airflow distribution change in roadway
分别对低、中、高速条件下所示积分方程进行计算，并将计算结果单位转换为国际单位，得到结果见表3。
表3 巷道风量计算结果Table 3 Calculation results of roadway air volume风速档位计算风量/(m3·s-1)低速0.484 3中速0.679 9高速0.957 1
2.5 巷道风量计算值与常规测定值对比
若常规巷道测定风量采用分格定点法，忽略靠近巷道壁的低风速区域，依排除低速区其他风速平均值作为巷道的平均风速，则此时低、中、高速档位对应的测定风量及其与巷道实际风量计算值间的数值对比如表4所示。
由表4可知，实际巷道通风测定若不考虑巷道边界区域风速低速区，则实际测定风量将高于计算风量值，约高出15%。巷道边界低速区的存在是由于风流受到巷道壁的黏滞作用，这种黏滞作用与巷道的粗糙度有关系，实际风速测时应考虑此黏滞作用对风速、风量的影响。有关巷道粗糙度造成巷道风流黏滞与巷道边界低速区域大小的关系需做进一步研究。
表4 巷道风量测定值与计算值对比Table 4 Calculated and measured results of roadway air volume风速档位测定风量/(m3·s-1)计算风量/(m3·s-1)相差百分比/%低速0.553 80.484 314.35中速0.783 00.679 915.17高速1.099 00.957 114.83
3 结语
通过对模拟巷道风流分布规律进行研究，以巷道风流分布规律为基础采用曲线回归拟合分析与积分计算相结合的方式，对巷道实际风量精确计算，此次试验研究结果总结分列如下：
1)三种风速条件下(高、中、低档)，巷道风流分布变化趋势大体一致，风流在巷道中的变化趋势与风速无关。
2)以风速为目标值V(x)，以距巷道壁距离为变量X，对三种风速条件下试验数据进行曲线拟合分析，得出相应的回归方程。根据回归方程可信度分析，回归方程可靠性极高，可作为该类型巷道风流变化规律数学表征模型，并可作为巷道风速值连续变化计算依据。
3)利用回归方程V(x)与距巷道壁x间的连续关系，构造相应的积分方程计算出巷道的实际风量值，结果显示风量计算值与实际测定值之间存在约15%的误差，实际通风测定时应考虑该误差。
4)由于巷道壁粗糙度而形成对风流的黏滞作用，使风流在靠近边壁处，风流速度较小，越靠近巷道中心区域，风流受此黏滞作用越不明显，风速也越大。有关巷道粗糙度造成巷道风流黏滞与巷道边界低速区域大小的关系需要做进一步研究。
[1] 李亚俊，李印洪，姚银佩，等.圆形巷道风流分布规律的数学模型表征[J].湖南有色金属，2020，36(3)：1-3.
LI Yajun，LI Yinhong，YAO Yinpei，et al.The study of mathematical model characterization of air flow distribution in roadway[J].Hunan Nonferrous metals，2020，36(3)：1-3.
[2] 李亚俊，李印洪，吴洁葵，等.巷道断面风流分布规律试验研究[J].有色金属(矿山部分)，2019，71(5)：102-104，110.
LI Yajun，LI Yinhong，WU Jiekui，et study on air flow distribution law of roadway[J].Nonferrous metals(Mining Section)，2019，71(5)：102-104，110.
[3] 李亚俊，李印洪，吴洁葵，等.热敏式风表风速测量参数优化[J].有色金属(矿山部分)，2019，71(4)：111-113.
LI Yajun，LI Yinhong，WU Jiekui，et al.The wind speed measurement optimization of thermosensitive anemometer[J].Nonferrous metals(Mining Section)，2019，71(4)：111-113.
[4] 杨宇.不同形状巷道断面风流—瓦斯耦合特性研究[D].太原：太原理工大学，2017.
YANG Yu.Study on the coupling characteristics of airflow and gas in tunnel section with different shapes[D].Taiyuan：Taiyuan University of Technology，2017.
[5] 罗永豪.巷道断面风速分布与煤矿通风系统实时诊断理论研究[D].太原：太原理工大学，2015.
LUO study on wind velocity distribution in the section of roadway and real time[D].Taiyuan：Taiyuan University of Technology，2015.
[6] 赵丹，黄福军，陈帅，等.点风速与平均风速在圆形巷道中的关系[J].辽宁工程技术大学学报(自然科学版)，2014，33(12)：1654-1659.
ZHAO Dan，HUANG Fujun，CHEN Shuai，et between point air velocity and average air velocity in circular roadway[J].Journal of Liaoning Technical University(Natural Science)，2014，33(12)：1654-1659.
[7] 周西华，孟乐，李诚玉，等.圆形管道风速测定与校正方法实验[J]验[J].辽宁工程技术大学学报(自然科学)，2012，31(6)：801-804.
ZHOU Xihua，MENG Le，LI Chengyu，et study on determination and correction method of wind speed in circular pipe[J].Journal of Liaoning Technical University(Natural Science)，2012，31(6)：801-804.
[8] 丁翠.矿井巷道风流状态“关键环”实验与数值研究[D].北京：中国矿业大学(北京)，2013.
DING and numerical study of “key ring” of airflow status in the mine tunnel[D].Beijing：China University of Mining and Technology(Beijing)，2013.
[9] 葛启发，朱维根，朱瑞军，等.深井难采矿体开采与通风控制技术研究[J].有色金属工程，2017，7(6)：86-89，99.
GE Qifa，ZHU Weigen，ZHU Ruijun，et of mining method and ventilation control technology for difficult-to-mining orebodies in deep mine[J].Nonferrous metals Engineering，2017，7(6)：86-89，99.
[10] 周志杨，王海宁，晏江波，等.改扩建矿井通风系统优化及动态管理模型创建[J].有色金属工程，2017，7(5)：75-79.
ZHOU Zhiyang，WANG Haining，YAN Jiangbo，et system optimization of mining reconstruction and extension and construction Its dynamic management model[J].Nonferrous metals Engineering，2017，7(5)：75-79.
[11] 支学艺，刘建华，赖春明.风流经转弯后稳定段长度数值模拟及应用研究[J].有色金属科学与工程，2016，7(1)：96-101.
ZHI Xueyi，LIU Jianhua，LAI simulation and application of airflow’s stable length after turning[J].Nonferrous metals Science and Engineering，2016，7(1)：96-101.
文章来源：风流一代 网址: http://fenglyd.400nongye.com/lunwen/itemid-51178.shtml
上一篇： 暂无
下一篇： 林业论文_咬定青山不放松
