寧波市地面沉降現狀分析及風險評價

高 峰，趙團芝，王文軍，章旬立，鄭凌逶

（1.寧波市自然資源生態修復和海洋管理服務中心，浙江 寧波 315000；2.浙大寧波理工學院土木建筑工程學院，浙江 寧波 315100）

0 引言

城市地面沉降表現為自然與人類活動作用下，地表土層壓縮變形、標高降低，是一種危害較大的緩變型環境地質災害[1-2]，威脅著低海拔城市運行安全與可持續發展。

一直以來，國內外對城市區域性地面沉降的關注度很高，在地面沉降調查與監測、地面沉降特征及分區、地面沉降發展趨勢預測、地面沉降影響因素及控制機理研究、地面沉降緩解與治理措施、地面沉降對社會經濟影響及法規制訂等方面[3-6]開展了廣泛研究。在地面沉降風險評價方面，趙團芝等[7-8]研究表明InSAR 技術在監測區域沉降方面具有可靠性與準確性。工程建設活動及第一軟土層壓縮變形是寧波市工程性地面沉降的主要原因。選取累計沉降量、沉降速率、第一軟土層厚度、建筑容積率作為地面沉降風險性評價指標，繪制寧波市地面沉降風險區劃圖并提出地面沉降防治對策。

風險評價被認為是一種能夠有效管控滑坡、地面沉降、泥石流等災害的理論方法。天津[9-10]、江蘇沿海地區（連云港、鹽城、南通）[11]、南京[12]、滄州[13]、北京[14]和阜陽[15]等地，結合各自地質環境與經濟發展特點，建立反映當地地面沉降特征的風險評價體系并進行區劃分析，目的是減少地面沉降給社會、經濟帶來損失[16]。InSAR 技術能夠在較大的范圍監測地面沉降[17]，但精度低于水準監測點；水準監測點適用于局部埋設，地面沉降結果更加準確[18]，但由于成本原因布設數量有限、覆蓋范圍較小。

本文利用寧波市2017—2020年已有InSAR 監測、水準監測、基巖標等監測數據，從而在較大的區域范圍內相對可靠地獲取地面沉降監測數據，系統分析近年地面沉降分布特征和發展趨勢，建立寧波市地面沉降風險評價體系，編制地面沉降風險區劃圖，從而為加強城市地質安全風險管控提供參考。

1 寧波市地面監測情況及沉降特征

1.1 地面沉降監測情況

寧波市歷來重視地下水監測和地面沉降監測。寧波市地面沉降監測中心始建于1983年，曾是國內第三處較為先進的地質環境監測中心；1999年寧波市啟動了首期“寧波市區地面沉降水準監測”項目；2014年在原地面沉降分層標的基礎上新建地面沉降監測中心，監測范圍覆蓋寧波繞城高速以內的中心城區；同年，寧波啟動了基于InSAR 的地面沉降監測工作。本文采用5 m分辨率Radarsat-2 C InSAR 遙感數據，使用衛星定位校正圖像數據，以提高雷達圖像的準確性和精度；InSAR數據處理主要包括預處理、差分干涉計算、時間/空間域形變估算、形變量計算等幾個步驟，最終獲得平均沉降速率和累積沉降量。

1.2 地面沉降特征

對InSAR 遙感監測數據與全市400 余處沉降監測點數據統一分析，全市2020年平均沉降速率約為5.67 mm/a，年沉降量低于5 mm 的區域面積為1 457.9 km2。寧波市形成的10 余個地面沉降帶，大多位于圍墾區域、城市開發建設熱點區域，地面沉降速率如圖1 所示，其中10 余處沉降帶分布見表1。

表1 地面沉降風險區劃影響因素數據來源Table 1 Data sources of influencing factors for land subsidence risk zoning

圖1 寧波市地面沉降速率概況Fig.1 Overview of land subsidence rates in Ningbo City

對寧波市2020年不同地面沉降速率范圍的區域面積占比進行統計（圖2）。

圖2 寧波市地面沉降速率圖Fig.2 Land subsidence rate map of Ningbo City

可以看到：除了水系、地面沉降不易發區和基巖山區以外，沉降速率低于10 mm/a 的區域面積占比接近88%，全市大部分區域的地面沉降情況穩定、相對安全。沉降速率在10～20 mm/a 的區域面積占比為8.63%，沉降速率在20～30 mm/a 的區域面積占比為2.90%，沉降速率高于30 mm/a 的區域面積占比為0.42%。鄞州中心區和東部新城的代表性監測點沉降量歷年變化情況顯示，隨著城市開發建設的逐步完成，2 個區域的沉降速率逐漸減緩，近3年沉降速率達到了歷年的最低水平，但今年的沉降速率相較于全市平均沉降速率水平仍然偏高。

由圖1 和圖2 總結寧波市地面沉降速率分布特征如下：

（1）濱海區域沉降速率相對較大。鄞州瞻岐、北侖梅山島等多處濱海區域沉降中心的平均沉降速率大于10 mm/a，其中最大值出現在象山新橋鎮圍填海區域，沉降中心平均沉降速率大于30 mm/a。

（2）城市開發建設強度較大的區域沉降速率相對較大。南部商務區、東錢湖旅游度假區等區域沉降中心的平均沉降速率多在20～40 mm/a。

（3）中心城區沉降速率相對較大。海曙高橋、江北洪塘等主要沉降中心的沉降速率多在10～30 mm/a。

綜上所述，工程建設引發的地面沉降在城市規劃建設集聚區日益凸顯，比如東部新城、鄞州中心區、奉化主城區和沿海產業帶等開發強度大、地質環境條件脆弱的地區年度沉降速率較大；大規模工程建設活動是寧波市地面沉降發生發展的主要因素，在后續地面沉降風險評價體系構建中需加以考慮。

2 地面沉降風險評價體系構建

2.1 評價體系說明

地面沉降風險評價屬于對復雜系統的定量化分析過程。一般采用線性加權模型，而權系數的確定又可通過層次分析、專家打分或二者綜合得出[19]。地面沉降的風險評價模型通?？杀硎緸椋篒=I危險性·I易損性，分析災害因子的危險性、承災體的脆弱性，根據地面沉降危險性評價指標、易損性評價指標結合風險評價模型開展寧波市地面沉降風險評價研究。前人研究[3,7,19]表明地面沉降危險性評價指標中具有代表性的因子主要包括：地面累計沉降量、沉降速率、軟土層厚度、填土厚度、地面建筑荷載、地下水位變化速率、地殼運動速率等；地面沉降易損性評價指標中具有代表性的因子主要包括：地面高程、沉降易發程度、交通線路分布、人口密度、單位面積GDP、建成區面積比重、重要水系密度等。

除了地面高程、地面累計沉降量、沉降速率等常規地面沉降風險因子，對于人口凈流入為正、城市建設高速發展的寧波市，城市人口密度、建設用地比重也應被重點考慮。從寧波市不同區域的地面沉降歷史數據來看，地面沉降速率較高的區域如市區、濱海圍墾區域、城市開發熱點區域等，與上述因子密切相關。

寧波市位于中國東南沿海的潮汐平原區，具有深厚軟土層和廣泛分布的松散層，含水層較多；地下水開采強度對地面沉降的影響已深入人心，且多年平均沉降速率預示著沉降易發程度，因此對于寧波市而言，除了地面沉降危險性評價因子、易損性評價因子，地面沉降易發程度作為1 個綜合性評價因子十分必要。

因此，本文構建的地面沉降風險評價體系如圖3 所示。首先考慮地質條件、水文地質條件、人為活動等的影響，對地面沉降易發程度進行分區，進而將地面沉降易發程度作為地面沉降風險評價體系中的1 個綜合性評價因子，并聯合地面高程、地面累計沉降量、沉降速率、城市人口密度、建設用地比重等5 個評價因子進行評價。

圖3 寧波市地面沉降風險評價體系Fig.3 Land subsidence risk assessment system for Ningbo City

2.2 評價體系權重及分級

浙江省自然資源廳印發的《浙江省城市地質安全風險“一張圖”編制技術指南》根據浙江省沿海平原地面沉降形成機理，結合長期地面沉降監測成果，通過“經驗總結+專家打分”相結合的方法，采用如德爾菲法（delphi method）或層次分析法（analytic hierarchy process，AHP）來匯總專家的意見、通過多輪迭代和反饋，使專家意見逐漸收斂，綜合確定了各評價指標權重。各項指標權重及等級劃定見表2。

表2 地面沉降風險區劃影響因素權重及分級Table 2 Weight and classification of factors influencing land subsidence risk zoning

通過逐一確定各評價因素的等級，再將各因素等級分值與其權重相乘得到單因素分值，最后將各因素分值累加，從而得到地面沉降綜合風險指數（W），按表3 劃分風險區等級。軌道交通建設過程中地下盾構隧道開挖也容易引起地面沉降發生發展[20-22]，根據浙江省自然資源廳印發的《浙江省城市地質安全風險“一張圖”編制技術指南》，將高鐵、軌道交通等重大工程保護區兩側200 m 范圍劃為地面沉降風險防控帶。

表3 地面沉降風險區等級劃分表Table 3 Classification table of land subsidence risk zone levels

綜合風險指數（W）計算公式如下：

式中：W——地面沉降綜合風險指數；

aj——j類影響因素的權重；

bj——j類影響因素的分值。

2.3 地面沉降易發程度分區

以區內第四紀松散層分布厚度、土層力學性質及空間分布、含水層發育層數及富水性、淺部主要壓縮層發育等為主要因素，開展寧波市地面沉降易發程度評價；根據評價結果進行易發程度分區，分為高易發、中易發、低易發和不易發區等四個級別。根據《地質災害危險性評估規范》（DB33/T 881—2012）[23]，用于地面沉降易發程度等級劃分的綜合分值計算公式為：

式中：A——綜合分值；

σi、τi——分別為地面沉降i項影響因素相應的量化指標賦分及權重，參考《地質災害危險性評估規范》（DB33/T 881—2012）[23]取值。

地面沉降易發程度分區結果表明：寧波市地面沉降無高易發區，中、低易發區主要與區內全新世軟土層厚大、歷史上大量開采地下水、局部高強度城市建設以及沿海圍墾工程等因素有關。

3 風險區劃結果與管控建議

3.1 風險區劃結果

風險區劃采用的影響因素包括地面高程、易發程度、地面累積沉降量、沉降速率、城市人口密度和建設用地比重，各項數據的來源情況見表4。

表4 地面沉降風險區劃影響因素數據來源Table 4 Data sources of influencing factors for land subsidence risk zoning

地面沉降風險等級計算和分析采用柵格計算方式進行。首先將表4 中6 類影響因素的源數據導入至GIS軟件并進行整理；然后分別對各影響因素的數據結果進行柵格劃分，即由矢量形式轉為柵格形式（柵格像元大小40 m×40 m），以城市人口密度為例，可得如圖4（a）所示結果；進一步根據表2 中的分值標準對各影響因素的柵格數據進行重新賦值得到如圖4（b）所示結果；接著利用GIS 軟件中的柵格計算工具將各影響因素的柵格分值數據乘以相應的權重并相加，獲得結果后根據表3進行地面沉降風險區等級劃分；最后對劃分結果進行矢量化即可得到最終的地面沉降風險區劃圖。本文中地理底圖采用國家最新版 CGCS 2000 大地坐標系，成圖的比例尺為 1∶50 000，圖中僅保留必要水系、交通、行政界線、地名等主要地理要素，以不同顏色代表不同風險區劃結果。

圖4 寧波市城市人口密度柵格數據圖Fig.4 Grid data map of urban population density in Ningbo City

寧波市地面沉降風險區劃結果如圖5 所示，寧波市全域無高風險區；中風險區面積132.6 km2，占比3.2%；低風險區面積2 214.8 km2，占比53.0%；中風險區、低風險區、水系、基巖山區以外的區域，為不易發區，面積1 827.9 km2，占比43.8%。相應地，本文1.2 節也對寧波市2020年不同地面沉降速率區域的面積占比進行了統計，沉降速率>20 mm/a 面積占比3.3%、5～20 mm/a 面積占比34.5%、0～5 mm/a 面積占比62.1%。地面沉降速率等監測結果作為地面沉降風險評價的關鍵指標之一，與地面沉降風險區劃結果存在著較強的相關性。需要注意的是，地面沉降風險也受到軟土層厚度、含水層數量、人口密度、建設用地比重等因素的影響，因此風險區劃分結果與地面沉降監測結果又存在一定的差別。

3.2 管控建議

地面沉降會對建（構）筑物和各類基礎設施帶來影響，造成建筑物地基下沉、房屋開裂、道路沉降、地下管線受損，形成的地面沉降漏斗會降低防洪排澇工程效能，臺風暴雨時會造成大面積積水、加劇洪澇災害。大面積均勻/不均勻沉降也會影響城市建設和資源開發，不利于城市基礎設施的運營維護；對于濱海城市，造成的高程損失會帶來海水倒灌等風險。因此，針對寧波市不同的地面沉降風險區分別提出管控建議（表5）。

表5 寧波市地面沉降風險區管控建議一覽表Table 5 Summary of control recommendations for ground subsidence risk zones in Ningbo City

4 結論

（1）結合2017—2020年寧波市InSAR 遙感監測數據與水準監測點數據，利用InSAR 技術在較大范圍內監測地面沉降、利用水準監測點更為準確地在局部測量地面沉降，從而在較大區域范圍內相對可靠地獲取地面沉降數據。結果表明全市形成的十余個沉降帶，大多位于圍墾區域、城市開發建設熱點區域。

（2）選取地面高程、沉降易發程度、地面累計沉降量、沉降速率、城市人口密度、建設用地比重等6 個評價因子建立地面沉降風險評價體系；其中沉降易發程度為綜合考慮地質條件、水文地質條件、人為活動影響后的綜合性評價因子。從而根據寧波市沿海平原地面沉降形成機理，結合長期地面沉降監測成果，采取“經驗總結+專家打分”相結合的方法綜合確定各評價指標、權重、分級，形成評價體系。

（3）地面沉降風險區劃結果表明，寧波市全域無高風險區，中風險區面積占比3.2%，低風險區面積占比53.0%，不易發區面積占比43.8%；設置高鐵、軌道交通等重大工程保護區兩側200 m 范圍為地面沉降風險防控帶。在此基礎上，提出了相應的管控建議如嚴格執行地下水開采管理、優化地面沉降監測網絡、針對重大工程開展監測及管理等。