專利名稱：燃?xì)夤艿佬孤z測和定位方法
技術(shù)領(lǐng)域：
本發(fā)明屬于氣體檢測和監(jiān)控技術(shù)，具體涉及一種基于地理信息和數(shù)據(jù)采集與監(jiān)控系 統(tǒng)技術(shù)的對燃?xì)夤芫W(wǎng)泄漏自動(dòng)檢測和定位的方法。
背景技術(shù)：
隨著城市燃?xì)鈶?yīng)用范圍的擴(kuò)大，天然氣管線系統(tǒng)的鋪設(shè)也趨于密集化。眾所周知， 天然氣是一種易燃易爆氣體，其儲(chǔ)存和輸送的安全性要求非常高。尤其對于長度在幾百 公里以上的燃?xì)廨斶\(yùn)管道，大部分都要經(jīng)過邊遠(yuǎn)及條件復(fù)雜的地區(qū)，這些地區(qū)除人員稀 少外，燃?xì)夤艿来蠖际请[蔽工程，點(diǎn)多、線長、面廣，發(fā)生泄漏事故之后難以及時(shí)發(fā)現(xiàn) 或確定泄漏地點(diǎn)，可能釀成較大的事故。而一旦發(fā)生事故將會(huì)造成巨大的生命財(cái)產(chǎn)損失 和環(huán)境污染，這是任何國家都必須面對的安全問題。如果能夠及時(shí)發(fā)現(xiàn)官網(wǎng)泄漏并確定 具體的泄漏地點(diǎn)，就能有效地避免或減輕泄漏事故所造成的危害，因此，對長距離燃?xì)?管網(wǎng)的泄漏自動(dòng)進(jìn)行檢測和定位，其經(jīng)濟(jì)和社會(huì)效益是不言而喻的。目前國內(nèi)的城市天 然氣管道仍然采用技術(shù)較為落后的人工巡檢檢漏方式，即巡檢員定期攜帶相關(guān)的檢漏儀 器沿管線進(jìn)行巡視，這樣有可能要等到泄漏發(fā)生很長時(shí)間后才能被發(fā)現(xiàn)，而對于一些地 形復(fù)雜的情況甚至檢測不到泄漏。
SCADA系統(tǒng)(Supervisory Control And Data Acquisition)系統(tǒng)，即數(shù)據(jù)采集與監(jiān)視控 制系統(tǒng)，是以計(jì)算機(jī)為基礎(chǔ)的生產(chǎn)過程控制與調(diào)度自動(dòng)化系統(tǒng)。它可以對現(xiàn)場的運(yùn)行設(shè) 備進(jìn)行監(jiān)視和控制，以實(shí)現(xiàn)數(shù)據(jù)采集、設(shè)備控制、測量、參數(shù)調(diào)節(jié)以及各類信號報(bào)警等 各項(xiàng)功能。由于各個(gè)應(yīng)用領(lǐng)域?qū)CADA的要求不同，所以不同應(yīng)用領(lǐng)域的SCADA系統(tǒng)發(fā) 展也不完全相同。
地理信息系統(tǒng)(GIS， geographic information system)是隨著地理科學(xué)、計(jì)算機(jī)技術(shù)、 遙感技術(shù)和信息科學(xué)的發(fā)展而發(fā)展起來的一個(gè)學(xué)科。它能夠把圖形管理系統(tǒng)和數(shù)據(jù)管理 系統(tǒng)有機(jī)的結(jié)合起來，從而克服了數(shù)據(jù)庫和圖形系統(tǒng)各自固有的局限性，使二者的優(yōu)勢 互補(bǔ)，功能更加齊全。
發(fā)明內(nèi)容
本發(fā)明的目的是，提供一種對長距離燃?xì)廨斶\(yùn)管線實(shí)時(shí)進(jìn)行泄漏監(jiān)測及泄漏點(diǎn)精確 定位的方法。數(shù)據(jù)采集與監(jiān)視控制系統(tǒng)(Supervisory Control And Data Acquisition)以下簡稱 SCADA;地理信息系統(tǒng)(Geographic Information System)以下簡稱GIS;全球定位系統(tǒng)-以下簡稱GPS。
燃?xì)夤艿佬孤z測和定位方法的原理是，采用壓力傳感器用來采集在時(shí)間序列上所 對應(yīng)的壓力信號，或者說是各個(gè)時(shí)刻的壓力值；采用GPS是強(qiáng)化采集的各壓力信號與 時(shí)間脈沖信號的同步關(guān)系；采用A/D模數(shù)轉(zhuǎn)換卡是將壓力傳感器接受到的模擬信號轉(zhuǎn)換 成數(shù)字信號。具體方法為在一段燃?xì)夤艿赖膬啥嗽O(shè)置壓力傳感器、GPS、 Modem和下 位機(jī)，通過網(wǎng)絡(luò)傳輸將采集的數(shù)據(jù)傳至上位機(jī)監(jiān)控系統(tǒng)中，由SCADA系統(tǒng)實(shí)時(shí)采集管道 中燃?xì)獾牧髁?、壓力，由GIS系統(tǒng)儲(chǔ)存燃?xì)夤艿赖膶傩詳?shù)據(jù)，包括管線、管長、管徑、 管壁厚，利用管網(wǎng)仿真計(jì)算進(jìn)行管道的泄漏自動(dòng)檢測和管道泄漏點(diǎn)的定位。
管網(wǎng)仿真計(jì)算的方法是將燃?xì)夤芫€按長度方向劃分網(wǎng)格(認(rèn)為管線截面方向燃?xì)?參數(shù)一致)，將管網(wǎng)劃分成時(shí)間和長度的二維網(wǎng)絡(luò)，通過離散網(wǎng)格，列取節(jié)點(diǎn)方程，并 將輸入管網(wǎng)的邊界條件也離散成一系列的節(jié)點(diǎn)方程，通過相應(yīng)的數(shù)值計(jì)算，對正常工況 下管網(wǎng)的壓力、溫度、流量參數(shù)信息進(jìn)行實(shí)時(shí)仿真計(jì)算
I0 = A0b Ab。 (B0bEb + el0) / Rb + A。bCb Ab0 /+ A0b AbcBcbEb / Rb
Ic = A; Ab0 (B0bEb + el0) / Rb + AcbCb Ab0 (c/el0 /刮+ A0b AbcBcbEb / Rb
e2。 = B0bLbBb0玍+ B0bLbBbc工 將式1° = A°bAb。 (B。bEb + ell() / Rb + A°bCbAb°(^。 / A) + A°bAbcBcbEb / Rb轉(zhuǎn)換為
,=[A。bCbAb。]-1 [I0 — A。bAb。(B。bEb +0/Rb - A。bAbcBcbEb /Rb]
式中A°b、 Ab°、 Abe、 Aeb、 Bb。、 B。b、 Bcb、 Bbc——轉(zhuǎn)換系數(shù)；
L——管道的流感，定義為管道兩端引起的壓力變化與流量變化率之比值，
R——管道的流阻，定義為管道壓差與流量之比值，
C——管道的流容，定義為流體質(zhì)量變化與引起變化的壓力變化之比值，
E——協(xié)變張量，表示壓力源，
e——協(xié)變張量，表示管道壓降，
et。和e2。表示電路斷開時(shí)左右兩端的管道壓降，
式中上、下角標(biāo)符號o和c分別表示電路斷開和閉合兩種情況；b表示該張量為原始張量。
通過龍格一庫塔(Runge-Kutta)法對上述方程進(jìn)行求解。給定管網(wǎng)結(jié)構(gòu)，即可確 定式中的轉(zhuǎn)換系數(shù)A。b Ab° A^A:;和Bb。B。bBebBbe，以及與時(shí)間無關(guān)的電感系數(shù)L 。動(dòng)態(tài) 仿真計(jì)算時(shí)的初始條件由穩(wěn)態(tài)過程的計(jì)算結(jié)果確定。由得到的e,。、 1°求得1、再由I。和 得到e,。的值。根據(jù)求得的r，以及e,。和e,。的值，通過轉(zhuǎn)化即可得到管道的氣體流量M
及壓力p，詳見圖1。
管道的泄漏自動(dòng)檢測和管道泄漏點(diǎn)的定位，具體方法如以下步驟 I .建立燃?xì)夤芫W(wǎng)地理信息系統(tǒng)，讀取和存儲(chǔ)燃?xì)夤芫W(wǎng)圖、管網(wǎng)屬性數(shù)據(jù)及空間地理
數(shù)據(jù)；
n.由數(shù)據(jù)采集與監(jiān)視控制系統(tǒng)實(shí)時(shí)采集每段燃?xì)夤艿赖膲毫?、溫度、流量參?shù)信息，
對于壓力信息采集的頻率高于1HZ;
ni.對采集到的壓力、溫度、流量參數(shù)信息通過網(wǎng)絡(luò)傳輸至上位機(jī)監(jiān)控系統(tǒng)中，與 地理信息系統(tǒng)預(yù)先儲(chǔ)存的相應(yīng)管網(wǎng)屬性數(shù)據(jù)進(jìn)行比較，對管網(wǎng)進(jìn)行實(shí)時(shí)仿真計(jì)算，得到 管網(wǎng)安全運(yùn)行條件下各個(gè)管段節(jié)點(diǎn)處燃?xì)獾牧髁?、壓力值，并將?shí)測值和仿真計(jì)算值進(jìn)
行比較；
IV. 計(jì)算各管道首末端處測量值和計(jì)算參數(shù)值之間的偏差，當(dāng)流量和壓力的偏差大于 預(yù)定的閥值，由上位機(jī)監(jiān)控系統(tǒng)報(bào)告顯示出現(xiàn)異常的管道信息；
V. 搜索極值點(diǎn)確定壓力突變點(diǎn)，用小波變換對異常管道的首末端壓力信息進(jìn)行分析 和處理，得到首末端的壓力突變點(diǎn)，具體為
通過求取小波系數(shù)在各尺度下的最大、最小值，獲取所對應(yīng)的時(shí)間上的采樣點(diǎn)數(shù)nu 和n12;
采集的壓力信號是以時(shí)間為橫坐標(biāo)，壓力為縱坐標(biāo)的曲線。nu和n!2指的是時(shí)間上 的采樣點(diǎn)數(shù)。比如采樣頻率是5Hz，即每0.2秒采樣一次。若1111=250， nl2=50，則 nn和ru2之間有(250—50) 0.2 =40s，就是說時(shí)間寬度是40s。 ^T(a，6)—指各尺度下
的小波變化系數(shù)。
②搜索在nn和ru2之間，由等效帶通濾波器帶寬比較大的第一個(gè)尺度分解得到高 頻信號尺度下的全部極值點(diǎn)。③確定各管道中的壓力突變點(diǎn)；
找出首末端奇異點(diǎn)對應(yīng)的時(shí)刻，計(jì)算時(shí)間差。
VI.根據(jù)發(fā)生泄漏管段的燃?xì)饬髁?、溫度、壓力參?shù)，通過泄漏定位公式計(jì)算泄漏 點(diǎn)位置，泄漏定位公式為
—丄(V-i^) + (V — Mj)(V + M2)X
2v — Wj + "2
x——泄漏點(diǎn)距異常管首段的距離，m L——異常管長度，m
v——管道傳輸介質(zhì)中負(fù)壓波的傳播速度，m/s W——異常管首段到泄漏點(diǎn)的天然氣流速，m/s u2——異常管首段到泄漏點(diǎn)的天然氣流速，m/s △t——負(fù)壓波傳播到上、下游傳感器的時(shí)間差，s。 對管道中燃?xì)饬魉俨捎梅侄纹骄魉賮磉M(jìn)行計(jì)算，具體公式為
"CP—燃?xì)馄骄魉?，m/s
—燃?xì)赓|(zhì)量流量，kg/s
一平均壓縮因子，
一氣體參數(shù)，
一燃?xì)馄骄鶞囟?K
—燃?xì)馄骄鶋毫Γ琍a
A~"異常管段截面積，m2。
圖1是管網(wǎng)仿真計(jì)算的方法流程圖。 圖2是本發(fā)明實(shí)施例奇異點(diǎn)搜索的算法流程圖。
圖2中山指的是等效帶通濾波器帶寬比較大的第一尺度分解得到的高頻信號；dW 指得是Daubechies母小波函數(shù)。
具體實(shí)施例方式
以下結(jié)合燃?xì)夤芫W(wǎng)泄漏自動(dòng)檢測和定位的實(shí)施例，對本發(fā)明做進(jìn)一步的說明。本實(shí)施例以天津市長度為62.1公里管線中的一段，即永清(首站)至王慶陀(末站) 之間36.2km長的一段燃?xì)夤苈愤M(jìn)行泄漏定位分析。
I. 利用燃?xì)夤芫W(wǎng)GIS系統(tǒng)，讀取和存儲(chǔ)燃?xì)夤芫W(wǎng)圖、管網(wǎng)屬性數(shù)據(jù)及空間地理數(shù) 據(jù)。 '
II. 由SCADA系統(tǒng)實(shí)時(shí)釆集每段燃?xì)夤艿赖膲毫?、溫度、流量參?shù)信息，壓力信息 采集的頻率為5HZ。
III. 對管網(wǎng)進(jìn)行實(shí)時(shí)仿真計(jì)算，得到管網(wǎng)安全運(yùn)行條件下各個(gè)管段節(jié)點(diǎn)處各個(gè)時(shí)刻燃 氣的流量、壓力值。上位機(jī)監(jiān)控系統(tǒng)于某日下午2:05檢測到首站流量為2.42kg/s，壓力 為1.428MP4末站流量為1.85kg/s，壓力為1.283MPa。而根據(jù)管網(wǎng)仿真軟件計(jì)算結(jié)果顯 示，首站流量為2.43kg/s，壓力為1.454MPa,末站流量為2.35kg/s，壓力為1.413MPa。
IV. 設(shè)定流量閥值為5X，壓力閥值為5%，末端流量相對偏差為
x 100% = 21.3% > 5%
2.35
末端壓力相對偏差為-
L413 —L283x100% = 9.2o/o>5% 1.413
流量壓力均偏差大于預(yù)定的閥值，主監(jiān)控室報(bào)警有異常情況，由上位機(jī)監(jiān)控系統(tǒng)報(bào) 告顯示永清(首站)至王慶陀(末站)的燃?xì)夤苈樊惓！?br>
V. 搜索極值點(diǎn)確定壓力突變點(diǎn)，用小波變換對異常管道的首末端壓力信息進(jìn)行分析
和處理，得到首末端的壓力突變點(diǎn)，具體為
① 通過求取小波系數(shù)在各尺度下的最大、最小值獲取它們對應(yīng)的nu， n12; 通過計(jì)算首站的nn、 012分別為
=應(yīng)(》PfT(fl,柳=1762
巧2 = min(S(『r(a,柳=1739
末站的nu、 2分別為
= max(2(肝(a,6))) = 1605 巧2 = min(S (rr(a,柳=1563
② 找出首末站在nn、 1112之間，在等效帶通濾波器帶寬比較大的第一個(gè)尺度分解所 得到的高頻信號尺度下的全部極值點(diǎn)；
8結(jié)果表明首站的所有極值點(diǎn)分別為n產(chǎn)1762，n^l754，nf1750， ri4=1739，末站的所 有極值點(diǎn)分別為nfl605, n2=1580， n3=1574， n4=1563。 ③確定各管道中的壓力突變點(diǎn)；
根據(jù)壓力突變點(diǎn)在小波變化各尺度下的系數(shù)應(yīng)該保持不變，而偽突變點(diǎn)的小波變換 系數(shù)隨尺度的增加迅速衰減的特性，對首末站所有極值點(diǎn)逐一進(jìn)行分析，得到首末站的 壓力突變點(diǎn)分別為n3和ri2。
找出首末端奇異點(diǎn)對應(yīng)的時(shí)刻，計(jì)算時(shí)間差；
由于泄漏會(huì)使異常管段產(chǎn)生一個(gè)瞬時(shí)的壓力突變，在壓力信號中表現(xiàn)為壓力突然下 降。但由于實(shí)際燃?xì)夤艿来嬖诠I(yè)噪聲以及泵的起停等因素的干擾，這些干擾也會(huì)對壓 力信號產(chǎn)生很大影響，因此采用小波變換對異常管道的首末站壓力信號進(jìn)行分析和處 理，從而得到實(shí)際的壓力突變點(diǎn)。
通過計(jì)算得到首末站檢測到負(fù)壓波的采樣數(shù)分別為1750和1580，則首末站檢測到 負(fù)壓波時(shí)間差為At二 (1750 — 1580) X0.2 = 34s。每0.2秒采樣一次。
VI.根據(jù)發(fā)生泄漏管段的燃?xì)饬髁?、溫度、壓力參?shù)，通過泄漏定位公式計(jì)算泄漏 點(diǎn)位置
2v — w+ "2
由SCADA系統(tǒng)采集的數(shù)據(jù)和GIS系統(tǒng)儲(chǔ)存的管道屬性信息顯示，L = 62.1km， v=320m/s， & =7.62m/s, m2 二7.53m/s,將以上數(shù)據(jù)全部代入泄漏定位公式進(jìn)行計(jì)算
丄=丄(Hi) + (v-^)(v + A)xA, —35751m 2v — ^ + w2
即泄漏點(diǎn)距首站35.751km，絕對誤差為36.20—35.751 =449m，相對誤差為1.2%。 本發(fā)明的特點(diǎn)及有益效果是，通過軟件系統(tǒng)如SCADA、 GIS、管網(wǎng)仿真計(jì)算等與硬件 設(shè)施的的配合，可以實(shí)現(xiàn)燃?xì)夤芫€的自動(dòng)監(jiān)測及定位功能，完成泄漏的快速檢測和報(bào)警， 迅速發(fā)現(xiàn)燃?xì)庑孤┪恢?，減少由泄漏引起的損失。具有檢測速度快，檢測精度高，漏報(bào) 率/誤報(bào)率低等優(yōu)點(diǎn)。同時(shí)本發(fā)明也可用于供油、供水等管道，具有較大的經(jīng)濟(jì)效益和深 遠(yuǎn)的社會(huì)效益。
權(quán)利要求
1.燃?xì)夤艿佬孤z測和定位方法，具有壓力傳感器、數(shù)據(jù)采集與監(jiān)視控制系統(tǒng)、全球定位系統(tǒng)、A/D模數(shù)轉(zhuǎn)換卡、下位機(jī)、上位機(jī)監(jiān)控系統(tǒng)，其特征是在一段燃?xì)夤艿赖膬啥嗽O(shè)置壓力傳感器、全球定位系統(tǒng)、A/D模數(shù)轉(zhuǎn)換卡和下位機(jī)，通過網(wǎng)絡(luò)傳輸將采集的數(shù)據(jù)傳至上位機(jī)監(jiān)控系統(tǒng)中，由數(shù)據(jù)采集與監(jiān)視控制系統(tǒng)實(shí)時(shí)采集管道中燃?xì)獾牧髁?、壓力，由地理信息系統(tǒng)儲(chǔ)存燃?xì)夤艿赖膶傩詳?shù)據(jù)，包括管線、管長、管徑、管壁厚，通過管網(wǎng)仿真計(jì)算進(jìn)行管道的泄漏自動(dòng)檢測和管道泄漏點(diǎn)的定位，具體方法如以下步驟I.建立燃?xì)夤芫W(wǎng)地理信息系統(tǒng)，讀取和存儲(chǔ)燃?xì)夤芫W(wǎng)圖、管網(wǎng)屬性數(shù)據(jù)及空間地理數(shù)據(jù)；II.由數(shù)據(jù)采集與監(jiān)視控制系統(tǒng)實(shí)時(shí)采集每段燃?xì)夤艿赖膲毫?、溫度、流量參?shù)信息，對于壓力信息采集的頻率高于1HZ；III.對采集到的壓力、溫度、流量參數(shù)信息通過網(wǎng)絡(luò)傳輸至上位機(jī)監(jiān)控系統(tǒng)中，與地理信息系統(tǒng)預(yù)先儲(chǔ)存的相應(yīng)管網(wǎng)屬性數(shù)據(jù)進(jìn)行比較，對管網(wǎng)進(jìn)行實(shí)時(shí)仿真計(jì)算，得到管網(wǎng)安全運(yùn)行條件下各個(gè)管段節(jié)點(diǎn)處燃?xì)獾牧髁俊毫χ?，并將?shí)測值和仿真計(jì)算值進(jìn)行比較；IV.計(jì)算各管道首末端處測量值和計(jì)算參數(shù)值之間的偏差，當(dāng)流量和壓力的偏差大于預(yù)定的閥值，由上位機(jī)監(jiān)控系統(tǒng)報(bào)告顯示出現(xiàn)異常的管道信息；V.搜索極值點(diǎn)確定壓力突變點(diǎn)，用小波變換對異常管道的首末端壓力信息進(jìn)行分析和處理，得到首末端的壓力突變點(diǎn)，具體為①通過求取小波系數(shù)在各尺度下的最大、最小值，獲取所對應(yīng)的時(shí)間上的頻數(shù)n11和n12；n11＝max(∑(WT(a，b)))n12＝min(∑(WT(a，b)))WT(a，b)——各尺度下的小波變化系數(shù)；②搜索在n11和n12之間，由等效帶通濾波器帶寬比較大的第一個(gè)尺度分解得到高頻信號尺度下的全部極值點(diǎn)；③確定各管道中的壓力突變點(diǎn)；④找出首末端奇異點(diǎn)對應(yīng)的時(shí)刻，計(jì)算時(shí)間差；VI.根據(jù)發(fā)生泄漏管段的燃?xì)饬髁?、溫度、壓力參?shù)，通過泄漏定位公式計(jì)算泄漏點(diǎn)位置，泄漏定位公式為<maths id="math0001" num="0001" ><math><![CDATA[ <mrow><mi>x</mi><mo>=</mo><mfrac> <mrow><mi>L</mi><mrow> <mo>(</mo> <mi>v</mi> <mo>-</mo> <msub><mi>u</mi><mn>1</mn> </msub> <mo>)</mo></mrow><mo>+</mo><mrow> <mo>(</mo> <mi>v</mi> <mo>-</mo> <msub><mi>u</mi><mn>1</mn> </msub> <mo>)</mo></mrow><mrow> <mo>(</mo> <mi>v</mi> <mo>+</mo> <msub><mi>u</mi><mn>2</mn> </msub> <mo>)</mo></mrow><mo>×</mo><mi>Δt</mi> </mrow> <mrow><mn>2</mn><mi>v</mi><mo>-</mo><msub> <mi>u</mi> <mn>1</mn></msub><mo>+</mo><msub> <mi>u</mi> <mn>2</mn></msub> </mrow></mfrac> </mrow>]]></math></maths>x——泄漏點(diǎn)距異常管首段的距離，mL——異常管長度，mv——管道傳輸介質(zhì)中負(fù)壓波的傳播速度，m/su1——異常管首段到泄漏點(diǎn)的天然氣流速，m/su2——異常管首段到泄漏點(diǎn)的天然氣流速，m/sΔt——負(fù)壓波傳播到上、下游傳感器的時(shí)間差，s。
2.根據(jù)權(quán)利要求1所述的燃?xì)夤艿佬孤z測和定位方法，其特征是對管道中燃?xì)饬魉俨捎梅侄纹骄魉賮磉M(jìn)行計(jì)算，具體公式為—W W; ""一 ^^ ——燃?xì)馄骄魉?，m/s 《 ——燃?xì)赓|(zhì)量流量，kg/s Zw——平均壓縮因子， ——?dú)怏w參數(shù)，r￡/7——燃?xì)馄骄鶞囟?，k——燃?xì)馄骄鶋毫Γ琍a A~—異常管段截面積，m2。
全文摘要
本發(fā)明公開了一種基于地理信息(GIS)和數(shù)據(jù)采集與監(jiān)控(SCADA)技術(shù)的對燃?xì)夤芫W(wǎng)泄漏自動(dòng)檢測和定位的方法。具體方法如以下步驟建立燃?xì)夤芫W(wǎng)GIS系統(tǒng)，讀取和存儲(chǔ)管網(wǎng)圖、管網(wǎng)屬性數(shù)據(jù)等；由SCADA系統(tǒng)實(shí)時(shí)采集每段燃?xì)夤艿赖膲毫?、溫度、流量參?shù)；對采集參數(shù)與儲(chǔ)存數(shù)據(jù)比較計(jì)算；計(jì)算管道首末端處測量值和計(jì)算參數(shù)值之間的偏差，顯示出現(xiàn)異常的管道信息；搜索極值點(diǎn)確定壓力突變點(diǎn)，根據(jù)發(fā)生泄漏管段的燃?xì)饬髁?、溫度、壓力參?shù)，通過泄漏定位公式計(jì)算泄漏點(diǎn)位置。本發(fā)明通過SCADA、GIS以及仿真計(jì)算，可以實(shí)現(xiàn)燃?xì)夤芫€的自動(dòng)監(jiān)測及定位，具有檢測速度快，檢測精度高，漏報(bào)率/誤報(bào)率低等特點(diǎn)。
文檔編號F17D5/00GK101625071SQ20091007010
公開日2010年1月13日 申請日期2009年8月7日 優(yōu)先權(quán)日2009年8月7日
發(fā)明者昭 楊, 壯 熊, 敏 邵 申請人:天津大學(xué)