1.一种基于实时监测风量的矿井通风系统异常的诊断方法，其特征在于，包括以下步骤：步骤一、确定矿井通风系统中通风机异常对分支风量的影响通风网路图中从进风井口沿着风流方向到达回风井口的节点和分支的交替序列称为有向通路；

包含通风机i的所有有向通路的集合称之为通风机i子网路，记为GI；

当通风系统中只有通风机i发生异常、其他通风机和通风网路的风阻无异常时，子网路GI中的分支风量随通风机i风量的增加而增加、减少而减少，即：(q′fi-q0fi)/(q′j-q0j)＞0(i＝1,2,…,f；j∈GI)；

当通风系统中只有一台通风机时，f＝1，通风机风量的变化率与通风网路中各分支风量的变化率相同，即：(q′fi-q0fi)/q0fi＝(q′j-q0j)/q0j＝常数；

其中，

i：通风机编号；

GI：包含通风机i的子网路；

f：通风机数量；

j：分支编号；

q0fi：通风正常时通风机i的风量；

q′fi：通风异常时通风机i的风量；

q0j：通风正常时分支j的风量；

q′j：通风异常时分支j的风量；

步骤二、确定矿井通风系统中分支异常对矿井通风系统风量分布的影响通风系统中分支j异常时，即通风系统中分支j风阻发生变化时，分支j本身的风量变化最大，离分支j越远的分支风量变化越小；

两条分支的距离等于两条分支间的分支数；

步骤三、矿井通风系统异常诊断约定矿井通风系统中仅有一处异常；

步骤3.1、风速传感器的实时监测风量正常通风时的分支风阻R0＝(r01,r02,…,r0j…,r0b)，设置风速传感器的分支S＝(s1,s2,…,sj,…,sv)，风速传感器的分支的实时监测风量和通风机的实时监测风量分别为Q′s＝(q′s1,q′s2,…,q′sj,…,q′sv)和Q′f＝(q′f1,q′f2,…,q′fi,…,q′ff)；

其中，r0j：分支j的正常风阻；

b：通风网路中的分支数；

sj：安设传感器j的分支编号；

v：设置风速传感器的分支数；

q′sj：风速传感器的分支sj的实时监测风量；

q′fi：通风机i的实时监测风量；

步骤3.2、通风系统正常

如果所有风速传感器的分支的实时监测风量都等于正常风量，即q′sj＝q0sj(j＝1,2,…,v)；

则通风系统正常；否则通风系统异常，转步骤3.3；

其中，

q0sj：风速传感器的分支S中sj的正常风量；

步骤3.3、通风机异常诊断

步骤3.3.1

将风速传感器的分支作为固定风量分支、风速传感器的分支的实时监测风量Q′s＝(q′s1,q′s2,…,q′sv)作为固定风量、风速传感器的分支之外的其他分支以r0j为风阻，进行通风网路解算以求通风网路的风量分布Q′＝(q′1,q′2,…,q′b)；

步骤3.3.2

求|q′fi-q0fi|(i＝1,2,…,f)最大的通风机编号，并将该通风机编号记为w；

步骤3.3.3

如果

(q′fw-q0fw)/(q′j-q0j)＞0(j∈GW)；

则通风机w异常；否则为通风网路中的分支异常，转步骤3.4；

其中，

w：发生异常的通风机；

GW：包含通风机w的子网路；

q′fw：通风机w的实时监测风量；

q0fw：通风机w的正常风量；

步骤3.4、分支异常诊断

步骤3.4.1

求(q′j-q0j)/(q′fw-q0fw)(j∈GW)最大的分支编号,并将该分支编号记为k；

步骤3.4.2

求离分支k最近的风速传感器的分支，并将该风速传感器的分支编号记为sd；

步骤3.4.3

以分支k的风量q′k为固定风量、其他分支以r0j为风阻求风量分布Q＝(q1,q2,…,qb)和分支k的风阻r′k；

步骤3.4.4

如果|q′sj-qsj|≤ε(j＝1,2,…,v)，则分支k为异常分支，r′k为异常分支k的异常风阻，计算结束；

步骤3.4.5

如果|q′sd-qsd|＞ε，则以q′sd为固定风量、r′k作为分支k的风阻求分支k的风量q′k，转步骤3.4.3；

步骤3.4.6

如果|q′sd-qsd|≤ε，但|q′sj-qsj|＞ε(j∈S；j≠sd)，则分支k不是异常分支，从与分支k连接的分支中选取(q′j-q0j)/(q′fw-q0fw)(j∈GW；j≠k)最大的分支，并将该分支的编号记为k，转步骤3.4.3；

其中，

k：异常分支编号；

qsj:分支k的风量q′k为固定风量时风速传感器的分支sj的计算风量值；

q′k：异常分支k的异常风量；

r′k：异常分支k的异常风阻；

sd：离分支k最近的风速传感器的分支编号；

q′sd：风速传感器实时监测的分支sd的风量；

qsd:分支k的风量q′k为固定风量时风速传感器的分支sd的计算风量值；

ε：设定的计算精度。