延时摄影因可靠、高效和低成本的优势,在冰川监测中应用广泛,特别是对于获取冰川表面连续变化信息而言。本文基于2020年3月—2021年9月物候相机拍摄的梅里雪山明永冰川末端照片及多期无人机影像,利用地面摄影测量技术和互相关算法,提取了日尺度冰川表面运动速度。结果表明:通过物候图像获取的冰川表面运动速度分辨率高,从海拔2 880~3 150 m a.s.l.,冰川总位移介于(129.38±7.76)~(669.95±247.88) m,年均表面运动速度达(79.14±4.74)~(412.86±152.75) m·a-1,呈从中间向两侧减缓的空间分布特征。冰川表面运动速度随季节变化,夏季流速[(0.13±0.06)~(1.99±0.37) m·d-1]快于冬季流速[(0.07±0.06)~(1.35±0.37) m·d-1]。与冬季流速相比,夏季流速受降水和气温升高的影响不稳定。根据流速分离结果,明永冰川末端底部全年处于融化或压融状态,底部滑动对冰川表面运动速度的贡献介于76%~93%。冬季底部滑动占表面流速高达82%,夏季底部滑动对冰川运动起绝对主导作用。本文采用的技术为进一步研究季风海洋型冰川的运动机制提供了参考方案。
利用1979—2016年ERA-Interim的月平均再分析资料,通过统计分析、诊断研究等方法,分析了青藏高原水汽收支的气候变化特征及其成因。结果表明:青藏高原水汽四个季节均是从南、西、北边界流入高原,从东边界流出高原;青藏高原春、夏、秋季都表现为一个明显的水汽汇,夏季是高原水汽输送最活跃的季节,而冬季水汽从高原向外输出;高原东部(南部)水汽输送和收支量远大于高原西部(北部)。青藏高原年和夏季有逐渐变湿的趋势,主要是由于东边界水汽流出量的显著减少;高原东、西部的水汽净收入均呈增加趋势,高原西部的增幅明显大于高原东部;东边界水汽流出量的减少是高原东部水汽净收入增加的原因,而南边界水汽流入量的明显增加是高原西部水汽净收入增加的原因;高原北部水汽收支的年际变化比高原南部更明显,高原北部有明显增湿的趋势,而高原南部增湿趋势不明显。夏季青藏高原南边界水汽收支变化主要受青藏高原以南、以东的季风控制区偏南风水汽输送的影响,与西太副高的变化也有密切联系。西边界水汽收支的变化主要受中纬度西风带水汽输送的影响,北边界水汽收支的变化主要受到“丝绸之路”遥相关波列的影响,东边界的水汽输出偏少与高原地区东北侧至贝加尔湖附近上空的反气旋异常和华南地区上空的气旋异常有关。
近相变区冻土导热系数是影响寒区岩土工程热稳定性的关键物性参数,针对该温度区间冻土导热系数难以精确测量的问题,提出了一种稳态热流计法试验获得导热系数、量子粒子群寻优算法近似回归定性温度的近相变区冻土导热系数测试方法,采用核磁共振技术对细砂土近相变区未冻水含量进行测试并探究了其与导热系数的关系,分析了负温区细砂土导热系数的分布特征和参数影响规律,建立了不同温度区间的导热系数预测模型。结果表明:负温区细砂土未冻水含量变化曲线与其导热系数呈现出极为相似的变化规律,温度-导热系数曲线可划分为-3.2 ℃为界的近相变区和稳定相变区;冻结细砂土导热系数的增幅区间主要发生在近相变区,增长幅度约占整个负温区导热系数总增幅的50%,且随温度降低,导热系数的增幅逐渐减小;预测模型平均相对误差为4.14%,±10%误差内占比为94.7%,稳定相变区确定系数R2为0.892,近相变区确定系数R2为0.883,证明了该模型具有较好的准确性。研究成果可为寒区工程温度场的计算与分析提供基础参数,亦可为寒区重大基础设施精细化设计和冻土风险精细化防控提供数据支持。
基于质子转移反应质谱仪(PTR-MS)和颗粒物(PM)采集和分析技术,对我国西北甘肃省和新疆维吾尔自治区(甘新地区)十三个城市大气中的挥发性有机物(VOCs)、PM(PM10、PM2.5和PM1)及其组分(水溶性无机离子、碳质气溶胶和无机元素)进行了采样分析,阐明了甘新地区城市大气中VOCs和PM的污染特征和健康风险。研究表明:污染特征分析表明平均总VOCs(TVOCs)浓度为(41.84±7.56) ppbv(ppbv为十亿分之一的体积混合比),氧化性VOCs(OVOCs)是VOCs的重要组分,VOCs组分总浓度高于国内外其他城市,主要为甲醇的浓度较高,而芳香烃则低于其他城市。PM10、PM2.5和PM1的平均浓度分别为(139.39±32.63) μg·m-3、(77.66±25.39) μg·m-3和(44.76±17.59) μg·m-3,水溶性无机离子(WSIIs)是该地区PM的重要组分。健康风险评估表明该地区VOCs的非致癌风险显著,致癌风险处于可接受水平;PM中重金属的非致癌风险和致癌风险均显著。乙醛的非致癌风险和萘的致癌风险较高,Mn的非致癌风险和As的致癌风险较高,因此应加强乙醛、萘、Mn和As元素的管控,以减少大气污染对人体健康的危害。
北极海冰范围退缩已对区域和全球气候变化、北极通航性乃至地缘政治格局产生了深远影响,开展北极海冰预测具有重要意义。本文使用美国冰雪数据中心(NSIDC)发布的海冰密集度、海冰冰龄遥感资料以及ERA5再分析资料,结合前人研究和海气耦合模式结果选择预测参量方案,开展了月尺度的海冰冰情预测。比较了支持向量机(SVR)、深度森林(DF)、LightGBM (LGB)、XGBoost (XGB)和CatBoost (CAT)等5种机器学习算法和以树模型LGB、XGB和CAT作为基模型,以贝叶斯回归、岭回归、套索回归和深度森林作为元模型的4种堆叠式集成学习模型,以及深度神经网络(DNN)、卷积神经网络(CNN)、时空卷积网络(ConvLSTM)3种深度学习模型在2000年测试集上对海冰范围和密集度空间分布的预测效果。结果表明:在海冰密集度预测中,ConvLSTM表现最优,套索堆叠集成学习模型预测效果次之。集成学习模型相较于三种单一树模型在预测效果上有约1%~4%的提升。在海冰范围预测中,堆叠式集成学习模型的预测效果最好。本研究为开展机器学习海冰预测奠定了重要基础。
气候变暖对北极多年冻土和植被产生了重要的影响。CLM(Community Land Model)是应用最广泛的陆面过程模式之一,但其中复杂的边界条件和参数化过程导致模式模拟结果存在一定的不确定性。本研究评估了CLM5.0对阿拉斯加多年冻土区表层土壤温度和碳循环的模拟能力,结果表明,CLM5.0可以捕捉到表层土壤温度的季节变化。在苔原和针叶林站点,CLM5.0在日尺度和月尺度都可以很好地模拟出总初级生产力(GPP)随时间的变化,但对净生态系统碳交换(NEE)的模拟结果存在一定的不确定性。CLM5.0可以较为合理地模拟高纬度多年冻土区的土壤温度季节变化,在未来的研究中可能还需要从结构、参数化方案等过程进行改进,从而进一步提升高纬度多年冻土区碳循环的模拟精度。
多年冻土区活动层是地表水和地下水相互转化中十分重要的交换通道,活动层土壤含水量是多年冻土区水文循环中重要的组成部分,其动态变化与寒区生态环境密切相关。在气候变化背景下,深入了解活动层土壤含水量的动态变化特征具有重要意义。本文利用ELM(Extreme Learning Machine)模型对青藏高原腹地不同海拔高度多年冻土区土壤含水量进行模拟分析,结果表明:与BP神经网络模型相比,二输入变量ELM模型的模拟精度更高;ELM模型模拟后1天土壤含水量的NSE值在0.69~0.87之间,其中坡下20 cm深度处模拟NSE取得最大值(0.87),并且模拟精度随着推后时间的增加有所提升,模拟后3天和后7天的NSE值分别在0.76~0.92和0.75~0.93之间;坡下各深度含水量的模拟效果优于坡上。在此基础上,通过设置不同的气候变化情景,研究土壤含水量在气候变化背景下的动态变化规律及响应特征。研究发现,升温导致冻结初期以及融化初期不同深度的土壤含水量均出现增大的趋势,在完全冻结期和完全融化期变化不明显。且随着气温增幅的加大,冻结初期以及融化初期的土壤含水量变化也逐渐增大,深层土壤含水量较浅层土壤含水量的增加更加显著。在降水增加的情景下,降水增加越大,土壤含水量的增加趋势越明显,但整体变化幅度较小;坡上各深度土壤含水量的增加主要发生在融化初期和完全融化期,坡下则主要集中在融化初期,相比于深层土壤,浅层土壤对降水增加的响应更加强烈。
青藏高原的许多湖泊发育有湖岸线、侵蚀阶地和高湖面沉积物,这些地貌和沉积单元是反映地质历史时期湖面变化的直接地貌和岩性证据。它们不仅反映了湖泊水文变化历史,还可作为揭示区域乃至全球环境变化的重要载体。在青藏高原众多湖泊中,纳木错最具典型性和代表性。纳木错湖岸最高古湖相沉积物位于现代湖面之上约150 m,高出当今溢流垭口约123 m。至今,有关纳木错异常高湖面的发育时代、湖侵-湖退过程、异常高湖面的形成机制、与其他湖泊之间的水文联系等科学问题仍不清楚或存在很大争议。本文以纳木错异常高湖面的特征、时代、期次和成因问题为主题,综述了一个多世纪以来国内外学者对纳木错湖面变化重建的研究成果,尤其对目前纳木错古湖面重建研究中存在的问题进行了总结,并对未来的工作重点进行了展望。此外,基于初步的调查结果,认为纳木错高湖面沉积可能形成于全新世中晚期,并指示了一次湖侵-湖退过程。南支槽的异常振荡及其导致的冬半年降水增加可能是引起湖面变化的主要驱动因素。
土地利用变化是引起碳排放增加的重要原因之一。探索区域土地利用变化引起的碳排放量和碳排放效率问题,对于把握各土地利用类型的碳排放特征、制定土地利用减排策略具有重要意义。本研究以集中国生态系统重要性与脆弱性于一体的典型区域的河西内陆河流域为研究区,利用1985—2020年的土地利用现状遥感数据,通过土地利用动态变化模型、排放系数法深入分析了河西内陆河流域的土地利用变化及其引起的碳排放量和碳排放效率的时空变化特征。结果显示:(1)1985—2020年绿洲规模呈持续扩张趋势,与1985年相比,2020年的绿洲面积相对增加了13.17%;建设用地以单向转入为主,其余土地利用类型之间的双向转变频繁。(2)各土地利用类型的净碳排放量总体呈增加趋势,2020年净碳排放量(113.17×104 t)相对1985年(48.57×104 t)增加了1.33倍,净碳排放量与绿洲面积显著正相关,耕地和建设用地扩张引起的碳排放量是净碳排放量增加的主要原因。碳排放弹性系数整体呈增加趋势,由1990年的-1.83增至2020年的205.91。(3)1985—2020年的耕地和建设用地碳排放强度、绿洲碳排放强度均呈增加趋势,分别由1985年的32.17 t?km-2和5.09 t?km-2增至2020年的60.74 t?km-2和10.48 t?km-2,但绿洲碳排放强度增加速率小于耕地和建设用地,同期的绿洲碳吸收强度普遍呈递减趋势。(4)工业城市嘉峪关市、金昌市的多年平均绿洲碳排放强度最大,以灌溉农业为主的行政单元相对较低,以畜牧业为主的祁连县最低。本研究通过定量分析河西内陆河流域土地利用的碳排放量及其变化,有助于宏观把握河西内陆河流域土地利用的碳排放特征,对于寻找解决土地利用的减排方案、加强土地管理和推动土地利用的低碳发展模式提供了科学认识基础和决策参考。
在青藏高原多年冻土退化趋势明显及地震活动频发的背景下,为研究多年冻土退化对青藏铁路桥梁桩基础地震易损性的影响,建立了桥梁桩基础-冻土相互作用有限元计算模型,并通过拟静力模型试验验证了模型的合理性。考虑多年冻土退化效应的影响,以桩基础曲率延性比为损伤指标,以地面峰值加速度PGA为地震动强度参数,并选取了符合场地特征的80条地震波,对不同多年冻土活动层厚度和轴压比下的原型桥梁桩基础进行了非线性地震损伤分析。拟静力试验和数值模拟结果表明:地震作用下桥梁桩基础顶部破坏严重,是整个桥墩-桩-土体系的薄弱部位;随着多年冻土活动层厚度的增加(退化程度加剧),周围土体对桩基础的约束减弱,桩身的最大曲率呈减小趋势,从而降低了地震作用下桥梁桩基础的破坏概率,但墩顶位移却有所增大,导致地震过程中落梁风险增加;而随着桥墩轴压比的增大,桩身最大曲率明显增加,多年冻土区桥梁桩基础的地震易损性显著增加。考虑到桥梁桩基础破坏后修复难度较高,因此多年冻土退化背景下青藏铁路桥梁抗震性能评估还需考虑高桥墩轴压比变化对桩基础地震易损性产生的不利影响。
大量工程实践表明,在寒冷山区公路建设过程中,路基阴阳坡效应对冻土区路基的稳定性具有重要影响。以甘肃省南部山区季节冻土路基为例,基于新建宕昌—迭部二级公路试验段的现场监测资料,分析了季节冻土山区公路路基阴阳面的地温和变形差异,提出针对性的两种新型防治措施——措施A(阴坡侧路基下半幅铺设10 cm厚XPS板,阳坡侧路基下半幅铺设6 cm厚XPS板)和措施B(路基下满幅+阴坡坡面铺设6 cm厚XPS板),并通过数值模拟分析了这两种措施控制路基阴阳坡效应的效果。结果表明,在甘肃南部季节冻土山区,山体遮挡作用使公路路基阴阳坡发生转换,对路基横向地温的影响不容忽视,试验段K18+180段阴坡路肩与阳坡路肩地温最大相差6 ℃左右,公路路基横向地温存在显著的阴阳坡效应。地温的横向差异导致公路路基土体冻深和冻胀变形的横向差异,两个试验段路基中心与阳坡路肩的最大季节冻深在第一个、第二个冻结期内均相差约0.8 m和0.9 m,K18+180段阴阳坡路肩冻胀量差值的最大值为2.8 mm,出现在春融初期。路堤高度为1.0 m、2.0 m时,措施A和措施B都能显著地减小公路路基阴阳坡温度差异,并且控制效果明显优于普通措施(满幅铺设8 cm厚XPS板、满幅铺设6 cm厚XPS板)。建议季节冻土山区公路设计时应充分考虑山体遮挡这一局地因素对路基土体冻融状态的影响,并采取阴阳坡两侧差异设计的结构型式,以减少因阴阳坡现象导致的路基病害。研究结果可为季节冻土山区路基阴阳坡效应及其防治研究以及冻害研究提供参考和借鉴。
在中国煤炭是一种重要的能源,大规模的开采活动严重破坏了周边的生态环境,而高原高寒地区煤矿开采对地表温度的时空影响规律有待深入研究。地表温度作为一种重要的地表参数对地表能量和水分平衡研究均具有重要意义,也可以反映区域生态环境的变化。本文以Landsat 5/8影像为数据源,利用面向对象的方法和普适性单通道算法分别得到了7期土地利用/覆盖数据和地表温度数据,研究了1990—2021年木里煤田聚乎更矿区从开采到生态修复过程的地表温度变化特征和影响因素。结果表明:利用面向对象的方法得到7期土地利用/覆盖数据很好地反映了聚乎更矿区的几个阶段,分别为缓慢开采期(1990年、2002年),快速开采期(2005年、2009年),缓慢修复期(2015年、2019年)和快速修复期(2021年);矿坑、建成区、排土场、堆煤场等煤矿土地类型地表温度最高,高寒草甸次之,高寒沼泽草甸地表温度较低,水体温度最低,裸地则在特高温和低温均有分布;草甸区地表温度影响因素的解释力大小排序为NDVI>反照率>坡向>海拔>坡度,矿区地表温度影响因素的解释力大小排序为坡向>反照率>NDVI>海拔>坡度,修复区地表温度影响因素的解释力大小排序为NDVI>反照率>海拔>坡向>坡度。研究工作有助于揭示高原矿区地表温度的变化规律,对于评估木里煤田生态修复效果具有重要的参考价值。
块石作为一种高孔隙率、大粒径的多孔介质,因其可通过调节内部冷暖季节对流方式实现多年冻土路基降温,而被广泛应用于青藏铁路工程以减弱气候变暖及工程扰动对多年冻土路基的影响。因此,精准描述块石路基结构换热机理并模拟其长期热防护效应十分重要。在一些领域内,已经证明使用局部非热平衡(LTNE)可以提升多孔介质的计算精度,但在块石结构孔隙尺度下的传热效能计算中仍较少使用。其中,多孔介质内部流固换热系数(hsf )是LTNE能量方程中重要的参数,该系数如何选取也决定着模型计算的准确性。本文以块石基底路基为例,通过数值模拟方法,将传统路基使用的局部热平衡(LTE)模型与LTNE模型的模拟结果与实测数据进行对比,并分析了不同经典hsf 对块石结构换热的影响。研究结果表明,块石基底路基局部非热平衡效应明显:模拟结果显示LTNE模型中,同一时间内块石与空气温度不相等,空气与块石在夏季和冬季的最大温差高达1.54 ℃和0.80 ℃;LTNE的模拟结果与实测数据更为接近,且LTE模型的降温幅度略高于预期,二者在2 m、5 m、10 m处的年平均误差百分数分别相差6%、1%和8%;不同研究者(Achenbach、Wakao、Dixon等、Pallares等)得出的hsf 在块石路基中变化趋势相近,峰值差异明显,但对块石路基下部土体温度的影响较小。
碎石层的传热机理研究对寒区路基改善技术的应用以及自然界中的“异常”冻土分布的解释有着重要的意义。很多研究者曾对碎石层进行了试验研究,但试验中的温度设置并无统一标准。相较于细颗粒土,粗颗粒材料中水分的影响较小,因此在试验中的温度设置对试验的结果起到了关键作用。为了对比不同的温度设置情况对碎石传热过程的影响,以不同的顶板温度、底板温度以及环境温度进行复合,并对不同组合情况应用COMSOL软件,在考虑碎石层多孔介质传热的情况下进行模拟。结果表明:即使在保温层包裹的情况下,环境温度依然能对温度剖面以及垂直方向上的热通量产生明显的影响。当环境温度与碎石层平均温度一致时,垂直方向上的热通量达到最小值,此时换热效率最低。因此,在此后的室内碎石层试验设计中,若需减少水平方向上的热量交换带来的影响,应尽量将环境温度设置在碎石层平均温度附近。工程中可通过降低碎石层侧边界温度来增强碎石层内外的热量交换,从而增强碎石层的换热效率。
冻融作用会对土体结构和性质造成影响,在冻融作用下,不同类型土壤的抗剪强度变化呈现出不同规律。本文选取青藏高原具有代表性的黏土,开展不同初始含水率和不同冻融循环次数下的重塑土抗剪强度变化研究,结合压电陶瓷监测试验,探究冻融循环下青藏黏土抗剪强度指标的变化模式和规律。结果表明:在反复冻融作用下,青藏黏土抗剪强度、黏聚力和内摩擦角随含水率增大呈阶梯下降;同一含水率的抗剪强度和黏聚力劣化曲线可以分为快速劣化和稳定两个阶段,内摩擦角呈先减小后增大的变化趋势;冻融循环过程中第1次循环对青藏黏土抗剪强度影响最为明显;采用压电陶瓷监测的应力波信号通过小波包分析可以得到真实信号振幅和应力波能量,二者均能反映青藏黏土抗剪强度的劣化趋势。本文提出的强度劣化指标能很好地反映青藏黏土抗剪强度的劣化程度,这在未来工程实际中有着一定的应用前景。
为研究寒区碳酸盐渍土强度特征,建立冻融循环作用下碳酸盐渍土的修正邓肯-张模型,指导碳酸盐渍土地区的工程建设,利用GDS三轴对取自黑龙江省某地的碳酸盐渍土进行不同围压、不同冻融循环次数下的UU试验。通过曲线拟合,获得了土样强度指标c和φ、模型参数Κ和n、破坏比Rf与冻融循环次数N的数学模型,根据模型理论值与实测值对比分析修正了邓肯-张模型切线模量Ei,从而建立了考虑碳酸盐渍土冻融循环次数N的修正邓肯-张模型。结果表明:碳酸盐渍土经冻融后,其不排水应力-应变曲线呈应变硬化型,且碳酸盐渍土在经历三次冻融循环后,强度降低幅度最大,占整个强度劣化过程的71%以上;冻融后,碳酸盐渍土内摩擦角φ、黏聚力c和参数K逐渐降低,破坏比Rf呈现出先降低后增加,而参数n表现出先增加后减小的趋势。本文提出了冻融循环作用下的修正邓肯-张模型,且模型预测精度高。
冻土动力学特性与冻土区基础设施工程的设计、运营、维护有着密切的联系。冻土的动力学特性是从冻土的滞回曲线中分析得到的,因此全面分析冻土滞回曲线形态特征及其影响因素对明确冻土动力学特性有着重要的意义。冻土滞回曲线的形态特征可用滞回曲线的轴向动应力幅值、不闭合程度、宽窄程度及面积4种参数来表征。冻土的轴向动应力幅值和动应变幅值所围成的骨干曲线反映了冻土等效变形的大小、非线性效应以及冻土的动弹性特性,冻土滞回曲线的不闭合程度反映了冻土的塑性变形性能,冻土滞回曲线的宽窄程度反映了冻土的黏滞特性,冻土滞回曲线的面积反映了冻土的耗能能力。为此,以我国西北季节冻土区粉质黏土为研究对象开展了负温动三轴试验,试验包括4种土温(-0.3 ℃、-1 ℃、-3 ℃、-5 ℃),3种土体含水率(14%、16%、18%),3种围压(0.1 MPa、0.2 MPa、0.3 MPa),以及3种加载频率(1 Hz、2 Hz、4 Hz)交叉组成的10种工况,探究了冻结粉质黏土滞回曲线形态特征及其影响因素。根据试验结果,主要分析了4种描述滞回曲线形态特征参数随着振次、围压、土温、含水率、加载频率以及动应变幅值等因素的变化规律。结果表明,振次对冻结粉质黏土滞回曲线形态特征的影响很小。随着动应变幅值增大,冻结粉质黏土的骨干曲线增长模式逐渐由近似线性模式转为明显的非线性模式;滞回曲线的不闭合程度和面积随动应变幅值的增大呈凹形曲线形态增长;宽窄程度呈凸形曲线形态增长。随着围压、含水率、加载频率的增大以及土温的降低,冻结粉质黏土的骨干曲线变得越陡峭,非线性效应表现得越弱;滞回曲线的不闭合程度、宽窄程度和面积都有不同程度的降低趋势。当围压在0.2 MPa以下、含水率在16%以下、土温在-3 ℃以下时,冻结粉质黏土滞回曲线形态特征随这3种因素的变化更明显。冻结粉质黏土滞回曲线的不闭合程度和宽窄程度受加载频率的影响较为均衡,并没有出现最优影响范围。与冻结粉质黏土滞回曲线的不闭合程度和宽窄程度随各影响因素变化程度相比,滞回曲线的面积受各因素的影响程度相对较小。在4种影响因素中,加载频率和土温对滞回曲线的面积影响较大,含水率的影响次之,围压对滞回曲线面积影响很小。当围压为0.1~0.3 MPa时,冻结粉质黏土的滞回曲线面积-动应变幅值关系曲线随围压变化非常小。
土石混合体-混凝土界面在寒区工程广泛分布,其界面处颗粒骨架受冻融影响塌陷、重组,极易诱发界面强度劣化,是影响寒区工程安全运营的难题。解析冻融影响下界面的强度行为,揭示界面强度与孔隙结构和颗粒骨架的互馈效应,是解决此难题的关键。开展土石混合体-混凝土界面室内直剪试验,借助NMR分层测试及PFC数值模拟,获取界面区抗剪强度参数及孔隙结构演化特性,明晰界面颗粒旋转运移特征;基于分形维数理论,定量评价界面孔隙结构变化,阐释冻融下界面强度劣化机理。即强度骤降阶段界面处土石颗粒受冻胀挤压,其有序排列被打乱,联结、互锁作用弱化,此时孔隙体积、分形维数均增大,界面整体性下降;反翘阶段界面处大块土颗粒分裂重组、骨架塌落,形成块石为核、黏土在外覆盖的包裹体结构,孔隙内部复杂程度降低导致抗剪强度小幅提升;随冻融次数增加,碎石外土颗粒剥落,界面区孔隙体积及颗粒旋转量缓慢增大,界面逐渐脱黏劣化。研究成果对工程构筑物孕灾机制具有指导价值。
冻土本构模型对于分析寒区岩土工程的安全和稳定性至关重要。低温吸力理论可以在物理层面解释温度对冻土力学特性产生影响的机理,且基于该理论所建立的冻土本构模型能实现融土和冻土本构关系的统一,所以近年来低温吸力理论受到越来越多的关注。本文首先通过分析饱和冻土中冰-水界面的表面张力作用,在低温吸力理论的基础上提出了低温粒间吸力的概念,并基于Clapeyron方程建立了低温粒间吸力与温度的关系式以及饱和冻土的有效应力公式。其次,针对多种类型的饱和冻土开展了温控侧限压缩试验和三轴剪切试验,发现低温粒间吸力对冻土力学特性的影响可以归结为其对冻土表观超固结性的影响,即冻土的表观超固结程度随低温粒间吸力的增大而增强。根据试验规律,建立了冻土表观前期固结压力和抗剪强度计算公式。最后,将以上所建立的冻土有效应力公式、前期固结压力公式和强度公式与超固结融土的统一硬化本构理论(即UH模型理论)结合,提出了一个新的冻土热力耦合弹塑性本构模型。使用该模型模拟冻土侧限压缩试验、单轴压缩试验和三轴剪切试验,发现其能合理描述冻土在不同温度和围压下的应力-应变关系,且计算结果与试验数据符合良好。
为了揭示路基冻结过程中地下水和土性对水分迁移规律的影响,针对开放体系和封闭体系的粉质黏土和砂土进行了单向冻结条件下的水分迁移试验。通过土柱上层位置设置碎石层,阻断液态水迁移路径,监测冻结过程中土柱的水热变化,结合土柱冻结深度、冻结速率曲线、含水率分布曲线和补水时程曲线,分析仅水汽补给时对土柱顶部水分聚集和冻结特征的影响。试验结果发现,无论是封闭体系还是开放体系,粉质黏土和砂土土柱都会在冻结区中形成两处水分聚集区:第一水分聚集区为控温板底部,以霜的形式聚集,主要是由土柱顶部土体的水汽迁移并凝华相变形成;第二水分聚集区为冻结区中液态水和气态水共同迁移形成,随着冻结锋面的向下推移,形成不连通孔隙的界面,液态水向0 ℃冰锋线迁移聚集并相变成冰,水汽迁移路径受阻而凝华成冰,致使该处含水率显著增加。相较于封闭体系,开放体系使两处水分聚集区产生更大的水分增量。相比于粉质黏土,砂土介质孔隙较大,在试验时间内水汽补给对水分聚集区的影响更明显,但由于砂土持水能力减弱,水汽补给速率随时间逐渐减小。
冻土强度是冻结法施工及冻土构筑物安全评估中的关键力学参数。然而,冻土强度特性十分复杂,受到环境温度、土体类型、含水量、应变速率等诸多因素的影响,目前还缺少统一的结论。鉴于此,文章对冻土强度的相关试验成果进行了全面回顾,并总结了各因素对冻土强度的影响规律及其作用机理。文献调研结果表明,冻土强度随温度的降低而增大,在工程常见温度范围内二者近似呈线性关系。温度降低所引起的冰强度增加和未冻水含量降低是冻土强度增大的主要原因;当温度低于一定阈值后(如-80 ℃),冻土强度达到最大值并不再变化。非饱和冻土的抗压强度随含水量的增加而增大,原因在于其饱冰度会随着含水量的增加而增加,使得土颗粒与冰晶体的黏结作用增强。对于饱和冻土,含水量的增加会导致冻土密实度降低,其强度随含水量的增加而降低。冻土抗拉强度与抗压强度具有较强的相关性,其压拉比随土体类型的不同存在较大差异,在3~12之间。整体来看,影响冻土强度的因素众多,目前的研究工作多基于单一影响因素,针对多种因素耦合作用、应变速率影响和抗拉强度的研究仍较少,相关工作有待进一步深入。
应变局部化问题是冻土材料破坏的前兆,目前相关研究还较为薄弱。在微极连续体理论框架下,通过考虑温度对相关模型参数的影响,建立了可以模拟冻土应变局部化问题的微极弹塑性数值框架。通过对冻结砂土平面应变试验进行数值分析,发现所建立的数值框架可实现对冻土剪切带产生、发展和形成的全过程模拟,且模拟获得的应力-应变曲线及剪切带破坏形态与试验结果吻合较好;进而探究了温度、端部约束及初始缺陷对剪切带发展规律的影响,发现温度对剪切带倾角和宽度影响较为显著,倾角随温度的降低而逐渐增大,宽度随温度的降低而逐渐加宽,具体剪切带的破坏形态和破坏位置需要综合考虑端部约束及初始缺陷的共同作用;通过对不同有限元网格条件下的冻土应变局部化进行分析,发现粗网格和细网格获得的模拟结果基本一致,说明所建立的数值框架可以基本克服网格依赖性问题。
为探究岩石在冻融循环作用下的裂隙演化特征,利用CT扫描技术,获取不同冻融次数岩石的三维图像,并通过孔隙网络模型参数进行定量表征,研究岩石随冻融循环的裂隙形貌演化过程,以及孔隙结构随冻融次数的变化。结果表明:完整岩石冻融微裂隙主要集中在岩石外围,呈环形片状;而裂隙岩石冻融产生的微裂隙主要存在于预制裂隙周围和岩桥区域,且分布范围更广。冻融循环作用后,完整与裂隙岩石的PNM中各类孔隙变化特征不一。前者孔隙数量减少,小孔占比逐渐减小,中大孔占比增大,孔隙分布位置较为均匀;后者孔隙数量增加,但中孔占比减小,大孔占比迅速增大且主要集中在岩心上下端面。在冻融过程中,岩石各类孔喉结构主要参数均呈增长趋势,裂隙岩石孔喉网络发育更剧烈,冻融结束后其最大喉道半径增长幅度约为完整岩石的2倍。在冻融前期,冻融循环作用主要促进岩石孔隙相互连通,后期主要表现为孔喉加速扩张,且裂隙岩石孔隙率增长速率均大于完整岩石。研究结果可为寒区岩体破坏机理及灾害防治提供参考。
山地冰川因具有高反照率、冰川风、逆温层及高值降水等特征而形成了独有的局地微气候,尤其是作为高值降水中心,对径流变化具有重要影响。本文基于祁连山中段北坡摆浪河21号冰川末端(海拔4 350 m)2020年9月2日—2021年8月28日的气象观测资料,开展了微气象特征分析,与临近不同海拔、下垫面的同期降水以及祁连山典型冰川区降水进行了比较,并对最大降水事件过程开展环流成因分析。研究发现:摆浪河21号冰川区气温超过0 ℃的天数有84 d,集中在5—9月;冰川风盛行,不同于其他冰川区的山谷风循环;天气主要以多云为主;入射与反射短波辐射月最大值分别出现在5月和4月。摆浪河21号冰川降水主要集中在4—8月,降水频次和强度均随着云量的增加而增加。观测年最大降水事件(2021年7月25—27日)属于局地对流降水,中高纬西北-东南向水汽输送为降水区提供了大量水汽;低层辐合和高层辐散、层结不稳定造成了强烈的暖空气上升,加之降水区位于槽后脊前不断有冷空气输入,冷暖交汇促使降水发生。
积雪作为大气固态降水过程的主要形式,是反映大气化学成分的良好媒介。然而,积雪样品在采集后的不同储存和处理方式可能会影响其原始信息的获取。本研究采集青藏高原腹地纳木错地区的积雪样品,使用不同前处理方法探究其对主要阳离子浓度的影响。结果表明:积雪样品融化后第3天及更长时间,使用不过滤的前处理方式会使测试结果明显偏离实际值。过滤操作能规避碳酸盐矿物溶解对Ca2+、Mg2+和Sr2+的富集,减少对大陆粉尘来源贡献的高估;对K+和Na+的影响则主要通过黏土矿物与水体间的离子交换实现,使得前者浓度降低、后者升高,从而可能低估生物质燃烧源的贡献、高估海盐输入。就不同孔径的滤膜而言,0.20 μm滤膜过滤效果优于0.45 μm滤膜。过滤后溶液的大部分离子在冷藏或常温存储的前10天变化较小(<10%)。综上,本文建议积雪样品在融化后当天使用孔径为0.20 μm的滤膜过滤,滤后溶液存储方式可选择冷藏或常温保存,并在10天或更短时间内完成测试。
积雪消融是气候变化和水资源管理的重要影响因素。本文联合Sentinel-1与Sentinel-2时序数据,提出了一种联合SAR与光学遥感数据的积雪消融识别方法。以祁连山区八宝河流域为研究区,利用ESTARFM时空数据融合模型得到与Sentinel-1同时相的Sentinel-2模拟数据,提取流域积雪覆盖范围;基于SAR多时相多极化变化检测算法提取湿雪;结合光学遥感与DEM数据对SAR提取的湿雪进行校正,最终获得八宝河流域消融期内18个时相干湿雪分布。利用GF-2影像和消融期的Sentinel-2影像分别对积雪和湿雪面积进行精度验证,结果表明该方法能够快速识别积雪消融时空变化,总体分类精度OA高达99%,Kappa系数高达0.86。同时,利用同一天早晚过境相差约12小时的升、降轨SAR数据对比分析了积雪消融变化特征。根据实验数据集分析表明,八宝河流域内干湿雪分布随时间剧烈变化,消融初期湿雪主要集中在河谷低海拔区域且消融速度快,随着气温的升高积雪消融至高海拔区域,而干雪主要分布在四周高海拔山区。且由于过境时间的差异,在整个消融期,降轨数据提取的湿雪面积小于升轨数据提取的湿雪面积。
MODIS V006版本数据仅提供了归一化积雪指数(NDSI),而用户往往关心的是直观的积雪分类,包括积雪范围或积雪覆盖率。美国国家冰雪数据中心推荐全球积雪范围最佳的NDSI阈值为0.4,但是青藏高原地形复杂多样,积雪斑块化特征明显,单一阈值并不能精确地判识不同下垫面上的积雪。青藏高原被称为地球的第三极,是中国三大稳定积雪区之一,蕴藏了大量的淡水资源。随着全球气候变暖,青藏高原地区积雪融化时间提前,冰川融水增加,影响河流水量,造成洪涝灾害,进而影响人类正常生产生活,因此通过确定不同下垫面阈值,改善传统阈值的积雪高估低估现象,提高积雪识别精度,进而更准确地探究青藏高原积雪状况,显得尤为迫切。本文以青藏高原为研究对象,首先生成MODIS逐日无云NDSI序列并进行验证;其次对应站点雪深数据与NDSI序列,证实在下垫面为林地和非林地的区域,去云NDSI序列与站点雪深均有良好的对应关系,确定不同下垫面最优阈值范围;最后在最优阈值范围内通过混淆矩阵确定最优阈值。计算得出,林地NDSI=0.03时,总体精度最高为94.02%,在该NDSI之下,高估误差OE和低估误差UE分别为1.21%和4.60%;非林地NDSI=0.26时,总体精度OA最高为94.27%,在该NDSI之下,高估误差OE和低估误差UE分别为0.51%和5.03%。因此选取优化后林地阈值为NDSI=0.03,非林地阈值为NDSI=0.26。为避免地面常规观测资料尺度上的局限性,本文采用高精度的Landsat 8 OLI卫星数据识别结果,作为“真值”对优化后阈值的判别结果进行“像元—像元”级别的验证。在定量验证中,优化后 NDSI阈值对MOD10A1 V006积雪判别结果的总体精度OA为84.21%,高估误差OE为5.33%,低估误差UE为10.46%;传统阈值对MOD10A1 V006积雪判别结果的总体精度OA为82.86%,高估误差OE为1.48%,低估误差UE为15.66%。可以看出在定量验证中,优化后阈值的积雪判别精度更高。同时在定性验证中,积雪大面积集中的区域,新的阈值与传统阈值提取效果均相对较好;积雪相对分散破碎的区域,优化后阈值能提取出大量积雪,传统阈值则不能。这表明考虑不同土地覆盖类型下的NDSI阈值优化可以有效地提高青藏高原积雪判别精度,为NDSI在积雪识别中的应用提供有力的支撑,有助于更准确地了解该地区积雪分布状况。
国土面积2/3为寒区的中国拥有广泛的路网、机场和城镇等人工地面系统,除冰雪技术研究与发展对道路安全与保障起到重要作用。为全面了解国内道路除冰雪领域研究现状和趋势,通过中外常用数据库核心集获取2010年以来我国学者在道路除冰雪领域的研究文献,借助CiteSpace软件对获取的文献信息中的发文量、文献来源和关键词进行分析。结果显示,我国道路工程领域对除冰雪技术愈发重视,自2010年以来我国在该领域的发文量逐年增加。同时发现,我国道路除冰雪领域学者在中文期刊的发文量占比在近年来则有所降低,更倾向于将研究成果发表于英文期刊。结果表明,长安大学、哈尔滨工业大学等高校是该领域主要的研究机构,且国内科研院所在国际合作方面还有待提升。在研究内容方面,我国道路工程学者对融雪剂的研究在减少,更加关注流体加热路面、电热路面等新型除冰雪技术。通过对2010年以来的高共被引文献分析,电热除冰雪技术、微波加热除冰雪技术、融雪抑冰路面等技术优势突出,更易引起学者关注,但在工程应用方面存在不足。通过文献综述发现,我国道路工程除冰雪领域急需解决上述新型除冰雪技术的能耗、工艺、成本、融雪效果及路用性能平衡等几个关键方面的问题。
从20世纪80年代初第一次阅读拉肯布鲁克博士的论文开始,我一直被他的冻土物理热学的研究方法吸引着,也在很长一段时间内影响着我对科研课题的选择。我曾把他的热学方程应用到自己的研究。1993年在北京举行的第六届国际冻土大会期间,我与拉肯布鲁克博士相见,并在会后一起参加了青藏高原调研旅行。尽管他个头高大,但说话却很温和,是一位非常耐心的良师益友。我曾通过信件联系他,并向他请求寄送其发表的文章给我,另外,他也曾作为专家为我的文章审稿。拉肯布鲁克博士的去世是地球物理学和冻土学界的一大损失,他是冻土物理热学研究的巨匠和先驱。他在美国地质调查局的专业研究生涯长达60多年,他的研究深入到地球物理、地热热流、多年冻土、地貌,以及大地构造和地震学,并在这些领域作出了杰出的贡献。他对深层冻土温度剖面的测量和研究使他比较早地提供了地球变暖及全球变化的证据。他对热油管道和多年冻土相互作用的研究避免了掩埋热油管道对冻土环境灾难性的破坏。拉肯布鲁克博士是一位多产的科学家,出版了大量的著作,他的许多文章已成为冻土学研究的基准和经典,这是他留给后人的宝贵遗产和财富。在我们缅怀这位大师的同时,也为他对地球物理及冻土学所作出的巨大贡献致以敬意。拉肯布鲁克博士不仅是一位科学家,他还是一位城市规划能手,他曾志愿担任所居城市市政规划委员会主任长达10年之久。为了纪念拉肯布鲁克博士,市政府最近决定用他的名字来命名一条人行道。