劉 斌，唐雅玲，馬晨燕，白少云

(1.云南省水利水電勘測(cè)設(shè)計(jì)研究院，云南 昆明 650021；2.武漢大學(xué)資源與環(huán)境科學(xué)學(xué)院，湖北 武漢 430079)

在日益復(fù)雜的全球氣候變化、海平面上升與不斷深入的城鎮(zhèn)化過(guò)程的耦合作用下，城市的防洪安全正面臨著嚴(yán)重的威脅。在分析國(guó)內(nèi)外城市洪澇災(zāi)害研究過(guò)程中發(fā)現(xiàn)，基于情景模擬分析的災(zāi)情模擬、動(dòng)態(tài)顯示及預(yù)警研究將成為今后災(zāi)害風(fēng)險(xiǎn)評(píng)估研究的必然趨勢(shì)[1]?？偟膩?lái)說(shuō)，基于情景模擬分析的洪澇風(fēng)險(xiǎn)評(píng)估方法至今已形成了較為完善的模型體系[2-6]。但該方法對(duì)研究區(qū)地理數(shù)據(jù)的完整性及精度要求較高，包括地形圖、河網(wǎng)數(shù)據(jù)、地下管網(wǎng)數(shù)據(jù)、高精度DEM、地表覆蓋數(shù)據(jù)等，若要進(jìn)一步提高風(fēng)險(xiǎn)評(píng)估的準(zhǔn)確性，從提高城市雨洪模型的精度出發(fā)是一個(gè)很好的思路。

目前，數(shù)字城市三維建模的數(shù)據(jù)主要來(lái)源于航空影像、衛(wèi)星影像和近景影像3種類(lèi)型。建模的方法主要為測(cè)量建模，而無(wú)論是基于傳統(tǒng)的航空攝影測(cè)量還是機(jī)載LiDAR的三維建模都需要耗費(fèi)大量的生產(chǎn)成本，且構(gòu)建的三維模型存在建筑物側(cè)面紋理信息不足、遮擋嚴(yán)重或覆蓋范圍有限的問(wèn)題，均難以生成滿(mǎn)足精度需求的城市三維模型[7-10]。利用傾斜攝影方式，獲取更高分辨率更豐富細(xì)節(jié)表現(xiàn)的多視影像，制作實(shí)景三維模型，既可以為城市洪澇模型提供高精度地理信息支持，又能增強(qiáng)情景模擬分析的動(dòng)態(tài)演示[11-15]。本文提出基于無(wú)人機(jī)傾斜攝影測(cè)量技術(shù)，結(jié)合GIS空間分析能力及情景模擬分析方法，從提高城市實(shí)景模型精度方面對(duì)城市雨洪進(jìn)行模擬分析。

1 傾斜攝影測(cè)量技術(shù)輔助雨澇災(zāi)情模擬

傾斜攝影測(cè)量技術(shù)是充分吸收了近幾年來(lái)視覺(jué)測(cè)量領(lǐng)域的最新成果而逐步發(fā)展起來(lái)的一項(xiàng)新的數(shù)據(jù)采集技術(shù)，傾斜方式可獲取到豐富的地物頂部及側(cè)面的高分辨率紋理。它不僅能夠真實(shí)地反映地物現(xiàn)狀，高精度地獲取物方紋理信息，還可通過(guò)先進(jìn)的定位、融合、建模等技術(shù)，生成真實(shí)的地表三維模型[9]。與傳統(tǒng)攝影測(cè)量產(chǎn)生的數(shù)字地表模型不同，傾斜攝影測(cè)量獲得的真三維模型采用空間三角網(wǎng)的方式描述，能描述地表復(fù)雜的地物，如房屋、街道、橋梁、陡坡、溝壑、公園、施工現(xiàn)場(chǎng)等。無(wú)人機(jī)搭載輕型傾斜攝影測(cè)量系統(tǒng)可在城市低空200～500 m的高度進(jìn)行快速的多視角影像采集，一架無(wú)人機(jī)可在幾天之內(nèi)完成幾十平方千米的城市多視角影像采集。使用地面控制點(diǎn)，通過(guò)軟件大規(guī)模并行計(jì)算，實(shí)現(xiàn)基于海量影像的區(qū)域網(wǎng)平差解算，獲取影像外方位元素，進(jìn)而以多機(jī)GPU并行加速為核心，實(shí)現(xiàn)城市三維模型的重建。重建的實(shí)景模型結(jié)構(gòu)完整，紋理清晰，在準(zhǔn)確性和效率上都優(yōu)于人工建模。

基于無(wú)人機(jī)傾斜攝影測(cè)量的方式已經(jīng)逐漸成為數(shù)字城市三維建模的主要技術(shù)手段，帶有地理信息的實(shí)景模型為城市洪澇分析提供重要的基礎(chǔ)信息，很好地彌補(bǔ)了傳統(tǒng)建模手段受限于數(shù)據(jù)分辨率過(guò)低而導(dǎo)致三維模型的精度較低的不足?；谥亟ㄈS模型，結(jié)合GIS空間分析及情景模擬，能夠增強(qiáng)災(zāi)情模擬動(dòng)態(tài)展示，為雨澇災(zāi)害防治提供服務(wù)。無(wú)人機(jī)傾斜攝影輔助雨澇災(zāi)害分析技術(shù)路線(xiàn)如圖1所示。

2 基于Smart3D的城市三維重建

(1)工程準(zhǔn)備。將外業(yè)采集回來(lái)的原始數(shù)據(jù)根據(jù)不同視角的相機(jī)分別存儲(chǔ)，完成Smart3D工作集群建立后新建Smart3D工程，將準(zhǔn)備好的影像文件、相機(jī)文件和POS文件導(dǎo)入Smart3D。

(2)空中三角測(cè)量。準(zhǔn)確的相機(jī)外方位元素是三維重建的必要條件，通過(guò)自動(dòng)連接點(diǎn)提取和地面控制點(diǎn)控制，完成空中三角測(cè)量。

(3)三維重建。在本文的試驗(yàn)環(huán)境下，慮及計(jì)算器集群的硬件設(shè)置及處理效率，進(jìn)行了分塊處理。對(duì)若干個(gè)瓦片進(jìn)行單獨(dú)的重建，按照密集點(diǎn)云生成、TIN模型構(gòu)建和紋理自動(dòng)映射3個(gè)步驟完成帶有LOD金字塔結(jié)構(gòu)的三維模型的構(gòu)建，得到如圖2所示的三維實(shí)景模型。

通過(guò)建立的三維實(shí)景模型可實(shí)現(xiàn)研究區(qū)的量測(cè)，包括任意點(diǎn)的三維坐標(biāo)信息、建筑高度面積信息等。然而，僅僅進(jìn)行簡(jiǎn)單的量測(cè)并不能實(shí)現(xiàn)城市雨洪的風(fēng)險(xiǎn)分析，借助于Smart3D強(qiáng)大的數(shù)據(jù)生產(chǎn)功能，生成0.1 m精度的DSM和DOM數(shù)據(jù)產(chǎn)品。通過(guò)無(wú)人機(jī)傾斜實(shí)景建模手段生產(chǎn)的數(shù)據(jù)產(chǎn)品在影像色彩、紋理精細(xì)程度上都優(yōu)于常規(guī)方式生產(chǎn)的產(chǎn)品效果，更有利于進(jìn)一步的研究分析。

3 基于情景模擬的城市雨洪風(fēng)險(xiǎn)評(píng)估

3.1 研究區(qū)處理

由于原始數(shù)據(jù)產(chǎn)品為分塊影像，處理前先進(jìn)行影像鑲嵌，大致勾畫(huà)出研究區(qū)邊界，得到矢量面數(shù)據(jù)，然后利用研究區(qū)邊界裁剪DOM、DSM，剔除不相關(guān)的數(shù)據(jù)。

結(jié)合高精度正射影像及Acute3D Viewer，對(duì)研究區(qū)域內(nèi)的房屋屋頂輪廓進(jìn)行矢量化處理。盡管結(jié)合Acute3D Viewer的三維觀察，可以盡可能地提高輪廓的精度，但是由于是人工手動(dòng)矢量化處理，加上在三維視角下不可避免的房屋遮擋的因素，仍會(huì)存在一定的誤差。圖3分別為研究區(qū)域DOM、DSM及矢量化圖。

3.2 城市降雨情景模擬

3.2.1 城市暴雨模型

降雨是內(nèi)澇災(zāi)害發(fā)生的主要致災(zāi)源，降雨量是風(fēng)險(xiǎn)評(píng)估中的主要驅(qū)動(dòng)變量。城市的降雨表現(xiàn)出在時(shí)間和空間分布上的不均勻性。本文城市暴雨模型利用研究區(qū)本地的暴雨強(qiáng)度公式，計(jì)算出降雨過(guò)程中任意時(shí)間段內(nèi)的平均降雨強(qiáng)度，模擬暴雨在時(shí)間上的變化特征。因?yàn)楸疚倪x取的研究區(qū)域較小，且地理環(huán)境特征較為一致，在本文中假設(shè)研究區(qū)是一個(gè)封閉的區(qū)域，不考慮與相鄰區(qū)域的連通性，并假設(shè)降雨在空間上均勻分布，根據(jù)研究區(qū)現(xiàn)行的暴雨強(qiáng)度公式，按不同歷時(shí)，計(jì)算不同重現(xiàn)期內(nèi)的暴雨強(qiáng)度，該暴雨強(qiáng)度乘以降雨歷時(shí)即可得到整個(gè)研究區(qū)內(nèi)的降雨總量，然后將降雨量均勻分布于整個(gè)研究區(qū)范圍內(nèi)，根據(jù)暴雨強(qiáng)度公式求出的暴雨強(qiáng)度單位為L(zhǎng)/(S·ha)，為了統(tǒng)一量綱(暴雨強(qiáng)度單位為mm/h)，對(duì)原公式進(jìn)行修正后，其表達(dá)式為

(1)

式中，T表示暴雨重現(xiàn)期,單位為a；q表示暴雨強(qiáng)度；t表示降雨歷時(shí)。

在城市降雨過(guò)程中，會(huì)產(chǎn)生各種不同程度的徑流損失，如由于城市地表的不平整產(chǎn)生的植被截留、蒸散發(fā)、下滲、地下徑流等，實(shí)際應(yīng)用中不可忽略這些因素對(duì)徑流量的影響。而在本文中，由于數(shù)據(jù)的缺失，則不考慮這些徑流損失的影響，將研究區(qū)內(nèi)的降雨總量作為徑流量。在城市產(chǎn)匯流過(guò)程中，還需要考慮城市的排水能力，在本文中，研究區(qū)域的排水標(biāo)準(zhǔn)按照一年一遇(48 mm/h)的降雨強(qiáng)度設(shè)計(jì)，在排水管網(wǎng)數(shù)據(jù)無(wú)法獲取的情況下，假設(shè)研究區(qū)內(nèi)排水管網(wǎng)的空間分布與排水能力一致。因此，最終研究區(qū)內(nèi)的內(nèi)澇積水量通過(guò)徑流量與排水量的差，乘以研究區(qū)面積計(jì)算得出。其公式表示為

W=(Q-V)×S

(2)

式中，W表示內(nèi)澇積水量；Q表示徑流量；V表示排水量(48 mm/h)；S表示研究區(qū)面積。

3.2.2 城市暴雨內(nèi)澇情景模擬

城市暴雨降水往往集中在前1～2 h內(nèi)，且氣象和市政工程暴雨計(jì)算標(biāo)準(zhǔn)也多參考1 h降雨量，因此，本文模擬了研究區(qū)內(nèi)1 h降雨歷時(shí)下重現(xiàn)期為5、10、20、50、200、500和1000 a的城市暴雨內(nèi)澇情景。根據(jù)計(jì)算得出研究區(qū)的總積水量，假設(shè)積水處于相對(duì)靜止?fàn)顟B(tài)，采用簡(jiǎn)化的處理方式，基于研究區(qū)高精度DSM，模擬內(nèi)澇淹沒(méi)區(qū)域和淹沒(méi)深度。將暴雨內(nèi)澇看作一種無(wú)源淹沒(méi)狀態(tài)，凡是高程值低于水位值的區(qū)域均計(jì)入淹沒(méi)區(qū)，借助Python的Arcpy模塊實(shí)現(xiàn)一定積水量下研究區(qū)域的淹沒(méi)面積和淹沒(méi)水面的高程值的計(jì)算，注意該淹沒(méi)水位高程為水準(zhǔn)高程。其編程思想為：假設(shè)洪水的淹沒(méi)水面永遠(yuǎn)是一個(gè)水平面，首先設(shè)定一個(gè)水位值，用該水位值作為高程屬性值的一個(gè)水平面，去切研究區(qū)DSM，得到該水位值對(duì)應(yīng)的水平面與DSM之間的庫(kù)容量；然后比較庫(kù)容量與積水量，若積水量小于庫(kù)容量，則減小設(shè)定的水位值，若積水量大于庫(kù)容量，則增加水位值，通過(guò)不斷調(diào)整水位值，讓庫(kù)容量漸漸逼近積水量，若庫(kù)容量等于積水量，則對(duì)應(yīng)的水位值就是洪水最終淹沒(méi)水面的高程值；最后根據(jù)上述方法，計(jì)算出研究區(qū)內(nèi)不同降雨情景下的淹沒(méi)水位、淹沒(méi)水位以及淹沒(méi)區(qū)面積，由于在重現(xiàn)期為5、10時(shí)，研究區(qū)內(nèi)除了低洼的水塘外幾乎沒(méi)有淹沒(méi)情況，則最終顯示重現(xiàn)期分別為20、50、200、500和1000 a的淹沒(méi)情況。

在實(shí)際應(yīng)用時(shí)，要分布式地考慮地物的蓄滯屬性和透水性，特別是在本次試驗(yàn)中，在影像圖上，可以看到研究區(qū)內(nèi)分布著較多植被，在每個(gè)柵格點(diǎn)的高程值里，每個(gè)植被都有其高程屬性，對(duì)研究區(qū)進(jìn)行淹沒(méi)分析時(shí)，會(huì)發(fā)現(xiàn)在一定高程下，植被會(huì)被歸為未被淹沒(méi)的區(qū)域，因此，在作實(shí)際淹沒(méi)分析時(shí)需要剔除掉植被對(duì)淹沒(méi)分布的影響。在本次試驗(yàn)中，結(jié)合正射影像將研究區(qū)內(nèi)的植被勾畫(huà)出來(lái)，并利用三維模型對(duì)植被高程進(jìn)行量測(cè)，給植被賦予高度值，以近似地剔除掉植被對(duì)淹沒(méi)結(jié)果的影響，修正公式為

Vij=Dij-Hij(i=1,2，…,n;j=1,2，…,m)

(3)

式中，i表示行數(shù)；j表示列數(shù)；Vij表示修改后的柵格單元的值；Dij代表原地形柵格單元值；Hij表示植被高度柵格單元值。根據(jù)經(jīng)植被修正后的研究區(qū)進(jìn)行計(jì)算分析，分別統(tǒng)計(jì)不同降雨情景下的研究區(qū)淹沒(méi)情況(見(jiàn)表1)。

表1 不同重現(xiàn)期下的淹沒(méi)情況

3.3 城市雨洪風(fēng)險(xiǎn)等級(jí)評(píng)估

根據(jù)3.2節(jié)通過(guò)模擬不同重現(xiàn)期下研究區(qū)的淹沒(méi)情況，用不同的淹沒(méi)水位裁切淹沒(méi)區(qū)域則可以得到淹沒(méi)區(qū)不同的淹沒(méi)深度，以淹沒(méi)深度作為風(fēng)險(xiǎn)評(píng)估的關(guān)鍵指標(biāo)，借鑒其他文獻(xiàn)的風(fēng)險(xiǎn)評(píng)價(jià)方法，將本文研究區(qū)暴雨洪澇災(zāi)害風(fēng)險(xiǎn)程度劃分為6個(gè)等級(jí)，見(jiàn)表2。

表2 城市雨洪風(fēng)險(xiǎn)等級(jí)

不同重現(xiàn)期下研究區(qū)的淹沒(méi)深度分布如圖4所示，并利用高精度的影像和高程信息，使用Terra Builder制作MPT文件，然后使用Terra Explorer Pro進(jìn)行淹沒(méi)場(chǎng)景的三維展示，可以很直觀地觀測(cè)到研究區(qū)在不同降雨情景下各個(gè)區(qū)域的淹沒(méi)情況，包括淹沒(méi)范圍、淹沒(méi)深度情況等。

通過(guò)模擬研究區(qū)域不同重現(xiàn)期下的暴雨情況，可以得到不同重現(xiàn)期下研究區(qū)的淹沒(méi)情況，通過(guò)對(duì)淹沒(méi)水深的提取及不同淹沒(méi)水深下研究區(qū)淹沒(méi)情況的三維顯示，可以很直觀地看到研究區(qū)內(nèi)每個(gè)地方的淹沒(méi)情況，為城市建筑的規(guī)劃布局提供一定的參考價(jià)值。

4 結(jié) 語(yǔ)

傾斜攝影技術(shù)基于多視影像，有效解決了傳統(tǒng)航測(cè)建模方式建筑物側(cè)面信息缺失的問(wèn)題，且其產(chǎn)品更加豐富，在生產(chǎn)實(shí)景三維模型的同時(shí)也能生產(chǎn)DOM、DSM等多種產(chǎn)品，大大增強(qiáng)了產(chǎn)品使用范圍，滿(mǎn)足不同工程需求。本文以淹沒(méi)深度為關(guān)鍵指標(biāo)進(jìn)行研究區(qū)洪澇風(fēng)險(xiǎn)的評(píng)估，基于Terra Explorer Pro的不同降雨情景下的淹沒(méi)三維模擬，直觀表達(dá)了淹沒(méi)區(qū)的分布隨淹沒(méi)水位的抬升的變化情況，為城市防洪減災(zāi)決策、制定城市發(fā)展規(guī)劃和土地利用規(guī)劃，提高城市對(duì)洪澇災(zāi)害的適應(yīng)能力提供了可借鑒的方法。