矿井通风阻力测定数据平差处理方法及应用

1、第15卷第10期2006年10月中国矿业CHINA MINING MAG AZINEOctober 2006矿井通风阻力测定数据平差处理方法及应用陈宙1赵恩平2蒋仲安1林建广1(11北京科技大学土木与环境工程学院教育部金属矿山高效开采与安全重点实验室北京100083;21邯邢矿山管理局北铭河铁矿邯郸056300摘要:本文论述了矿井通风阻力测定数据平差处理的两种方法:条件平差和间接平差,并基于通风网络和平差理论编制了相应的通风阻力测量数据平差软件,最后以西石门铁矿中区通风阻力实测数据为例,验证了该种分析方法的可靠性和及实用性。关键词:通风阻力测量条件平差间接平差中图分类号:TD722文献标识码:

2、B 文章编号:1004-4051(200610-0105-04DATA ADJUSTMENT METH OD OF RESISTANCE MEASUREMENT FOR MINE VENTI LATIONAN D ITS APPL ICATIONChen Zhou 1Zhao Enping 2Jiang Zhongan 1Lin Jianguang 1(11Civil and Environment Engineering School ,USTB Beijing 100083;21Beiminghe iron ore ,Hanxing Metallurgical Mine administr

3、ation Handan 056300Abstract :Two kinds of mine ventilation resistance measurement data adjustment method :ConditionAdjustment and Indirect Adjustment ,are discussed in the article.The application about mine ventilation re 2sistance measurement data adjustment is developed on the basis of ventilation

4、 network and adjustment theo 2ry.Then ,taking Xishimen Iron ore ventilation resistance measurement data as an example ,the article veri 2fies the reliability and practicability of these analytical methods.K ey w ords :Mine ventilation Resistance Measurement Condition adjustment Indirect adjustment收稿

5、日期:2006-04-17作者简介:陈宙(1982-男硕士研究方向:安全技术及工程1概述矿井通风阻力测定是矿井通风技术管理工作的重要内容之一,测量数据的好坏主要是由测量数据的误差大小来判断。尽管目前有较多的通风阻力测定方法,但无论采用哪种方法都会产生一定的误差,同时受观测条件和环境、观测者、测量仪器的影响,其测量结果不可避免地存在着或多或少误差。因此,如何处理带有误差的观测值,找出待求量的最佳估值,是测量平差所研究的内容。测量平差即是测量数据调整的意思,其基本的定义是:依据某种最优化准则,把一系列带有观测误差的测量数据,求定未知量的最佳估值及精度的理论和方法。在测绘工程和其它工程领域中,只带有

6、偶然误差的观测占有大多数,是比较普遍的情形,这是测量平差研究的基础。2数据平差处理方法211条件平差一般测量时为了检查错误和提高最后成果的精度,常常做多余观测。如按理论,在一回路中,如果想得到所有巷道的风压,只需知道其中任意n 21条分支的风压就可以了。观测这样的n 21条分支风压,称为必要观测。但是通常观测n 条分支的风压,观测的第n 条分支风压为多余观测。由于观测值带有误差,致使闭合回路风压之和不等于0,产生了不符值。如果有r 个多余观测,就产生r 个平差值条件方程。采用平差值条件方程来达到消除不符值,同时运用求条件极值的原理来解出改正数,从而求得平差值,这样以条件方程为函数模型的平差方法

7、,就是条件平差。确立方程等式个数。条件方程的个数与多余观测的个数相同,即r =n 2t ;其中n 为观测的数据数,t 为必要观测数。中国矿业第15卷列立方程。确定了方程的个数后,列出的方程之间要求线性无关。对于一个网络回路来说,其包含的分支可能又被另外的回路所包含,所列出的条件方程多种多样。但是所有列出的方程中必然只有r个线性无关平差值方程。一般一个独立回路对应一个线性无关的方程。条件平差一般方程如下所示:a11L31+a12L32+a1n L3n+c1=0a21L31+a22L32+a2n L3n+c2=0a r1L31+a r2L32+a rn L3n+c r=0(1式中a ij(i=1,

8、2,r;j=1,2,n为条件方程系数;c ij(i=1,2,r为条件方程常数项;L3未知平差值。以L3i=L i+v i(i=1,2n(i=1,2n代入(1式中得条件方程:a11v1+a12v2+a1n v n+w1=0a21v1+a22v2+a2n v n+w2=0a r1v1+a r2v2+a rn v n+w r=0(2式中v i(i=1,2,n为实测值与平差值的差;w i为条件方程闭合差,或称为条件方程的不符值;即:a11L1+a12L2+a1n L n+c1=w1a21L1+a22L2+a2n L n+c2=w2a r1L1+a r2L2+a rn L n+c2=w r(3(1(2和

9、(3式矩阵表达式为:AL3+C=O(4AV+W=O(5AL+C=W(6按照最小二乘法准则,要求V T PV取极小值,因此通过条件极值的拉格朗日乘数法,设其乘数为K=k1,k2k rT,K称为联系数向量。组成新函数:=V T PV-2K T(AV+W将对V求一阶导,并令其为0,可得:V=P-1A T K(7代入(5得:A P-1A T K+W=O(8联立解方程(7(8求出V值,代入L3i= L i+v i便可获得平差值L3i 212间接平差间接平差是通过选定足以确定某个平差问题的t个未知数来消除观测值之间的不符值,并且求自由极值的方法解出未知参数的最或然值,从而求得各个观测值的平差值。间接平差方

10、程中未知数的个数等于必要观测的个数t,方程的个数为观测值的个数。其一般方程为:a11x1+a12x2+a1t x t+c1=L1+v1a21x1+a22x2+a2t x t+c2=L2+v2a n1x1+a n2x2+a nt x t+c n=L n+v n(9令l i=c i-L i则得误差方程的矩阵表达式为:V=A X+l(10同上要求V T PV取极小值,按照数学上自由极值理论,即:V T PVx=VT PA=O当L是独立观测值时,P是对角阵。转置得A T PV=O(11联立式(10,(11可求得X,代入(10式求出改正向量V。由此得观测值的平差值L3i=L i +v i。3通风阻力测定

11、数据平差及软件的编制311通风阻力测定数据平差通风阻力测量数据的平差分为风量平差和阻力平差。通过平差后的阻力和风量从而求出巷道的校正风阻;下面以分支数n,节点数m的风网为例,分别讲述风量平差和阻力平差。(1风量平差若采用条件平差,根据节点风量平衡原理,只需测量n-m+1条分支的风量(分支必须为余支,就可以求出其它分支风量。其余多于余支数目的测量为多余观测。一般我们可以认为多余观测数越多,平差值越接近真值。因此,假设我们测得全部分支风量为n个,则多余观测数即条件方程的个数为m-1个。此时条件平差方程(4的系数矩阵A(M×n为节点独立关联矩阵(M=m-1,W(M×1为节点风量不

12、符值矩阵,P(n×n为权对角线矩阵。若采用间接平差,同样必须测量n-m+1条分支风量,即风量间接平差方程数最少应大于n-m+1。假设测得全部分支风量n个,则系数矩阵A(M×n为节点独立关联矩阵,P(n×n为权对角线矩阵,向量X的维数为n-m+1。601第10期陈宙等:矿井通风阻力测定数据平差处理方法及应用(2阻力平差若采用条件平差,根据回路风压平衡理论,在每一个独立的回路中,如果回路包含K 条分支,只要测得回路的巷道数k -1,则可以计算出其余一条巷道的风压。那么在整个风网中只需测得m -1条分支风压,而计算出其余n -m +1条分支风压。同样假设我们测得所有分支

13、阻力,则条件平差方程(4的系数矩阵A (M ×n 为独立回路矩阵(M =n -m +1,W (M ×1为回路风压不符值矩阵,P (n ×n 为权对角线矩阵。若采用间接平差,则风压间接平差方程数最少应大于m -1。如果测得全部分支风压n 个,则系数矩阵A (M ×n 为独立回路矩阵,P (n ×n 为权对角线矩阵,向量X 的维数为m -1。312程序的编制根据条件平差和间接平差方程及回路矩阵构造算法,本文采用条件平差处理方法模型编制了相应的计算程序V EN T O P T ,本程序采用基于win 2dows 操作系统的面向对象语言delp hi

14、语言编写,并在PC 机上成功运行。程序结构框图如图1所示: 图1程序结构框图该计算程序适用于各种测量方法的结果平差处理。程序运行初始需要输入网络的拓扑值、测得的分支风量和阻力值。程序运行后能够计算出平差后的风量和阻力值,并得出相应的风阻。程序提供打印数据功能,并可导出数据到excel 中保存。4应用实例 邯邢矿山局西石门铁矿中区通风系统由1#,2#副井,402120m 水平回风井,中区回风斜井组成。现在402120m 水平回风井安装K45262No19型200kW 离心式风机1台,1002120m 斜坡道废弃油库中安装K40242No15风机110kW1台。通风阻力测定分2组,采用双测点同步法

15、测量。数据处理通过编制的V EN T O P T 软件运行计算而得。计算中选用的权矩阵是由测量值的倒数构成的。西石门铁矿中区通风系统网络图如图2所示,通过V EN T O P T 软件计算平差值如表1所示。图2西石门铁矿中区通风系统网络图表1矿井通风网络实测值与平差值对比表分支号巷道名称始点末点长度/m断面/m2风量/m 3s -1阻力/N 风阻/N s 2m 8实测平差实测平差实测平差120m 下盘7205101512301000100001000010000100001000701中国矿业第15卷续表分支号巷道名称始点末点长度/m断面/m2风量/m3s-1阻力/N风阻/Ns2m8实测平差实

16、测平差实测平差由表1可知,各条巷道的风阻值均在合理范围内,风量、阻力均为0的巷道表示该巷道已安装风门或封堵,风量为负值表明风流方向与设定的始点2末点方向相反。同时可以看出实测数据和平差结果存在差异,计算机平差后的结果符合通风网络风流流动基本规律:即流入与流出各节点的空气质量流量达到平衡,各独立回路的风压代数和等于零。合理地消除了测量过程中的随机误差,为后续的风网解算和优化调节奠定了合理的基础。5结论(1平差是以测量数据无粗差,只有随机误差,随机误差成正态分布为基础,其平差计算公式由最小二乘法准则导出,对于具有较大误差的数据应首先人为剔除,避免影响数据平差精度。(2在平差处理中的条件平差和间接平差是可以相互转换的,不过就计算机编程实现难易程度而言,使用条件平差方法模型比较容易实现。(3对于没有进行全风网测量的情况下,在建立平差方程的时候,关键是要求确定平差方程的个数,和各类矩阵所代表的实际意义。(4平差计算中测量精度权矩阵有多种取法。不同的权取值获得的平差结果有差异,通常取风量、阻力或风阻测量值的倒数为精度权矩阵。但是不论取何种权矩阵,其平差后的结果均会满足通风网络风流流动基本规律。利用风量平差结果和风阻测量值倒数为精度权矩阵的阻力平差结果,可计算出风阻平差值。(5通过程