基于有限元方法的冕寧跨斷層水準(zhǔn)大幅變化機(jī)理

岳沖 趙靜 牛安福 吉平 李曉帆 蘇琴 楊凡

摘要：分析計(jì)算安寧河斷裂帶上冕寧跨斷層水準(zhǔn)場(chǎng)地原始觀測(cè)數(shù)據(jù)及年均速率結(jié)果，結(jié)合1-1C測(cè)段加密觀測(cè)數(shù)據(jù)的差異性特征，重點(diǎn)對(duì)滑坡、抽水2個(gè)主要影響因素，利用有限元方法構(gòu)建該場(chǎng)地三維Mohr-Coulomb彈塑性滑坡體模型及二維抽水影響的單向流固耦合模型，模擬獲得2種主要因素的影響模式。結(jié)果顯示：冕寧水準(zhǔn)1-1C測(cè)線中距離1C點(diǎn)約300 m范圍為古滑坡體所處邊坡的坡腳區(qū)域，古滑坡體對(duì)坡腳區(qū)域長(zhǎng)期的累積作用使得該區(qū)域存在隆升變形特征，進(jìn)而使1-1C測(cè)段在抽水造成的大幅下降過(guò)程中出現(xiàn)沉降量逐漸減小后再次增大的特征;地下水抽取造成沿水井形成了較為明顯的抽水漏斗沉降變形特征，在含水系統(tǒng)平衡被打破后，1號(hào)點(diǎn)所處土體邊坡更易受到抽水帶來(lái)的沉降變形影響，出現(xiàn)與井口處相當(dāng)甚至比井口附近更大的沉降變形。

關(guān)鍵詞：冕寧跨斷層水準(zhǔn);有限元法;Mohr-Coulomb;流固耦合模型

中圖分類號(hào)：P315.725文獻(xiàn)標(biāo)識(shí)碼：A文章編號(hào)：1000-0666（2020）04-0666-08

0引言

斷層作為地殼形變的薄弱地帶，強(qiáng)地震的孕育及發(fā)生往往伴隨著顯著的斷層變化（王永安等，2003，2004;岳沖等，2019）。定點(diǎn)跨斷層測(cè)量作為斷層形變監(jiān)測(cè)的有效手段，能夠直接反映測(cè)點(diǎn)處斷層段落的運(yùn)動(dòng)變化，但受布設(shè)場(chǎng)地等環(huán)境因素的影響，并不能保證所有的觀測(cè)樁均布設(shè)在基巖上，特別是位于川滇地區(qū)的跨斷層場(chǎng)地，由于地形起伏變化較大，部分觀測(cè)樁布設(shè)于山體底部土體內(nèi)，觀測(cè)樁穩(wěn)定性的差異會(huì)影響測(cè)項(xiàng)異常變化的可信程度。在實(shí)際的異常分析過(guò)程中，除了分析跨斷層場(chǎng)地所跨斷層本身調(diào)整運(yùn)動(dòng)帶來(lái)的變化以外，更多需要結(jié)合跨斷層水準(zhǔn)場(chǎng)地條件分析觀測(cè)樁周邊地下水位升降、人為載荷干擾、坡體穩(wěn)定性變化等因素的影響（馬伶俐等，2019），只有通對(duì)比分析不同影響模式下變形特征，才能結(jié)合實(shí)際異常測(cè)項(xiàng)變化，進(jìn)一步對(duì)異常信度等做出評(píng)判。

位于安寧河斷裂帶上的冕寧跨斷層水準(zhǔn)場(chǎng)地3-1測(cè)段水準(zhǔn)自2016年出現(xiàn)大幅下降變化，國(guó)內(nèi)學(xué)者進(jìn)行了大量研究工作，趙強(qiáng)和蘇琴（2018）利用BABS-InSAR技術(shù)獲取冕寧水準(zhǔn)場(chǎng)地地面沉降結(jié)果，提出該場(chǎng)地受溫泉開采影響的結(jié)論，但受測(cè)點(diǎn)所處山坡植被的影響，未能給出測(cè)線變化較大的1號(hào)點(diǎn)及所處邊坡的位移變化，其計(jì)算結(jié)果顯示沉降漏斗中心最大的沉降速率為27 mm/a，與1號(hào)點(diǎn)附近過(guò)渡點(diǎn)沉降量值相當(dāng)，與沉降漏斗沉降量向邊緣遞減的認(rèn)識(shí)并不一致;馬伶俐等（2019）針對(duì)該測(cè)項(xiàng)周邊滑坡體及抽水井等影響因素的現(xiàn)場(chǎng)落實(shí)以及數(shù)據(jù)的進(jìn)一步跟蹤分析提出了場(chǎng)地受溫泉單一因素影響的質(zhì)疑。因此針對(duì)冕寧跨斷層水準(zhǔn)場(chǎng)地大幅變化的成因及機(jī)理，本文重點(diǎn)討論滑坡及抽水2個(gè)影響因素，通過(guò)有限元方法構(gòu)建冕寧場(chǎng)地三維Mohr-Coulomb彈塑性滑坡體模型及二維抽水影響的單向流固耦合模型，獲得滑坡及抽水2種因素的影響模式，進(jìn)而對(duì)跨斷層水準(zhǔn)異常變化的成因機(jī)理進(jìn)行研究。

1研究區(qū)概況及異常分析

1.1研究區(qū)概況

本文主要研究位于安寧河斷裂帶上的冕寧跨斷層水準(zhǔn)場(chǎng)地（28.56°N，102.21°E），據(jù)衛(wèi)星影像和實(shí)地觀察，安寧河斷裂帶在該場(chǎng)地主要由安寧河?xùn)|支斷裂（主斷裂）和安寧河西支斷裂組成（圖1），安寧河西支斷裂在冶勒盆地大壩子切割全新世洪積扇，形成邊坡脊、斷塞塘等顯著的斷錯(cuò)地貌，呈近NS向延伸;安寧河?xùn)|支斷裂從紫馬垮向南經(jīng)彝海、小鹽井至馬尿河邊，總體走向近NS，沿?cái)嗔炎呦蚍植加写闋顢嗳?、斷層崖、反向陡坎等斷錯(cuò)地貌，表明其具有強(qiáng)烈的全新世活動(dòng)性。東支斷裂具有規(guī)模較大的破碎帶，冕寧水準(zhǔn)場(chǎng)地3號(hào)位于基巖破碎帶西側(cè)，其北側(cè)存在2～3條分支斷裂，向南在馬尿河南岸至秧柴溝一線，形成顯著的槽谷地貌。地震研究43卷第4期岳沖等：基于有限元方法的冕寧跨斷層水準(zhǔn)大幅變化機(jī)理

1.2冕寧跨斷層水準(zhǔn)異常分析

冕寧跨斷層水準(zhǔn)測(cè)線由3個(gè)觀測(cè)端點(diǎn)構(gòu)成（圖1c）（趙強(qiáng)，蘇琴，2018;馬伶俐等，2019），其中1號(hào)點(diǎn)位于斷層下盤，其臺(tái)基位于山坡底部，主要由松散的沖洪積塊碎石層組成，充填物以砂為主，含少量泥質(zhì)成分;2號(hào)點(diǎn)位于河流西側(cè)，臺(tái)基為昔格達(dá)組（Q1x）深灰色半固結(jié)的泥巖、粉砂巖;3號(hào)點(diǎn)位于斷層上盤，其臺(tái)基處于馬尿河右岸階地上，階地堆積物由大塊的漂石、卵石及砂組成，厚度約2 m，階地基座為花崗巖。整個(gè)場(chǎng)地長(zhǎng)約3 km（其中1-2測(cè)段測(cè)線長(zhǎng)約1.67 km，2-3測(cè)段測(cè)線長(zhǎng)約1.33 km），為單線水準(zhǔn)場(chǎng)地。整個(gè)觀測(cè)測(cè)線建于1980年11月，于1981年9月開始跨斷層短水準(zhǔn)觀測(cè)，觀測(cè)之初，每年觀測(cè)3～4個(gè)周期;1996—2012年，每年觀測(cè)6個(gè)周期;2013年4月以后，每年觀測(cè)12期，于2017年12月停測(cè)。

由圖2a原始曲線可以看出：自20世紀(jì)80年代開始觀測(cè)以來(lái)，3-2測(cè)段數(shù)據(jù)變化較為穩(wěn)定，證明2號(hào)、3號(hào)點(diǎn)穩(wěn)定性較強(qiáng)，觀測(cè)期間未發(fā)生顯著相對(duì)高差變化;相較之下，3-1測(cè)段于2006年之后整體運(yùn)動(dòng)速率有所增強(qiáng)，證明1號(hào)點(diǎn)自2006年以來(lái)出現(xiàn)下沉變化，自2016年6月開始，出現(xiàn)顯著壓縮變化，達(dá)19 mm，而3-2測(cè)段未出現(xiàn)同步異常變化，證明該測(cè)段變化可能主要由1號(hào)點(diǎn)引起。

為進(jìn)一步定量分析跨斷層水準(zhǔn)變化，排除降雨、氣溫等年度影響因素，本文通過(guò)計(jì)算斷層年均速率對(duì)比2個(gè)測(cè)段的活動(dòng)強(qiáng)度變化（江在森等，2001;岳沖等，2017）。其中，斷層的年均速率計(jì)算公式為：vt=1n∑ni=1（hti-ht-1i）（i=1，2，…，n）（1）式中：n為一年中觀測(cè)期數(shù);hti為第t年的第i個(gè)觀測(cè)值。

圖2b計(jì)算結(jié)果顯示：2016年前，3-1，3-2測(cè)段年均速率分別為-0.27，0.003 mm/a，3-1測(cè)段年均速率高于3-2測(cè)段，也證明1號(hào)測(cè)點(diǎn)的活動(dòng)性較強(qiáng)，但2個(gè)測(cè)段年均速率均在0值附近變化，也反映冕寧跨斷層水準(zhǔn)場(chǎng)地所在區(qū)域較弱的斷層垂向活動(dòng)特征。

為進(jìn)一步獲得1號(hào)點(diǎn)周圍區(qū)域內(nèi)位移下降變化，收集2016年1月以來(lái)1-1C測(cè)段（圖1c ）內(nèi)觀測(cè)過(guò)渡點(diǎn)結(jié)果，測(cè)量結(jié)果（以2016年1月測(cè)量結(jié)果為對(duì)比數(shù)據(jù)，其它月份均與此月結(jié)果對(duì)比）如圖3所示。通過(guò)差分?jǐn)?shù)據(jù)可以看出：2016年1—6月（圖3a），距離1號(hào)點(diǎn)800 m范圍內(nèi)，點(diǎn)位變化相對(duì)穩(wěn)定，而在距離1號(hào)點(diǎn)約800～1 300 m范圍內(nèi)出現(xiàn)小幅隆升-下降變化;自2016年7月起（圖3b），距離1號(hào)點(diǎn)800 m范圍內(nèi)陸續(xù)出現(xiàn)大幅沉降變化，沉降變化時(shí)段主要集中于2016年7—10月，2016年11月后整個(gè)測(cè)段沉降速率降低，其中最大的沉降點(diǎn)位于距離1號(hào)點(diǎn)64.5 m區(qū)域，最大沉降幅度達(dá)到27.5 mm;而在距離1號(hào)點(diǎn)約800～1 300 m范圍出現(xiàn)先隆升后下降變化，也證明在該區(qū)域內(nèi)可能存在其它影響因素，使得其變化與0～800 m范圍內(nèi)存在較大差異。

2計(jì)算方法及模型設(shè)計(jì)

2.1強(qiáng)度折減法

針對(duì)邊坡穩(wěn)定性的分析方法主要分為極限平衡法和有限元（或有限差分）法，而有限元法中強(qiáng)度折減法通過(guò)分析獲得安全系數(shù)，既能保持有限元模擬復(fù)雜問(wèn)題的優(yōu)勢(shì)，又能給出明確的安全系數(shù)，在邊坡穩(wěn)定性模擬計(jì)算中得到了廣泛的應(yīng)用（鄭穎人等，2005;黃盛銓等，2008;呂慶等，2008;江權(quán)等，2009;張強(qiáng)等，2011;Griffiths，Lane，1999）。

強(qiáng)度折減法最早由Zienkiewic等（1975）提出，抗剪強(qiáng)度折減系數(shù)定義為：在外載荷保持不變的情況下，邊坡內(nèi)土體所能提供的最大抗剪強(qiáng)度與外載荷在邊坡內(nèi)所產(chǎn)生的實(shí)際剪應(yīng)力之比。折減后的抗剪強(qiáng)度參數(shù)為：c′=c/Fcφ′=arctan（tanφ/Fφ）（2）式中：c，φ分別為材料的峰值粘聚力和內(nèi)摩擦角;c′，φ′分別為臨界破壞狀態(tài)下材料的粘聚力和內(nèi)摩擦角;Fc和Fφ為粘聚力和抗滑動(dòng)折減系數(shù)，通過(guò)不斷迭代模擬計(jì)算，逐漸增大折減系數(shù)，直至邊坡達(dá)到發(fā)生失穩(wěn)破裂的臨界狀態(tài)（計(jì)算不收斂的邊界）。

2.2滑坡體模擬

本文以1號(hào)點(diǎn)所處邊坡為主要研究區(qū)域設(shè)計(jì)三維滑坡體模型，模型尺寸為800 m×1 000 m，整個(gè)邊坡高度參考場(chǎng)地等高線分布為300 m（馬伶俐等，2019），邊坡角度為30°，模型示意圖見(jiàn)圖4a。本文使用有限元軟件ABAQUS完成整個(gè)模型構(gòu)建及模擬計(jì)算工作，其中采用線性六面體單位（C3D8R）進(jìn)行網(wǎng)格剖分，網(wǎng)格間距為10 m，整個(gè)模型共獲得159 200個(gè)單元，171 296個(gè)節(jié)點(diǎn)。

1號(hào)點(diǎn)臺(tái)基主要由松散的沖積碎石組成，內(nèi)含砂、泥質(zhì)等（馬伶俐等，2019），土體采用理想線彈塑性Mohr-Coulomb模型，模型的基本參數(shù)設(shè)置如下：密度為1 500 kg/m3，彈性模量為10 MPa，泊松比為0.35，粘聚力為20 kPa（陳正漢，郭楠，2019），內(nèi)摩擦角為30°（李寧等，2018;Cai，Ugal，2004）。因沉積土深度的不確定性，分別考慮沉積層厚度分別為20，40，60，80 m和全尺寸沉積土共5個(gè)模型，以探索不同厚度沉積層對(duì)地表位移的影響，其中基巖考慮為彈性體，密度為2 630 kg/m3，彈性模量為300 MPa，泊松比0.36（陳衛(wèi)兵等，2007;呂慶等，2008;許瑞寧，2016;陳正漢，郭楠，2019）。

在整個(gè)模擬計(jì)算過(guò)程中，模型的邊界條件設(shè)置如下：模型底部固定，其西側(cè)和東側(cè)邊界沿X方向固定，其北側(cè)和南側(cè)邊界沿Z方向固定，模型表面保持自由（邊界條件按照不同坐標(biāo)軸顏色進(jìn)行區(qū)分）。計(jì)算時(shí)考慮重力作用影響，重力加速度為9.8 m/s2。計(jì)算步驟如下：首先對(duì)模型施加重力載荷作用，建立模型應(yīng)力平衡狀態(tài);在應(yīng)力平衡狀態(tài)下逐漸增加強(qiáng)度折減系數(shù)（折減系數(shù)的峰值為2），迭代計(jì)算，直至數(shù)值不收斂，以此作為土坡達(dá)到臨界破壞的標(biāo)準(zhǔn)。

2.3抽水影響模擬

本文結(jié)合冕寧跨斷層水準(zhǔn)場(chǎng)地周邊遙感影像（圖1a）、等高線分布及異常落實(shí)結(jié)果（馬伶俐等，2019），構(gòu)建跨1號(hào)點(diǎn)至水井的二維剖面模型（圖4b）。根據(jù)InSAR等監(jiān)測(cè)結(jié)果顯示，沿短軸方向沉降范圍約0.4 km（趙強(qiáng)，蘇琴，2018），因此，整個(gè)模型設(shè)計(jì)長(zhǎng)600 m（水井中心距左側(cè)邊界400 m），井口處高200 m，1號(hào)點(diǎn)距井口215 m，1號(hào)點(diǎn)所處邊坡比井口處高50 m，水井深度為70 m。根據(jù)場(chǎng)地周邊大于10 m的沉積層結(jié)果（馬伶俐等，2019），井口以下20 m設(shè)置為沉積層（彈性模量5 MPa，泊松比0.4）（金小榮等，2005;曾超峰，2014）;馬伶俐等（2019）認(rèn)為該水井可能將55 m處含水細(xì)層打通，因此，井下20～55 m處設(shè)置為砂巖（彈性模型10 MPa，泊松比0.32），井口55 m以下為昔格達(dá)組地層（楊氏模量80 MPa，泊松比0.4）（孫長(zhǎng)升，2015;許瑞寧，2016;張偉松，盧海峰，2017;楊寶，2018;杜宇翔等，2020）。整個(gè)模型兩側(cè)沿水平方向固定，底部沿垂直方向固定，其上側(cè)地表均勻施加抽水引起的地表壓強(qiáng)變化，水井井底和邊坡表面設(shè)置滲透系數(shù)為0，主要模擬沉降變化最大時(shí)段2016年7—10月共計(jì)4個(gè)月抽水量引起的地表變形。

3模擬結(jié)果及分析

3.1累積位移變化分析

本文重點(diǎn)研究邊坡失穩(wěn)對(duì)跨斷層水準(zhǔn)的影響，共模擬5組不同厚度沉積層模型失穩(wěn)破壞前累積位移結(jié)果。通過(guò)在模型地表構(gòu)建剖面線（圖5a中黑線所示，從模型左側(cè)至右側(cè)點(diǎn)號(hào)依次為1～86號(hào)，其中，坡腳點(diǎn)為34號(hào)點(diǎn)，坡頂點(diǎn)為75號(hào)點(diǎn)），模擬結(jié)果顯示折減系數(shù)最大達(dá)1.16時(shí)，模型達(dá)到失穩(wěn)破壞的臨界狀態(tài)，我們獲得接近失穩(wěn)前折減系數(shù)為1.1時(shí)，5組模型地表垂向累積位移結(jié)果如圖5b所示。從模擬結(jié)果可以看出：松散沉積層厚度由20 m增加至全尺寸的沉積土模型，在滑坡體邊坡處的變形區(qū)（位移絕對(duì)值大于0 m）逐漸向坡腳處（34號(hào)點(diǎn)）擴(kuò)展，在坡頂處變形區(qū)逐步向模型右側(cè)（86號(hào)點(diǎn)）擴(kuò)展，即滑坡失穩(wěn)造成的變形范圍以及沉降位移量均顯著增大，最大位移量變化主要集中在坡頂點(diǎn)附近區(qū)域;而隨著沉積層厚度的不斷增加，失穩(wěn)帶來(lái)的坡體運(yùn)動(dòng)逐步在坡腳區(qū)域累積，進(jìn)而造成邊坡坡腳區(qū)域地表的隆升變形。

為了進(jìn)一步展示邊坡失穩(wěn)初期模型內(nèi)部位移特征，圖5a為強(qiáng)度折減系數(shù)為1.01時(shí)的全尺寸沉積土模型位移結(jié)果。圖5a顯示，隨著土邊坡整體強(qiáng)度不斷降低，坡頂區(qū)域土體逐漸失穩(wěn)并向邊坡中部區(qū)域遷移運(yùn)動(dòng)，當(dāng)強(qiáng)度折減系數(shù)達(dá)到1.16時(shí)，在邊坡中部區(qū)域形成較大垮塌位移，進(jìn)而造成滑坡運(yùn)動(dòng)，如圖1a中紫色框即為在山坡中上部區(qū)域發(fā)生的古滑坡跡象，而滑坡發(fā)生后整個(gè)坡頂區(qū)域土體向坡體中部及坡腳區(qū)域堆積，進(jìn)一步帶來(lái)坡腳區(qū)域的隆升。而本文所研究的跨斷層水準(zhǔn)1-1C測(cè)線距離1C點(diǎn)約300 m范圍即為該古滑坡體所處邊坡的坡腳區(qū)域，古滑坡體對(duì)坡腳區(qū)域的長(zhǎng)期累積作用使得該區(qū)域存在隆升變形特征，進(jìn)而使1-1C測(cè)段在大幅沉降過(guò)程中出現(xiàn)沉降量先減小后增大的跡象。而冕寧水準(zhǔn)2號(hào)點(diǎn)跨過(guò)了該滑坡區(qū)域，且2號(hào)點(diǎn)布設(shè)于半固結(jié)的昔格達(dá)組地層上，相對(duì)較為穩(wěn)定，因此滑坡體并未影響2號(hào)點(diǎn)變化，冕寧3-2測(cè)線變化相對(duì)較為穩(wěn)定。3.2抽水影響模擬結(jié)果

施加4個(gè)月抽水量模擬地表沉降的結(jié)果如圖6所示。為對(duì)比抽水量變化對(duì)地表沉降量的影響，模型模擬抽水量由10 m2/h（圖6b中模型1）增加至20 m2/h（圖6b中模型2），圖6a展示抽水量為20 m2/h時(shí)的地表沉降位移結(jié)果。從圖6a可以看出：伴隨著不斷抽水，原來(lái)平衡狀態(tài)下的含水系統(tǒng)被打破，孔隙水壓力減小，導(dǎo)致有效應(yīng)力變化，進(jìn)而造成土體變形，沿水井形成了較為明顯的抽水漏斗沉降變形，而整個(gè)左側(cè)邊坡在抽水過(guò)程中作為自然排水面，加之較為松散的沉積層結(jié)構(gòu)，存在放大沉降位移量現(xiàn)象;從具體沉降量（圖6b）可以看出，模型1中井口處沉降量達(dá)到13.3 mm，井口左側(cè)邊坡區(qū)域作為自由排水面，在整個(gè)抽水過(guò)程中坡腳區(qū)域隨水體逐步滲透排出、沉積土體固結(jié)并向坡腳區(qū)域累積，沉降量出現(xiàn)先減小后逐漸增加趨勢(shì)，最大達(dá)到15.8 mm后逐步穩(wěn)定;而當(dāng)抽水量達(dá)到20 m2/h時(shí)，整個(gè)區(qū)域沉降量明顯增加，同時(shí)整個(gè)邊坡沉降變化速率也顯著增大，其中井口處沉降量達(dá)到26.6 mm，坡腳最小沉降量為24.7 mm，1號(hào)點(diǎn)區(qū)域沉降量達(dá)到31.6 mm。

因此，作為土體邊坡上的1號(hào)點(diǎn)，因受邊坡結(jié)構(gòu)本身以及沉積層的影響，在整個(gè)抽水過(guò)程中，更易受到抽水帶來(lái)的沉降變形，出現(xiàn)與井口處相當(dāng)甚至比井口附近更大的沉降變形，但后期隨著抽水量的穩(wěn)定以及地表水、地下水的不斷補(bǔ)給，地下含水系統(tǒng)將再次達(dá)到平衡狀態(tài)，地下水水位下降一定幅度后將趨于穩(wěn)定，進(jìn)而地表沉降位移也將趨于穩(wěn)定。

4結(jié)論

本文針對(duì)冕寧跨斷層水準(zhǔn)3-1測(cè)線出現(xiàn)的大幅沉降變化異常，計(jì)算并分析了冕寧場(chǎng)地年均變化速率、測(cè)段加密數(shù)據(jù)，結(jié)合異常落實(shí)過(guò)程中的滑坡及抽水2個(gè)主要因素，利用有限元方法構(gòu)建冕寧場(chǎng)地Mohr-Coulomb彈塑性滑坡體模型及抽水影響的單向流固耦合模型，通過(guò)模擬獲得2種主要因素的影響模式，結(jié)果顯示：

（1）斷層年均速率結(jié)果顯示3-1，3-2測(cè)段年均速率分別為-0.27和0.003 mm/a，3-1測(cè)段運(yùn)動(dòng)速率高于3-2測(cè)段，也證明1號(hào)點(diǎn)活動(dòng)性較強(qiáng);1-1C測(cè)段加密觀測(cè)結(jié)果顯示沉降變形主要集中于2016年7—10月，最大的沉降點(diǎn)位于距離冕寧1號(hào)點(diǎn)64.5 m點(diǎn)位區(qū)域，最大沉降幅度達(dá)到27.5 mm，1-1C測(cè)段沉降量呈“增加—減小—增加—減小”的變化特征。

（2）模擬5組不同厚度沉積層模型失穩(wěn)破壞前累積位移結(jié)果顯示：隨松散沉積層厚度的不斷增加，失穩(wěn)變形的范圍逐步擴(kuò)大，邊坡沉降位移也逐漸增大，最大位移變化主要集中在坡頂點(diǎn)附近區(qū)域，失穩(wěn)帶來(lái)的坡體運(yùn)動(dòng)逐步在坡腳區(qū)域形成累積，進(jìn)而造成邊坡坡腳區(qū)域地表的隆升變形。冕寧水準(zhǔn)1-1C測(cè)線中距離1C點(diǎn)約300 m范圍為古滑坡體所處邊坡的坡腳區(qū)域，古滑坡體對(duì)坡腳區(qū)域長(zhǎng)期的累積作用使得該區(qū)域存在隆升變形特征，進(jìn)而使1-1C測(cè)段在抽水造成的大幅下降過(guò)程中出現(xiàn)沉降量逐漸減小后再次增大的跡象。

（3）模擬4個(gè)月抽水造成地表沉降的結(jié)果顯示：沿水井形成了較為明顯的抽水漏斗沉降變形特征，1號(hào)點(diǎn)所處土體邊坡，因本身自由排水面以及松散沉積層的影響，在含水系統(tǒng)平衡被打破后，更易受到抽水帶來(lái)的沉降變形影響，出現(xiàn)與井口處相當(dāng)甚至比井口附近更大的沉降變形。

受該地區(qū)地下水文資料、水井鉆孔資料等限制，本文對(duì)地下介質(zhì)分布、物性參數(shù)等進(jìn)行必要的簡(jiǎn)化，模擬得到的造成1號(hào)點(diǎn)沉降變化的抽水量介于10～20 m2/h，與實(shí)地異常落實(shí)的10 m2/h結(jié)果（馬伶俐等，2019）基本一致，通過(guò)定量模擬滑坡與抽水兩種影響因素的影響模式，為跨斷層測(cè)量出現(xiàn)大幅變化后的判定提供定性分析的依據(jù)，未來(lái)隨著該區(qū)域研究資料的進(jìn)一步完善，作者也將進(jìn)一步跟蹤補(bǔ)充該區(qū)域的研究結(jié)果。

兩位審稿人對(duì)本文進(jìn)行認(rèn)真的審閱并提出了寶貴建議，四川省地震局測(cè)繪工程院提供高質(zhì)量的測(cè)量數(shù)據(jù)，在此表示衷心感謝。

參考文獻(xiàn)：

陳衛(wèi)兵，鄭穎人，雷文杰.2007.沉埋樁加固滑坡體模型試驗(yàn)的三維有限元模擬[J].巖土力學(xué)，28（增刊1）：273-276.

陳正漢，郭楠.2019.非飽和土與特殊土力學(xué)及工程應(yīng)用研究的新進(jìn)展[J].巖土力學(xué)，40（1）：1-54.

杜宇翔，盛謙，王帥，等.2020.昔格達(dá)組半成巖微觀結(jié)構(gòu)與力學(xué)性質(zhì)研究[J].巖土力學(xué)，41（4）：1-13.

黃盛銓，劉君，孔憲京.2008.強(qiáng)度折減DDA法及其在邊坡穩(wěn)定分析中的應(yīng)用[J].巖石力學(xué)與工程學(xué)報(bào)，（增刊1）：2799-2806.

江權(quán)，馮夏庭，向天兵.2009.基于強(qiáng)度折減原理的地下洞室群整體安全系數(shù)計(jì)算方法探討[J].巖土力學(xué)，30（8）：2483-2488.

江在森，丁平，王雙緒，等.2001.中國(guó)西部大地形變監(jiān)測(cè)與地震預(yù)報(bào)[M].北京：地震出版社.

金小榮，俞建霖，祝哨晨，等.2005.基坑降水引起周圍土體沉降性狀分析[J].巖土力學(xué)，26（10）：54-60.