王正維，陳迎輝，楚國華

（1.河北省張家口地質(zhì)環(huán)境監(jiān)測站，河北張家口 075000；2.石家莊經(jīng)濟(jì)學(xué)院，河北石家莊 050031）

0 引言

自2003年6月中國氣象局和國土資源部合作開展全國地質(zhì)災(zāi)害氣象預(yù)警工作以來，各省市陸續(xù)開展了此項(xiàng)工作，在很大程度上減少了人員傷亡和財(cái)產(chǎn)損失，為防災(zāi)減災(zāi)工作開辟了一條新的途徑，已成為保護(hù)廣大人民群眾生命財(cái)產(chǎn)安全的重要手段。但地質(zhì)災(zāi)害氣象預(yù)警工作的理論基礎(chǔ)和技術(shù)方法仍在進(jìn)一步研究完善過程中。在理論研究方面劉傳正等對我國地質(zhì)災(zāi)害的預(yù)警方法開展了深入系統(tǒng)的研究，提出了全面考慮地質(zhì)災(zāi)害潛勢度、危險(xiǎn)度和危害度的“三度”概念［1、2］。目前國內(nèi)開展此項(xiàng)工作的主要技術(shù)方法有考慮前期過程雨量的綜合判據(jù)圖法、網(wǎng)格剖分概率量化綜合預(yù)警模型法、考慮降雨特征量的預(yù)警指標(biāo)法以及基于GIS和神經(jīng)網(wǎng)絡(luò)的地質(zhì)災(zāi)害預(yù)警系統(tǒng)等［3］。由于地質(zhì)災(zāi)害本身的復(fù)雜性，預(yù)警理論方法還處于研究探索階段，各種方法應(yīng)根據(jù)實(shí)際應(yīng)用效果加以調(diào)整。本文以河北省張家口市為例，以改進(jìn)的綜合判據(jù)圖法為理論依據(jù)，利用Delaunay三角網(wǎng)和邊界限定自動(dòng)剖分三角網(wǎng)進(jìn)行地質(zhì)災(zāi)害預(yù)警的自動(dòng)計(jì)算、判別，并開發(fā)了預(yù)警軟件。

1 張家口市地質(zhì)災(zāi)害氣象預(yù)警工作現(xiàn)狀

張家口市位于河北省西北部，面積3.7×104km2，既有山區(qū)、盆地，又有高原、丘陵，地貌及地質(zhì)構(gòu)造復(fù)雜，地質(zhì)環(huán)境脆弱，地質(zhì)災(zāi)害具有種類多、分布范圍廣、突發(fā)性強(qiáng)、危害性大等特點(diǎn)，是地質(zhì)災(zāi)害多發(fā)區(qū)。截止2008年底，已發(fā)現(xiàn)各類地質(zhì)災(zāi)害及隱患708處，主要威脅著村莊、居民、農(nóng)田、礦山、水庫、學(xué)校、交通干線、鄉(xiāng)間公路和水利設(shè)施等，在很大程度上影響著經(jīng)濟(jì)建設(shè)的快速發(fā)展。因此，在張家口市開展地質(zhì)災(zāi)害氣象預(yù)警工作意義重大。

從2005年6月開始，張家口市國土資源局和張家口市氣象局聯(lián)合開展了汛期地質(zhì)災(zāi)害氣象預(yù)警工作。市氣象臺(tái)提供雨量實(shí)測和預(yù)測資料，張家口市地質(zhì)環(huán)境監(jiān)測站根據(jù)氣象臺(tái)提供的實(shí)測預(yù)測雨量結(jié)合相關(guān)地質(zhì)資料進(jìn)行綜合研究分析，在與各有關(guān)方面充分會(huì)商后做出地質(zhì)災(zāi)害預(yù)警預(yù)報(bào)。根據(jù)掌握和收集的大量資料，將全市劃分了16個(gè)地質(zhì)災(zāi)害預(yù)警區(qū)。根據(jù)已經(jīng)發(fā)生的地質(zhì)災(zāi)害單體或群體前期累計(jì)過程雨量，參考相鄰區(qū)域和地質(zhì)環(huán)境條件相似區(qū)域的預(yù)警判據(jù)，確定了各個(gè)預(yù)警區(qū)的地質(zhì)災(zāi)害發(fā)生的臨界過程雨量判據(jù)。在2005～2006年的地質(zhì)災(zāi)害氣象預(yù)警工作中，采用各雨量站點(diǎn)預(yù)警日前期過程雨量累加，與雨量站點(diǎn)所在地質(zhì)災(zāi)害預(yù)警區(qū)預(yù)警判據(jù)對比，分析判斷后得出地質(zhì)災(zāi)害預(yù)警級別。由于采用的網(wǎng)絡(luò)數(shù)據(jù)傳輸、計(jì)算機(jī)軟件輔助計(jì)算、預(yù)警級別人工判別方法，速度慢、效率低，計(jì)算和判別過程誤差大、自動(dòng)化程度低，因此，需要開發(fā)一套預(yù)警全過程自動(dòng)計(jì)算判別軟件，實(shí)現(xiàn)預(yù)警過程的自動(dòng)化，以提高計(jì)算速度、準(zhǔn)確度，消除人工判別誤差，提升地質(zhì)災(zāi)害氣象預(yù)警工作效率。

2 預(yù)警理論方法

張家口市地質(zhì)災(zāi)害方面研究程度低，各方面資料較少，因此，實(shí)現(xiàn)按網(wǎng)格剖分的地質(zhì)災(zāi)害潛勢度預(yù)警方法比較困難，所以，考慮當(dāng)日雨量與累積雨量的綜合判別表方法作為預(yù)警的基礎(chǔ)理論依據(jù)［4］。根據(jù)以往的研究成果，將張家口市劃分為16個(gè)地質(zhì)災(zāi)害預(yù)警區(qū)，進(jìn)而確定了各個(gè)預(yù)警區(qū)地質(zhì)災(zāi)害發(fā)生的臨界雨量判據(jù)。

目前，國內(nèi)從國家級到省級的地質(zhì)災(zāi)害氣象預(yù)警業(yè)務(wù)普遍采用前期累計(jì)雨量、有效累計(jì)雨量與最新日雨量等預(yù)警指標(biāo)，張家口市地質(zhì)災(zāi)害氣象預(yù)警分析雨量采用預(yù)警日當(dāng)日預(yù)測雨量和2、4、7、10、15日前期累計(jì)雨量進(jìn)行綜合預(yù)警分析。按照國土資源部和中國氣象局的聯(lián)合規(guī)定，地質(zhì)災(zāi)害氣象預(yù)警分為5級。

國內(nèi)最新研究成果表明［5-7］，地質(zhì)災(zāi)害發(fā)生的規(guī)模和類型與降雨量有明顯的正相關(guān)，而不同降雨類型下災(zāi)害的發(fā)生情況也存在明顯差異，地質(zhì)災(zāi)害的發(fā)生存在明顯的滯后性，一般集中在降雨后72h內(nèi)發(fā)生，從理論和國內(nèi)有關(guān)觀測事實(shí)來說，地質(zhì)災(zāi)害發(fā)生的當(dāng)日激發(fā)雨量對地質(zhì)災(zāi)害的發(fā)生有著重要影響，5天以內(nèi)降水情況對地質(zhì)災(zāi)害的影響隨時(shí)間延續(xù)而明顯減弱，超過5天以上的前期降雨情況對地質(zhì)災(zāi)害發(fā)生的影響很小。因此，作者采用的6個(gè)時(shí)間段根據(jù)預(yù)警判據(jù)計(jì)算的預(yù)警結(jié)果對預(yù)警日綜合預(yù)警結(jié)果取值影響程度各不相同，距預(yù)警日時(shí)間越長，影響越小。本文使用權(quán)重系數(shù)來表示前期6個(gè)時(shí)間段的預(yù)警計(jì)算結(jié)果對綜合預(yù)警結(jié)果的影響程度，將計(jì)算取得的各預(yù)警級別各時(shí)間段的權(quán)重系數(shù)之和稱為預(yù)警指數(shù)。只有當(dāng)預(yù)警指數(shù)達(dá)到或超過某一設(shè)定值時(shí)，即達(dá)到預(yù)警指數(shù)臨界值，該預(yù)警級別計(jì)算結(jié)果才有效，這個(gè)計(jì)算過程可以有效避免以往綜合預(yù)警級別取各時(shí)間段預(yù)警結(jié)果極大值時(shí)距離預(yù)警日較遠(yuǎn)時(shí)間段的降雨情況對綜合預(yù)警結(jié)果的過大影響，體現(xiàn)距離預(yù)警日越近的降雨情況對綜合預(yù)警結(jié)果影響越大的特點(diǎn)。經(jīng)過對相關(guān)資料的研究成果的分析，初步為6個(gè)時(shí)間段預(yù)警結(jié)果分別設(shè)定不同的權(quán)重系數(shù)，見表1：

綜合預(yù)警結(jié)果取值時(shí)，按預(yù)警級別將6個(gè)時(shí)間段的預(yù)警計(jì)算結(jié)果權(quán)重進(jìn)行相加，如果6個(gè)時(shí)間段預(yù)警結(jié)果有一項(xiàng)以上為某一預(yù)警級別時(shí)，則要在該級別的權(quán)重之和上加上基本權(quán)重，最后取各個(gè)預(yù)警級別中預(yù)警指數(shù)（即總權(quán)重）達(dá)到預(yù)警指數(shù)臨界值的級別中最高級別為預(yù)警日綜合預(yù)警級別。根據(jù)6個(gè)時(shí)間段預(yù)警結(jié)果權(quán)重和預(yù)警指數(shù)臨界值取值可以看出，當(dāng)日計(jì)算結(jié)果對綜合預(yù)警級別取值起決定性作用，5日以后時(shí)段預(yù)警結(jié)果對綜合預(yù)警級別取值作用很小。

表1 預(yù)警指數(shù)臨界值與權(quán)重系數(shù)Table 1 Critical value of early warning index and weight coefficients

Delaunay［8］三角形網(wǎng)格具有算法簡單、生成容易、有較強(qiáng)的控制性、單元具有較好的形狀和較強(qiáng)的邊界適應(yīng)性，在二維空間計(jì)算方面能夠表現(xiàn)出良好的性能，因此，作者在對張家口市地質(zhì)災(zāi)害氣象預(yù)警采用Delaunay三角形網(wǎng)絡(luò)剖分方法，充分利用Delaunay三角形網(wǎng)絡(luò)剖分方法，充分利用Delaunay三角形網(wǎng)絡(luò)的性能優(yōu)勢進(jìn)行等值線一成和圖形計(jì)算。

以張家口市界為邊界限定條件，對預(yù)警區(qū)域進(jìn)行Delaunay三角形網(wǎng)格邊界限定自動(dòng)剖分，剖分取得的等邊三角形網(wǎng)絡(luò)與各預(yù)警區(qū)多邊形進(jìn)行區(qū)對區(qū)空間判別計(jì)算，用自動(dòng)剖分三角形的形心代表整個(gè)三角形區(qū)域。通過上述計(jì)算，將預(yù)警區(qū)的數(shù)據(jù)信息轉(zhuǎn)換到自動(dòng)剖分三角形的形心上，可以得出每一個(gè)計(jì)算三角形的6個(gè)時(shí)間段雨量預(yù)警判據(jù)。

對張家口市域內(nèi)的210個(gè)自動(dòng)雨量站點(diǎn)進(jìn)行凸域點(diǎn)集Delaunay三角形網(wǎng)格自動(dòng)構(gòu)網(wǎng)，建立預(yù)警區(qū)的雨量站點(diǎn)三角網(wǎng)。在具體計(jì)算時(shí)，從雨量數(shù)據(jù)庫中取出前期過程雨量，利用三角網(wǎng)的網(wǎng)對網(wǎng)計(jì)算進(jìn)行預(yù)警三角網(wǎng)空間插值，計(jì)算出每一個(gè)預(yù)警三角形的過程雨量，進(jìn)而計(jì)算出各時(shí)間段預(yù)警級別和綜合預(yù)警級別。

對預(yù)警日綜合預(yù)警級別計(jì)算完成后，預(yù)警結(jié)果按等級進(jìn)行三角形區(qū)域拓?fù)淙诤嫌?jì)算，確定各個(gè)預(yù)警區(qū)域級別、面積、邊界后，進(jìn)行計(jì)算機(jī)輸出和保存［9］。

3 預(yù)警方法的實(shí)現(xiàn)與軟件開發(fā)

3.1 雨量數(shù)據(jù)庫的設(shè)計(jì)與實(shí)現(xiàn)

雨量數(shù)據(jù)庫用于存儲(chǔ)各自動(dòng)雨量站點(diǎn)的實(shí)測與預(yù)測雨量，將其按日期順序設(shè)計(jì)為固定長度的數(shù)據(jù)存儲(chǔ)區(qū)塊，以預(yù)警計(jì)算日為關(guān)鍵字段，實(shí)現(xiàn)雨量數(shù)據(jù)按日期快速隨機(jī)定位、查詢、存取，數(shù)據(jù)庫設(shè)計(jì)為固定長度，設(shè)置雨量數(shù)據(jù)初始值，可以在讀取過程雨量時(shí)判斷數(shù)據(jù)存儲(chǔ)狀態(tài)。

3.2 預(yù)警區(qū)的建立

將張家口市現(xiàn)有16個(gè)預(yù)警區(qū)作為基本預(yù)警區(qū)（圖1中實(shí)線表示），基本預(yù)警區(qū)的外部包容邊界為市界，為了便于進(jìn)行市界邊界處的自動(dòng)剖分三角網(wǎng)對預(yù)警區(qū)的區(qū)對區(qū)判別計(jì)算，將市界外側(cè)一定距離范圍內(nèi)的區(qū)域作為擴(kuò)展預(yù)警區(qū)，擴(kuò)展預(yù)警區(qū)各項(xiàng)數(shù)據(jù)與相鄰的基本預(yù)警區(qū)相同。

圖1 預(yù)警區(qū)與雨量站點(diǎn)三角網(wǎng)Fig.1 Early warning area and triangular networking of rainfall stands

3.3 雨量站點(diǎn)三角網(wǎng)的建立

將全市232個(gè)自動(dòng)雨量站點(diǎn)按凸域二維點(diǎn)集合的形式進(jìn)行Delaunay三角形自動(dòng)構(gòu)網(wǎng)（圖1中虛線表示）。利用雨量站點(diǎn)三角網(wǎng)可以實(shí)現(xiàn)對自動(dòng)剖分三角網(wǎng)的三角形雨量插值。

3.4 自動(dòng)剖分三角網(wǎng)的實(shí)現(xiàn)

考慮到計(jì)算機(jī)性能和實(shí)際預(yù)警精度需要，設(shè)定三角形邊長為2.5km，每個(gè)三角形的面積2.7km2，以張家口市界為限定邊界，建立自動(dòng)剖分等邊三角網(wǎng)，對所有自動(dòng)生成的三角形進(jìn)行空間多邊形（三角形）對市界多邊形自動(dòng)檢測判別，保證所得三角形與基本預(yù)警區(qū)的相關(guān)性，實(shí)際運(yùn)算得到1.6萬多個(gè)三角形。

圖2 程序設(shè)計(jì)主要流程圖Fig.2 Process of the program design

3.5 預(yù)警三角網(wǎng)的實(shí)現(xiàn)

預(yù)警三角網(wǎng)是在自動(dòng)剖分三角網(wǎng)的基礎(chǔ)上構(gòu)建的，將自動(dòng)剖分三角網(wǎng)與所有預(yù)警區(qū)多邊形進(jìn)行區(qū)對區(qū)相交判別分析，得到預(yù)警三角網(wǎng)。設(shè)計(jì)預(yù)警三角網(wǎng)的數(shù)據(jù)結(jié)構(gòu)不存儲(chǔ)自動(dòng)剖分三角網(wǎng)數(shù)據(jù)項(xiàng)中的位置與空間關(guān)系數(shù)據(jù)，而是建立自動(dòng)剖分三角形與預(yù)警三角形的相互指針鏈接，可以快速實(shí)現(xiàn)兩個(gè)三角網(wǎng)中相對應(yīng)三角形的數(shù)據(jù)交換，減少內(nèi)存占用，提高程序運(yùn)行速度。

3.6 軟件設(shè)計(jì)的主要流程

軟件開發(fā)的基本思路是以雨量數(shù)據(jù)庫、雨量站點(diǎn)三角網(wǎng)和地質(zhì)災(zāi)害預(yù)警三角網(wǎng)為基礎(chǔ)，雨量數(shù)據(jù)存儲(chǔ)于雨量數(shù)據(jù)庫中，利用雨量站點(diǎn)三角網(wǎng)計(jì)算出各節(jié)點(diǎn)各時(shí)段雨量值，進(jìn)一步自動(dòng)生成等值線；利用雨量站點(diǎn)三角網(wǎng)對地質(zhì)災(zāi)害預(yù)警三角網(wǎng)形心的自動(dòng)插值，計(jì)算出每個(gè)地質(zhì)災(zāi)害預(yù)警三角形各時(shí)段的累計(jì)雨量值，通過與預(yù)警判據(jù)對比計(jì)算出各時(shí)段預(yù)警級別，進(jìn)而計(jì)算出綜合預(yù)警級別，軟件設(shè)計(jì)的主要流程見圖2。

3.7 軟件開發(fā)與操作

以Microsoft Visual C++6.0為開發(fā)工具，借助MFC類庫的強(qiáng)大功能進(jìn)行軟件開發(fā)，為加快運(yùn)行速度，提高程序靈活性，擺脫商品軟件二次開發(fā)的使用限制，全部代碼由軟件自身實(shí)現(xiàn)，數(shù)據(jù)庫格式使用自定義格式。

圖3 程序主要功能模塊Fig.3 Main function modules of the program

地質(zhì)災(zāi)害氣象預(yù)警程序由建立預(yù)警模型、導(dǎo)入雨量數(shù)據(jù)、預(yù)警計(jì)算和顯示輸出四個(gè)功能模塊組成（圖3），建立預(yù)警模型和導(dǎo)入雨量數(shù)據(jù)是地質(zhì)災(zāi)害預(yù)警的基礎(chǔ)，預(yù)警計(jì)算模塊實(shí)現(xiàn)了該程序的核心功能顯示保存模塊實(shí)現(xiàn)了預(yù)警計(jì)算結(jié)果的對外輸出保存功能。

圖4 2008年7月29日實(shí)測15日降雨等值線Fig.4 Contours of 15-day rainfall measured on July 29，2008

軟件操作主要分為五類：基礎(chǔ)數(shù)據(jù)庫操作、雨量數(shù)據(jù)庫操作、預(yù)警計(jì)算、預(yù)警結(jié)果的顯示存儲(chǔ)、輔助工具選項(xiàng)。

基礎(chǔ)數(shù)據(jù)庫操作主要包括雨量數(shù)據(jù)庫的建立、預(yù)警區(qū)數(shù)據(jù)文件的建立、邊界限定自動(dòng)剖分三角網(wǎng)的生成、預(yù)警三角網(wǎng)的構(gòu)建、雨量站點(diǎn)三角網(wǎng)的構(gòu)建。通過菜單操作，控制參數(shù)從文本文件輸入，只要輸入相關(guān)數(shù)據(jù)文件名稱，便可完成操作。雨量數(shù)據(jù)庫操作主要包括成批導(dǎo)入實(shí)測雨量、成批導(dǎo)入預(yù)測雨量、導(dǎo)入實(shí)測雨量和導(dǎo)入預(yù)測雨量，這此操作將文本文件的數(shù)據(jù)導(dǎo)入雨量數(shù)據(jù)庫中。

圖5 2008年8月8日預(yù)測10日降雨等值線Fig.5 Contours of 10-day rainfall estimated on August 8，2008

預(yù)警計(jì)算部分可以實(shí)現(xiàn)指定預(yù)警日期的地質(zhì)災(zāi)害自動(dòng)預(yù)警，還可以實(shí)現(xiàn)歷史日期的預(yù)警結(jié)果顯示。

顯示編輯和預(yù)警/等值線部分主要是對預(yù)警計(jì)算結(jié)果的顯示輸出，具有放大、縮小、原始尺寸顯示、設(shè)置縮放系數(shù)、存儲(chǔ)圖像、顯示圖例等功能，顯示內(nèi)容可以選擇6個(gè)時(shí)間段的實(shí)測或預(yù)測雨量等值線、預(yù)警結(jié)果，可以同時(shí)選擇顯示基礎(chǔ)界線、雨量站點(diǎn)三角網(wǎng)、預(yù)警三角網(wǎng)。顯示的內(nèi)容均可以保存為JPEG格式的圖像文件。

所有的參數(shù)可以通過配置文件修改，軟件設(shè)計(jì)可以適應(yīng)不同地區(qū)、不同條件的地質(zhì)災(zāi)害預(yù)警工作。

4 軟件應(yīng)用

使用該軟件系統(tǒng)進(jìn)行了張家口市2007年和2008年各136天的汛期地質(zhì)災(zāi)害氣象預(yù)警，軟件運(yùn)行正常，預(yù)警結(jié)果與人工計(jì)算預(yù)測結(jié)果基本相同，共預(yù)警9次，其中氣象部門實(shí)測降雨量與預(yù)測降雨量基本相同的6次均成功預(yù)警，避免了人員傷亡和經(jīng)濟(jì)損失。在自動(dòng)剖分三角形邊長2.5km、三角形數(shù)量1.6萬多個(gè)的條件下，使用P4 3.0GHz雙核處理器、512M內(nèi)存的計(jì)算機(jī)進(jìn)行計(jì)算，預(yù)警計(jì)算時(shí)占用內(nèi)存43M，計(jì)算時(shí)間1分05秒。實(shí)踐證明，該軟件可以滿足市級地質(zhì)災(zāi)害預(yù)警工作的時(shí)間要求和精度要求（圖4～圖7）。

圖6 2008年9月2日預(yù)警等級計(jì)算結(jié)果Fig.6 Rand of forewarning by calculation on September 2，2008

圖7 2008年9月8日預(yù)警等級計(jì)算結(jié)果Fig.7 Rank of forewarning by calculation on September 8，2008

5 結(jié)論

本文在地質(zhì)災(zāi)害氣象預(yù)警的綜合判據(jù)圖法的基礎(chǔ)上，對地質(zhì)災(zāi)害預(yù)警方法進(jìn)行了改進(jìn)，對不同時(shí)段按綜合判據(jù)圖法計(jì)算的預(yù)警等級給予不同的權(quán)重系數(shù)，并引入預(yù)警指數(shù)和預(yù)警指數(shù)臨界值進(jìn)行綜合預(yù)警等級計(jì)算，體現(xiàn)了距離預(yù)警日越近的降雨情況對綜合預(yù)警結(jié)果影響越大的特點(diǎn)。

使用Delaunay三角網(wǎng)理論方法進(jìn)行地質(zhì)災(zāi)害預(yù)警區(qū)邊界限定自動(dòng)剖分、構(gòu)建雨量站點(diǎn)三角網(wǎng)，利用Delaunay三角網(wǎng)的特性進(jìn)行雨量自動(dòng)插值、等值線追蹤、預(yù)警結(jié)果計(jì)算和輸出，開發(fā)了預(yù)警軟件，采用與數(shù)據(jù)存儲(chǔ)結(jié)構(gòu)相適應(yīng)的優(yōu)化算法，提高了程序運(yùn)行速度，實(shí)現(xiàn)了各項(xiàng)基礎(chǔ)數(shù)據(jù)庫建立、雨量數(shù)據(jù)庫存取、預(yù)警計(jì)算全過程自動(dòng)化，參數(shù)使用配置文件讀取，軟件模型應(yīng)用靈活、適應(yīng)性好，在基礎(chǔ)資料充分的情況下、軟件可以擴(kuò)展構(gòu)建成三角形單元網(wǎng)格剖分概率量化綜合預(yù)警模型。

軟件未進(jìn)行不同地區(qū)和不同條件下的測試，各時(shí)間段的預(yù)警結(jié)果權(quán)重系數(shù)、預(yù)警指數(shù)、預(yù)警指數(shù)臨界值的取值方法和大小還需要在應(yīng)用中進(jìn)一步探討。

受氣象部門降雨量預(yù)報(bào)精度和準(zhǔn)確度限制，預(yù)警日預(yù)測雨量只能取預(yù)報(bào)的可能降雨量值的最大值，在一定程度上降低了地質(zhì)災(zāi)害氣象預(yù)警精度和準(zhǔn)確度。在得到氣象部門等相關(guān)方面技術(shù)支持的條件下，軟件還可以實(shí)現(xiàn)自動(dòng)導(dǎo)入各雨量站點(diǎn)實(shí)測和預(yù)測雨量、預(yù)警結(jié)果網(wǎng)絡(luò)傳輸?shù)裙δ?，進(jìn)一步提高地質(zhì)災(zāi)害氣象預(yù)警工作的效率和準(zhǔn)確度，更好的為地質(zhì)災(zāi)害防治工作服務(wù)。

［1］劉傳正.地質(zhì)災(zāi)害預(yù)警工程體系探討［J］.水文地質(zhì)工程地質(zhì)，2000，27（4）：1-4.

［2］劉傳正.突發(fā)性地質(zhì)災(zāi)害的監(jiān)測預(yù)警問題［J］.水文地質(zhì)工程地質(zhì)，2001，28（4）：1-4.

［3］李曉.重慶地區(qū)的強(qiáng)降雨過程與地質(zhì)災(zāi)害的相關(guān)分析［J］.中國地質(zhì)災(zāi)害與防治學(xué)報(bào)，1995，6（3）：39-42.

［4］陳百煉，楊勝元，楊森林，等.基于GIS的地質(zhì)災(zāi)害氣象預(yù)警方法初探［J］.中國地質(zhì)災(zāi)害與防治學(xué)報(bào)，2005，16（4）：93-96.

［5］鐘蔭乾.滑坡與降雨關(guān)系及其預(yù)報(bào)［J］.中國地質(zhì)災(zāi)害與防治學(xué)報(bào)，1998，9（4）：81-86.

［6］潘懋，李鐵鋒.災(zāi)害地質(zhì)學(xué)［M］.北京：北京大學(xué)出版社，2002，10-11.

［7］傅鶴林，周中，卜翠松，等.地質(zhì)災(zāi)害預(yù)測預(yù)報(bào)國內(nèi)外現(xiàn)狀［J］.湘南學(xué)院學(xué)報(bào)，2006，27（2）：43-48.

［8］鄔吉明，沈隆均，張景琳.Delaunay三角網(wǎng)格的一種快速生成法［J］.數(shù)值計(jì)算與計(jì)算機(jī)應(yīng)用，2001，22（4）：267-275.

［9］章孝燦，黃智才，戴企成，等.GIS中基于拓?fù)浣Y(jié)構(gòu)和凸殼技術(shù)的快速TIN生成算法［J］.計(jì)算機(jī)學(xué)報(bào)，2002，25（11）：1212-1218.

［10］周培德.計(jì)算幾何——算法分析與設(shè)計(jì)［M］.北京：清華大學(xué)出版社，2000，103-109.

［11］盧朝陽，吳成柯.任意多邊形內(nèi)帶特征約束的散列數(shù)據(jù)的最優(yōu)三角剖分.計(jì)算機(jī)設(shè)計(jì)與圖形學(xué)學(xué)報(bào)，1997，9（4）：302-308.