基于多传感器信息融合的长输管线泄漏检测方法与流程

本发明涉及管线泄漏检测技术领域，具体涉及一种基于多传感器信息融合的长输管线泄漏检测方法。

背景技术：

长输管线不仅包含对石油和天然气等石化产品的输送，还包含对液氨、低温甲醇、低温乙二醇等化学品的输送。无论长输管线经过偏远的山区、人烟稀少的平原，还是穿过人口密集的生活区，管线内介质的泄漏不仅会造成环境污染，还带来经济损失，甚至影响到周边居民的人身财产安全。因此，对长输管线泄漏的可靠检测是一项极具意义的系统工程。

长输管线泄漏的原因，既有不法份子盗油盗气、施工单位野蛮施工等第三方破坏，也有地震、滑坡、泥石流等自然灾害，还有防护层老化、焊接质量缺陷等管道自身缺陷。目前常用的泄漏检测方法有负压波法、次声波法和分布式光纤检测法等。其中负压波法基于压力对泄漏位置进行定位，针对第三方破坏引起的泄漏，压力变化明显，适用性较好，但是当管道输送介质为气体，或管道自身存在渗漏的情况下，压力变化微小，负压波法无法有效检测。次声波法通过检测泄漏介质与泄漏处管壁的高速摩擦产生的次声波信号进行定位检测，其定位方法同负压波法类似，由于输气管线泄漏时同样产生次声波信号，因此在输气管线同样适用。但是，现场发生多点泄漏时，不同泄漏点产生的次声波会相互干扰，引起传感器无法分辨，从而导致定位不准。分布式光纤测温法沿管线敷设感温光纤，感温光纤上每隔0.5m有一个测温点，当管网有泄漏时，泄漏介质引起周边温度场的改变，通过测温报警点即可对泄漏处进行定位。分布式光纤测温法不仅能解决管道渗漏引起的泄漏检测，而且施工方便、易维护。由于长输管线敷设路径长(一般有几十至上百公里不等)，同一时间途径地域的温度不尽相同；且有些地域本身的气候变化显著，昼夜温差较大。因此，单纯的分布式光纤测温法难以应对复杂多变的气候变化，会产生一定的误报警。综上所述，目前的测量方法难以满足对架空敷设长输管线健康状态的有效监测与泄漏定位。

技术实现要素：

本发明的目的就是针对现有技术的缺陷，提供一种基于多传感器信息融合的长输管线泄漏检测方法，其基于流量、温度时间梯度和温度空间梯度三参量异质多传感器信息模糊融合，可降低环境温度和恶劣气候变化等对测量结果造成的影响，检测结果具有较高的可靠性。

本发明一种基于多传感器信息融合的长输管线泄漏检测方法，其技术方案为：在长输管线的首端和末端设置流量计，以及长度方向上均匀设置温度传感器，基于所述流量计和温度传感器进行泄露检测的方法包括

步骤1，实时采集流量数据及温度数据，且每次采集的数据为一帧；

步骤2，基于采集到的当前帧流量数据计算累积流量差值δf，并对所述累积流量差值δf进行模糊处理；

步骤3，基于采集到的当前帧温度数据计算温度时间梯度δti^(t)，并对所述温度时间梯度δti^(t)进行模糊处理；

步骤4，基于采集到的当前帧温度数据计算温度空间梯度δti^(s)，并对所述温度空间梯度δti^(s)进行模糊处理；

步骤5，根据累积流量差值δf、温度时间梯度δti^(t)、温度空间梯度δti^(s)的模糊输入语言变量，结合预先建立的模糊对照库，输出管线发生泄漏的可能性p，对所述输出管线发生泄漏的可能性p进行模糊处理后的管线泄漏概率p(k)进行判断，当p(k)超过阈值pth时，判断为管线泄漏。

较为优选的，还包括

步骤6，当判断为管线泄漏时，泄漏点集合c由泄漏位置集合a和泄漏位置集合b共同确定，且c＝a∩b；

其中，泄漏位置集合a和泄漏位置集合b分别根据温度时间梯度和温度空间梯度的模糊语言值进行确定。

较为优选的，所述基于采集到的当前帧流量数据计算累积流量差值δf包括：

计算首站流量计后n个温度测点的平均值和末站流量计前n个温度测点的平均值；

利用首站流量计后n个温度测点的平均值和末站流量计前n个温度测点的平均值分别对首站和末站瞬时流量进行温度补偿，得到补偿流量；

根据所述补偿流量分别计算首站累积流量和末站累积流量；

根据首站累积流量和末站累积流量计算累积流量差值δf。

较为优选的，所述基于采集到的当前帧温度数据计算温度时间梯度δti^(t)为：

根据公式δti^(t)(k)＝|ti(k)-ti(k-t)|计算温度时间梯度δti^(t)；

其中，ti(k)为管线上第i个测点在第k帧时的温度值，ti(k-t)为管线上第i个测点在第k-t帧时的温度值，δti^(t)(k)为第i个测点在第k帧时的温度时间梯度值。

较为优选的，所述基于采集到的当前帧温度数据计算温度空间梯度δti^(s)为：

根据公式δti^(s)(k)＝|ti(k)-ti+1(k)|计算温度空间梯度δti^(s)；

其中，ti(k)为管线上第i个测点在第k帧时的温度值，ti+1(k)为管线上第i+1个测点在第k帧时的温度值，δti^(s)(k)为第i个测点在第k帧时的温度空间梯度值。

较为优选的，所述模糊对照库为：

当δf＝nb时，p＝s；

当δf＝ns时

若δti^(t)和δti^(s)均取值b，p＝b；

若δti^(t)和δti^(s)一个取值b，一个取值m，p＝m；

其余情况下p＝s；

当δf＝o时

若δti^(t)和δti^(s)均取值b，p＝b；

若δti^(t)和δti^(s)之间有一个或两个取值m，p＝m；

其余情况下p＝s；

当δf＝ps时

若δti^(t)和δti^(s)均取值b，p＝b；

若δti^(t)和δti^(s)之间有一个取值s，p＝s；

其余情况下p＝m；

当δf＝pb时

若δti^(t)和δti^(s)之间有1个取值b，p＝b；

若δti^(t)和δti^(s)之间有一个取值s，p＝s；

其余情况下p＝m；

其中，s、m、b为语言值的大、中、小，nb、ns、o、ps、pb分别为负大、负小、零、正小、正大。

较为优选的，所述步骤2中，基于采集到的当前帧流量数据计算累积流量差值δf之前，对采集到的当前帧流量数据进行去噪处理。

较为优选的，所述累积流量差值δf的模糊处理包括：

设置累积流量差值δf的基本论域、语言值的离散论域、量化因子、语言变量的语言值和对应每个语言值的隶属度函数μ(δf)。

较为优选的，所述温度时间梯度δti^(t)的模糊处理包括：

设置温度时间梯度δti^(t)的基本论域、语言值的离散论域、量化因子、语言变量的语言值和对应每个语言值的隶属度函数μ(δti^(t))。

较为优选的，所述温度空间梯度δti^(s)的模糊处理包括：

设置温度空间梯度δti^(s)的基本论域、语言值的离散论域、量化因子、语言变量的语言值和对应每个语言值的隶属度函数μ(δti^(s))。

本发明的有益效果为：通过流量、温度时间梯度和温度空间梯度三个参量信息进行异质多传感器模糊融合分析，继而进行长输管线的泄漏与定位检测，检测结果可靠。通过温度时间梯度、温度空间梯度等进行泄漏分析，降低环境温度和恶劣气候变化等对测量结果造成的影响；通过流量计前后n个温度测点的平均值分别补偿首站与末站的瞬时流量，减小流量采集误差，提升流量参量的准确性。

附图说明

图1为本发明用于实施基于多传感器信息融合的长输管线泄漏检测方法的系统设置示意图；

图2为本发明基于多传感器信息融合的长输管线泄漏检测方法流程示意图。

具体实施方式

下面结合附图和具体实施例对本发明作进一步的详细说明，便于清楚地了解本发明，但它们不对本发明构成限定。

如图1所示，用于实施基于多传感器信息融合的长输管线泄漏检测方法的系统包括控制器1，容积式流量计2，分布式光纤温度传感器3和dts测温主机4等模块组成；分布式测温光纤4随管线敷设，每间隔0.5m有1个温度测点，通过光通率变化进行温度测量，测温精度可达1℃，测温范围为-180～300℃；dts测温主机3布置在现场，采集分布式测温光纤4每个测点的原始温度，在主机内进行平滑、累加平均和小波去噪等数据处理，然后通过modbus协议或tcp/ip协议将测点温度返回给控制器1；容积式流量计2分布在管线的首、末两端，实时监测流入和流出管线介质的体积流量，流量信号以4～20ma+hart信号形式，通过阻燃多股铜芯聚乙烯绝缘聚氯乙烯护套铜丝编织屏蔽计算机电缆传送至控制器1。控制器1为逻辑控制模块，对采集的温度和流量数据进行分析、预警和存储。

本发明为独立于用户基本过程控制系统之外的第三方子系统，不直接参与工艺过程的操控，通过modbus和泄漏联锁干节点接口为操作人员提供长输管线的健康状态信息。

如图2所示，基于多传感器信息融合的长输管线泄漏检测方法的流程如下：

步骤1，控制器从容积流量计和dts主机采集一帧流量(首站和末站)及温度(每个测温点)数据，设此次信号采集为第k帧，fa(k)为首站在第k帧检测的流量，fb(k)为末站在第k帧检测的流量，ti(k)为第i个测点在第k帧检测到的温度(温度单位为℃)。

步骤2，基于采集到的当前帧流量数据计算累积流量差值δf，并对所述累积流量差值δf进行模糊处理。对首站和末站的原始流量数据分别进行去噪处理；然后分别计算首站/末站流量计前后n个温度测点的平均值，并分别对首站和末站瞬时流量进行温度补偿；之后，分别计算首站与末站的累积流量，继而得到首站和末站的累积流量差值δf；最后，确定流量参量的论域、模糊语言值以及每个模糊语言值对应的隶属度函数，进行流量参量的模糊化处理。

其具体方法如下：

根据采集的第k帧温度信息，计算首站流量计后n个温度测点温度的平均值设fa′(k)为去噪后的瞬时流量，fa″(k)为补偿后瞬时流量，温度补偿按下式进行，式中t0为被测管线的设计温度(温度单位为℃)。末站瞬时流量补偿方法与首站相同，此处不过多赘述。

之后，根据补偿流量分别计算首站和末站的累积流量fa∑(k)和fb∑(k)，并通过公式(2)计算得到第k帧累积流量差值δf(k)：

δf(k)＝fa∑(k)-fb∑(k)(2)

最后，进行流量参量的模糊化。设δf参量的基本论域为[-900m3/h，900m3/h]，该语言值的离散论域为{-3，-2，-1，0，1，2，3}，则量化因子为1/300，语言变量的语言值设为nb,ns,o,ps,pb等5个，对应每个语言值的隶属度函数μ(δf)如下表1所示。

表1δf输入参量隶属度函数

步骤3：基于采集到的当前帧温度数据计算温度时间梯度δti^(t)，并对所述温度时间梯度δti^(t)进行模糊处理。dts主机已对每个测点的温度信号进行平滑、累加平均和小波去噪等分析，控制器在此基础上根据公式(3)依次计算每个测点在周期t的温度时间梯度值δti^(t)(k)；

δti^(t)(k)＝|ti(k)-ti(k-t)|(3)

然后，进行温度时间梯度参量的模糊化。设δti^(t)参量的基本论域为[0℃，10℃]，该语言值的离散论域为{0，1，2，3}，则量化因子为3/10，语言变量的语言值设为s,m,b等3个，对应每个语言值的隶属度函数μ(δti^(t))如下表2所示。

表2δti^(t)和δti^(s)输入参量隶属度函数

步骤4，基于采集到的当前帧温度数据计算温度空间梯度δti^(s)，并对所述温度空间梯度δti^(s)进行模糊处理。控制器根据公式(4)依次计算每个测点在本帧的温度空间梯度值δti^(s)(k)。

δti^(s)(k)＝|ti(k)-ti+1(k)|(4)

然后，进行温度空间梯度参量的模糊化。设δti^(s)参量的基本论域为[0℃，10℃]，该语言值的离散论域为{0，1，2，3}，则量化因子为3/10，语言变量的语言值设为s,m,b等3个，对应每个语言值的隶属度函数μ(δti^(s))如表2所示。

步骤5，根据累积流量差值δf、温度时间梯度δti^(t)、温度空间梯度δti^(s)的模糊输入语言变量，结合预先建立的模糊对照库，输出管线发生泄漏的可能性p，对所述输出管线发生泄漏的可能性p进行模糊处理后的管线泄漏概率p(k)进行判断，当p(k)超过阈值pth时，判断为管线泄漏。

本例模糊规则库如下：

(1)当δf＝nb时，无论δti^(t)和δti^(s)如何取值，p＝s。

(2)当δf＝ns时，当且仅当δti^(t)和δti^(s)均取值b时，p＝b；当δti^(t)和δti^(s)一个取值b，一个取值m时，p＝m；其余情况下p＝s。

模糊规则如表3-1所示

表3-1模糊规则-1

(3)当δf＝o时，当且仅当δti^(t)和δti^(s)均取值b时，p＝b；当δti^(t)和δti^(s)之间有一个或两个取值m时，p＝m；其余情况下p＝s。

模糊规则如表3-2所示

表3-2模糊规则-2

(4)当δf＝ps时，当且仅当δti^(t)和δti^(s)均取值b时，p＝b；当δti^(t)和δti^(s)之间有一个取值s时，p＝s；其余情况下p＝m。

模糊规则如表3-3所示

表3-3模糊规则-3

(5)当δf＝pb时，当δti^(t)和δti^(s)之间有1个取值b时，p＝b；当δti^(t)和δti^(s)之间有一个取值s时，p＝s；其余情况下p＝m。

模糊规则如表3-4所示

表3-4模糊规则-4

由以上模糊规则可以看出，本例采用小(s)、中(m)、大(b)三个语言值描述输出模糊变量p，每个语言值的隶属度函数μ(p)如表4所示。输出变量p的基本论域为[0,1]，对应的离散论域为{0，1，2，3，4}，则其量化因子为4。

表4输出参量p隶属度函数

采用最大隶属度法去模糊化，去模糊化后即可得到第k帧数据分析中管线泄漏概率p(k)。本例规定，当p(k)>0.7时，系统即判定管线泄漏。

步骤6，泄漏点集合c确定。泄漏点集合c由温度时间梯度确定的泄漏集合a和温度空间梯度确定的泄漏集合b共同确定，确定方法如下：

c＝a∩b(5)

对于温度时间梯度参量，其泄漏点集合a根据公式(6)确定，且当所有点的温度时间梯度模糊变量落入离散论域0时，集合a为全集。

a＝{i|δti^(t)＝1或δti^(t)＝2或δti^(t)＝3}(6)

对于温度空间梯度参量，其泄漏点集合b根据公式(7)确定，且当所有点的温度空间梯度模糊变量落入离散论域0时，集合b为全集。

b＝{i|δti^(s)＝1或δti^(s)＝2或δti^(s)＝3}(7)

步骤7，控制器综合步骤5和步骤6的结果作为模糊判定输出，并进行下一帧数据采集，继续执行步骤1至步骤6。

在一个更加优选的实施例中，可根据项目具体实施情况，将管线分割为若干个相对独立的物理区域，对每一个区域分别进行本发明所述方法的检测，以提高流量参量的响应精度。

本说明书未作详细描述的内容属于本领域专业技术人员公知的现有技术。