熊成利，金洪芳，鄭永忠

(浙江省第二測繪院，浙江 杭州 310012)

目前，中國19個省份中超過50個城市發(fā)生了不同程度的地面沉降，重災區(qū)主要為長江三角洲地區(qū)、華北平原和汾渭盆地這3個區(qū)域[1]。地面沉降作為一種累進性的地質災害,影響范圍廣,防治難度大,已成為困擾城市建設、影響居民安全生活的一種重要地質災害，并且嚴重影響著城市地區(qū)的經(jīng)濟增長和可持續(xù)發(fā)展[2]。為了有效地控制城市地面沉降，杭嘉湖地區(qū)沿?？h市政府城市規(guī)劃、地質等管理部門每年都投入大量人力和財力對地面沉降及地下水位進行不間斷的動態(tài)監(jiān)測，并積累了大量豐富的觀測數(shù)據(jù)，為城市規(guī)劃和建設發(fā)展提供了重要的科學依據(jù)。然而由于歷史條件或技術原因的限制，這些大量的寶貴資料都是通過手工整理，工作效率低，查閱資料也極為不便，大批的數(shù)據(jù)沒有得到科學的管理、分析和有效的利用。如何有效管理這些數(shù)據(jù)，并進行合理分析，進而掌握沿?？h市的地面沉降變化規(guī)律,分析地面沉降機理，對防震減災、城市規(guī)劃建設等都具有十分重要的意義。

1 需求分析

為科學地管理多年歷史沉降監(jiān)測資料，建立一套合理科學的基于GIS的沉降信息管理和共享服務平臺，實現(xiàn)不同時期、不同觀測點的多源多時相沉降觀測數(shù)據(jù)的自動入庫、數(shù)值分析、可視化顯示、趨勢分析與預測、各類報表、圖形的自動生成與輸出，及時快速地掌握地面沉降變化規(guī)律，提供處理措施，為城市規(guī)劃、城市建設與管理及社會公眾信息需求提供服務，提高城市規(guī)劃和城市空間的管理開發(fā)與利用水平，筆者所在單位組織相關技術力量進行了沿?？h市級城區(qū)地面沉降信息管理系統(tǒng)的研究與開發(fā)，并推廣至杭嘉湖多個縣市，取得了很好的用戶評價。

2 系統(tǒng)設計

2.1 系統(tǒng)總體框架

本系統(tǒng)采用B/S三層架構，PLSQL Developer進行Oracle 11g數(shù)據(jù)庫管理，IIS發(fā)布網(wǎng)站，實現(xiàn)對地面沉降數(shù)據(jù)的管理、統(tǒng)計查詢分析和預測分析等。系統(tǒng)的技術框架如圖1所示。

2.2 關鍵技術

2.2.1 數(shù)據(jù)庫設計

基于B/S框架的沉降信息的管理與可視化系統(tǒng)設計中，數(shù)據(jù)庫設計是整個設計中較為重要的部分之一，其設計的好壞對系統(tǒng)性能有很大影響。本系統(tǒng)運用Oracle實現(xiàn)與沉降相關數(shù)據(jù)交換與存儲。與沉降相關的數(shù)據(jù)主要有站點的基本數(shù)據(jù)及站點的多期時間序列沉降數(shù)據(jù)，在該系統(tǒng)的數(shù)據(jù)庫設計中主要分沉降站點基本數(shù)據(jù)表和多期歷史沉降數(shù)據(jù)表，相關E-R如圖2所示。系統(tǒng)根據(jù)該E-R圖進行了沉降數(shù)據(jù)相關表(沉降觀測點表、歷年高程記錄表)、底圖配置表、用戶表、圖片或文件與觀測點對應表等詳細的數(shù)據(jù)表設計。

2.2.2 沉降數(shù)據(jù)預測方法研究

利用數(shù)學模型將多期觀測的地面沉降量數(shù)據(jù)擬合成曲線，進而計算后期沉降的方法在國內(nèi)外得到了廣泛的應用。這類方法包括雙曲線法、三點法、Asaoka法、泊松曲線法、灰色理論和人工神經(jīng)網(wǎng)絡法等。各種預測模型的適用性不同，如雙曲線法適合長期預測，時間越長，準確度越高[3]；而Asaoka法的優(yōu)點在于它可利用較短時間內(nèi)的觀測資料，就能得到較為可靠的預測結果[4]；灰色模型要求原始數(shù)據(jù)變化較為平緩[5]。城市范圍內(nèi)地面觀測點較多，分布范圍廣，地質條件不盡相同；任何一種地面沉降模型都是在一定地質條件下推導出來的，在城市沉降預測分析中，不可能預先為觀測點選擇適用的模型，也不可能設計一種通用的地面沉降模型[6]。因此，沉降預測模型計算出結果后，還需對結果進行精度驗證，可以選用平均絕對百分誤差(MAPE)進行模型精度評價[1]。通過這種精度評價可以得出的結論是：越多的原始序列期數(shù)越能反映觀測點的沉降趨勢；在實際應用中，應優(yōu)先信任短期預測結果；可在觀測期數(shù)5、6期以上時優(yōu)先采用灰度模型；在缺少數(shù)據(jù)樣本的支持下，應當優(yōu)先采用Asaoka法該法進行預測；三點法在當前樣本下比較適合長期預測；雙曲線法適應性較強，其預測精度不高；觀測點在不同預測模型的預測精度也不相同，可以對模型進行改進。

2.2.3 天地圖與沉降數(shù)據(jù)融合

地面沉降信息管理系統(tǒng)包括多期沉降觀測數(shù)據(jù)的整理與和建庫、沉降觀測數(shù)據(jù)的空間可視化與查詢，以及沉降觀測數(shù)據(jù)時空分析，為城市地面沉降的預測、防控、減災提供支持。沉降觀測數(shù)據(jù)通常為單一離散觀測點，因此，對沉降觀測信息的管理、分析和預測需要城市級基礎地理信息數(shù)據(jù)支撐。本系統(tǒng)提出以采用智慧城市地理信息共享平臺架構如縣市天地圖作為城市離散沉降信息管理的空間展示基礎。

2.2.4 基于交叉驗證的空間插值方法

在城市沉降監(jiān)測中，常用的單點測量手段獲取的單站觀測數(shù)據(jù)是離散的、單一的。由于觀測站點個數(shù)受到限制，獲取的觀測資料只能反映局部沉降監(jiān)測信息，無法從整體上對觀測區(qū)進行監(jiān)測。為了彌補這一不足，近年來，多個專家學者利用空間插值的方法進行了地面沉降由點到面的監(jiān)測研究，并在城市地面沉降觀測中得到了較多應用。常見的空間插值方法有：IDW法、Kriging和樣條函數(shù)方法[7]。對于不同的觀測點分布、不同的地形環(huán)境和觀測數(shù)據(jù)量，其插值計算結果在一定范圍內(nèi)存在差異。IDW插值把估算值和實測點間的距離作為權重因子，權重值為距離反比[8-9]；Kriging插值是通過一組具有z值的分散點生成估計表面的高級地統(tǒng)計過程[10-11]；樣條函數(shù)插值是通過控制估計方差，利用一些特征節(jié)點，用多項式擬合的方法來產(chǎn)生平滑的插值曲線[12-14]。在城市地面沉降插值計算時，3種空間插值模型表現(xiàn)出各自的特點[15]:IDW法簡單易行，效率高，但是易受極值影響；Kriging法考慮區(qū)域空間特性，但是算法復雜度高并且耗時；樣條函數(shù)法保凸性、逼真性和平滑性好，但是誤差估計困難要求點數(shù)量足夠且分布均勻。

進行空間插值后需要對插值結果進行精度驗證。首先假定某些觀測點高程值未知，通過周圍站點實測值插值來估算，然后計算驗證站點的實測值與估算值之間的誤差值。被估算高程值的站點稱為驗證站點，用于插值的站點稱為采樣站點。常用的評估方法有平均絕對誤差(MAE)和平均誤差平方的平方根(RMSIE)。采用這兩種評估標準能很好地反映不同插值方法的相對精度，以此選擇最適宜的空間插值方法進行城市沉降計算。

2.2.5 沉降數(shù)據(jù)導出工具設計

導出工具用于將沉降趨勢預測結果轉換為天地圖支持的格式。它是基于Visual Studio 2015集成開發(fā)環(huán)境，采用C#語言，利用二次開發(fā)平臺ArcGIS Engine 10.1進行開發(fā)的。導出工具的地圖瀏覽功能基于桐鄉(xiāng)地圖的MXD格式的地圖文檔；數(shù)據(jù)庫使用沉降信息管理系統(tǒng)的數(shù)據(jù)庫，通過服務器發(fā)布的數(shù)據(jù)庫服務進行通信。導出工具根據(jù)用戶選擇的條件對觀測點進行沉降趨勢預測并將結果以專題圖的形式展現(xiàn)在地圖中，且可以生成相應的圖例。最后可以將生成的展示結果以MXD的格式保存至本地硬盤。其工作流程如圖3所示。

2.3 功能模塊

本系統(tǒng)具體的功能模塊包括：

(1) 觀測點與高程信息管理。存放觀測點的位置信息和高程信息，是本系統(tǒng)進行沉降統(tǒng)計與分析的基礎數(shù)據(jù)。

(2) 文檔管理。存放與觀測點相關的信息，如觀測點的照片信息和點之記等。文檔管理可以按照文檔內(nèi)容的性質對文檔進行分門別類，提高檢索的便捷度。支持上傳、編輯、下載、標注、條件查詢等文件管理操作。

(3) 底圖管理。用戶可以增加、刪除、修改底圖信息，切換使用不同的底圖。

(4) 數(shù)據(jù)導入與導出。數(shù)據(jù)導入可以隨時更新系統(tǒng)基礎數(shù)據(jù)，使得系統(tǒng)能夠存儲更多觀測點的歷年高程信息，導入沉降數(shù)據(jù)時還可以驗證其有效性、正確性；數(shù)據(jù)導出可以選擇一段數(shù)據(jù)，得到觀測點的位置信息、高程信息及沉降統(tǒng)計信息等。

(5) 用戶管理。系統(tǒng)采用的用戶分為一般用戶、數(shù)據(jù)管理員及超級管理員等3種角色。不同類型的用戶分別有不同的權限：一般用戶只有查看預覽與檢索查詢沉降信息的功能；數(shù)據(jù)管理員不但有預覽沉降信息的功能，而且還有增加、編輯、底圖配置、統(tǒng)計分析、預測分析、InSAR數(shù)據(jù)分析等權限；超級管理員具有系統(tǒng)的最高權限，除了有數(shù)據(jù)管理員的權限之外，還有對用戶管理的權限。

(6) 報告導出。用戶可以根據(jù)需要，選擇任一時間段的控制點信息，生成包含有趨勢預測信息、沉降等值線信息、DEM變化信息、淹沒信息、沉降量統(tǒng)計信息的報告并可以打印。

(7) 空間查詢和模糊查詢?？臻g查詢指用戶可以在地圖上框選觀測點，查詢其高程信息和屬性信息。模糊查詢主要是指用戶可以根據(jù)觀測點中名字的關鍵字查詢到觀測點，并在底圖上高亮顯示。

(8) 圖層管理。用戶可以根據(jù)自己的需要選擇哪些觀測點可以出現(xiàn)在底圖上，哪些工具需要被隱藏。

(9) 趨勢預測。將雙曲線法、三點法、Asaoka法及灰度模型等4種方法分別應用于已知觀測年份的沉降數(shù)據(jù)預測未知年份的沉降數(shù)據(jù)，并將預測結果以折線圖或沉降圖的形式展示出來。

(10) 沉降統(tǒng)計。根據(jù)數(shù)據(jù)庫中已經(jīng)存在的觀測點數(shù)據(jù)，對某一時間段或所有觀測點的沉降量進行統(tǒng)計，最后以折線圖和柱狀圖的形式展現(xiàn)出來。

(11) 沉降等值線。根據(jù)已知高程信息，得到沉降等值線，并在地圖上展現(xiàn)出來。

(12) 影像分析。影像分析主要是根據(jù)InSAR數(shù)據(jù)進行分析，分為沉降分析和淹沒分析。將分析的結果進行渲染，并在地圖上顯示出來。

(13) 決策分析。已知某地區(qū)內(nèi)一些點的高程信息和沉降信息，提供一組預備建設的規(guī)劃圖進行檢測(可以是點、線、面)，不同的建設物對其所占區(qū)域的高程信息和沉降信息等有不同的要求，要求對其中不適合進行建設的地方進行標注提示，以輔助最終決策的制定。

3 系統(tǒng)實現(xiàn)與推廣應用

本系統(tǒng)開發(fā)完成后在杭嘉湖沿海地區(qū)產(chǎn)生了積極的響應，如在桐鄉(xiāng)、海寧和海鹽、嘉善地區(qū)都得到了很好的應用。系統(tǒng)中沉降等值線功能、沉降點預測、沉降面預測分別如圖4和圖5所示。

4 結 語

本文利用現(xiàn)代的數(shù)據(jù)庫技術、GIS技術及網(wǎng)絡技術等，在網(wǎng)絡化的計算機硬件基礎上，緊密結合沿?？h市級城區(qū)地面沉降監(jiān)測的業(yè)務流程，通過對現(xiàn)有的多期地面沉降動態(tài)監(jiān)測數(shù)據(jù)分析和時空處理，構建地面沉降監(jiān)測數(shù)據(jù)庫，并設計了沉降監(jiān)測數(shù)據(jù)管理及更新維護、數(shù)據(jù)后處理、數(shù)據(jù)查詢統(tǒng)計、沉降趨勢分析、沉降模擬預測和可視化表達與分析等功能模塊，開發(fā)建立了比較完善的地面沉降信息管理系統(tǒng)，并加以推廣應用，對沿?？h市級城市輕量級地面沉降信息管理、預警、防震減災、城市規(guī)劃建設等工作都具有積極的意義。