An in-situ method for assessing soil aggregate stability in burned landscapes
-
摘要: 森林火灾后因火烧迹地土壤斥水性,导致坡面径流和土壤可蚀性增强,提高了火后泥石流易发性,而土壤团聚体稳定性是影响土壤入渗能力和侵蚀敏感性的关键指标。目前常用于火烧迹地土壤团聚体稳定性测定的水滴冲击测定方法(counting the number of water drop impacts,CND),不适用于原位测定且耗时较长(滴定一组团聚体需要数小时)。因此文章提出一种基于冲击振荡破坏效应的团聚体稳定性测定方法(shock and vibration damage method,SVD)。充分考虑容重、有机质含量和斥水性对土壤团聚体稳定性的影响,通过室内火烧模拟试验,制备了13种类型的土壤团聚体。采用自制的试验仪器进行SVD法正交试验测定土壤团聚体质量损失率,并与传统CND法测得的破坏团聚体的水滴数量进行对比。结果表明:SVD法的测定MT-6方案(冲击高度1 m、容器容水量40%、冲击5次、测定团聚体20颗)与CND法的测定结果具有很强的一致性(Kendall系数=0.797)和相关性(R2=0.634),测定时间较短(测定一组团聚体约5 min),且测定结果区分度较好(约62%的团聚体MLR位于区分度良好的40%~60%区间),将其作为SVD法的最优测定方案。此外,SVD法试验装置结构简单、便携易拆卸,可用于原位快速且定量地区分火烧迹地不同火烈度下土壤团聚体稳定性水平,对火烧迹地土壤侵蚀、水土流失治理以及火后泥石流起动机理研究具有重要指导意义。Abstract: Due to soil repellency in burned areas, slope runoff and soil erodibility escalates following forest fires, increasing the vulnerability to post-fire debris flows. Soil aggregate stability is a critical determinant of soil infiltration capacity and erosion susceptibility. The prevalent method of assessing soil aggregate stability in burned areas, the counting the number of water drop impacts (CND) method, is time-intensive and impractical for in-situ measurements. In response, this study introduces a novel technique based on the shock and vibration damage (SVD) effect for evaluating soil aggregate stability in burned areas. Thirteen distinct soil aggregate types were meticulously prepared for indoor simulated fire testing, with due consideration to factors such as bulk weight, organic matter content, and water repellency, which influence stability of soil aggregates. Employing a custom-built test apparatus, the mass loss rate (MLR) of soil aggregates was determined through orthogonal experiments using the SVD method and compared against the standard CND technique's quantification of water droplet-induced aggregate destruction. The findings demonstrated that SVD method, employing Test Scheme 6 (testing 20 aggregates, 1-meter impact height, 40% water content, and five impacts), exhibits excellent agreement (Kendall coefficient = 0.797) and correlation (R2 = 0.634) with CND method outcomes. This testing scheme, characterized by rapid determination and effective discrimination, is identified as the optimal testing approach. The SVD testing apparatus is straightforward, portable, and easily disassembled, rendering it suitable for on-site use. It can be used to distinguish the stability level of soil aggregates swiftly and quantitatively under various fire intensities in burned areas in situ, which is an important guiding significance for the study of soil erosion, erosion control, and post-fire debris flow initiation mechanism in burned areas.
-
0. 引言
西南山地存在众多的不同成因的堆积体斜坡,坡体上常聚居村落、交通和水利水电等基础设施。在多种因素的作用下,堆积体斜坡常易形成大规模的滑坡,威胁人民生命财产安全。这类滑坡地处高山峡谷地貌,地质环境条件复杂多变,诱发滑坡变形的因素众多,其中地震和强降雨是主要的诱发因素,特别是两因素常相互作用防治难度大[1-2]。
炉霍县旦都乡马居滑坡位于鲜水河断裂带附近,受近年来地震活动[3]和降雨影响,从2015年开始明显变形,且逐年加剧,现已演化为一大型滑坡。滑坡的破坏模式和滑坡机理在区域内具有典型性,本次拟对此滑坡的变形破坏机理进行研究,以期理清这类滑坡变形发展的过程和各因素相互作用的影响,进一步优化防治方案,评价方案的有效性和安全性,为地质灾害防治和风险管控提供依据。
1. 研究区地质环境条件
炉霍县地处青藏高原东南边缘,高山峡谷地貌区,海拔高度3300~3600 m,区域内河流为鲜水河支流。平均年降雨量700 ~920 mm。近10年降雨量较平稳。区域内地质构造复杂,新近构造活动强烈(图1)。历史上遭受多次地震的影响,1973年炉霍县鲜水河断裂带发生7.6级地震,2014年11月22日鲜水河断裂带的康定地区发生6.3级地震及多次余震。研究区基本地震动峰值加速度为0.3 g,反应谱特征周期为0.4 s[4]。
图 1 研究区地震断裂及历史地震分布图[4]Figure 1. Distribution map of seismic faults and historical earthquakes in the study area研究区地层主要为三叠系上统雅江组上段(
)含粉砂炭质绢云板岩和变质长石石英砂岩、第四系残坡积层( )、滑坡堆积层( )。2. 滑坡基本特征
2.1 滑坡基本概况
滑坡整体呈圈椅状,滑坡后缘以寺庙(马居庙)外侧围墙前拉张裂缝为界,前缘以俄觉龙日沟沟床为界,左侧和右侧以水沟为界,滑坡面积4.6×104 m2,厚10~13 m,体积60×104 m3,属中型土质滑坡。滑坡主滑方向为100°,整体坡度约30°,呈前陡后缓的趋势,前缘坡度约35°~40°,后缘坡度约20°~25°。
2.2 变形特征
根据滑坡的变形破坏特征将该滑坡划分为2个变形区,其中Ⅰ区为寺庙及下方坡体,Ⅱ区为寺庙外北侧区域。变形Ⅰ区主要变形特征有:前缘沟床侵蚀垮塌下错台坎,横向拉张裂缝,后缘拉裂缝、寺庙房屋裂缝等。变形Ⅱ区主要变形特征有:后缘拉裂缝、地面开裂、电杆歪斜、前缘多级下错陡坎群等。
调查表明,该滑坡区在2015年以前未出现可见变形,受2014年11月康定地震及余震影响,该区震感明显,2015年雨季时滑坡前缘陆续变形,以后每年雨季滑坡体变形呈现出了逐年加剧的趋势,滑坡体后缘也陆续出现了变形,2020年发生了一次短历时强降雨,雨强约为50年一遇降雨,诱发了俄觉龙日沟泥石流,对坡脚冲刷淘蚀,滑坡前缘出现了大规模的滑坡变形,并进一步牵引后缘变形加剧。
2.3 坡体结构特征
滑体主要由第四系全新统滑坡堆积层组成。厚度为10~13 m,主要成份为含黏碎石土,所充填黏性土褐黄色,稍湿,无摇振反应,干强度较高,切面较光滑,碎石成分主要为变质砂质、板岩,常见钙质胶结,粒径一般2~10 cm,含量大于50%。滑带位于滑体或基覆界面,含水量较大,地下水活动痕迹明显。滑床为三叠系上统雅江组上段含粉砂炭质绢云板岩(T3y3)和变质长石石英砂岩互层(表1)。滑坡区工程地质平面图见图2,滑坡体结构特征见图3。
表 1 滑坡区主要岩土参数表Table 1. Geotechnical parameters of landslide area岩土体 体积模量/MPa 剪变模量
/MPa天然状态 饱和状态 重度/(t·m−3) 内聚力/kPa 内摩擦角/(°) 重度/(t·m−3) 内聚力/kPa 内摩擦角/(°) 滑坡体 15.60 5.18 2.01 24.3 30.8 2.05 22.9 27.8 强风化 15.20 7.80 22 40 35 - - - 中风化基岩 20.00 12.00 27 1900 42.8 - - - 3. 滑坡变形破坏机理分析
马居滑坡是位于地震区的一堆积体滑坡,勘查表明,该滑坡坡度较陡,滑坡体结构为含黏碎石土层。下覆基岩为含粉砂炭质绢云板岩具有较好的隔水性,基覆界面地下水活动较强。受2014年地震作用滑坡体结构完整性和稳定性变差,随后受降雨和坡脚下切作用坡体变形持续加大,滑坡体处于不稳定—欠稳定状态。
该滑坡变形过程及所受影响因素明确,但总体来说,该滑坡地质环境条件复杂多变,诱发因素和作用方式多变,通过对滑坡变形特征及变形历史分析,基本理清了近期滑坡变形的时空关系,为更好地查清该滑坡变形破坏机理,更好地指导滑坡防治方案的制定,这里对滑坡三个重要的影响因素进行数值分析。选取的分析剖面为滑坡主滑方向的2-2′剖面,模型水平方向左右采用位移约束,模型的底部采用竖向位移约束。共划分14070个结点,6811个单元。参数见表1。
3.1 地震作用下坡体变形机理分析
对马居滑坡变形历史调查表明,2014年发生在鲜水河断裂带上的康定地震,对该滑坡稳定性影响大,为了进一步查明地震对该滑坡具体影响程度,这里选取该地区多遇地震加速度及特征周期人工合成地震波进行时程分析[5]。人工合成的地震波见图4。
时程分析获得了斜坡在地震作用下坡体的塑性区分布图及最大剪应变增量云图[6](图5、图6)。
模拟结果表明地震作用对斜坡稳定性影响大,斜坡体内塑性区基本上成片相连,特别是斜坡顶部出现了大面积的塑性区,该区域也是地表裂缝发现最多的区域,剪应变云图表明在基覆界面处剪应变增量最大,且出现了两个较大的区域,特别是沟道向上约10~30 m的位置,勘查表明该区裂缝发育密集。在地震荷载作用下,坡体稳定性下降,斜坡土体结构遭到破坏,地震后不久,在坡体缓坡平台前缘部分相继出现开裂,该点随即纳入监控。
3.2 坡脚沟道下切诱发坡体变形破坏分析
2020年7月该区域暴发了短历时强降雨,俄觉龙日沟暴发了50年一遇泥石流,坡脚沟道淘蚀严重,滑坡体前缘变形加剧。为了查明沟道下切对滑坡稳定性影响,以本次勘查沟道内的基覆界面作为最大冲刷深度进行数值分析,采用参数折减法进行稳定系数计算。获得下切后稳定系数为1.00,最大剪应变区域延伸至沟道底部,表明沟道下切后牵引作用明显,直接导致坡体变形加剧(图7)。
3.3 不同降雨频率条件下稳定性分析
该滑坡不但受到短历时强降雨冲刷坡脚的影响,变形历史表明,长历时降雨也对滑坡的稳定性影响较大。研究表明[7-19],滑坡体采用Mein-Larson降雨入渗模型能较好地模拟出降雨对坡体含水率影响,可计算出不同时刻坡体湿润锋深度,降雨结束后,再采用渗透模型计算各时刻地下水的分布特征,对于特征时间结点采用强度折减法计算出各时刻滑坡稳定系数,以此评价降雨全历时滑坡的稳定性变化。滑坡体的饱和渗透系数为0.42 m/h,100年一遇降雨强度为0.034 m/h,50年一遇降雨强度为0.031 m/h,20年一遇降雨强度为0.026 m/h,10年一遇降雨强度为0.023 m/h,常年降雨强度为0.015 m/h,滑坡体的饱和渗透系数大于小时降雨强度,则垂直湿润锋深度
可采用以下公式计算:(1) 式中:M——饱和含水率与初始含水率的差值(本次取0.1);
p——降雨强度m/h;
t——降雨历时/h。
本次降雨历时综合近年来该地区长历时强降雨并综合类似文献[6]后综合确定计算时长为10 d。强降雨入渗坡体后采用地下水渗透模型进行模拟分析,地表及沟道表明设置为透水边界,其余设置为不透水边界条件。
分析结果可以看出强降雨的雨强越大对滑坡的稳定性影响越大,强降雨在降雨过程中坡体浅表饱和,湿润锋下移,但是滑坡稳定下降不明显,随着湿润锋向坡体转移过中滑坡体的稳定性有较大幅度的下降,随后随着地下水沿着隔水界面活动时,滑坡稳定性长期处于不稳定状态,待滑坡体内的地下水基本被疏干后滑坡体整体稳定性逐渐恢复到初始状态。说明滑坡风险管控中不但要关注降雨时刻会诱发滑坡灾害,在降雨过后较长时间(可能5 d左右)内照样会出现滑坡灾害(图8)。
本次对于降雨入渗研究表明厚层堆积体滑坡的稳定性较降雨条件下浅表层滑坡稳定性的变化在降雨阶段有所不同,主要体现在浅表层无限边坡和岩质边坡[20]的稳定性在降雨期间下降很快[6],而厚层堆积体滑坡的整体稳定性虽有向下的趋势,但是整体稳定性下降不大(图9)。
3.4 滑坡变形破坏机理小结
文中对马居滑培的机理分析遵循该滑坡变形历史过程进行模拟,力求获得该滑坡地震、坡脚下切、强降雨入渗等因素的影响。分析结果表明:马居滑坡先是地震作用对坡体损伤作用明显,促使坡体裂隙发育,位于基覆界面的滑面基本形成,滑坡体整体处于欠稳定状态,后期随着降雨入渗的影响,滑坡体处于缓慢变形,2020年的短历时强降雨诱发泥石流冲刷下切坡脚,诱发前缘大规模滑坡,牵引作用使得滑面进一步贯通,变形加剧,坡体整体稳定性进一步下降。模拟结果很好地拟合了滑坡变形破坏过程,进一步揭示了滑坡变形破坏机理,这为工程防治和风险管控打下了基础。
4. 滑坡防治方案分析
4.1 防治方案
马居滑坡坡体上常住人口多,周边斜坡失稳后将直接威胁寺庙安全,其后果严重,属大型地质灾害。该滑坡地质环境条件复杂,滑坡诱发因素多,地震、降雨等叠加后易形成特大型灾害。
滑坡体及周边地形坡度大,从寺庙至沟道斜坡基本上不具备设置抗滑支挡的工程施工条件,但该段斜坡稳定性最差,在天然状态下其稳定系数介于1.03~1.05,沟道下切后坡体的稳定性将下降到1.0附近,若继续冲刷将整体失稳。寺庙建在缓坡平台上,寺庙前缘有一较缓坡平台,可作为抗滑桩施工场地,若在此设置抗滑桩将直接针寺庙建筑物进行保护。
所以,在综合考虑滑坡变形破坏机理和施工条件的基础上,确定了本滑坡治理方案:在滑坡体中部(寺庙前缘)抗滑支挡,坡脚沟道固床防止下切。
4.2 坡脚防护
坡脚沟道下切持续牵引坡体整体变形,坡脚沟道的防护措施是必要的。坡脚防护有两种方案,一种是沟道固床,保持现有状态,沟道不继续下切,前缘坡体在自然条件下逐步稳定。另外一种方案就是回淤压脚,提高滑坡前缘稳定性。后一种方案属于主动防治方案,这里选取回淤高度2 m、4 m、6 m、8 m、10 m进行回淤压脚效果分析,分析采用参数折减法。
分析结果表明,不同回淤高度的压脚提升滑坡稳定性差距较大,回淤高度8 m以内对滑坡前缘稳定性基本上无明显改善,只有当回淤高度超过8 m时才有较大幅度的提高。通过对滑坡前缘非稳定状态下的最大剪应变云图分析,滑坡前缘最不利滑面的剪出口距沟底的高度大于8 m,若回淤高度小于8 m,对稳定性的贡献基本上很小,所以稳定系数一直处于徘徊于1.03附近,若要大幅度提高滑坡整体稳定性,回淤压脚的高度至少要大于10 m,考虑到该滑坡体坡脚沟道纵坡较陡、沟道较长,必须采用多道高坝才能达到治理的效果(图10)。
4.3 极端条件下防治方案安全性评价
滑坡区地处中国最活跃的地震断裂带上,发生地震可能性大,对在地震工况及强降雨条件下防治方案的安全性进行评价是必要的[21-22]。本次考虑到在极端条件下,如地震和长历时强降雨来袭后,滑坡前缘失稳,或者在坡脚淘蚀作用长期作用下滑坡前缘坡体累进破坏导致前缘整体失稳,使抗滑桩前缘失去支挡,在这种条件下对上部抗滑桩方案的安全性进行评价。滑坡区基本地震动峰值加速度为0.3 g,多遇地震动峰值加速度为0.1 g,罕遇地震0.57g[4](通过数值分析获得了各工况条件下滑坡稳定性分区云图图11-图15)。
数值分析表明前缘失去支挡后,抗滑桩方案对上部保护对象能有效的保护,稳定性差的部分就是上部滑体及后缘陡坡部分,天然工况及三种地震工况下稳定系数依次降低,分别为1.4,1.3,1.2,1.0。数值分析结果说明在极端情况发生后,上部滑坡不会随着前缘失去支挡后随即破坏且还具有一定安全余度,但是该斜坡坡度陡,易在余震和强降雨作用下再次发生变形破坏,故在极端情况发生后,抗滑桩方案能够为斜坡上的居民赢得撤离时间,提高抗风险的能力。
5. 结论
通过对马居堆积体滑坡变形破坏机理及防治方案研究,得到如下结论:
(1)马居滑坡受地震作用和降雨作用双重影响下,地震作用下坡体结构损伤严重,滑坡体稳定性下降,后期降雨作用下滑坡变形进一步加剧,滑坡目前处于欠稳定和不稳定状态。
(2)对降雨入渗坡体的研究表明,短历时降雨对滑坡的稳定性影响较小,长历时降雨在土体的湿润锋到达一定深度后滑坡的稳定性才有较大幅度的下降,若降雨历时较短,对滑坡的稳定性影响有限,这相对于无限斜坡整体稳定性要好。
(3)短历时强降雨能形成较强的地表径流,形成洪水或泥石流冲刷和淘蚀滑坡坡脚,对滑坡的稳定性的影响较大,易形成牵引式滑坡,进而堵塞沟道形成泥石流的堵溃点。
(4)推荐采用抗滑桩加固床防治方案,固床回淤高度应综合考虑滑坡剪出口位置,抗滑支档方案在极端情况下有一定抗风险能力。
-
图 3 SVD法与CND法测试结果
注:①(a)为强能组,对应方案MT-4、9、10、15、18、19、20、22、24;(b)为中能组,对应方案MT-1、6、12、13、21、23;(c)为低能组,对应方案MT-2、3、5、7、8、11、14、16、17、25。②高温亲水:ST-1~ST-3,轻度斥水:ST-7~ST-9,强烈斥水:ST-10~ST-11,严重斥水:ST-12,极端斥水:ST-13。③同一土壤层不同小写字母代表显著性差异(P < 0.05);同一斥水性不同大写字母代表显著性差异(P < 0.05),相同字母表示差异不显著,误差条表示测量结果平均值的标准误差。
Figure 3. Test results of the SVD method and CND method
表 1 供试土壤主要物理化学性质
Table 1 Key physical and chemical properties of the studied soils
指标 高温亲水组 天然对照组 轻度斥水组 强烈−极端斥水组 ST-1 ST-2 ST-3 ST-4 ST-5 ST-6 ST-7 ST-8 ST-9 ST-10 ST-11 ST-12 ST-13 温度/(°C) 600 600 600 常温 常温 常温 200 200 400 200 400 200 400 时间/min 120 120 120 \ \ \ 5 60 5 30 5 60 5 深度/cm 0~2 2~4 4~6 0~2 2~4 4~6 0~2 2~4 4~6 0~2 2~4 0~2 0~2 斥水性 亲水 亲水 亲水 强烈
斥水轻度
斥水亲水 轻度
斥水轻度
斥水轻度
斥水强烈
斥水强烈
斥水严重
斥水极端
斥水干容重/(g·cm−3) 0.95 1.04 1.16 0.95 1.04 1.16 0.95 1.04 1.16 0.95 1.04 0.95 0.95 初始有机质含量/% 14.0 9.48 6.93 14.0 9.48 6.93 14.0 9.48 6.93 14.0 9.48 14.0 14.0 注:亲水(WDPT≤5 s),轻度斥水(WDPT:5~60 s),强烈斥水(WDPT:60~600 s),严重斥水(WDPT:600~3600 s),极端斥水(WDPT:>3600 s)[23],ST-1~13代表团聚体类型编号。 表 2 团聚体稳定性测定正交试验设计表
Table 2 Table for orthogonal test design for determining aggregation stability
测定方法编号 A/% B/cm C/次 D/颗 测定方法编号 A/% B/cm C/次 D/颗 MT-01 40 80 10 10 MT-14 80 100 1 10 MT-02 100 50 3 10 MT-15 20 100 10 25 MT-03 80 50 5 15 MT-16 20 50 1 5 MT-04 20 80 15 15 MT-17 60 80 5 5 MT-05 60 50 10 20 MT-18 100 200 5 25 MT-06 40 100 5 20 MT-19 20 150 5 10 MT-07 80 80 3 25 MT-20 60 200 15 10 MT-08 100 80 1 20 MT-21 40 200 1 15 MT-09 100 150 10 15 MT-22 80 200 10 5 MT-10 80 150 15 20 MT-23 100 100 15 5 MT-11 60 150 1 25 MT-24 20 200 3 20 MT-12 40 150 3 5 MT-25 60 100 3 15 MT-13 40 50 15 25 注:影响因素A是容器容水量,B是冲击高度,C是冲击次数,D单次测试所需土壤团聚体颗数,MT-1~25代表SVD法试验方案。 表 3 SVD法各方案与CND法的一致性检验结果
Table 3 Consistency test results between various schemes of the SVD method and the CND method
方案 Kendall系数 P 方案 Kendall系数 P 方案 Kendall系数 P 方案 Kendall系数 P MT-1 0.827 0.071 MT-8 0.687 0.170 MT-15 0.791 0.089 MT-22 0.726 0.134 MT-2 0.299 0.845 MT-9 0.709 0.149 MT-16 0.918 0.037 MT-23 0.684 0.173 MT-3 0.495 0.456 MT-10 0.723 0.137 MT-17 0.581 0.304 MT-24 0.679 0.178 MT-4 0.654 0.206 MT-11 0.890 0.045 MT-18 0.874 0.051 MT-25 0.604 0.270 MT-5 0.761 0.108 MT-12 0.874 0.051 MT-19 0.659 0.199 MT-6、
MT-11、
MT-160.797 0.000 MT-6 0.904 0.041 MT-13 0.780 0.095 MT-20 0.815 0.076 MT-7 0.728 0.133 MT-14 0.659 0.199 MT-21 0.794 0.087 注:表中P为显著性,P值小于0.05,说明具有显著一致性;Kendall系数值代表一致性程度(较差:<0.2,一般:0.2~0.4,中等:0.4~0.6,较强:0.6~0.8,很强:0.8~1.0)[25]。 -
[1] 胡卸文,金涛,殷万清,等. 西昌市经久乡森林火灾火烧区特点及火后泥石流易发性评价[J]. 工程地质学报,2020,28(4):762 − 771. [HU Xiewen,JIN Tao,YIN Wanqing,et al. The characteristics of forest fire burned area and susceptibility assessment of post-fire debris flow in jingjiu township,xichang City[J]. Journal of Engineering Geology,2020,28(4):762 − 771. (in Chinese with English abstract) DOI: 10.13544/j.cnki.jeg.2020-224 HU Xiewen, JIN Tao, YIN Wanqing, et al . The characteristics of forest fire burned area and susceptibility assessment of post-fire debris flow in jingjiu township, xichang City[J]. Journal of Engineering Geology,2020 ,28 (4 ):762 −771 . (in Chinese with English abstract) DOI: 10.13544/j.cnki.jeg.2020-224[2] BADÍA-VILLAS D,GONZÁLEZ-PÉREZ J A,AZNAR J M,et al. Changes in water repellency,aggregation and organic matter of a mollic horizon burned in laboratory:Soil depth affected by fire[J]. Geoderma,2014,213:400 − 407. DOI: 10.1016/j.geoderma.2013.08.038
[3] VARELA M E,BENITO E,KEIZER J J. Effects of wildfire and laboratory heating on soil aggregate stability of pine forests in Galicia:The role of lithology,soil organic matter content and water repellency[J]. CATENA,2010,83(2/3):127 − 134.
[4] ZAVALA L M,GRANGED A J P,JORDÁN A,et al. Effect of burning temperature on water repellency and aggregate stability in forest soils under laboratory conditions[J]. Geoderma,2010,158(3/4):366 − 374.
[5] JORDÁN A,ZAVALA L M,MATAIX-SOLERA J,et al. Effect of fire severity on water repellency and aggregate stability on Mexican volcanic soils[J]. CATENA,2011,84(3):136 − 147. DOI: 10.1016/j.catena.2010.10.007
[6] WIETING C,EBEL B A,SINGHA K. Quantifying the effects of wildfire on changes in soil properties by surface burning of soils from the Boulder Creek Critical Zone Observatory[J]. Journal of Hydrology:Regional Studies,2017,13:43 − 57.
[7] JIMÉNEZ-PINILLA P,MATAIX-SOLERA J,ARCENEGUI V,et al. Advances in the knowledge of how heating can affect aggregate stability in Mediterranean soils:a XDR and SEM-EDX approach[J]. CATENA,2016,147:315 − 324. DOI: 10.1016/j.catena.2016.07.036
[8] 胡卸文,王严,杨瀛. 火后泥石流成灾特点及研究现状[J]. 工程地质学报,2018,26(6):1562 − 1573. [HU Xiewen,WANG Yan,YANG Ying. Research actuality and evolution mechanism of post-fire debris flow[J]. Journal of Engineering Geology,2018,26(6):1562 − 1573. (in Chinese with English abstract) DOI: 10.13544/j.cnki.jeg.2018-073 HU Xiewen, WANG Yan, YANG Ying . Research actuality and evolution mechanism of post-fire debris flow[J]. Journal of Engineering Geology,2018 ,26 (6 ):1562 −1573 . (in Chinese with English abstract) DOI: 10.13544/j.cnki.jeg.2018-073[9] 郭军玲,王虹艳,卢升高. 亚热带土壤团聚体测定方法的比较研究[J]. 土壤通报,2010,41(3):542 − 546. [GUO Junling,WANG Hongyan,LU Shenggao. Comparative studies on measurement methods for aggregate stability of subtropical soils[J]. Chinese Journal of Soil Science,2010,41(3):542 − 546. (in Chinese with English abstract) DOI: 10.19336/j.cnki.trtb.2010.03.007 GUO Junling, WANG Hongyan, LU Shenggao . Comparative studies on measurement methods for aggregate stability of subtropical soils[J]. Chinese Journal of Soil Science,2010 ,41 (3 ):542 −546 . (in Chinese with English abstract) DOI: 10.19336/j.cnki.trtb.2010.03.007[10] SCHOMAKERS, J., et al . Measurement of soil aggregate stability using low intesity ultrasonic vibration[J]. Spanish Journal of Soil Science,2011 ,1 (1 ):8 −19 .SCHOMAKERS, J., et al. Measurement of soil aggregate stability using low intesity ultrasonic vibration[J]. Spanish Journal of Soil Science,2011,1(1):8 − 19.
[11] AKSAKAL E L,ANGIN I,SARI S. A new approach for calculating aggregate stability:mean weight aggregate stability (MWAS)[J]. CATENA,2020,194:104708. DOI: 10.1016/j.catena.2020.104708
[12] AMÉZKETA E. Soil aggregate stability:a review[J]. Journal of Sustainable Agriculture,1999,14(2/3):83 − 151.
[13] 黄悦等,基于高能水分特性法的土壤团聚体结构稳定性研究进展[J]. 水土保持研究,2022 :1 − 8. [HUANG Yue,et al. Progress in Research on Structural Stability of Soil Aggregates Based on High Energy Moisture Characteristic Method[J]. Research of Soil and Water Conservation,2022 :1 − 8. (in Chinese with English abstract) HUANG Yue, et al. Progress in Research on Structural Stability of Soil Aggregates Based on High Energy Moisture Characteristic Method[J]. Research of Soil and Water Conservation, 2022 : 1 − 8. (in Chinese with English abstract)
[14] THOMAZ E L. Effects of fire on the aggregate stability of clayey soils:A meta-analysis[J]. Earth-Science Reviews,2021,221:103802. DOI: 10.1016/j.earscirev.2021.103802
[15] 成小彬. 水锤振荡器的关键技术研究与结构设计[D]. 北京:中国石油大学(北京),2017. [CHENG Xiaobin. Research on key technology of water hammer oscillator and its structural design[D]. Beijing:china university of petroleum (bei jing),2017. (in Chinese with English abstract) CHENG Xiaobin. Research on key technology of water hammer oscillator and its structural design[D]. Beijing: china university of petroleum (bei jing), 2017. (in Chinese with English abstract)
[16] IMESON A C,VIS M. Assessing soil aggregate stability by water-drop impact and ultrasonic dispersion[J]. Geoderma,1984,34(3/4):185 − 200.
[17] MCCALLA T M. Water-drop method of determining stability of soil structure[J]. Soil Science,1944,58(2):117 − 122. DOI: 10.1097/00010694-194408000-00002
[18] ALMAJMAIE A,HARDIE M,ACUNA T,et al. Evaluation of methods for determining soil aggregate stability[J]. Soil and Tillage Research,2017,167:39 − 45. DOI: 10.1016/j.still.2016.11.003
[19] RAWLINS B G,TURNER G,WRAGG J,et al. An improved method for measurement of soil aggregate stability using laser granulometry applied at regional scale[J]. European Journal of Soil Science,2015,66(3):604 − 614. DOI: 10.1111/ejss.12250
[20] 张绍科,胡卸文,王严,等. 四川省冕宁县华岩子沟火后泥石流成灾机理[J]. 中国地质灾害与防治学报,2021,32(5):79 − 85. [ZHANG Shaoke,HU Xiewen,WANG Yan,et al. Disaster mechanism of post-fire debris flow in Huayanzi gully,Mianning County,Sichuan Province[J]. The Chinese Journal of Geological Hazard and Control,2021,32(5):79 − 85. (in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.2021.05-09 ZHANG Shaoke, HU Xiewen, WANG Yan, et al . Disaster mechanism of post-fire debris flow in Huayanzi gully, Mianning County, Sichuan Province[J]. The Chinese Journal of Geological Hazard and Control,2021 ,32 (5 ):79 −85 . (in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.2021.05-09[21] 胡彩莉,马玉贞,郭超,等. 烧失量法测定土壤有机质含量的实验条件探究[J]. 地球与环境,2016,44(1):110 − 118. [HU Caili,MA Yuzhen,GUO Chao,et al. Optimization of the experiment conditions for estimating organic matter content with loss-on-ignition method[J]. Earth and Environment,2016,44(1):110 − 118. (in Chinese with English abstract) DOI: 10.14050/j.cnki.1672-9250.2016.01.015 HU Caili, MA Yuzhen, GUO Chao, et al . Optimization of the experiment conditions for estimating organic matter content with loss-on-ignition method[J]. Earth and Environment,2016 ,44 (1 ):110 −118 . (in Chinese with English abstract) DOI: 10.14050/j.cnki.1672-9250.2016.01.015[22] 孔凡伟. 如何精测土壤容重[J]. 黑龙江农业科学,2010,(10):178. [KONG Fanwei. How to accurately measure soil bulk density[J]. Heilongjiang Agricultural Sciences,2010,(10):178. (in Chinese) KONG Fanwei . How to accurately measure soil bulk density[J]. Heilongjiang Agricultural Sciences,2010 , (10 ):178 . (in Chinese)[23] BISDOM E B A,DEKKER L W,TH SCHOUTE J F. Water repellency of sieve fractions from sandy soils and relationships with organic material and soil structure[J]. Geoderma,1993,56(1/2/3/4):105 − 118.
[24] 程佳会,范成伟,徐晓兵,等. 正交试验下黏土导热系数影响分析[J]. 江西建材,2022(3):34 − 36. [CHENG Jiahui,FAN Chengwei,XU Xiaobing,et al. Influence analysis of thermal conductivity of clay under orthogonal test[J]. Jiangxi Building Materials,2022(3):34 − 36. (in Chinese with English abstract) CHENG Jiahui, FAN Chengwei, XU Xiaobing, et al . Influence analysis of thermal conductivity of clay under orthogonal test[J]. Jiangxi Building Materials,2022 (3 ):34 −36 . (in Chinese with English abstract)[25] 侯勋,方刚. 基于肯德尔系数的改进ID3算法[J]. 科学技术创新,2021,(22):21 − 22. [HOU Xun,FANG Gang. Improved ID3 algorithm based on Kendall coefficient[J]. Scientific and Technological Innovation,2021,(22):21 − 22. (in Chinese) HOU Xun, FANG Gang . Improved ID3 algorithm based on Kendall coefficient[J]. Scientific and Technological Innovation,2021 , (22 ):21 −22 . (in Chinese)[26] LI Qiwen,AHN S,KIM T,et al. Post-fire impacts of vegetation burning on soil properties and water repellency in a pine forest,South Korea a[J]. Forests,2021,12(6):708,2021
[27] 李萌. 热声制冷机板叠内工质流动与换热规律研究[D]. 包头:内蒙古科技大学,2021. [LI Meng. Study on the flow and heat transfer laws of working medium in the stack of thermoacoustic refrigerator[D]. Baotou:Inner Mongolia University of Science & Technology,2021. (in Chinese with English abstract) LI Meng. Study on the flow and heat transfer laws of working medium in the stack of thermoacoustic refrigerator[D]. Baotou: Inner Mongolia University of Science & Technology, 2021. (in Chinese with English abstract)