郭 萌，陳 剛，黃 驍,3，王文霞，鄭小燕

（1.北京市地質礦產勘查開發(fā)局，北京 100195; 2.北京市地質工程勘察院，北京 100037; 3.中國地質大學工程技術學院，北京 100083）

北京大興規(guī)劃新城地面沉降研究

郭 萌1，陳 剛2，黃 驍2,3，王文霞2，鄭小燕2

（1.北京市地質礦產勘查開發(fā)局，北京 100195; 2.北京市地質工程勘察院，北京 100037; 3.中國地質大學工程技術學院，北京 100083）

本文分析研究了大興規(guī)劃新城地面沉降發(fā)育特點、成因機理及影響因素，建立模型預測了地面沉降發(fā)展趨勢，為未來新城建設規(guī)劃以及地質災害防治，提供了科學依據。

地面沉降；含水層；壓縮層；沉降模型

北京市的地面沉降主要發(fā)生在平原地區(qū)，目前仍處于快速發(fā)展階段，沉降速率和范圍比較大[1]，已經形成了多個沉降區(qū)[2]。地面沉降對北京城市建設的影響在一些地區(qū)已有顯現(xiàn)，對北京城市建設未來發(fā)展規(guī)劃存在著潛在影響，是地區(qū)經濟可持續(xù)發(fā)展的重要制約因素之一[1]。

北京大興規(guī)劃新城位于榆垡 禮賢沉降區(qū)的西北部，目前的沉降量大于50mm的地區(qū)占到規(guī)劃區(qū)的65%以上，最大累計沉降量達400mm[3]，地面沉降已經成為影響大興新城規(guī)劃與經濟建設的主要地質災害之一。本次研究以規(guī)劃的大興新城（面積約163.03 km2）為中心向四周擴展，研究區(qū)面積達445km2。對研究區(qū)地面沉降發(fā)育特點、成因機理及影響因素進行了分析，預測了地面沉降發(fā)展趨勢，為未來新城規(guī)劃建設以及地質災害防治具有實際意義。

1 地面沉降現(xiàn)狀分析

（1）地面沉降狀況

大興規(guī)劃新城位于榆垡 禮賢沉降中心西北部。本次研究區(qū)南部龐各莊距該沉降中心約2.3km，截止2009年，研究區(qū)累計地面沉降量為40～560mm（見圖1）。其中西北部狼垡、蘆城一帶沉降量較小，為40～50 mm；向南沉降量逐步增大，在北臧村、龐各莊一帶沉降量達400～550 mm。研究區(qū)往南靠近大興禮賢 榆垡沉降中心，地面沉降量較大，北部地面沉降量較小。榆垡 禮賢沉降中心對研究區(qū)南部區(qū)域的影響是明顯的。

（2）地面沉降發(fā)展趨勢

大興榆垡 禮賢沉降中心形成于上世紀70年代后期；1987～1999年，榆垡 禮賢沉降區(qū)面積迅速擴大，以20～35mm/a的沉降速率快速發(fā)展。1999～2005年為地面沉降快速發(fā)展階段，以37mm/a的沉降速率快速發(fā)展[1]。研究區(qū)近幾年地面沉降以每年10～20mm的速率往東北擴展（見圖2），地面沉降速率大致呈逐漸減緩的態(tài)勢（見圖3）。

圖2 研究區(qū)近4年地面沉降量等值線圖

圖3 研究區(qū)部分地段地面沉降速率曲線圖

2 地面沉降成因分析

地面沉降產生的原因是多方面的，一般情況下可分為自然因素和人為因素。自然因素包括地質構造運動和土的次固結沉降等；人為因素包括超量開采地下水、工程建設等。研究區(qū)地面沉降主要與開采地下水和巖土體結構特征有關[2]、[4]、[5]。

（1）過量開采地下水和連續(xù)多年降水偏少導致地面沉降

據水務部門資料，大興區(qū)多年平均地下水可開采量26230.8萬m3；豐水年可開采量32323.1萬m3；平水年可開采量25221.9萬m3；枯水年可開采量20325.5萬m3；特枯水年可開采量15370.7萬m3。按多年平均（29144.2 萬m3）的開采水平，大興區(qū)地下水超采量平均年近3000萬m3，要使地下水達到平衡的降水量為591.2 mm；在豐水年份地下水略有富余，而平水年份超采地下水3500萬m3，枯水年份超采近8500萬m3，特枯水年份超采超過13500萬m3，使得地下水水位不斷下降。

北京地區(qū)枯水年份已經連續(xù)了12年。為了滿足大興區(qū)不斷增加的供水需求，在地下水富水性較好的中北部地區(qū)建立了兩座水廠，每年大量開采地下水，供水能力7.7萬m3/d，供水人口18萬人。從20 世紀80 年代起，隨著降水量的逐年減少，地下水開采規(guī)模不斷擴大，導致地下水位逐年下降，在黃村、龐各莊、榆垡、禮賢等地形成地下水 降落漏斗 。根據1980～2008年地下水年平均埋深值（見圖4），研究區(qū)地下水位總體呈現(xiàn)持續(xù)下降態(tài)勢，地下水埋藏深度從3.07m下降到16.78m。從1984年以來，該區(qū)潛水水位下降了約17m，承壓水水位下降了約10.0～15.0m左右。其中1991～2000年全區(qū)地下水開采量由24900萬m3/a增加到35800萬m3/a，地下水位下降了3m，呈區(qū)域性下降。

圖4 研究區(qū)地下水多年埋深變化過程線

大興區(qū)長期超量開采地下水，造成地下水位的大幅下降導致含水層上覆土層孔隙水壓力降低，使土層固結失水，土層壓縮，形成地面沉降。

地面沉降速率和沉降的大小受地下水位下降速率和下降的幅度控制。由地下水位標高與沉降量關系曲線圖（圖5及圖6）可以看出，地面沉降與地下水位變化密切相關。當?shù)叵滤幌陆邓俾蕼p小時，沉降速率減緩；當?shù)叵滤幌陆邓俾首兇髸r，沉降速率增大。因此，研究區(qū)地面沉降隨地下水開采強度而變化，地面沉降量與地下水位下降量呈正相關關系，超量開采和連續(xù)多年的降雨量偏少導致地下水位下降是引起大興新城工作區(qū)地面沉降的主要原因。

圖5 黃村水位標高與沉降量變化曲線圖

圖6 龐各莊水位標高與沉降量變化曲線圖

（2）地層特性導致南北地面沉降差異

研究區(qū)地面沉降南北差異主要是受其地層巖性空間分布特征影響。工作區(qū)北部土體以單一結構的砂卵石層為主，夾少量粘性土層，含水層顆粒粗，抗壓強度高，可壓縮性低，不易產生壓縮變形，隨著地下水位的升降，其主要表現(xiàn)為彈性變形特征。盡管北部地區(qū)超采地下水引起地下水位大幅下降，但砂卵石層釋水壓縮量很小，且可壓縮的第四系松散層厚度相對較?。?0～70m）（見圖7），因此地面沉降表現(xiàn)不明顯，其1955～2009年地面累計沉降量小于50mm。中部至南部地區(qū)土層逐漸過渡為多層中細砂、薄層砂卵石層與多層粘性土，并以壓縮性較高的粘性土為主的多層土體結構類型。隨著土層粘土含量和厚度的增加，地下水開采引起的土層釋水壓縮量也逐漸增大，且可壓縮的第四系松散層厚度相對較厚（可達250～300m）（見圖7），沉降量逐漸增大。故研究區(qū)地面沉降在空間分布上呈現(xiàn)出南部大于北部的明顯態(tài)勢。

研究區(qū)由北向南整體上土層的壓縮模量逐漸減小，而沉降量逐漸增大（見圖8）。土層壓縮模量越小，表明土層的固結程度低，在荷載作用下更易于壓縮變形。因此當超采地下水引起地下水位大幅下降，土層的孔隙水壓力減小，有效應力增大，南部土層的壓縮變形較大，地面沉降也大。此外研究區(qū)由北向南第四系松散層厚度逐漸增大，可供壓縮的土層厚度增加，沉降量隨之逐漸增大（見圖9）。

總之，地層巖性及結構特征是大興新城區(qū)產生地面沉降的重要地質基礎條件，地下水開采和連續(xù)的枯水年是地面沉降產生的直接誘因。研究區(qū)地面沉降發(fā)生發(fā)展規(guī)律與超量抽取地下水、地層巖性及其結構特征密切相關。

圖7 研究區(qū)第四系厚度等值線圖

圖8 研究區(qū)土層壓縮模量與累積沉降量變化曲線圖

圖9 研究區(qū)第四系松散層厚度與累積沉降量變化曲線圖

3 地面沉降預測

許多學者致力于地面沉降的研究，提出許多計算模型[6-10]，目前對地面沉降的研究主要集中在以地下水開采為主要原因的地面沉降模擬和預測[11]。本次依據地下水動力場 土應力場耦合機理，在充分考慮地層特性、水文地質條件和沉降特點的基礎上，利用Modfow軟件，建立起一個概化為三層結構的非均質各向異性、準三維非穩(wěn)定流地下水系統(tǒng)和地面沉降模型，模擬研究區(qū)2014年地面沉降變化趨勢。

（1）地層概化

由于第四系地層巖性、成因類型、埋藏條件、地下水開采層位等與地面沉降密切相關，而且地下水開采是形成地面沉降的主要原因，因此必須對研究區(qū)的含水層組和壓縮層進行結構劃分[12]。本次將研究區(qū)第四系含水層在垂向上劃分為3層，各層特征見表1。

表1 研究區(qū)概化含水層特征表

經概化后，模型中的第二含水層、第三含水層分別為研究區(qū)的第一、第二壓縮層。各壓縮層特征見表2。

表2 研究區(qū)概化壓縮層特征表

（2）沉降模型擬合與預測

本次研究利用地下水流數(shù)值模擬軟件MODFLOW和IBS軟件包進行模擬求解。根據模型范圍、含水層結構特征、地下水流動特征以及沉降區(qū)特征，將研究區(qū)剖分為100 119個單元間距為200m的矩形網格（見圖10）。

圖10 模型結構及網格剖分示意圖

模型的擬合主要包括水流模型和沉降模型的擬合。水流模型中考慮區(qū)域流場以及單點動態(tài)變化的擬合，沉降模型中主要考慮區(qū)域沉降擬合。單點動態(tài)曲線擬合結果表明，各孔水位計算值與觀測值的變化趨勢一致，淺井的擬合效果略好于深井的擬合效果。區(qū)域沉降使用1955～2009年累計沉降曲線進行擬合，能夠反映沉降趨勢，體現(xiàn)出沉降由南向北逐漸變小的趨勢。

通過該模型計算了未來5年的地面沉降（見圖11），在未來5年內，研究區(qū)地面沉降呈平緩發(fā)展的態(tài)勢，沉降速率大概在8～16mm/a。累計沉降量大的地區(qū)主要位于南部的北高各莊、龐各莊一帶，累積沉降量將大于500mm，成為地面沉降危害較嚴重區(qū)；在規(guī)劃新城范圍內，南部的羅奇營、北臧村、大臧村一帶，累積沉降為300～500mm，為地面沉降危害一般區(qū)；其他區(qū)域處于輕微區(qū)。

圖11 研究區(qū)2014年預測累計沉降量及危險性分區(qū)圖

4 結論

（1）大興規(guī)劃新城區(qū)位于榆垡 禮賢沉降中心西北部。地面沉降將影響新城的規(guī)劃和建設。

（2）研究區(qū)近幾年地面沉降以每年10～20mm的速率往東北擴展，地面沉降速率呈逐漸減緩的態(tài)勢。

（3）地層巖性及結構特征是研究區(qū)產生地面沉降的重要地質背景，地下水開采是地面沉降產生的直接誘因。研究區(qū)地面沉降發(fā)生發(fā)展規(guī)律與超量抽取地下水及地層巖性及其結構特征密切相關。

（4）通過地面沉降模型對地下水動態(tài)以及區(qū)域地面沉降測繪數(shù)據進行擬合；預測得出至2014年，研究區(qū)地面沉降將進一步增大，大興規(guī)劃新城遭受地面沉降的威脅將愈發(fā)嚴重。

[1] 賈三滿等.北京市地面沉降發(fā)展及對城市建設的影響[J].城市地質, 2007.

[2]賈三滿等.北京地面沉降機理研究初探[J].城市地質, 2007.

[3]王翊虹等.北京大興規(guī)劃新城前期區(qū)域工程地質勘查報告.北京市地質工程勘察院,2010.

[4]陳正松.上海地區(qū)地面沉降原因分析[J].大地測量與地球動力學, 2009.

[5]董得茂等.北京地面沉降形成現(xiàn)狀機理與控制對策.全國地面沉降學術研究討會論文集,2002.

[6] 董國鳳．地面沉降預測模型及其應用研究[D],天津大學, 2006.

[7]李成柱,周志芳.地面沉降的數(shù)值計算模型與流固禍合理論綜述[J],勘察科學技術, 2006.

[8] 佘文芳．滄州地區(qū)地面沉降模型研究[D],中國地質大學, 2007.

[9] 周俊等．基于情景分析的天津市濱海新區(qū)地面沉降預測[J],地質災害與環(huán)境保護. 2010.

[10]焉建國等.基于AR模型的上海地區(qū)地面沉降預測分析,大地測量與地球動力學, 2009.

[11]Y Yun and D Kim. Integritymonitoring algorithms using filte2 ring app roaches for higher navigation performance: Consider2 ation of the non - Gaussian GNSS Measurements [ R ]. IONGNSS 20 th ITM of the Satellite Division, Sep tember 2007.

[12] 劉予.北京市平原地面沉降區(qū)含水巖組和可壓縮層劃分[J],城市地質, 2007.

Study on the Surface Subsidence in Beijing Metro Planning Daxing

GUO Meng1, CHEN Gang2, HUANG Xiao2,3, WANG Wenxia2, ZHENG Xiaoyan2
(1.Beijing Geology Prospecting & Developing Bureau, Beijing 100195; 2. Beijing Institute of Geological Engineering, Beijing 100037; 3. China University of Geosciences, Beijing 100083)

This paper studies the growth characteristic, the Mechanism and the infuence factor of the surface subsidence in Beijing Metro Planning Daxing， and has established the subsidence model and forecasted the developing trend of the surface subsidence. This study can provide actual guiding signifcance to the construction and the disaster prevention.

Land subsidence；Aquifer; Compression layer；Subsidence model

X43

A

1007-1903(2011)01-0017-05

郭 萌（1958- ），男，北京市地質礦產勘查開發(fā)局地環(huán)處處長，高級工程師，主要從事地質災害監(jiān)測與評估工作。電子信箱：dkjdhc@163.com