发明名称 基于广义模糊双曲正切模型的输油管网泄漏检测方法
摘要 一种基于广义模糊双曲正切模型的输油管网泄漏检测方法,属于管道检测技术领域。按如下步骤:1、拾取到负压波信号,传递给信号调理板,计算负压波产生的起始位置;2、利用广义模糊双曲正切模型对负压波来源进行分类,判断负压波的产生是由于泄漏、压力越站还是工况调整;3、若负压波来源于压力越站,则检测前一段管道上有没有泄漏,回到步骤1;若负压波来源于泄漏则执行步骤4;若负压波来源于工况调整则执行步骤5;4、给出泄漏警报;本发明的优点:采用广义模糊双曲正切模型来区分负压波的来源。将阀门开度,泵的状态,流量,温度,压力,密度作为广义模糊双曲正切模型的输入变量,并通过输出值来判断是否发生泄漏,防止误报警。
申请公布号 CN102242872B 申请公布日期 2013.01.30
申请号 CN201110169827.1 申请日期 2011.06.22
申请人 东北大学 发明人 冯健;刘金海;张化光;马大中;魏向向;李健
分类号 F17D5/06(2006.01)I 主分类号 F17D5/06(2006.01)I
代理机构 沈阳东大专利代理有限公司 21109 代理人 朱光林
主权项 1.一种基于广义模糊双曲正切模型的输油管网泄漏检测方法,其特征在于,按如下步骤:步骤1、安装在管网中一段管道两端的数据采集装置中的压力传感器拾取到负压波信号,传递给信号调理板,根据负压波到达两端的时间差计算负压波产生的起始位置;根据负压波到达单条管道两端的时间差计算负压波产生的起始位置;公式如下,<maths num="0001"><![CDATA[<math><mrow><mi>X</mi><mo>=</mo><mfrac><mrow><mi>L</mi><mo>-</mo><mi>&alpha;</mi><mrow><mo>(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>-</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow></mrow><mn>2</mn></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中X-输油管道泄漏点首端的距离,m;L-输油管道长度,m;α-负压波在输油管道中的传播速度,m/s;t<sub>1</sub>-负压波到达首端时间,s;t<sub>2</sub>-负压波到达末端时间,s;管道中存在柴油、汽油混输的情况,此时波速在两种液体中的速度是不同的,在这种情况下检测到负压波存在时计算负压波产生的起始位置公式转换为:<maths num="0002"><![CDATA[<math><mrow><mi>X</mi><mo>=</mo><mfrac><mrow><mrow><mo>(</mo><msub><mi>t</mi><mn>1</mn></msub><mo>-</mo><msub><mi>t</mi><mn>2</mn></msub><mo>-</mo><mfrac><mi>S</mi><msub><mi>&alpha;</mi><mn>1</mn></msub></mfrac><mo>)</mo></mrow><msub><mi>&alpha;</mi><mn>2</mn></msub><mo>+</mo><mi>S</mi><mo>+</mo><mi>L</mi></mrow><mn>2</mn></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中t<sub>1</sub>-负压波到达首端时间,s;t<sub>2</sub>-负压波到达末端时间,s;S-输油管道中汽油长度,m;α<sub>1</sub>-负压波在汽油中传播速度,m;α<sub>2</sub>-负压波在柴油中传播速度,m;L-输油管道长度,m;X-输油管道泄漏点距首端的距离,m;步骤2、某一管段上负压波是由本段上某点泄漏、其他管段泄漏产生的压力越站或工况调整产生,利用广义模糊双曲正切模型对负压波来源进行分类,判断负压波的产生是由于泄漏、压力越站还是工况调整;所述的利用广义模糊双曲正切模型对负压波来源进行分类方法如下:1)、模型输出值为1时,代表泄漏发生,负压波来源于泄漏;2)、模型输出值为0.5时,代表负压波来源于压力越站;3)、模型输出值为0时,代表负压波来源于工况调整;广义模糊双曲正切模型第l条模糊规则的形式为:<maths num="0003"><![CDATA[<math><mrow><msup><mi>R</mi><mi>l</mi></msup><mo>:</mo><mi>IF</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>-</mo><msub><mi>d</mi><mn>11</mn></msub><mo>)</mo></mrow><mi>is</mi><msub><mi>F</mi><msub><mi>x</mi><mn>11</mn></msub></msub><mi>and</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>-</mo><msub><mi>d</mi><mn>12</mn></msub><mo>)</mo></mrow><mi>is</mi><msub><mi>F</mi><msub><mi>x</mi><mn>12</mn></msub></msub><mi>and</mi><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mi>and</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>-</mo><msub><mi>d</mi><mrow><mn>1</mn><msub><mi>w</mi><mn>1</mn></msub></mrow></msub><mo>)</mo></mrow><mi>is</mi><msub><mi>F</mi><msub><mi>x</mi><mrow><mn>1</mn><msub><mi>w</mi><mn>1</mn></msub></mrow></msub></msub><mi>and</mi></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><mrow><mo>(</mo><msub><mi>x</mi><mn>2</mn></msub><mo>-</mo><msub><mi>d</mi><mn>21</mn></msub><mo>)</mo></mrow><mi>is</mi><msub><mi>F</mi><msub><mi>x</mi><mn>21</mn></msub></msub><mi>and</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>2</mn></msub><mo>-</mo><msub><mi>d</mi><mn>22</mn></msub><mo>)</mo></mrow><mi>is</mi><msub><mi>F</mi><msub><mi>x</mi><mn>22</mn></msub></msub><mi>and</mi><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mi>and</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>2</mn></msub><mo>-</mo><msub><mi>d</mi><mrow><mn>2</mn><msub><mi>w</mi><mn>2</mn></msub></mrow></msub><mo>)</mo></mrow><mi>is</mi><msub><mi>F</mi><msub><mi>x</mi><mrow><mn>2</mn><msub><mi>w</mi><mn>2</mn></msub></mrow></msub></msub><mi>and</mi><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><mi>and</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>n</mi></msub><mo>-</mo><msub><mi>d</mi><mrow><mi>n</mi><mn>1</mn></mrow></msub><mo>)</mo></mrow><mi>is</mi><msub><mi>F</mi><msub><mi>x</mi><mrow><mi>n</mi><mn>1</mn></mrow></msub><mtext></mtext></msub><mi>and</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>n</mi></msub><mo>-</mo><msub><mi>d</mi><mrow><mi>n</mi><mn>2</mn></mrow></msub><mo>)</mo></mrow><mi>is</mi><msub><mi>F</mi><msub><mi>x</mi><mrow><mi>n</mi><mn>2</mn></mrow></msub></msub><mi>and</mi><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mi>and</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>n</mi></msub><mo>-</mo><msub><mi>d</mi><mrow><mi>n</mi><msub><mi>w</mi><mi>n</mi></msub></mrow></msub><mo>)</mo></mrow><mi>is</mi><msub><mi>F</mi><msub><mi>x</mi><mrow><mi>n</mi><msub><mi>w</mi><mi>n</mi></msub></mrow></msub></msub></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><mi>THEN</mi><msup><mi>y</mi><mi>l</mi></msup><mo>=</mo><msub><mi>c</mi><msub><mi>F</mi><mn>11</mn></msub></msub><mo>+</mo><msub><mi>c</mi><msub><mi>F</mi><mn>12</mn></msub></msub><mo>+</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><msub><mi>c</mi><msub><mi>F</mi><mrow><mn>1</mn><msub><mi>w</mi><mn>1</mn></msub></mrow></msub></msub><mo>+</mo><msub><mi>c</mi><msub><mi>F</mi><mn>21</mn></msub></msub><mo>+</mo><msub><mi>c</mi><msub><mi>F</mi><mn>22</mn></msub></msub><mo>+</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><msub><mi>c</mi><msub><mi>F</mi><mrow><mn>2</mn><msub><mi>w</mi><mn>2</mn></msub></mrow></msub></msub><mo>+</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><mo>+</mo><msub><mi>c</mi><msub><mi>F</mi><mrow><mi>n</mi><mn>1</mn></mrow></msub></msub><mo>+</mo><msub><mi>c</mi><msub><mi>F</mi><mrow><mi>n</mi><mn>2</mn></mrow></msub></msub><mo>+</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><msub><mi>c</mi><msub><mi>F</mi><mrow><mi>n</mi><msub><mi>w</mi><mi>n</mi></msub></mrow></msub></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,w<sub>i</sub>为将x<sub>i</sub>线性变换的个数,i=1,…,n;d<sub>ij</sub>为x<sub>i</sub>线性变换点,i=1,…,n;j=1,…w<sub>i</sub>;<img file="FDA00002053154000026.GIF" wi="76" he="65" />为与<img file="FDA00002053154000027.GIF" wi="156" he="61" />对应的模糊子集,包括正P和负N两个语言值,当<img file="FDA00002053154000028.GIF" wi="76" he="65" />为正P时,<img file="FDA00002053154000029.GIF" wi="75" he="65" />为<img file="FDA000020531540000210.GIF" wi="98" he="65" />当<img file="FDA000020531540000211.GIF" wi="76" he="65" />为负N时,<img file="FDA000020531540000212.GIF" wi="76" he="65" />为<img file="FDA000020531540000213.GIF" wi="229" he="65" />是与<img file="FDA000020531540000214.GIF" wi="75" he="64" />对应的输出常数,IF中输入变量以及THEN中输出常数项都是可选的,但是输出项<img file="FDA000020531540000215.GIF" wi="70" he="53" />与输入变量是一一对应的,即如果在IF部分包括<img file="FDA000020531540000216.GIF" wi="106" he="65" />则在THEN部分应包括<img file="FDA000020531540000217.GIF" wi="71" he="54" />项;相反,如果IF部分不包括<img file="FDA000020531540000218.GIF" wi="108" he="65" />则在THEN部分也不包括<img file="FDA000020531540000219.GIF" wi="71" he="55" />项;给定一组广义双曲正切型模糊规则基,定义广义输入变量x<sub>i</sub>=x<sub>z</sub>-d<sub>zj</sub>           (4)d<sub>zj</sub>为x<sub>z</sub>线性变换点,j=1,......,w<sub>z</sub>;设<img file="FDA000020531540000220.GIF" wi="232" he="113" />式中w<sub>z</sub>为将x<sub>z</sub>线性变换的个数,其中z=1,......,n,取广义输入变量对应的模糊集合<img file="FDA000020531540000221.GIF" wi="48" he="56" />和<img file="FDA000020531540000222.GIF" wi="61" he="56" />的隶属函数为<img file="FDA000020531540000223.GIF" wi="146" he="62" />和<img file="FDA000020531540000224.GIF" wi="178" he="62" /><maths num="0008"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>&mu;</mi><msub><mi>P</mi><msub><mi>x</mi><mi>i</mi></msub></msub></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><msup><mi>e</mi><mrow><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>k</mi><msub><mi>x</mi><mi>i</mi></msub></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msup></mtd></mtr><mtr><mtd><msub><mi>&mu;</mi><msub><mi>N</mi><msub><mi>x</mi><mi>i</mi></msub></msub></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><msup><mi>e</mi><mrow><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>+</mo><msub><mi>k</mi><msub><mi>x</mi><mi>i</mi></msub></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msup></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,<img file="FDA000020531540000226.GIF" wi="48" he="56" />为常数,将<img file="FDA000020531540000227.GIF" wi="155" he="62" />简写为<img file="FDA000020531540000228.GIF" wi="281" he="62" />简写为<img file="FDA000020531540000229.GIF" wi="182" he="56" />简写为k<sub>i</sub>,根据模糊规则基得出如广义模糊双曲正切模型:<maths num="0009"><![CDATA[<math><mrow><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mfrac><mrow><msub><mi>c</mi><msub><mi>P</mi><mi>i</mi></msub></msub><msup><mi>e</mi><mrow><msub><mi>k</mi><mi>i</mi></msub><msub><mi>x</mi><mi>i</mi></msub></mrow></msup><mo>+</mo><msub><mi>c</mi><msub><mi>N</mi><mi>i</mi></msub></msub><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>k</mi><mi>i</mi></msub><msub><mi>x</mi><mi>i</mi></msub></mrow></msup></mrow><mrow><msup><mi>e</mi><mrow><msub><mi>k</mi><mi>i</mi></msub><msub><mi>x</mi><mi>i</mi></msub></mrow></msup><mo>+</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>k</mi><mi>i</mi></msub><msub><mi>x</mi><mi>i</mi></msub></mrow></msup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,<img file="FDA00002053154000031.GIF" wi="198" he="113" />为合计的规则数目,<img file="FDA00002053154000032.GIF" wi="44" he="45" />是与<img file="FDA00002053154000033.GIF" wi="48" he="56" />对应的输出常量,<img file="FDA00002053154000034.GIF" wi="51" he="45" />与<img file="FDA00002053154000035.GIF" wi="61" he="56" />对应的输出常量;x<sub>1</sub>,…,x<sub>n</sub>为从管道中采集来的数据,压力、流量、温度、阀门开度、泵的状态等的改变量,f(x)为模型输出,根据f(x)判断管网中负压波来源;步骤3、若负压波来源于压力越站,则检测前一段管道上有没有泄漏,回到步骤1;若负压波来源于泄漏则执行步骤4;若负压波来源于工况调整则执行步骤5;步骤4、给出泄漏警报;步骤5、结束。
地址 110819 辽宁省沈阳市和平区文化路3号巷11号