（岩土工程专业论文）冻土路基温度场的试验研究及仿真分析.pdf

摘要 摘要 本文以青藏铁路清水河段为背景，根据该试验段高寒、高海拔的自然环境 以及多年冻土的不良地质条件，对多年冻土温度场受路基施工的扰动规律进行 了分析。 在综合考虑了太阳辐射、气温、风速、风向及蒸发等真实气象条件下，提出 将诸多气象因素叠加为第1 类边界条件的建模方法，建立了路基非稳态温度场 有限元分析模型：应用a n s y s 9 0 大型有限元软件，对不同断面、不同地基高度、 不同季节的冻土路基温度场进行有限元计算。 通过有限元计算结果与实测结果的比较，验证了模型的可行性。分析了冻 土路基温度分布对其沉降变形的影响规律。研究表明： 1 ) 青藏高原多年冻土区的冻土环境十分脆弱，铁路路基的修筑，改变了原 来天然地表与外界的热交换关系。填土热量的下传，使冻土路基上限受到扰动， 在路基地基面下一定深度范围内形成了融化夹层，它是控制着路基自身稳定后 沉降变形的主要因素，也使得多年冻土在垂向上呈现出不衔接状的特点。 2 ) 由于路堤走向和冻土区气候等原因，造成了路堤阴阳两面热量不平衡， 阳面一侧吸收的热量远大于阴面- - n ，阳面一侧多年冻土的融化深度大于阴面 一侧，使得融化夹层的中心向阳面一侧偏移，进而造成路基两侧的不均匀沉降。 3 ) 由于路堤填土热量随时间的增长慢慢消散，直至最后消失，预测了远期 路基是稳定的，并计算出了路基阴阳两面温度的差异。 保护冻土环境是保证路基稳定的前提，由此提出了一些青藏铁路基施工期 和运营期采取的保护冻土环境措施。 关键词：冻土路基，温度场，冻土上限，变形，冻土保护 a b s t r a c t a b s t r a c t b a s e do nt h eq i n g s h u i h ep r o j e c t ，ap a r to fq i n g h a i t i b e tr a i l w a y , a c c o r d i n gt o t h eb a dg e o l o g i c a lc o n d i t i o n ( h i g ha l t i t u d e ，h i 曲c o l dr e g i o n ，p e r m a f r o s 0 ，t h i sp a p e r a n a l y z e st h ei n f l u e n c er o l e so fe m b a n k m e n tc o n s t r u c t i o no nt e m p e r a t u r ed i s t r i b u t i o n o ft h ep e r m a f r o s t i nt h i sp a p e r ，m a n ya c t u a lm e t e o r o l o g i c a lf a c t o r ss u c ha ss o l a rr a d i a t i o n ，a i r t e m p e r a t u r e ，w i n dv e l o c i t y , w i n dd i r e c t i o na n d s o i lm o i s t u r ee v a p o r a t i o nw e r e c o n s i d e r e d ，a n da n e wm o d e l i n gm e t h o dt oa d dt h e ma st h ef i r s tb o u n d a r yc o n d i t i o n s w a sp u tf o r w a r d a c c o r d i n g l y , f e mm o d a l so fe m b a n k m e n tt h e r m a lr e g l m ew e r es e t u p t h et r a n s i e n tt h e r m a lr e g i m ei ne m b a n k m e n tw e r en u m e r i c a l l ya n a l y z e db yf e m s o f t w a r ea n s y s 9 0w i t hd i f f e r e n tc r o s s s e c t i o n ，d i f f e r e n te m b a n k m e n td e p t h ， d i f f e r e n ts e a s o n t h ef e a s i b i l i t yo ft h em o d e lw a sv a l i d a t e db yc o m p a r i n gt h ef i n i t ee l e m e n tr e s u l t a n dt h em e a s u r e dr e s u l t m o r e o v e rt h ei n f l u e n c eo ft h et e m p e r a t u r ed i s t r i b u t i o no ft h e p e r m a f r o s te m b a n k m e n to ni t sd e f o r m a t i o nw a sa n a l y z e d r e s e a r c hi n d i c a t e ： 1 ) t h ef r o z e ns o i le n v i r o n m e n th a sb e e nv e r yf r a g i l ei nq i n g h a i t i b e tp l a t e a u p e r m a f r o s ta r e a s ，e m b a n k m e n tc o n s t r u c t i o no np e r m a f r o s tw o u l dg i v es u b s t a n t i a l d i s t u r b a n c ef o rt h eh e a tt r a n s f e rb a l a n c eb e t w e e nt h eg r o u n ds u r f a c ea n de n v i r o n m e n t t h eh e a to fe m b a n k m e n tp a s s i n gd o w nt ot h eg r o u n dm a k e st h eu p p e rl i m i to ft h e p e r m a f r o s tr o a df a l l ，a n dl i t t l eb yl i a l eat h a w i n gi n t e r l a y e ri sf o r m e d ，w h i c hi st h e m a i nc o n t r o l l i n g f a c t o r _ o fr o a d b e ds t a b i l i z a t i o na n ds e t t l e m e n t ，a n da l s oc a u s e sa n u n l i n k e dc o n d i t i o no fp e r m a f r o s ti nv e r t i c a l 2 ) b e c a u s eo fd i r e c t i o no fs u b - g r a d ea n dp e r m a f r o s ta r e a sc l i m a t e t h eh e a to f e m b a n k m e n t st w os i d e si si m b a l a n c e ，a n dt h et e m p e r a t u r ep o s i t i v es i d ei sf a rg r e a t e r t h a n n e g a t i v e s ，t h ep o s i t i v et h a wd e p t hi sg r e a t e rt h a nt h en e g a t i v e s ，t h e nt h e t h a w i n gc o r es h i f t st o w a r dt h ep o s i t i v es i d e e v e n t u a l l yt h et h a w i n gs e t t l e m e n ti nt w o s i d e so fe m b a n k m e n ti sn o te v e n “3 ) t h eh e a to fe m b a n k m e n ti sg r a d u a l l yd i s a p p e a r e dw i t ht i m ei n c r e a s i n g ，a n d a b s t r a c t t h ee m b a n k m e n tw i l lb e c o m em o r ea n dm o r es t a b i l i t yi nf u t u r e ，a n dt h ed i f f e r e n to f t h et w os i d e s t e m p e r a t u r ei sc o m p u t e d p r o t e c t i n gt h ee n v i r o n m e n to ff r o z e ns o i li st h ep r e c o n d i t i o nf o re m b a n k m e n t s t a b i l i t y , t h e r e b ys o m em e a s u r e sa r ep r o v i d e dt op r o t e c tt h ef r o z e ns o i le n v i r o n m e n t d u r i n ge n g i n e e rc o n s t r u c t i o na n do p e r a t i o n k e yw o r d s ：f r o z e ns o i le m b a n k m e n t ，t e m p e r a t u r ef i e l d ，p e r m a f r o s tt a b l e ，d e f o r m a t i o n ， f r o z e ns o i lp r o t e c t i n g 独创性声明 本人声明所呈交的论文是我个人在导师指导下进行的研究工作 及取得的研究成果。尽我所知，除了文中特别加以标注和致谢的地方 外，论文中不包含其他人已经发表或撰写过的研究成果，也不包含为 获得石家庄铁道学院或其它教育机构的学位或证书而使用过的材料。 与我一同工作的同志对本研究所做的任何贡献均已在论文中作了明 确的说明并表示了谢意。 关于论文使用授权的说明 2 岬、岁、j 6 本人完全了解石家庄铁道学院有关保留、使用学位论文的规定， 即：学院有权保留送交论文的复印件，允许论文被查阅和借阅；学校 可以公布论文的全部或部分内容，可以采用影印、缩印或其他复制手 段保存论文。 ( 保密的论文在解密后应遵守此规定) 签名：蛀导师签名： 幽日期： 第一章绪论 1 1 国内外研究现状 1 1 1 冻土研究概况 第一章绪论 冻土是一种由固体矿物颗粒、粘塑性冰包裹体和液相水( 未冻水和强结合水) 和气态包裹体( 水汽和空气) 组成的特殊土体，是广泛分布在地球表层的一种 低温地质体，其特殊性主要表现在它的工程性质与温度密切相关。常规类土体 的性质主要受其颗粒的矿物和机械成分、密度和含水量的控制，只要这些因素 确定，土的性质就基本稳定，因此多表现为静态特性。冻土与普通土不同，冻 土的特性除与上述因素有关外，还受含冰量的控制，而含冰量直接与温度相联 系。在人类生产活动的范围内( 一般不超过2 0 m ) ，由于气候季节变化，引起土 温的变化是不可避免的，因此冻土的性质随时都在变化，表现为动态特性。 冻土区有丰富的土地、森林和矿藏资源，它的存在及其演变对人类的生存 环境、生产活动和可持续发展具有重要影响。根据冻土存在时间的长短，地球 上主要分布着两种冻土：一种为多年冻土，两年以上处于冻结状态，只有表层 几米的土层处于夏融冬冻的状态；另一种为季节冻土，只在地表几米范围内冬 季冻结夏季消融。 冻土的分布面积约占全球陆地面积的2 3 ，主要分布在俄罗斯、加拿大、 中国和美国的阿拉斯加等地，其中我国的多年冻土分布面积约占世界多年冻土 分布面积的1 0 ，占我国国土面积的2 1 5 ，季节冻土分布面积占国土面积的 5 3 5 n 我国冻土研究起步较晚，最早见于1 9 5 8 年出版的多年冻土的工程地质和 铁路建筑一书，是对东北多年冻土地区铁路建设的总结。从1 9 6 0 年起，铁道 部科学研究院西北研究所与铁道部第一勘测设计院、中国科学院寒区旱区环境 与工程研究所等有关部门，开始对青藏公路沿线冻土进行考察和定点观测研究， 尤其在铁路科学研究大会战期间，对青藏铁路沿线进行了多学科的综合性研究， 取得了大量的研究成果。2 0 世纪7 0 年代中期到2 0 世纪9 0 年代中期，随着中央 决定青藏公路修建工程的实施，取得了丰富的科研成果。在高原多年冻土地区 第一章绪论 公路工程地质研究中，采用探地雷达与钻、坑探相结合的综合勘探方法，取得 了公路沿线详细、连续的冻土工程地质资料，根据青藏公路冻土温度与气温、 海拔、纬度的关系，提出了路基工程地质分区、分类，对冻土环境和多年冻土 地温带的变化趋势进行预测。我国冻土科学研究工作者结合工程实践开展了大 量野外调查分析研究、室内实验研究和计算机数值分析，取得了丰硕的成果， 突出的表现在土体冻胀和盐胀机理研究及冻害盐害防治等方面。研究成果除以 论文的形式发表在国际多年冻土会议论文集和国内专业性杂志上外，比较集中 的反映在土的冻胀与建筑物冻害防治( 童长江，1 9 8 5 ) 、冻结凿井冻土壁的 工程性质( 吴紫汪，1 9 8 8 ) 、冻土的温度水分应力及其相互作用( 安维东， 1 9 8 9 ) 、冻土中水分迁移的实验研究( 徐学祖，1 9 9 1 ) 、冻土强度和蠕变( 吴 紫汪，1 9 9 4 ) 、冻融土中水热输运问题( 李述训，程国栋，1 9 9 5 ) 、多年冻土 退化与道路工程( 藏恩穆，1 9 9 9 ) 、冻土物理学( 徐学祖，2 0 0 1 ) 等专著中。 国外俄罗斯是世界上最早注意到冻土上的建筑条件的国家。在1 9 0 4 年时修 建的阿穆尔和外贝加尔铁路时遇到了大量的冻土问题，1 9 1 2 年俄罗斯学者出版 了永久冻土与永久冻土上的建筑物一书，阐述了在永久冻土条件下建造建 筑物的方法，为冻土学的发展做出了巨大的贡献，促进了俄罗斯的冻土研究工 作。1 9 2 7 年m i 苏姆金发表的苏联境内的多年冻土专著，标志着冻土成为 一门专门的学科。1 9 5 8 年，n i a 崔托维奇出版的专著冻土地基与基础比较 系统的阐明了冻土的基本物理性质、力学性质、冻土中的物理力学过程及冻土 上基础设计、施工方法。1 9 7 0 年原苏联科学院出版了房屋及构筑物与永久冻 土的相互热作用一书，论述了冻土融化深度和正融土沉降预报。 美国陆军部寒区研究与工程实验室( c r r e l ) 从二战开始，为解决阿拉斯 加等寒区的道路、机场及其它军事设施的建筑问题，开展调研积累多年冻土区 特殊工程的设计标准与施工准则的素材。1 9 8 0 年c r r e l 出版了深季节冻土区 和多年冻土区基础设计和施工一书，提出了深季节冻土区和多年冻土区特殊 工程设计标准和施工准则【1 4 l 。 在加拿大，冻土研究的蓬勃发展也主要源于极地多年冻土区石油天然气等 资源的开采，在工程实际建设当中，与采用融化的方法比较起来，偏重保持地 基土冻结状态的方法。当然在公路运营期间，保持地基土的冻结状态，对保证 路基的稳定性而言，将形成最有利的条件。 一2 一 第一章绪论 1 1 2 冻土温度场的发展状况 冻土温度场已有1 7 0 多年的研究历史，但直到1 9 世纪前期，冻土及其温度 场的研究基本仍为初步探索过程。1 8 9 0 年俄国成立了冻土研究委员会，此后开 展了较广泛的研究。进入2 0 世纪，苏联又经历了1 9 4 5 1 9 6 0 年和1 9 6 1 1 9 7 1 年两个较快的发展时期。开展了与温度有关的热力学、热物理学、土壤水热改 良、工程建筑稳定性以及地球表面和岩石圈温度场的形成等方面的实验研究和 以解析解为主的理论计算研究。2 0 世纪7 0 年代后，计算机和数值方法在苏联冻 土领域广泛应用。使以前难以解决的复杂几何形状和地质条件，考虑热质交换 的非线性问题的深度和广度都有了新的发展。真正开始理论性研究并作为这门 学科理论的奠基人为c y m t h h m m 。在温度场等热物理研究方面出现了 l 渐p 且b l i e b b a 等代表人物。 北美、西北欧等其它国家和地区，与苏联一样主要是自然资源的开发推动 了这些国家冻土温度场及其有关学科研究的发展。2 0 世纪阿拉斯加金矿的开采、 1 9 4 2 年北美战备公路的严重冻害，开始了本课题从理论上较全面的研究。在加 拿大，这项研究的蓬勃发展也主要源于极地多年冻土区石油天然气等资源的开 采。2 0 世纪7 0 年代这些国家相继进入研究的高潮，除自然资源开发的需要外， 其主要原因之一是现代监测手段和计算技术在冻土研究中得到广泛应用，并求 得了一维非线性温度场的数值解。同期还发展了与温度有关学科的研究。 r e y n o l d s 和e m a n u e l 对表面效应与温差的特性进行了研究，随后的研究主 要是预测不同桥型的温度场分布。1 9 7 5 年，k e h l b e c k 假定热流动是一维的，给 出了这类问题的理论解。e m e r s o n 提出了有限差分模型，h u n t 和c o o k e 对其进 行了改进。许多学者都研究过管道表面结冰的问题。e l m e r 提出了与此问题相关 的准静态解；p e k e r i s e t a l 给出了一阶正确准静态解；h w a n g 考察了由c a r s l a w e t a l 提出的埋设管道冻结的二维数值解的正确性；正冻和正融土热物理学一书 ( h o p x a e b ，1 9 6 4 ) 论述了正冻、正融土中的热质迁移。随后，冻土温度状况 计算方法( 西c j i b 删a h ，1 9 7 3 ，徐学祖等译，1 9 8 2 ) 一书讨论了关于冻土温度场 的多种计算方法，为以后人们计算、研究提供了理论依据。 我国冻土温度场的研究起步较晚，正式开始温度场理论研究仅有5 0 余年历 史。2 0 世纪5 0 年代余力教授开始对人工凿井温度场的研究；徐攸( 1 9 6 2 ) 研究 了天然细沙地基的温度场。在此期间，主要开展室内外观测和实验的经验方程 3 第一章绪论 的计算研究。用稳定、非稳定一维线性问题的解析法和半解析法进行研究。2 0 世纪6 0 年代，z n k 开始研究公路桥的温度行为通过对几座桥梁的调查，认为桥 梁的温度场分布大气温度、风、湿度、太阳辐射、材料。 2 0 世纪8 0 年代以后，世界各先进寒区国家的冻土研究利用其它学科成就推 动本学科向前发展。主要向着多维多相非线性问题和多场相互作用问题理论模 型的建立和求解发展：向着研制和应用现代化高效能高精度的试验技术发展。 在温度场计算方面国内外学者专家们做了大量的工作，如2 0 世纪8 0 年代杨诗 秀等人的土壤冻结条件下的水热耦合迁移的数值模拟；安维东等人的渠道冻结 时热质迁移的数值模拟；1 9 8 2 年日本北海道大学低温科学研究所木下诚一论述 了土体冻胀机理、冻结和冻胀的数学解析及人工冻土利用。国内土的冻胀与 建筑物冻害防治( 童长江、管枫年，1 9 8 5 ) 、冻土的温度水分应力及其相互作 用( 安维东等，1 9 8 9 ) 、冻结凿井冻土壁的工程性质( 吴紫汪、丁德文等， 1 9 8 8 ) ；2 0 世纪9 0 年代，李述训的立井冻结法凿井工程中的热工计算，对不同 的温度场进行了计算。何春雄等建立了隧道内空气与围岩对流换热及固体传热 影响的综合模型，分析了隧道内空气分别为层流和紊流的情况时，隧道内气温 及围岩冻结、融化状况；赖远明等根据传热学和渗流理论导出了带相变的温度 场和渗流场耦合问题的控制微分方程，应用伽辽金法导出了这一问题的有限元 计算公式，以大坂山隧道为例考察了渗流和保温层对隧道围岩冻深的影响规律。 陈尚桥、黄润秋用有限元反演的方法对深埋隧洞地温场进行了数值模拟研究， 并对进一步掘进时的地温场变化作了预测。何春雄等曾就壁面温度随气温周期 性变化的情况，分析计算了隧道围岩的温度场。冻土中水分迁移的实验研究 ( 徐学祖、邓友生，1 9 9 1 ) 、冻融土中的水热输运问题( 李述训、程国栋，1 9 9 5 ) 等几本专著从不同侧面论述了有关冻土温度场的问题，在冻土温度场研究方面 做出了重大贡献忙1 。 长期的工程实践和科学研究，使人们认识到冻土温度场是冻土生长条件、 冻土工程、冻土环境、冻土改造和利用以及冻土预报等问题研究的基础。它将 普通冻土学和工程冻土学密切联系，成为深入研究冻土学的纽带。 自从g a l i e l o 系统地引进实验方法和n e w t o n 奠定物理学的理论基础以来， 近代科学方法一直分为实验研究和理论分析两大方面。实验研究观察实体本身， 理论分析求出由数学模型所代表的结果。因此，一个物理过程的最可靠数据资 料应该由实验测量得到。然而，采用实物进行实验研究是极其昂贵的，有些条 4 第一章绪论 件下甚至是不可能实现的；而模型实验并不是总能模拟实体的各方面特征，测 量在许多情况下存在着相当大的困难，测量仪表本身也存在着误差，且所得结 果推广到实体时无一般规律可循，所以，实验研究效能不高。必须与其它方面 相结合。理论分析基于根据实际的物理模型和实验结果建立的数学模型，即描 述物理过程的一组微分方程。尽管数学模型解析解的结果具有普遍意义，但由 于实际问题中只有很少一部分方程组可以获得解析解，而得到的解析解往往含 有无穷级数、特殊函数以及关于特征值的超越方程，应用价值不高。因而，经 典的数学方法无疑不能为求解现代科学、技术、工程中的大量数学模型提供足 够的实际可行的途径。 就理论研究文献而言，冻土温度场研究已经从早期的定性描述发展到今天 利用计算机进行数值模拟分析。期间经历了观测实验经验方程、大为简化的计 算近似处理、均质一维、二维线性稳定问题和一维非稳定线性问题的解析计算， 以及近期的数值模拟研究等发展过程。 经验方程一般为观测实验结果的统计表达式。它主要考虑某些重要影响因 素，对其它的影响因素则以修正系数的形式表示。近似计算对问题做出抽象后， 在数学上大为简化了计算公式。因此他们均是以满足一股工程计算为目的的近 似计算方法。在一般的工程课题计算中，不可能要求计算公式所得的精度比给 定的基本参数更高。当然要考虑到计算对所计算问题的影响程度和是否满足实 践要求。此外，在冻土学科发展前提下，计算手段落后以及野外考察条件的限 制，往往以手工计算为主。为了计算简便起见，对在工程精度允许范围内的近 似公式一直得到普遍的推广和应用。 线性稳定温度场理论主要研究冻结层和融土层内单向和平面的温度分布问 题。如近似的研究最大冻结深度和最大融化深度内拟定稳定阶段的温度场问题， 以及多年冻土上限以上基本无相变的冻土温度场问题。稳定温度场忽略了时间 变化产生的影响。由于根据稳定温度场建立起来的公式一般比较简单，使用也 比较方便，所以凡是能用稳定温度场公式解决问题的地方尽量使用这一公式。 但要严格把握应用范围。 为了解决冻土温度场随时间变化和相变影响问题，发展了一维线性相变问 题的解析解。自从非稳定相变温度场理论问世以来，冻土温度场的解析方法研 究有了新的突破。非稳定温度场理论比稳定温度场理论更趋完善，更符合实际。 它最大的优越性是可以解决整个过程的发展变化情况，迄今已广泛应用于生产 第一章绪论 实践中。但一般来说，解析法只适用于几何形状简单，均质各项同性的情况。 目前最大时间尺度和平面、空间等非稳定温度场问题方面解析计算尚很困难， 因而限制了它的应用范围。这对深入理论研究提出了解决计算方法的课题。 随着电子计算机的出现和飞跃发展以及伴随的计算方法的不断改进，正在 把数值模拟方法推向人类科学研究活动的前沿。借助计算机和有限元、有限差 分、边界元等数值计算方法，人们几乎能对任何一个描述实际问题的数学模型 求出数值解，数值模拟方法己上升为一种对科学的定量研究起重要作用的、与 实验研究和理论分析相并列的第三种科学方法。在物理过程可以用确定的数学 模型近似描述的前提下，数值模拟方法具有其它方法无可比拟的一些优点：可 以模拟真实条件和理想条件的状态；成本低；计算速度快；获得的数据资料完 备。因而数值模拟方法被认为是一种大有前途和卓有成效的方法，采用数值模 拟方法可以较少实验的工作量，拓宽研究的范围，深化研究的内容和提高研究 成果的质量。从某种意义上说，在特定参数下进行一次数值模拟相当于进行了 一次试验，历史上就曾有过首先通过数值模拟发现新现象而后又由实验验证的 例子。2 0 世纪7 0 年代以来，数值计算方法在冻土温度场研究中的采用，使以前 由于几何形状复杂、地质条件特殊而难以解决的有关冻土温度场的一些难题， 比如热质交换的非线性问题的研究在深度和广度上都有了较大的发展。目前， 数值模拟方法已被认为是冻土预报中最有效、最基本的方法之一。因此，在冻 土温度场的研究中有必要采用数值模拟方法。事实上。即使在进行实验研究时， 人们也希望同时获得数值模拟求得的数值解，以补充资料和进行理论分析。可 见，在冻土温度场的研究中，数值模拟方法是有力的工具，进行数值模拟研究 是完全必要的。 为适应寒区科技的发展和面向国民经济建设主战场的需要，我国冻土学研 究正在实现从工程到资源、环境，从描述到预测、预报，从认识到改造、利用 的三大转变，我国寒区的开发和道路工程的修建、养护与提高道路等级工作的 发展，促使冻土热状况动态特征数值模拟的研究不断深入。 6 一 第一章绪论 1 2 冻土的基本理论 1 2 1冻土的成分、构造及其力学性质 我国北方和西部与西南高山地区，分布着广大的冻土地带。在这些地区的 一切生产活动和经济活动、自然环境的利用和保护，都受它的影响。地基土的 冻结融化等地质作用，给建筑工程设施造成的危害，已越来越引起重视。通过 长期的工程实践和科学研究+ 人们逐步认识了冻土冻结融化的变化规律及其有 关特性，了解冻土的有关基本原理和特性，对工程设计、施工与管理都是十分 必要的。 1 ，2 1 1 冻土的成分与构造特征 冻土是较复杂的四相体系。其基本成分是固体矿物颗粒、粘塑性的冰包裹 体、液相水( 未冻水和强结合水) 和气体包裹体( 水气和空气) 。 ( 1 ) 固体颗粒是冻土的骨架成分，对冻土性质( 如冻土矿物颗粒尺寸、形状、 分散度、矿物成分等1 有重要影响，而这些特征对冻土的物理力学、物理化学性 质又有直接影响。 ( 2 ) 冰是构成冻土的必需组分，它是胶结土粒骨架的固态水。冰是一种物理 化学性质极为特殊的单矿物低温水化岩石，结构很不稳定，受温度影响极大。 当温度降低达到一定负温时，极为致密、坚硬的有序结构，当温度升高时，分 子间的连结作用削弱，冰的强度降低。冰有强烈的各向异性，当力作用在不同 方向时表现出冰的粘塑性变形与脆性。冰在很小的荷载作用下产生流动变形， 只有在瞬时荷载作用下才有弹性。冻土温度或作用荷载变化时，还发生冰一水 相转变。冻土中的冰相成分制约了冻土的结构特征和物理力学性质。 ( 3 ) 冻土中的液相水主要是未冻水。在一般负温条件下，冻土中总含有未冻 结的液态水，它主要存在于冻结细粒土中，是土粒表面吸附的强结合水。冻土 中未冻水与土的矿物成分、可溶盐含量、冻结温度、外部荷载作用有关。冻土 中的未冻水含量及其在外界作用下的变化，多方面制约着冻土的物理、热学、 电学及力学性质。 ( 4 ) 冻土中的气态成分是水气和空气。一般仅在冻土未被水和冰完全饱和的 情况下才存在水气，土中吸附的气体与土的矿物成分、分散度、孔隙率及有机 物含量有关。 一1 一 第一章绪论 上述冻土中的各相成分相互作用，影响和决定了冻土的各种性质。根据总 含水率、土的类别，多年冻土可分为少冰、多冰、富冰、饱冰冻土和含土冰层， 青藏铁路高原多年冻土区工程设计及施工暂行规定中，对多年冻土分类如表卜1 所示。 表1 - 1 多年冻土分类 冻土类型土的类别总含水率w a 融化后的潮湿程度 粉黏粒含量= 1 5 的粗颗粒土，a 1 5 的粗颗粒土w a 1 2 = 细砂、粉砂w a 1 4稍湿 ”?粉土w a 1 7 一 黏性士w a w p坚硬 + 粉黏粒含量= 1 5 的粗颗粒土 1 0 = w a 1 5 的粗颗粒土1 2 = w 1 5 =细砂、粉砂1 4 = w a 1 8潮湿 ”? 粉土 1 7 = w a 2 1 一 黏性土 w p - - - - w a ：w p 4 - 4 硬塑 富 粉黏粒含量= 1 5 的粗颗粒土 1 5 = w a 2 5 饱和出水 ( 出水量小于1 0 ) 冻 细砂、粉砂 1 8 = w a 2 8 饱和 士 粉土 2 1 = w a 3 2 塾堂圭! ! ! ! 苎丛三监! ! !墼望 饱 粉黏粒含量- 1 5 的粗颗粒土2 5 吡 4 4 f 出水嚣譬豪2 0 ) 鋈黼粒冀圣扩蹴土2 5 = w ：a 7 4 2 8 = w 1 5 的粗颗粒土， 粉黏粒质量 1 0 的细砂 不考虑不考虑 w = 1 2 1 0 ， i 级 口= l 不冻胀 粉砂 1 2 1 0 粉土 w = 1 9 1 5 黏性土 w = w p + 2 2 0 粉黏粒含量 1 5 的粗颗粒土，w = 1 2= 1 0 粉黏粒质量 1 0 的细砂1 2 1 0 粉砂 粉土 黏性土 w = 1 4 1 4 w = 1 9 w = 1 9 1 9 1 5 = 2 0 2 0 ， 1 5 的粗颗粒土 粉黏粒质量 1 0 的细砂 粉砂 粉土 1 2 1 8 1 4 w = 1 9 1 9 w = 2 3 1 9 1 o = 1 5 1 5 3 5 r = 6 燃 黏性土 w r + 2 一 ，8：0 t 5 粉黏粒质量 1 0 的细砂 ”，一1 “。 粉砂1 9 w = 2 3= 1 0 粉土 2 2 1 2 特强冻 胀 一1 3 一 第一章绪论 ，7 ：等1 0 0 ( 1 - 7 ) 二d z 。= h 一z( 1 - 8 ) 式中，a z 地表冻胀量； z 。设计冻深( m m ) ； h 冻层厚度( m m ) 。 1 2 2 4 冻土的主要力学特性 ( 1 ) 分类 冻土对工程建筑物稳定性产生影响的主要原因是冻土的融沉作用，故其工 程上的分类主要是根据冻土层的总含水量和融化后的沉降性评价而进行融沉分 类的，按融沉分类为：不融沉、弱融沉、融沉、强融沉和融陷等五类。 ( 2 ) 主要物理力学性质 冻土是由矿物颗粒、固体水、液相水和气体组成的四项体系，其物理力学 性质主要取决于冻土中的含水量( 包括固相和液相含水量) 、冻土中的未冻水的含 量、原状结构冻土的重度、固体矿物颗粒相对密度等四个方面特征。 借鉴国内外多年来冻土物理力学性质的研究，对影响工程质量的冻土的物 理力学性质及其相关因素进行分析总结如下： 冻土的抗压强度与负温度、矿物成分、及颗粒组成的关系。冻土随着负 温度的降低，其极限抗压强度急剧增加；从一5 起随着负温度的降低其极限 抗压强度基本呈现直线变化。 冻土抗压强度与总含水量的关系。当负温度恒定时，随着冻土含水量的 增加其抗压强度也增加：孔隙全充满冰时抗压强度最大，超过这个限度反而降 低。 冻土的抗压强度与外力作用时间的关系。冻土瞬时抗压强度很大，但在 长期荷载的作用下就会降低很多，当应力小于某一界限时，荷载长期作用下也 不破坏。 冻土抗剪强度与负温度、外压力、荷载作用时间以及含水量的关系。冻 土的抗剪强度长期随负温度的降低而显著增大；在一定范围内与法向压力呈直 线关系；冻土的抗剪强度比瞬时抗压强度小的多；在冻土含冰量饱和之前，其 含冰量的增加将增加其抗剪强度。 第一章绪论 冻土的变形性质。多年冻土在短期荷载作用下，压缩变形很低，类似岩 石，一般不计变形，但在长期作用下，其变形相当大。 融化下沉。在热力作用下，冻土逐渐融化，导致冻土的一系列物理力学 性质的变化，在自重和外加荷载的作用下产生融化下沉和压密现象，这将使地 基产生不均匀沉降。 冻土物理力学性质表明冻土的力学性质是不稳定的，它赢接取决于其外部 环境是否发生变化以及变化的方向。在施工过程中，不可避免的会引起多年冻 土外部环境( 温度、作用力、含水量等) 的变化，从而造成不稳定因素的产生。因 此，如何解决冻土地区的施工技术是工程成败的关键。 1 3 选题目的及意义 随着我国改革开放的深入，经济的迅速发展以及西部大开发战略的实施， 使得青藏铁路的修筑迫在眉睫。青藏高原是世界上面积最大、海拔最高的高原， 地理位置独特，自然环境恶劣，日温差大，地质条件复杂，1d 中出现正负温的日 数每年有1 8 0d 左右( 比东北严寒地区多1 倍左右) 。青藏铁路穿越多年冻土地 段，冻土问题是青藏铁路的一大难题，其特殊性和复杂性在世界上是独一无二 的，素有“世界屋脊”、“地球第三极”之称。 冻土问题是青藏铁路的一大难题，主要是指土中固液转化以及冷生构造作 用，给常规工程措施和施工带来的困难。其特殊性和复杂性在世界上独一无二。 在青藏高原多年冻土地区修建铁路是从未实践过的领域，随着几十年来自然条 件和气候的变化、科学技术的突飞猛进、科研成果和工程实践经验的积累，我 们对自然和冻土的认识也在不断加深。 冻土与融土的最重要差别在于冻土中含有冰体，并且正冻土和已冻土中存 在着温度梯度作用下的水分迁移现象。由于冻土中冰的胶结连接作用，保持负 温不变条件下冻土具有较高的坚硬度和强度，但当温度升高或降低时( 包括在 负温范围内) ，冻土的性质将发生重大变化，从而决定了冻土的强度具有不稳定 性。当冻土融化时，结构中的冰胶结作用随即产生破坏，冻土显著变形，压缩 性和透水性显著增大，强度和稳定性明显降低，融化处发生快速的局部下沉， 并常伴随有已融土的稀释体从承重面挤出；而当土壤冻结时，伴随土中孔隙水 和外界水源补给水结晶形成多晶体、透镜体和冰夹层等形成的冰侵入体，土体 第一章绪论 体积增大，导致地表产生不均匀冻胀，由此将对修筑于冻土上的各类建筑物产 生根本无法允许的变形和破坏。 长期的工程实践表明，青藏铁路成败的关键在于路基工程。青藏铁路绝大 部分是路基工程。而路基的强度和稳定性，是确保铁路正常运营的首要条件。 多年冻土地区道路的修筑，改变了地表和大气之间的热交换关系。根据多年的 观测资料，虽然青藏高原多年冻土区的气温冻结指数比融化指数大的多，但是 在填土路提的高度小于一定高度时，地基的冻结系数较低，从而使多年冻土区 的冻土融化，冻土上限下降，导致路基沉降，路面毁坏，严重影响道路运营。 因此，高原冻土区铁路路基是按照以下三个基本原则进行的，应该单独应用或 综合考虑。一是保持冻土处于冻结状态，及保护冻土原则。对于路基而言，是 指应用该原则设计的年限内，能保持热稳定性。即人为上限始终控制在指定的 范围内，保持其下卧多年冻土的冻结状态。二是控制多年冻土逐渐融化或局部 融化，即延缓融化。对路基而言延缓融化原则是指在设计使用年限内，允许设 计路基基底的多年冻土逐渐完全融化或产生局部融化，可以将融化速率和深度 控制在路基稳定性所允许的变形范围内。三是预先融化或清除多年冻土，即破 坏多年冻土原则。即将路基基底多年冻土融化或清除，并将融化后的水分疏干。 多年冻土区的工程实践和科学研究使人们逐渐深刻地认识到，冻土温度场 是冻土生成条件、冻土工程、冻土环境、冻土改造和冻土预报的基础，是联系 普通冻土学和工程冻土学的纽带。与土壤冻结和融化有关的水分迁移过程、成 冰规律和冻胀融沉现象等在一定条件下实际上完全由过程的热状态决定。无论 是研究冻土的区域分布规律、物理力学性质，还是选择冻土地区工程建筑物设 计的原则、确定具体的设计参数，都需要了解并掌握冻土的温度场分布情况， 因此，对冻土温度场的研究是工程冻土学研究的主要内容之一。 研究围岩冻土温度场的根本目的就在于预测冻土温度场的变化规律及其可 能导致的融沉与冻胀破坏，提出调控冻融过程的方法和治理各种冻害的措施， 确保高寒区、冻土隧道工程具有稳定性、耐久性和经济合理性。尽管这一问题 的研究已经取得了很大进展，但在全球范围内，冻土的融沉和冻胀问题迄今为 止仍然是严重影响冻土隧道工程稳定性的突出问题，许多方面的研究仍处于空 白或需要进一步深化。 长期以来，我国多年冻土区的道路工程主要是按保护冻土原则设计和修筑 的，在气候变暖、冻土退化和人类各种工程与经济活动对冻土温度状况的影响 一1 6 一 第一章绪论 不断扩大的背景下，冻土温度场的研究正面临新的挑战：如何确保寒区现有隧 道工程的稳定性，冻土区新的隧道工程应怎样设计、施工和养护，这需要冻土 温度场的研究尽快作出准确的回答，并提出解决问题的原则与具体措施。 与气候变化对冻土温度场的影响相比，人类的各类工程活动和经济活动的 影响更直接、更大。人类工程经济活动极大地影响了寒区土的冻融作用，造成 了脆弱冻土环境及生态环境变化。由于人类的工程活动改变了多年冻土的热量 平衡，引起多年冻土地区路基变形，特别是引起高温使多年冻土地区路基热融 下沉问题。这是至今还未彻底解决的一大难题。因此为保证冻土地基的稳定性， 研究其内在的温度场分布规律就变得十分迫切。建立符合实际的冻土路基温度 场，确保冻土地基的稳定性，研究路基沉降受冻土上限移动变化和地温的影响 规律，从而达到控制路基不均匀变形的目的，同时也在为在该地区进行路基施 工及其其他建筑工程提供重要的参考。1 9 9 9 年1 月5 7 日在北京召开的中国铁 路重大工程地质问题研讨会上，“铁道工程作用下高原冻土的融沉及其危害”问题 已被我国专家列为铁道工程地质研究的重大课题之一，足见研究冻土温度场的 重要性和紧迫性， 在多年路基施工过程中，为了确保地基的稳定性，多年冻土的融化深度越 小越好。因此，计算冻土地基温度场的变化规律，得出多年冻土的融化圈，进 一步计算地基应力和位移，提出有效的旌工方案是必要的，为以后多年冻土区 工程的设计和施工提供理论依据和参考。 1 4 课题研究内容和方案 1 4 1 研究内容 ( 1 ) 研究冻土的有关原理和特性，及冻土的成分、构造、力学性质和工程特 性。 ( 2 ) 通过研究热力学、热传导学，确定冻土温度场的模型理论。 ( 3 ) 通过大量现场实测数据基础上，研究多年冻土地基温度分布规律及冻土 冻融特性和地温分布规律。 ( 4 ) 运用a n s y s 有限元分析软件，建立冻土路基基温度场的模型。在 a n s y s 模型上，通过数值分析模拟冻土路基温度分布规律及其对路基沉降变形 1 7 第一章绪论 的影响。 ( 5 ) 通过对模型模拟数据和实测数据的对比，预测远期路基稳定性状况，研 究控制路基变形的措施。 1 4 2 研究方案 通过结合冻土、热力学和热传导学的基本理论，建立冻土路基温度场的微 分方程，得出冻土路基的理论模型。 分析现场实测数据，确定冻土类型，得出冻土地基温度分布规律。找出影响 地温场变化的主要因素及一般导热微分方程与定解条件，确定路基边界条件， 得到围岩温度场计算程序及参数，从而创建a n s y s 计算模型。通过模型得出模 拟数值结果。 分析模拟结果。通过对模型模拟数据和实测数据的对比，预测远期路基稳定 性状况，研究控制路基变形的措施。 一1 8 第二章试验段的工程地质特征 第二章试验段工程地质特征 青藏铁路沿线多年冻土位于中低纬度地区，由于其海拔高，气温低，降水 少，导致多年冻土比较发育。其中格尔木至拉萨段将穿越5 0 0k m 的多年冻土地 段，由于沿线地质( 包括构成地层的岩性、含水量和地中热流) 、地理( 地形、 植被、雪盖、地面条件等) 条件的不同，在空间上冻土的状态可能出现巨大差 异。清水河试验段位于青藏铁路格拉段昆仑山口至风火山之间的广阔的楚玛尔 河高原准平原区，起讫里程为d k l 0 2 4 + 9 0 0 d k l 0 2 5 + 0 0 0 ，长度1 0 0m 。具有 下列工程地质、水文及气候特征。 2 1试验段地质、水文及气候概况 2 1 1 地貌及水文特征 该试验段位于北纬3 2 。3 6 。的( 长) 江( 黄) 河源多年冻土区的片状多年 冻土带中，平均海拔45 0 0m ，地势开阔平坦，属典型的高原河流相沉积地貌， 植被覆盖率为1 0 ，主要以高寒草原和垫状植被为主：热融塘湖分布