【正文】 參考文獻 [1]唐紅梅，陳洪凱，關(guān)明芳《邊坡滲流研究現(xiàn)狀及趨勢分析》， 重慶交通學院學報， [2]沈 燁《邊坡穩(wěn)定分析方法綜述》，科技資訊， [3]廖珊珊，張玉成，胡海英《邊坡穩(wěn)定性影響因素的探討》，廣東水利水電， [4]李春忠《地下水滲流作用下基坑邊坡穩(wěn)定性研究》， [5]陳五一，羅玉龍《滲流對土質(zhì)邊坡穩(wěn)定性影響分析》，武漢大學學報， [6]鄔攀《滲流作用下邊坡穩(wěn)定性分析》， [7]謝羅峰《滲流作用下邊坡穩(wěn)定性研究》，南京水利科學研究院， [8]毛昶熙、陳平等《滲流作用下的壩坡穩(wěn)定有限單元分析》，巖土工程學報 ， 。 3) 進行巖質(zhì)岸坡裂隙滲流機理及滲流 應(yīng)力耦合作用研究 ， 量化裂隙水壓力的計算理論及變化規(guī)律 。 2) 庫位變動帶飽和滲流、非飽和滲流相互轉(zhuǎn)化過程研究 ， 明確岸坡巖土體在不同時間、不同空間部位的滲透系數(shù)取值規(guī)律 。 據(jù)此可見 ， 水庫每年水位變幅 30m 左右 。 而對于三峽水庫而言 ， 庫岸邊坡滲流問題顯得極為重要 。 主要側(cè)重于理論的引用 ， 鮮見與野外鉆孔、油水實驗資料結(jié)合對松散介質(zhì)滲透系數(shù)進行空間估值的實例 。 我國學者胡尊國 ( 1992) 研究指出多孔介質(zhì)的孔隙形狀具有分形的特征 ， 并認為分維數(shù)是刻畫或描述不同地質(zhì)成因多孔介質(zhì)的一個定量指標 ，但未進一步 說明如何利用分維值估計滲透系數(shù)等參數(shù) 。 Ababou( 1988) 根據(jù) Mount sinon 含水層中得到的數(shù)據(jù)分析 ， 指出滲透系數(shù)的對數(shù)表現(xiàn)出分形規(guī)律 。 間接解法 ： 給定參數(shù)試探值 ， 求解正問題 ， 得到計算水頭值 ， 再與已知觀測水頭值比較 ， 求得水頭誤差最小條件下的參數(shù) 。 雖不全等于實際巖體水力學參數(shù) ， 但反映了整個研究區(qū)域巖體水力學特征 ， 其關(guān)鍵是動態(tài)數(shù)據(jù)的可靠性和模型選擇的合理性 。 數(shù)學模型反演求解法 運用地下水動態(tài)信息 ， 用建立的數(shù)學模型來反演求解 。 找出了它的主變異方向和變異程度 ， 并對滲透系數(shù)進行了量優(yōu)估計 ， 進而給出了它的空間變異分布圖和參數(shù)分區(qū)圖 。 Demarsily、 Kitanidis、 Vomvoris、 Amleto 和 Murashige、雷志棟和楊詩秀把它引入水文及水文地質(zhì)領(lǐng)域 ， 研究了含水層滲透系數(shù)的空間變化特征 。 大多數(shù)研究者選擇隨機模型作為研究手段 ， 是假定含水介質(zhì)滲透系數(shù)滿足對數(shù)正態(tài)分布 ， 但缺乏完備的理論和充分的實踐資料支持 ， 因此其應(yīng)用受到很大的限制 。 因此 ， 該方法的最大特征是考慮了參數(shù)的結(jié)構(gòu)性和各向異性特征 。 這一方法的思路是 ： 將研究區(qū)劃分為方形網(wǎng)格 ，每一網(wǎng)格被認為一致均勻 ， 某一網(wǎng)格參數(shù)被認為與其周圍緊鄰的網(wǎng)格多數(shù)相關(guān) ，井建立參數(shù)自回歸方程 。 ? 連續(xù)相關(guān)參數(shù)法 美國學者 Celhar 最先考慮了參數(shù)之間的相關(guān)性 ， 用連續(xù)相關(guān)參數(shù)法研究了含水介質(zhì)參數(shù)和地下水流問題 ， 該方法的思路可歸納為 ： 滲透系數(shù) K 可表示為其均值 K和擾動值 K’之和 ， 用 Fourierstielt ijes 積分變換表示各類擾動值及滲透系數(shù)的自協(xié)方差函數(shù) ； 由 Fourierst ieltijes 逆變換求得滲透系數(shù)的譜密度函數(shù) ， 并由水流微分方程建立滲透系數(shù)和水位譜密度函數(shù)之間的關(guān)系 ， 通過研究滲透系數(shù)的譜密度函數(shù) ， 達到研究參數(shù)空間變異的目的 。 最早提出離散獨立參數(shù)法的學者是 Warren 和 Price( 1961)， 運用該方法研究了滲透系數(shù)的不同頻率分布對等效滲透系數(shù)的作用 。 其基本思想是將研究區(qū)劃分為不同的塊段 。Hoeksema 和 kitcmidis( 1985) 研究指出 ， 除極少數(shù)外 ， 導水系數(shù)、滲透系數(shù)和貯水系數(shù)等參數(shù)都可以認為是服從對數(shù)正態(tài)分布的 。 2） 隨機理論 Law( 1944) 研究油田勘探的巖芯資料時 ， 首先發(fā)現(xiàn)滲透系數(shù)服從對數(shù)正態(tài)分布 。Manath N. Panda 等在研究用粒度分布參數(shù)計算松散介質(zhì)的滲透率中指出 ， 對非均質(zhì)介質(zhì)而言 ， 介質(zhì)粒度均值的變化對滲透率的貢獻比孔隙度或分選程度要大 。 王釗等通過一些簡化假定 ， 推導出了適用于潛水完整單井的滲透系數(shù)表達式 ， 經(jīng)修正得出了非完整潛水單井天然來水過程測定滲透系數(shù)的方法 。 滲透系數(shù) 確定型模型結(jié)合實驗方法 包括室內(nèi)和野外實驗方法 ， 其中野外實驗方法包括野外幾何測量法和野外抽( 壓 ) 水試驗法 。凌道盛等提出了有自由面滲流分析的虛節(jié)點法 ， 其特點屬于一種真正的固定網(wǎng)格法 ， 對網(wǎng)格不作任何改動 ， 并能精確地描述跨越自由面單元的滲透矩陣 ， 可同時適用于平面滲流問題和空間問題 。 Desai 法的缺點是每一迭代步需要計算出自由面的具體位置、自由面的面積以及法向流量且工作量大 ， 而且這一方法的全部調(diào)整過程均依賴于第一次有限元分析結(jié)果 。 王賢能則根據(jù)高斯點的位置確定高斯點 對應(yīng)的滲透系數(shù) ， 計算單個高斯點對滲透矩陣的貢獻 ， 對所有高斯點求和得到單元的滲透矩陣 ； 第二類是調(diào)整流量法 ， 如 Desai 提出的剩余流量法 ， 張有天的初流量法等 。 Bathe 法在迭代過程中把跨越自由面和自由面以上單元的滲透系數(shù)一律按實際滲透系數(shù)的千分之一計算 ， 與實際不吻合 。 第二是固定網(wǎng)格法 : 自 1973 年 Ne uman 提出用不變網(wǎng)格分析有自由面滲流的 Galerkin 法以來 ， 出現(xiàn)了多種固定網(wǎng)格法 。 滲流自由面 第一是移動網(wǎng)格法 : 這類方法把自由面當作可移動邊界處理 ， 在迭代過程中不斷改變自由面的位置 ， 使網(wǎng)格發(fā)生變化 ， 直到自由面位置穩(wěn)定為止 。 王媛等采用四自由度全耦合法 ， 進行了基于增量理論的復雜裂隙巖體滲流與應(yīng)力彈塑性全耦合法分析 。 朱珍德利用流體擴散能量迭加原理推求滲透張量公式 ， 提出了滲透張量隨裂隙損傷加劇而改變以及裂隙巖體非穩(wěn)態(tài)滲流場與損傷場合模型 ,并以此編成有限元程序應(yīng)用到三峽永久船閘高邊坡穩(wěn)定性之中 ； 盛金昌建立能反映裂隙剪脹特性的非線性彈塑性本構(gòu)模型 ， 考慮裂隙非線性法向變形和剪脹對裂隙隙寬的影響 ， 比較實際地反映了裂隙變形對巖體滲透性的影響 。 建立邊坡巖體滲流場與應(yīng)力場耦合分析數(shù)學模型也是一個重要的問題 。 關(guān)于剪應(yīng)力和剪應(yīng)變對滲透場的影響 ， 目前研究成果較少 。 法國專家 Louis認為在低應(yīng)力狀態(tài)下滲透系數(shù)與有效應(yīng)力之間呈負指數(shù)關(guān)系即 K= K 0 e aw ( 9) 式中 ， K裂隙法線方向的有效應(yīng)力 ； K 0初始滲透系數(shù)， a常系數(shù) Carlsson 等基于室內(nèi)試驗認為滲透系數(shù)和法向應(yīng)力之間的存在高度是非線性 ， 并得到滲透系數(shù)與自重應(yīng)力之間的關(guān)系 ； 郭雪莽考慮到巖體滲流 應(yīng)力耦合作用是變形引起滲透系數(shù)的改變 ， 推導了滲透系數(shù)與應(yīng)變的本構(gòu)關(guān)系 。 應(yīng)力場和滲流場的耦合作用 巖體邊坡滲流分析中 ， 應(yīng)力場和滲流場的耦合問題十分重要 ， 著名的法國拱壩的失事為一典型實例 。其最大的特點是可以直接從實驗中直接觀察到滲流現(xiàn)象 。 介玉新等將適體坐標變換方法用于滲流計算 ， 即通過求Poisson 方程自動生成實際物理邊界的計算網(wǎng)絡(luò) ， 并將滲流問題變換直角網(wǎng)絡(luò)系統(tǒng)下進行有限差分離散和數(shù)值計算 ， 對不同的幾何邊界實現(xiàn)統(tǒng)一的數(shù)值求解算法 。 數(shù)