发明名称 一种校核矿井巷道摩擦风阻值的方法
摘要 本发明属于采矿工程矿井通风领域,具体涉及一种校核矿井巷道摩擦风阻值的方法。本发明技术方案是首先测试得到矿井巷道的摩擦风阻,将其录入到本发明的矿井通风智能管理系统中,通过矩阵演算得到矿井通风网络中的关键通路以及关键通路中的高信任分支,以这一分支的摩擦风阻值为基准,反过来校核其他巷道的摩擦风阻,进而得到整个矿井通风网络各巷道最接近真值的摩擦风阻,然后在系统中验证性工业试验,以确保本发明技术方案的准确性。本发明方发依托PC计算机,实现了矿井通风系统在计算机屏幕上的反演与模拟,与矿井现场测试值最大误差小于4%,为矿井通风管理与技术改造提供了科学依据,为煤炭矿井的安全生产提供了有力保障。
申请公布号 CN102323998B 申请公布日期 2014.02.05
申请号 CN201110299643.7 申请日期 2011.09.29
申请人 辽宁工程技术大学 发明人 贾廷贵;曲国娜;贾宝山
分类号 G06F19/00(2011.01)I;E21F1/02(2006.01)I 主分类号 G06F19/00(2011.01)I
代理机构 沈阳东大专利代理有限公司 21109 代理人 李在川
主权项 1.一种校核矿井巷道摩擦风阻值的方法,其特征在于按照以下步骤进行:首先测定矿井各通风巷道的通风阻力h和风量Q,并将实际的矿井各通风巷道网络转化为矿井通风网络数据结构录入计算机中,根据通风阻力定律式(1),求得R<sub>j</sub>;<maths num="0001"><![CDATA[<math><mrow><msub><mi>h</mi><mi>j</mi></msub><mo>=</mo><msub><mi>R</mi><mi>j</mi></msub><mo>&CenterDot;</mo><msubsup><mi>Q</mi><mi>j</mi><mn>2</mn></msubsup><mo>,</mo><mrow><mo>(</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中:h<sub>j</sub>代表矿井通风系统中的第j巷道的通风阻力,单位Pa;R<sub>j</sub>代表矿井通风系统中的第j巷道的摩擦风阻,单位Ns<sup>2</sup>/m<sup>8</sup>;Q<sub>j</sub>代表矿井通风系统中的第j巷道的风量,单位m<sup>3</sup>/s;将R<sub>j</sub>录入矿井通风智能管理系统中,完成数字化系统数据的初始赋值;然后,采用深度优先搜索法遍历矿井通风网络,得到矿井通风网络独立通路矩阵P,枚举并对比各通路阻力H<sub>i</sub>与矿井通风总阻力H<sub>t</sub>的相对误差ε<sub>i</sub>;P=(p<sub>ij</sub>)<sub>s×n</sub>,i=1,2,…,s;j=1,2,…,n    (2)<maths num="0002"><![CDATA[<math><mrow><msub><mi>p</mi><mi>ij</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn><mo>,</mo><mrow><mo>(</mo><msub><mi>e</mi><mi>j</mi></msub><mo>&Element;</mo><msub><mi>P</mi><mi>i</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mn>0</mn><mo>,</mo><mrow><mo>(</mo><msub><mi>e</mi><mi>j</mi></msub><mo>&NotElement;</mo><msub><mi>P</mi><mi>i</mi></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msub><mi>H</mi><mi>i</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>p</mi><mi>ij</mi></msub><mo>&CenterDot;</mo><msub><mi>h</mi><mi>j</mi></msub><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mi>s</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>ε<sub>i</sub>=|H<sub>i</sub>-H<sub>t</sub>|/H<sub>t</sub>×100%    (5)其中:s代表矿井通风系统中构成独立通路的个数;n代表矿井通风系统中巷道个数;P<sub>ij</sub>代表当矿井通风系统中的第j巷道在第i个通路上时取1值,当矿井通风系统中的第j巷道不在第i个通路上时取0值;H<sub>i</sub>代表第i条通路上的通风阻力,单位Pa;H<sub>t</sub>代表矿井通风系统中矿井通风总阻力,单位Pa;ε<sub>i</sub>代表通路通风阻力与矿井通风系统总阻力的相对误差,单位%;选择ε<sub>i</sub>≤8%的通路作为关键通路集,在关键通路集中选择各巷道的测定通风阻力h<sub>j</sub>为基础,设定W<sub>j</sub>为第j条边或第j条分支在其所在回路中的逼近系数,根据回路压力平衡定律式(6),有:<maths num="0004"><![CDATA[<math><mrow><msub><mi>W</mi><mi>j</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>l</mi></munderover><msub><mi>c</mi><mi>ij</mi></msub><msub><mi>h</mi><mi>j</mi></msub><mo>,</mo><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths><img file="FDA0000363014490000022.GIF" wi="1356" he="224" />其中:C为独立回路矩阵;c<sub>ij</sub>为回路矩阵中的元素,即当矿井通风系统中的第j巷道在第i个回路上且风流方向与回路方向一致时取1值,即当矿井通风系统中的第j巷道在第i个回路上且风流方向与回路方向相反时取-1值,当矿井通风系统中的第j巷道不在第i个回路上时取0值;选取出关键通路中各巷道的逼近系数W<sub>j</sub>,使其组成最佳逼近系数集,在数值集中设k个数据筛选区间,每个区间的域由式(8)确定,W<sub>j</sub>在域X<sub>t</sub>的频率由式(9)得到,<maths num="0005"><![CDATA[<math><mrow><msub><mi>X</mi><mi>t</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>X</mi><mn>0</mn></msub><mo>=</mo><msub><mi>W</mi><mn>1</mn></msub><mo>+</mo><mfrac><mn>1</mn><mi>k</mi></mfrac><mrow><mo>(</mo><msub><mi>W</mi><mi>max</mi></msub><mo>-</mo><msub><mi>W</mi><mi>min</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>X</mi><mi>t</mi></msub><mo>=</mo><msub><mi>X</mi><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><mfrac><mn>1</mn><mi>k</mi></mfrac><mrow><mo>(</mo><msub><mi>W</mi><mi>max</mi></msub><mo>-</mo><msub><mi>W</mi><mi>min</mi></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>,</mo><mrow><mo>(</mo><mi>t</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><msub><mi>b</mi><mi>j</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mi>a</mi><mo>=</mo><mi>a</mi><mo>+</mo><mn>1</mn><mo>,</mo></mtd><mtd><msub><mi>W</mi><mi>j</mi></msub><mo>&Element;</mo><mrow><mo>(</mo><msub><mi>X</mi><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>X</mi><mi>t</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>a</mi><mo>=</mo><mi>a</mi></mtd><mtd><msub><mi>W</mi><mi>j</mi></msub><mo>&NotElement;</mo><mrow><mo>(</mo><msub><mi>X</mi><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>X</mi><mi>t</mi></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>,</mo><mrow><mo>(</mo><mi>t</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo><mrow><mo>(</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mi>l</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中:W<sub>max</sub>代表关键通路中各巷道的最逼近数值;W<sub>min</sub>代表关键通路中各巷道的最远离数值;X<sub>t</sub>代表数据区间函数;l代表独立回路个数,即矿井通风系统中构成独立通路的个数;k代表数值滤波的数据区间个数;b<sub>j</sub>代表W<sub>j</sub>在X<sub>t</sub>区间的频率函数;a代表统计W<sub>j</sub>在域X<sub>t</sub>出现的频数;根据频率函数b<sub>j</sub>,由高频率确定通路集中关键巷道的原则,获取最高信任测试巷道摩擦风阻R<sub>双-i</sub>,从而确定了矿井通风网络的关键通路及关键通路中高信任分支,以此为测试基准,将R<sub>双-i</sub>代入公式(1)中,得到本方法测得的巷道通风阻力h<sub>双-i</sub>;以h<sub>双-i</sub>为测试基准,将其代入公式(1)中,计算得到其它巷道的摩擦风阻R;对矿井通风系统的每一条巷道重复进行上述测定和计算,对于每一条巷道共得到i个摩擦风阻值,选取其中出现频率最高的风阻值,确定其为整个矿井通风网络各巷道最接近真实值的摩擦风阻值;将校核后的各巷道的摩擦风阻值重新录入矿井通风智能管理系统数据库,生成校核后的矿井通风智能管理系统;在校核后的矿井通风智能管理系统中对矿井通风系统进行验证性工业试验,打开风门,改变调窗面积,与矿井实际工业试验得到的摩擦风阻真实值进行比较,验证其校核处理的精度与准度,在工业试验验证符合实际情况的前提下,便可使用该系统进行矿井通风的设计与管理,当发生与实际情况不符的情况,则返回上述步骤重新进行校核。
地址 110043 辽宁省阜新市中华路47号
您可能感兴趣的专利