基于3S集成的突發水污染模擬研究

楊明祥，解建倉，李建勛，謝玉峰

（西安理工大學，710048，西安）

基于3S集成的突發水污染模擬研究

楊明祥，解建倉，李建勛，謝玉峰

（西安理工大學，710048，西安）

針對當前水污染研究多局限于模型開發，可視化工作較為薄弱的問題，提出了一種基于3S集成的突發水污染模擬方案，并討論了河流三角化、離散點濃度擬合、基于視域的數據動態加載等關鍵技術，對松花江某取水口上游突發石油污染事件進行了動態模擬。結論表明，該方法可以形象地模擬污染物的擴散和運移過程，并提供定點的污染監測功能，為污染模型的檢驗提供了依據，為應急方案的制定提供了全方位的信息支持。

3S集成；三維可視化；水污染擴散

近年，我國GDP保持了較快的增長勢頭，經濟建設突飛猛進，但在經濟發展的同時，各種環境問題日益突出，特別是突發水污染更是頻繁發生。由于突發水污染具有影響范圍大、不確定性因素多、往往距離人口稠密區較近等特點，對人民生命財產和生態平衡造成了巨大威脅，制約了經濟社會發展。因此，如何有效地分析突發水污染事件，制定有效地應急措施成為社會關注的熱點問題。目前為止，國內對于突發水污染的研究多局限于污染物遷移擴散的模型開發，對其運移過程的可視化模擬還非常薄弱，特別是三維的可視化模擬方面，還沒有成熟的系統出現。這樣，就導致了現有模型多處于孤立狀態，沒有與其配套的驗證和展示系統，不能很好地指導應急方案的制定。本文提出了基于3S集成的突發水污染模擬方案，利用3S集成平臺的坐標和影像系統來實現污染物的定位及其周邊地形的展示服務，為應急方案的制定提供全方位的信息支持，達到提高方案制定速度和質量的目的。

一、3S集成下的水污染模擬體系結構

1.3S集成平臺

3S技術是指遙感技術（RS，remote sensing）、地理信息系統（GIS，geographic information system）和全球定位 系 統 （GPS，global positioning system）。3S集成則是三者的有機結合，其集成可以提供一些特定的功能，遙感（RS）和全球定位系統（GPS）的集成，可實現河流周邊地形的展示。而地理信息系統（GIS）和其他各數據庫的聯合集成將構成一個功能強大的數據存儲和分析平臺，可以為河流水污染模擬提供強有力的技術支撐和保障。

水利行業自20世紀80年代初開始應用遙感（RS）技術，即通過對地觀測獲取信息。對GIS的使用則始于20世紀80年代后期，在經歷了認識了解和初步應用兩個階段后，現已步入深入應用的階段，且很快就與生產實際緊密地結合起來。全球定位系統（GPS）在水利行業的應用始于20世紀90年代初，但發展非常迅速，在地面及水下地形測繪中使用已很普遍。可見，水利行業經過多年的發展，已經積累了大量的空間信息資源，完全有條件開發和利用3S集成平臺解決相關問題，特別是突發水污染問題，由于其涉及地域面積廣，影響時間長，有明顯的時空分布特點，僅靠模型計算難以實現對全局信息的掌握，很有必要引入3S集成平臺，對污染數據進行全方位的信息展示和分析。

2.三維水質模擬體系結構

3S集成平臺上三維水質模擬包含三個方面的內容：河面三維仿真、影像數據展示和污染物擴散運移可視化，三者通過地理空間坐標相聯系形成污染物擴散運移可視化展示平臺，在此基礎上進行各種統計分析。系統結構如圖1所示。

二、關鍵技術

1.河流三角化

為了實現河流的三維仿真，需要對河流進行三角化，本文提出了一種基于經緯方格的剖分算法來實現三角化。其步驟如下所述：

①解析河道二維GIS數據，獲取待處理區域的經緯方格，如圖2所示。經線（緯線）之間的間距越小，則單個經緯方格所圍成的面積越小，代表河面被剖分得越為精細。

②求出經緯方格與河岸的所有交點，以及位于兩岸之間的方格頂點，如圖2中黑色圓點即為所求的點。該過程的核心是求經緯方格與河岸的交點。

③將相鄰兩經線上的點依次連接成三角形，如圖2所示。

圖1 突發水污染模擬系統結構圖

圖2 河面三角網

經過以上三個步驟，實現了河面的三角化，但這僅僅是二維平面的三角化，要實現河流的三維空間模擬，還要利用3S集成平臺中的DEM模型，根據空間經緯度信息獲取高程，這樣，每個三角形的三個頂點都被賦予了不同的高程信息，模擬出的河流可以表現出實際的地貌信息。

2.離散點濃度擬合

河流的三角化為水質仿真提供了結構基礎，通過將網格點賦予不同的濃度信息，可以實現不同污染情況的表現，然而污染物擴散模型計算出來的結果是離散點的濃度信息，這些點在空間上并不與網格點一一對應，需要將其擬合到網格點上才能實現三維的水質模擬。由于河道狹長且蜿蜒曲折，傳統單一的擬合方法不能準確反映濃度信息，本文借鑒分段的思想，每隔一定距離將河道分割成一段，在每一個河段內使用最小二乘法擬合出網格點的濃度。需要說明的是，在河段較長的情況下，這種方法在斷面處會形成較為明顯的濃度跳變現象，所以在斷面處需要使用線性插值進行濃度平滑。

3.水質仿真可視化

通過擬合獲取網格點的污染物濃度后，再進行水質的仿真可視化。水質仿真可視化是指通過圖形圖像的方式對水質模型仿真過程進行跟蹤、駕馭和結果的顯示，同時實現整個過程的可視化，提供迅速、直觀、高效和形象的模型展示效果。本文所述仿真可視化過程借助3S集成平臺來實現，用戶可以通過鼠標在屏幕上直觀形象地進行操作 （開始、暫停、加速、減速等），根據實際需要控制整個仿真過程。污染物濃度與展示顏色之間的映射是可視化過程中的核心，其含義是把計算機數值模擬的數據轉換為可供繪制的幾何圖素和屬性，它決定在最后的圖像中應該看到什么，又如何將其表現出來。

本文中，水質按被污染程度劃分為5個等級，分別為Ⅰ類水、Ⅱ類水、Ⅲ類水、Ⅳ類水、Ⅴ類水。Ⅰ類水水質高錳酸鹽指數(COD)濃度0～2.5，Ⅱ類水2.5～4.0，Ⅲ類水4.0～6.0，Ⅳ類水6.0～10.0，Ⅴ類水 10.0～15.0，通過不同的顏色表示被污染的程度，最好的Ⅰ類水質對應色帶中藍色部分，愈往右走，表示被污染的程度愈高，右端紅色表示污染程度最嚴重，為Ⅴ類水。根據顏色和濃度的映射關系，將網格點著色，然后以三角形為單位，利用三個頂點的顏色信息進行平滑渲染，實現濃度信息向圖像信息的轉化，多組濃度值映射成一個圖像序列，根據間隔的時間將這些真彩色圖像順序播放成動畫，即實現了水污染過程中污染物的擴散演進模擬。

4.分析統計

突發水污染影響面廣，危害大，發展瞬息萬變，在污染發生后，需要根據污染物發展狀態制定應急方案，最大限度降低損失。使用三維圖形圖像技術模擬污染物運移和擴散可以為應急方案的制定提供直觀、形象和全景式的信息支持，但污染過程中的某些定量信息也很重要，例如對某個特定取水口的保護措施制定中，決策者需要了解的信息是取水口開始被污染的時間、污染物濃度達到峰值的時間以及水質達標的時間，所以還需要提供污染信息的分析統計功能，主要包括定點污染統計和重度污染區影響統計分析。

定點污染統計：水面上某一確定點在污染事件發生后其污染物濃度的變化。這里可以直接對污染擴散模型計算出來的某個時間點的成果文件進行分析。實際計算時，在成果文件中找出定點所在的三角形，然后根據三角形頂點差值計算出定點的濃度值。利用這種方法計算出各個時間點的濃度值，按照時間序列繪制濃度變化曲線，從而可以直觀地看到定點水質的變化過程，為應急方案制定提供定量的數據支持。

重度污染區影響統計分析：在特定時間點對Ⅳ類、Ⅴ類水影響區域進行分析統計，主要包括影響到的取水口、灌區和生態保護區等。首先，將濃度達到Ⅳ類水及其以上的網格點找到，獲取重度污染區域（不止一個），然后獲取各個區域的邊界 （由網格點組成），最后根據這些重度污染區域的邊界對保護對象進行篩選，即可得到該時間點被重度污染的保護對象。這些信息可以幫助決策者判斷污染事件的危害和影響程度，正確決策，盡可能降低重要保護目標的受污染程度。

5.基于視域的數據動態加載

水質模擬過程中，網格點對內存的耗費是不可忽視的。一個待模擬的網格點包括以下信息：經度、緯度和顏色，經度和緯度使用浮點數表示，顏色使用RGB值表示，故一個點占用28個字節，如果緯線和經線的間隔都取為5m，對50km長、平均300m寬的河道進行剖分，一個時間點的數據需要占用16MB的空間，這樣大的數據量在兩幅圖像更替的時候必然導致畫面的延遲，然而水質仿真需要達到平滑的顏色漸變效果，對實時性要求較高，是不允許出現延遲現象的。在系統運行過程中，人們能觀察到的是一個視域范圍內的圖形圖像信息，視域外的信息對用戶沒有任何意義，于是考慮對網格點數據進行基于視域的動態加載和卸載：在圖像刷新的時候，判斷當前視域經緯度的最大值和最小值，只有處于這個范圍內的信息才被加載，反之則被卸載。通過這種數據動態加載的方法，實現了內存使用量的降低，從而達到了畫面的流暢播放。獲取視域經緯度范圍的關鍵代碼如下所示：

dc是3S集成平臺運行中的上下文對象，包含平臺當前的基本信息。

三、應用實例

松花江流域跨黑龍江省、吉林省和內蒙古自治區，匯入黑龍江后，經俄羅斯注入太平洋。流域全長2214.3 km，流域面積 55.68萬 km2，流域總人口 5669 萬（2000年），位居長江、黃河之后，為我國第三大河，是我國重要的工業、農業、能源基地。隨著我國振興東北老工業基地戰略的實施，松花江流域的工業經濟得到了再次發展，但與此同時，突發水污染事故也頻頻發生，嚴重威脅流域內人民生命財產安全，給國民經濟造成了巨大的損失。如何及時應對和處理流域內的突發水污染事件，保護取水口等重要目標，成為當前治理松花江流域污染的關鍵。

松花江突發水污染水質時空分布模擬系統是利用現代網絡通信、計算機、信息管理技術，并借助3S集成平臺對水質的時空分布信息進行展示，為應急決策的制定提供直觀且全方位的信息支持。系統主要模塊及功能見表1。

系統以JAVA語言開發，采用C/S模式，實現了對松花江干流突發水污染水質時空分布情況的動態可視化模擬。假設位于“××”取水口上游5 km處，發生了油罐車翻車事故，有10m3油料進入河道，在水質數學模型支持下模擬油膜的擴散、分散過程，經系統模擬在翻車事故發生后的1小時9分鐘取水口附近水質變為Ⅳ類水，在2小時42分鐘后水質達到Ⅲ類水標準，從模擬圖像可以看出，污染物呈帶狀分布，濃團位于河流中間部分，污染區域逐漸變大，且有逐漸稀釋的趨勢，與實際情況相符。

表1 系統主要功能

四、結 論

本文分析了突發水污染的特點以及我國的研究現狀，闡述了傳統水污染仿真手段的不足——缺乏空間定位信息、脫離實際地貌，提出了基于3S集成平臺對污染物運移和擴散進行可視化展示的方法。實例研究表明，該方法形象直觀地展示了污染物遷移擴散過程，為檢驗污染模型的準確性提供了參考和依據，也為快速有效處置突發水污染事件提供了全新的思路和有力的技術支持。

[1]韓曉剛,黃廷林.我國突發性水污染事件統計分析[J].水資源保護,2010,26(1)：84-90.

[2]尹剛.基于3S技術的河流水污染監測信息系統的構建 [J].裝備環境工程,2009,6(4)：87-91.

[3]張波,王橋,等.基于SD-GIS的突發水污染事故水質時空模擬[J].武漢大學學報(信息科學版),2009,34(3)：348-351.

[4]曹增節.數字地球空間疊加與信息可視化[J].浙江工業大學學報,2002,30(3)：256-260.

[5]葉松,譚德寶,等.水污染擴散模擬三維可視化研究 [J].系統仿真學報,2008,20(16)：4451-4453.

[6]董志穎,李兵,等.水質預警理論初探[J].水資源研究,2002,23(1)：36-38.

[7]丁賢榮,徐健,等.GIS與數模集成的水污染突發事故時空模擬[J].河海大學學報(自然科學版),2003,31(2)：203-206.

[8]陳志云.松花江區水資源質量評價綜述[J].水資源保護,2007,23(1)：6-8.

[9]金春久,王超,等.松花江干流水質模型在流域水資源保護管理中的應用[J].水利學報,2010,41(1)：86-92.

責任編輯 鄧淑珍

The research of sudden water pollution simulation based on 3 S integration

Yang Mingxiang,Xie Jiancang,LiJianxun，XieYufeng

This paper pointed out the problems that current water pollution research more confined to model development,and visualization work was weak.Then,it put forward a method of sudden water pollution simulation that based on 3S integration and discussed the river triangle melt,the discrete points concentration fitting,the data dynamic loading based on the horizon,and some key techniques.Simulated a sudden water pollution broke upstream a certain intake of songhua river.The conclusion show that this method can simulate the diffusion and transport process of pollutants better and provide the function of fixed-pointmonitoring,the basis for pollution model test and the comprehensive information support for emergencyplan formulation.

3Sintegration；3 Dvisuali zation；water pollution diffusion

國家863計劃資助項目(2006AA01A126)，國家自然科學基金資助項目（50979088）。

X522

A

1000-1123（2011）11-0012-03

2011-02-01

楊明祥，碩士研究生。