鄭艷華 潘新平 林榮列 王婉芝 陸奕海 黃清輝 陳銳賢

摘 要： 為了快速準(zhǔn)確地確定氣體泄漏源的位置，提出基于加權(quán)質(zhì)心算法的氣體源定位法。采用氣敏傳感器對目標(biāo)區(qū)域的氣體濃度進行檢測，并通過無線傳感器網(wǎng)絡(luò)把收集的數(shù)據(jù)反饋給協(xié)調(diào)器，協(xié)調(diào)器把打包好的數(shù)據(jù)傳遞給主處理器，處理器依靠測量的結(jié)果計算出氣體源位置，并可視化顯示。以乙醇蒸汽為實驗?zāi)繕?biāo)氣體，平均定位相對誤差小于5.9%，驗證了該方法的準(zhǔn)確性。

關(guān)鍵詞： 氣敏傳感器； 加權(quán)質(zhì)心算法； STM32處理器； ZigBee無線傳感器網(wǎng)絡(luò)； 氣體源定位； 相對誤差

中圖分類號： TN99?34； TP391 文獻標(biāo)識碼： A 文章編號： 1004?373X（2018）11?0073?04

Gas source location system based on weighted centroid algorithm

ZHENG Yanhua， PAN Xinping， LIN Ronglie， WANG Wanzhi， LU Yihai， HUANG Qinghui， CHEN Ruixian

（School of Physics and Electronic Engineering， Guangzhou University， Guangzhou 510006， China）

Abstract： The gas source location method based on weighted centroid algorithm is proposed to determine the position of gas leakage source quickly and accurately. The gas sensor is used to detect the gas concentration of the target area. The wireless sensor network is used to feed the collected data back to the coordinator. The data after packaging is transmit to the main processor by means of the coordinator. The processor can calculate the gas source location for visual display according to the measured results. The experiment of taking the ethanol vapor as the target gas is performed， which can prove that the average relative error of location is less than 5.9%， and the gas source position algorithm is correct.

Keywords： gas sensor； weighted centroid algorithm； STM32 processor； ZigBee wireless sensor network； gas source position； relative error

隨著工業(yè)的迅速發(fā)展，存儲氣體的倉庫越來越多，當(dāng)發(fā)生意外事件導(dǎo)致氣體發(fā)生泄漏，或者一些工廠的有毒危險氣體發(fā)生偷排時，如不能及時尋找泄漏的地點，將給人們的生活和安全帶來極大的隱患[1?4]。從20世紀(jì)60年代開始就有眾多研究人員對危險氣體的擴散規(guī)律進行了研究?，F(xiàn)有的氣體擴散模型主要有高斯模型、BM 模型、FEM3模型、箱及相似模型、三維現(xiàn)象傳遞模型、淺層模型以及板塊模型等[5?9]。上述模型都有各自的局限性，不同的擴散模型適用于不同類型的氣體。只有通過實驗測試才能找出特定氣體泄漏規(guī)律適用的模型[9?10]。本文利用傳感器網(wǎng)絡(luò)和加權(quán)質(zhì)心算法對氣體進行定位研究，并進行了實驗驗證。

1 氣體擴散模型

在無風(fēng)連續(xù)泄漏點源的擴散模型中，由菲克擴散定律可知[11]，擴散過程中濃度隨時間的變化率等于該處的擴散通量隨距離變化率的負(fù)值。進而得經(jīng)典的擴散方程：

[?C?t+u?C?x=Kx?2C?x2+Ky?2C?y2+Kz?2C?z2] （1）

式中：[C]為[C（x，y，z，t）]，表示[t]時刻[（x，y，z）]處氣體濃度，單位為mg/m3；[u]表示風(fēng)速，單位為m/s；[t]表示泄漏氣體擴散時間，單位為s；[Kx，Ky，Kz]分別表示[x，y，z]軸方向上的湍流擴散系數(shù)，單位為m2/s。

當(dāng)氣體恒以泄漏速度[Q]在[t0]時刻開始擴散，[Kx=Ky=Kz=K]，風(fēng)速[u=]0時，并設(shè)氣體泄漏源與坐標(biāo)系原點重合，則無風(fēng)連續(xù)泄漏點源的擴散模型可簡化為微分方程：

[?C?t=K?2C?x2+?2C?y2+?2C?z2] （2）

在初始條件[x=y=z=0]，[C→∞]和邊界條件[x=y=z=]0，[C→∞]下的解為：

[C（x，y，z，t）=Q4Kπd2πd2K（t-t0）∞e-?2d?] （3）

這里用[d]表示測量點到源點的距離。

又因研究的是穩(wěn)態(tài)下的模型，即達到平衡狀態(tài)，等價于[t→∞]時的情況，這時有：

[2πd2K（t-t0）∞e-?2d?=1] （4）

而在穩(wěn)態(tài)下，濃度僅是位置的函數(shù)，與時間[t]無關(guān)，定義[P=QK]為環(huán)境場源參數(shù)，所以式（3）轉(zhuǎn)化為：

[C（x，y，z）=P4πd] （5）

2 系統(tǒng)硬件設(shè)計

系統(tǒng)總體的硬件結(jié)構(gòu)框圖如圖1所示。在終端處包括ZigBee節(jié)點和氣敏傳感器；在客戶端處包括STM32F407核心處理器、SD卡模塊、溫度傳感器和TFT觸摸顯示屏。每個節(jié)點配置2個氣敏傳感器，通過節(jié)點處自身的CC2530處理器使用內(nèi)部ADC采集氣敏傳感器的電壓，把數(shù)據(jù)傳遞到協(xié)調(diào)器，協(xié)調(diào)器把每個節(jié)點的數(shù)據(jù)按順序打包壓縮，通過串口通信發(fā)送到STM32F407處理器。STM32通過中斷接收到串口傳遞的數(shù)據(jù)，再根據(jù)定位算法進行信號處理，將狀態(tài)信息和處理的結(jié)果保存在SD卡內(nèi)并在TFT屏幕上實時顯示。

本實驗采用FIGARO公司設(shè)計的酒精傳感器TGS2620。為了使輸入阻抗同其匹配，因而在傳感器模擬量輸出口加入跟隨器MCP6022，高阻抗信號輸入跟隨器阻抗變換后，輸出和A/D相匹配的低阻抗信號，將輸入對A/D的影響降至最低，提高對負(fù)載電阻電壓的采樣精度。

ZigBee節(jié)點的控制芯片使用TI公司的CC2530，內(nèi)部集成了以8051內(nèi)核為處理器和高性能的射頻收發(fā)器，內(nèi)置8通道12位A/D轉(zhuǎn)換器、看門狗和睡眠定時器等。其UART0引腳是P0_2，P0_3，用于和外部芯片進行通信，本實驗用其與核心處理器STM32串口通信。

中央控制器采用STM32F407單片機，可以運行μC/OS等實時操作系統(tǒng)和GUI圖形庫。而且其自帶的FSMC復(fù)用接口、SDIO接口和USB HOST接口有利于本設(shè)計高效快速地刷新TFT顯示和讀取SD卡/U盤文件。

3 系統(tǒng)的軟件設(shè)計

系統(tǒng)軟件設(shè)計分為無線傳感網(wǎng)絡(luò)的搭建、數(shù)據(jù)的采集與接收、數(shù)據(jù)處理方法和用戶界面設(shè)計4個部分。

3.1 無線傳感網(wǎng)絡(luò)的搭建

系統(tǒng)上電后，協(xié)調(diào)器開始工作，首先硬件以及軟件架構(gòu)的各模塊初始化，建立傳感器網(wǎng)絡(luò)。協(xié)調(diào)器完成初始化工作后就處于等待狀態(tài)，當(dāng)收到其他節(jié)點請求入網(wǎng)時，允許該節(jié)點入網(wǎng)并且為該節(jié)點分配網(wǎng)絡(luò)地址，其流程圖如圖2所示。一般來說，協(xié)調(diào)器的初始化工作包括硬件初始化、接口初始化以及軟件架構(gòu)各模塊的初始化。

3.2 數(shù)據(jù)的采集、傳送和接收

終端傳感器節(jié)點最重要的功能便是采集現(xiàn)場環(huán)境數(shù)據(jù)并且發(fā)送給協(xié)調(diào)器。根據(jù)終端傳感器節(jié)點的功能，可以將其軟件設(shè)計分為五個部分：系統(tǒng)初始化、時間同步、數(shù)據(jù)采集、傳輸以及睡眠。終端傳感器節(jié)點工作流程圖如圖3所示。

終端節(jié)點在網(wǎng)絡(luò)中除了選擇數(shù)據(jù)傳輸路徑、轉(zhuǎn)發(fā)數(shù)據(jù)包、延長通信距離外還有一個功能是網(wǎng)絡(luò)管理，維護無線傳感網(wǎng)絡(luò)的正常工作。終端節(jié)點上電后，發(fā)現(xiàn)網(wǎng)絡(luò)便請求入網(wǎng)并與第一個響應(yīng)它的協(xié)調(diào)器節(jié)點綁定，然后作為網(wǎng)絡(luò)中的路由節(jié)點等待其他的節(jié)點加入。

協(xié)調(diào)器節(jié)點收到終端傳感器節(jié)點與路由節(jié)點的數(shù)據(jù)信息后，通過數(shù)據(jù)包解析確定該數(shù)據(jù)為有效數(shù)據(jù)信息后，通過串口上傳給核心處理器。

3.3 基于加權(quán)質(zhì)心定位算法

定位算法步驟如下：

1） 設(shè)傳感器[i]的位置已知，為[（xi，yi，zi）]，測量濃度為[Ci，]則傳感器的權(quán)重為[qi=Cii=1nCi]，[n]為參與定位傳感器的個數(shù)。

可得預(yù)測的源點坐標(biāo)[（x，y，z）]為：

[x=i=1nxiqi=i=1nxiCii=1nCiy=i=1nyiqi=i=1nyiCii=1nCiz=i=1nziqi=i=1nziCii=1nCi] （6）

2） 求出各傳感器節(jié)點[（xi，yi，zi）]到預(yù)測源點[（x，y，z）]的距離[di]，由式（5）可得氣體環(huán)境場源參數(shù)：

[Pi=4πdiCi] （7）

最終，氣體環(huán)境場源參數(shù)可取為：

[P0=1ni=1nPi] （8）

3） 以初步定位坐標(biāo)[（x，y，z）]為中心，以定位區(qū)域邊長的10%為邊長畫出一個矩形區(qū)域作為精細(xì)定位區(qū)域，將此區(qū)域劃分為邊長為1 cm的網(wǎng)格，格點坐標(biāo)記為[（x，y，z）]。

4） 將由初步定位中預(yù)測所得的環(huán)境場源參數(shù)[P0]及測得的濃度值[Ci]代入式（6）分別計算理想源點[（x，y，z）]和每個格點[（x，y，z）]到各傳感器節(jié)點[（xi，yi，zi）]的距離[di]和[d′i]，記誤差為[E=i=1nd′i-di]。

5） 遍歷格點[（x，y，z）]，尋找使得[E]取最小值的格點，其坐標(biāo)[（x，y，z）]即為精細(xì)定位結(jié)果。

3.4 用戶界面設(shè)計

本系統(tǒng)采用ST_emWin圖形庫實現(xiàn)GUI界面。操作界面一共有5個應(yīng)用窗口，分別是實時數(shù)據(jù)采集窗口、濃度變化曲線窗口、定位結(jié)果二維平面顯示窗口、系統(tǒng)信息設(shè)置窗口和SD卡文件管理窗口。

4 氣體源定位實驗結(jié)果及分析

氣體源定位實驗在25 m×10 m的實驗室內(nèi)進行，密封實驗室門窗，形成相對穩(wěn)定的密閉實驗環(huán)境。在直徑為9 cm的容器內(nèi)倒入25 ml的酒精形成乙醇蒸汽源。傳感器陣列等距離放置在正方形四邊。檢測區(qū)域面積分別為2 m×2 m，3 m×3 m和4 m×4 m。在相同的檢測區(qū)域下改變氣體源坐標(biāo)，實驗結(jié)果與誤差如表1所示（[誤差=定位結(jié)果坐標(biāo)到實際坐標(biāo)的距離檢測區(qū)域的邊長]）。

系統(tǒng)每5 s采集數(shù)據(jù)并定位計算，結(jié)果以二維平面坐標(biāo)顯示在用戶界面上。定位結(jié)果如圖4所示，檢測區(qū)域為3 m×3 m，氣體源坐標(biāo)為（2.50，1.20） m，2次計算結(jié)果分別為（2.47，1.34） m和（2.65，1.23） m。

實驗結(jié)束后，由保存在SD卡的坐標(biāo)數(shù)據(jù)得到的定位坐標(biāo)效果如圖5所示。

系統(tǒng)在開始定位階段，氣體源還未真正擴散，每個傳感器測量的濃度都很低且相近，根據(jù)定位算法特性，系統(tǒng)定位的結(jié)果為檢測區(qū)域的正中心周圍。隨著氣體源的泄漏，靠近氣體源的傳感器濃度慢慢增大，計算的定位結(jié)果往氣體源點坐標(biāo)靠近。在氣體擴散過程的后期階段，源頭的酒精量減少，場源減弱，濃度下降，則定位結(jié)果靠近傳感器陣列的幾何中心。

在圖4中，靠右較大的二維平面圖是整個實驗過程中定位計算結(jié)果顯示，左上角小圖是第1 143次定位計算的結(jié)果坐標(biāo)顯示。在整個定位過程中，因為開始階段和后期階段測量濃度很低，導(dǎo)致定位結(jié)果靠近檢測區(qū)域的中心，不是真正定位結(jié)果。圖5坐標(biāo)顯示的是去除開始和后期階段數(shù)據(jù)，紅色小方點是整個實驗過程定位計算的結(jié)果，綠色菱形點是整個過程定位計算結(jié)果的平均值，藍(lán)色大方格點是氣體源坐標(biāo)點。

5 結(jié) 語

本文基于加權(quán)質(zhì)心算法和無線傳感器網(wǎng)絡(luò)對氣體源坐標(biāo)進行定位。同時設(shè)計了分布式氣敏傳感器陣列以及無線傳感器網(wǎng)絡(luò)采集數(shù)據(jù)，并使用ST?ewMin圖形庫設(shè)計用戶界面，讓定位結(jié)果以直觀的二維平面顯示。本實驗以酒精蒸汽為目標(biāo)氣體在密閉實驗室下進行氣體源定位實驗。結(jié)果證明本系統(tǒng)在一定檢測區(qū)域內(nèi)能夠進行定位，說明了本系統(tǒng)的可靠性。

參考文獻

[1] 陳平，石磊，趙建勇，等.氣體變電站故障泄漏預(yù)警回收技術(shù)研究與應(yīng)用[J].高壓電器，2016，52（12）：244?248.

CHEN Ping， SHI Lei， ZHAO Jianyong， et al. Research on early warning and recovery technology of gas leakage [J]. High voltage apparatus， 2016， 52（12）： 244?248.

[2] 史騰飛，管義鋒.基于PHAST的CNG運輸船裝卸載過程氣體泄漏后果評估[J].船舶工程，2016，38（12）：83?87.

SHI Tengfei， GUAN Yifeng. Consequence evaluation of gas leakage in CNG carrier loading and unloading process based on PHAST [J]. Ship engineering， 2016， 38（12）： 83?87.

[3] 孫凱.危險氣體無線監(jiān)測系統(tǒng)的設(shè)計及氣體源定位算法的研究[D].成都：西南交通大學(xué)，2016.

SUN Kai. Design of wireless monitoring system for dangerous gas and study on gas source localization [D]. Chengdu： Southwest Jiaotong University， 2016.

[4] CAO M L， MENG Q H， JING Y Q， et al. Distributed sequential location estimation of a gas source via convex combination in WSNs [J]. IEEE transactions on instrumentation & measurement， 2016， 65（6）： 1484?1494.

[5] 朱伯齡，於孝春，李育娟.氣體泄漏擴散過程及影響因素研究[J].石油與天然氣化工，2009，38（4）：354?358.

ZHU Boling， YU Xiaochun， LI Yujuan. Study on the diffusion process and influence factors of gas leakage [J]. Chemical engineering of oil and gas， 2009， 38（4）： 354?358.

[6] 林元棣，廖瑞金，夏桓桓，等.油紙絕緣微水二次擴散過渡過程的分析[J].電工技術(shù)學(xué)報，2016，31（10）：16?25.

LIN Yuandi， LIAO Ruijin， XIA Huanhuan， et al. Analysis of transition process during secondary diffusion of moisture in oil?paper insulation [J]. Transactions of China electrotechnical society， 2016， 31（10）： 16?25.

[7] 潘新平，鄭艷華，林榮列，等.基于氣體擴散模型的氣體源定位系統(tǒng)[J].國外電子測量技術(shù)，2017，36（9）：72?76.

PAN Xinping， ZHENG Yanhua， LIN Ronglie， et al. Gas source location system based on gas diffusion model [J]. Foreign electronic measurement technology， 2017， 36（9）： 72?76.

[8] CAO M L， MENG Q H， JING Y Q， et al. Distributed sequential location estimation of a gas source via convex combination in WSNs [J]. IEEE transactions on instrumentation & measurement， 2016， 65（6）： 1484?1494.

[9] WADA Y， MATSUKURA H， ISHIDA H. Estimation of gas source location from fluctuating readings of gas sensors and anemometer on mobile robot in outdoor environment [J]. ECS transactions， 2016， 75（16）： 99?106.

[10] 馮志華，聶百勝.危險氣體泄漏擴散的實驗方法研究[J].中國安全科學(xué)學(xué)報，2006，16（7）：18?23.

FENG Zhihua， NIE Baisheng. Experimental methodology for study diffusion of leaked hazardous gas [J]. China safety science journal， 2006， 16（7）： 18?23.

[11] 陳立偉，楊建華，孫亮，等.基于分布式傳感器陣列的靜態(tài)氣體源定位方法[J].電子科技大學(xué)學(xué)報，2014，43（2）：212?215.

CHEN Liwei， YANG Jianhua， SUN Liang， et al. Odor source localization algorithm based on spatially distributed sensors array [J]. Journal of University of Electronic Science and Technology of China， 2014， 43（2）： 212?215.