一种基于模糊推理的管道泄漏定位方法
【专利摘要】本发明提供一种基于模糊推理的管道泄漏定位方法，包括：实时采集管道首末两端的压力信号和温度信号；若实时采集的管道两端的压力信号中存在非工况调整引起的压力突变信号，则当前管道存在泄漏点，计算压力突变信号到达管道首末两端的时间差；利用负压波法、广义递归神经网络定位方法、利用温度梯度法，定位管道泄漏点到管道首端的距离；根据利用负压波法、广义递归神经网络定位方法、温度梯度法定位得到的管道泄漏点到管道首端的距离，利用模糊推理算法对泄漏点进行定位。采用三种不同的定位方法对管道泄漏点进行定位，经过模糊推理，鱼群算法隶属度函数参数自适应优化，赋予三种定位结果不同的置信度，从而实现泄漏点的精确定位。
【专利说明】
一种基于模糊推理的管道泄漏定位方法
【技术领域】
[0001]本发明属于管道检测【技术领域】，具体涉及一种基于模糊推理的管道泄漏定位方法。
【背景技术】
[0002]随着输油管道的增多和管道服役时间的增长，管道运输的安全状态不容乐观。很多管道的老化程度比较严重，且已进入泄漏事故多发期。此外，人为破坏管道窃取石油资源的现象也相当严重，带来重大安全隐患。因此，管道泄漏检测成为管道安全生产管理的重要工作内容，对管道泄漏点进行定位也成了工作中的重中之重。
[0003]目前，已经出现了温度梯度法、负压波法、管道模型法等多种方法实现对管道泄漏点的定位，但是，在这些方法中，有些方法简单方便、实用性强，但是定位精度较差；而另外一些方法定位精度较准确，但是方法复杂，需要考虑多种因素。因此，设计一种定位精度高、实用性强的管道泄漏定位装置具有重要的现实意义。
【发明内容】
[0004]针对现有技术存在的问题，本发明提供一种基于模糊推理的管道泄漏定位方法。
[0005]一种基于模糊推理的管道泄漏定位方法，包括以下步骤:
[0006]步骤1:实时采集管道首末两端的压力信号和温度信号；
[0007]步骤2:判断实时采集的管道两端的压力信号中是否存在压力突变信号，是，则执行步骤3，否则返回步骤I ;
[0008]步骤3:判断当前管道是否存在工况调整，是，则返回步骤1，否则，当前管道存在泄漏点，计算压力突变信号到达管道首末两端的时间差；
[0009]步骤4:利用负压波法，根据压力突变信号到达管道首末两端的时间差，定位管道泄漏点到管道首端的距离；
[0010]步骤5:利用广义递归神经网络定位方法，根据压力突变信号到达管道首末两端的时间差，定位管道泄漏点到管道首端的距离；
[0011]步骤6:利用温度梯度法，根据采集到的管道首末两端的温度信号，定位管道泄漏点到管道首端的距离；
[0012]步骤7:根据利用负压波法、广义递归神经网络定位方法、温度梯度法定位得到的管道泄漏点到管道首端的距离，利用模糊推理算法对泄漏点进行定位；
[0013]步骤7.1:根据历史数据建立模糊推理规则，用来描述管道的决策条件空间与其对应的决策结果空间的关系；
[0014]所述管道的决策条件空间包括管道首末两端的温度差和管道首末两端压力差；
[0015]所述决策结果空间包括管道运行正常、管道发生泄漏和管道工况调整，其中，管道发生泄漏决策结果的三个程度模糊子集分别为严重泄漏、较大泄漏和小泄漏；
[0016]步骤7.2:对各管道的决策条件空间和各决策结果分别建立模糊隶属度函数；
[0017]步骤7.3:利用鱼群算法优化各模糊隶属度函数的顶点位置和模糊隶属度函数宽度；
[0018]步骤7.3.1:将各决策条件的模糊隶属度函数的顶点位置和模糊隶属度函数宽度、决策结果的模糊隶属度函数的顶点位置和模糊隶属度函数宽度作为基因，并将各基因任意排列构成人工鱼个体形式；
[0019]步骤7.3.2:按照构成的人工鱼个体形式随机产生N组数据，即N个人工鱼个体；
[0020]步骤7.3.3:分别计算N个人工鱼个体的食物浓度，记录食物浓度最大值及其对应的人工鱼个体；
[0021]步骤7.3.4:设定人工鱼个体的感知距离、人工鱼个体移动的最大步长、拥挤度因子、寻优最大次数、食物浓度阈值；
[0022]步骤7.3.5:对人工鱼个体进行寻优；
[0023]步骤7.3.5.1:确定与当前人工鱼个体间距离小于感知距离的所有人工鱼个体，并确定这些人工鱼个体对应的模糊隶属度函数的顶点位置的平均值及模糊隶属度函数宽度的平均值，同时确定这些人工鱼个体中食物浓度最大的人工鱼个体及其食物浓度；
[0024]步骤7.3.5.2:根据聚群、追尾和觅食确定新的人工鱼个体；
[0025]步骤7.3.5.3:计算新的人工鱼个体的食物浓度，判断该新的人工鱼个体的食物浓度是否大于记录的食物浓度最大值，是，则记录该新的人工鱼个体的食物浓度及其对应的人工鱼个体；否则，维持记录的食物浓度最大值；
[0026]步骤7.3.5.4:判断是否所有人工鱼个体均完成步骤7.3.5.1?步骤7.3.5.3:是，则执行步骤7.3.5.5，否则对下一条人工鱼个体执行步骤7.3.5.1?步骤7.3.5.3 ；
[0027]步骤7.3.5.5:判断是否达到寻优最大次数或者当前记录的食物浓度最大值大于食物浓度阈值，则当前记录的食物浓度最大值对应的人工鱼个体为最优人工鱼个体，进而确定优化的各模糊隶属度函数的顶点位置和模糊隶属度函数宽度，否则，返回步骤7.3.5.1 ；
[0028]步骤7.4:根据确定的优化的各模糊隶属度函数的顶点位置和模糊隶属度函数宽度，得到管道发生泄漏决策结果对应的模糊隶属度函数，进而得到严重泄漏的置信度、较大泄漏的置信度和小泄漏的置信度；
[0029]步骤7.5:将严重泄漏的置信度、较大泄漏的置信度和小泄漏的置信度分别作为利用负压波法、温度梯度法、广义递归神经网络定位方法定位得到的管道泄漏点到管道首端的距离的置信度，计算管道泄漏点到管道首端的距离&:
COfirK, +coJ\x,+COfl-W
[0030]xL =-广 i f 丄 r~
?)/丨)十 α)/丨-TCOf2
[0031]其中，cof^、Cof1和Cof2分别为利用负压波法、广义递归神经网络定位方法、温度梯度法定位得到的管道泄漏点到管道首端的距离的置信度，X1 > X2和X3分别为利用负压波法、广义递归神经网络定位方法、温度梯度法定位得到的管道泄漏点到管道首端的距离。
[0032]有益效果:
[0033]采用三种不同的定位方法对管道泄漏点进行定位，利用压力信号到达管道首末端的时间差实现了基于管道机理模型的负压波定位和基于数据驱动方法的广义递归神经网络模型定位，同时，利用管道首末端的温度信号实现了基于温度梯度的定位；对上述三种定位结果进行二次定位，经过模糊推理，鱼群算法隶属度函数参数自适应优化，赋予三种定位结果不同的置信度，从而实现泄漏点的精确定位。
【专利附图】
【附图说明】
[0034]图1为本发明一种实施方式的基于模糊推理的管道泄漏定位装置结构框图；
[0035]图2为本发明一种实施方式的第一信号调理模块的电路原理图；
[0036]图3为本发明一种实施方式的A/D转换芯片与FPGA的接口电路图；
[0037]图4为本发明一种实施方式的基于模糊推理的管道泄漏定位方法流程图；
[0038]图5为本发明一种实施方式模糊推理方法流程图；
[0039]图6为本发明一种实施方式优化运行状态和决策结果的隶属度函数参数优化流程图。
【具体实施方式】
[0040]下面结合附图对本发明的【具体实施方式】做详细说明。
[0041]本实施方式中，实现基于模糊推理的管道泄漏定位方法所采用的基于模糊推理的管道泄漏定位装置，如图1所示，该装置包括压力传感器模块、温度传感器模块、第一信号调理模块、第二信号调理模块、第一 A/D转换模块、第二 A/D转换模块和FPGA中央处理单元，其中FPGA中央处理单元模块包括时序控制模块和基于模糊推理的泄漏定位模块，时序控制模块用于产生A/D转换模块的工作时序，基于模糊推理的泄漏定位模块包括负压波定位模块、广义递归神经网络定位模块、温度梯度定位模块和模糊推理模块。
[0042]压力传感器模块的型号为PT500-502，温度传感器模块的型号为PT100，第一 A/D转换模块和第二 A/D转换模块的型号均为ADS7844，FPGA中央处理单元模块的型号为EP3C25Q240C8。
[0043]压力传感器和温度传感器均安装在管道首末两端，分别采集管道首端和末端的压力和温度信号，压力传感器输出端连接第一信号调理模块的输入端，温度传感器输出端连接第二信号调理模块的输入端，第一信号调理模块的输出端、第二信号调理模块的输出端分别连接第一 A/D转换模块的输入端和第一 A/D转换模块的输入端，第一 A/D转换模块的输出端、第二 A/D转换模块的输出端分别连接FPGA中央处理单元的输入端，FPGA中央处理单元的时序控制模块分别连接第一 A/D转换模块和第二 A/D转换模块。
[0044]第一信号调理模块的电路原理图如图2所示，该模块实现信号的滤波和放大，压力传感器的输出首先经过滤波电路滤波，然后经一个1K的电阻R2连接到运算放大器的反相输入端，同相输入端接2.5V的参考电压，运算放大器的输出端连接电阻R3的一端、电阻Rl的一端及电容C2的一端，电阻R3的另一端作为信号调理模块的输出端连接A/D转换芯片的输入端，电阻Rl的另一端连接运算放大器的反相输入端，电容C2的另一端接地。本实施方式中运算放大器的型号为AD824 ο第二信号调理模块的原理与第一信号调理模块相同。
[0045]第一 A/D转换模块/第二 A/D转换模块与FPGA的接口电路图如图3所示，第一A/D转换模块、第二 A/D转换模块均采用ADS7844的A/D转换芯片，A/D转换芯片将电压信号转换为数字信号，A/D转换芯片的6个不同的输出端分别连接FPGA时序控制模块的自定义I/O 口，即A/D转换芯片的DCLK端连接FPGA的I/0.23端，A/D转换芯片的CS端连接I/0.24端、A/D转换芯片的Din端连接I/0.25端、A/D转换芯片的Busy端连接I/0.26端，A/D转换芯片的Dout端连接I/0.27端，FPGA的型号为EP3C25Q240C8。
[0046]本实施方式中，基于模糊推理的管道泄漏定位方法，如图4所示，包括以下步骤:
[0047]步骤1:实时采集管道首末两端的压力信号和温度信号；
[0048]步骤2:判断实时采集的管道两端的压力信号中是否存在压力突变信号，是，则执行步骤3，否则返回步骤I ;
[0049]步骤3:判断当前管道是否存在工况调整，是，则返回步骤1，否则，当前管道存在泄漏点，计算压力突变信号到达管道首末两端的时间差；
[0050]步骤4:利用负压波法，根据压力突变信号到达管道首末两端的时间差，定位管道泄漏点到管道首端的距离X1 ;
[0051]X1=I(Zjv-Af)
2
[0052]其中，L为管道的长度，V为管道输油介质中压力波的传播速度，AC为管道实时压力突变信号到达管道首末两端的时间差；
[0053]步骤5:利用广义递归神经网络定位方法，根据压力突变信号到达管道首末两端的时间差，定位管道泄漏点到管道首端的距离X2 ;
[0054]基于广义递归神经网络模型建立的理论时间差Λ t与泄漏位置X的对应关系Z =F(At),利用压力突变信号到达管道首、末端的时间差，定位出管道泄漏点到管道首端的距离X2 ；
[0055]步骤5-1:从管段首端开始，到末端为止(全长为L)，构造模拟泄漏位置序列Z ={O, Z1, Z21,…，ZJ，I为单位步长，I越小，定位精度越高；
[0056]步骤5-2:计算负压波在管道介质中的传播速度；
[0057]V (Z) = {v (O)，V (Z1)，V (Z21)，…，V (Zl) }
[0058]和管道内传输介质的流速
[0059]u (Z) = {u (O), u (Z1), u (Z21),..., u (Zl) }；
[0060]步骤5-3:构建模拟泄漏位置对应的理论时间差序列Λ t = { Λ ^ Λ Λ t21,…，AtL}，其中
?>ΧIi./.1
「00611 Δ/ = -dx — -dx
L 」Ji) V(Z) -U(Z)^Z)-hu(Z)
[0062]步骤5-4:建立基于模拟泄漏位置对应的理论时间差序列的广义递归神经网络模型，利用历史数据进行训练；
[0063]步骤5-5:将管道实时压力突变信号到达管道首末两端的时间差作为广义递归神经网络模型的输入，广义递归神经网络模型的输出为管道泄漏点到管道首端的距离x2。
[0064]步骤6:利用温度梯度法，根据采集到的管道首末两端的温度信号，定位管道泄漏点到管道首端的距离X3 ;
[0065]根据管道首端的温度数据、管道的流量信息和管道固有特性构建管道泄漏温度下降曲线1，再根据管道末端的温度数据、管道的流量信息和管道固有特性构建管道泄漏温度上升曲线2，利用曲线I和曲线2的交点定位出管道泄漏点到管道首端的距离X3 ；
[0066]具体步骤如下:
[0067]步骤6-1:构建管道泄漏温度下降曲线1，公式如下:
KinD ‘Ki π?)、 / KnDX、
[0068]t2 ^t0--—— + (广丨 ~tQ + —exp(----)
CC.Kj C.
[0069]其中，t2表示到首端距离为X处的温度，t0表示X处周围介质的温度，K表示介质传热系数，i表示管道的水力坡降，D表示管道直径，C表示比热容，h表示管道首端温度，G表示流量。
[0070]步骤6-2:构建管道泄漏温度上升曲线2，公式如下:
r nKinD tKhiD、 t KnDX ,
[0071]ι4 = l() + ■ (I W1--^^)CXP(———)
[0072]其中，t4表示到末端距离为X处的温度，t3表示管道末端温度。
[0073]步骤6-3:利用管道泄漏温度下降曲线I和管道泄漏温度上升曲线2的交点，定位出管道泄漏点距离管道首段的距离x3。
[0074]步骤7:根据利用负压波法、广义递归神经网络定位方法、温度梯度法定位得到的管道泄漏点到管道首端的距离，利用模糊推理算法对泄漏点进行定位；
[0075]如图5所示，具体步骤如下:
[0076]步骤7.1:根据历史数据建立模糊推理规则，用来描述管道的决策条件空间与其对应的决策结果间的关系；
[0077]管道的决策条件空间包括管道首末两端的压力差和温度差；
[0078]决策结果空间包括管道运行正常、管道发生泄漏和管道工况调整，其中，管道发生泄漏决策结果的三个泄漏程度模糊子集分别为严重泄漏、较大泄漏和小泄漏；
[0079]模糊推理规则Ri的文字描述形式如下:
[0080]Ri:If运行状态属于决策条件空间的子集，
[0081]then当前管道运行状况的决策结果是...。
[0082]管道的决策条件空间C = {Cl，c2}，其中C1代表首末两端温度差，C2代表首末两端压力差，D表示决策结果空间，决策条件C」(j = I, 2)表示管道运行状态,D= {(^，‘屯},其中Cl1代表运行正常，d2代表发生泄漏，d3代表调节工况。
[0083]步骤7.2:对各管道的决策条件空间和各决策结果分别建立模糊隶属度函数；
[0084]决策条件空间中每一类决策条件表示的运行状态被划分为5个程度的模糊子集，当前运行状态属于每一个模糊子集的隶属度为E = {e0, e1；..., e4}, ek(k = 0, I,...4)代表由该决策条件表示的运行状态偏离其正常值的程度；每一类决策结果同样被划分为3个程度模糊子集，每一类决策结果属于其对应的模糊子集的隶属度为F = {f0, f1； f2}，fk(k =0，I, 2)表示该类决策结果可能性的程度。
[0085]步骤7.3:利用鱼群算法优化各模糊隶属度函数的顶点位置和模糊隶属度函数宽度；
[0086]如图6所示，具体步骤如下:
[0087]步骤7.3.1:将各决策条件的模糊隶属度函数的顶点位置和模糊隶属度函数宽度、决策结果的模糊隶属度函数的顶点位置和模糊隶属度函数宽度作为基因，并将各基因任意排列构成人工鱼个体形式；
[0088]将运行状态Cj的模糊隶属度函数进行实数编码为Wj1Wf Vj5Wj5作为一个基因,Vji=(i = I, 2…5)表示运行状态Cj中第i个隶属度函数的顶点位置,Wji表示运行状态Cj中第i个隶属度函数宽度；将决策结果火的模糊隶属度函数进行实数编码为:CrtZrtCrfZrfCriZri作为一个基因，cri (i = I, 2，3)表示决策结果火中第i个隶属度函数的顶点位置，Zri表示决策结果4中第i个隶属度函数宽度。将由运行状态得到的2个基因同由决策结果得到的3个基因任意排列构成一个人工鱼个体形式X。
[0089]步骤7.3.2:按照构成的人工鱼个体形式随机产生N组数据，即N个人工鱼个体，N=50 ；
[0090]每组数据均在其所对应的隶属度函数论域内，则每一组数据可视为一个人工鱼Xi, i = 1,2, - ,50 ；
[0091]步骤7.3.3:分别计算N个人工鱼个体Xi的食物浓度，记录食物浓度最大值及其对应的人工鱼个体；
[0092]食物浓度Y (Xi)公式如下:
[0093]./ I 1-/- , 5Hllll(VV)1., VV..),.1mill(U;'、, U;.,)
[0094]步骤7.3.4:设定人工鱼个体的感知距离T、人工鱼个体移动的最大步长St印、拥挤度因子P、寻优最大次数NUM为200、食物浓度阈值Tm ；
[0095]步骤7.3.5:对人工鱼个体进行寻优；
[0096]步骤7.3.5.1:确定与当前人工鱼个体间距离小于感知距离的所有人工鱼个体，并确定这些人工鱼个体对应的模糊隶属度函数的顶点位置的平均值及模糊隶属度函数宽度的平均值，同时确定这些人工鱼个体中食物浓度最大的人工鱼个体及其食物浓度；
[0097]令Clij = I I X1-Xj I I表示人工鱼Xi和Xj之间的距离；
[0098]确定人工鱼Xi的邻域(即Clij < T)内的人工鱼条数s (包括Xi本身)，即确定与当前人工鱼个体间距离小于感知距离的所有人工鱼个体，将s条人工鱼任意排列，表示为:
[0099]X”X2，…，Xi，…Xs，取P 为:
Γη1ηη1Σν^/1 Σ wM Σν^5Σ zJki Σ6>3 Σ Zfti甘
LU I υυ」f=i/=ι /=i/=i /=i/=i /=i共
?.，-- * * 595 * ' ?95***55* * * 9
SSSSSSSS
中，P表示此邻域Clij < T鱼群的中心位置，Vfji (i = 1,2,…5)表示第f条人工鱼Xf的运行状态h中第i个隶属度函数的顶点位置，表示第f条人工鱼Xf的运行状态h中第i个隶属度函数宽度，Cfki表示第f条人工鱼Xf的决策结果dk中第i个隶属度函数的顶点位置，Zfki表示第f条人工鱼Xf的决策结果dk中第i个隶属度函数宽度。选取X1J2,…，Xi，...Xs中食物浓度最大的人工鱼，标记为Xmax,分别计算Xmax的食物浓度Y(Xmax)和中心位置的食物浓度Y (P)。
[0101]步骤7.3.5.2:根据聚群、追尾和觅食确定新的人工鱼个体；
[0102](I)聚群
[0103]如果满足Y (P)/(s-Ι) > P.Y(Xi)且Y(Xi) < Y(P)的条件时，则根据
P-X
[0104]^ ' =X, +rand()-step-^p_x^ 得到新的个体鱼 Xinext。(rand()为随机产生的
O?I范围内的随机数)。
[0105](2)追尾
[0106]如果满足Y (Xmax) / (s-Ι) > P.Y (Xmax)且 Y (Xi) < Y (Xmax)的条件时，则根据\ = Xi + rand0.step.: ^ ^得到的新人工鱼Xinext。如果聚群和追尾的条件同时满足，则比较Y (Xmax)和Y (P)，按照值较大的方式进行移动。
[0107](3)觅食
[0108]如果⑴和⑵两个条件都不满足，则在人工鱼Xi的邻域内(即Clij < T)中随机
X — X
选择一个人工鱼Xj，若Y(Xi) <Y(Xj)，根据An = Xi +randO-step.^得到新的人工鱼Ximxt ;否则重新随机选择重复判断，若反复5次后仍不满足条件，则按照
[0109]Xinext = Xj+rand ().step 随机得到 Xinexto
[0110]计算新的人工鱼Xinext的食物浓度Y (Xinext)。
[0111]步骤7.3.5.3:计算新的人工鱼个体的食物浓度，判断该新的人工鱼个体Xinrait的食物浓度Y(Xinrart)是否大于记录的食物浓度最大值，是，则记录该新的人工鱼个体的食物浓度及其对应的人工鱼个体；否则，维持记录的食物浓度最大值；
[0112]步骤7.3.5.4:判断是否所有50条人工鱼个体均完成步骤7.3.5.1?步骤7.3.5.3:是，则执行步骤7.3.5.5，否则对下一条人工鱼个体执行步骤7.3.5.1?步骤
7.3.5.3 ;
[0113]步骤7.3.5.5:判断是否达到寻优最大次数200或者当前记录的食物浓度最大值大于食物浓度阈值Tni，则当前记录的食物浓度最大值对应的人工鱼个体为最优人工鱼个体，进而确定优化的各模糊隶属度函数的顶点位置和模糊隶属度函数宽度，否则，将新个体鱼作为Xi，返回步骤7.3.5.1 ；
[0114]步骤7.4:根据确定的优化的各模糊隶属度函数的顶点位置和模糊隶属度函数宽度，得到管道发生泄漏决策结果对应的模糊隶属度函数，进而得到严重泄漏的置信度、较大泄漏的置信度和小泄漏的置信度；
[0115]步骤7.5:将严重泄漏的置信度、较大泄漏的置信度和小泄漏的置信度分别作为利用负压波法、温度梯度法、广义递归神经网络定位方法定位得到的管道泄漏点到管道首端的距离的置信度，计算管道泄漏点到管道首端的距离&:
cof'x, + cof\x, + cof,.\\
[0116]Xr————^
?.ο/0 + Cf；/, +COf2
[0117]其中，cof^、Cof1和Cof2分别为利用负压波法、广义递归神经网络定位方法、温度梯度法定位得到的管道泄漏点到管道首端的距离的置信度，X1 > X2和X3分别为利用负压波法、广义递归神经网络定位方法、温度梯度法定位得到的管道泄漏点到管道首端的距离。
[0118]决策结果空间中的d2代表发生泄漏，其3个泄漏程度模糊子集为:严重泄漏(&)、较大泄漏(f\)和小泄漏(f2)，其相应的置信度分别为COf^COf1和COf2，分别将该三个置信度等同于利用负压波法、广义递归神经网络定位方法定位、温度梯度法得到的管道泄漏点到管道首端的距离的置信度，进而完成泄漏点定位。
【权利要求】
1.一种基于模糊推理的管道泄漏定位方法，其特征在于:包括以下步骤: 步骤1:实时采集管道首末两端的压力信号和温度信号； 步骤2:判断实时采集的管道两端的压力信号中是否存在压力突变信号，是，则执行步骤3，否则返回步骤I ; 步骤3:判断当前管道是否存在工况调整，是，则返回步骤1，否则，当前管道存在泄漏点，计算压力突变信号到达管道首末两端的时间差； 步骤4:利用负压波法，根据压力突变信号到达管道首末两端的时间差，定位管道泄漏点到管道首端的距离； 步骤5:利用广义递归神经网络定位方法，根据压力突变信号到达管道首末两端的时间差，定位管道泄漏点到管道首端的距离； 步骤6:利用温度梯度法，根据采集到的管道首末两端的温度信号，定位管道泄漏点到管道首端的距离； 步骤7:根据利用负压波法、广义递归神经网络定位方法、温度梯度法定位得到的管道泄漏点到管道首端的距离，利用模糊推理算法对泄漏点进行定位； 步骤7.1:根据历史数据建立模糊推理规则，用来描述管道的决策条件空间与其对应的决策结果空间的关系； 所述管道的决策条件空间包括管道首末两端的温度差和管道首末两端压力差； 所述决策结果空间包括管道运行正常、管道发生泄漏和管道工况调整，其中，管道发生泄漏决策结果的三个程度模糊子集分别为严重泄漏、较大泄漏和小泄漏； 步骤7.2:对各管道的决策条件空间和各决策结果分别建立模糊隶属度函数； 步骤7.3:利用鱼群算法优化各模糊隶属度函数的顶点位置和模糊隶属度函数宽度；步骤7.3.1:将各决策条件的模糊隶属度函数的顶点位置和模糊隶属度函数宽度、决策结果的模糊隶属度函数的顶点位置和模糊隶属度函数宽度作为基因，并将各基因任意排列构成人工鱼个体形式； 步骤7.3.2:按照构成的人工鱼个体形式随机产生N组数据，即N个人工鱼个体； 步骤7.3.3:分别计算N个人工鱼个体的食物浓度，记录食物浓度最大值及其对应的人工鱼个体； 步骤7.3.4:设定人工鱼个体的感知距离、人工鱼个体移动的最大步长、拥挤度因子、寻优最大次数、食物浓度阈值； 步骤7.3.5:对人工鱼个体进行寻优； 步骤7.3.5.1:确定与当前人工鱼个体间距离小于感知距离的所有人工鱼个体，并确定这些人工鱼个体对应的模糊隶属度函数的顶点位置的平均值及模糊隶属度函数宽度的平均值，同时确定这些人工鱼个体中食物浓度最大的人工鱼个体及其食物浓度； 步骤7.3.5.2:根据聚群、追尾和觅食确定新的人工鱼个体； 步骤7.3.5.3:计算新的人工鱼个体的食物浓度，判断该新的人工鱼个体的食物浓度是否大于记录的食物浓度最大值，是，则记录该新的人工鱼个体的食物浓度及其对应的人工鱼个体；否则，维持记录的食物浓度最大值； 步骤7.3.5.4:判断是否所有人工鱼个体均完成步骤7.3.5.1?步骤7.3.5.3:是，则执行步骤7.3.5.5，否则对下一条人工鱼个体执行步骤7.3.5.1?步骤7.3.5.3 ； 步骤7.3.5.5:判断是否达到寻优最大次数或者当前记录的食物浓度最大值大于食物浓度阈值，则当前记录的食物浓度最大值对应的人工鱼个体为最优人工鱼个体，进而确定优化的各模糊隶属度函数的顶点位置和模糊隶属度函数宽度，否则，返回步骤7.3.5.1 ；步骤7.4:根据确定的优化的各模糊隶属度函数的顶点位置和模糊隶属度函数宽度，得到管道发生泄漏决策结果对应的模糊隶属度函数，进而得到严重泄漏的置信度、较大泄漏的置信度和小泄漏的置信度； 步骤7.5:将严重泄漏的置信度、较大泄漏的置信度和小泄漏的置信度分别作为利用负压波法、温度梯度法、广义递归神经网络定位方法定位得到的管道泄漏点到管道首端的距离的置信度，计算管道泄漏点到管道首端的距离&:
τ _ cof,r\ + COf2Xs 十 cyj/;.\:2
L CQfn+cqf\+COf2 其中，COf^、COf1和COf2分别为利用负压波法、广义递归神经网络定位方法、温度梯度法定位得到的管道泄漏点到管道首端的距离的置信度，X1^X2和X3分别为利用负压波法、广义递归神经网络定位方法、温度梯度法定位得到的管道泄漏点到管道首端的距离。
【文档编号】F17D5/02GK104197203SQ201410334406
【公开日】2014年12月10日 申请日期:2014年7月14日 优先权日:2014年7月14日
【发明者】张化光, 吴振宁, 汪刚, 刘金海, 冯健, 马大中, 王一, 宁迪, 李芳明, 卢森骧, 许相凯, 刘喆, 周坤, 赵重阳 申请人:东北大学