【正文】 根據(jù)三峽水庫調(diào)度運(yùn)營規(guī)劃 ， 水庫建造完成投入運(yùn)行后 ， 洪期水位保持在 145m 附近 ， 枯水季節(jié)水位保持在 175m附近 ， 每年雨季來臨前降水空庫確保防洪 。 朝倫巴根和 劉延璽等用泛克里格法 ， 對(duì)通遼河西水資源評(píng)價(jià)與管理研究區(qū)內(nèi)含水層系統(tǒng)的滲透系數(shù)進(jìn)行了結(jié)構(gòu)分析 。 Bulness( 1946) 和 Warren( 1961) 研究油田巖芯資料時(shí) ， 證實(shí)了 Law 的發(fā)現(xiàn) 。 固定網(wǎng)格法又可以分成兩種 ： 一種是調(diào)整單元滲透矩陣法 ， 如 Bathe 的單元滲透矩陣調(diào)整法 ， 李春華的改進(jìn)單元滲透矩陣調(diào)整法 ， 王賢能等的高斯點(diǎn)法 。 常用的方法有砂槽模型法、縫隙水槽法、水力網(wǎng)模型法、電擬試驗(yàn)法和電阻網(wǎng)模型法 。 雙重介質(zhì)模型的基本假定是裂隙 巖體中孔隙、裂隙分布均勻 ， 即在微元體中 ，孔隙與裂隙體積之比保持不變 ， 認(rèn)為裂隙巖體由孔隙性差而導(dǎo)水性強(qiáng)的裂隙系統(tǒng)和孔隙性好而導(dǎo)水性差的巖塊孔隙系統(tǒng)共同構(gòu)成的連續(xù)介質(zhì) 。之后考慮滲流作用的邊坡應(yīng)力 變形關(guān)系的有限元數(shù)值模擬技術(shù)發(fā)展起來 ， 并出現(xiàn)了一些新方法 ， 如張孟喜 (1991)提出了有限元追蹤法 、 邵龍?zhí)?(1997)提出了水流邊坡穩(wěn)定法 、 毛趙熙 (20xx)提出了圓弧有限元法等 ； 童富果 (20xx)提出在一定的滲流場(chǎng)下 ,可以將滲流場(chǎng)產(chǎn)生的滲透力作為荷載 ， 分析 坡體的應(yīng)力場(chǎng) ， 進(jìn)而評(píng)價(jià)坡體穩(wěn)定性 ； 田杰 (20xx)認(rèn)為由于坡體應(yīng)力場(chǎng)的存在和改變導(dǎo)致多孔介質(zhì)的孔隙率發(fā)生改變 ， 從而導(dǎo)致其滲透系數(shù)變化 ， 進(jìn)而對(duì)滲流場(chǎng)產(chǎn)生影響 ， 而滲流場(chǎng)產(chǎn)生的滲透壓力又直接影響了應(yīng)力場(chǎng) ， 因此基于比奧固結(jié)理論分析了滲流場(chǎng) 應(yīng)力場(chǎng)的 耦合 作用 。 非確定性分析方法是根據(jù)滑坡預(yù)測(cè)理論的廣義系統(tǒng)科 學(xué)原理 ， 在類比法的基礎(chǔ)上應(yīng)運(yùn)而生的一類研究方法 。 發(fā)展至今主要有以下幾種方法 ： 有限單元法 (FEM)： 該方法未將坡體視為剛體 ， 考慮了邊坡體的非均質(zhì)性 、不連續(xù)性等特征 ， 并綜合考慮了巖土體的應(yīng)力 、 應(yīng)變發(fā)展規(guī)律 ,能夠準(zhǔn)確地以應(yīng)力 、 應(yīng)變?yōu)橐蛩胤治鲞吰伦冃纹茐牡陌l(fā)生機(jī)制 ,對(duì)了解邊坡的應(yīng)力分布 、 變形發(fā)展有很大幫助 。 本文對(duì)國內(nèi)外滲流及邊坡穩(wěn)定的研究現(xiàn)狀進(jìn)行了歸納總結(jié)。這樣的滲流場(chǎng)相互作用，將對(duì)邊坡 穩(wěn)定性 產(chǎn)生 影響。 從本質(zhì)上講 ,塑性極限分析法和極限平衡法一樣 ， 都只是靜態(tài)解 ，只能反映某個(gè)特定時(shí)空下邊坡的穩(wěn)定狀況 ， 不能反應(yīng)邊坡巖土體的應(yīng)力 變形發(fā)展變化規(guī)律 ， 更不能模擬邊坡失穩(wěn)破壞的動(dòng)態(tài)過程 。 在邊坡穩(wěn)定性分析中數(shù)值方法與傳統(tǒng)的極限平衡法相比雖然具有諸多優(yōu)勢(shì) ，但由于 地質(zhì)條件的復(fù)雜性和認(rèn)識(shí)的局限性 ， 建立的地質(zhì)模型和采用的巖土體本構(gòu)關(guān)系往往與實(shí)際情況存在著一定的差異 ， 且計(jì)算參數(shù)的選取也可能是以某種簡化為基礎(chǔ)的 ， 因而與實(shí)際情況會(huì)有誤差 ， 從而影響了計(jì)算結(jié)果的精度 。 數(shù)值方法發(fā)展到今天 ， 其常用的方法有 有限元法 、 邊界元法 、 有限差分法等等 。 離散模型是應(yīng)用最廣的一種模型 ， 但難于用于進(jìn)行復(fù)雜的三維滲流分析 。 介玉新等將適體坐標(biāo)變換方法用于滲流計(jì)算 ， 即通過求Poisson 方程自動(dòng)生成實(shí)際物理邊界的計(jì)算網(wǎng)絡(luò) ， 并將滲流問題變換直角網(wǎng)絡(luò)系統(tǒng)下進(jìn)行有限差分離散和數(shù)值計(jì)算 ， 對(duì)不同的幾何邊界實(shí)現(xiàn)統(tǒng)一的數(shù)值求解算法 。 滲流自由面 第一是移動(dòng)網(wǎng)格法 : 這類方法把自由面當(dāng)作可移動(dòng)邊界處理 ， 在迭代過程中不斷改變自由面的位置 ， 使網(wǎng)格發(fā)生變化 ， 直到自由面位置穩(wěn)定為止 。Manath N. Panda 等在研究用粒度分布參數(shù)計(jì)算松散介質(zhì)的滲透率中指出 ， 對(duì)非均質(zhì)介質(zhì)而言 ， 介質(zhì)粒度均值的變化對(duì)滲透率的貢獻(xiàn)比孔隙度或分選程度要大 。 大多數(shù)研究者選擇隨機(jī)模型作為研究手段 ， 是假定含水介質(zhì)滲透系數(shù)滿足對(duì)數(shù)正態(tài)分布 ， 但缺乏完備的理論和充分的實(shí)踐資料支持 ， 因此其應(yīng)用受到很大的限制 。 主要側(cè)重于理論的引用 ， 鮮見與野外鉆孔、油水實(shí)驗(yàn)資料結(jié)合對(duì)松散介質(zhì)滲透系數(shù)進(jìn)行空間估值的實(shí)例 。 2) 庫位變動(dòng)帶飽和滲流、非飽和滲流相互轉(zhuǎn)化過程研究 ， 明確岸坡巖土體在不同時(shí)間、不同空間部位的滲透系數(shù)取值規(guī)律 。 數(shù)學(xué)模型反演求解法 運(yùn)用地下水動(dòng)態(tài)信息 ， 用建立的數(shù)學(xué)模型來反演求解 。 其基本思想是將研究區(qū)劃分為不同的塊段 。 王賢能則根據(jù)高斯點(diǎn)的位置確定高斯點(diǎn) 對(duì)應(yīng)的滲透系數(shù) ， 計(jì)算單個(gè)高斯點(diǎn)對(duì)滲透矩陣的貢獻(xiàn) ， 對(duì)所有高斯點(diǎn)求和得到單元的滲透矩陣 ； 第二類是調(diào)整流量法 ， 如 Desai 提出的剩余流量法 ， 張有天的初流量法等 。 法國專家 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)系 。 建立了由管狀線單元、縫狀面單元和帶狀體單元組合而成的三維裂隙網(wǎng)絡(luò)滲流數(shù)值模型 。 Sammoril(1991)等人對(duì)滲流方程采用 Richard 方程 ， 提出了穩(wěn)定性分析采用的簡化 Bishop 法 。 當(dāng)前應(yīng)用較多的非確定性定量分析方法主要有 ： 經(jīng)驗(yàn)方法 ， 數(shù)理統(tǒng)計(jì)方法 ， 信息模型法 ， 模糊數(shù)學(xué)評(píng)判法 ， 灰色系統(tǒng)方法 ， 模式識(shí)別方法 ， 非線性模型預(yù) 測(cè)法 ， 人工智能法等 。 此外 ，裂隙塊體介質(zhì)的大變形及破壞問題的分析也適用于該法 。該方法的優(yōu)點(diǎn)是綜合考慮了影響邊坡穩(wěn)定的各種因素，能快速的對(duì)其穩(wěn)定性做出評(píng)價(jià)和預(yù)測(cè)。 滲流場(chǎng)和應(yīng)力場(chǎng)作為邊坡穩(wěn)定分析的重要組成部分，二者相互作用、相互聯(lián)系。該方法要經(jīng)過多次試算，并且最后得出的滑動(dòng)面也并不一定是邊坡實(shí)際發(fā)生滑動(dòng)的滑動(dòng)面；第二，僅滿足 力和力矩的平衡，在滑動(dòng)面上滿足莫爾 庫倫破壞準(zhǔn)則和應(yīng)力的邊界條件，無法反映巖土體自身的應(yīng)力 應(yīng)變關(guān)系；第三，難以分析較復(fù)雜的土層，土體填筑、開挖等施工活動(dòng)對(duì)邊坡的影響也難以分析；第四，只能分析靜態(tài)條件下邊坡體的穩(wěn)定狀態(tài)，不能夠模擬邊坡失穩(wěn)的整個(gè)動(dòng)態(tài)過程；第五，分析中沒有把水作為一種變量考慮其中，因而不能反映水對(duì)邊坡體穩(wěn)定的影響。 該方法較有限元方法能更全面地考慮巖土體的不連續(xù)性 、 各向異性及大變形特征 ，且求解速度較快 。 以解析法為主要研究手段的發(fā)展階段 隨著地下水開采量的越來越大 ， 人們發(fā)現(xiàn)地下水的運(yùn)動(dòng)狀態(tài)表現(xiàn)出明顯的隨時(shí)間變化的特征 ， (1928)研究了地下水運(yùn)動(dòng)的不穩(wěn)定性和承壓含 水層貯水性質(zhì) ， (1935)在此基礎(chǔ)上提出了地下水流向承壓水井的非穩(wěn)定流公式 。 等效非連續(xù)介質(zhì)模型亦稱離散模型 ， 忽略了巖塊的孔隙系統(tǒng) ， 把巖體視為簡單的裂隙介質(zhì) ， 用裂隙水力參數(shù)或幾何參數(shù)來表征裂隙巖體內(nèi)滲透空間結(jié)構(gòu)的具體布局 。 數(shù)值法 數(shù)值方法是把整個(gè)滲流區(qū)分割成若干個(gè)形狀規(guī)則的單元 ， 各個(gè)單元根據(jù)需要選擇合適的水文地質(zhì)參數(shù) ， 把所有單元合在一起就能表現(xiàn)出滲流區(qū)域在幾何上的不規(guī)則性和在水文地質(zhì)上的非均質(zhì)性 ， 代表實(shí)際的滲流區(qū) 。 迄今 ， 滲流場(chǎng)和應(yīng)力場(chǎng)耦合分析方法的研究較少 ， 一般都將兩場(chǎng)分開計(jì)算