基于稳态模型的气体管道泄漏定位方法的研究 .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 气体管道泄漏定位的实验验证 对中国石油大学 (
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 稳态 模型 气体 管道 泄漏 定位 方法 研究
限制150内