基于稳态模型的气体管道泄漏定位方法的研究 .pdf

《基于稳态模型的气体管道泄漏定位方法的研究 .pdf》由会员分享，可在线阅读，更多相关《基于稳态模型的气体管道泄漏定位方法的研究 .pdf（6页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、第3 1 卷第 3 期 2 0 1 0年 3 月 仪 器 仪 表 学 报 C h i n e s e J o u rna l o f S c i e n t i fi c I n s t r u me n t V0 1 31 No 3 Ma r 2 01 0 基于稳态模型的气体管道泄漏定 位方法 的研究 术 孙良。 ，王建林 ( 1 北京化工大学信息科学与技术学院北京1 0 0 0 2 9 ； 2 中国石油大学信息与控制工程学院东营2 5 7 0 6 1 ) 摘要： 基于管道中气体流动的稳态模型， 根据对温度的处理方式， 分别提出气体管道泄漏定位的等温定位法和变步长龙格一 库塔法。通过实际气体

2、管道泄漏的实验验证对模型法的应用场合进行了分析， 并对两种方法的定位效果进行了对比。通过非 等温天然气管道泄漏的仿真试验 ， 对温差大小 、 气体特性估算的准确度和算法中模型的简化程度对定位效果的影响进行了研 究。仿真研究表明， 管道两端气体特性估算的准确程度决定了算法的定位精度， 而算法中的管道模型可以进行适度的简化。此 外 ， 变步长的龙格一 库塔法在非等温管道上的泄漏定位效果要优 于等温定位法 。 关键词： 气体管道； 泄漏定位 ； 模型； 龙格一 库塔法 中图分类号 ： T E 9 7 3 T P 2 7 4 文献标识码 ： A 国家标 准学科分类代码 ： 4 4 0 5 5 Re s

3、 e a r c h o n t h e m e t h o d f o r l e a k l o c a l i z a t i o n i n g a s p i p e l i ne ba s e d o n s t e a dy - s t a t e mo de l S u n L i a n g 一，Wa n g J i a n l i n ( J S c h o o l o f I n f o r m a t i o n S c i e n c e a n d T e c h n o l o g y ， B e q i n g U n i v e r s i t y of C

4、h e m i c a l T e c h n o l o g y ， B e q i n g 1 0 0 0 2 9 ，C h i n a ； 2 C o l l e g e of I n f o r m a t i o n a n d C o n t r o l E n g i n e e r i n g， C h i n a U n i v e r s i t y ofP e t r o l e u m， D o n g y i n g 2 5 7 0 6 1 ， C h i na ) Ab s t r a c t ：B a s e d o n t h e s t e a d y s t

5、 a t e mo d e l o f fl o w i n a g a s p i p e l i n e ， t w o me t h o d s f o r l e a k l o c a l i z a t i o n i n g as p i p e l i n e a r e p r o p o s e d，wh i c h u s e d i f f e r e n t t e mp e r a t u r e t r e a t me n t s r e s p e c t i v e l y On e me t h o d i s t h e i s o t h e r ma

6、 l l e a k l o c a t i o n me t h o d， wh i c h i s b a s e d o n t h e i s o t h e rm al a s s u mp t i o n；t h e o t h e r o n e i s v a ria b l e - - s t e p Ru n g e - Ku t t a me t h o d for n o n - i s o t h e rm al s i t u a _ t i o nTh r o u g h e x p e rime n t s c o n d u c t e d o n a r e

7、 al g a s p i p e l i n e，t h e a p p l i c a t i o n o c c a s i o n s o f t h e mo d e l b ase d l e a k l o c a t i o n me t h od s a r e a n aly z e d，a n d t h e p e r f o rm a n c e c o mp a ris o n b e t we e n t h e t wo me t h o d s i s als o ma d e Ac c o r d i n g t o t h e s i mula t i o

8、 n s o f l e a k a g e i n a n a t u r al g as p i pe l i n e，t h e e f f e c t s o f t e mp e r a t u r e d i f f e r e n c e，t h e a c c u r a c y o f e s t i ma t i o n i n g a s p r o p e r t i e s a n d t h e s i mp l i fi c a t i o n o f t h e mo d e l o n t h e a c c ura c y o f l e a k l o c

9、 a t i o n a r e s t u d i e dS i mu l a t i o n s r e s u l t s s h o w t h a t t h e e s t i ma t i o n a c c u r a c y i n g as p r o p rie t i e s a t e n d s d e t e r mi n e s t h e l e a k l o c a t i o n p r e c i s i o n o f t h e t wo me t h o d s ，wh i l e t h e mo d e l c a n mo d e r a

10、t e l y b e s i mp l i fi e d Co mp a r i s o n b e t we e n t h e t wo me t h o d s als o c o n fir ms t h a t t h e a c c u r a c y o f t h e v a r i a b l e s t e p Ru n g e Ku t t a me t h o d i s b e t t e r t h a n t h a t o f t h e i s o t h e rm al l e a k l o c a t i o n me t h o d for l e

11、ak l o c a t i o n i n n o n - i s o t h e r ma l p i pe l i n e Ke y wo r ds：g a s p i p e l i n e；l e a k l o c a t i o n；mo d e l ；Ru n g e Ku t t a me t h o d 1 引 言 管道的泄漏检测 与定位是管道安 全运行 的重 要保 收稿 日期 ： 2 0 0 9 - 0 7 R e c e i v e d D a t e ： 2 0 0 9 -0 7 基金项 目： 北京市 自然科学基金( N o 4 0 8 2 0 2 2 ) 资助项 目

12、 障， 目前有多种管道泄漏检测与定位 方法 ， 但 每种方 法都有其局限性 。其中基于模型的管道泄漏检测与定位 方法以管道输送 的稳态模 型或瞬态模型为基础， 实现对 管道沿线未知参数的预测， 在管道的泄漏检测 中具有较 学兔兔 w w w .x u e t u t u .c o m5 6 6 仪器仪表学报 第 3 1卷 大 的潜力 。 基于模型的泄漏定位方法分为基于瞬态模型和基于 稳态模型的泄漏定位方法。管道瞬态模型的求解只能采 用数值方法 ， 或者以这些数值方法为基础的状态空间 法 ， 计算误差大 ， 且计算 复杂 ， 耗时长 ， 难 以满足实时 定位的需要。基于稳态模 型的管道泄漏定位方

13、法计算简 单 ， 在液体管道 上的应 用 。 已经较为成熟。而气体管 道的流动过程受温度和气体可压缩性的影响较大 ， 相 比 液体管道来说气体管道的模型更为复杂。为提高计算速 度 ， 基于模型的气体管道 泄漏定位方法 ” 对于管道模 型不可避免地要进行各种简化 ， 但 目前还没有相关文献 对模型法的泄漏定位效果 、 模型 的简化对泄漏 定位 的影 响等进行研究和分析 。 本文根据气体管流的稳态模 型， 针对温度 的不 同处 理方式分别提 出气体管道泄漏定位的等温定位法和变步 长龙格一 库塔法 ， 可应用于等温和非等温气体管道 的泄漏 定位 ， 并通过实际管道的实验验证和管道的仿 真研究对 这两

14、种方法的应用场合 、 应用效果和模 型简化 的影响等 进行探讨 2 气体管道基本模型 忽略高程变化 ， 根据质量 、 动量和能量的守恒可以得 到气体 的稳态管流方程 ： 击 ( ) = 0 O x + = 一 D 2 ( 1 ) ( 2 ) p w A ( h+ 一 ) + ( T一 )：0 ( 3 ) 厶 式中： P是气体密度， k g m； 是气体流速 ， m s ； A 是水力摩阻系数； 是管道的轴向长度， 11 1 ； P为压力， P a ； D 是 管 道 内径 ， m； h为 焓 ， J- k g ； K 为 传 热 系 数 W 13 1 - K； 是气体 的温度 ， r， n 为

15、 土壤温度 ， K 。A 为管道截面积 ， m ， 假设在整个管长范 围内管道截 面积 保持不变。 气体管道模型的求解必须结合气体状态方程： P = f ( P ， T ) ( 4 ) 即可得到管道沿线压力 、 流量和温度 的分布。气体 状态方程的选取决定了模型求解的复杂程度。 3 基于模型的气体管道泄漏定位算法 3 1 等温定位法 该方法将管道中的气体流动过程视为等温过 程， 因 此可以不考虑能量方程 ， 对式 ( 1 ) 和( 2 ) 沿管长进行积分 求解 ， 可得到管 内稳态流动关系： ( JP ； 一 ) G Z R T 一 2 1 n ( P a P )=A L D ( 5 ) 式中

16、： 尺是气体常数 ， k J K ； A为管长范围内的 平均摩阻系数； z为管道 内气体的平均压缩系数； L为管 道总长 ， m； T是管道 内平均温度 ； P 为管 道人 口压力； 为管道 H 口的压力 ； G为 管道 内气体 的质 量 流速， G p w ， k g 1 1s。 当管道发生泄漏时 ， 可以将管道 以泄漏点为界分成 两个管段 ， 如图 1 所示。 图 1 泄漏管道示意 图 F i g I S c h e ma t i c r e p r e s e n t a t i o n s f o r l e a k i n g p i p e l i n e 由图 1 可见 ， 泄漏

17、点上游管段的质量流速为 G ， 泄 漏点处的压力 为 ， 管道 始端压力 为 ， 根据式 ( 5 ) 则有 ( P 2pP ) G Z R T一2 1 n P 0 P L=A L r D ( 6 ) 式中： 为泄漏点到管道始端的距离。泄漏点下游管段 的质量流速为 ， 管道末端压力为 ， 根据式( 5 ) 则有 ( P P ； ) G z Z R T一2 1 n P P =A ( 一 r ) D ( 7 ) 综合式( 6 ) 和( 7 ) ， 可以得到 一 ! 竺 ： 塑! 二 ： 塑! L 1 G z 2 Z RT一1 c a2 Z RT ， ( P ；一 P ) G q Z R T一 2 1

18、 n ( P Q P L ) L r = 一 A D ( 8 ) ( 9 ) 式中： 平均摩阻系数 A需要在管道无泄漏时通过测得管 道始末端压力和流量 ， 根据式 ( 5 ) 进行辨识得到。 文献 1 1 、 1 2 和 1 3 所提 出的泄漏定位方法都属 于等温定位法 。其 中文献 1 1 和 1 2 所提 出算法忽略 式 ( 6 ) 和( 7 ) 中的 2 1 n P 。 P 和 2 1 n P， P 项 ， 这两项在泄 漏量较大时会引入计算误差。而文献 1 3 为便于采用 状态空间的形式在文献 1 1 和 1 2 的基础上对流量进 行线性化 ， 从而引入更大的计算误差。 等温定位法根据泄

19、漏后管道两端 的压力和流量即可 确定泄漏位置， 计算简单， 实施方便。但对温度的处理上 采用平均温度进行计算 ， 对于非等温管道会引入一定 的 计算误差 。 3 2 变步长龙格- 库塔法 对于非等温 的流动过程 ， 由于式 ( 1 ) 、 ( 2 ) 和( 3 ) 的 求解没有解析表达式 ， 一般多采用精度较高的四阶龙格一 学兔兔 w w w .x u e t u t u .c o m第 3 期 孙 良等 ： 基于稳态模型的气体管道泄漏定位方法的研究 5 6 7 库塔法 ， 求 解 时需要将 管道分 为 段 ， 每段 管长 为 ，根据管道人 口条件下的温度 、 压力 和流速等参数 以 为步长可

20、求得管道从始端到末端各截面的参数。 通过对 4 1 节中的气体管道进行泄漏仿真 ， 泄漏前 后管道沿线压力 分 布见图 2所 示 ， 其 中曲线 1和 曲线 A O B分别为泄漏前后管道压力梯度 曲线。由图 2可见 ， 无泄漏时 ， 管道 的压力梯度 曲线是 光滑的 ， 当泄漏发 生 后 ， 管道的压力梯度曲线 在泄漏点处 出现拐点。分别 按 照管道始端和末端条件下 的温度 、 压力和流速 等参数求 解正压力梯度 曲线 A A 和逆压力梯 度曲线 B B ， 两者 的 交点 O所对应的位置即为泄漏点。其中 A A 的求解仍按 式( 1 ) 、 ( 2 ) 和( 3 ) 进行， 而 B B 的求

21、解是从管道末端向管 道始端进行 ， 因此需要对摩擦阻力项 和热损失项 的方 向 性加 以调整 ， 即式( 2 ) 调整为 ： 箬+ = 去 ( 1 0 ) 式( 3 ) 调整为： 2 p w A ( h+ ) 一K r r D( T一7 “0 )=0 ( 1 1 ) 其余各项保持不变 。 管道长度 图 2 管道泄漏定位示 意图 F i g 2 Sc h e ma t i c r e p r e s e nt a t i o ns f o r l e a k l o c a t i o n 根据定位的原理可知， 计算管段 A x 的大小决定 了定 位精度。 越小 ， 定位精度越高， 但计算 耗时

22、越长 ， 难 以 完成实时泄漏定位的需要 ； A x 越大 ， 计算速度越快 ， 但定 位精度越低。为此本 文提 出变步 长的龙 格一 库塔法进 行 泄漏定位， 即首先选取大的计算步长， 在整个管长范围内 进行泄漏位置的初选 ， 然后逐步缩小泄漏定位 的搜索 范 围并逐步递减计算步长 ， 直到最终满足泄漏定位精度 的 需要。计算流程如图 3 所示。 图 3中 l c 为设定 的定位精度。 c为计 算步长递减 的比例 ， k c 1 。初始步长和 Ji c的选取影 响计算精度 和 计算速度。初始步长的选取应 首先保证龙格一 库塔法 的 计算不发散， 其次由于气体管道压力分布 的非线性， 为保 证

23、计算精度， 步长不能选取过大。可事先逐步调整步长 的大小 ， 以 和 2分别作 为步长对模 型进 行计 算 ， 若两者的计算结果相差在误差允许范围内， 可将此时 的 作为初始步长。k 选取过大和过小都会影响计算 的速度 ， 本文取 c为 0 5 。 变步长龙格一 库塔法的计算量要比等温定位法大 ， 但 可用于非等温气体管道的泄漏定位 。 图3 变步长龙格一 库塔法计算流程 Fi g 3 Pr o c e d ur e o f v a ria b l e s t e p Ru n g e Ku t t a me t h o d 4 实验验证 4 1 气体管道泄漏定位的实验验证 对中国石油大学 (

24、 华 东) 储 运实验室的一条气体管 道进行 不 同 泄漏 点 和 泄 漏 量 的 现场 试 验 ， 管 道 总 长 2 5 2 4 m， 管径 2 5 m m， 介质为空气 。分别在距离管道始端 3 8 5 m和 2 1 1 m处设置泄漏点 ， 和 ， 并在各泄漏点安 装手动阀控制泄漏量的大小 ， 管道的始端和末端分别安装 有压力传感器和流量计。实验设备流程图如图 4 所示。 图 4 实验设备流程 图 F i g 4 S c h e ma t i c d i a g r a m o f e x p e ri me n t f a c i l i t y 学兔兔 w w w .x u e t

25、u t u .c o m5 6 8 仪器仪表学报 第 3 1 卷 图4中 P 。、 Q m q为管道始端 的压力 和流量测量 ， P 、 Q m z 为管道末端的压力和流量测量 。 图5为泄漏点 厶 的一次泄漏测试所测到的管道两 端参数变化以及泄漏定位曲线 ， 泄漏率为 1 1 。 、 趟 、 R 幽 时间 s ( a ) 泄漏前后管道两端压力变化曲线 ( a ) V a ri a t i o n o f t h e p r e s s u r e a l b o t h e n d s b e f o r e a n d a f t e r l e a k a g e 时 间 ， s ( b

26、 ) 泄漏前后管道两端流量变化皓线 ( b ) V a ri a t i o n of t h e fl o w r a t e a t b o t h e n d s b e f o r e a n d a f t e r l e a k a g e 时 间 ， ( C ) 泄漏定位曲线 ( C ) E s t i m a t i o n o f l e a k l o c a t i o n 图 5 泄漏前后管道两端参数变化及泄漏定位 F i g 5 Va r i a t i o n o f p a r a me t e r s a t b o t h e n d s of t h e p

27、 i p e l i n e b e f o r e a n d a f t e r l e a k a g e a n d e s t i ma t i o n of l e a k l o c a t i o n 由图 5的( a ) 和( b ) 可见 ， 管道发生泄漏 时， 管道两 端参数的变化经历 一个稳 态一 快 瞬态一 慢 瞬态一 稳态 的过 程。其 中的快瞬态过程是流动参数的变化在泄漏点和管 道两端之间来回传播反射 的过程 ， 由于泄漏点 的工况未 知， 因此无论瞬态模型还是稳态模型都无法对这个过程 进行描述， 因此基于模型 的泄漏定位 方法在这个过程是 无法进行定位的。只有进

28、入慢瞬态过程时， 泄漏点的工 况基本稳定 ， 基于模型的泄漏定位算法才能稳定定位 ， 见 图 5 ( C ) 。由图 5可见 ， 在 2 2 S 到 3 0 S的时间内， 由于泄 漏导致管道两端参数变化剧烈 ， 而 3 0 S以后变化趋缓， 泄漏定位算法才能够输出稳定的定位结果 。 根据在泄漏点 。 和 进行的泄漏测试 ， 分别按等温 定位法和变步长龙格一 库塔法计算泄漏位置 ， 计算值与实 际泄漏位置的对 比见表 1 。 表 1 泄漏定位结果与定位误差 Tab l e l Le a k l oc at i o n r e s ul t s a nd I o c a fiO n e r r o

29、r s 由表 l 可知， 两种定位算法都在小泄漏时定位误差 较大 ， 而大泄漏时定位误差较小。对该等温管道来说 ， 两 种算法都可以较为准确的进行泄漏定位 ， 且定位误差基 本相当。但变步长 龙格 库塔法计算量要 比等稳定位法 大 ， 因此等温定 位法更 适用 于等温管道 或温差 不 大的 管道。 4 2 气体管道泄漏定位的仿真验证 对于非等温的天然气管道， 管道工况较为复杂， 一是 管道沿程温度的变化 ， 管道两端存在温差 ； 二是天然气的 组分复杂 ， 因此对气体的密度 、 压缩系数等气体特性进行 估算时往往不可避免地存在误差。采用仿真的气体管道 进行泄漏的模拟 ， 可以充分考察介质特性的

30、估算 、 管道两 端温差大小以及模型的简化程度对两种算法的定位精度 的影响。 本文采用的管道为普光气 田的一段外输管道， 管道 口径为 6 3 2 m m， 长度为 1 0 0 k m， 介质为天然气 ， 各组分所 占百分比见表 2所示 。正 常工况 下， 管道 始端 压力为 9 M P a ， 末端压力为 8 0 6 M P a ， 流量为 4 1 4 1 0 N m h 。 管道始端温度为 3 1 0 K， 土壤温度为 2 8 0 K 。分别在距离 管道始端 1 0 k m， 5 0 k m和 9 0 k m处按分支管道模拟泄 漏 ， 泄漏率分别为 0 5 ， 1 0 和 5 。 学兔兔

31、w w w .x u e t u t u .c o m第 3期 孙 良等 ： 基 于稳态模型的气体管道泄漏定位方法的研究 表 2 仿真管道 中天然气 组分 Ta b l e 2 Co mpo s i t i on o f na t u r a l g a s i n t h e s i mul a ting pi pe l i n e CH4 C2 H6 C3 H8 CO2 H2 S N2 H2 HE 7 6 5 2 0 0 1 2 0 0 0 0 8 8 6 3 0 1 4 1 4 0 0 5 5 2 0 0 2 0 0 01 0 气体管道的仿真模型采用 B WR S 状态方程 ， 可以比

32、较准确地描述气体 的特性 。但若在变步长龙格一 库塔法 中仍采用该状态方程 ， 需要在每个管段都进行气体特性 的计算 ， 会使得管道模型求解复杂 ， 计 算量 大， 则会 因计 算耗时多而难 以满足实时泄漏定位的需要， 因此在定位 算法中必须对管道模型进行简化以加快计算速度。本文 通过 3个方案对管道两端气体特性估算的准确度和算法 中模型简化对泄漏定位 的影响进行研究 ： 方案 1 ： 采用理想的空气近似实际 的天然气 ， 变步长 龙格一 库塔法中管道模 型计算采用理想气体的状态方程 ； 等温定位法的平均压缩系数取为 1 。 方案 2 ： 对 管道两端气体质量流量 的估算采用带压 缩 因子的实

33、际气体的状态方程 ， 即 P = ， 其中 z的 计算采用前苏联气体研究所针对石油伴生气的压缩系数 计算公式 ， 即： z = 1 0 0 0 1 + 5 P ( 1 2 ) 式 中： P为气体压力 ， l O P a 。变步长龙格 库塔法中管道 模型计算引入始末两端平均压缩系数。等温定位法 中的 平均压缩系数也同样取值。 方案 3 ： 对管道两端气体质量流量的估算采用 B WR S 状态方程 ， 变步长龙格一 库塔法 中管道模型的计算仍引人 始末两端的平均压缩系数 ， 只是压缩系数的计算也采用 B WR S状态方程 ； 等温定位法中的平均压缩 系数也 同样 取值。 各方案 中的等温定位法均以

34、管道始末端的平均温度 作为管道温度。三种方案的定位结果对 比见表 3所示 。 表 3 三个方案的泄漏定位结果对 比 Tabl e 3 Compa r i s o n o f l e ak l o c a t i o n a mon g t hr e e s c he me s 由表 3可见 ， 各方案在小泄漏率时， 定位误差大， 随 着泄漏率的增大， 定位误差越来越小。三种方 案分别对 应不同精度的介质特性的估算 ， 随着对 管道两端介质特 性估算精度 的提高 ， 定位误差也越来 越小。其 中方 案 3 中两种定位方法的定位精度都较高 ， 因此可以说对管道 两端气体特性估算越准确 ， 定位精度

35、越高 ， 而算法中的管 道模型只需引入平均压缩系数进行补偿即可， 无需在每 个管段都进行复杂的气体特性计算。从而大大简化了管 道模型 ， 提高了变步长龙格一 库塔法的计算速度 。此外 由 表 3 可见 ， 各方案中变步长龙格一 库塔法 的定位误差要小 于等温定位算法。 为考察管道两端温差大小 对算法定位精度 的影 响， 本文通过仿真 ， 提高管道入 口温度 ， 维持土壤 温度不变 ， 实现管道两端温差 的增大。针对 1 0 k m的泄漏位置 ， 取 泄漏率为 1 ， 两种算法均采用方案 3 ， 分别 考察温差大 小对两种算法定位精度的影响 ， 见表 4所示。 表 4 两端温差变化对 定位结果的

36、影响 Ta bl e 4 Ef f k t o f t e mpe r a t ur e di ffe r e nc e be t we e n bo t h e n ds o n l e ak l o c a tio n r e s ult 始端温度 K 始末端温差 K等温定位法 k m 变步长龙一 库法 k m 由表 4可见 ， 随着始末两端温差的增大 ， 两种算法的 定位误差均增加 ， 但变步长龙格一 库塔法的定位误差变化 要远小于等温定位法。因此变步长龙格 库塔法更适用 于非等温气体管道的泄漏定位 。 学兔兔 w w w .x u e t u t u .c o m5 7 0 仪器仪表学

37、报 第 3 1卷 5 结 论 本文基于管道 内气体流动的稳态模 型， 针对温度的 不同处理方式 ， 分别提 出了管道泄漏定位 的等温定位法 和变步长龙格一 库塔法。 通过对实际气体管道泄漏定位的实验分析可知基于 模 型的泄漏定位法可以在泄漏发生后的慢瞬态过程输出 稳定的定位结果。对于等温管道， 两种定位方法的定位 精度基本相当， 但等温定位法计算量小 ， 更适合于等温管 道或两端温差不大的管道 的泄漏定位 。 通过对天然气管道泄漏的仿真研究表明， 对于非等温管 道， 管道两端介质特性估算的准确性对泄漏定位精度影响较 大， 而算法中的管道模型可以引入管道两端平均压缩因子进行 补偿即可达到较好的定

38、位效果。从而简化了模型的计算， 加快 了算法的计算速度。此外， 变步长龙格一 库塔法在非等温气体 管道上的定位精度要略高于等温定位法， 且受管道两端温差变 化的影响较小， 因此更适用于非等温气体管道的泄漏定位。 参考文献 1 刘恩斌 ， 彭善碧 ， 等基于瞬态模 型的油气 管道泄漏检 测 J 天然气工业 ， 2 0 0 5 ， 2 5 ( 6 ) ： 1 0 2 1 0 4 LI U E B，PENG SH B，e t a 1 Le a ka g e d e t e c t i o n o f o i l a nd g a s p i p e l i n e s b y t r a n s i

39、 e n t m o d e l J N a t u r a l G a s I n d u s t r y ， 2 0 0 5 ， 2 5 ( 6 ) ：1 0 2 1 0 4 J 2 J A N D R E W F C ， P E D R O L ， B R Y A N W K A s e l e c t i v el i t e r a t u r e r e v i e w o f t r a n s i e n t b a s e d l e a k d e t e c t i o n m e t h ods J J o m n al o f H y d r o E n v i r

40、o n m e n t R e s e a r c h ， 2 0 O 9 ， 4 ( 2 ) ： 2 1 2 2 2 7 3 G E C H H， Y A N G H Y ， Y E H， e t a 1 A f a s t l e a k l o c a ri n g m e t h o d b a s e d o n w a v e l e t t r a n s f o r m J T s i n g h u a S c i e n c eT e c h n o l o g y ， 2 0 0 9 ，1 4 ( 5 ) ： 5 5 1 5 5 5 4 Y A N G J ，WE N Y

41、 M，L I P L e a k l o c a t i o n u s i n g b l i n d s y s t e m i d e n t i f i c a t i o n i n w a t e r d i s t r i b u t i o n p i p e l i n e s J J o r u n a l o f S o u n d a n d V i b r a t i o n 2 0 0 83 1 0：1 3 4 1 4 8 5 刘恩斌 ， 彭善 碧 基于瞬态模型的油气管道泄漏检测 J 天然气工业 ， 2 0 0 5 ， 2 5 ( 6 ) ： 1 0 2 1 0 4

42、 L I U EN B P E NG S H B L e a k a g e d e t e c t i o n o f o i l a n d g a s p i p e l i n e s b y t r a n s i en t mo d e 1 Na t u r a l Ga s I n d u s t r y J 2 0 0 5 ， 2 5 ( 6 ) ： 1 0 2 1 0 4 6P H Y S I K E R R T M o d e l b a s e d p i p e l i n e l e a k d e t e c t i o n a n d l o e a l i z

43、a t i o n J 3 R I n t e r n a t i o n a l ， 2 o 0 3 ， 4 2 ( 7 ) ： 4 5 5 - 4 6 0 7 B I L L MA N N L ，I S E R M A N N R L e a k d e t e c t i o n m e t h o d s f o r p i p e l i n e s J A u t o m a t i e a ， 1 9 8 7 ， 2 3 ( 3 ) ： 3 8 1 3 8 5 I 8 I B E N K H E R O U F A，A L L I D I N A Y A Lea k d e t

44、e c t i o n a n d l o c a t i o n i n g a s p i p e l i n e s J I E E P r o c e e d i n g s ，1 9 8 8 ， 1 3 5 ( 2 ) ： 1 4 2 1 4 8 9 Z H A N G X J S t a t i s t i c a l l e a k d e t e c t i o n i n g a s a n d l i q u i d p i p e l i n e s J P i p e s&P i p e l i n e s i n t e r n a t i o n a l ，1 9

45、9 3 ， ( 8) ：2 6 2 9 1 1 0 T u Y Q，Q I U D Q， Z H A N G D A r e v i s e d s t a t i c a l d e d u c t i o n a p p r o a c h t o p i p e l i n e l e a k d e t e c t i o n 2 0 1 E E E I ns t r ume n t a t i o n a n d Me a s u r e me n t Te C hn o l o g y Co n f e r - e n e e ， v a i l ，C o l o r a d o

46、， U S A， 2 0 0 3 ( 2 ) ：1 6 3 9 1 6 4 2 1 1 孟令稚基于瞬态模型法的输气管道泄漏监测与定位 技术 J 北京交通大学学报， 2 0 0 8 ， 3 2 ( 3 ) ： 7 3 - 7 7 MENG L YRe s e a r c h o n l e a k de t e c t i o n a nd po s i t i o n f o r n a t ur a l g a s p i p e l i ne b a s e d o n t r a ns i e n t mo d e l me t h o d J J o u r n a l o f B e i j i n g J i a o T o n g U n i v e r s i t y ，2 0 0 8 ，3 2 ( 3 )：7 3 - 7