李勝天，聶 飛，王 萃

(江西省地質(zhì)局地理信息工程大隊(duì)，330001，南昌)

0 引言

地面沉降是指地面高程的降低, 又稱地面下沉或地沉, 是以緩慢、難以察覺的向下垂直運(yùn)動(dòng)為主，是指在自然和人為因素作用下，由于地殼表層土體壓縮而導(dǎo)致區(qū)域性地面標(biāo)高降低的一種環(huán)境現(xiàn)象。目前，地面沉降己成為城市化進(jìn)程中普遍存在的生態(tài)環(huán)境問題, 成為世界范圍內(nèi)制約社會(huì)經(jīng)濟(jì)可持續(xù)發(fā)展的重要地質(zhì)災(zāi)害之一。監(jiān)測方面,在傳統(tǒng)D-InSAR技術(shù)基礎(chǔ)上發(fā)展的永久散射體干涉測量技術(shù)(PS)能很好地解決相關(guān)問題,有助于提高形變的時(shí)空分辨率及數(shù)據(jù)處理的精度。Ferretti等首先提出基于幅度變化特征識(shí)別PS點(diǎn)的算法，并把這種技術(shù)稱為永久散射體技術(shù)。隨后在國際上眾多學(xué)者展開了深入研究,如Crosetto等'不斷完善PS干涉測量的技術(shù)方法，并在監(jiān)測地表形變方面取得了很多成功的應(yīng)用成果；而對(duì)地面沉降災(zāi)害風(fēng)險(xiǎn)研究較少，主要集中在地面沉降危險(xiǎn)性分級(jí)標(biāo)準(zhǔn)制定和危險(xiǎn)性評(píng)價(jià)。

1 研究區(qū)概況

杭州市濱江區(qū)，位于浙江省杭州市南，錢塘江下游南岸，為錢塘江泥沙淤積而成的沙土平原(圖1)。地質(zhì)屬錢塘江沖積平原，地勢平坦，地面自然標(biāo)高為5.2～6.2 m(黃海高程)，地表以下5～14 m范圍內(nèi)為粉砂、粉細(xì)砂，地耐力為100～120 kPa。由于城市建設(shè)用地的不斷擴(kuò)張，尤其是城市核心區(qū)的建筑群密度的不斷增大，地面下陷等沉降災(zāi)害頻發(fā)，對(duì)公共基礎(chǔ)設(shè)施和人民群眾的財(cái)產(chǎn)安全造成嚴(yán)重威脅。因此，為了減少沉降災(zāi)害的發(fā)生頻率，在城市建設(shè)中有效規(guī)避沉降災(zāi)害的風(fēng)險(xiǎn)，開展地質(zhì)災(zāi)害風(fēng)險(xiǎn)評(píng)價(jià)尤為迫切。根據(jù)《地質(zhì)災(zāi)害風(fēng)險(xiǎn)調(diào)查評(píng)價(jià)技術(shù)要求(試行)》對(duì)地質(zhì)災(zāi)害風(fēng)險(xiǎn)評(píng)價(jià)的總體要求，需采用定性與定量相結(jié)合的方法，在易發(fā)性、危險(xiǎn)性、易損性評(píng)價(jià)基礎(chǔ)上，劃分為極高、高、中、低4個(gè)等級(jí)，將危險(xiǎn)性和易損性評(píng)價(jià)結(jié)果疊加分析，形成風(fēng)險(xiǎn)評(píng)價(jià)與區(qū)劃結(jié)果。傳統(tǒng)的沉降風(fēng)險(xiǎn)評(píng)價(jià)依賴于實(shí)地測量監(jiān)測技術(shù)，即使是一系列的觀測點(diǎn)數(shù)據(jù)也難以直觀地表現(xiàn)一定區(qū)域內(nèi)地面沉降的空間分布，而有著大范圍、高密度、時(shí)效性優(yōu)勢的基于合成孔徑雷達(dá)差分干涉技術(shù)(InSAR)的地面沉降監(jiān)測手段可以更加方便地提取范圍內(nèi)地表形變速率或形變量，為定量的沉降易發(fā)性評(píng)估提供了依據(jù)；利用層次分析法(AHP)以用地類型、建筑高程、GDP和人口分布4個(gè)權(quán)重因子評(píng)估沉降易損性；在易發(fā)性的基礎(chǔ)上，根據(jù)降雨對(duì)土壤孔隙水壓力的影響確定降雨工況下的沉降失穩(wěn)概率，依靠危險(xiǎn)性指數(shù)法結(jié)合10年一遇，20年一遇，50年一遇，100年一遇的降雨工況分別進(jìn)行沉降災(zāi)害危險(xiǎn)性評(píng)價(jià)。將危險(xiǎn)性和易損性評(píng)價(jià)結(jié)果疊加運(yùn)算，通過矩陣分析法的方式形成風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果，劃分地質(zhì)災(zāi)害風(fēng)險(xiǎn)區(qū)段。

圖1 濱江區(qū)區(qū)位與土地利用圖

2 研究流程與方法

本文分析思路及分析方法主要是通過D-InSAR和PS-InSAR技術(shù)分析sentinel-1A數(shù)據(jù)分別得到地面累積沉降量和年沉降速率信息，運(yùn)用模糊層次分析法建立易發(fā)性指標(biāo),結(jié)合Python代碼設(shè)計(jì)的ROC曲線利用其ACU值進(jìn)行檢驗(yàn)易發(fā)性模型的有效性；通過危險(xiǎn)性指數(shù)法利用不同程度的降水情況獲取失穩(wěn)概率，得到不同降雨工況下的危險(xiǎn)性指標(biāo)；將獲取的土地利用類型數(shù)據(jù)、GDP和人口密度數(shù)據(jù)分權(quán)疊加計(jì)算得到沉降災(zāi)害的易損性指標(biāo)。對(duì)危險(xiǎn)性和易損性指標(biāo)做無量綱化處理，利用矩陣分析法得到沉降風(fēng)險(xiǎn)性結(jié)果。按照《地質(zhì)災(zāi)害風(fēng)險(xiǎn)調(diào)查評(píng)價(jià)技術(shù)要求(試行)》中的分級(jí)要求制作濱江區(qū)沉降風(fēng)險(xiǎn)性評(píng)價(jià)。

3 研究結(jié)果分析

3.1 D-InSAR與PS-InSAR分析結(jié)果

濱江區(qū)研究數(shù)據(jù)基于PS-InSAR的2019年7月至2021年4月21景哨兵1號(hào)數(shù)據(jù)獲取的點(diǎn)數(shù)據(jù)文件，利用空間克里金插值對(duì)于屬性表中的velocity字段進(jìn)行掩膜濱江區(qū)制作的沉降速率柵格圖(單位:mm/a)(圖2)。濱江區(qū)研究數(shù)據(jù)基于D-InSAR的2019年12月至2020年12月哨兵1號(hào)數(shù)據(jù)獲取的累計(jì)沉降量柵格圖(單位：m)(圖3)。

圖2 濱江區(qū)PS_InSAR沉降柵格數(shù)據(jù)

圖3 濱江區(qū)D_InSAR沉降柵格數(shù)據(jù)

3.2 灰色預(yù)測

GM(1,1)模型是常用的灰色預(yù)測模型，其原理是：對(duì)某一數(shù)據(jù)序列用累加的方式生成一組趨勢明顯的新數(shù)據(jù)序列，按照新的數(shù)據(jù)序列的增長趨勢建立模型進(jìn)行預(yù)測，然后再用累減的方法進(jìn)行逆向計(jì)算，恢復(fù)原始數(shù)據(jù)序列，進(jìn)而得到預(yù)測結(jié)果。由此本次使用灰色預(yù)測的GM(1,1)模型對(duì)濱江區(qū)高易發(fā)區(qū)(易發(fā)性綜合得分最高的區(qū)域)當(dāng)中選擇4個(gè)灰色預(yù)測點(diǎn)進(jìn)行未來的升降情況(圖4、表1、表2、圖5、圖6)。

圖4 灰色預(yù)測專題圖

圖5 點(diǎn)號(hào)28851的預(yù)測曲線圖

圖6 點(diǎn)號(hào)119447預(yù)測曲線圖

3.3 易發(fā)性評(píng)價(jià)

濱江區(qū)地形、地質(zhì)構(gòu)造、土質(zhì)情況等自然要素作用，位于杭州灣沖擊平原，河流、錢塘江潮沖擊，堆積大量松散物，地處西湖復(fù)向斜。濱江城市的發(fā)展等人類活動(dòng)；城市建筑建設(shè)施工，高樓林立；采用模糊層次分析法，根據(jù)地面沉降易發(fā)性評(píng)價(jià)指標(biāo)體系，依據(jù)土地利用類型、地面沉降速率、地面沉降累計(jì)量3個(gè)指標(biāo)對(duì)地面易發(fā)性沉降的影響，構(gòu)建模糊一致矩陣；易發(fā)性綜合得分=地面沉降速率*0.533+累計(jì)沉降量*0.333+土地利用類型*0.134(圖7、圖8)。

表1 灰色預(yù)測屬性表

表2 灰色預(yù)測等級(jí)對(duì)照表

圖7 易發(fā)性評(píng)價(jià)模型圖

圖8 易發(fā)性分區(qū)結(jié)果圖

ROC曲線精度驗(yàn)證對(duì)濱江區(qū)易發(fā)性得分利用柵格計(jì)算器進(jìn)行歸一化取值，在其中獲取檢驗(yàn)點(diǎn)來進(jìn)行ROC概率驗(yàn)證(圖9)。ROC曲線縱軸代表真陽性率(True Positive Rate，TPR)，在地質(zhì)災(zāi)害易發(fā)性評(píng)價(jià)中代表災(zāi)害面積頻率；其橫軸代表假陽性率(FalsePositiveRate，F(xiàn)PR)，在地質(zhì)災(zāi)害易發(fā)性評(píng)價(jià)中代表易發(fā)區(qū)面積頻率。通常用ROC曲線下面積AUC作為一個(gè)客觀定量的評(píng)價(jià)指標(biāo)來衡量模型預(yù)測的準(zhǔn)確程度。AUC值介于0～1之間，越接近1模型預(yù)測準(zhǔn)確性越高。AUC值結(jié)果是0.71預(yù)測結(jié)果較好，易發(fā)性模型精度較好。

圖9 ROC曲線結(jié)果圖

3.4 危險(xiǎn)性評(píng)價(jià)

根據(jù)浙江省自然資源廳關(guān)于印發(fā)《浙江省縣(市、區(qū))地質(zhì)災(zāi)害風(fēng)險(xiǎn)普查與鄉(xiāng)鎮(zhèn)(街道)地質(zhì)災(zāi)害風(fēng)險(xiǎn)調(diào)查評(píng)價(jià)工作方案》的通知，確定危險(xiǎn)性評(píng)價(jià)方法為危險(xiǎn)性指數(shù)法，結(jié)合易發(fā)性綜合得分(圖10、圖11)。計(jì)算方法如下：

Hi=Yi/Ymax×Pi

式中：Hi為某種工況下第i個(gè)評(píng)價(jià)單元危險(xiǎn)性指數(shù)(危險(xiǎn)性概率)；Yi為第i個(gè)評(píng)價(jià)斜坡單元易發(fā)性指數(shù)；Ymax為最大易發(fā)性指數(shù)；Pi為某種工況下第i個(gè)評(píng)價(jià)單元的給定時(shí)間段內(nèi)的失穩(wěn)概率；失穩(wěn)概率：Pi=L/Lmax/day(基于極值降雨假設(shè)的Pi確定方法求出4種工況下的失穩(wěn)概率)。

圖10 危險(xiǎn)性評(píng)價(jià)模型

3.5 易損性評(píng)價(jià)

易損性是指受災(zāi)體受到災(zāi)害破壞機(jī)會(huì)的多少與發(fā)生損毀的難易程度。主要由受災(zāi)體自身?xiàng)l件和社會(huì)經(jīng)濟(jì)條件所決定，其目的是分析現(xiàn)有經(jīng)濟(jì)技術(shù)條件下人類社會(huì)對(duì)地質(zhì)災(zāi)害的抗御能力，確定不同社會(huì)經(jīng)濟(jì)要素的易損性參數(shù)。承載體易損性：道路，建筑；社會(huì)經(jīng)濟(jì)易損性：GDP，年末常住人口；采用AHP法對(duì)沉降易損性指標(biāo)賦予權(quán)重，利用python腳本進(jìn)行分析出圖(圖12、圖13)。

圖11 危險(xiǎn)性不同工況分級(jí)圖

圖12 易損性評(píng)價(jià)模型

圖13 易損性分級(jí)圖

3.6 風(fēng)險(xiǎn)性綜合評(píng)價(jià)

根據(jù)浙江省自然資源廳關(guān)于印發(fā)《浙江省縣(市、區(qū))地質(zhì)災(zāi)害風(fēng)險(xiǎn)普查與鄉(xiāng)鎮(zhèn)(街道)地質(zhì)災(zāi)害風(fēng)險(xiǎn)調(diào)查評(píng)價(jià)工作方案》，利用矩陣分析法，結(jié)合易損性和危險(xiǎn)性分析進(jìn)行風(fēng)險(xiǎn)性的評(píng)價(jià)分析，通過二者進(jìn)行歸一化賦值，進(jìn)行風(fēng)險(xiǎn)性分級(jí)劃分(表3)。

3.6.1 風(fēng)險(xiǎn)性指標(biāo)確定 地面沉降地質(zhì)災(zāi)害的風(fēng)險(xiǎn)性主要取決于地面沉降發(fā)生的危險(xiǎn)性及承災(zāi)體的易損性2個(gè)方面。根據(jù)濱江區(qū)易損性分區(qū)以及不同工況下危險(xiǎn)性分區(qū)的評(píng)價(jià)結(jié)果，計(jì)算濱江區(qū)的沉降風(fēng)險(xiǎn)值，得到濱江區(qū)風(fēng)險(xiǎn)性結(jié)果評(píng)價(jià)圖，計(jì)算方法：

Ri=Hi×Vi

式中：Ri為某種工況下第i個(gè)評(píng)價(jià)單元風(fēng)險(xiǎn)性指數(shù)；Hi為某種工況下第i個(gè)評(píng)價(jià)單元危險(xiǎn)性指數(shù)；Vi為第i個(gè)評(píng)價(jià)斜坡單元易損性指數(shù)。

3.6.2 沉降風(fēng)險(xiǎn)性分析 杭州市濱江區(qū)地面沉降的主要因素是建筑施工的土塊挖方導(dǎo)致的地面沉降，以及錢塘江等水域的對(duì)土質(zhì)滲透作用的影響。由于建筑施工會(huì)導(dǎo)致地基沉降進(jìn)而在周圍區(qū)域形成沉降漏斗，因此在水域附近和建筑工地附近區(qū)域都是重點(diǎn)區(qū)域，相關(guān)部門應(yīng)重點(diǎn)關(guān)注，積極采取相關(guān)措施，如控制建筑物密度、在修建道路鐵路等設(shè)施時(shí)應(yīng)避開沉降區(qū)域或者加固道路等，以減少損失。

圖14 暴雨工況下的沉降風(fēng)險(xiǎn)性評(píng)價(jià)圖

4 結(jié)論與展望

本研究在應(yīng)用地質(zhì)災(zāi)害評(píng)價(jià)規(guī)范的基礎(chǔ)上建立了多層次的評(píng)價(jià)模型，整個(gè)過程基本上用模型來實(shí)現(xiàn)數(shù)據(jù)處理與分析，部分評(píng)價(jià)過程使用基于arcpy的python代碼來實(shí)現(xiàn)，可移植性和擴(kuò)展性得到提升，空間適用性廣。該評(píng)價(jià)模型不僅可以用于濱江區(qū)的沉降監(jiān)測，還可以用于其他城市區(qū)域的沉降監(jiān)測。其中易損性評(píng)價(jià)將承災(zāi)體易損性與社會(huì)經(jīng)濟(jì)易損性相結(jié)合，前者通過土地利用類型和建筑高度對(duì)區(qū)域開展建筑易損性分級(jí)，通過道路等級(jí)對(duì)道路緩沖區(qū)進(jìn)行賦值，較好地?cái)M合了城市建筑密集區(qū)域特別是高層建筑密集區(qū)域的易損性特點(diǎn)。時(shí)效性強(qiáng)。哨兵1號(hào)數(shù)據(jù)免費(fèi)，更新快且容易獲取，在此基礎(chǔ)上開展沉降監(jiān)測可省去外業(yè)大范圍監(jiān)測的工作量。使用PS-Insar和D-Insar獲取沉降數(shù)據(jù)可直接得到沉降量空間(時(shí)間)上的變化規(guī)律，并用ROC曲線進(jìn)行精度校驗(yàn)，便于后續(xù)的易發(fā)性評(píng)價(jià)。使用灰色預(yù)測模型對(duì)多個(gè)點(diǎn)位進(jìn)行沉降速率預(yù)測，便于綜合歷史數(shù)據(jù)開展決策。隨著城市化的建設(shè)和開發(fā)，地面沉降漸漸演化成眾多城市的通病，本次項(xiàng)目研究的目的就是為了發(fā)現(xiàn)這一潛在的緩慢下沉的災(zāi)害，通過考察研究，技術(shù)設(shè)計(jì)規(guī)劃，形成一套評(píng)價(jià)體系來應(yīng)對(duì)地面沉降這一問題。