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