賈超++張國榮++王嘉斌++張永偉

摘要：地下水開采誘發(fā)的不良環(huán)境及地質問題日益凸顯，其中地面沉降致災問題已廣受關注。華北平原山東區(qū)是若干地面沉降顯著區(qū)之一，地下水開采是其主要誘發(fā)因素。針對地面沉降問題，以華北平原德州區(qū)為典型工程背景，采用水文地質工程地質與巖土力學相結合的方法，通過現(xiàn)場及室內試驗，在獲得研究區(qū)巖土力學參數(shù)的基礎上，建立流固耦合數(shù)值模型，對研究區(qū)地面沉降的歷史進行重演，對地面沉降的發(fā)展趨勢進行預測。在獲得研究區(qū)流場及地面沉降基本規(guī)律的基礎上，提出減少地下水開采量、定期回灌地下水、調整開采井開采層位及開采井平面布局等工程措施進行地面沉降控制，并采用建立的數(shù)值模型從定量分析角度評價這些措施的效果。

關鍵詞：地下水；地面沉降；流固耦合；水文地質；工程地質；巖土力學；華北平原

中圖分類號：P641；P642.26文獻標志碼：A

0引言

地面沉降是一種自然因素或人類工程活動引發(fā)的地下松散地層固結壓縮，從而導致一定區(qū)域范圍內地面高程降低的緩變地質災害現(xiàn)象。在眾多誘發(fā)因素中，過量抽取地下水是造成地面沉降最主要的原因之一[1]。根據(jù)相關資料統(tǒng)計報道，中國目前已有超過50個城市發(fā)生了明顯的地面沉降現(xiàn)象，其中最典型區(qū)域有長江三角洲、華北平原和汾渭盆地[2]。針對這些地區(qū)的地面沉降災害問題，已有薛禹群、吳吉春、駱祖江、彭建兵、殷躍平、邵景力等眾多專家學者進行了長期持續(xù)的工作，并取得了顯著的研究成果[118]。相比于其他地區(qū)，華北平原山東地區(qū)地面沉降的機理研究工作尚顯不足。在華北平原，山東地面沉降尤以魯北地區(qū)最為嚴重，德州是魯北地區(qū)地面沉降的典型代表。德州地區(qū)地面沉降已對當?shù)丨h(huán)境及生產(chǎn)造成嚴重影響，制約當?shù)厣鐣?jīng)濟的可持續(xù)發(fā)展[3]，并且根據(jù)該區(qū)地面沉降長期監(jiān)測資料顯示，地面沉降尚無明顯停止的跡象[48]。

已有研究成果大多針對德州地區(qū)地下水開采與地面沉降的關系，應用統(tǒng)計學方法進行數(shù)理統(tǒng)計回歸分析，但僅有這些成果并不能準確認識到該區(qū)地面沉降的深層機理及沉降本質。地下水開采誘發(fā)地面沉降的本質原因是：地下水開采導致各地層應力調整，引發(fā)各地層壓縮總和在地面上顯現(xiàn)。因此，采用水文地質工程地質與巖土力學相結合的方法更能揭示德州地區(qū)地面沉降的本質規(guī)律。該方法已被廣泛應用于其他地區(qū)地面沉降研究中，但在德州地區(qū)尚無這方面的研究成果。

基于此，本文通過現(xiàn)場工作獲取進尺達800 m的全取芯鉆孔試樣，進行試樣的實驗室測試，取得翔實的物理力學參數(shù)；分析流場分布規(guī)律，建立地面沉降流固耦合分析模型，重演歷史沉降情景，預測沉降發(fā)展趨勢；針對沉降發(fā)展趨勢及規(guī)律，提出若干控制地面沉降的工程措施，并對提出的措施進行綜合對比分析[1920] ，為山東德州地區(qū)防治地面沉降理論及分析計算奠定基礎。

1研究區(qū)概況

山東省德州市德城區(qū)沉降資料及前期工作相對較好。該區(qū)東與德州市陵城區(qū)接壤，南與山東省平原縣、武城縣毗鄰，西、西北、東北分別與河北省故城縣、景縣、吳橋縣連接。其經(jīng)度范圍為116°10′15″E～116°30′40″E，緯度范圍為 37°10′15″N～37°30′30″N，面積為539 km2。

研究區(qū)地下水多年處于超量開采狀態(tài)，地下水位呈持續(xù)下降趨勢，并伴隨嚴重的地面沉降問題。該區(qū)自1989年開始進行地面精密測量，2006年400 mm沉降等值線圈定沉降漏斗面積為344 km2，歷年沉降發(fā)展過程見圖1。

圖1山東省德州市德城區(qū)地面沉降發(fā)展過程

Fig.1Land Subsidence Development Process in Decheng District of Dezhou City， Shandong Province

2現(xiàn)場及室內試驗

山東省德州市德城區(qū)在工程地質深層鉆孔資料方面比較欠缺。山東省國土資源部門在該區(qū)進行了2處800 m深層全取芯鉆孔鉆探工作，山東大學承擔了京臺高速公路德州服務區(qū)鉆位所得的巖芯參數(shù)測試工作。該處鉆孔共取173個原狀土樣、15個平行樣和37個擾動砂樣（圖2）。

根據(jù)現(xiàn)場地質鉆探和室內試驗結果，對德城區(qū)地層以工程地質及水文地質相近的原則進行地層概化。概化后德城區(qū)工程地質特征層分布見圖3；實驗室測試所得各地層物理力學參數(shù)取值范圍見表1。

3數(shù)值模型

3.1數(shù)值模型的建立

以圖3劃分的工程地質特征層為所考慮的地層分布，以德州市德城區(qū)行政邊界為模型外邊界，建立研究區(qū)數(shù)值模型，總面積為539 km2。

數(shù)值模型共納入開采深井263口，按各井實測流量賦值，部分井點信息見表2。井位分布見圖4。

德州地面沉降觀測井2號孔396～400 m深度處

圖2現(xiàn)場巖芯鉆孔地質編錄

Fig.2Geological Record of Field Core Hole

開采井直徑相比于整個區(qū)域面積非常微小，不利于數(shù)值模型的建立，因此，借鑒Oh等的研究成果[15]，將開采井概化為線，開采井邊界條件采用質量流量進行表示。該方法既不會影響計算結果的精確度，又可使大量開采井實際位置及實際開采量

圖3工程地質特征層劃分

Fig.3Division of Engineering Geological Feature Layers

加入大區(qū)域模型得以順利實現(xiàn)；經(jīng)過驗證，該方法正確有效。

表1地層物理力學參數(shù)

Tab.1Physical Mechanics Parameters of Strata

土層厚度/m密度/（kg·m-3）彈性模量/MPa泊松比滲透系數(shù)/（m·d-1）貯水系數(shù)/m-1孔隙度

第1層2001 900～2 190100～1500.30～0.451.800 0～2.400 00.010 00.40～0.50

第2層1002 000～2 20040～600.25～0.040.005 0～0.008 00.005 00.24～0.35

第3層1502 000～2 300110～1600.25～0.301.500 0～1.900 00.008 00.30～0.40

第4層1002 050～2 25050～650.20～0.330.004 0～0.005 50.000 10.24～0.27

第5層2502 100～2 300190～2100.17～0.221.400 0～1.600 00.007 00.30～0.40

表2開采井開采流量

Tab.2Flows of Groundwater Withdrawal from Extraction Wells

開采井編號開采井經(jīng)緯度實測流量/（m3·d-1）質量流量/（kg·（m·s）-1）

K1（37°28′10.86″N，116°17′55.53″E）-82.19-0.006 341 821

K2（37°27′3.56″N，116°18′20.69″E）-82.19-0.006 341 821

K3（37°26′43.02″N，116°21′37.20″E）-109.59-0.008 456 019

K4（37°26′43.02″N，116°21′37.21″E）-136.99-0.010 570 216

K5（37°27′58.89″N，116°19′12.09″E）-82.19-0.006 341 821

K7（37°28′8.18″N，116°18′34.75″E）-493.15-0.038 051 698

K8（37°28′3.26″N，116°18′44.72″E）-301.37-0.023 253 858

K9（37°29′39.51″N，116°19′1.54″E）-104.11-0.008 033 179

注：開采井K3已出邊界。

圖4開采井分布

Fig.4Distribution of Extraction Wells

建立的數(shù)值模型見圖5。三維網(wǎng)格共劃分為202 964個單元，計算平臺采用Comsol Multiphisices多物理場仿真分析平臺。

圖5數(shù)值模型三維網(wǎng)格剖分

Fig.53D Mesh Subdivision of Numerical Model

3.2流固耦合的實現(xiàn)及其支配方程

德州市德城區(qū)中深層承壓水的大量抽取超過承壓水的補給能力，引起承壓水位降低，也就是土層的孔隙水壓力降低導致了承壓含水層土體的變形。這是一個包含滲流與固體力學相互作用的過程，其壓縮變形具有一定的時間效應。

流固耦合支配方程主要由流體控制方程和多孔介質控制方程組成。本文采用比奧固結理論，其基本控制方程有流體控制方程和多孔介質控制方程。

（1）流體控制方程為

SymbolQC@ σ=ρg（2）

式中：σ為總應力張量；ρ為總密度；g為重力加速度。

流固耦合的實現(xiàn)過程就是在初始條件和邊界條件的基礎上對流體控制方程和多孔介質控制方程的耦合求解過程。

3.3初始條件和邊界條件

德州市德城區(qū)地面沉降監(jiān)測資料自1991年開始，取為計算初始時間點。整個模型外邊界有水頭監(jiān)控點，采用水頭邊界；模型內所含263口開采井取為流量邊界，按實際開采量施加，考慮數(shù)值模型建立及計算簡便，均轉化為質量流量（圖6）。

3.4模型識別與校驗

采用實測巖土力學參數(shù)、初始條件及邊界條件對德州市德城區(qū)地面沉降進行分析計算。圖7為1991～2007年德州市德城區(qū)計算地面沉降量分布。將圖7與圖1的實測地面沉降進行對比，計算所得整體沉降趨勢及沉降量與實測資料較為吻合。

圖6三維開采井點位置

Fig.6Location of 3D Extraction Well Points

圖71999～2007年累計地面沉降（單位：m）

Fig.7Accumulative Land Subsidence from 1999 to 2007 （Unit：m）

為進一步說明建立的數(shù)值模型合理性及正確性，從圖7所得流固耦合三維數(shù)值模擬計算結果中提取德州市德城區(qū)8個點位（平均覆蓋整個計算區(qū)域），與該8個點位的實測沉降量進行對比分析。

提取的8個點位逐年沉降模擬值及實測值對比曲線見圖8。從圖8可知，沉降量模擬值與實測值基本吻合，滿足精度需求。建立的數(shù)值模型可以較好地重現(xiàn)德州市德城區(qū)地下水開采誘發(fā)地面沉降的歷史過程。

圖8沉降監(jiān)測點實測值與模擬值對比曲線

Fig.8Comparison Curves Between Measured and Simulated Land Subsidences in the Monitoring Points

4地面沉降發(fā)展趨勢

應用建立的數(shù)值模型預測德州市德城區(qū)地面沉降在現(xiàn)狀地下水開采量情況下的地面沉降發(fā)展趨勢。

以2010年作為地面沉降預測初始時刻，地下水開采量采用同期各開采井統(tǒng)計量。圖9展示了計算得到的德州市德城區(qū)1 m沉降范圍演化。由圖9可知，研究區(qū)若保持現(xiàn)狀的地下水開采條件，其沉降范圍具有不斷擴展的趨勢。沉降中心位置基本保持穩(wěn)定，但對同一位置沉降量隨時間的推移不斷加大，由此形成的沉降漏斗范圍不斷擴大，漏斗中心高程不斷下降，因此，沉降災害具有加劇趨勢，應采取切實有效的控制措施加以處置。

圖9不同年份1 m沉降范圍演化預測（單位：m）

Fig.9Evolution Prediction of 1 m Subsidence Range in Different Years （Unit：m）

表3列出了德州市德城區(qū)代表性點位的沉降量預測值。根據(jù)數(shù)值模型預測結果，在現(xiàn)狀條件下開

表3現(xiàn)狀地下水開采條件下不同沉降監(jiān)測點地面沉降預測結果

Tab.3Predictions of Land Subsidence in Different Monitoring Points Under the Condition of Current Groundwater Exploitation

年份監(jiān)測點D62地面沉降量/mm監(jiān)測點Q1地面沉降量/mm監(jiān)測點Q6地面沉降量/mm監(jiān)測點D101地面沉降量/mm監(jiān)測點D14地面沉降量/mm監(jiān)測點DJ9地面沉降量/mm監(jiān)測點D110地面沉降量/mm監(jiān)測點D111地面沉降量/mm

20151 316.21 269.5910.8962.8503.3801.6687.5910.8

20161 345.01 297.3931.8986.5514.3823.7705.5931.8

20171 373.71 324.9952.31 009.2525.7844.6723.0952.3

20181 402.51 352.5972.91 031.7537.2865.4740.4972.9

20191 431.41 380.1993.51 053.7548.8885.6757.7993.5

20201 460.61 407.91 014.21 075.6560.5905.7775.01 014.2

20211 490.41 436.11 035.01 097.2572.3925.6792.51 035.0

20221 520.21 464.21 055.71 118.9584.2945.5809.91 055.7

20231 550.01 492.41 076.51 140.5596.0965.3827.31 076.5

20241 580.21 521.01 097.41 162.3608.0985.3845.01 097.4

20251 611.41 550.31 118.51 184.3620.21 005.5863.11 118.5

采地下水10年，總體地面沉降速率相對于2007年之前有所減緩，地面沉降中心位置基本穩(wěn)定，但累積沉降量不斷加大。尤其是在開采井分布較為集中的地段，累計沉降量增大的趨勢尤為嚴重。從開采井分布位置分析，沉降漏斗基本處于開采井集中分布區(qū)，因此，集中開采地下水將會加劇地面沉降量，導致漏斗高程加速下降，若不及時采取有效措施，地下水開采誘發(fā)的地面沉降將越發(fā)難以控制。

5地面沉降控制措施

德州市德城區(qū)非常迫切地需要進行地面沉降控制，本文提出減小地下水開采量、定期回灌地下水、調整開采井開采層位及開采井平面布局等工程措施，并采用建立的數(shù)值模型從定量分析角度評價這些措施的效果。

5.1減小地下水開采量

根據(jù)《山東省地面沉降防治規(guī)劃（2012～2020年）》，從2015年起減少德州市德城區(qū)第一和第二承壓含水層的地下水開采量。根據(jù)建立的數(shù)值模型預測未來20年地面沉降的變化情況，計算方案見表4。根據(jù)數(shù)值模擬結果，隨著地下水開采量的減少，德州市德城區(qū)地面沉降速率明顯降低（圖10）。當?shù)叵滤_采量減少30%及以上時，部分地面沉降監(jiān)測點在減少地下水開采量的最初幾年內出現(xiàn)了地面沉降趨于停止或少量回彈的現(xiàn)象，但是隨著時間推移，回彈量減少，地面沉降繼續(xù)發(fā)展。總體來說，地面沉降速率較未減少地下水開采量前明顯降低，因此，減少地下水開采量在一定程度上可以有效控制德州市德城區(qū)地面沉降。

表4減少地下水開采量計算方案

Tab.4Schemes for Reducing the Groundwater Exploitation

方案1方案2方案3方案4方案5

開采量減少10%開采量減少20%開采量減少30%開采量減少40%開采量減少50%

圖10減少地下水開采量不同計算方案下沉降監(jiān)測點累計地面沉降量對比

Fig.10Comparisons of Cumulative Land Subsidence Under Different Schemes for Reducing the Groundwater Exploitation

5.2定期回灌地下水

人工回灌是指利用工程設施將地表水注入地下含水層，以增加地下水儲量的措施。根據(jù)德州市德城區(qū)地面沉降中心分布狀況，在其中心布設15口回灌井，這些井年回灌量相同，根據(jù)總回灌量的不同，設置4種回灌方案（表5）。假設回灌的地下水均達到承壓含水層，根據(jù)計算結果，當?shù)叵滤毓嗔枯^少時，對地面沉降速率的影響較?。换毓嗑浇孛娉两盗渴艿叵滤毓嘤绊戄^大，遠離回灌井的地面效果不明顯（圖11）。

圖11不同回灌方案下沉降監(jiān)測點累計地面沉降量對比

Fig.11Comparisons of Cumulative Land Subsidence Under Various Recharge Schemes

表5地下水回灌方案

Tab.5Schemes for Groundwater Recharge

方案方案1方案2方案3方案4

總回灌量/（104m3·年-1）200400600800

5.3調整開采井開采層位

根據(jù)《德州市地面沉降地質調查與評價》，德州市德城區(qū)淺層微咸水、咸水盈余量每年為395.3×104 m3，尚處于未開發(fā)階段。由此設計了調整開采井開采層位的方案：從2015年起，潛水層每年多開采390×104 m3的水量，將第一或第二承壓含水層的地下水開采量每年減少390×104 m3。由于潛水層水質較差，可將開采的潛水層地下水作為城市綠化、消防和對水質要求不高的工業(yè)用水水源。

通過計算，調整地下水開采層位后的累計沉降量與未調整開采層位時相比有大幅度降低；同時，減少第二承壓含水層的地下水開采量比減少等量的第一承壓含水層的開采量對地面沉降的控制效果要好，而且沉降中心有所減緩。因此，在條件允許的情況下，地下水盡量從潛水層或上部第一承壓含水層開采。

5.4調整開采井平面布局

在開采量一定的情況下，為緩解地面沉降梯度大的問題，可通過調整開采井及開采量的分布來實現(xiàn)對地面沉降的控制。結合德州市的實際情況，重點針對地面沉降嚴重的沉降中心區(qū)進行開采井平面布局研究。調整開采井平面布局后，沉降中心區(qū)的開采井間距至少達到1 000 m（圖12）。假設沉降中心區(qū)開采井的開采流量重新調整，即各開采井的年開采量相同。

圖12開采井平面布局的調整

Fig.12Adjustment of Plane Layout for Extraction Wells

根據(jù)布局調整進行計算，結果顯示在總開采量不變的條件下，重新調整開采井和開采流量的布局，可以在很大程度上減小或減緩德州市德城區(qū)地面沉降速率，使沉降中心沉降梯度減小。因此，在開采量一定的情況下，為最大程度減少地面沉降，可以根據(jù)德州市德城區(qū)的實際情況，適當調整開采井的分布和各開采井的年開采量。

5.5小結

根據(jù)各工程措施所進行的計算分析結果表明，各工程措施均對控制研究區(qū)地面沉降產(chǎn)生效果，但也存在差異性。因此，可以考慮將多種方案進行綜合利用：根據(jù)當?shù)厣鐣?jīng)濟發(fā)展趨勢，合理規(guī)劃城市發(fā)展布局，這可以與地下水開采井平面布局相結合；在不同階段可以設置不同的地下水開采量上限，在枯水季和豐水期可以設置不同的量值；在科學合理評估各開采層位不同地下水開采量情況下，可以隨時間的推移適當調整不同開采層位，盡量保持各層位具有回灌回彈的特性，但在回灌過程中一定要注意水環(huán)境污染控制問題，因為地下水一旦遭到污染，治理將會非常困難。

6結語

地下水開采已導致華北平原德州區(qū)出現(xiàn)了大范圍的地面沉降，且有加劇的趨勢。已有研究成果已不能從根本上解決地面沉降控制的深層次問題。本文從地面沉降的本質出發(fā)，采用水文地質工程地質及巖土力學相結合的方法，建立山東省德州市德城區(qū)地面沉降流固耦合數(shù)值模型，成功反演該區(qū)地面沉降歷史發(fā)展情景，并對該區(qū)未來地面沉降發(fā)展趨勢進行了定量化預測；在此基礎上，提出研究區(qū)地面沉降若干控制措施，并對提出的工程措施進行評價。本次研究成果是山東省德州市地面沉降研究總結，為該區(qū)在地面沉降機理方面的進一步研究奠定了基礎。

參考文獻：

References：

[1]薛禹群，吳吉春，張云，等.長江三角洲（南部）區(qū)域地面沉降模擬研究[J].中國科學：D輯，地球科學，2008，38（4）：477492.

XUE Yuqun，WU Jichun，ZHANG Yun，et al.Simulation of Land Subsidence in the Southern Yangtze River Delta Plain[J].Science in China：Series D，Earth Sciences，2008，38（4）：477492.

[2]殷躍平，張作辰，張開軍.我國地面沉降現(xiàn)狀及防治對策研究[J].中國地質災害與防治學報，2005，16（2）：18.

YIN Yueping，ZHANG Zuochen，ZHANG Kaijun.Land Subsidence and Countermeasures for Its Prevention in China[J].The Chinese Journal of Geological Hazard and Control，2005，16（2）：18.

[3]何國榮，馮在敏，馮克印.水準測量網(wǎng)監(jiān)測在德州市地面沉降中的應用[J].山東國土資源，2012，28（6）：2830.

HE Guorong，F(xiàn)ENG Zaimin，F(xiàn)ENG Keyin，et al.The Monitoring of Land Subsidence Based on the Leveling Surveying Network in Dezhou City[J].Land and Resources in Shandong Province，2012，28（6）：2830.

[4]王貞國，王彥俊，馮守濤，等.德州市深層地下水人工回灌試驗淺探[J].山東國土資源，2005，21（6/7）：8285.

WANG Zhenguo，WANG Yanjun，F(xiàn)ENG Shoutao，et al.Study on Artificial Recharge Experiment of Deep Underground Water in Dezhou City[J].Land and Resources in Shandong Province，2005，21（6/7）：8285.

[5]孟慶峰，戴魯旗.德州市地面沉降現(xiàn)狀及防治對策[J].山東國土資源，2008，24（11）：89.

MENG Qingfeng，DAI Luqi.Land Subsidence and Countermeasures for Its Prevention in Dezhou City[J].Land and Resources in Shandong Province，2008，24（11）：89.

[6]王彥俊，劉桂儀.德州深層地下水降落漏斗相關模型建立及其應用[J].山東地質，1998，14（2）：3845.

WANG Yanjun，LIU Guiyi.Correlation Models Establishment and Its Application in Dezhou Funnel of Internal Underground Water[J].Geology of Shandong，1998，14（2）：3845.

[7]楊麗芝.德州深層地下水位降落漏斗演變機制與可調控性研究[D].北京：中國地質科學院，2009.

YANG Lizhi.From Princple and Controllingadjusting Research About Deep Groundwater Depression Cone in Dezhou[D].Beijing：Chinese Academy of Geological Sciences，2009.

[8]蔡文曉.德州市深層地下水開采與地面沉降關系研究[D].長春：吉林大學，2009.

CAI Wenxiao.Research on Relationship Between Deep Groundwater Exploitation and Land Subsidence in Dezhou[D].Changchun：Jilin University，2009.

[9]姬亞東，柴學周，劉其聲，等.大區(qū)域地下水流數(shù)值模擬研究現(xiàn)狀及存在問題[J].煤田地質與勘探，2009，37（5）：3236.

JI Yadong，CHAI Xuezhou，LIU Qisheng，et al.Current Situation and Existing Problems of Numerical Simulation of Groundwater Flow in a Largescale Area[J].Coal Geology and Exploration，2009，37（5）：3236.

[10] ZHAO C Y，ZHANG Q，DING X L，et al.Monitoring of Land Subsidence and Ground Fissures in Xian，China 20052006：Mapped by SAR Interferometry[J].Environmental Geology，2009，58（7）：15331540.

[11]駱祖江，曾峰，李穎.地下水開采與地面沉降控制三維全耦合模型研究[J].吉林大學學報：地球科學版，2009，39（6）：10801086.

LUO Zujiang，ZENG Feng，LI Ying.Study on Three Dimensional Full Coupling Model of Groundwater Exploitation and Landsubsidence Control[J].Journal of Jilin University：Earth Science Edition，2009，39（6）：10801086.

[12]劉聰，彭建兵，陳立偉，等.西安地裂縫場地地震效應分析[J].工程地質學報，2012，20（6）：979985.

LIU Cong，PENG Jianbing，CHEN Liwei，et al.Finite Element Analysis of Seismic Effects on Ground Fissures in Xian[J].Journal of Engineering Geology，2012，20（6）：979985.

[13]魏加華，崔亞莉，邵景力，等.濟寧市地下水與地面沉降三維有限元模擬[J].長春科技大學學報，2000，30（4）：376380.

WEI Jiahua，CUI Yali，SHAO Jingli，et al.3Dfinite Element Method Simulation of Groundwater and Land Subsidence in Jining City[J].Journal of Changchun University of Science and Technology，2000，30（4）：376380.

[14]萬偉峰.西安市地下水開采地面沉降數(shù)值模擬及防治方案研究[D].西安：長安大學，2008.

WAN Weifeng.Numerical Simulation of Land Subsidence Due to Groundwater Withdraw in Xian and Its Prevention Schemes[D].Xian：Changan University，2008.

[15]OH H J，AHN S C，CHOI J K，et al.Sensitivity Analysis for the GISbased Mapping of the Ground Subsidence Hazard near Abandoned Underground Coal Mines[J].Environmental Earth Sciences，2011，64（2）：347358.

[16]SUN Y H，ZHANG X D，MAO W F，et al.Mechanism and Stability Evaluation of Goaf Ground Subsidence in the Third Mining Area in Gong Changling District，China[J].Arabian Journal of Geosciences，2015，8（2）：639646.

[17]KIM K D，LEE S，OH H J.Prediction of Ground Subsidence in Samcheok City，Korea Using Artificial Neural Networks and GIS[J].Environmental Geology，2009，58（1）：7172.

[18]OUYANG Z H.The Role of Potential Soil Cavity on Ground Subsidence and Collapse in Coal Mining Area[J].Journal of Coal Science and Engineering，2010，16（3）：240245.

[19]張永偉.華北平原德州地面沉降成生機理、監(jiān)測預警與可控性研究[D].濟南：山東大學，2014.

ZHANG Yongwei.Formation Mechanism，Monitoring and Warning，Controlling Research of Subsidence of Dezhou in North China Plain[D].Jinan：Shandong University，2014.

[20]賈超，張國榮，何奔，等.華北平原德州地面沉降監(jiān)測預警系統(tǒng)建設報告[R].濟南：山東大學，2014.

JIA Chao，ZHANG Guorong，HE Ben，et al.Land Subsidence Monitoring and Warning System Construction Report in Dezhou District of North China Plain[R].Jinan：Shandong University，2014.