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

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

關(guān)鍵詞：地下水；地面沉降；流固耦合；水文地質(zhì)；工程地質(zhì)；巖土力學(xué)；華北平原

中圖分類號：P641；P642.26文獻(xiàn)標(biāo)志碼：A

0引言

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

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

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

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ù)下降趨勢，并伴隨嚴(yán)重的地面沉降問題。該區(qū)自1989年開始進(jìn)行地面精密測量，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)場及室內(nèi)試驗(yàn)

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

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

3數(shù)值模型

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

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

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

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

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

Fig.2Geological Record of Field Core Hole

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

圖3工程地質(zhì)特征層劃分

Fig.3Division of Engineering Geological Feature Layers

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

表1地層物理力學(xué)參數(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)緯度實(shí)測流量/（m3·d-1）質(zhì)量流量/（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流固耦合的實(shí)現(xiàn)及其支配方程

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

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

（1）流體控制方程為

SymbolQC@ σ=ρg（2）

式中：σ為總應(yīng)力張量；ρ為總密度；g為重力加速度。

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

3.3初始條件和邊界條件

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

3.4模型識別與校驗(yàn)

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

圖6三維開采井點(diǎn)位置

Fig.6Location of 3D Extraction Well Points

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

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

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

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

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

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

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

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

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

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

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

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

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

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

年份監(jiān)測點(diǎn)D62地面沉降量/mm監(jiān)測點(diǎn)Q1地面沉降量/mm監(jiān)測點(diǎn)Q6地面沉降量/mm監(jiān)測點(diǎn)D101地面沉降量/mm監(jiān)測點(diǎn)D14地面沉降量/mm監(jiān)測點(diǎn)DJ9地面沉降量/mm監(jiān)測點(diǎn)D110地面沉降量/mm監(jiān)測點(diǎ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)定，但累積沉降量不斷加大。尤其是在開采井分布較為集中的地段，累計沉降量增大的趨勢尤為嚴(yán)重。從開采井分布位置分析，沉降漏斗基本處于開采井集中分布區(qū)，因此，集中開采地下水將會加劇地面沉降量，導(dǎo)致漏斗高程加速下降，若不及時采取有效措施，地下水開采誘發(fā)的地面沉降將越發(fā)難以控制。

5地面沉降控制措施

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

5.1減小地下水開采量

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

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

Tab.4Schemes for Reducing the Groundwater Exploitation

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

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

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

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

5.2定期回灌地下水

人工回灌是指利用工程設(shè)施將地表水注入地下含水層，以增加地下水儲量的措施。根據(jù)德州市德城區(qū)地面沉降中心分布狀況，在其中心布設(shè)15口回灌井，這些井年回灌量相同，根據(jù)總回灌量的不同，設(shè)置4種回灌方案（表5）。假設(shè)回灌的地下水均達(dá)到承壓含水層，根據(jù)計算結(jié)果，當(dāng)?shù)叵滤毓嗔枯^少時，對地面沉降速率的影響較??；回灌井附近地面沉降量受地下水回灌影響較大，遠(yuǎn)離回灌井的地面效果不明顯（圖11）。

圖11不同回灌方案下沉降監(jiān)測點(diǎ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調(diào)整開采井開采層位

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

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

5.4調(diào)整開采井平面布局

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

圖12開采井平面布局的調(diào)整

Fig.12Adjustment of Plane Layout for Extraction Wells

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

5.5小結(jié)

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

6結(jié)語

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

參考文獻(xiàn)：

References：

[1]薛禹群，吳吉春，張云，等.長江三角洲（南部）區(qū)域地面沉降模擬研究[J].中國科學(xué)：D輯，地球科學(xué)，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].中國地質(zhì)災(zāi)害與防治學(xué)報，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]何國榮，馮在敏，馮克印.水準(zhǔn)測量網(wǎng)監(jiān)測在德州市地面沉降中的應(yīng)用[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]王貞國，王彥俊，馮守濤，等.德州市深層地下水人工回灌試驗(yàn)淺探[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]王彥俊，劉桂儀.德州深層地下水降落漏斗相關(guān)模型建立及其應(yīng)用[J].山東地質(zhì)，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]楊麗芝.德州深層地下水位降落漏斗演變機(jī)制與可調(diào)控性研究[D].北京：中國地質(zhì)科學(xué)院，2009.

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

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

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

[9]姬亞東，柴學(xué)周，劉其聲，等.大區(qū)域地下水流數(shù)值模擬研究現(xiàn)狀及存在問題[J].煤田地質(zhì)與勘探，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].吉林大學(xué)學(xué)報：地球科學(xué)版，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]劉聰，彭建兵，陳立偉，等.西安地裂縫場地地震效應(yīng)分析[J].工程地質(zhì)學(xué)報，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ì)寧市地下水與地面沉降三維有限元模擬[J].長春科技大學(xué)學(xué)報，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].西安：長安大學(xué)，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]張永偉.華北平原德州地面沉降成生機(jī)理、監(jiān)測預(yù)警與可控性研究[D].濟(jì)南：山東大學(xué)，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)測預(yù)警系統(tǒng)建設(shè)報告[R].濟(jì)南：山東大學(xué)，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.