摘要
降雨是诱发自然地质灾害发生以及诸多岩土体工程破坏失稳的主要原因之一。因此研究降雨非饱和入渗机理及其计算方法具有重要的意义。目前对于降雨非饱和入渗的机理认识以及计算方法的探讨都存在一些不足之处。首先无论试验研究还是模型理论的建立都没有对小降雨强度的自由入渗阶段给与足够的重视;其次传统建立在入渗量基础上的降雨入渗模型无法很好解决岩土体工程的渗透稳定性问题,即工程土体内渗透荷载以及负压变化规律;再次,现将降雨入渗作为一种边界条件来处理的计算方法不符合降雨的物理过程且会带来较大误差。因此,本文针从水势的角度对降雨非饱和入渗展开了试验研究,并在此基础上进行了理论模型和计算方法的研究,希望能为降雨条件下的岩土工程问题性问题提供新的思路和计算方法。
本文以试验为基础,试验结果分析、机理研究、理论模型推导和数值计算方法相结合。首先设计了室内可精确控制雨强的降雨非饱和入渗试验装置,分别以粉土和粉砂土进行了包括自由入渗和压力入渗阶段的降雨试验,以及考虑了非饱和入渗参数影响的试验;基于试验结果从水势(负压)的角度对降雨非饱和入渗过程进行了描述,并进一步分析了降雨强度、非饱和入渗系数以及负压之间相互关系。现有的降雨非饱和入渗计算将降雨入渗的过程作为一种边界条件来处理的方法未能体现降雨入渗过程中几个关键物理量的变化关系及其随时间的变化关系。本文从负压的角度推导了降雨自由入渗和压力入渗模型,且应用降雨边界上的入渗模型联立计算方法来替代传统的边界条件处理方法。该方法符合降雨入渗的物理过程,且误差较小,具有较为明显的物理意义和很好的工程应用价值。
本文发展了模型相应的数值计算方法,并结合已开发的降雨非饱和计算软件Unsat-Seepage,编制了相应的程序。应用本文的计算方法,对交河故城土遗址崖体在考虑降雨入渗条件下的稳定性问题进行了计算分析,验证了成果的正确有效性。
Rainfall always induced catastrophic disasters and landslide, so it is a great significance of the study of rainfall infiltration mechanism and it is calculation methods. There are some shortages in current study about rainfall infiltration, the free seepage stage has always been neglected in study of rainfall experiments and theoretical models although the stage is very important to engineering seepage stability problems; it is reasonless that current calculation methods treat rainfall infiltration as a simple boundary, which not only brings more errors but also has not reflected the practical rainfall course. So the studies of rainfall unsaturated infiltration experiments、theoretical models and calculation methods have accomplished in this paper, hope to provide new research ways to rainfall infiltration stability problems of geotechnical engineering.
The main works include tests of rainfall infiltration equipments、results analysis、theoretical models study and numerical calculation. The test equipment of rainfall infiltration has been designed firstly, which can exactly control rainfall intensity. Twelve tests has been done using different soils as specimen which to illustrate the mechanism of free and pressure infiltration stages from a point of view of suction, and the rainfall course has been definitely described. Furthermore, we analyzed the relationship of rainfall intensity soil unsaturated parameters and suction. Based on suction, we conducted the rainfall free and pressure infiltration models, they can tell the essence of rainfall infiltration, and the research developed relative calculation methods to take place of the current method treating the rainfall infiltration as a boundary. All these works not only have clearly practical significance and good precision, but also is well adapt to engineering application.
The paper design the relative numerical calculation method based on the rainfall infiltration theoretical models we conducted, combining with rainfall infiltration compute program Unsat-Seepage, the computer codes have been developed to study the rainfall infiltration and its’engineering application. As an elucidation, the stability problems of Ruins of Jiaohe slopes under rainfall infiltration have been simulated, using the method developed the paper. The result showed that all works of the paper are more applicable and reasonable.
引文
[1]王恭先.面向21世纪我国滑坡灾害防治的思考.兰州滑坡泥石流学术研讨会文集.兰州:兰州大学出版社,1981:1-8.
[2]李广信.岩土工程的经济、安全与可持续发展--国家科技中长期发展规划建议.岩土工程界,2004,7(7):1-7.
[3]潘懋,李铁锋.灾害地质学.北京:北京大学出版社,2002.
[4]钟立勋,中国重大地质灾害事例分析,中国地质灾害与防治学报,1999.9
[5]王瑞久,水文条件对河谷岩石边坡的影响,水文地质工程地质,1999.6.
[6]孙广忠,中国典型滑坡,北京:科学出版社,1998.
[7]郑明新王金才,降雨滑坡的时空预测新方法,中国地质灾害与防治学报, 1997.9: 8 (3):25-29.
[8]中华人民共和国水利部. SL264-2001.中华人民共和国行业标准-水利水电工程岩石试验规程.北京:中国水利水电出版社, 2001.
[9] Malone A W,黄润秋,边坡安全与滑坡风险管理-香港的经验,地质科技管理,1999,第5期:6-18.
[10] U.S. Department of the interior, U.S. Geological Survey. Landslide types and processes [OL]. [2004].
[11] Bishop A W, Alpan I, Blight G E, and B. Factors contraling the shear strength of partly saturated cohesive soils. ACSE Research Conference on Shear Strength of Cohesive Soil, University of Colorado, Boulder, CO. 1960:503-532.
[12] W. R. Pucell. Capillary pressures their measurement using mercury and the calculation of permeability. AIME, 1949, V.186, P39-48.
[13]荣冠,王思敬,王恩志等.强降雨下元磨典型边坡工程稳定性分析[J].岩土力学, 2008, 27(4):704-711.
[14]王思敬.中国岩石力学与工程世纪成就.南京:河海大学出版社, 2004:21-35.
[15]陈祖煜.土质边坡稳定分析.北京:中国水利水电出版社,2003.
[16]孙广忠,中国典型滑坡,北京:科学出版社,1998.
[17]钟荫乾,滑坡与降雨关系及其预报,中国地质灾害与防治学报,1998.12:9 (4).
[18]李悼芬汪敏,港渝两地滑坡灾害的对比研究,岩石力学与工程学报,2000.7:19 (4).
[19] Rutqvist J, Wu Y S, Tsang C F, et al.A modeling approach for analysis of coupled multiphase fluid flow,heat transfer,and deformation in fractured porous rock.International Journal of Rock Mechanics and Mining Sciences, 2002,39(4):429-442.
[20] Dershowitz W.S., A new three dimensional model for flow in fractured rock. Int.Assoc.Hydrogeol., 1985, (17):441-448.
[21]朱伟,山村和也.雨水、洪水渗透时河堤的稳定性,岩土工程学报,1999.7:21 (4).
[22]陈学东.浅层非饱和带降雨入渗规律的试验与数值研究[硕士学位论文].南京:河海大学,2005.
[23] (英)E.HoeK J.WBray,岩石边坡工程,冶金工业出版社,1983.
[24]殷坤龙韩再生李志中,国际滑坡研究的新进展,水文地质工程地质,2000.5.
[25]谢新宇,朱向荣,谢康和,等.饱和土体一维大变形固结理论新进展.岩土工程学报, 1997, 19(4):30-38.
[26]龚壁卫詹良通刘艳华王斌,非饱和膨胀土的抗剪强度特性研究,长江科学院院报,2000.8:17(5).
[27]龚壁卫,非饱和击实膨胀土总应力强度探讨,长江科学院院报,1998.6:15(3).
[28] Fredlund D G,Morgenstern N R,Widger R A. Shear strength of unsaturated soils[J].Canadian Geotechnical Journal, 1978, 15(3):313–321.
[29] Fredlund D G,Rahardio H.非饱和土土力学[M].陈仲颐译.北京:中国建筑工业出版社,1997.
[30] Richards L A. Capillary conduction of liquids in porous mediums [J]. Physics, 1931,1: 318–333.
[31]雷志栋杨诗秀谢森传.土壤水动力学,清华大学出版社,1988.
[32] Green W H, and Ampt G A. Studies on soil physics:1,Flow of air and water though soil[J].J Agric Sci,1911,4(1):1~24.
[33] Van Genuchten M Th. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils[J]. Soil Science Society of American Journal, 1980, 44(5): 892-898.
[34]朱伟,程南军,陈学东等.浅谈非饱和渗流的几个基本问题[J].岩土工程学报,2006,28(2):235–240.
[35]李秀珍,许强,黄润秋,汤明.高滑坡预报判据研究中国地质灾害与防治学报,2003,14(4):173-184.
[36] Horton R E. An approach toward a physical interpretation of infiltration capacity[J].Soil Sci, Soc Am Proc,1940,3:399~417.
[37]陈正汉.重塑非饱和土的变形、强度、屈服和水量变化特性[J].岩土工程学报,1999.21(1):82-90.
[38] Philip J.R. Numerical Solution of Quations of the Diffusion Type With DiffusivityConcentration-Dependent. Australian J. Physics Vol 10, Nol ,1957.
[39] Childs E.C. the permeability of porous materials. Proceedings of the Royal Society of London. Series A. 1948, 201; 392-405.
[40]胡波.非饱和土-水特征曲线的研究.武汉大学.2005.
[41]詹良通,吴宏伟,包承纲,龚壁卫.降雨入渗条件下非饱和膨胀土边坡原位监测.岩土力学,2003,2:151-158.
[42]张文杰,詹良通,凌道盛,陈云敏.水位升降对库区非饱和土质岸坡稳定性的影响.浙江大学学报:工学版,2006,40(8):1365-1370.
[43] Khalili N, and Khabbaz M H, A unique relationship for the determination of the shear strength of unsaturated soil. Geotechnique, 1998,48(5):681-687.
[44] Vanapalli S K, Fredlund D G, Pufahl D E, Clfton A W. model for the prediction of shear strength with respect to soil suction. Canadian Geotechnical Journal,1996,33:379-392.
[45]缪林昌,殷宗泽非饱和土的剪切强度.岩土力学,1999,3:1-6.
[46] Richards L A. Capillary conduction of liquids through porous medium.Physics,1931,1:318-333.
[47] Gardner W R. Some steady-state solution of liquids through porous mediums. Physics,1931,1:228-232.
[48] Brooks R H, Corey A T. Hydraulic prosperities of porous medium. Fort Collins, Colorado. Hydrology Paper, 1964.
[49]姚海林,郑少河,李文斌等.降雨入渗对非饱和膨胀土边坡稳定性影响的参数研究[J].岩石力学与工程学报。2002,21(7):1034-1039.
[50]栾茂田,李顺群等.非饱和土的理论土水特征曲线[J].岩土工程学报.2005,27(6):611-615.
[51] Leong E C , Rahardjo H. Review of soil water characteristic curve equations. Journal of Geotechnical , 1994 31:521-532.
[52] Leong E C , Rahardjo H. Permeability functions for unsaturated soils. Journal of Geotechnical and Geo-enviromental Engineering ASCE, 1997,123(12) :1118-1126.
[53]陈文亮,王占礼.人工模拟降雨特性的试验研究.水土保持通报,1991,11(2):55-62.
[54]赵志进,李桂英.人工模拟降雨机具和方法的发展研究与展望.中国水土保持.1991,11(2):1-4.
[55]吴林高.渗流力学.上海:上海科技文献出版社, 1996.
[56]顾慰慈.渗流计算原理及应用.北京:中国建材工业出版社, 2000.
[57]刘宁.渗流力学基础.北京:石油工业出版社, 1985.
[58]吴林高.渗流力学.上海:上海科学技术文献出版社, 1996.
[59]冯康.数值计算方法.北京:国防工业出版社, 1978.
[60]薛禹群.地下水动力学.北京:地质出版社, 1997.
[61]王思敬.地球内外动力耦合作用与重大地质灾害的成因初探.工程地质学报, 2002, 10(2):115-117.
[62] Uyakon P S. Pinder G F. Computation Methods in Subsurface Flow. New York : Academic Press, 1983.
[63]周佩华,张学栋,唐克丽.黄土高原土壤侵蚀与旱地农业国家重点实验室土壤侵蚀模拟实验大厅降雨装置.水土保持通报. 2000, 20(4):27-31.
[64]刘震,郝振纯.人工模拟降雨系统总体设计.水利水电技术. 1998,29(8):1-4.
[65]徐向舟,刘大庆,张红武,董占地,朱明东.室内人工模拟降雨试验研究,北京林业大学学报,2006,28(5):52-58.
[66]孙大松,刘鹏,夏小和,王建华非饱和土的渗透系数.水利学报,2004,3:71-75.
[67]刘海宁,姜彤,刘汉东.非饱和土渗透函数的间接确定.岩土力学,25(11):1795-1799.
[68]徐永福.非饱和土本构模型研究综述。水利水电科技进展.1996,16(5):4-9.
[69]张有天陈平,有自由面渗流分析的初流量法,水利学报,1988.8.
[70]朱岳明,Darcy渗流量计算的等效结点流量法,河海大学学报,1997.7.
[71]梁业国熊文林,有自由面渗流分析的子单元法,水利学报,1997.8 :34-38.
[72]陈万吉汪自力,瞬态有自由面渗流分析的不动网格一高斯点有限元法,大连理工大学学报,1991.9.
[73]速宝玉沈振中,用变分不等式理论求解渗流问题的截止负压法,水利学报,1996.3
[74]陈守义,考虑入渗和蒸发影响的土坡稳定性分析方法,岩土力学,1997.6.
[75] R.A.弗里泽J.A.彻里,地下水,地震出版社,1987.12.
[76]庄一鸽林三益,水文预报,水利电力出版社,1998.
[77]姚贤良程云生,土壤物理学,农业出版社,1986.
[78] [苏] K.C.巴斯宁耶夫, A.M.费拉西莫夫等编著.地下流体力学.北京:石油工业出版社, 1990.
[79]郭禄光,樊功渝.最小二乘法与测量平差.上海:同济大学出版社, 1985.
[80]陈宝林.最优化理论与方法.北京:清华大学出版社, 1989.
[81]谢政,李建华,汤泽滢.非线性最优化.长沙:国防科技大学出版社, 2006.
[82]席少霖.非线性最优化方法.北京:高等教育出版社, 1992.
[83]张光澄,王文娟,韩会磊,等.非线性最优化计算方法.北京:高等教育出版社, 2005.
[84]潘祖梁.非线性问题的数学方法及其应用.杭州:浙江大学出版社, 1998.
[85]张韵华,奚梅成,陈长松.数值计算方法和算法.北京:科学出版社, 2000.
[86]蔡大用,白峰杉.高等数值分析.北京:清华大学出版社, 1997.
[87]李庆扬,王能超,易大义.数值分析.北京:清华大学出版社, 2008.
[88] Brand E W,Premchitt J, Phillipson H B. Relationship between Rainfall and Landsildes. Proceedings of the fourth International Symposium on Landslides, Toronto,1,Ontario Bitech, Vancourver, Canada, 1984:377-384.
[89]黄康乐.求解非饱和土壤水流问题的一种数值方法[J].水利学报,1987,(9):9-16.
[90]朱学愚,谢春红,钱孝星.非饱和流动问题的SUPG有限元法[J].水利学报,1994, (6): 37-42.
[91] Neuman S P. Saturated-unsaturated seepage by finite element. J.Hydraul. Div, Am, Soc. Civ. Eng. 1973, 99(12):2233-2250.
[92] Gallipoli D. Constitutive and numerical modeling of unsaturated soils. Ph.D. thesis, University of Glassgow, UK. 2000.
[93] Petterson K E. The early history of circular sliding surface. Geotechnique,1955,5:275-296.
[94] Deng G, Shen Z J, Chen T L, Yang D Q. Numerical simulation of slope failure during/after rainfall. In: Proceeding of Int. Conf. on Slope Engineering, HongKong, 2003:207-212.
[95]邵明安王全九,推求土壤水分运动参数的简单入渗法/理论分析,土壤学报, 2000.2:37(1).
[96]张家发李思慎叶自桐.高边坡山体饱和非饱和渗流场的初步分析,人民长江,1998.1:20(1).
[97]吴梦喜,高莲士.饱和-非饱和土体非稳定渗流数值分析[J].水利学报, 1999, 12:38-42.
[98]彭华,饱和-非饱和岩土非稳定渗流有限元分析研究[J],水动力学研究与进展。2002,17(2):254-259.
[99]周剑平. Origin实用教程7.5版.西安:西安交通大学出版社, 2007.
[100]叶卫平. Origin7.5科技绘图及数据分析.北京:机械工业出版社, 2006.
[101]张昌松. Origin7.5入门与提高.北京:化学工业出版社, 2008.
[102]孙冬梅,朱岳明,张明进.非饱和带水–气二相流数值模拟研究[J].岩土工程学报,2007,29(4):560–565
[103]朱伟,程南军,陈学东等.浅谈非饱和渗流的几个基本问题[J].岩土工程学报,2006,28(2):235–240.
[104]林鸿洲.降雨诱发土质边坡失稳的试验与数值分析研究[博士学位论文].北京:清华大学水利水电工程系,2004.
[105]阎岩.渗流作用下岩石蠕变试验与变参数方程的研究[博士学位论文].北京:清华大学水利水电工程系,2009.
[106]刘晓丽.水岩耦合过程及其多尺度行为理论与应用研究[博士学位论文].北京:清华大学水利水电工程系,2008.
[107]梁爱民.非饱和土壤渗透特性及饱和入渗机理试验研究[博士学位论文].大连:大连理工大学,2008.
[108] Pruess K,Narasimhan J S Y &Tsang Y W On thermchydralogic conditions near high-level nuclear wastes emplaced in partially saturated fractured tuff,2 Effective continuum approximation Water Resour Res,1990,26(6):1249~1261.
[109]薛禹群,朱学愚.地下水动力学,地质出版社,1979.7.
[110]刘蕴华,张乃良。数学物理方程及其应用。南京:河海大学出版社,1994.
[111]郑颖人,赵尚毅.有限元强度折减法在土坡与岩坡的应用[J].岩石力学与工程学报,2004,23(19):3 381–3 388.
[112]连镇营,韩国城,孔宪京.强度折减有限元法研究开挖边坡的稳定性[J].岩土工程学报,2001,23(4):406–411.
[113] Kueper B H, C S Haase,and H King, Leakage of DNAPL from waste impoundments constructed in fractured rock and clay,. Theory and case History. Can. Geotech. J, 1992. 29.
[114]邓建辉,张嘉翔,闵弘等.基于强度折减概念的滑坡稳定性三维分析方法(II):加固安全系数计算[J].岩土力学,2004,25(6):871–875.
[115]郑宏,李春光,李焯芬,等.求解安全系数的有限元法[J].岩土工程学报,2002,24(5):626–628.
[116]栾茂田,武亚军,年廷凯.强度折减有限元法中边坡失稳的塑性区判据及其应用[J].防灾减灾工程学报,2003,23(3):1–8.
[117]赵尚毅,郑颖人,时卫民等.用有限元强度折减法求边坡稳定安全系数[J].岩土工程学报,2002,24(3):333–336.
[118]刘金龙,栾茂田,赵少飞,袁凡凡,王吉利.关于强度折减有限元方法中边坡失稳判据的讨论.岩土力学,2005,26(8):1345-1348.
[119] Mckee C R, Bumb A C, Koenig R A. Stress-dependent permeability and porosity of coal and other geologic formations. SPE Formation Evaluation, 1988,3(1):81-91.
[120] Thomas L K, Chin L Y, Pierson R G, et al. Coupled geomechanics and reservoir simulation. SPE Journal, 2003,8(4):350-358.
[121] Rosenblum L, Earnshaw R A, Encarnacao J, et al. Scientific Visualization: Advances and Challenges. New York: Academic Press, 1994.
[122] Dershowitz W S,Einstein H H. Three-dimensional flow modeling in jointed rock masses[A].// Proc. the 6th Int. Cong.,ISRM[C], 1987. 87-92.
[123] Wang J S Y, and T N Narasimhan. Hydrologic mechanisms governing fluid flowin a partially saturated, fractured, porous medium,. Water Resour. Res.,1985, 21 (12).
[124]李荣建,于玉贞,李广信.有限元强度折减法在非饱和土坡稳定分析中的应用.水利水电技术,2006, 37(3):42-45.
[125] Itasca Consulting Group Inc.FLAC3D user guide(version 3.0)[R].Minneapolis,Minnesota,USA: Itasca Consulting Group Inc.2005.
[126]毛昶熙等。渗流数值计算与程序应用。南京:河海大学出版社,1999.
[127]荣冠,张伟,周创兵.降雨入渗条件下边坡岩体饱和非饱和渗流计算[J].岩土力学,2005,26(10):1 545–1 550.
[128]朱伟,秦建设,高玉峰.求解非饱和土堤初始浸润线的一种解析法[J].河海大学学报(自然科学版),2003,31(3):314–317.
[129]吴梦喜,张学勤,有自由面渗流分析的虚单元法,水利学报,1994.8.
[130] Louis C A. Study on groundwater flow in jointed rock and its influence on the stability of rock masses[R]. London:Imperial College,1969.
[131] Louis C. Rock hydraulics. // Muller L. Rock mechanics. 1974: 299-387.
[132] Jones F O. A laboratory study of the effects of confining pressure on fracture flow and storage capacity in carbonate tocks. J.Petrol.Tech., 1975,21.
[133] Snow D T, Rock fracture spacings, openings and porosities. J. Soil. Mech. Found. Div., Proc. ASCE, 1968, 94(SMI):73-91
[134] Gale J E., The effects of fracture type ( Induced versus natural) on the stress-fracture closure permeability relationships. In: Proc. 23th Symp. On Rock Mech., Berkeley, California, 1982.
[135]吴宏伟,陈守义,庞宇威.雨水入渗对非饱和土坡稳定性影响的参数研究,岩土力学,1999.3,20 (1).
[136]王金生,杨志峰,陈家军,王志明.包气带土壤水分滞留特征研究,水利学报, 2000.2 (2).
[137] Peters R R and E A Klavetter. A continuum model for water movement in unsaturated fractured rock mass. Water Resour. Res, 1998, 24 (3).
[138]张有天,陈祖煌,刘中等.岩体渗流对三峡船闸高边坡稳定的影响[R],北京:清华大学中国水利水电科学研究院中国科学院地质研究所,1995.
[139] Istok J. Groundwater Modeling by the Finite Element Method. Washington DC: American Geophysical Union, 1989.
[140] Bach L B..Soil water movement in response to temperature gradients: Experimental measurements and model evoluation. Soil Sci.Soc.Am.J. 1992, 56:37-46.
[141]周创兵,叶自桐,韩冰.岩石节理非饱和渗流特性初步研究,岩土工程学报,1998.11:20 (6).
[142]李锡夔,范益群.非饱和土变形及渗流过程的有限元分析.岩土工程学报, 1998, 20(4):20-24.
[143]孔祥言.高等渗流力学.合肥:中国科学技术大学出版社, 1999.
[144]刘建军,薛强.岩土工程中的若干流-固耦合问题.岩土工程界, 2004, 7(11): 27-30.
[145]叶金汉.岩石力学参数手册.北京:水利电力出版社, 1991.
[146] Zaman M M, Abdulraheem A, Roegiers J C. Reservior compaction and surface subsidence in the North Sea area. Journal of Geology. 1995, (5):373-379.
[147]赵慧丽,张弥,李兆平.含水量对北京地区非饱和土抗剪强度影响的试验研究,石家庄铁道学院学报2001,Vol.14(4):30-33.
[148]戴福初,李悼芬,王思敬.松散击实火山岩坡残积土的应力应变特性及其对滑坡的意义,岩土工程学报,1999.5:21 (3).
[149] Gerke H H,and M T Van Genuchten. A dual-porosity model for simulating thepreferential movement of water and solute in structured porous media. WaterResour. Res.,1993,29(2).
[150] Aarnes J, Hou T Y. Multiscale Domain Decomposition Methods for Elliptic Problems with High Aspect Ratios. Acta Mathematicae Applicatae Sinica, 2002, 18(1): 63-76.
[151] Efendiev Y, Hou T Y, Luo W. Preconditioning of Markov Chain Monte Carlo Simulations Using Coarse-Scale Models. SIAM Journal of Scientific Computing, 2006, 28(2): 776-803.
[152] Efendiev Y, Hou T Y, Strinopoulos T. Multiscale Simulations of Porous Media Flows in Flow-Based Coordinate System. Comput. Geoscience , 2008, 12:257-272.
[153]韩冰,叶自桐,周创兵.岩体裂隙饱和/非饱和渗流机理初步研究,水科学进展,1999, 12:10 (4).
[154] Bai Y L, Xia M F, Wei Y J, et al. Non-Equilibrium Evolution of Collective Microdamage and Its Coupling with Meso-scopic Heterogeneities and Stress Fluctuations. // Horie Y, Davisona L, Thadhani N N. Shock Dynamics and Non-Equilibrium Mesoscopic Fluctuations in Solids. New York: Springer-Verlag, 2002: 255-278.
[155]王思敬.岩体工程地质力学中的数值分析.岩体工程地质力学问题(八),科学出版社, 1987.
[156]王思敬.工程地质学的任务与未来.工程地质学报. 1999, 7(3):195-199.
[157]王思敬.环境工程地质评价、预测与对策分析.地质灾害与环境保护, 1997,8(1):27-34.
[158]王思敬.水库地区的水岩作用及其地质环境影响.工程地质学报, 1996, 4(3):1-9.
[159] Hou T Y, Wu X H. A Multiscale Finite Element Method for Elliptic Problems in Composite Materials and Porous Media. J. Comput Phys. , 1997, 134:169-189.
[160] Milly C D, Advances in modeling of water in the unsaturated zone, Transport in Porous Media,1998,3:491-514
[161] Pcck A J., Moisture profile development and air compression during water uptake by bounded porous bodices:2. Horizontal columns, Soil Sci. 1965a,100:333-340
[162]张辉.裂隙孔隙岩体非饱和渗透率的探讨.云南水力发电, 2002(1):54-57.
[163] Harrison J P, Hudson J A. Engineering rock mechanics. Part 2:illustrative workable examples. In: S.arkk.a P, Eloranta P, editors. Oxford: Pergamon, 2000.
[164] Hudson J A. Rock engineering case histories: key factors, mechanisms and problems. Sarkka, E. Rock Mechanics-a challenge for society. Rotterdam: Balkema, 2001:13-20.
[165]王恩志,孙役,邓旭东.单裂隙非饱和渗透荷载分布规律,清华大学学报(自然科学版),2000,40 (4).
[166] (美)YANG H.HUANG著,包承纲王清友宋昆仑涂帆等译,土坡稳定分析,清华大学出版社,1988.10.
[167]王恩志,杨成田.裂隙网络及地下水网络数值方法的研究,勘察科学技术,1991(4):13~16.
[168]王恩志.岩体裂隙的网络分析及渗流模型.岩石力学与工程学报. 1993,12(3):214~221.
[169]王恩志,王洪涛,孙役.三维裂隙网络渗流数值模型研究,工程力学,1997(增2):520~525.
[170]王恩志,王洪涛,孙役.双重裂隙系统渗流模型研究,岩石力学与工程学报,1998,8,Vol17(4):400~406.
[171]王恩志,王洪涛,邓旭东.‘以管代孔’-岩体排水孔模拟方法探讨,岩石力学与工程学报,2001年,21(3):346~349.
[172]王恩志,王洪涛,王慧明.‘以缝代井列’-岩体排水孔幕模拟方法探讨,岩石力学与工程学报,2002年,21(1):98~101.
[173]王恩志,孙役,黄远智.三维离散裂隙网络渗流模型与实验模拟,水利学报,2002年,307(5):37~40.
[174]王恩志,邓旭东,王远亮,魏寿松.天生桥一级水电站三维岩体渗流分析,水力发电学报,2000年,71(4):25~32.
[175]王恩志,李昂,毛文然,吴海亮.回龙抽水蓄能电站上库岩体三维渗流分析,水力发电学报,2002年,76(1):96~101.
[176] Farhat C, Roux F X. A method of finite element tearing and interconnecting and its parallel solution algorithm. Int. J. Num. Meth. Eng., 1991,32(6):1205-1227.
[177] Huhes T J R, Levit J, Winget J. Eelement-by-element implicit algorithms for heat conduction analysis. Journal of Engineering Mechanics, 1983, 109(6):576-585.
[178] Kazemi H,Merrill J R, Porterfield K L, et al. Numerical simulation of water-oil flow in naturally fractured reservoir. Trans. AIME, 1976, 261:317-326.
[179]王洪涛,王恩志.各向异性裂隙岩体渗透系数计算方法探讨,武汉水利电力大学学报,1997,30(2),49~53.
[180]王洪涛,王恩志,岩体主干裂隙系统三维非稳定渗流分析模型,水动力学研究与进展,1998,A13(2),206~213.
[181]王洪涛,王恩志,裂隙岩体渗流耦合数值分析方法及其工程应用,水利水运科学研究,1998年,第4期,35~40
[182]孙满利,王旭东,李最雄等。交河故城瞭望台保护加固技术[J].岩土力学, 2007,28(1):163-168.
[183]解耀华.交河故城的历史与保护修缮工程[A].
[184]田琳.交河故城生土建筑遗存保护加固实验[A].
[185]解耀华编.交河故城保护与研究[C].乌鲁木齐,新疆人民出版社,1999.
[186]刘波,韩彦辉. FLAC原理、实例与应用指南.北京:人民交通出版社, 2005.