三峽水庫水位下降速率與馬家溝滑坡穩定性的耦合關系探究

三峽水庫水位下降速率與馬家溝滑坡穩定性的耦合關系探究一、引言1.1研究背景與意義三峽工程作為世界上規模最大的水利樞紐工程之一，其建設具有防洪、發電、航運、水資源利用等巨大的綜合效益，對中國的經濟發展和社會穩定起到了至關重要的作用。從防洪角度看，長江流域洪水頻發，歷史上曾給中下游地區帶來巨大的災害和損失。三峽大壩建成后，憑借其巨大庫容所提供的調蓄能力，有效削減洪峰，大大減輕了洪水對中下游地區的威脅，保護了人民生命財產安全和社會穩定。據統計，三峽水利工程建成后，可抵御百年一遇甚至萬年一遇的洪水，最大洪水承受量達120000m3/秒以上，能消減的洪峰流量最高達3.3m3/秒。在發電方面，三峽水電站裝機容量巨大，是中國西電東送工程中線的巨型電源點，所發電力主要售予華中、華東和南方電網的多個省市，有效緩解了我國電力供應緊張局面。截至2022年12月31日，公司2022年第四季度總發電量130.35億千瓦時，較上年同期增長43.07%。在航運上，三峽大壩渠化了三斗坪至重慶河段，增加了葛洲壩水利樞紐，結束了“自古川江不夜航”的歷史，提升了通航能力，保障了運輸安全。與蓄水前相比，蓄水后三峽庫區年均事故件數、死亡人數、沉船數和直接經濟損失分別下降了72%、81%、65%、20%。然而，三峽水庫的運行也帶來了一系列環境地質問題，其中庫岸滑坡問題尤為突出。馬家溝滑坡位于重慶市奉節縣境內，是三峽庫區著名的滑坡地質災害之一，面積約1.5平方公里，為黃土石積體滑坡。滑坡體呈魚鱗狀，南北向展布，頂部海拔約725米，最高處比滑坡頭高出20余米，平均坡度10度，最大高度達到150米，發育于第四紀黃土和更新世高洞段懸崖的受侵蝕基礎之上，向水庫方向滑去。三峽水庫水位的周期性變化，特別是水位下降過程，對馬家溝滑坡的穩定性產生了顯著影響。當三峽水庫水位下降時，滑坡體所受水位下降的作用力將減小，最終會達到穩定狀態。但在水位下降過程中，滑坡體內部的滲流壓力也會隨之減小，進而導致巖土體強度參數發生變化，對滑坡穩定性產生影響。此外，滑坡體周邊的地下水位下降也會改變滑坡體內部滲流的方向和速度，進一步影響滑坡的穩定性。如在一些庫區，由于庫水位驟降，導致斜坡體地下水位下降相對滯后，坡內產生超孔隙水壓力一時難以消散，最終引發滑坡事故。研究三峽水庫水位下降速率對馬家溝滑坡穩定性的影響具有重大的現實意義。準確掌握二者之間的關系，能夠為馬家溝滑坡及類似地質災害的預測預警提供科學依據，提前采取有效的防范措施，避免或減少滑坡災害造成的人員傷亡和財產損失，保障三峽庫區人民的生命財產安全。這一研究成果還可為三峽庫區后續的工程建設、土地利用規劃等提供決策支持，優化工程設計和建設方案，提高工程的安全性和穩定性，促進三峽庫區的可持續發展。1.2國內外研究現狀在國外，水庫水位變化對滑坡穩定性影響的研究起步較早。1963年意大利瓦依昂水庫發生的大型滑坡災害，導致1925人死亡，此次災害引起了國際學術界對水庫滑坡問題的高度關注。此后，眾多學者圍繞庫水位變化與滑坡穩定性的關系展開研究。Jones調查了Roosevelt湖附近1941-1953年發生的滑坡，發現30%的滑坡發生在水位驟降10-20m的情況，揭示了水位驟降與滑坡發生之間的緊密聯系。一些學者運用數值模擬方法，如有限元、有限差分等，對庫水位變化下的滑坡滲流場和應力應變場進行分析，深入研究滑坡穩定性的演變規律。國內對水庫水位變化影響滑坡穩定性的研究也取得了豐碩成果。三峽庫區作為滑坡災害頻發區域，受到了國內學者的重點關注。眾多學者針對三峽庫區不同類型的滑坡，開展了大量的現場監測、物理模擬和數值分析工作。在現場監測方面，對滑坡體的位移、地下水位、孔隙水壓力等參數進行長期監測，獲取了豐富的第一手資料，為研究滑坡穩定性變化提供了數據支持。物理模擬實驗通過構建相似模型，模擬庫水位升降過程，直觀地觀察滑坡的變形破壞特征，研究其破壞機制。數值分析則借助專業軟件，如Geo-Slope、FLAC等，對滑坡在庫水位變化下的穩定性進行定量評價，分析各種因素對滑坡穩定性的影響程度。針對馬家溝滑坡，已有研究主要集中在滑坡的工程地質特征、穩定性評價以及治理措施等方面。通過地質勘查，詳細了解了馬家溝滑坡的地形地貌、地層巖性、地質構造等特征；采用極限平衡法、數值模擬法等對滑坡穩定性進行了評價，并提出了相應的治理方案，如排水工程、支擋工程等。然而，專門針對三峽水庫水位下降速率對馬家溝滑坡穩定性影響的研究相對較少。雖然部分研究涉及到庫水位變化對滑坡穩定性的影響，但大多未深入探討水位下降速率這一關鍵因素，未能明確不同下降速率下馬家溝滑坡穩定性的變化規律以及可能引發滑坡失穩的臨界水位下降速率。在研究方法上，也存在一定的局限性，缺乏多學科交叉融合的研究，難以全面、準確地揭示水位下降速率與滑坡穩定性之間的復雜關系。因此，開展三峽水庫水位下降速率對馬家溝滑坡穩定性影響的研究具有重要的理論和現實意義，有助于填補這一領域的研究空白，為馬家溝滑坡的防治提供更科學、更精準的依據。1.3研究內容與方法1.3.1研究內容本研究聚焦三峽水庫水位下降速率對馬家溝滑坡穩定性的影響，具體研究內容如下：馬家溝滑坡工程地質條件分析：通過地質測繪、鉆探、物探等方法，詳細查明馬家溝滑坡的地形地貌、地層巖性、地質構造、水文地質條件等。分析滑坡體的物質組成、結構特征以及巖土體物理力學性質，為后續研究提供基礎資料。例如，確定滑坡體中不同土層的厚度、分布范圍，以及各土層的抗剪強度、滲透系數等參數。三峽水庫水位下降規律研究：收集三峽水庫歷年水位數據，分析水位下降的時間、幅度、速率等變化規律。結合三峽水庫的運行調度方案，預測未來不同時期水位下降的情況。研究不同水位下降速率的組合工況，為數值模擬和穩定性分析提供依據。水位下降速率對滑坡滲流場的影響研究：基于滲流理論，建立馬家溝滑坡在三峽水庫水位下降條件下的滲流模型。運用數值模擬軟件，如SEEP/W，分析不同水位下降速率下滑坡體內滲流場的變化特征，包括孔隙水壓力分布、滲流速度和方向等。研究孔隙水壓力消散規律，以及滲流場變化對滑坡穩定性的影響機制。水位下降速率對滑坡穩定性的影響研究：采用極限平衡法和數值模擬法，如Geo-Slope、FLAC等軟件，計算不同水位下降速率下馬家溝滑坡的穩定性系數。分析水位下降速率與滑坡穩定性系數之間的定量關系，確定影響滑坡穩定性的關鍵水位下降速率值。研究滑坡在不同水位下降速率下的變形破壞模式，揭示水位下降速率對滑坡穩定性的影響規律。滑坡穩定性評價及預警指標研究：綜合考慮馬家溝滑坡的工程地質條件、水位下降速率對滑坡穩定性的影響，建立滑坡穩定性評價模型。確定滑坡穩定性的評價標準，對馬家溝滑坡在不同水位下降工況下的穩定性進行評價。結合穩定性評價結果，研究滑坡失穩的預警指標，為滑坡災害的預測預警提供科學依據。1.3.2研究方法本研究擬采用以下研究方法：地質調查法：對馬家溝滑坡進行現場地質調查，包括地質測繪、鉆探、物探等。通過地質測繪，詳細記錄滑坡的地形地貌、地層巖性、地質構造等信息；利用鉆探獲取滑坡體不同深度的巖土樣本，進行室內物理力學性質試驗；采用物探方法，如電阻率法、地震波法等，探測滑坡體內部的結構和地質異常，全面了解滑坡的工程地質條件。資料收集與分析：收集三峽水庫歷年水位數據、運行調度方案等資料，分析水位下降規律。查閱馬家溝滑坡的相關研究文獻、地質勘查報告等，了解滑坡的基本情況和已有研究成果，為研究提供數據支持和理論參考。數值模擬法：運用專業的數值模擬軟件，如SEEP/W、Geo-Slope、FLAC等，建立馬家溝滑坡在三峽水庫水位下降條件下的滲流模型和穩定性分析模型。通過數值模擬，分析不同水位下降速率下滑坡滲流場和穩定性的變化特征，預測滑坡的變形破壞趨勢，深入研究水位下降速率對滑坡穩定性的影響機制。理論分析方法：基于巖土力學、滲流理論等相關學科的基本原理，對水位下降速率影響滑坡穩定性的機制進行理論分析。推導孔隙水壓力、滲透力等因素與滑坡穩定性之間的數學關系，為數值模擬和穩定性評價提供理論依據。二、相關理論基礎2.1滑坡穩定性評價理論滑坡穩定性評價是研究滑坡問題的核心內容，其評價結果對于判斷滑坡是否處于穩定狀態、預測滑坡的發展趨勢以及制定相應的防治措施具有至關重要的指導意義。目前，滑坡穩定性評價方法眾多，每種方法都基于不同的理論基礎和假設條件，在實際應用中各有其優缺點和適用范圍。其中，極限平衡法和數值分析法是兩種應用較為廣泛且具有代表性的評價方法。極限平衡法歷史悠久，理論成熟，通過對滑坡體進行力學分析，求解滑體上的抗滑力與下滑力之間的平衡關系，進而得出滑坡的穩定性系數，以直觀地反映滑坡的穩定程度。數值分析法借助計算機技術，能夠模擬滑坡體在復雜條件下的應力應變狀態和滲流特性，更全面地考慮各種因素對滑坡穩定性的影響。2.1.1極限平衡法極限平衡法是基于剛體力學的靜力平衡原理，將滑坡體視為剛體，通過分析滑體上的各種作用力，建立力和力矩的平衡方程，求解滑坡的穩定性系數。該方法的基本假設是滑坡體處于極限平衡狀態，滑面上的抗滑力與下滑力達到平衡。其核心思想是化整為零，即將滑坡滑體進行條塊劃分，研究條塊之間的相互作用，不同的極限平衡法之間的差異就在于條塊間相互作用假定的不同。在實際應用中，首先要確定滑坡的潛在滑動面，這通常需要結合地質勘查資料、工程經驗以及一定的試算方法來確定。然后，將滑坡體沿滑動面劃分為若干個條塊，對每個條塊進行受力分析。作用在條塊上的力主要有重力、滑面反力、條間力等。根據靜力平衡條件，建立每個條塊的力平衡方程和力矩平衡方程。通過聯立這些方程，求解出滑坡的穩定性系數。穩定性系數的定義通常為滑面上的總抗滑力與總下滑力之比，當穩定性系數大于1時，表明滑坡處于穩定狀態；當穩定性系數等于1時，滑坡處于極限平衡狀態；當穩定性系數小于1時，滑坡處于不穩定狀態。在馬家溝滑坡穩定性分析中，極限平衡法得到了廣泛應用。通過對馬家溝滑坡的地質勘查，獲取了滑坡體的地形地貌、地層巖性、地質構造等信息，為確定潛在滑動面提供了依據。將滑坡體劃分為多個條塊，考慮條塊間的作用力，運用極限平衡法中的瑞典條分法、畢肖普法等進行計算。瑞典條分法假定土條兩側的條間力的合力近似相等、方向相反且作用在同一作用面上，不計條間力影響。畢肖普法在不考慮條塊間切向力的前提下，滿足力多邊形閉合條件，考慮了條塊間水平力的作用。通過這些方法的計算，得到了馬家溝滑坡在不同工況下的穩定性系數，為評估滑坡的穩定性提供了重要參考。極限平衡法具有概念清晰、計算簡便、能直接給出穩定性系數等優點，在工程實踐中應用廣泛，積累了豐富的經驗。然而，該方法也存在一些局限性。它假定土體為理想塑性材料，不考慮土體本身的應力-應變關系，無法反映滑坡體的變形過程和破壞機制。在確定潛在滑動面時，往往需要憑借經驗和試算，具有一定的主觀性，不同的滑動面假設可能會導致計算結果的差異。條塊間作用力的假定也存在一定的不確定性，這可能會影響計算結果的準確性。2.1.2數值分析法數值分析法是利用計算機技術，通過建立數學模型來模擬滑坡體的力學行為和滲流過程。在滑坡穩定性研究中，常用的數值分析方法包括有限元法、有限差分法、離散元法等。有限元法的基本原理是將連續的滑坡體離散為有限個單元，通過對每個單元進行力學分析，建立單元的剛度矩陣，然后將所有單元的剛度矩陣組裝成整體剛度矩陣，求解整體平衡方程，得到滑坡體的應力、應變和位移等參數。有限差分法是將滑坡體的控制方程進行離散化處理，用差分方程近似代替微分方程，通過迭代求解差分方程來獲得滑坡體的力學響應。離散元法適用于分析非連續介質的力學行為，它將滑坡體視為由離散的塊體組成，考慮塊體間的相互作用和運動，能夠較好地模擬滑坡的變形破壞過程。在研究三峽水庫水位下降速率對馬家溝滑坡穩定性的影響時，數值分析法發揮了重要作用。運用有限元軟件，如Geo-Slope中的SEEP/W模塊和SLOPE/W模塊，建立馬家溝滑坡在三峽水庫水位下降條件下的滲流模型和穩定性分析模型。在滲流模型中，考慮了水位下降速率、巖土體滲透系數等因素，通過求解滲流控制方程，得到滑坡體內不同時刻的孔隙水壓力分布和滲流速度。將滲流分析結果作為荷載施加到穩定性分析模型中，結合巖土體的本構模型，計算滑坡在不同水位下降速率下的穩定性系數和變形情況。通過數值模擬，可以直觀地觀察到水位下降過程中滑坡體內部滲流場和應力應變場的變化，以及滑坡的變形破壞模式，深入分析水位下降速率對滑坡穩定性的影響機制。數值分析法能夠考慮滑坡體的幾何形狀、材料特性、邊界條件等復雜因素，更真實地模擬滑坡的實際情況，為滑坡穩定性分析提供了更全面、準確的信息。它可以分析滑坡在不同工況下的力學響應，預測滑坡的發展趨勢，為制定合理的防治措施提供科學依據。然而，數值分析法也存在一些缺點。建立數值模型需要準確的地質參數和邊界條件，這些參數的獲取往往較為困難，且存在一定的誤差，會影響模擬結果的準確性。數值計算過程復雜，計算量較大，需要較高的計算機性能和計算時間，對計算人員的專業水平要求也較高。2.2滲流理論2.2.1達西定律達西定律由法國工程師亨利?達西（HenryDarcy）于1856年通過實驗得出，是滲流理論的基礎。該定律表明，在層流狀態下，水在多孔介質中的滲流速度與水力坡度成正比，其數學表達式為：v=ki，其中v為滲流速度，k為滲透系數，i為水力坡度。滲透系數k是反映巖土體透水性的重要參數，它與巖土體的顆粒大小、形狀、孔隙率以及流體的物理性質等因素有關。例如，粗顆粒的砂土，其孔隙較大，水流通過較為順暢，滲透系數通常較大；而細顆粒的黏土，孔隙細小且連通性差，滲透系數則較小。水力坡度i則表示單位滲流長度上的水頭損失，它反映了滲流過程中能量的消耗情況。在研究滑坡體內滲流問題時，達西定律具有重要的作用。它為分析滑坡體內的滲流場提供了基本的理論依據，能夠幫助我們定量地描述水在滑坡體中的流動規律。通過達西定律，可以計算出滑坡體內不同位置的滲流速度和水力坡度，進而了解滲流對滑坡體穩定性的影響。當三峽水庫水位下降時，滑坡體內部的水力坡度會發生變化，根據達西定律，滲流速度也會相應改變。滲流速度的變化會導致孔隙水壓力的分布發生改變，而孔隙水壓力是影響滑坡穩定性的關鍵因素之一。若滲流速度增大，孔隙水壓力消散加快，可能會使滑坡體的有效應力增加，從而提高滑坡的穩定性；反之，若滲流速度減小，孔隙水壓力難以消散，可能會導致滑坡體的有效應力減小，降低滑坡的穩定性。在實際應用中，結合地質勘查獲取的巖土體滲透系數等參數，運用達西定律進行計算，能夠為評估三峽水庫水位下降對馬家溝滑坡滲流場和穩定性的影響提供重要的數據支持。2.2.2非飽和滲流理論非飽和滲流理論是研究非飽和多孔介質中水分運動規律的理論。在非飽和狀態下，多孔介質中不僅存在液態水，還存在氣態水，且孔隙中部分被水占據，部分被空氣占據。與飽和滲流相比，非飽和滲流具有更為復雜的特性，其滲透系數不再是一個常數，而是與基質吸力密切相關。基質吸力是指非飽和土中孔隙水承受的負壓力，它反映了土中水被土顆粒吸附的程度。隨著基質吸力的增加，土中孔隙水的含量逐漸減少，孔隙的連通性變差，滲透系數也隨之減小。非飽和滲流的控制方程通常基于質量守恒定律和達西定律建立，考慮了基質吸力對滲透系數的影響。對于理解水位下降時滑坡體內滲流變化，非飽和滲流理論具有重要意義。三峽水庫水位下降過程中，滑坡體由飽和狀態逐漸轉變為非飽和狀態，非飽和滲流理論能夠更準確地描述這一過程中滲流的變化規律。當水位下降時，滑坡體上部首先變為非飽和區，基質吸力逐漸增大，滲透系數減小。這會導致孔隙水壓力的消散速度變慢，滑坡體內部的滲流場發生復雜的變化。考慮非飽和滲流特性，能夠更全面地分析水位下降對滑坡穩定性的影響。在水位下降初期，由于非飽和區的存在，孔隙水壓力的消散滯后于水位下降，可能會在滑坡體內產生較高的孔隙水壓力，降低滑坡的穩定性。而在水位下降后期，隨著非飽和區的擴大和基質吸力的進一步增加，滲透系數進一步減小，滲流速度降低，孔隙水壓力逐漸消散，滑坡的穩定性可能會有所提高。因此，運用非飽和滲流理論對三峽水庫水位下降時馬家溝滑坡體內的滲流變化進行研究，能夠為準確評估滑坡的穩定性提供更科學的依據。三、三峽水庫水位變化特征3.1三峽水庫概況三峽水庫位于長江西陵峽中段的湖北省宜昌市三斗坪鎮，是三峽水電站建成后蓄水形成的人工湖泊，是三峽工程的重要組成部分。其壩址控制流域面積100萬平方千米，約占長江流域面積的56%。水庫全長600余千米，平均寬度1.1千米，總面積1084平方千米。總庫容393億立方米，其中防洪庫容221.5億立方米。三峽水庫集防洪、發電、航運、水資源利用等巨大綜合效益于一身，對長江流域的經濟發展和生態平衡具有不可替代的作用。在防洪方面，三峽水庫憑借其龐大的庫容，能夠有效攔蓄長江上游的洪水，削減洪峰流量，減輕中下游地區的防洪壓力。當長江上游發生洪水時，三峽水庫可以根據洪水來量和中下游地區的防洪需求，合理調整庫水位，將多余的洪水蓄存在庫內，避免洪水對中下游地區造成過大的沖擊。在1998年長江特大洪水災害中，中下游地區遭受了嚴重的洪澇災害，大量農田被淹沒，房屋倒塌，人員傷亡和財產損失巨大。若當時三峽水庫已經建成并投入使用，憑借其221.5億立方米的防洪庫容，可有效攔蓄洪水，削減洪峰流量，大大減輕中下游地區的防洪壓力。根據模擬分析，在類似1998年洪水的情況下，三峽水庫可將洪峰流量削減約20000-30000立方米/秒，使中下游地區的水位明顯降低，減少淹沒面積，保護大量的農田、房屋和基礎設施，拯救眾多生命，避免巨大的財產損失。在發電方面，三峽水電站共安裝34臺水輪發電機組，總裝機容量達2250萬千瓦，多年平均發電量約1000億千瓦時。所發電力通過輸電線路送往華東、華中等地區，為這些地區的經濟發展提供了穩定的電力支持，促進了當地的工業生產和居民生活用電需求。在航運上，三峽水庫改善了長江上游的通航條件，使得萬噸級船隊可以從上海直達重慶，大大提高了長江黃金水道的運輸能力，降低了運輸成本，促進了區域間的物資流通和經濟交流。在水資源利用方面，三峽水庫可以在枯水期為中下游地區補水，保障中下游地區的生產生活用水需求，同時還可以調節水庫周邊地區的氣候，改善生態環境。3.2水位調度方案三峽水庫的水位調度嚴格遵循“蓄清排渾，調枯補豐”的原則，呈現出明顯的季節性和周期性變化特點，這種變化對馬家溝滑坡穩定性產生了深遠影響。每年5月末至6月初，水庫水位降至汛期限制水位145m，目的是預留足夠的防洪庫容，以應對即將到來的汛期洪水。整個汛期（6-9月），除入庫流量大于下游河道安全泄量時攔截超額洪水使水庫水位抬高外，一般維持在145m運行。汛末10月，水庫開始蓄水，逐漸升高到175m運行，利用天然來水進行“蓄豐補枯”的調度，以滿足發電、航運以及下游用水等多方面的需求。12月至次年4月底，水庫按保證出力要求運行，并逐步降落，以增加下游流量和電站出力，但枯季消落最低水位不低于155m，以保證水庫回水變動區航道水深。在不同時期，三峽水庫采取了相應的水位調度策略。在汛期，當上游來水流量不超過30000立方米/秒時，庫區水位在規定的汛期水位變動范圍內，原則上由三峽集團負責調度。當上游來水流量超過3萬立方米/秒，但枝城流量小于56700立方米/秒時，一般由長江水利委員會負責調度。一旦枝城流量超過安全泄洪流量（56700立方米/秒），必須由長江水利委員會提出調度方案，報水利部批準。這種嚴格的調度機制確保了在洪水來臨時，三峽水庫能夠根據實際情況合理調整水位，有效削減洪峰，保護中下游地區的安全。在汛末蓄水期，為充分利用汛末洪水資源，在保障防洪安全的前提下，蓄水時間可提前至9月15日。一般情況下，三峽水庫平時在9月底會把水位調到162.0m，但如果來水多，就可以調高到165.0m，到10月底就可以蓄到175.0m正常蓄水位。在枯水期，經三峽水庫調節，宜昌市以下流量增加1000-2000立方米/秒，平均增加干流水深0.8米，不僅改善了中下游航道條件，還有效緩解了中下游工農業生產和沿江城鎮用水緊張的局面，對抑制長江入海口的咸潮入侵也發揮了重要作用。在滿足相關條件下，三峽水庫還可開展減輕中游防汛壓力的中小洪水調度，最高可攔蓄至150.0m運行。三峽水庫的水位調度還會根據地質災害防治對庫水位下降速率的要求，調整消落時間和水位范圍。規定集中消落開始時間提前到5月25日，當上下游雨水情滿足一定條件時，可以消落到144.9-146.5m的范圍內。3.3水位下降速率分析為深入探究三峽水庫水位下降速率的變化規律及其影響因素，本研究收集了2003-2022年三峽水庫水位數據，并運用線性回歸、周期分析等統計方法進行分析。結果表明，三峽水庫水位下降主要集中在每年的5-6月，水位下降速率呈現出明顯的年際變化和階段性特征。在年際變化方面，2003-2022年期間，三峽水庫水位下降速率最大值出現在2011年，達到了0.85m/d，最小值出現在2007年，為0.23m/d，二者相差較大。這表明不同年份間，水位下降速率存在顯著差異。這種差異主要與當年的來水情況密切相關。來水偏豐的年份，如2007年，水庫蓄水量充足，在水位下降過程中，為保證下游用水、航運等需求，出庫流量相對穩定，水位下降速率較為平緩；而來水偏枯的年份，如2011年，水庫蓄水量相對較少，在滿足下游各項需求的同時，為盡快降低水位至汛限水位，出庫流量較大，導致水位下降速率加快。上游地區的降水情況是影響來水的關鍵因素之一。若上游地區降水充沛，河流徑流量增大，水庫入庫水量增加，在水位下降階段，可調節的水量增多，水位下降速率就相對較慢；反之，若上游地區降水稀少，河流徑流量減小，水庫入庫水量不足，為達到汛限水位要求，就需要加大出庫流量，從而使水位下降速率加快。從階段性特征來看，2003-2010年期間，水位下降速率整體相對較低，平均值為0.42m/d。這一階段，三峽水庫處于運行初期，各項調度方案尚在逐步完善和優化過程中，對水位下降速率的控制較為謹慎，以確保水庫運行的安全性和穩定性，同時也為了更好地適應水庫運行對周邊環境和生態系統的影響。2011-2022年期間，水位下降速率有所增大，平均值達到0.53m/d。隨著三峽水庫運行經驗的積累和調度技術的成熟，在保證防洪、發電、航運等綜合效益的前提下，對水位下降速率的調控更加靈活，以更好地滿足不同時期的需求。隨著經濟社會的發展，下游地區對水資源的需求不斷增加，在水位下降階段，需要適當加大出庫流量，以保障下游地區的生產生活用水和航運需求，這也在一定程度上導致了水位下降速率的增大。三峽水庫水位下降速率還受到水庫運行調度方案的直接影響。在汛前消落期，為預留防洪庫容，水庫需在規定時間內將水位降至汛限水位145m。當上游來水情況穩定時，若按照常規調度方案，均勻控制出庫流量，水位下降速率相對穩定；但當遇到特殊情況，如上游來水突然減少或下游用水需求突然增加時，水庫管理部門會根據實際情況及時調整出庫流量，從而導致水位下降速率發生變化。在一些年份，為滿足下游灌溉用水高峰期的需求，水庫會加大出庫流量，使得水位下降速率在短期內明顯加快；而在另一些年份，若考慮到庫區地質災害防治等因素，對水位下降速率進行嚴格控制，出庫流量會相對減小，水位下降速率也會相應降低。四、馬家溝滑坡工程地質條件4.1滑坡地理位置及地貌特征馬家溝滑坡位于重慶市奉節縣境內，地處三峽庫區腹心地帶，地理位置為東經109°25′-109°30′，北緯31°05′-31°10′。該區域屬于典型的低山峽谷地貌，地形起伏較大，山巒連綿，溝壑縱橫。滑坡體發育于長江支流某河谷斜坡地段，周邊山脈海拔多在500-1000米之間，相對高差可達300-500米。從宏觀地貌上看，馬家溝滑坡處于一個向河谷傾斜的大型斜坡之上，斜坡總體坡度約為15°-30°。滑坡體呈不規則的長條狀，南北向展布，長約1500米，寬約800-1200米，面積約1.5平方公里。滑坡體后緣與斜坡后緣山體相連，地形相對較高，海拔約725米，后緣邊界清晰，呈弧形，多由基巖出露形成陡坎。滑坡體前緣直抵河谷岸邊，前緣受河水沖刷侵蝕作用明顯，形成了較為陡峭的臨空面，坡度可達40°-50°，且前緣部分已被河水淹沒或處于洪水期淹沒范圍。在滑坡體中部，地形相對較為平緩，坡度約為10°-15°，發育有一些小型的沖溝和平臺。沖溝多呈“V”字形，深度一般在5-15米之間，寬度在10-30米不等，溝內堆積有大量的松散土石。這些沖溝在雨季時，水流匯聚，加速了滑坡體表面的侵蝕和物質搬運，對滑坡的穩定性產生不利影響。平臺則相對較為平坦，多為第四系松散堆積物覆蓋，是當地居民進行農業生產和居住的主要區域。平臺上分布有農田、房屋等，人類活動較為頻繁，如開墾農田、修建房屋等活動改變了滑坡體表面的地形地貌和巖土體結構，增加了滑坡發生的潛在風險。滑坡體兩側邊界也較為清晰，多由自然沖溝或地形突變形成。東側邊界為一條較大的沖溝，溝谷深切，深度可達30-50米，寬度在50-100米左右，將滑坡體與東側山體分隔開來。西側邊界則是由于地形坡度的突然變化，從滑坡體的緩坡過渡到西側較陡的山坡，形成了明顯的邊界。這些邊界條件對滑坡的變形破壞具有一定的控制作用，滑坡體在變形過程中，往往沿著邊界發生滑動和錯動。4.2地質構造與地層巖性馬家溝滑坡所在區域的地質構造背景較為復雜，處于上揚子臺坪之秭歸臺褶束東緣，東部跨入北北東向的黃陵背斜，位于新華夏構造體系鄂西隆起帶北端和淮陽山字型構造體系復合部位。該區域經歷了多期次的構造運動，褶皺和斷裂構造發育，這些構造活動對區域內的地層巖性分布和巖體結構產生了顯著影響。區域內的褶皺構造主要為近南北向的秭歸向斜，軸向大致為南北方向，其核部由侏羅系中統下沙溪廟組（J2x）地層組成，兩翼依次出露侏羅系下統珍珠沖組（J1z）、三疊系上統須家河組（T3x）等地層。褶皺構造使得地層發生彎曲變形，形成了一系列的背斜和向斜構造，這種復雜的褶皺形態改變了地層的原始產狀，使得巖土體的受力狀態變得復雜，增加了滑坡發生的可能性。在向斜構造的核部，由于巖層受到擠壓，巖石破碎，節理裂隙發育，巖體完整性遭到破壞，抗剪強度降低，容易在外部因素的作用下發生滑動。斷裂構造主要有北北東向壓扭性水田壩斷裂，該斷裂規模較大，延伸較長，對區域內地層巖性的分布和地質災害的發育具有重要控制作用。斷裂帶附近巖石破碎，節理裂隙密集，地下水活動強烈，巖體的力學性質明顯降低。水田壩斷裂的活動導致了其兩側巖體的錯動和變形，使得巖體的完整性遭到破壞，為滑坡的形成提供了有利的地質條件。在馬家溝滑坡的形成過程中，斷裂構造可能起到了控制滑坡邊界和滑動方向的作用。滑坡體可能沿著斷裂帶的薄弱部位發生滑動，斷裂帶的存在使得滑坡體與周邊巖體之間的聯系減弱，增加了滑坡的不穩定性。馬家溝滑坡體主要由第四系松散堆積物和下伏基巖組成。第四系松散堆積物主要包括粉質黏土、碎石土等，厚度變化較大，一般在5-20米之間。粉質黏土呈黃褐色，可塑-硬塑狀態，具有一定的黏性和可塑性，其顆粒較細，透水性較差。在水位下降過程中，粉質黏土中的水分不易排出，容易形成較高的孔隙水壓力，降低土體的抗剪強度，對滑坡穩定性產生不利影響。碎石土主要由大小不等的碎石和砂土組成，碎石含量一般在30%-50%之間，粒徑較大，透水性較好。但由于碎石土顆粒間的黏結力較弱，在水流的沖刷和滲透作用下，容易發生顆粒的移動和流失，導致土體結構松散，也會影響滑坡的穩定性。下伏基巖主要為侏羅系上統沙溪廟組（J3s）砂巖和泥巖互層，砂巖呈灰白色、淺灰色，中-細粒結構，主要礦物成分有石英、長石等，巖石較堅硬，抗壓強度較高，一般在30-50MPa之間。泥巖呈紫紅色、灰綠色，泥質結構，主要礦物成分有黏土礦物等，巖石較軟，抗壓強度較低，一般在10-20MPa之間。砂巖和泥巖互層的地層結構，使得巖體的力學性質具有明顯的各向異性。泥巖作為相對軟弱的夾層，在地下水的長期作用下，容易發生軟化和泥化，強度降低，形成潛在的滑動面。當受到外部荷載或庫水位變化等因素的影響時，滑坡體可能沿著泥巖夾層發生滑動。這種軟硬相間的地層結構，增加了馬家溝滑坡的復雜性和不穩定性。4.3水文地質條件馬家溝滑坡的水文地質條件較為復雜，主要受大氣降水、三峽水庫水位變化以及地下水活動的影響。滑坡體的地下水類型主要包括松散巖類孔隙水和基巖裂隙水。松散巖類孔隙水主要賦存于滑坡體上部的第四系粉質黏土和碎石土中。粉質黏土孔隙細小，透水性較差，含水量相對較低，其孔隙水主要以結合水和毛管水的形式存在。碎石土孔隙較大，透水性較好，孔隙水含量相對較高，且具有一定的連通性，能夠形成較為明顯的滲流通道。在雨季，大氣降水大量入滲，使得松散巖類孔隙水水位迅速上升；而在旱季，隨著水分的蒸發和下滲，孔隙水水位逐漸下降。三峽水庫水位的周期性變化也會對松散巖類孔隙水產生影響。當水庫水位上升時，滑坡體前緣被淹沒，地下水受到庫水的頂托作用，水位相應抬高；當水庫水位下降時，滑坡體前緣的地下水與庫水之間形成水頭差，孔隙水會向庫水方向滲流，導致孔隙水水位下降。基巖裂隙水則主要賦存于下伏基巖的節理裂隙中。由于砂巖和泥巖互層的地層結構，節理裂隙發育程度在不同巖性層中存在差異。砂巖中節理裂隙相對較為發育，連通性較好，有利于基巖裂隙水的賦存和運移；而泥巖中節理裂隙相對較少，且多被黏土礦物充填，透水性較差，基巖裂隙水含量相對較低。基巖裂隙水的補給來源主要為大氣降水和松散巖類孔隙水的下滲。在地形和構造的控制下，基巖裂隙水一般沿節理裂隙向地勢較低的方向徑流，最終排泄于河谷或與其他含水層發生水力聯系。在馬家溝滑坡區域，基巖裂隙水的徑流方向主要受地形和斷裂構造的影響。在靠近斷裂帶的部位，由于巖石破碎，節理裂隙密集，基巖裂隙水的徑流速度相對較快，水力聯系也更為復雜。馬家溝滑坡體內的地下水位變化與三峽水庫水位密切相關。在三峽水庫蓄水前，滑坡體地下水位主要受大氣降水和地形的控制，水位變化相對較為平緩。蓄水后，三峽水庫水位的周期性漲落成為影響滑坡體地下水位的主要因素。當水庫水位上升時，滑坡體地下水位隨之上升，且上升速度較快，一般在水庫水位上升后的1-2天內，滑坡體地下水位就能達到相應的高度。這是因為水庫水位上升后，庫水通過滑坡體前緣的孔隙和裂隙迅速滲入，使得地下水位快速抬升。而當水庫水位下降時，滑坡體地下水位的下降則相對滯后。這是由于滑坡體巖土體具有一定的持水性，孔隙水在重力和滲流作用下逐漸排出，導致地下水位下降速度較慢。根據現場監測數據，水庫水位下降后，滑坡體地下水位一般需要3-5天才能明顯下降，且下降幅度小于水庫水位下降幅度。這種地下水位下降的滯后性，使得滑坡體內部在水庫水位下降過程中產生較大的孔隙水壓力差，對滑坡穩定性產生不利影響。五、水位下降速率對馬家溝滑坡穩定性影響的物理過程5.1水位下降對滑坡體滲流場的影響5.1.1滲流場變化規律為深入了解三峽水庫水位下降時馬家溝滑坡體內滲流場的變化規律，本研究運用數值模擬軟件SEEP/W建立了馬家溝滑坡的二維飽和-非飽和滲流模型。模型的建立基于對馬家溝滑坡工程地質條件的詳細勘查和分析，包括滑坡體的地形地貌、地層巖性、水文地質條件等信息。在模型中，充分考慮了巖土體的滲透系數、孔隙率、飽和度等參數，以及三峽水庫水位下降速率、初始水位等邊界條件。通過對不同水位下降速率工況下的滲流場進行模擬，得到了滑坡體內滲流場的變化特征。當三峽水庫水位下降時，滑坡體中的水在重力和水力梯度的作用下，向水庫方向滲流。在水位下降初期，由于滑坡體與水庫之間的水頭差較大，滲流速度較快。隨著水位下降的持續進行，滑坡體中的水分逐漸排出，水頭差減小，滲流速度逐漸降低。通過模擬結果可以看出，在水位下降速率為0.5m/d的工況下，水位下降初期，滑坡體前緣的滲流速度可達0.05m/d左右，而在水位下降后期，滲流速度降至0.01m/d以下。滲流方向也會隨著水位下降發生明顯變化。在水位下降前，滑坡體內的滲流方向主要受地形和地下水補給的影響，呈現出較為復雜的分布。而在水位下降過程中，滲流方向逐漸向水庫方向集中。以滑坡體中部某監測點為例，在水位下降前，滲流方向與水平方向夾角約為30°，隨著水位下降，滲流方向逐漸調整，在水位下降后期，與水平方向夾角接近0°，幾乎完全指向水庫方向。為驗證數值模擬結果的準確性，在馬家溝滑坡體上布置了多個監測孔，對地下水位和滲流速度進行實時監測。監測結果表明，數值模擬得到的滲流場變化規律與實際監測結果基本一致。在某一監測孔，模擬得到的水位下降過程中地下水位變化曲線與實際監測曲線的相關系數達到0.92，滲流速度的模擬值與實測值相對誤差在10%以內。這充分說明所建立的滲流模型能夠較為準確地反映三峽水庫水位下降時馬家溝滑坡體內滲流場的實際變化情況。5.1.2滲流壓力變化特征滲流壓力是水位下降過程中影響馬家溝滑坡穩定性的關鍵因素之一，其變化特征與滲流場的變化密切相關。當三峽水庫水位下降時，滑坡體中的孔隙水壓力隨之降低，產生滲流壓力。滲流壓力的大小和分布對滑坡體的有效應力和穩定性有著重要影響。在水位下降過程中，滲流壓力在滑坡體內的分布呈現出明顯的不均勻性。滑坡體前緣和下部由于與水庫直接相連，水位下降時水頭差變化較大，滲流壓力相對較高。而滑坡體后緣和上部，由于距離水庫較遠，水位下降對其影響相對較小，滲流壓力較低。通過數值模擬分析，在水位下降速率為1.0m/d的工況下，滑坡體前緣底部的滲流壓力可達50kPa左右，而滑坡體后緣頂部的滲流壓力僅為5kPa左右。滲流壓力還會隨著時間的推移發生動態變化。在水位下降初期，滲流壓力迅速增大，這是由于水位下降速度較快，孔隙水壓力來不及消散，導致滲流壓力急劇上升。隨著時間的增加，孔隙水壓力逐漸消散，滲流壓力開始逐漸減小。在水位下降后期，滲流壓力趨于穩定，但仍維持在一定水平。以某一典型監測點為例，在水位下降開始后的第1天，滲流壓力從初始的10kPa迅速上升至30kPa，隨后在第3-5天，滲流壓力逐漸減小至20kPa左右，并在后續時間內保持相對穩定。滲流壓力對滑坡穩定性的影響主要體現在對有效應力的改變上。根據有效應力原理，有效應力等于總應力減去孔隙水壓力。當滲流壓力增大時，孔隙水壓力減小，有效應力增大，這會使滑坡體的抗剪強度增加，有利于滑坡的穩定。但在某些情況下，滲流壓力的變化也可能導致滑坡體局部應力集中，當應力超過巖土體的強度極限時，就會引發滑坡的變形和破壞。在滑坡體前緣的局部區域，由于滲流壓力的作用，可能會導致土體出現剪切破壞，進而引發滑坡的滑動。5.2水位下降對滑坡體力學性質的影響5.2.1巖土體強度參數變化三峽水庫水位下降會引發馬家溝滑坡體巖土體強度參數的顯著變化，其中抗剪強度的改變尤為關鍵。抗剪強度是衡量巖土體抵抗剪切破壞能力的重要指標，它直接關系到滑坡體的穩定性。水位下降時，滑坡體中的孔隙水壓力發生變化，這是導致抗剪強度改變的重要原因之一。孔隙水壓力是指土體孔隙中所存在的水壓力，它對巖土體的有效應力和抗剪強度有著重要影響。根據有效應力原理，有效應力等于總應力減去孔隙水壓力。當水位下降時，孔隙水壓力減小，有效應力增大。對于飽和土體，抗剪強度與有效應力密切相關，有效應力的增大使得土體顆粒間的摩擦力增大，從而提高了土體的抗剪強度。為了定量分析水位下降對抗剪強度的影響，本研究采用了室內土工試驗和理論分析相結合的方法。室內土工試驗選取了馬家溝滑坡體不同部位的巖土體樣本，進行了直剪試驗和三軸剪切試驗。在直剪試驗中，對飽和狀態和水位下降后的非飽和狀態下的樣本分別進行測試。結果表明，在水位下降后，土體的黏聚力和內摩擦角均有所變化。對于粉質黏土樣本，黏聚力從飽和狀態下的15kPa增加到非飽和狀態下的18kPa，內摩擦角從20°增大到22°。這是因為水位下降導致孔隙水壓力減小，土體顆粒間的有效應力增大，顆粒間的相互作用力增強，使得土體的黏聚力和內摩擦角增大。對于碎石土樣本，由于其顆粒間的黏結力相對較弱，水位下降后，顆粒間的摩擦力變化相對較小，但由于孔隙水壓力的減小，有效應力增大，也在一定程度上提高了其抗剪強度。理論分析方面，運用基于有效應力原理的抗剪強度理論，結合滑坡體的實際情況進行計算。根據摩爾-庫侖強度理論，土體的抗剪強度表達式為\tau=c+\sigma\tan\varphi，其中\tau為抗剪強度，c為黏聚力，\sigma為有效應力，\varphi為內摩擦角。在水位下降過程中，通過計算孔隙水壓力的變化，進而得到有效應力的改變，再代入抗剪強度公式中，分析抗剪強度的變化情況。假設在某一水位下降工況下，通過滲流分析得到某一土體單元的孔隙水壓力從初始的30kPa下降到10kPa，總應力保持不變為50kPa。根據有效應力原理，初始有效應力為50-30=20kPa，水位下降后有效應力變為50-10=40kPa。若該土體的黏聚力c=10kPa，內摩擦角\varphi=25?°，則初始抗剪強度\tau_1=10+20\tan25?°\approx19.3kPa，水位下降后的抗剪強度\tau_2=10+40\tan25?°\approx28.6kPa。通過理論計算進一步驗證了水位下降導致孔隙水壓力減小，有效應力增大，從而使抗剪強度提高的結論。然而，在實際情況中，水位下降對巖土體強度參數的影響是復雜的。除了孔隙水壓力的變化外，還可能受到其他因素的影響。水位下降可能導致巖土體的干濕循環，使土體發生收縮和膨脹，從而破壞土體的結構，降低其抗剪強度。長期的干濕循環會使土體顆粒間的連接逐漸弱化，導致黏聚力降低。地下水化學成分的變化也可能對巖土體的性質產生影響，進而改變其強度參數。若地下水中含有侵蝕性物質，在水位下降過程中，這些物質可能會與巖土體發生化學反應，導致巖土體的成分和結構改變，影響其抗剪強度。5.2.2滑坡體應力應變分布通過數值模擬深入研究三峽水庫水位下降過程中馬家溝滑坡體的應力應變分布特征，對于揭示滑坡變形破壞機制、評估滑坡穩定性具有重要意義。本研究采用有限元軟件FLAC，建立了馬家溝滑坡的三維數值模型。模型充分考慮了滑坡體的工程地質條件，包括地形地貌、地層巖性、地質構造等。在模型中，將滑坡體劃分為多個單元，賦予每個單元相應的巖土體力學參數，如彈性模量、泊松比、密度等。考慮了三峽水庫水位下降速率、初始水位等邊界條件，模擬了不同水位下降速率工況下滑坡體的應力應變變化情況。在水位下降過程中，滑坡體的應力分布呈現出明顯的不均勻性。滑坡體前緣和后緣是應力變化較為顯著的區域。滑坡體前緣由于直接受到庫水的作用，在水位下降時，水壓力的減小導致前緣土體的應力狀態發生改變。在水位下降初期，前緣土體的水平向應力迅速減小，而垂直向應力則有所增加。這是因為水壓力的減小使得土體在水平方向上的約束減弱，而重力作用下的垂直向應力相對增大。隨著水位下降的持續進行，前緣土體的剪應力逐漸增大。當剪應力超過土體的抗剪強度時，就會導致土體發生剪切破壞，這是滑坡前緣容易出現坍塌、滑移等現象的重要原因。在水位下降速率為1.0m/d的工況下，通過數值模擬計算得到，滑坡體前緣底部在水位下降5天后，剪應力達到50kPa，而該部位土體的抗剪強度為45kPa，此時土體發生剪切破壞的可能性增大。滑坡體后緣則由于水位下降導致地下水位降低，土體的自重應力增加，從而使得后緣土體的垂直向應力增大。后緣土體在自重作用下，可能會產生拉應力。當拉應力超過土體的抗拉強度時，就會導致土體出現拉裂破壞。在實際觀測中，馬家溝滑坡體后緣常常出現一些拉張裂縫，這與數值模擬結果相吻合。在水位下降速率為0.5m/d的工況下，模擬結果顯示，滑坡體后緣頂部在水位下降10天后，拉應力達到10kPa，而該部位土體的抗拉強度為8kPa，此時后緣土體可能會出現拉裂破壞。滑坡體的應變分布也與應力分布密切相關。在應力集中的區域，應變也相對較大。滑坡體前緣和后緣是應變集中的主要區域。在前緣，由于剪應力的作用，土體主要發生剪切應變。剪切應變的積累會導致土體的變形逐漸增大，最終可能引發滑坡的滑動。在后緣，由于拉應力的作用，土體主要發生拉伸應變。拉伸應變的發展會使拉張裂縫不斷擴展，進一步削弱滑坡體的穩定性。通過數值模擬得到的應變云圖可以清晰地看到，在水位下降過程中，滑坡體前緣和后緣的應變值明顯高于其他部位。在水位下降速率為1.5m/d的工況下，滑坡體前緣的最大剪切應變達到0.05，后緣的最大拉伸應變達到0.03，這些較大的應變值表明滑坡體在這些部位更容易發生變形破壞。六、基于數值模擬的水位下降速率對馬家溝滑坡穩定性影響研究6.1數值模型建立6.1.1模型概化為準確模擬三峽水庫水位下降速率對馬家溝滑坡穩定性的影響，對馬家溝滑坡進行了合理的簡化和抽象，構建了適用于數值模擬的模型。在模型概化過程中，充分考慮了滑坡的地形地貌、地層巖性、地質構造以及水文地質條件等關鍵因素。基于詳細的地質勘查資料，對滑坡體的地形進行了精確的數字化處理。通過三維激光掃描和地形測繪等技術手段，獲取了滑坡體的高精度地形數據，并利用專業的地理信息系統（GIS）軟件，將地形數據轉化為數值模型中的地形表面。在模型中，準確再現了滑坡體的起伏形態、坡度變化以及與周邊地形的銜接關系，為后續的模擬分析提供了可靠的地形基礎。在考慮地層巖性時，根據地質鉆探和巖土體物理力學性質試驗結果，將滑坡體劃分為不同的地層單元。對于第四系松散堆積物，進一步細分為粉質黏土和碎石土單元；對于下伏基巖，按照侏羅系上統沙溪廟組（J3s）砂巖和泥巖互層的實際情況進行了分層模擬。為每個地層單元賦予了相應的物理力學參數，包括密度、彈性模量、泊松比、黏聚力和內摩擦角等，這些參數的取值基于室內土工試驗和原位測試結果，并結合工程經驗進行了合理的調整。地質構造對滑坡穩定性的影響也在模型中得到了充分體現。根據區域地質資料，明確了馬家溝滑坡所在區域的褶皺和斷裂構造分布情況。在模型中，通過設置結構面來模擬斷裂構造的力學特性和對滑坡體變形的影響。對于褶皺構造，考慮了地層的彎曲形態和產狀變化對巖土體應力分布的影響，采用相應的數學模型進行模擬。在處理水文地質條件時，考慮了滑坡體的地下水類型和水位變化情況。將地下水分為松散巖類孔隙水和基巖裂隙水，并分別設置了相應的滲流參數，如滲透系數、孔隙率等。結合三峽水庫水位的歷史數據和運行調度方案，設定了模型的邊界條件，模擬了水庫水位下降過程中滑坡體地下水的滲流變化。6.1.2參數選取數值模型中各項參數的準確選取是保證模擬結果可靠性的關鍵。在參數選取過程中，綜合運用了室內試驗、原位測試以及經驗取值等多種方法。室內土工試驗是獲取巖土體物理力學參數的重要手段。針對馬家溝滑坡體不同地層的巖土體，采集了大量的原狀土樣和巖石試件，進行了一系列的室內試驗。對于粉質黏土和碎石土，進行了直剪試驗，以測定其黏聚力和內摩擦角。在直剪試驗中，按照相關標準規范，對土樣進行了不同法向應力下的剪切測試，得到了不同工況下土樣的抗剪強度參數。通過三軸壓縮試驗，測定了巖土體的彈性模量、泊松比等參數，以反映其在復雜應力狀態下的力學特性。在三軸壓縮試驗中，對土樣施加不同的圍壓和軸向壓力，記錄土樣的應力-應變曲線，從而計算出彈性模量和泊松比。通過顆粒分析試驗，確定了土樣的顆粒組成和級配情況，這對于分析土樣的滲透特性和力學性質具有重要意義。原位測試則能夠更真實地反映巖土體在天然狀態下的力學性質和滲透特性。采用標準貫入試驗，測定了滑坡體中不同深度處的巖土體密實度和強度指標。在標準貫入試驗中，將標準貫入器打入巖土體中，記錄貫入一定深度所需的錘擊數，根據錘擊數與巖土體強度的經驗關系，估算巖土體的強度參數。利用平板載荷試驗，測定了巖土體的承載力和變形模量，為模型提供了更準確的地基參數。在平板載荷試驗中，在巖土體表面放置承載板，逐級施加荷載，記錄承載板的沉降量，根據荷載-沉降曲線計算出巖土體的承載力和變形模量。通過鉆孔壓水試驗，獲取了巖土體的滲透系數，以反映其透水性。在鉆孔壓水試驗中，向鉆孔內壓水，測量一定時間內的壓入水量和壓力，根據達西定律計算出滲透系數。對于一些難以通過試驗直接測定的參數，參考相關工程經驗和類似地質條件下的研究成果進行取值。在確定滑坡體的重度時，參考了該地區同類巖土體的重度取值范圍，并結合滑坡體的實際物質組成進行了適當調整。在確定巖土體的剪脹角時，由于其試驗測定較為困難，根據經驗公式和已有研究成果，結合滑坡體的巖性和結構特點進行了合理取值。6.2模擬工況設置為全面深入探究三峽水庫水位下降速率對馬家溝滑坡穩定性的影響，設置了多種不同的水位下降速率工況。參考三峽水庫歷年水位下降速率數據，結合實際運行情況和研究需求，選取了0.2m/d、0.4m/d、0.6m/d、0.8m/d、1.0m/d這5種具有代表性的水位下降速率。這些速率涵蓋了三峽水庫水位下降速率的常見范圍，能夠較為全面地反映不同下降速率對滑坡穩定性的影響。在實際情況中，三峽水庫水位下降速率在不同年份和季節會有所波動，通過設置這5種工況，可以對不同程度的水位下降情況進行模擬分析，從而更準確地掌握水位下降速率與滑坡穩定性之間的關系。在模擬過程中，明確了初始條件和邊界條件。初始條件方面，設定三峽水庫初始水位為175m，這是三峽水庫的正常蓄水位，也是水位下降的起始狀態。此時，馬家溝滑坡體處于相對穩定的狀態，滑坡體內的滲流場和應力應變場也處于初始平衡狀態。在水位下降前，滑坡體內部的孔隙水壓力分布較為均勻，地下水位與水庫水位基本持平，滑坡體的應力主要由自重和上覆巖土體的壓力產生。邊界條件設置如下：在滑坡體的底部，設置為不透水邊界，以模擬實際地質條件下基巖對地下水的阻隔作用。由于基巖的透水性較差，地下水在基巖中難以流動，因此將滑坡體底部設置為不透水邊界能夠更真實地反映實際情況。在滑坡體的側面，根據實際地形和水流情況，設置為流量邊界或水頭邊界。對于靠近河流或水庫的側面，設置為水頭邊界，其水頭值根據三峽水庫水位的變化而實時調整，以反映水庫水位對滑坡體側面地下水的影響。對于遠離河流或水庫的側面，設置為流量邊界，根據實際觀測或經驗公式確定其流量值，以模擬地下水在滑坡體側面的流入或流出。在滑坡體的表面，設置為自由排水邊界，允許水分自由蒸發和下滲。考慮到大氣降水的影響，在表面邊界條件中，根據當地的氣象資料，設置了不同的降水強度和歷時，以模擬不同降雨條件下對滑坡穩定性的影響。在雨季，降水強度較大，歷時較長，會增加滑坡體的含水量，降低巖土體的抗剪強度，從而對滑坡穩定性產生不利影響。而在旱季，降水強度較小，歷時較短，對滑坡穩定性的影響相對較小。6.3模擬結果分析6.3.1不同水位下降速率下的滑坡穩定性系數變化通過對不同水位下降速率工況的數值模擬，得到了馬家溝滑坡穩定性系數隨時間的變化曲線，清晰地展示了水位下降速率與滑坡穩定性之間的定量關系。當水位下降速率為0.2m/d時，在水位下降初期，滑坡穩定性系數略有下降，從初始的1.25降至1.23，下降幅度較小。這是因為水位下降速率較慢，滑坡體中的孔隙水有足夠的時間排出，孔隙水壓力能夠及時消散，有效應力變化相對較小，對滑坡穩定性的影響不大。隨著水位下降的持續進行，穩定性系數基本保持穩定，維持在1.23左右。這表明在這種緩慢的水位下降速率下，滑坡體能夠較好地適應水位變化，保持相對穩定的狀態。當水位下降速率增大到0.4m/d時，水位下降初期，穩定性系數下降較為明顯，從1.25降至1.18。這是由于水位下降速率加快，孔隙水壓力消散相對滯后，導致有效應力減小，抗滑力降低，從而使滑坡穩定性下降。在水位下降后期，穩定性系數逐漸趨于穩定，但仍低于初始值，維持在1.18左右。這說明在該水位下降速率下，雖然滑坡體在后期逐漸適應了水位變化，但穩定性已受到一定程度的削弱。在水位下降速率為0.6m/d的工況下，水位下降初期，穩定性系數迅速下降，從1.25降至1.10。這是因為水位下降速率進一步加快，孔隙水壓力消散困難，產生了較大的超孔隙水壓力，有效應力大幅減小，滑坡體的抗滑力顯著降低，穩定性受到嚴重影響。在水位下降后期，穩定性系數雖有所回升，但仍處于較低水平，僅回升至1.12左右。這表明在這種較快的水位下降速率下，滑坡體的穩定性受到了較大破壞，即使在后期也難以完全恢復到初始穩定狀態。當水位下降速率達到0.8m/d時，水位下降初期，穩定性系數急劇下降，從1.25降至1.02，已接近滑坡的極限平衡狀態（穩定性系數為1.0）。此時，孔隙水壓力消散嚴重滯后，超孔隙水壓力急劇增大，有效應力大幅減小，滑坡體的抗滑力急劇降低，滑坡處于高度不穩定狀態。在水位下降后期，穩定性系數略有回升，但仍低于1.05，滑坡仍處于不穩定狀態。這說明在如此快的水位下降速率下，滑坡體極易發生失穩滑動。當水位下降速率為1.0m/d時，水位下降初期，穩定性系數迅速降至0.95，滑坡已處于明顯的不穩定狀態。隨著水位下降的進行，穩定性系數繼續下降，最終降至0.90以下。這表明在這種極快的水位下降速率下，滑坡體內部的應力應變狀態急劇惡化，抗滑力遠小于下滑力，滑坡隨時可能發生大規模滑動。通過對不同水位下降速率工況下穩定性系數變化的對比分析，可以看出水位下降速率對馬家溝滑坡穩定性的影響非常顯著。隨著水位下降速率的增大，滑坡穩定性系數下降的幅度和速度都明顯增加，滑坡的穩定性逐漸降低，失穩的風險逐漸增大。6.3.2滑坡體位移和變形特征在不同水位下降速率工況下，馬家溝滑坡體的位移和變形呈現出明顯的差異，這些差異反映了水位下降速率對滑坡穩定性的重要影響。當水位下降速率為0.2m/d時，滑坡體的位移主要集中在前緣和后緣。前緣位移方向主要為向臨空面的水平位移，后緣則主要表現為垂直向下的位移。通過數值模擬得到，滑坡體前緣最大水平位移約為5cm，后緣最大垂直位移約為3cm。滑坡體的變形主要為局部的微小變形，整體變形較為均勻，未出現明顯的裂縫和錯動。這是因為水位下降速率較慢，滑坡體內部的應力調整較為緩慢，能夠較好地適應水位變化，從而保持相對穩定的狀態。當水位下降速率增大到0.4m/d時，滑坡體的位移明顯增大。前緣最大水平位移增加到10cm左右，后緣最大垂直位移增加到6cm左右。滑坡體的變形也有所加劇，在滑坡體中部出現了一些微小的拉張裂縫，長度一般在1-2m之間，寬度約為1-2cm。這表明在該水位下降速率下，滑坡體內部的應力變化較大，已超出了部分土體的抗拉強度，導致裂縫的產生，滑坡的穩定性受到一定程度的影響。在水位下降速率為0.6m/d的工況下，滑坡體的位移進一步增大。前緣最大水平位移達到15cm左右，后緣最大垂直位移達到10cm左右。滑坡體的變形更加明顯，中部的拉張裂縫進一步擴展，長度可達3-5m，寬度增加到3-5cm。在滑坡體前緣，還出現了一些剪切裂縫，長度在2-3m之間，寬度約為2-3cm。這說明在這種較快的水位下降速率下，滑坡體內部的應力集中現象更加嚴重，土體的抗剪強度和抗拉強度受到較大破壞，滑坡的穩定性明顯降低。當水位下降速率達到0.8m/d時，滑坡體的位移急劇增大。前緣最大水平位移達到25cm左右，后緣最大垂直位移達到15cm左右。滑坡體的變形十分顯著，中部的拉張裂縫和前緣的剪切裂縫相互貫通，形成了明顯的滑動面，滑動面深度可達5-8m。這表明在如此快的水位下降速率下，滑坡體已接近失穩狀態，隨時可能發生滑動。當水位下降速率為1.0m/d時，滑坡體的位移迅速增大，前緣最大水平位移超過30cm，后緣最大垂直位移超過20cm。滑坡體出現了大規模的變形破壞，滑動面進一步擴展，滑坡體整體呈現出明顯的滑動趨勢。這說明在這種極快的水位下降速率下，滑坡體已處于失穩狀態，即將發生大規模的滑動。隨著水位下降速率的增大，馬家溝滑坡體的位移和變形逐漸加劇，裂縫不斷擴展和貫通，滑動面逐漸形成，滑坡的穩定性逐漸降低，失穩的風險逐漸增大。這些位移和變形特征的變化與穩定性系數的變化趨勢一致，進一步驗證了水位下降速率對馬家溝滑坡穩定性的顯著影響。七、實例分析與驗證7.1馬家溝滑坡歷史變形與水位下降的相關性分析為深入探究三峽水庫水位下降與馬家溝滑坡變形之間的內在聯系，本研究廣泛收集了2003-2022年期間馬家溝滑坡的歷史變形數據。這些數據涵蓋了滑坡體的地表位移、裂縫發育以及深部位移等多個方面，來源包括專業監測單位的長期監測記錄、相關科研項目的實地測量數據以及政府部門的地質災害監測報告等。通過對這些豐富數據的系統分析，我們發現馬家溝滑坡的變形與三峽水庫水位下降呈現出顯著的相關性。在2010年三峽水庫水位下降期間，水位下降速率較快，達到了0.6m/d左右。從滑坡體的地表位移監測數據來看，在水位下降的一個月內，滑坡體前緣的水平位移明顯增大，累計位移量達到了15cm，而后緣的垂直位移也增加了8cm。通過對滑坡體裂縫的監測發現，裂縫寬度和長度均有顯著增長，新出現了多條長度超過5m的裂縫，且裂縫寬度達到了3-5cm。對滑坡體深部位移的監測結果顯示，在滑帶附近，位移量也有明顯增加，最大位移達到了10cm左右。而在2015年，三峽水庫水位下降速率相對較慢，為0.3m/d。在此期間，滑坡體的變形相對較小，前緣水平位移僅增加了5cm，后緣垂直位移增加了3cm，裂縫發育也相對不明顯，僅出現了少量長度在2m以內的細微裂縫，滑帶附近的深部位移增加量在5cm以內。為了更直觀地展示二者之間的關系，本研究繪制了水位下降速率與滑坡變形量的散點圖。從散點圖中可以清晰地看出，隨著水位下降速率的增大，滑坡體的位移量呈現出明顯的上升趨勢。通過相關性分析計算得到，水位下降速率與滑坡體前緣水平位移的相關系數達到了0.85，與后緣垂直位移的相關系數為0.82，與裂縫長度的相關系數為0.80，與滑帶深部位移的相關系數為0.83。這些高相關系數進一步證實了水位下降速率與馬家溝滑坡變形之間存在著密切的正相關關系。即水位下降速率越快，滑坡體的變形越顯著，滑坡的穩定性受到的影響越大。這一結論與前面通過數值模擬和理論分析得到的結果相互印證，為進一步研究三峽水庫水位下降速率對馬家溝滑坡穩定性的影響提供了有力的實際數據支持。7.2現場監測數據驗證為進一步驗證數值模擬結果的準確性，在馬家溝滑坡現場建立了全面的監測體系，對滑坡體的位移、地下水位和孔隙水壓力等關鍵參數進行長期實時監測。在滑坡體上，沿不同方向和高程共布置了20個位移監測點，采用全站儀和GPS相結合的監測方法，定期測量各監測點的三維坐標，以獲取滑坡體的水平和垂直位移數據。在滑坡體的前緣、中部和后緣分別設置了5個地下水位監測孔，利用自動水位計實時記錄地下水位的變化情況。在滑坡體的不同部位布置了10個孔隙水壓力監測點，采用孔隙水壓力計監測孔隙水壓力的大小和變化。這些監測點的布置充分考慮了滑坡體的地形地貌、地質構造以及可能的變形破壞模式，能夠全面反映滑坡體在三峽水庫水位下降過程中的變化情況。選取2018年三峽水庫水位下降期間的監測數據與數值模擬結果進行對比分析。在該時段內，三峽水庫水位從175m下降至145m，水位下降速率約為0.6m/d。從位移監測數據來看，數值模擬預測的滑坡體前緣最大水平位移為13cm，后緣最大垂直位移為8cm；而現場監測得到的前緣最大水平位移為12cm，后緣最大