大氣污染物擴(kuò)散高斯模型模擬

1、92.2大氣污染物擴(kuò)散的高斯模型模擬：可視化模擬點(diǎn)源大氣污染的擴(kuò)散9.2.2 Gaussian Atmospheric Dispersion Model突發(fā)性大氣污染事故時(shí)有發(fā)生，對(duì)大氣污染擴(kuò)散進(jìn)行模擬和分析，有利于減小事故的危害，減輕人員傷亡和財(cái)產(chǎn)損失。高斯擴(kuò)散模型是國(guó)際原子能機(jī)構(gòu)(IAEA)推薦使用于重氣云擴(kuò)散模擬的數(shù)學(xué)模型，該模型在非重氣云擴(kuò)散的應(yīng)用日益廣泛。高斯擴(kuò)散模型是描述大氣對(duì)有 害氣體的輸移、擴(kuò)散和稀釋作用的物理或數(shù)學(xué)模型，是進(jìn)行災(zāi)害預(yù)測(cè)和救援指揮的有力手段之一。高斯擴(kuò)散模型高斯模型又分為高斯煙團(tuán)模型和高斯煙羽模型。大氣污染物泄漏分為瞬時(shí)泄漏和連續(xù)泄漏，瞬時(shí)泄漏是指污染物泄放的時(shí)

2、間相對(duì)于污染物擴(kuò)散的時(shí)間較短如突發(fā)泄漏等的情形，連續(xù)泄漏則是指污染物泄放的時(shí)間較長(zhǎng)的情形。瞬時(shí)泄漏采用高斯煙團(tuán)模型模擬，而連續(xù)泄漏采用高斯模型煙羽模型模擬。高斯模型適用于非重氣云氣體，包括輕氣云和中性氣云氣體。 要求氣體在擴(kuò)散過(guò)程中，風(fēng)速均勻穩(wěn)定。在高斯煙團(tuán)模型中，選擇風(fēng)向建立坐標(biāo)系統(tǒng)，即取泄漏源為坐標(biāo)原點(diǎn)，x軸指向風(fēng)向，y軸表示在水平面內(nèi)與風(fēng)向垂直的方向，z軸則指向與水平面垂直的方向，具體公式見(jiàn)式(9.1):(x ut)2C(x,y,z, t)y2e花(z_H)2(z H)2,(eP +eF)(9.1)其中：|C(x, y,z,t) |為泄漏介質(zhì)在某位置某時(shí)刻的濃度值；Q為污染物單位時(shí)間排放

3、量(mg/s);匕、已、二分別x、y、z軸上的擴(kuò)散系數(shù)，需根據(jù)大氣穩(wěn)定度選擇參數(shù)計(jì)算得到(m); x、y、z表示x、y、z上的坐標(biāo)值(m); u表示平均風(fēng)速(m/s); t表示擴(kuò)散時(shí)間(s); H 表示泄漏源的高度(m)。同理，高斯煙羽模型的表達(dá)式如：C(x,y,z,t)= Q2 y 2占 e2二 u二 z(z H)2(z H )2(9.2)技術(shù)方法若用高斯模型算出空間每一個(gè)點(diǎn)在一個(gè)時(shí)刻的污染濃度，這個(gè)計(jì)算量是很大的。因此所設(shè)計(jì)的系統(tǒng)一般都是采用先進(jìn)行圖層網(wǎng)格化，由高斯模型計(jì)算出有限個(gè)網(wǎng)格點(diǎn)的上的污染物濃度，在進(jìn)行空間內(nèi)插得到面上每一個(gè)點(diǎn)的污染物濃度，并由此得到污染物濃度的等值線。整個(gè)過(guò)程的示

4、意圖如圖 所示幵始I加載地圖圖921大氣污染擴(kuò)散的高斯模擬的步驟1）圖層網(wǎng)格化圖層網(wǎng)格格式分為結(jié)構(gòu)化網(wǎng)格、非結(jié)構(gòu)化網(wǎng)格。結(jié)構(gòu)化網(wǎng)格是指網(wǎng)格中每個(gè)結(jié)點(diǎn)都有數(shù)量相同的相鄰點(diǎn)，如正方形格網(wǎng)，而非結(jié)構(gòu)化網(wǎng)格則不同。由于結(jié)構(gòu)化網(wǎng)格易于實(shí)現(xiàn)，便于進(jìn)行插值處理，故多用于實(shí)際應(yīng)用。 考慮到氣體污染物質(zhì)量濃度的空間變化頻繁，采用固定的結(jié)構(gòu)化網(wǎng)格，以事故發(fā)生中心地為整個(gè)區(qū)域的幾何中心，采用等間距條件將圖層的二維空間離散化，計(jì)算每個(gè)網(wǎng)格點(diǎn)上面的污染物濃度，自動(dòng)生成反映大氣污染物質(zhì)量濃度分布的等間距網(wǎng)格。2）空間插值空間插值常用于將離散點(diǎn)的測(cè)量數(shù)據(jù)轉(zhuǎn)換為連續(xù)的數(shù)據(jù)曲面，以便與其他空間現(xiàn)象的分布模式進(jìn)行比較，它包括了空

5、間內(nèi)插和外推兩種算法?？臻g插值的理論假設(shè)是空間位置上越靠近的點(diǎn)，越可能具有相似的特征值；而距離越遠(yuǎn)的點(diǎn)，其特征值相似的可能性越小?？臻g插值的方法可以分為整體插值和局部插值兩類。整體插值方法用研究區(qū)所有采樣點(diǎn)的數(shù)據(jù)進(jìn)行全區(qū)特征的擬合，局部插值方法僅僅用鄰近的數(shù)據(jù)點(diǎn)來(lái)估計(jì)未知的值。整體插值方法一般包括邊界內(nèi)插方法、趨勢(shì)面分析和變化函數(shù)插值。局部插值方法只使用臨近的數(shù)據(jù)點(diǎn)來(lái)估計(jì)未知點(diǎn)的數(shù)值，一般包括最鄰近點(diǎn)法（泰森多邊形法 ） 、移動(dòng)平均插值法 （距離倒數(shù)插值法 ）、樣條函數(shù)插值法和空間自協(xié)方差最佳法（克里金插值法 ）。在 ArcGIS 中，有如下插值方法可供使用：全局多項(xiàng)式插值法，反距離權(quán)重法，含

6、障礙 的擴(kuò)散插值法， 含障礙的核插值法， 局部多項(xiàng)式插值法， 徑向基函數(shù)插值法和移動(dòng)窗口克里 金法。3） 系統(tǒng)實(shí)現(xiàn)由于ArcGIS等軟件無(wú)法直接實(shí)現(xiàn)大氣污染物的高斯模擬，需要借助諸如MapX、ArcGISEngine的二次開(kāi)發(fā)實(shí)現(xiàn)氣體濃度值在 GIS平臺(tái)上動(dòng)態(tài)展示的技術(shù)方法。這里以基于 C#.NET 和GIS組件ArcGIS Engine的二次開(kāi)發(fā)為例，介紹實(shí)現(xiàn)二三維模擬大氣污染物高斯擴(kuò)散的方案。ArcGIS Engine是用于構(gòu)建定制應(yīng)用的一個(gè)完整的嵌入式 GIS組件庫(kù)。本案例模擬共包括 五個(gè)功能模塊： 空間數(shù)據(jù)庫(kù)模塊、地圖基本操作模塊、事故參數(shù)設(shè)置模塊、 事故模擬分析模 塊和文檔保存輸出模

7、塊。各個(gè)模塊的主要功能具體如下：a. 空間數(shù)據(jù)庫(kù)模塊主要是指地圖中各個(gè)圖層中的空間數(shù)據(jù)庫(kù)，包含安全數(shù)據(jù)、 加氣站、消火栓、避難場(chǎng)所等應(yīng)急設(shè)施信息以及道路、鐵路建筑物，居民點(diǎn)、橋梁、水系、湖泊等基礎(chǔ)數(shù)據(jù)。數(shù)據(jù)庫(kù)采用OraclelOg作為有關(guān)數(shù)據(jù)庫(kù)的開(kāi)發(fā)與管理工具，在ArcGIS中通過(guò)Add Data添加數(shù)據(jù)庫(kù)中相關(guān)數(shù)據(jù)，并按一定的圖層順序存放，形成 .mxd 格式文檔。b. 地圖基本操作模塊中，可以對(duì)地圖進(jìn)行打開(kāi)圖形、保存圖形、漫游、放大、縮小、 全屏顯示、地圖窗口刷新、屬性選擇、距離量算、面積量算等地理信息系統(tǒng)中的基本操作， 在程序中是通過(guò) ArcGIS Engine組件的IToolbarCo

8、ntrol接口中AddItem 來(lái)實(shí)現(xiàn)。c. 事故參數(shù)設(shè)置模塊， 包括事發(fā)地點(diǎn)查詢和定位、 事發(fā)時(shí)間參數(shù)設(shè)置、 泄漏參數(shù)設(shè)置、 氣象參數(shù)設(shè)置等信息。d. 事故模擬分析模塊， 是在高斯煙羽模型和高斯煙團(tuán)模型的基礎(chǔ)上， 利用 ArcGIS Engine組件在電子地圖上模擬毒氣連續(xù)泄漏或持續(xù)泄漏的擴(kuò)散過(guò)程，主要功能有動(dòng)態(tài)模擬擴(kuò)散過(guò)程模擬、整個(gè)擴(kuò)散過(guò)程的軌跡回放、 某個(gè)時(shí)刻或者某個(gè)時(shí)間段氣體的擴(kuò)散范圍并分析該時(shí)刻或 者時(shí)間段受影響單位的統(tǒng)計(jì)信息。e. 三維分析模塊，利用ArcGIS 3D Analyst擴(kuò)展模塊提供的三維數(shù)據(jù)查看環(huán)境實(shí)現(xiàn)系統(tǒng)的三維可視化。在 ArcSce ne中導(dǎo)入城市三維數(shù)據(jù)，通過(guò)3

9、D View Tools工具條實(shí)現(xiàn)對(duì)圖層的放大、縮小、旋轉(zhuǎn)等操作。調(diào)用IScene-graph的Locate（）方法實(shí)現(xiàn)建筑物的查詢。加載三維大氣擴(kuò)散表面圖層，顯示 TIN 要素中的三角面，設(shè)置要素圖層的顏色、透明度等屬性，顯示大 氣擴(kuò)散的三維區(qū)域。f. 文檔保存輸出模塊， 是指將事故基本信息、 警戒區(qū)域的擴(kuò)散圖形、 受影響的重要單位、城鄉(xiāng)道路、居民建筑、水域橋梁、行政區(qū)域、消防監(jiān)控等信息保存輸出在Word 文檔，為應(yīng)急預(yù)案提供信息支持。 實(shí)例模擬假設(shè)一輛裝載液氯的槽罐車在城市郊區(qū)發(fā)生了車禍， 造成了氯氣泄露事故。 事故發(fā)生時(shí) 間為2012年07月01日12時(shí)，氣體泄漏方式為瞬時(shí)泄漏，泄漏氣體總量為5 000 kg,云量為 3 成，晴天， 3 級(jí)西北風(fēng)。將上述信息輸入模擬系統(tǒng), 包括時(shí)間參數(shù)、 泄露參數(shù)和天氣參數(shù)。 設(shè)置參數(shù)完畢后進(jìn)入 計(jì)算分析過(guò)程, 將某時(shí)刻的計(jì)算結(jié)果導(dǎo)入基礎(chǔ)地理圖中進(jìn)行污染影響范圍繪制。設(shè)定事故發(fā)生一小時(shí)后,對(duì)事故影響的范圍進(jìn)行模擬,得到一小時(shí)后影響范圍