青藏高原多年冻土区斜坡稳定性研究
|
摘要
青藏高原地区冻土分布面积达1.5×10~6 km~2，占我国多年冻土分布面积的70％，在全球气候转暖的背景下，多年冻土有着退化的趋势。而青藏高原地区的各类工程活动则进一步加剧了多年冻土退化的趋势，由此引发的冻融灾害不仅威胁着公路、铁路、桥梁、输油管道等线性工程的建设，而且导致寒区生态环境恶化，使多年冻土丧失恢复能力，热融滑塌就是最为典型的斜坡冻融灾害。
本文从青藏高原多年冻土区冻土斜坡失稳类型分析、典型冻土斜坡地质结构分析入手，通过对青藏公路K3035里程典型斜坡热融滑塌模型试验研究、稳定性分析、现场监测、温度场和变形场数值分析，开展了青藏高原多年冻土区斜坡稳定性研究，取得了如下成果：
① 在全面总结青藏公路沿线有关斜坡研究成果的基础上，通过调查对青藏高原多年冻土区冻土斜坡失稳类型进行了划分，根据其形成原因和冻土赋存环境，正融土滑坡可以划分为两种类型：融冻泥流型和热融滑塌型；正冻土滑坡可划分为两种类型，即蠕变型滑坡和崩塌型滑坡。
② 在青藏高原多年冻土区斜坡地带活动层土体沿确定的滑动面(冻融界面、厚层地下冰面)在融化季节滑动形成融冻泥流和热融滑塌两类典型的滑坡灾害。热融滑塌的失稳因素可以归结为：a．厚层地下冰的存在且埋藏深度较浅易受地表气温变化的影响，形成冰土分界面相对隔水且抗剪强度极低的滑坡条件，这是热融滑塌形成的内在因素；b．不合理取土或自然营力侵蚀引起的厚层地下冰的暴露且使坡体形成临空面而使其失去支承力，这是热融滑塌形成的外在因素(诱发因素)。
③ 青藏公路K3035热融滑塌的显著特点是滑塌体不是整体一次性滑动，目前的滑坡范围扩展是从后缘原坡体融化、开裂、坍塌，塌落的土体沿滑面滑动。扩展发育时期为每年气温正值期，且主要集中在0～2m深度土体范围内，处于融化状态的7-10月。综合2000～2002年3年内滑坡范围扩展情况，估算现今滑坡范围向后扩展的速度为5m／a。
④ 任何滑坡的产生，水分条件的变化是其重要的影响因素。热融滑塌更是如此，现场监测结果表明：融化季节在青藏公路K3035滑体2m以上范围的土体含水量高于冰层之上的粉质粘土液限值；在青藏公路K3057处大于同期对未滑动活动层土体的含水量的测量结果，即滑体水分显著高于原土体。冰面水分的聚集主要是由地下冰的融化所引起，正是由于大量融化水的存在，使得土体强度显著降低，并沿地下冰面产生了滑动。
⑤ 根据相似理论第一定律，应用积分类比法推导并建立了冻土斜坡模型试验的相似指标和相似判据，得出以原状土作为模型介质时，6个相似常数可以减少为2个相似常数，即c_τ和c_ι。模型与原型的时间比例尺是由几何比例尺决定的，即c_τ=c_ι~2，仅有一个相似参数为自变量，另一个为因变量。应用自行设计的变形测量装置对斜坡土体表面冻胀和融沉进行了计算和分析，得出冻土斜坡的融沉变形远远大于冻胀变形的结论，与青藏公路沿线多年冻土地区公路路基变形观测结果一致。通过模型试验，得出变形场受温度场控制的结论，并在多个冻融循环中，存在着一个斜坡后继变形的“温度门限”，控制着斜坡变形
的进一步发展。这些成果填补了多年冻土区滑坡研究的空白，为青藏高原多年
冻土区道路工程建设中滑坡整治提供了依据。
⑥本文推导了相变热传导问题微分控制方程的一般形式，应用有限元法对
青藏公路K3035斜坡地温场和变形场进行了模拟，建立了冻土斜坡温度场和变
形场数值分析模型。模拟了滑坎后缘拉张区5年内位移变化情况，得出冻土滑
坡体每一年中位移由两个阶段即稳定阶段和非稳定阶段组成的结论，稳定阶段
对应台阶面，非稳定阶段对应台阶壁;每一个台阶面均向右微微“上翘”，可
以得出土体在冻结过程中还发生上坡位移，而不单单仅存在下坡位移，这一计
算结果与风火山冻土蠕变监测结果取得了很好的一致。从模拟坡面内不同深度
对应的下坡位移可以得出冻土斜坡的位移主要发生在浅部，深部基本不发生位
移的结论，与青藏公路K3035斜坡现场测斜监测资料一致。采用基于灰色系统
理论的GM(1，1)模型，对滑塌距离动态变化进行预测，计算滑塌距离为5.76m，
误差为12.94%，说明所建立的GM(1，1)模型可以作为青藏公路K3035滑塌距离
年内月平均变化的预测模型。
⑦应用有效应力原理推导了不同渗流条件下无限斜坡稳定性分析的计算
公式，对各种不同表达形式的安全系数进行了讨论和对比，得出在流线平行于
斜坡条件下，两类特殊条件(干土坡、饱水土坡)下的分析结果与相同情况下
安全系数统一式结果相一致的结论;当已=O及0(m(l时，无限斜坡安全系数的
通式为尸=[(1一m)+m
上]竺旦竺
ysat tana
流线为相互平行的水平线、流线与坡向一致两种
情况分别为孔隙压力和安全系数的两种特例。流线与水平线的夹角为。‘刀‘a，
二，.、二二/。、__刀_____:__:____翔。二。，，、。卜、、、。
孔隙水压力u(刀)=mzy*一二mzr，sin，a，刀一。和刀=a分别对应上述两种极
口
值情况(流线为相互平行的水平线、流线与坡向一致)。
⑧绘制出了不同坡度、不同含水条
Permafrost area of Qinghai-Tibet Platean is 1. 5x106 km2,70 percent of permafrost area in China. On the basis of global warming, there is a degradation tendency for permafrost. The degradation tendency of permafrost will become worse caused by all kinds of human engineering activities. Geological hazards from freeze-thaw action will not only threat linear engineering construction for example, highway, railway, bridge, oil pipeline, etc et, but result in deterioration for ecological environment of cold region to lose recovery ability for permafrost. Thermal-thaw sliding collapse is a typical example.
The following researches have been carried out in this paper, including slope instability patterns analysis in permafrost region of Qinghai-Tibet Platean, geological structure analysis for typical permafrost slope, model experiment of thermal-thaw sliding collapse of Qinghai-Tibet milestone K3035, stability analysis method, in-situ monitoring and numerical simulation of temperature field, displacement field. The following conclusions have been put forward.
(1)On the basis of concluding research results of permafrost slope along Qinghai-Tibet highway, permafrost instability patterns are divided into two categories, freezing soil sliding and thawing soil sliding.Gelifluction & thaw slumping are concluded in thawing soil sliding. Creep sliding & collapse sliding are concluded in freezing soil sliding, It is precondition that soils are saturated by water and thick ice takes as sliding face for gelifluction and thaw slumping.
(2)Factors caused by thaw slumping are as follows, a. thick layer ground ice and shallow embedding depth, b. thick layer ground ice uncovered
(3)Thaw slumping in milestone K3035 of Qinghai-Tibet highway is not sliding at a time. At present, movement of slope scale begin with melting, crack, collapse. During positive temperature every year, sliding collapse begin to develop with depth in 2m,between July and October, in the period of 3 years between 2000 and 2002 ,it is estimated velocity of thaw slumping is 5m/a.
(4)Water conditions are important for any slop sliding. It show that during thaw season, for collapse sliding of milestone K3035 of Qinghai-Tibet highway , water content of soil above 2m in depth is more than 30%, it is higher than liquid limit value situated on thick layer ice. Water content is more than 19% for collapse sliding of milestone K3057 of Qinghai-Tibet highway. At the same time, water content of no sliding soil is less than 18%, water content of sliding soil is obviously higher than that of no sliding soil. Water Congregation at ice surface is mainly caused by ground ice melting. It is melting water that soil strength is obviously decrease, and sliding along thick ice surface. 〢ccording to the first law of simulation theory, simulation index and criterion for permafrost slope model have be concluded and built by applying integration comparison method. It is concluded that taking in-situ soil as similar material, six similar parameters will
be decreased two parameters, this is ct and cl .Time scale factor between prototype and
scale model is determined by geometry scale factor, this is ,ct =cl2 . Only one similar factor
is self-variable, another is dependency variable. The following conclusions have been concluded, deformation field is controlled by temperature field. During a series of temperature cycles, there are critical temperature value successive displacement. Frost heaving and thaw settlement for slope soil surface are calculated and analyzed by applying self-design model triangle frame. It is concluded that frost heaving is much more than thaw settlement, it is in accordance with observation results along Qinghai-Tibet highway.
(6)Differential equation of heat transfer with phase change is derived in this paper. Simulation of temperature displacement field of permafrost slope at milestone K3035 of Qinghai-Tibet highway have been carried out by Marc finite element software. Displacement of 5 years is simulated behind sliding step. It is drawn out that displacemen
本文从青藏高原多年冻土区冻土斜坡失稳类型分析、典型冻土斜坡地质结构分析入手，通过对青藏公路K3035里程典型斜坡热融滑塌模型试验研究、稳定性分析、现场监测、温度场和变形场数值分析，开展了青藏高原多年冻土区斜坡稳定性研究，取得了如下成果：
① 在全面总结青藏公路沿线有关斜坡研究成果的基础上，通过调查对青藏高原多年冻土区冻土斜坡失稳类型进行了划分，根据其形成原因和冻土赋存环境，正融土滑坡可以划分为两种类型：融冻泥流型和热融滑塌型；正冻土滑坡可划分为两种类型，即蠕变型滑坡和崩塌型滑坡。
② 在青藏高原多年冻土区斜坡地带活动层土体沿确定的滑动面(冻融界面、厚层地下冰面)在融化季节滑动形成融冻泥流和热融滑塌两类典型的滑坡灾害。热融滑塌的失稳因素可以归结为：a．厚层地下冰的存在且埋藏深度较浅易受地表气温变化的影响，形成冰土分界面相对隔水且抗剪强度极低的滑坡条件，这是热融滑塌形成的内在因素；b．不合理取土或自然营力侵蚀引起的厚层地下冰的暴露且使坡体形成临空面而使其失去支承力，这是热融滑塌形成的外在因素(诱发因素)。
③ 青藏公路K3035热融滑塌的显著特点是滑塌体不是整体一次性滑动，目前的滑坡范围扩展是从后缘原坡体融化、开裂、坍塌，塌落的土体沿滑面滑动。扩展发育时期为每年气温正值期，且主要集中在0～2m深度土体范围内，处于融化状态的7-10月。综合2000～2002年3年内滑坡范围扩展情况，估算现今滑坡范围向后扩展的速度为5m／a。
④ 任何滑坡的产生，水分条件的变化是其重要的影响因素。热融滑塌更是如此，现场监测结果表明：融化季节在青藏公路K3035滑体2m以上范围的土体含水量高于冰层之上的粉质粘土液限值；在青藏公路K3057处大于同期对未滑动活动层土体的含水量的测量结果，即滑体水分显著高于原土体。冰面水分的聚集主要是由地下冰的融化所引起，正是由于大量融化水的存在，使得土体强度显著降低，并沿地下冰面产生了滑动。
⑤ 根据相似理论第一定律，应用积分类比法推导并建立了冻土斜坡模型试验的相似指标和相似判据，得出以原状土作为模型介质时，6个相似常数可以减少为2个相似常数，即c_τ和c_ι。模型与原型的时间比例尺是由几何比例尺决定的，即c_τ=c_ι~2，仅有一个相似参数为自变量，另一个为因变量。应用自行设计的变形测量装置对斜坡土体表面冻胀和融沉进行了计算和分析，得出冻土斜坡的融沉变形远远大于冻胀变形的结论，与青藏公路沿线多年冻土地区公路路基变形观测结果一致。通过模型试验，得出变形场受温度场控制的结论，并在多个冻融循环中，存在着一个斜坡后继变形的“温度门限”，控制着斜坡变形
的进一步发展。这些成果填补了多年冻土区滑坡研究的空白，为青藏高原多年
冻土区道路工程建设中滑坡整治提供了依据。
⑥本文推导了相变热传导问题微分控制方程的一般形式，应用有限元法对
青藏公路K3035斜坡地温场和变形场进行了模拟，建立了冻土斜坡温度场和变
形场数值分析模型。模拟了滑坎后缘拉张区5年内位移变化情况，得出冻土滑
坡体每一年中位移由两个阶段即稳定阶段和非稳定阶段组成的结论，稳定阶段
对应台阶面，非稳定阶段对应台阶壁;每一个台阶面均向右微微“上翘”，可
以得出土体在冻结过程中还发生上坡位移，而不单单仅存在下坡位移，这一计
算结果与风火山冻土蠕变监测结果取得了很好的一致。从模拟坡面内不同深度
对应的下坡位移可以得出冻土斜坡的位移主要发生在浅部，深部基本不发生位
移的结论，与青藏公路K3035斜坡现场测斜监测资料一致。采用基于灰色系统
理论的GM(1，1)模型，对滑塌距离动态变化进行预测，计算滑塌距离为5.76m，
误差为12.94%，说明所建立的GM(1，1)模型可以作为青藏公路K3035滑塌距离
年内月平均变化的预测模型。
⑦应用有效应力原理推导了不同渗流条件下无限斜坡稳定性分析的计算
公式，对各种不同表达形式的安全系数进行了讨论和对比，得出在流线平行于
斜坡条件下，两类特殊条件(干土坡、饱水土坡)下的分析结果与相同情况下
安全系数统一式结果相一致的结论;当已=O及0(m(l时，无限斜坡安全系数的
通式为尸=[(1一m)+m
上]竺旦竺
ysat tana
流线为相互平行的水平线、流线与坡向一致两种
情况分别为孔隙压力和安全系数的两种特例。流线与水平线的夹角为。‘刀‘a，
二，.、二二/。、__刀_____:__:____翔。二。，，、。卜、、、。
孔隙水压力u(刀)=mzy*一二mzr，sin，a，刀一。和刀=a分别对应上述两种极
口
值情况(流线为相互平行的水平线、流线与坡向一致)。
⑧绘制出了不同坡度、不同含水条
Permafrost area of Qinghai-Tibet Platean is 1. 5x106 km2,70 percent of permafrost area in China. On the basis of global warming, there is a degradation tendency for permafrost. The degradation tendency of permafrost will become worse caused by all kinds of human engineering activities. Geological hazards from freeze-thaw action will not only threat linear engineering construction for example, highway, railway, bridge, oil pipeline, etc et, but result in deterioration for ecological environment of cold region to lose recovery ability for permafrost. Thermal-thaw sliding collapse is a typical example.
The following researches have been carried out in this paper, including slope instability patterns analysis in permafrost region of Qinghai-Tibet Platean, geological structure analysis for typical permafrost slope, model experiment of thermal-thaw sliding collapse of Qinghai-Tibet milestone K3035, stability analysis method, in-situ monitoring and numerical simulation of temperature field, displacement field. The following conclusions have been put forward.
(1)On the basis of concluding research results of permafrost slope along Qinghai-Tibet highway, permafrost instability patterns are divided into two categories, freezing soil sliding and thawing soil sliding.Gelifluction & thaw slumping are concluded in thawing soil sliding. Creep sliding & collapse sliding are concluded in freezing soil sliding, It is precondition that soils are saturated by water and thick ice takes as sliding face for gelifluction and thaw slumping.
(2)Factors caused by thaw slumping are as follows, a. thick layer ground ice and shallow embedding depth, b. thick layer ground ice uncovered
(3)Thaw slumping in milestone K3035 of Qinghai-Tibet highway is not sliding at a time. At present, movement of slope scale begin with melting, crack, collapse. During positive temperature every year, sliding collapse begin to develop with depth in 2m,between July and October, in the period of 3 years between 2000 and 2002 ,it is estimated velocity of thaw slumping is 5m/a.
(4)Water conditions are important for any slop sliding. It show that during thaw season, for collapse sliding of milestone K3035 of Qinghai-Tibet highway , water content of soil above 2m in depth is more than 30%, it is higher than liquid limit value situated on thick layer ice. Water content is more than 19% for collapse sliding of milestone K3057 of Qinghai-Tibet highway. At the same time, water content of no sliding soil is less than 18%, water content of sliding soil is obviously higher than that of no sliding soil. Water Congregation at ice surface is mainly caused by ground ice melting. It is melting water that soil strength is obviously decrease, and sliding along thick ice surface. 〢ccording to the first law of simulation theory, simulation index and criterion for permafrost slope model have be concluded and built by applying integration comparison method. It is concluded that taking in-situ soil as similar material, six similar parameters will
be decreased two parameters, this is ct and cl .Time scale factor between prototype and
scale model is determined by geometry scale factor, this is ,ct =cl2 . Only one similar factor
is self-variable, another is dependency variable. The following conclusions have been concluded, deformation field is controlled by temperature field. During a series of temperature cycles, there are critical temperature value successive displacement. Frost heaving and thaw settlement for slope soil surface are calculated and analyzed by applying self-design model triangle frame. It is concluded that frost heaving is much more than thaw settlement, it is in accordance with observation results along Qinghai-Tibet highway.
(6)Differential equation of heat transfer with phase change is derived in this paper. Simulation of temperature displacement field of permafrost slope at milestone K3035 of Qinghai-Tibet highway have been carried out by Marc finite element software. Displacement of 5 years is simulated behind sliding step. It is drawn out that displacemen
引文
[1]周幼吾、程国栋等，冻土环境研究现状及任务，第四届全国冰川冻土学术会议论文选集(冻土学)科学出版社，1990，171～176
[2]徐学祖、丁德文等，冻土改造和利用的研究现状和展望，第四届全国冰川冻土学术会议论文选集(冻土学)科学出版社，1990，191～199
[3]徐学祖、王家澄，中国冻土分布及其地带性规律的初步探讨 第二届全国冻土学术会议论文选集 科学出版社，1982，3～12
[4]吴紫汪、赖远明等著，寒区隧道工程 北京：海洋出版社 2003年，1～
[5]周幼吾、杜榕桓 青藏高原冻土初步考察 科学通报，1963．2
[6]周幼吾，多年冻土 中国自然地理(地貌)，北京：科学出版社 1980
[7]周幼吾、郭东信 我国多年冻土的主要特征 冰川冻土 1982，4 (1)：1～19
[8]王绍令、米海珍 青藏公路铺筑沥青路面后路基下多年冻土的变化《冰川冻土》，1993，15(4)：566～573
[9]秦大河(总主编)，丁一江(主编) 中国西部环境演变评估(第二卷)中国西部环境变化的预测 2002，北京，科学出版社 174～178
[10]朱林楠、吴紫汪等 青藏高原东部的冻土退化《冰川冻土》，1995，17 (2)：120～124
[11]吴青柏 人类工程活动下冻土环境变化和工程适应性研究 中国科学院博士学位论文 2000年11月62～
[12]王绍令 青藏公路风火山地区的热融滑塌 《冰川冻土》，1990，12 (1)：63～68
[13]中国科学院兰州冰川冻土所 吴紫汪等《冻土路基工程》兰州大学出版社 1988 2～
[14]中国科学院兰州冰川冻土沙漠研究所 著《冻土》科学出版社 1975．6
[15]周幼吾，郭东信等．中国冻土[M]．北京：科学出版社，2000．403～404．
[16]Tarr R S. Rapidity of Weathering and Stream Erosion in the Arctic Latitudes[J]. American Geologist, 1897, 19: 131～136.
[17][苏]叶米里扬诺娃娃著 滑坡作用的基本规律 铁道部科学研究院西北所滑坡室译，重庆出版社 1986
[18]崔之久 1980初探青藏高原特殊的冰缘现象 科学通报 第7期
[19]王文宝 多年冻土地区铁路路基边坡滑塌及其整治 第三届全国冻土学术会议论文选集 中国科学院兰州冰川冻土研究所 科学出版社 1989
[20]张廷军、周幼吾 1989第三届全国冻土学术会议论文选集 48～53
[21]程国栋 青藏铁路工程与多年冻土相互作用及环境效应 中国科学院院刊 2002(1)21～25
[22]张长庆，朱林南等，中国南水北调西线工程地区的冻土与工程问题，1993，冰川冻土，Vol．15(1)：90～95
[23] McRoberts E .C,& Morgenstern N.R The stability of thawing slopes Can.Geotech.J. 1974,11:447～469
[24] Clark M J. Advance in Periglacial Geomorphology[M]. New York: John Wiley & Sons Ltd, 1988, 325～359.
[25] Niu Fujun, Cheng Guodong, Xie Qun. Study on Instability of Slopes in Permafrost Regions of Qinghai-Tibet High Plateau[C]. Proceedings of the 5th International Symposium on Permafrost Engineering. Yakutsk: Permafrost Inst. SB RAS Press, 2002.
192～197．
[26]殷坤龙等 国际滑坡研究的新进展 水文地质工程地质 2000，27(5)：1～3
[27]郭东信等 青藏公路风火山垭口盆地融冻泥流阶地初步研究 冰川冻土 1993，15 (1)：58～62
[28]王绍令 青藏公路风火山地区的热融滑塌 冰川冻土 1990 vol．12 (1)：63～70
[29]Wang B, French H M. In Situ Creep of Frozen Soil[J]. Tibet Pleateau, China. Can. Geotech J., 1995, 23: 545～552.
[30]Zhu Cheng. Zhang Jian-xin, CHENG Peng. Rock Glaciers in the Central Tianshan Mountain, China[J]. Permafrost and Periglacial Processes, 1996, 7:69～78.
[31]Weeks A G The stability of natural slopes in south-east England as affected by periglacial activity, Quarterly Journal of Engineering Geology, 1969 2:49～63
[32]Hutchinson,J.N. Periglacial solifluction:an approximate mechanism for clay soil,Geotechnique, 1974,Vo 1.24:438～443
[33]郑启浦 东北多年冻土地区不良地质现象对铁路工程的危害及防治 中国地质灾害与防治学报 1993，4(1)：56～61
[34]李永春、赵康年 青藏公路K2023反季滑坡实例 地质灾害与环境保护 2002，13(4)：33～35
[35]张建永 滑坡研究现状综述 中国岩溶 1999，18 (3)：280～286
[36]Niu F J, Cheng G D, Xie Q. Study on Instability of Slopes in Permafrost Regions of Qinghai-Tibet High Plateau[C]. Proceedings of the 5th International Symposium on Permafrost Engineering. Yakutsk: Permafrost Inst. SB RAS Press, 2002. 192～197.
[37]Hakata M, Yoshimasa K, Kazuomi H. Soil Moisture Conditions during Thawingon a Slope in the Daisetsu Mountains, Hokkaido, Japan[J]. Permafrost and Periglacial Processes, 2001, 12: 211～218.
[38]丁德文、罗学波 冻土热工模型的理论基础 科学通报 1979(8)：360～364
[39]李之光编著，热力设备实验研究基础，北京：国防工业出版社 1971：35～39
[40]李之光编著 相似与模化：(理论及应用)北京：国防工业出版社1982：54～94
[41]林韵梅编著．实验岩石力学—模拟研究[M]．北京：煤炭工业出版社，1984：3～17
[42]魏先祥、赖远明 相似方法的原理及应用 兰州：兰州大学出版社 2001：1～38]
[43]邱国庆主 编冻土学辞典 214～科学出版社 1994年5月
[44]青藏高原多年冻土地区公路路基变形 刘永智等 《冰川冻土》2002(1)：10～15
[45]李东庆等，边坡渗流对冻土地区路基稳定性的影响分析，兰州火学学报(自然科学版)，2000，36(3)：175～179
[46]薛禹群 主编 地下水动力学原理 地质出版社，1986年12月136～206
[47]柴军瑞等，岩体渗流场与温度场耦合的连续介质模型 地下水 1997，(2)
[48 朱伯芳等 水工混凝土结构的温度应力与温度控制 水利电力出版社，1976
[49]朱伯芳 有限单元法原理及应用 水利电力出版社 1978．8：251～254
[50]安维东主编 冻土的温度水分应力及其相互作用 兰州大学出版社，1989：212～217
[51]李东庆 吴紫汪等 多年冻土退化的数值模拟计算 第五届全国冰川冻土学大会论文集(下册)，甘肃文化出版社 1996：1082～1086
[52]李东庆 吴紫汪等 花石峡连续多年冻土区湿润性地段路基稳定性模拟分析 岩土工程学报 1999，21 (1)：17～20
[53]Bonacina, Comini,G, et al. Numerical solution of phase change problems, Int. J. Heat
Transfer,1973,16:1825～1832
[54]数学手册编写组 数学手册 1979．5 高等教育出版社：892～894
[55]郭宽良等编著 计算传热学 中国科学技术大学出版社 1988．12：59～69
[56]朱林楠 高原冻土区不同下垫面的附面层研究 冰川冻土 1998，10(1)：8～14
[57]吴紫汪 程国栋等 冻土路基工程 兰州大学出版社 1998年9月，48～60
[58]王补宣 工程传热传质学 (上册)科学出版社 1982年8月，185～
[59]李述训 程国栋 郭东信 气候持续变暖条件下青藏高原多年冻土变化趋势数值模拟中国科学 1996(4)：342～347
[60]刘慈群 郭白奇 二维路堤有相变的热传导的数值计算方法 计算机科学与应用 1981
[61]郭白奇 刘慈群 有相变热的热传导问题 应用数学学报 1977 (4)：51～57
[62]令锋 吴紫汪 渗流对多年冻土区路基温度场影响的数值模拟 冰川冻土 1999(2)：115～119
[63]王贤能 黄润秋 深埋长隧洞温度场的评价预测 水文地质工程地质 1996 (6)：6～10
[64]陈尚桥 黄润秋 深埋隧洞地温场的数值模拟研究 地质灾害与环境保护 1995，6(2)：30～36
[65]徐学祖 王家澄、张立新 著 冻土物理学，科学出版社，2001：39～98
[66]徐学祖 陶兆祥 傅素兰 典型融冻土的热学性质 中国科学院兰州冰川冻土研究所集刊(第2号)科学出版社，1981，55～71
[67]徐学祖 邓友生 冻土中水分迁移的实验研究 科学出版社 1991，71～94
[68]吴紫汪等 冻结凿井冻土壁的工程性质 兰州大学出版社 1988年9月，1～80
[69]邱国庆等主编 冻土学辞典 甘肃科学技术出版社 1994年5月，23，41～
[70]李述训 程国栋 编著 冻融土中的水热输运问题 兰州大学出版社 1995年7月，25～26
[71]王竹溪 热力学简程 人民教育出版社，1979
[72]陈火红 编著 Marc有限元实例分析教程 机械工业出版社，2002年4月，341～374
[73]王书法 朱维申 节理岩体热应力问题的非线性连续变形分析方法 岩石力学与工程学报 2000，19(5)：643～646
[74]邓聚龙著 灰色控制系统 华中工学院出版社，1985：318～336
[75]彭国伦编著Fortran 95程序设计 中国电力出版社 2002年9月
[76]徐明编著Fortran powerstation version 4.0基础教程 清华大学出版社 2000年6月
[77]徐士良编 计算机常用算法 清华大学出版社 1989年11月
[78]谭浩强 田淑清 FORTRAN语言—FORTRAN 77结构化程序设计 清华大学出版社 1990年3月
[79]张咸恭 地下水对工程和环境的作用 工程地质学报 1993，1(1)：1～6
[80]王旭升等．斜坡地下水渗透力计算与稳定性分析[J]．水文地质工程地质，2003，30 (2)：41～45
[81]李媛等，孟家山黄土—红层接触面滑坡破坏机理研究[J]．水文地质工程地质，2001，28 (1)：52～54
[82]牛富俊等．青藏高原多年冻土区斜坡类型及典型斜坡稳定性研究[J]．冰川冻土，2002，24 (5)：608～613
[83]陈祖煜著 土质边坡稳定分析—原理 方法 程序 中国水利水电出版社 2003年1月，15～21
[84]靳德武 牛富俊等 不同渗流条件下无限斜坡稳定性分析探讨 地质灾害与环境保护
2003，14(4)：63～67
[85]胡中雄编著，土力学与环境土工学 上海：同济大学出版社：273～278
[86]Braja M.Das Principles of Geotechnical Engineering USA: Maple-Vall Book Manufacturing Group :431～437
[87]靳德武 牛富俊等 不同渗透流条件下无限斜坡稳定性分析方法探讨 地质灾害与环境保护 2003，14(4)：63～67
[88]程国栋 何平 多年冻土地区线性工程建设 冰川冻土 2001，23(3)：213～217
[89]郭兴民 朱林楠 小型冻土模型试验装置的研制冻土 工程国家重点实验室年报(1993)：293～295
[90]程圣国等，土质滑坡相似材料试验设计理论及评价方法研究，水力发电，2002(4)：21～22
[91]何建、马景骏，潮湿路基温度场、湿度场耦合作用计算模型 哈尔滨工程大学学报，2001，22(3)：59～62
[92]刘鸿绪 对冻土学几个基本术语的商榷 第五届全国冰川冻土学大会论文集(上册)：541～544
[93]马巍 程国栋等 多年冻土地区主动冷却地基方法研究 冰川冻土 2002，24(5) 579～587
[94]C.BONCINA & G.COMINI On the solution of the nonlinear heat conduction equations by numerical methods Int.J.Heat Mass Transfer, 1973,Vol. 16:581～589
[95]雷志栋等 地下水浅埋条件下越冬期土壤水热迁移的数值模拟 冰川冻土，1998，20(1)：51～54
[96]何平 程国栋等 饱和正冻土中水、热力耦合模型 冰川冻土 2000，22(2)：135～138
[97]杨海蓉 多年冻土地区防治路基融化下沉及提高其稳定性的措施 冰川冻土，1985，7(1)：83～88
[98]冯文杰等 不同边界条件对多年冻土上限影响的模型试验研究 冰川冻土 2001，23(4)：353～359
[99]何平 程国栋 土体冻结过程中的热质迁移研究进展 冰川冻土，2001，23(1)：92～98
[100]夏才初 李永盛编著 地下工程测试理论与监测技术 同济大学出版社 1999年8月，94～118
[101]夏才初 潘国荣等编著 土木工程监测技术 中国建筑工业出版社 2001年7月 186～213
[102]刘玉海 陈志新 倪万魁 赵法锁等 大同城市地质研究 西安地图出版社 1995
[103]杨英华 主编 土力学 地质出版社 1987年5月157～178
[104]胡广韬 杨文元 主编 工程地质学 地质出版社 158～172
[105]杨林德等编著 岩土工程问题的反演理论与工程实践 科学出版社 1996年4月，1～63
[106]刘祖典编著 黄土力学与工程 陕西科学技术出版社 1997年4月 160～248
[107]尚岳全 黄润秋编著 工程地质研究中的数值分析方法 成都科技大学出版社，1991年8月，96～134
[108]黄润秋 张倬元 王士天著 黄河拉西瓦水电站高边坡稳定性的系统工程地质研究 成都科技大学出版社 1991年3月 52～67
[109]孙广忠著 工程地质与地质工程 地震出版社 1993年6月，155～182
[110]赵彭年 编著 松散介质力学 地震出版社 1995年6月，101～108
[111]崔之久 青藏高原冰缘地貌的基本特征 中国科学 1981(6)：724～733
[112]王绍令 青藏公路沿线现代冰缘作用类型及其特征 冰川冻土，1986，8(2)
[113]何平 朱元林等 冻土的变形性能与泊松比 地下空间 1999，19(5)：504～507
[114]王正中 沙际德等 正交各向异性冻土与建筑物相互作用的非线性有限元分析 土木工程学报 1999，32(3)：55～60
[115]陈志新 黄土高边坡稳定性影响因素分析 全国岩土与工程学术人会文集 人民交通出版社，601～605
[116]倪万魁 黄土高边坡层状结构与稳定性关系研究 工程地质学报 2002，10：279～284
[117]陈志新 倪万魁 著 延安滑坡及其灾害防治 西安地图出版社，2002年2月
[2]徐学祖、丁德文等，冻土改造和利用的研究现状和展望，第四届全国冰川冻土学术会议论文选集(冻土学)科学出版社，1990，191～199
[3]徐学祖、王家澄，中国冻土分布及其地带性规律的初步探讨 第二届全国冻土学术会议论文选集 科学出版社，1982，3～12
[4]吴紫汪、赖远明等著，寒区隧道工程 北京：海洋出版社 2003年，1～
[5]周幼吾、杜榕桓 青藏高原冻土初步考察 科学通报，1963．2
[6]周幼吾，多年冻土 中国自然地理(地貌)，北京：科学出版社 1980
[7]周幼吾、郭东信 我国多年冻土的主要特征 冰川冻土 1982，4 (1)：1～19
[8]王绍令、米海珍 青藏公路铺筑沥青路面后路基下多年冻土的变化《冰川冻土》，1993，15(4)：566～573
[9]秦大河(总主编)，丁一江(主编) 中国西部环境演变评估(第二卷)中国西部环境变化的预测 2002，北京，科学出版社 174～178
[10]朱林楠、吴紫汪等 青藏高原东部的冻土退化《冰川冻土》，1995，17 (2)：120～124
[11]吴青柏 人类工程活动下冻土环境变化和工程适应性研究 中国科学院博士学位论文 2000年11月62～
[12]王绍令 青藏公路风火山地区的热融滑塌 《冰川冻土》，1990，12 (1)：63～68
[13]中国科学院兰州冰川冻土所 吴紫汪等《冻土路基工程》兰州大学出版社 1988 2～
[14]中国科学院兰州冰川冻土沙漠研究所 著《冻土》科学出版社 1975．6
[15]周幼吾，郭东信等．中国冻土[M]．北京：科学出版社，2000．403～404．
[16]Tarr R S. Rapidity of Weathering and Stream Erosion in the Arctic Latitudes[J]. American Geologist, 1897, 19: 131～136.
[17][苏]叶米里扬诺娃娃著 滑坡作用的基本规律 铁道部科学研究院西北所滑坡室译，重庆出版社 1986
[18]崔之久 1980初探青藏高原特殊的冰缘现象 科学通报 第7期
[19]王文宝 多年冻土地区铁路路基边坡滑塌及其整治 第三届全国冻土学术会议论文选集 中国科学院兰州冰川冻土研究所 科学出版社 1989
[20]张廷军、周幼吾 1989第三届全国冻土学术会议论文选集 48～53
[21]程国栋 青藏铁路工程与多年冻土相互作用及环境效应 中国科学院院刊 2002(1)21～25
[22]张长庆，朱林南等，中国南水北调西线工程地区的冻土与工程问题，1993，冰川冻土，Vol．15(1)：90～95
[23] McRoberts E .C,& Morgenstern N.R The stability of thawing slopes Can.Geotech.J. 1974,11:447～469
[24] Clark M J. Advance in Periglacial Geomorphology[M]. New York: John Wiley & Sons Ltd, 1988, 325～359.
[25] Niu Fujun, Cheng Guodong, Xie Qun. Study on Instability of Slopes in Permafrost Regions of Qinghai-Tibet High Plateau[C]. Proceedings of the 5th International Symposium on Permafrost Engineering. Yakutsk: Permafrost Inst. SB RAS Press, 2002.
192～197．
[26]殷坤龙等 国际滑坡研究的新进展 水文地质工程地质 2000，27(5)：1～3
[27]郭东信等 青藏公路风火山垭口盆地融冻泥流阶地初步研究 冰川冻土 1993，15 (1)：58～62
[28]王绍令 青藏公路风火山地区的热融滑塌 冰川冻土 1990 vol．12 (1)：63～70
[29]Wang B, French H M. In Situ Creep of Frozen Soil[J]. Tibet Pleateau, China. Can. Geotech J., 1995, 23: 545～552.
[30]Zhu Cheng. Zhang Jian-xin, CHENG Peng. Rock Glaciers in the Central Tianshan Mountain, China[J]. Permafrost and Periglacial Processes, 1996, 7:69～78.
[31]Weeks A G The stability of natural slopes in south-east England as affected by periglacial activity, Quarterly Journal of Engineering Geology, 1969 2:49～63
[32]Hutchinson,J.N. Periglacial solifluction:an approximate mechanism for clay soil,Geotechnique, 1974,Vo 1.24:438～443
[33]郑启浦 东北多年冻土地区不良地质现象对铁路工程的危害及防治 中国地质灾害与防治学报 1993，4(1)：56～61
[34]李永春、赵康年 青藏公路K2023反季滑坡实例 地质灾害与环境保护 2002，13(4)：33～35
[35]张建永 滑坡研究现状综述 中国岩溶 1999，18 (3)：280～286
[36]Niu F J, Cheng G D, Xie Q. Study on Instability of Slopes in Permafrost Regions of Qinghai-Tibet High Plateau[C]. Proceedings of the 5th International Symposium on Permafrost Engineering. Yakutsk: Permafrost Inst. SB RAS Press, 2002. 192～197.
[37]Hakata M, Yoshimasa K, Kazuomi H. Soil Moisture Conditions during Thawingon a Slope in the Daisetsu Mountains, Hokkaido, Japan[J]. Permafrost and Periglacial Processes, 2001, 12: 211～218.
[38]丁德文、罗学波 冻土热工模型的理论基础 科学通报 1979(8)：360～364
[39]李之光编著，热力设备实验研究基础，北京：国防工业出版社 1971：35～39
[40]李之光编著 相似与模化：(理论及应用)北京：国防工业出版社1982：54～94
[41]林韵梅编著．实验岩石力学—模拟研究[M]．北京：煤炭工业出版社，1984：3～17
[42]魏先祥、赖远明 相似方法的原理及应用 兰州：兰州大学出版社 2001：1～38]
[43]邱国庆主 编冻土学辞典 214～科学出版社 1994年5月
[44]青藏高原多年冻土地区公路路基变形 刘永智等 《冰川冻土》2002(1)：10～15
[45]李东庆等，边坡渗流对冻土地区路基稳定性的影响分析，兰州火学学报(自然科学版)，2000，36(3)：175～179
[46]薛禹群 主编 地下水动力学原理 地质出版社，1986年12月136～206
[47]柴军瑞等，岩体渗流场与温度场耦合的连续介质模型 地下水 1997，(2)
[48 朱伯芳等 水工混凝土结构的温度应力与温度控制 水利电力出版社，1976
[49]朱伯芳 有限单元法原理及应用 水利电力出版社 1978．8：251～254
[50]安维东主编 冻土的温度水分应力及其相互作用 兰州大学出版社，1989：212～217
[51]李东庆 吴紫汪等 多年冻土退化的数值模拟计算 第五届全国冰川冻土学大会论文集(下册)，甘肃文化出版社 1996：1082～1086
[52]李东庆 吴紫汪等 花石峡连续多年冻土区湿润性地段路基稳定性模拟分析 岩土工程学报 1999，21 (1)：17～20
[53]Bonacina, Comini,G, et al. Numerical solution of phase change problems, Int. J. Heat
Transfer,1973,16:1825～1832
[54]数学手册编写组 数学手册 1979．5 高等教育出版社：892～894
[55]郭宽良等编著 计算传热学 中国科学技术大学出版社 1988．12：59～69
[56]朱林楠 高原冻土区不同下垫面的附面层研究 冰川冻土 1998，10(1)：8～14
[57]吴紫汪 程国栋等 冻土路基工程 兰州大学出版社 1998年9月，48～60
[58]王补宣 工程传热传质学 (上册)科学出版社 1982年8月，185～
[59]李述训 程国栋 郭东信 气候持续变暖条件下青藏高原多年冻土变化趋势数值模拟中国科学 1996(4)：342～347
[60]刘慈群 郭白奇 二维路堤有相变的热传导的数值计算方法 计算机科学与应用 1981
[61]郭白奇 刘慈群 有相变热的热传导问题 应用数学学报 1977 (4)：51～57
[62]令锋 吴紫汪 渗流对多年冻土区路基温度场影响的数值模拟 冰川冻土 1999(2)：115～119
[63]王贤能 黄润秋 深埋长隧洞温度场的评价预测 水文地质工程地质 1996 (6)：6～10
[64]陈尚桥 黄润秋 深埋隧洞地温场的数值模拟研究 地质灾害与环境保护 1995，6(2)：30～36
[65]徐学祖 王家澄、张立新 著 冻土物理学，科学出版社，2001：39～98
[66]徐学祖 陶兆祥 傅素兰 典型融冻土的热学性质 中国科学院兰州冰川冻土研究所集刊(第2号)科学出版社，1981，55～71
[67]徐学祖 邓友生 冻土中水分迁移的实验研究 科学出版社 1991，71～94
[68]吴紫汪等 冻结凿井冻土壁的工程性质 兰州大学出版社 1988年9月，1～80
[69]邱国庆等主编 冻土学辞典 甘肃科学技术出版社 1994年5月，23，41～
[70]李述训 程国栋 编著 冻融土中的水热输运问题 兰州大学出版社 1995年7月，25～26
[71]王竹溪 热力学简程 人民教育出版社，1979
[72]陈火红 编著 Marc有限元实例分析教程 机械工业出版社，2002年4月，341～374
[73]王书法 朱维申 节理岩体热应力问题的非线性连续变形分析方法 岩石力学与工程学报 2000，19(5)：643～646
[74]邓聚龙著 灰色控制系统 华中工学院出版社，1985：318～336
[75]彭国伦编著Fortran 95程序设计 中国电力出版社 2002年9月
[76]徐明编著Fortran powerstation version 4.0基础教程 清华大学出版社 2000年6月
[77]徐士良编 计算机常用算法 清华大学出版社 1989年11月
[78]谭浩强 田淑清 FORTRAN语言—FORTRAN 77结构化程序设计 清华大学出版社 1990年3月
[79]张咸恭 地下水对工程和环境的作用 工程地质学报 1993，1(1)：1～6
[80]王旭升等．斜坡地下水渗透力计算与稳定性分析[J]．水文地质工程地质，2003，30 (2)：41～45
[81]李媛等，孟家山黄土—红层接触面滑坡破坏机理研究[J]．水文地质工程地质，2001，28 (1)：52～54
[82]牛富俊等．青藏高原多年冻土区斜坡类型及典型斜坡稳定性研究[J]．冰川冻土，2002，24 (5)：608～613
[83]陈祖煜著 土质边坡稳定分析—原理 方法 程序 中国水利水电出版社 2003年1月，15～21
[84]靳德武 牛富俊等 不同渗流条件下无限斜坡稳定性分析探讨 地质灾害与环境保护
2003，14(4)：63～67
[85]胡中雄编著，土力学与环境土工学 上海：同济大学出版社：273～278
[86]Braja M.Das Principles of Geotechnical Engineering USA: Maple-Vall Book Manufacturing Group :431～437
[87]靳德武 牛富俊等 不同渗透流条件下无限斜坡稳定性分析方法探讨 地质灾害与环境保护 2003，14(4)：63～67
[88]程国栋 何平 多年冻土地区线性工程建设 冰川冻土 2001，23(3)：213～217
[89]郭兴民 朱林楠 小型冻土模型试验装置的研制冻土 工程国家重点实验室年报(1993)：293～295
[90]程圣国等，土质滑坡相似材料试验设计理论及评价方法研究，水力发电，2002(4)：21～22
[91]何建、马景骏，潮湿路基温度场、湿度场耦合作用计算模型 哈尔滨工程大学学报，2001，22(3)：59～62
[92]刘鸿绪 对冻土学几个基本术语的商榷 第五届全国冰川冻土学大会论文集(上册)：541～544
[93]马巍 程国栋等 多年冻土地区主动冷却地基方法研究 冰川冻土 2002，24(5) 579～587
[94]C.BONCINA & G.COMINI On the solution of the nonlinear heat conduction equations by numerical methods Int.J.Heat Mass Transfer, 1973,Vol. 16:581～589
[95]雷志栋等 地下水浅埋条件下越冬期土壤水热迁移的数值模拟 冰川冻土，1998，20(1)：51～54
[96]何平 程国栋等 饱和正冻土中水、热力耦合模型 冰川冻土 2000，22(2)：135～138
[97]杨海蓉 多年冻土地区防治路基融化下沉及提高其稳定性的措施 冰川冻土，1985，7(1)：83～88
[98]冯文杰等 不同边界条件对多年冻土上限影响的模型试验研究 冰川冻土 2001，23(4)：353～359
[99]何平 程国栋 土体冻结过程中的热质迁移研究进展 冰川冻土，2001，23(1)：92～98
[100]夏才初 李永盛编著 地下工程测试理论与监测技术 同济大学出版社 1999年8月，94～118
[101]夏才初 潘国荣等编著 土木工程监测技术 中国建筑工业出版社 2001年7月 186～213
[102]刘玉海 陈志新 倪万魁 赵法锁等 大同城市地质研究 西安地图出版社 1995
[103]杨英华 主编 土力学 地质出版社 1987年5月157～178
[104]胡广韬 杨文元 主编 工程地质学 地质出版社 158～172
[105]杨林德等编著 岩土工程问题的反演理论与工程实践 科学出版社 1996年4月，1～63
[106]刘祖典编著 黄土力学与工程 陕西科学技术出版社 1997年4月 160～248
[107]尚岳全 黄润秋编著 工程地质研究中的数值分析方法 成都科技大学出版社，1991年8月，96～134
[108]黄润秋 张倬元 王士天著 黄河拉西瓦水电站高边坡稳定性的系统工程地质研究 成都科技大学出版社 1991年3月 52～67
[109]孙广忠著 工程地质与地质工程 地震出版社 1993年6月，155～182
[110]赵彭年 编著 松散介质力学 地震出版社 1995年6月，101～108
[111]崔之久 青藏高原冰缘地貌的基本特征 中国科学 1981(6)：724～733
[112]王绍令 青藏公路沿线现代冰缘作用类型及其特征 冰川冻土，1986，8(2)
[113]何平 朱元林等 冻土的变形性能与泊松比 地下空间 1999，19(5)：504～507
[114]王正中 沙际德等 正交各向异性冻土与建筑物相互作用的非线性有限元分析 土木工程学报 1999，32(3)：55～60
[115]陈志新 黄土高边坡稳定性影响因素分析 全国岩土与工程学术人会文集 人民交通出版社，601～605
[116]倪万魁 黄土高边坡层状结构与稳定性关系研究 工程地质学报 2002，10：279～284
[117]陈志新 倪万魁 著 延安滑坡及其灾害防治 西安地图出版社，2002年2月