王城斌， 于清楊， 劉 晨， 張 超， 劉 偉， 宋澤宇

(吉林大學(xué)建設(shè)工程學(xué)院，長(zhǎng)春 130012)

地面沉降不僅危及人類(lèi)位于地表及地下的已建工程構(gòu)筑物，還會(huì)限制城市的進(jìn)一步開(kāi)發(fā)和建設(shè)。在全球面臨氣溫升高導(dǎo)致海平面上升問(wèn)題的當(dāng)下，地面沉降對(duì)沿海城市和地區(qū)影響尤甚。導(dǎo)致地面沉降的人為影響因素眾多，例如基坑、隧洞的開(kāi)挖，地表高層、超高層建筑建設(shè)，地下水的開(kāi)采等[1-7]，隨著人類(lèi)社會(huì)的日益發(fā)展，超大及特大城市不斷涌現(xiàn)，城市地面沉降日益顯著，特別是近幾十年，全世界范圍內(nèi)已有多個(gè)城市的發(fā)展受到由于過(guò)量開(kāi)采地下水所引起的地面沉降的困擾。蔡文曉認(rèn)為，在德州市，由深層地下水的超量開(kāi)采引起的地面沉降占總沉降量高達(dá)90%，且由地下水的開(kāi)采引起的地面沉降影響范圍較大[8]。

基于有效應(yīng)力原理可知，抽水引起的地面沉降是由孔隙靜水壓力降低，土體有效應(yīng)力增大所致[9]。對(duì)于抽水所引起的地面沉降，丁德民等認(rèn)為，基于比奧固結(jié)理論建立的滲流場(chǎng)和應(yīng)力場(chǎng)有限元耦合模型進(jìn)行的地面沉降計(jì)算結(jié)果與實(shí)際情況吻合度良好[10]。

大量研究表明地下構(gòu)筑物的存在會(huì)降低地下土體的滲透系數(shù)[11-14]。許燁霜分析地下?lián)跛畼?gòu)筑物對(duì)土體滲透性的影響時(shí)采用等效介質(zhì)理論計(jì)算等效滲透系數(shù)，進(jìn)而分析計(jì)算擋水構(gòu)筑物對(duì)地面沉降的影響，得到了擋水構(gòu)筑物的設(shè)置會(huì)增大地面沉降量的結(jié)論[15]。但是，這種等效介質(zhì)計(jì)算理論并未將地下構(gòu)筑物位置的影響考慮在內(nèi)，土體不同位置的水頭變化也很難準(zhǔn)確監(jiān)測(cè)或分析計(jì)算得到，且用目前廣泛采用的固結(jié)沉降理論和解析解方法計(jì)算抽水所引起的地面沉降往往沒(méi)有考慮土層的側(cè)摩阻力[16-17]，只能得到實(shí)際情況的近似解。為了分析地下構(gòu)筑物對(duì)抽水引起的地面沉降的影響和土體各參數(shù)敏感性變化規(guī)律，需得到沉降量的精確解，考慮通過(guò)數(shù)值計(jì)算方法來(lái)得到地面沉降精確解，分析擋水構(gòu)筑物對(duì)地面沉降的影響特征，并由 MATLAB編程和公式推導(dǎo)計(jì)算分析和驗(yàn)證擋水構(gòu)筑物的設(shè)置對(duì)土體參數(shù)相對(duì)于地面最大沉降量敏感性的影響。

1 抽水沉降模型的建立

1.1 基于滲透固結(jié)理論的抽水沉降力學(xué)模型

比奧從固結(jié)機(jī)理出發(fā)，考慮土體三向滲流和變形，且考慮土體變形與孔壓的相互作用，在太沙基固結(jié)理論的基礎(chǔ)上推導(dǎo)了準(zhǔn)確反應(yīng)孔隙壓力消散和土體固結(jié)關(guān)系的真三維固結(jié)方程[18](假設(shè)土體骨架微小彈性變形，變形符合壓縮定律；地下水不可壓縮，滲流符合達(dá)西定律)：

(1)

式(1)中，G為剪切模量，μ為泊松比，P為孔隙水壓力，γw為4 ℃下水的容重，u、v、w為土層x、y、z向的位移，Kx、Ky、Kz為土層x、y、z向的滲透系數(shù)，η2為拉普拉斯算子。采用有限元法，以比奧固結(jié)理論為基礎(chǔ)，通過(guò)構(gòu)建地下應(yīng)力場(chǎng)和滲流場(chǎng)耦合模型，研究擋水構(gòu)筑物對(duì)抽水引起的地面沉降影響征。

1.2 抽水沉降的幾何模型

圖1 計(jì)算模型示意圖Fig.1 Schematic diagram of the calculation model

考慮多層含水層和弱隔水層時(shí)的計(jì)算模型如圖1所示，由于一般情況下，地下土層一般為隔水層與含水層互層，由于出于規(guī)律性研究的考量，選取的土層不宜過(guò)少，但太多亦意義不大，結(jié)合文獻(xiàn)[10,15]，分析模型由含水層與弱隔水互層的11個(gè)土層組成。出于定性研究的目的，地下?lián)跛畼?gòu)筑物很少有深入到地下500 m的情況，且主要分布于從地表至地下200 m范圍內(nèi)，所以除地表土層外，其余地層選取為50 m厚，即包括1個(gè)35 m潛水層、5個(gè)厚度為50 m的承壓含水層、5個(gè)厚度為50 m的弱隔水層。由于以邊界不透水的工況為研究對(duì)象，研究對(duì)象應(yīng)選為大范圍，故模型長(zhǎng)、高分別為1 500 m和535 m，厚度1 m。初始地下水位取5 m。模型側(cè)面設(shè)為法向位移約束，底面設(shè)為全約束，表面為透水邊界，側(cè)面和底面均為不透水邊界。開(kāi)采井類(lèi)型為排井，研究范圍內(nèi)只有一個(gè)且位于模型中間，方向與寬方向垂直。地下構(gòu)筑物位于排水井左側(cè)，為了避免擋水構(gòu)筑物與周邊土體由于材質(zhì)不同而引起的應(yīng)力分布變化對(duì)擋水效應(yīng)的分析，本文計(jì)算過(guò)程中，地下構(gòu)筑物相對(duì)于土層設(shè)為不透水材料。

1.3 計(jì)算參數(shù)

現(xiàn)結(jié)合一般實(shí)際情況，以文獻(xiàn)[10,15]研究區(qū)域中各土層分布特征和相關(guān)土體參數(shù)大小為基礎(chǔ)，以達(dá)到了解地下?lián)跛畼?gòu)筑物對(duì)土體參數(shù)敏感性影響規(guī)律的研究目的。由于地表土層一般缺乏上覆荷載，故選為潛水含水層，巖性為粉砂。由于接近地表的土層，接收滲透雨水和地表水較多，第一隔水層選為黏土。起初，隨著深度增大，在承壓水的作用下，細(xì)粒顆粒被水流沖走，含水層構(gòu)架顆粒增大，故選為細(xì)砂和中細(xì)砂，深處含水層由于承壓加大，土層孔隙減小，水流變慢，含水層選為粉細(xì)砂。除接近地表的隔水層外，其余隔水層皆為粉質(zhì)黏土。由于選取模型土層較厚，為了使模型與一般實(shí)際更吻合，隔水層皆設(shè)置為弱隔水層。具體各物理力學(xué)參數(shù)見(jiàn)表1。

計(jì)算模型包括不考慮擋水構(gòu)筑物和考慮擋水構(gòu)筑物兩大類(lèi)，在不考慮擋水構(gòu)筑物時(shí)，計(jì)算開(kāi)采層分別為第1～5含水層時(shí)的地面沉降；考慮擋水構(gòu)筑物時(shí)，計(jì)算分析擋水構(gòu)筑物各種進(jìn)深工況時(shí)的影響，考慮到目前一般人類(lèi)活動(dòng)中，擋水構(gòu)筑物進(jìn)深、地下水開(kāi)采深度及開(kāi)采量等情況，設(shè)置開(kāi)采層為第2含水層，抽水速度設(shè)置為40 m3/d。

2 抽水引起的地面沉降數(shù)值計(jì)算結(jié)果分析

土的壓縮固結(jié)其實(shí)不僅與土體本身的性質(zhì)有關(guān)，也與邊界條件、抽水情況和人類(lèi)活動(dòng)等有關(guān)。條件不同，得到的結(jié)果往往有很大的差異。在地表處設(shè)置24個(gè)監(jiān)測(cè)點(diǎn)，其中抽水井附近監(jiān)測(cè)點(diǎn)較密。經(jīng)相關(guān)研究，人類(lèi)活動(dòng)10 a左右為一個(gè)周期，故以下計(jì)算結(jié)果為開(kāi)采含水層10 a后地面沉降的情況。

2.1 不設(shè)置擋水構(gòu)筑物時(shí)地面沉降數(shù)值計(jì)算結(jié)果分析

由圖2為不設(shè)置擋水構(gòu)筑物時(shí)土體沉降云圖，計(jì)算結(jié)果如圖3，在沒(méi)有擋水構(gòu)筑物的情況下，抽取地下水引起的地面沉降量隨著抽水層深度的增大而增大，開(kāi)采第五含水層10 a后的最大沉降量達(dá)0.415 5 m，是開(kāi)采第一含水層10 a后最大沉降量的2.61倍。地表沉降較為平緩，沉降差較小，最大沉降差發(fā)生在抽取第二含水層時(shí)，最大沉降差為0.041 3 m。分析認(rèn)為這是由于模型邊界不透水的設(shè)定和各含水層良好的透水性所致。

圖2 不設(shè)置擋水構(gòu)筑物時(shí)土體沉降云圖Fig.2 Cloud diagram of land subsidencewithout water retaining structure

2.2 設(shè)置擋水構(gòu)筑物條件下地面沉降數(shù)值計(jì)算結(jié)果分析

考慮擋水構(gòu)筑物深度的影響時(shí)，將擋水構(gòu)筑物設(shè)置在抽水井左側(cè)距離水井100 m處，抽取第二含水層。考慮到抽水模型的対稱(chēng)性，且抽水井左側(cè)擋水構(gòu)筑物的設(shè)置對(duì)右側(cè)的沉降量無(wú)影響，故只統(tǒng)計(jì)設(shè)置有地下?lián)跛畼?gòu)筑物的抽水井一側(cè)15個(gè)監(jiān)測(cè)點(diǎn)的沉降，分析隨著擋水構(gòu)筑物逐漸深入地下時(shí)的地表沉降特征。

表1 各土層物理力學(xué)性質(zhì)參數(shù)表Table 1 Parameters of physical and mechanical propertyof each soil layer

圖4為設(shè)置擋水構(gòu)筑物時(shí)土體沉降云圖，計(jì)算結(jié)果如圖5、圖6所示：在抽水井處，地面沉降量隨著擋水構(gòu)筑物的深入逐漸增大，當(dāng)擋水構(gòu)筑物由未完全深入第二含水層底部的225 m處深入到第二含水層底部的235 m時(shí)，即由部分阻擋到完全阻擋抽水層時(shí)，地面沉降增幅最大，增幅達(dá)21.73%，其他情況下隨著擋水構(gòu)筑物深入地面沉降量增幅都較小。擋水構(gòu)筑物另一側(cè)，距離抽水井足夠遠(yuǎn)的650 m處，如圖6，沉降量起初隨著擋水構(gòu)筑物的逐漸深入小幅增大，隨后深入到第二含水層底部時(shí)，沉降量出現(xiàn)較大幅度的減小，顯示沉降量小于未設(shè)置擋水構(gòu)筑物時(shí)，隨著擋水構(gòu)筑物的繼續(xù)深入地下，沉降量呈現(xiàn)小幅減小趨勢(shì)。結(jié)合前面介紹的抽取各含水層的沉降特點(diǎn)，結(jié)合圖3可知，在650 m處起初隨著擋水構(gòu)筑物的入深，沉降量有小幅增大是由于擋水構(gòu)筑物只擋住了淺部土層的透水，在沒(méi)有全擋住抽水層之前卻加大了深部土層的開(kāi)采量，進(jìn)而導(dǎo)致?lián)跛畼?gòu)筑物剛開(kāi)始深入地下時(shí)，650 m處沉降量小幅增大?？傮w來(lái)看，隨著擋水構(gòu)筑物的深入，地表沉降差增大。另外，在抽水井擋水構(gòu)筑物附近，地表的沉降量也是隨著擋水構(gòu)筑物的深入逐漸增大，擋水構(gòu)筑物深入到完全擋住抽水井另一側(cè)時(shí)，即擋水構(gòu)筑物深入地下535 m，在擋水構(gòu)筑物處監(jiān)測(cè)點(diǎn)，雖然在理論上地下水頭減小量是降低，但比起未設(shè)置擋水構(gòu)筑物時(shí)沉降量增大29.95%，分析表明這主要是由于抽水井一側(cè)發(fā)生的較大沉降量增大地層側(cè)摩阻力所致。

圖3 開(kāi)采不同含水層時(shí)地面沉降曲線Fig.3 Curves of land subsidence whenpumping different aquifer

圖4 擋水構(gòu)筑物深入到第二含水層底部時(shí)抽水土體沉降云圖Fig.4 Cloud diagram of land subsidence with water retainingstructure penetrating to the bottom of the second aquifer

圖5 設(shè)置擋水構(gòu)筑物一側(cè)不同擋水構(gòu)筑物入深條件下抽水引起的地面沉降曲線Fig.5 Curves of land subsidence caused by pumping water onthe side with water retaining structure which penetratinginto underground differently

圖6 隨著擋水構(gòu)筑物深入地下兩典型監(jiān)測(cè)點(diǎn)沉降曲線圖Fig.6 Curves of land subsidence of two typicalmonitoring points with water retaining structurepenetrating into underground

3 基于二次響應(yīng)面法的最大地面沉降土體參數(shù)相對(duì)敏感性和結(jié)構(gòu)可靠度分析

3.1 二次響應(yīng)面法理論

對(duì)大型復(fù)雜的結(jié)構(gòu)進(jìn)行計(jì)算分析時(shí)，采用傳統(tǒng)的解析法得到結(jié)構(gòu)的顯式函數(shù)往往很困難，而且推導(dǎo)得到的函數(shù)在應(yīng)用上一般均具有限制性和近似性的特征。這種情況下，不改變?cè)蟹匠痰暮x，在概率意義上又能代替顯示函數(shù)的多項(xiàng)式響應(yīng)面法應(yīng)運(yùn)而生[19]。多項(xiàng)式響應(yīng)面法是通過(guò)一系列數(shù)值模擬計(jì)算,采用多項(xiàng)式函數(shù)擬合未知的極限狀態(tài)面。其中分析結(jié)構(gòu)非線性函數(shù)的不含交叉項(xiàng)的二次多項(xiàng)式響應(yīng)面函數(shù)得到廣泛應(yīng)用，亦選取此種方法：

(2)

式(2)中:Xi為與函數(shù)值有關(guān)的參數(shù)樣本點(diǎn)，系數(shù)a、bi和ci由選取2n+1個(gè)樣本點(diǎn)來(lái)獲得。

3.2 結(jié)構(gòu)可靠度JC算法

地下土體的復(fù)雜性決定了土體參數(shù)往往不能準(zhǔn)確獲取，基于土體隨機(jī)變量分布參數(shù)來(lái)獲得的結(jié)構(gòu)可靠度具有較好的現(xiàn)實(shí)意義。

JC算法又稱(chēng)驗(yàn)算點(diǎn)法或改進(jìn)的一次二階矩法，是在中心點(diǎn)法的基礎(chǔ)上考慮了變量的概率分布，正態(tài)分布變量需轉(zhuǎn)化為等效正態(tài)分布[20]。采用結(jié)構(gòu)在驗(yàn)算點(diǎn)處的切平面代替極限狀態(tài)曲面，并利用改進(jìn)的一次二階矩法計(jì)算結(jié)構(gòu)可靠度。結(jié)構(gòu)極限狀態(tài)功能函數(shù)可表示為

Z=g(X1,X2,…,Xn)

(3)

(4)

(5)

(6)

由可靠度定義有

(7)

結(jié)合式(4)～式(7)可得

(8)

滿(mǎn)足這方程式的驗(yàn)算點(diǎn)為

(9)

二次響應(yīng)面法新驗(yàn)算點(diǎn)

(10)

式(10)中：Xi為JC算法初始迭代點(diǎn)。

3.3 土體參數(shù)相對(duì)敏感度計(jì)算方法

(11)

(12)

式中:UXi和σXi分別為隨機(jī)變量分布參數(shù)的均值和方差，X*為隨機(jī)變量分布參數(shù)X的驗(yàn)算點(diǎn)，CX=σXρXσX，為隨機(jī)變量分布參數(shù)X的協(xié)方差矩陣，ρX是相關(guān)系數(shù)矩陣。

3.4 最大地面沉降土體參數(shù)相對(duì)敏感度和結(jié)構(gòu)可靠度計(jì)算結(jié)果分析

選取無(wú)擋水構(gòu)筑物和設(shè)置擋水構(gòu)筑物深入地下235 m，即正好完全阻擋第二含水層兩種典型工況，分別開(kāi)采第一含水層、同時(shí)開(kāi)采第一、二含水層和開(kāi)采第二含水層的三組工況為對(duì)比研究對(duì)象。

工況一：無(wú)擋水構(gòu)筑物，開(kāi)采第一含水層。

工況二：擋水構(gòu)筑物深入地下235 m，開(kāi)采第一含水層。

工況三：無(wú)擋水構(gòu)筑物， 等量開(kāi)采第一、二含水層。

工況四：擋水構(gòu)筑物深入地下235 m， 等量開(kāi)采第一、二含水層。

工況五：無(wú)擋水構(gòu)筑物，開(kāi)采第二含水層。

工況六：擋水構(gòu)筑物深入地下235 m，開(kāi)采第二含水層。

各參數(shù)分布特征的選取結(jié)合文獻(xiàn)[22]而得，隨機(jī)變量分布參數(shù)的均值由各土層設(shè)計(jì)值加權(quán)得到，如表2所示。

表2 敏感性分析的土體參數(shù)Table 2 Soil parameters for sensitivity analysis

考慮彈性模量和滲透系數(shù)兩個(gè)最敏感參數(shù)間相關(guān)性，ρ(E,K)=-0.2,ρ(E,K)=0,ρ(E,K)=0.2，其他參數(shù)互相獨(dú)立。

由于選取抽水含水層較厚，計(jì)算結(jié)構(gòu)可靠度的最大沉降量選取不宜過(guò)小，結(jié)合規(guī)范可知，一般工程中允許最大地面沉降量為0.2 m左右,且實(shí)際情況中互層地層概化后一般為10～20 m，出于選取地層厚度及前面計(jì)算得出地表沉降量的考慮，計(jì)算參數(shù)相對(duì)敏感性和結(jié)構(gòu)度時(shí)，設(shè)置最大允許地面沉降為0.8 m。利用二次響應(yīng)面法、JC算法，MATLAB編程計(jì)算參數(shù)相對(duì)敏度和結(jié)構(gòu)可靠度的具體計(jì)算流程如圖7。

計(jì)算結(jié)果如表3所示：對(duì)比同一工況下的參數(shù)敏感度可知，影響地面沉降的參數(shù)按敏感性大小排序依次是E、K、μ、ρ、φ，其中C和e敏感性極小，E、K、μ對(duì)沉降量的敏感性占絕對(duì)地位。φ的敏感性大于C的原因推測(cè)是開(kāi)采的含水層是深層，側(cè)限土壓力較大；而e從無(wú)到有的變化可能與開(kāi)采深度有關(guān)。對(duì)比工況一和二，工況三和四，工況五和六，分析可知，當(dāng)開(kāi)采含水層相同時(shí)，E、ρ和e敏感性均隨擋水構(gòu)筑物的設(shè)置敏感性減小，而K、μ、φ和C增大。

地面沉降由剪應(yīng)變和正應(yīng)變兩項(xiàng)控制，而由剪切應(yīng)力公式：

τ=kρghtanφ+C

(13)

式(13)中:g是重力加速度，k為側(cè)壓系數(shù)，k=μ/(1-μ),由于密度ρ與正應(yīng)力相關(guān)：

圖7 土體參數(shù)相對(duì)敏感度和結(jié)構(gòu)可靠度計(jì)算流程圖Fig.7 Flow chart of parameter sensitivity andstructure reliability calculation process

表3 土體參數(shù)相對(duì)敏感度計(jì)算結(jié)果Table 3 Computing results of relative sensitivities of soil parameters

σ=ρgh

(14)

且由于發(fā)生變形沉降過(guò)程中S(σ)>S(τ),可知與剪應(yīng)變即側(cè)摩阻力控制地面沉降量正相關(guān)的參數(shù)有泊松比、摩擦角和內(nèi)聚力，由此可見(jiàn)，擋水構(gòu)筑物的設(shè)置，使得由土層側(cè)摩阻力控制的沉降比例提高，側(cè)摩阻力的重要性增大,與前述數(shù)值模擬地面沉降規(guī)律分析一致。

對(duì)比工況一、三和五、工況二、四和六無(wú)法得出明顯規(guī)律，即環(huán)境相同，開(kāi)采含水層不同時(shí)，未發(fā)現(xiàn)明顯規(guī)律。

為了驗(yàn)證敏感性參數(shù)分析的正確性，進(jìn)一步計(jì)算所有隨機(jī)變量分布參數(shù)變異和部分隨機(jī)變量分布參數(shù)變異對(duì)可靠指標(biāo)結(jié)果的影響。計(jì)算結(jié)果如表4所示。

表4 考慮不同變量的結(jié)構(gòu)可靠度Table 4 Structure reliability considering different variables

從表4可以看出，只考慮參數(shù)彈性模量、滲透系數(shù)、泊松比的變異與考慮全部參數(shù)變異得到的可靠度指標(biāo)基本一致，進(jìn)一步驗(yàn)證了這三個(gè)參數(shù)對(duì)地面沉降影響較大的分析結(jié)果。倘若工程級(jí)別較低，地面沉降可靠度的計(jì)算可視其他參數(shù)為常數(shù)，以便減小不必要的工作量。同時(shí)，由計(jì)算結(jié)果可知，考慮彈性模量與滲透系數(shù)的正相關(guān)性和設(shè)置擋水構(gòu)筑物，計(jì)算所得的可靠度指標(biāo)均減小，設(shè)置擋水構(gòu)筑物的影響與前述數(shù)值模擬結(jié)果一致。

4 考慮地面沉降側(cè)摩阻力的二維控制方程推導(dǎo)驗(yàn)證土體參數(shù)敏感性變化

地下土體是一個(gè)整體，不同位置的各點(diǎn)變形量具有關(guān)聯(lián)性和傳遞性，特別是存在沉降差的相鄰兩點(diǎn)。因而由開(kāi)采地下水而導(dǎo)致每一個(gè)土體微單元的固結(jié)壓縮量不僅與該單元孔隙水壓力消散導(dǎo)致有效應(yīng)力增加有關(guān)，當(dāng)與周?chē)鷨卧w存在沉降差時(shí)，還受到周?chē)馏w牽引而導(dǎo)致的有效應(yīng)力減小有關(guān)。

考慮側(cè)摩阻力，基于假定①各土層土體均質(zhì)、各向同性，和假定②土體只發(fā)生豎向位移、土體應(yīng)力應(yīng)變關(guān)系符合線性關(guān)系。進(jìn)行理論推導(dǎo)[23]，開(kāi)采地下水引起的土層壓縮公式為

(15)

式(15)中：n為水頭發(fā)生變化的土層數(shù)，Δdi, w是土層i水頭減小量，γw為地下水重度，mi為土層i厚度，Ei,s為彈性模量。

考慮側(cè)摩阻力徑向傳遞時(shí)，選取單位長(zhǎng)度微元體進(jìn)行分析：

(16)

式(16)中：τ、τrz為土層剪應(yīng)力，r是土體與抽水井之間的距離。

由剪應(yīng)力、剪應(yīng)變和位移的關(guān)系，可知：

τrz=Gγrz

(17)

(18)

式中：G為剪切模量，γrz為總剪應(yīng)變,ur和s分別為r和z方向上的位移。

當(dāng)忽略影響非常小的徑向位移，只考慮豎向沉降時(shí)，有：

(19)

聯(lián)立式(16)和式(19)可得：

(20)

由側(cè)摩阻力控制的沉降量為

(21)

式(21)中：M為抽水影響的地下深度。

增加的剪應(yīng)力，即側(cè)摩阻力對(duì)地面沉降有抑制作用

(22)

由于二次導(dǎo)項(xiàng)較小，此處可以不予考慮可得

(23)

所以控制方程的分析結(jié)果與敏感性參數(shù)分析結(jié)果一致，再一次驗(yàn)證了擋水構(gòu)筑物的設(shè)置會(huì)增大滲透系數(shù)K和與側(cè)摩阻力相關(guān)的參數(shù)泊松比μ、內(nèi)摩擦角φ和內(nèi)聚力C敏感性，相應(yīng)地其他參數(shù)敏感性會(huì)減小的結(jié)論。

5 結(jié)論

(1)在未設(shè)置擋水構(gòu)筑物時(shí)，開(kāi)采地下水時(shí)地表沉降差均較小，但是隨著開(kāi)采含水層深度增加，地表沉降量增大。

(2)隨著擋水構(gòu)筑物逐漸深入地下，擋水構(gòu)筑物兩側(cè)的沉降規(guī)律有所不同，開(kāi)采井一側(cè)地表沉降為逐漸增大，另一側(cè)為先增大后減小，沉降量變化幅值較?。坏徽摀跛畼?gòu)筑物深入地下多少米都會(huì)加大地表最大沉降量和沉降差，其中擋水構(gòu)筑物對(duì)開(kāi)采層的擋水效應(yīng)所引起的地面沉降最為顯著。

(3)對(duì)地面沉降有影響的參數(shù)按敏感性大小排序依次是彈性模量、滲透系數(shù)、泊松比、重度、內(nèi)摩擦角，其中內(nèi)聚力和孔隙比敏感性極小，彈性模量、滲透系數(shù)、泊松比對(duì)沉降量的敏感性占絕對(duì)地位。

(4)擋水構(gòu)筑物的設(shè)置會(huì)加大滲透系數(shù)和與側(cè)摩阻力相關(guān)參數(shù)的敏感性，與側(cè)摩阻力相關(guān)的參數(shù)有泊松比、內(nèi)摩擦角和內(nèi)聚力，相應(yīng)地其他與正應(yīng)變相關(guān)的物理力學(xué)參數(shù)，如彈性模量、密度和孔隙比敏感性減小。故當(dāng)計(jì)算分析地面沉降量時(shí)，應(yīng)該適當(dāng)適時(shí)考慮土層側(cè)摩阻力的影響。

(5)以最大沉降量為結(jié)構(gòu)可靠指標(biāo)，地下?lián)跛畼?gòu)筑物的設(shè)置和考慮彈性模量與滲透系數(shù)的正相關(guān)性計(jì)算得到的結(jié)構(gòu)可靠度均降低。土層由側(cè)摩阻力控制的沉降量大小主要取決于土層的沉降差和發(fā)生變形的土層厚度。