Application of multiple InSAR techniques and SAR data from multi-sources to landslide deformation monitoring: A case study of the Zhixincun landslide in Jilin Province
-
摘要:
为实现对吉林省治新村滑坡的有效监测，文章选取2017年27景Sentinel-1A数据，基于小基线雷达干涉测量技术（SBAS-InSAR）对治新村滑坡进行形变监测，分析了其时序演化态势。选用2016、2017年2景ALOS-2数据，采用差分雷达干涉测量技术（D-InSAR）监测了该滑坡形变体的特征。SBAS-InSAR对滑坡形变时序演化态势进行监测，而D-InSAR则主要对滑坡具体的形变体进行形变监测，且L波段的ALOS-2数据穿透性强于C波段的Sentinel-1A数据，可以获得更完整的干涉信息，两者监测结果可交叉验证，提高结果的可靠性。SBAS-InSAR监测结果表明：治新村滑坡汇水区滑坡后缘在监测期间发生了沉降，并且在2017年7月5—29日期间滑坡后缘地表沉降达12.47 mm，监测期间平均沉降速率为2.88 mm/a；位于山谷的受威胁居民区发生了抬升，至2017年12月8日平均累计抬升达19.59 mm，监测期间平均抬升速率19.99 mm/a。D-InSAR结果显示：治新村滑坡汇水区斜坡存在5处主要变形体，面积最大变形体17973 m2，位于西侧斜坡，最不稳定变形体位于斜坡东侧，监测期间平均累计形变量最大达49.9 mm。两种监测方法都表明，滑坡灾害威胁主要来自植被覆盖较差的西侧斜坡，雨季是治新村滑坡灾害防治的重点时期。
-
关键词:
- 滑坡 /
- 治新村 /
- SBAS-InSAR /
- D-InSAR /
- ALOS-2 /
- Sentinel-1A
Abstract:In order to realize effective monitoring of Zhixincun landslide, this paper selected 27 sentinel-1A data in 2017, and conducted deformation monitoring of Zhixincun landslide based on small baseline radar interferometry technology (SBAS-InSAR), and analyzed its temporal evolution situation. Using ALOS-2 data from 2016 and 2017, differential radar interferometry (D-InSAR) was used to monitor the characteristics of the landslide variant. SBAS-InSAR monitors the temporal evolution situation of landslide deformation, while D-InSAR mainly monitors the deformation of specific landslide shape and variation. Moreover, the penetration of L-band ALOS-2 data is stronger than that of C-band sentinel-1A data, which can obtain more complete interference information. The monitoring results of both can be cross-verified. Improve the reliability of the results. The SBAS-InSAR monitoring results showed that the slope end of the landslide catchment area in Zhixincun had subsidence during the monitoring period, and the surface subsidence at the landslide end reached 12.47mm from July 5 to July 29, with an average subsidence rate of 2.88mm/a during the monitoring period. Uplift occurred in the threatened residential areas in the valley, with an average cumulative uplift of 19.59mm on December 8 and an average uplift rate of 19.99mm/a during the monitoring period. The D-InSAR results showed that there were five major deformations on the slope of Zhixincun landslide catchment area. The largest deformations with an area of 17 973m2 were located on the west side of the slope, and the most unstable deformations were located on the east side of the slope. The average cumulative shape variable reached 49.9mm during the monitoring period. Both monitoring methods showed that the threat of landslide disaster mainly came from the west slope with poor vegetation cover, and the rainy season was the key period of landslide disaster prevention and control in Zhixincun.
-
Keywords:
- landslide /
- Zhixincun /
- SBAS-InSAR /
- D-InSAR /
- ALOS-2 /
- Sentinel-1A
-
0. 引言
青藏高原东北部位于中国东部季风半湿润区、西北内陆干旱区和青藏高原高寒区的交汇部，是气候与环境变化复杂而敏感的地区[1]。区内黄土广泛发育，且具有沉积厚度大、地层较完整等特征，特别是西宁盆地及其周边地区，是整个青藏高原东北部主要的黄土沉积区[2]。受新构造运动的影响，区内黄土侵蚀切割强烈，沟壑纵横，呈现梁、峁、塬和黄土丘陵等地貌景观[3]。由于西宁盆地及其周边地区降雨多为短时暴雨和阵雨，加之山麓斜坡坡度较大，植被覆盖度低，故降雨入渗深度相对较浅，导致盆地及其周边地区黄土浅层滑坡广泛发育，且具有分布规律性差、前期变形小、分布范围大、面小点多等特征，很难进行有效预测和提前防护，给区内黄土地区生态地质环境和工程安全带来严重威胁[4 − 6]。
植物通过力学效应和水文效应两方面作用，能够控制水土流失、浅层滑坡等地质灾害现象的发生[7 − 8]。植物固土护坡力学效应在于根系较土体而言具有显著较高的抗拉强度，当含有植物根系的土体受到剪切或张拉作用时，根系的存在能够将土体内部剪切滑动面上的剪应力和张裂面上的张应力转化为自身所承受的拉应力，并通过根-土界面剪应力将拉应力传递至稳定土体，从而起到增强土体抗剪强度和提高边坡稳定性的效果[9 − 12]。根系固土护坡力学效应可分为加筋和锚固两方面作用[13]。由于黄土浅层滑坡深度多为2～3 m[4, 6]，而草本植物和幼龄木本植物的根系主要分布在地表以下1～1.5 m深度范围内[14 − 16]，故对于大多数草本植物和幼龄木本植物而言，由于根系分布深度有限，因此主要起到加筋土体的作用，体现在增强根系分布层土体抗剪强度。而对于一些龄期较长的木本植物而言，根系能生长至深部土体，并足以穿过浅层滑坡滑动面，具有锚固浅层土体的作用[17 − 18]。因此，一般认为，大龄期木本植物根系锚固作用是实现区内黄土浅层滑坡生态防护的关键。
灌木林是青藏高原东北部黄土区主要的人工林和天然林类型，区内常见的灌木植物如柠条锦鸡儿（Caragana korshinskii Kom.）、柽柳（Tamarix chinensis Lour.）和沙棘（Hippophae rhamnoides Linn.）等[19 − 20]。这些灌木植物中，柠条锦鸡儿与同属的小叶锦鸡儿（Caragana microphylla Lam.）和甘蒙锦鸡儿（Caragana opulens Kom.）具有适应性强、抗寒耐旱、抗风蚀耐啃食和易繁殖等特性，是区内防风固沙、水土保持、荒山绿化以及高速公路护坡绿化的首选树种[21 − 23]。柠条锦鸡儿根系具有极强的向深层湿润土体延伸的特性[21 − 22]，是典型的垂直根型植物，其主根明显，侧根发达，属于主侧根均衡发育型灌木[22]。胡崇礼[24]对生长于西宁盆地西山的柠条锦鸡儿根系分布特征进行了调查研究，结果表明生长期为18 a的柠条锦鸡儿根系深度可达7.0 m。因此，从根系分布特征的角度而言，柠条锦鸡儿等垂直根型灌木具有锚固斜（边）坡浅层土体的潜力。然而，目前有关植物锚固机理方面的研究，多是以植物原位拉拔试验、根系模型拉拔试验或根系拉拔过程数值模拟为基础，着重讨论根系抗拔出特性及其影响因素[25 − 31]。且已有的根系抗拔出试验中，根系是由地表垂直[25 − 31]或水平荷载作用下被拔出土体[27, 32]。这与浅层滑坡过程中，受滑坡下滑力作用下根系的实际受力状态存在区别。因为分布于滑动面上部和下部的根系受力方向实际上是不同的，且两部分根系在锚固土体过程中所发挥的作用亦不尽相同[33]。因此，拉拔试验所获得的根系抗拔出力更适合用于评价植物在风力和地表径流作用下的立地稳定性，而无法直接用于定量评价植物根系对斜（边）坡浅层土体稳定性的贡献。故在探讨植物根系锚固效应时，有必要先着手分析滑动面上下两部分根系的根-土相互作用过程，在此基础上，合理确定根系锚固力大小和计算模型，实现定量评价植物锚固作用对斜（边）坡浅层土体稳定性的贡献。
基于此，该项研究以广泛分布于青藏高原东北部黄土区的灌木柠条锦鸡儿为研究对象。首先阐明黄土浅层滑坡滑动时，分布于滑动面上下两部分的柠条锦鸡儿根系根-土相互作用力学过程，明确根系锚固机理，并进一步提出根系对斜（边）浅层土体稳定性贡献定量评价方法。以此为基础，以生长期为11 a的柠条锦鸡儿为试验对象，通过根系挖掘试验、根系原位拉拔试验和理论计算相结合的方法，合理确定假定滑动面条件下根系锚固力取值，并进一步定量评价了降雨入渗条件下，柠条锦鸡儿根系对黄土斜（边）浅层土体稳定性的增强作用。该项研究成果对于深入探讨灌木植物护坡力学机制，进一步实现合理利用植物有效地防治黄土地区水土流失、浅层滑坡等斜坡地质灾害具有理论指导意义和实际应用价值。
1. 柠条锦鸡儿根系锚固理论
1.1 柠条锦鸡儿根系锚固机理分析
植物根系锚固效应作用的对象是斜（边）坡浅层土体[13]。对于黄土斜（边）坡而言，浅层滑坡主要是降雨入渗诱发的[6, 18]。已有研究表明，黄土高原地区每年降水最大入渗深度为1～3 m左右，半湿润地区大丰水年份方可达到5 m以下[34]。而浅层滑坡深度多为2～3 m[4, 6]。在降雨入渗深度范围内，土体含水量增大导致浅层土体自重应力增大。与此同时，含水量增大能够显著降低黄土抗剪强度，从而诱发斜（边）坡浅层土体失稳滑动。
当植物根系分布深度大于浅层滑坡滑动面深度时，根系能够发挥其锚固作用[17 − 18]。图1所示为柠条锦鸡儿根系锚固黄土斜坡浅层土体的实例，以及黄土浅层滑坡发生时，柠条锦鸡儿根系根-土相互作用过程示意图。由图1（b）可知，当滑动面贯通滑体下滑或潜在滑动面局部发生蠕滑变形时，滑体和滑床将沿着滑动面发生相对错动，穿过滑动面的柠条锦鸡儿根系将承受拉应力，并通过根-土相互作用，将拉应力传递至根周土体（图1（b）假定根-土界面剪应力与拉拔力满足静力平衡条件，且已达到极限平衡状态）。位于滑动面上下两侧的根系所受拉应力方向均朝向滑动面，并在根-土界面产生背离滑动面的抗拔出力。虽然穿过滑动面的根系作用机理类似于全长粘结式锚杆，然而，由于根系在地表没有“锚头”结构，因此，滑动面上下两侧根系与土体之间的抗拔出力，实际为一对作用力与反作用力。为了便于区别，该项研究将滑动面以上根系受拉产生的背离滑动面的抗拔出力称为锚固反力，用Fj表示；将滑动面以下根系抗拔出力用Fm表示。根系锚固反力对于根系锚固斜（边）坡浅层土体具有重要意义。因为滑动面以上根系一旦与滑体分离，则将不再提供锚固反力，无法继续将滑坡下滑力传递至滑动面以下根系，从而导致锚固作用失效。也就是说，如果视滑动面以下的根系为锚固段，则锚固段根系所能提供锚固力的大小，在一定程度上取决于滑动面以上根系锚固反力的大小，反之亦然。综上所述，在滑动面一定的情况下，根系所能提供的实际锚固力大小为锚固段根系最大抗拔出力和滑动面以上根系锚固反力之间的最小值。由于黄土地区浅层滑坡主要受降雨入渗影响，因此滑动面以上属于降雨入渗区。入渗区内土体和根-土界面抗剪强度降低，根-土相互作用产生的抗拔出力相对于降雨入渗之前将显著下降。因此，降雨条件下滑动面以上根系锚固反力的大小，是评价根系锚固效果时需着重分析和考虑的条件。
![]()
图 1 柠条锦鸡儿根系锚固黄土斜坡浅层土体及其锚固机理示意图Figure 1. Schematic diagram of the shallow soil body of the C. korshinskii roots anchoring loess and its anchoring mechanism1.2 柠条锦鸡儿根系对黄土斜（边）坡浅层土体稳定性贡献评价计算模型
该项研究在极限平衡理论基础上，参考《土钉支护技术规范》（GJB 5055—2006）[35]中边坡土钉支护普通条分法稳定性分析计算公式，提出考虑柠条锦鸡儿根系锚固作用条件下，斜（边）坡浅层土体稳定性计算式（1）。
(1) 式中：Fs——斜（边）稳定性系数；
Wi——土条i自重/kN；
Qi——作用于土条i的地面、地下载荷/kN；
αi——土条i圆弧破坏面切线与水平面的夹角/（°）；
φj——土条i圆弧破坏面处第j层土的内摩擦角/（°）；
Rk——穿过土条破坏面上第k排灌木根系所能提供的锚固力/kN，取锚固段根系抗拔出力Fm和滑动面以上根系锚固反力Fj的最小值；
Shk——第k排灌木根系的水平间距/m；
βk——第k排灌木根系锚固段轴线与该处破坏面切线之间的夹角/（°）；
cj——土条i圆弧破坏面处第j层土的黏聚力/kPa；
滑动面以上根系锚固反力可由式（2）计算[36]。
(2) 式中：Fj——锚固反力/kN；
τrs——根-土界面抗剪强度/kPa；
D——单根平均根径/m；
crs——根-土间黏聚力/kPa；
φrs——根-土间摩擦角/（°）；
K0——静止土压力系数，取0.6；
γ——土体重度/（kN·m−3）；
l——根系垂直埋置深度/m。
由此可知，在明确滑动面几何形态和滑动面以上柠条锦鸡儿根系特征的基础上，即可进行锚固反力计算。计算锚固反力时，结合柠条锦鸡儿侧根生长具有不确定性的实际情况，出于安全储备的角度，该项研究提出仅考虑由柠条锦鸡儿主根所贡献的锚固反力，而不考虑滑动面以上侧根对锚固反力的贡献。同时，以根系埋深代替主根实际长度进行锚固反力计算。因此，仅需明确柠条锦鸡儿根系主根在滑动面以上的平均根径和根-土界面抗剪强度指标，即可计算根系锚固反力。基于此，通过原位拉拔试验进一步确定滑动面以下根系抗拔出力数值，即可合理确定作用在滑动面上的根系锚固力取值，并通过式（1）计算出根系锚固作用条件下浅层滑坡稳定性系数。
2. 实例分析
2.1 假定滑动面以上柠条锦鸡儿根系分布特征与根周土体物理力学特性
试验区位于在青海省西宁市大有山，为天然黄土斜坡。试验区斜坡坡向向东，坡度22°～43°。试验区内生长的灌木柠条锦鸡儿属于人工灌木林。试验选取1株生长期为11 a（年龄通过年轮确定）的柠条锦鸡儿进行，采用原位挖掘法调查假定滑动面以上柠条锦鸡儿根系分布特征。考虑到黄土区浅层滑坡深度多为2～3 m。鉴于此，该项研究最初设定假定滑动面最大深度为地表以下2.5 m，该深度亦为此次原位根系挖掘调查的下限。但是，在实际挖掘过程中，受控于探坑施工空间和坑周排土空间，并考虑到坑周堆土稳定性，最终挖掘深度为2 m（假定滑动面最大深度亦改调整为2 m）。现场挖掘过程如图2所示。为了保证获得原位根系的准确信息，边挖掘边测量。现场测量的根系形态学指标包括根幅、主根长、主根根径和侧根根数、根长等。
![]()
图 2 假定滑动面以上根系原位挖掘试验Figure 2. In-situ excavation test of root system above the hypothetical sliding surface表1为试验株地上植株和0～2 m深度范围内根系形态学指标统计数据。由该表可知，在0～2 m深度范围内，柠条锦鸡儿总根数达到27根，其中，除主根外还包括一级侧根20根，二级侧根和三级侧根各3根，主根平均根径为0.027 m。累计总根长34 m，根表面积0.87 m2，根幅为1.99 m。
表 1 试验株地上植株和根系形态学指标统计Table 1. Statistical analysis of morphological indexes of above-ground plants and root system above the hypothetical sliding surface of the testing plants植株及根系指标 树龄/a 冠幅/m 株高/m 总根数/根 根幅/m 一级侧根/根 二级侧根/根 三级侧根/根 主根平均根径/ m 总根长/m 根表面积/m2 取值 11 1.76 2.20 27 1.99 20 3 3 0.027 34 0.87 在根系调查的同时，按0.3 m间隔用环刀采集根周原状土体试样，每层采集4个试样，同时用铝土盒采集3个天然土试样。带回实验室后分别用于天然密度、抗剪强度指标和天然含水量测试。试验点植株根周土体物理力学特性如表2所示。深度范围为0～2 m
表 2 试验株根周土体物理力学特性Table 2. Physical and mechanical properties of root-soil interface of the testing plants物理力学性质 天然密度/（g·cm−2） 天然含水量/% 黏聚力/kPa 内摩擦角/（°） 取值 1.40±0.08 7.56±1.21 15.83±7.97 13.56±4.25 2.2 假定滑动面以下锚固段柠条锦鸡儿根系原位拉拔试验
通过现场挖掘表明，试验株分布于地表2 m以下深度的主根和侧根均属于垂直根系（垂直根系指根轴线与水平线夹角为60°～90°的根系[37]。因此，该项研究通过原位拉拔试验获得其抗拔出力。原位拉拔试验仪器采用自行设计加工的便携式野外拉拔仪，拉拔仪由支架、导链、拉力传感器（最大量程30 kN，精度0.5%）、拉线式位移传感器（最大量程1 m，精度0.5%）、数据采集系统、笔记本电脑组成，由蓄电池提供电源。拉拔试验时，首先将根系固定在导链下端挂钩上，然后人工匀速上拉导链（平均拉拔速率为21.50±3.71 cm/min），直到根系拉拔力超过峰值并达到位移量程极限时结束试验。试验过程中，数据采集系统将连续采集根系抗拔出力和对应的相对位移并在笔记本电脑终端软件页面显示，试验结束后，软件会自动保存此次试验相关数据。由于除主根外，试验株在地表2 m深度仍然有4根根系垂直向下生长。因此，共开展5次原位拉拔试验。
原位拉拔试验数据如表3所示。由表3可知，在地表以下2 m位置处，试验株垂直根系根径为0.005～0.020 m，抗拔出力为0.206～1.666 kN。拉拔试验过程中，5根根系并未被完整拔出土体，而是以根皮和木质部同时拉断的形式断裂。
表 3 试验株锚固段根系原位拉拔试验结果Table 3. In-situ pullout test results of root system in anchoring section of testing plants序号 抗拔出力/kN 拉拔端根径/m 断裂方式 1 0.206 0.007 根皮和木质部同时拉断 2 1.243 0.020 3 0.519 0.005 4 1.666 0.012 5 1.168 0.011 2.3 柠条锦鸡儿根系对黄土边坡浅层土体稳定性贡献评价
2.3.1 计算模型建立
在GeoStudio软件Slope/W模块中建立边坡模型（图3）。设置边坡坡度为30°，坡面长度为15 m，坡高为7.5 m。由于该项研究中假定滑动面深度、野外根系挖掘调查深度以及根系原位拉拔试验中根系出露深度均为地表以下2 m。因此，为获得最大深度为2 m的潜在滑动面，以便于进行柠条锦鸡儿根系锚固作用对黄土边坡浅层土体稳定性贡献评价，模型中设置地表以下2 m深度范围为降雨入渗区，且假设降雨入渗深度范围内土体处于饱和状态。由表2可知0～2 m深度范围内，试验株根周土体天然抗剪强度平均值为15.83±7.97 kPa和13.56°±4.25°。韩帅[38]的研究指出，原状黄土在饱和状态下的黏聚力与内摩擦角相对于其在天然含水量状态下的黏聚力和内摩擦角分别降低了47.5%和7.7%，基于此，该项研究中，饱和土体抗剪强度指标分别取8.31 kPa和12.51°。饱和重度取17.7 kN/m3。模型中入渗层下部非饱和土体物理力学参数采用表2中所示的0～2 m深度范围内土体物理力学参数。采用Morgenstern-Price法获得边坡潜在滑动面位置和滑体各条块特征参数。由图3所示，边坡潜在滑动面为圆弧形，滑动面最大深度与降雨入渗深度一致。
2.3.2 根系锚固反力计算与锚固力取值
由于黄土处于饱和状态时，土体黏聚力和内摩擦角降低幅度较大，加之根和根周土体弹性模量存在较大差异，故此时根-土界面破坏基本为根周土体剪切破坏[39]。因此，在降雨工况下，根系锚固反力计算时，根-土界面黏聚力和内摩擦角亦选择饱和黄土抗剪强度指标。该项研究结合边坡模型潜在滑动面条块特征，分别计算了试验株生长在不同条块中心轴线位置处的锚固反力数值，具体计算参数和结果如表4所示。
表 4 试验株生长于不同条块时锚固反力计算参数与结果Table 4. Calculation parameters and results of anchoring reaction force of testing plants growing on different blocks条块号 根-土间黏聚力/kPa 根-土间摩擦角/（°） 单根平均根径/m 根系垂直埋置深度/m 静止土压力系数 土体重度/（kN·m−3） 锚固反力/kN 1 8.31 12.51 0.040 0.461 0.6 17.7 0.513 2 8.31 12.51 0.033 1.204 0.6 17.7 1.207 3 8.31 12.51 0.029 1.655 0.6 17.7 1.546 4 8.31 12.51 0.027 1.912 0.6 17.7 1.709 5 8.31 12.51 0.027 2.000 0.6 17.7 1.805 6 8.31 12.51 0.027 2.000 0.6 17.7 1.805 7 8.31 12.51 0.027 1.930 0.6 17.7 1.729 8 8.31 12.51 0.029 1.748 0.6 17.7 1.650 9 8.31 12.51 0.033 1.485 0.6 17.7 1.529 10 8.31 12.51 0.033 1.147 0.6 17.7 1.141 11 8.31 12.51 0.036 0.737 0.6 17.7 0.762 12 8.31 12.51 0.040 0.257 0.6 17.7 0.278 注：表中单根平均根径取对应条块滑动面深度范围内主根平均根径；根系垂直埋置深度取值为对应条块滑动面深度，即滑动面中间点至地表的距离。 如前文所述，根系实际作用于滑动面上的锚固力应取滑动面以下锚固段根系最大抗拔出力和滑动面以上根系锚固反力的最小值。由原位拉拔试验实际测量数据显示，在地表以下2 m位置处，试验株5根垂直单根之间最大水平间距为0.65 m，故可同时作用于边坡模型中的同一块条块滑动面上。因此，试验株根系在深度为2 m的滑动面上所能提供的抗拔出力合力为4.802 kN。由表4可知，边坡模型潜在滑动面第5和6号滑块对应的滑动面深度最大，均为2 m。因此，试验株根系作用于这两块条块上的锚固段根系抗拔出力为4.802 kN。与此同时，第5和6号滑块对应的根系锚固反力均为1.805 kN，小于滑动面下部锚固段根系抗拔出力。相比之下，其他各条块滑动面深度均有所降低，而对应的根系锚固反力也随之降低。与此同时，由于滑动面以下根系埋深增加，故其他各条块滑动面以下根系抗拔出力理论上会相应增大。因此，在该项研究中，单株柠条锦鸡儿根系作用于各条块时，滑动面以上根系锚固反力均低于锚固段根系抗拔出力。因此，在计算根系锚固作用下边坡潜在滑动面稳定性系数时，锚固力取各条块对应的根系锚固反力值。
2.3.3 柠条锦鸡儿根系锚固作用条件下边坡潜在滑动面稳定性系数
采用式（1）计算无根系锚固作用和有根系锚固作用条件下模型边坡潜在滑动面稳定性系数。其中，有根系锚固作用条件下潜在滑动面稳定性系数计算时，首先计算单株根系锚固力分别作用于1～12号条块时，边坡潜在滑动面稳定性系数。在此基础上，进一步计算多株根系同时作用条件下，边坡潜在滑动面稳定性系数。经过野外测量，柠条锦鸡儿沿斜坡倾向方向行间距约为3～4 m，因此，在计算多株根系锚固作用条件下边坡潜在滑动面稳定性系数时，设计三组植株分布情况，每组中植株间距2块条块（约3 m）。即A组：植物根系作用于1、4、7和10号条块、B组：作用于2、5、8和11号条块以及C组：作用于3、6、9和12号条块。
各条件下边坡潜在滑动面稳定性系数计算结果如表5所示。由表5可知，通过式（1）计算得到的无根系锚固作用条件下边坡潜在滑动面稳定性系数为1.3370。当单株根系锚固于潜在滑动面不同条块上时，潜在滑动面稳定性系数均有提高，但增幅有限，为0.018%～0.427%。其中，当根系锚固力作用于潜在滑动面第4块条块时，潜在滑动面稳定性系数最大，为1.3427。当根系锚固力作用于潜在滑动面中上部位置处（2～7号）的条块时，潜在滑动面稳定性系数相对高于根系锚固力作用于最顶部（1号）和下部位置处（8～12号）条块的稳定性系数。这主要与条块下滑力大小、对应的锚固反力大小和βk值大小存在关系。由表5亦可知，当4株柠条锦鸡儿根系等间距锚固于潜在滑动面时，不同组合条件下边坡潜在滑动面稳定性系数为1.3508～1.3519，较无根系锚固作用条件下增长1.035%～1.111%，且显著高于（P<0.05，ANOVA）单株根系作用时对应的稳定性系数。总体而言，试验株根系锚固作用能够提高降雨入渗条件下黄土斜（边）坡浅层土体稳定性，但是作用效果有限。
表 5 柠条锦鸡儿根系锚固作用下边坡潜在滑动面稳定性系数计算结果Table 5. Calculation results of stability coefficient of potential sliding surface in the C. korshinskii roots system anchored slope锚固力作用
条块号土条
自重
/kPa土条
宽度
/m圆弧破坏面切线与
水平面的夹角/（°）根系锚固段轴线与
该处破坏面切线
之间的夹角/（°）根系水平
间距/m圆弧破坏面
处土体内
摩擦角/（°）圆弧破坏面处
土体黏聚力
/kPa根系锚
固力/kN稳定性
系数增幅/% 1 7.374 0.903 57.966 32.034 2 12.510 8.310 0.513 1.3391 0.156 2 19.252 0.903 50.216 39.784 2 12.510 8.310 1.207 1.3417 0.349 3 26.465 0.903 43.601 46.399 2 12.510 8.310 1.546 1.3426 0.417 4 30.577 0.903 37.659 52.341 2 12.510 8.310 1.709 1.3427 0.427 5 29.967 0.839 32.352 57.648 2 12.510 8.310 1.805 1.3425 0.414 6 29.967 0.839 27.531 62.469 2 12.510 8.310 1.805 1.3420 0.378 7 30.434 0.891 22.776 67.224 2 12.510 8.310 1.729 1.3413 0.324 8 27.560 0.891 18.043 71.957 2 12.510 8.310 1.650 1.3406 0.272 9 23.414 0.891 13.435 76.565 2 12.510 8.310 1.529 1.3399 0.217 10 18.082 0.891 8.914 81.086 2 12.510 8.310 0.513 1.3388 0.134 11 11.619 0.891 4.449 85.551 2 12.510 8.310 1.207 1.3379 0.071 12 4.055 0.891 0.010 89.990 2 12.510 8.310 1.546 1.3372 0.018 A - - - - - - - - 1.3510 1.047 B - - - - - - - - 1.3519 1.111 C - - - - - - - - 1.3508 1.035 对照 - - - - - - - - 1.3370 - 注：该项研究根系所能提供的锚固力取值为试验株根系单株作用在各条块时对应的锚固反力，并假设锚固力方向铅垂向下，故灌木根系锚固段轴线与该处破坏面切线之间的夹角为滑块滑动面切线与铅垂方向的夹角；灌木根系的水平间距取值由野外实际测量数据确定。 值得一提的是，该项研究在计算根系锚固反力时，出于安全储备仅考虑了主根而未考虑侧根对锚固反力的贡献。且以根系埋深代替单根实际长度进行锚固反力计算。虽然计算结果在理论上低于试验株根系对模型边坡潜在滑动面稳定性的真实贡献，相对较为保守，但是其结果在一定程度上依然能够反映该年龄段柠条锦鸡儿根系锚固作用对降雨条件下黄土斜（边）坡浅层土体稳定性的作用效果。与此同时，该设定亦可以简化分析和计算过程，便于定量评价柠条锦鸡儿及其他垂直根型植物对降雨型黄土浅层滑坡稳定性的贡献。此外，该项研究仅以生长期为11 a，处于中龄期[40]的灌木柠条锦鸡儿作为研究对象。由于植物根系生长特征是随着其生长期的增长而不断变化的，随着生长期增大，植物垂直根系和水平根系数量、长度、根径以及根-土界面抗剪特性还会进一步变化。故随着龄期变化，根系实际作用于斜（边）坡浅层土体的锚固力和锚固效果势必存在动态变化的特性。因此，掌握植物在特定气候土壤条件区域内根系随龄期变化的生长特征，是明确其根系固土护坡效应动态规律，科学合理利用植物有效地防治斜（边）浅表层土体地质灾害的基础。
3. 结论
（1）由于根系在地表没有“锚头”结构，故在确定滑动面几何特征的情况下，根系所能提供的实际锚固力大小取滑动面以下锚固段根系最大抗拔出力和滑动面以上根系锚固反力之间的最小值较为合理。
（2）生长期为11 a的单株柠条锦鸡儿根系锚固于最大厚度为2 m的圆弧形滑动面不同条块上时，潜在滑动面稳定性系数增幅为0.018%～0.427%；当4株柠条锦鸡儿根系以2块条块的间距（约3 m）作用于潜在滑动面时，潜在滑动面稳定性系数可提高1.035%～1.111%，显著高于（P<0.05，ANOVA）单株根系作用时的稳定性系数。试验株根系锚固作用能够提高降雨入渗条件下黄土斜（边）坡浅层土体稳定性，但是作用效果有限。
该项研究成果对于研究区采用垂直根型灌木植物进行浅层滑坡等地质灾害现象的生态防护，以及对进一步探讨根-土相互作用机制具有实际指导意义和理论研究价值。
-
![]()
图 3 SBAS-InSAR时空基线图与像对连接图
Figure 3. SBAS-InSAR space-time baseline diagram and image pair connection diagram
![]()
图 5 SBAS-InSAR监测点时序特征
Figure 5. Temporal characteristics of SBAS-InSAR monitoring points
表 1 Sentinel-1A数据集
Table 1 Sentinel-1A data set
序号 成像日期 时间基线/d 空间基线/m 01 2017-01-06 60 43.8 02 2017-01-18 48 19.9 03 2017-01-30 36 75.9 04 2017-02-11 24 170.1 05 2017-02-23 12 109.8 06 2017-03-07 0 0 07 2017-03-19 12 −15.5 08 2017-03-31 24 32.9 09 2017-04-12 36 57.4 10 2017-04-24 48 128.3 11 2017-05-06 60 92.7 12 2017-05-18 72 47.9 13 2017-05-30 84 −70.3 14 2017-06-11 96 104.5 15 2017-06-23 108 95.2 16 2017-07-05 120 25.1 17 2017-07-29 144 43.2 18 2017-08-10 156 71.3 19 2017-08-22 168 27.8 20 2017-09-03 180 62.6 21 2017-09-15 192 107.1 22 2017-09-27 204 62.3 23 2017-10-09 216 70.1 24 2017-10-21 228 122.2 25 2017-11-02 240 145.3 26 2017-11-26 264 56.1 27 2017-12-08 276 61.5 
下载: 导出CSV
表 2 ALOS-2影像信息
Table 2 ALOS-2 image information
序号 成像日期 时间基线/d 空间基线/m 01 2016-07-26 392 −225.5 02 2017-08-22
下载: 导出CSV
表 3 斜坡变形体
Table 3 Deformation of slope
编号 面积/m² 最大沉降量/mm 最小沉降量/mm 平均沉降量/mm 1 5318 58.6 29.4 43.2 2 17973 50.6 24.2 36.4 3 4636 68.7 29.8 43.5 4 3043 71.8 31.3 49.9 5 11281 47.6 20.1 38.3
下载: 导出CSV
-
[1] 殷跃平,王文沛. 高位远程滑坡动力侵蚀犁切计算模型研究[J]. 岩石力学与工程学报,2020,39(8):1513 − 1521. [YIN Yueping,WANG Wenpei. Study on calculation model of dynamic erosion plowing of high-level remote landslide[J]. Chinese Journal of Rock Mechanics and Engineering,2020,39(8):1513 − 1521. (in Chinese with English abstract) Yin Yueping, Wang Wenpei. Study on calculation model of dynamic erosion plowing of high-level remote landslide[J]. Chinese Journal of Rock Mechanics and Engineering, 2020, 39（8）: 1513-1521. （in Chinese with English abstract）
[2] 黄海峰,林海玉,吕奕铭,等. 基于小型无人机遥感的单体地质灾害应急调查方法与实践[J]. 工程地质学报,2017,25(2):447 − 454. [HUANG Haifeng,LIN Haiyu,LYU Yiming,et al. Micro unmanned aerial vehicle based remote sensing method and application for emergency survey of individual geohazard[J]. Journal of Engineering Geology,2017,25(2):447 − 454. (in Chinese with English abstract) HUANG Haifeng, LIN Haiyu, LYU Yiming, et al. Micro unmanned aerial vehicle based remote sensing method and application for emergency survey of individual geohazard[J]. Journal of Engineering Geology, 2017, 25（2）: 447-454. （in Chinese with English abstract）
[3] 朱庆,曾浩炜,丁雨淋,等. 重大滑坡隐患分析方法综述[J]. 测绘学报,2019,48(12):1551 − 1561. [ZHU Qing,ZENG Haowei,DING Yulin,et al. A review of major potential landslide hazards analysis[J]. Acta Geodaetica et Cartographica Sinica,2019,48(12):1551 − 1561. (in Chinese with English abstract) ZHU Qing, ZENG Haowei, DING Yulin, et al. A review of major potential landslide hazards analysis[J]. Acta Geodaetica et Cartographica Sinica, 2019, 48（12）: 1551-1561. （in Chinese with English abstract）
[4] 张勤,黄观文,杨成生. 地质灾害监测预警中的精密空间对地观测技术[J]. 测绘学报,2017,46(10):1300 − 1307. [ZHANG Qin,HUANG Guanwen,YANG Chengsheng. Precision space observation technique for geological hazard monitoring and early warning[J]. Acta Geodaetica et Cartographica Sinica,2017,46(10):1300 − 1307. (in Chinese with English abstract) DOI: 10.11947/j.AGCS.2017.20170453 ZHANG Qin, HUANG Guanwen, YANG Chengsheng. Precision space observation technique for geological hazard monitoring and early warning[J]. Acta Geodaetica et Cartographica Sinica, 2017, 46（10）: 1300-1307. （in Chinese with English abstract） DOI: 10.11947/j.AGCS.2017.20170453
[5] 刘文,王猛,朱赛楠,等. 基于光学遥感技术的高山极高山区高位地质灾害链式特征分析—以金沙江上游典型堵江滑坡为例[J]. 中国地质灾害与防治学报,2021,32(5):29 − 39. [LIU Wen,WANG Meng,ZHU Sainan,et al. An analysis on chain characteristics of highstand geological disasters in high mountains and extremely high mountains based on optical remote sensing technology:A case study of representative large landslides in upper reach of Jinsha River[J]. The Chinese Journal of Geological Hazard and Control,2021,32(5):29 − 39. (in Chinese with English abstract) LIU Wen, WANG Meng, ZHU Sainan, et al. An analysis on chain characteristics of highstand geological disasters in high mountains and extremely high mountains based on optical remote sensing technology: a case study of representative large landslides in upper reach of Jinsha River[J]. The Chinese Journal of Geological Hazard and Control, 2021, 32（5）: 29-39. （in Chinese with English abstract）
[6] XIE Mingli,ZHAO Weihua,JU Nengpan,et al. Landslide evolution assessment based on InSAR and real-time monitoring of a large reactivated landslide,Wenchuan,China[J]. Engineering Geology,2020,277:105781. DOI: 10.1016/j.enggeo.2020.105781
[7] 王洵,周云,孙蒙,等. 青海玉树Mw6.9级地震震源破裂过程[J]. 地质通报,2014,33(4):517 − 523. [WANG Xun,ZHOU Yun,SUN Meng,et al. Rupture process of the 2010 Mw6.9 Yushu earthquake in Qinghai Province[J]. Geological Bulletin of China,2014,33(4):517 − 523. (in Chinese with English abstract) WANG Xun, ZHOU Yun, SUN Meng, et al. Rupture process of the 2010 Mw6.9 Yushu earthquake in Qinghai Province[J]. Geological Bulletin of China, 2014, 33(4): 517-523. （in Chinese with English abstract）
[8] 刘星洪,姚鑫,周振凯,等. 滑坡灾害InSAR应急排查技术方法研究[J]. 地质力学学报,2018,24(2):229 − 237. [LIU Xinghong,YAO Xin,ZHOU Zhenkai,et al. Study of the technique for landslide rapid recognition by InSAR[J]. Journal of Geomechanics,2018,24(2):229 − 237. (in Chinese with English abstract) LIU Xinghong, YAO Xin, ZHOU Zhenkai, et al. Study of the technique for landslide rapid recognition by insar[J]. Journal of Geomechanics, 2018, 24(2): 229-237. （in Chinese with English abstract）
[9] 张佳佳,田尤,陈龙,等. 澜沧江昌都段滑坡发育特征及形成机制[J]. 地质通报,2021,40(12):2024 − 2033. [ZHANG Jiajia,TIAN You,CHEN Long,et al. Development and formation mechanism of landslides along Changdu section of Lancang River[J]. Geological Bulletin of China,2021,40(12):2024 − 2033. (in Chinese with English abstract) ZHANG Jiajia, TIAN You, CHEN Long, et al. Development and formation mechanism of landslides along Changdu section of Lancang River[J]. Geological Bulletin of China, 2021, 40(12): 2024-2033. （in Chinese with English abstract）
[10] 许军强,马涛,卢意恺,等. 基于SBAS-InSAR技术的豫北平原地面沉降监测[J]. 吉林大学学报(地球科学版),2019,49(4):1182 − 1191. [XU Junqiang,MA Tao,LU Yikai,et al. Land subsidence monitoring in North Henan plain based on SBAS-InSAR technology[J]. Journal of Jilin University (Earth Science Edition),2019,49(4):1182 − 1191. (in Chinese with English abstract) XU Junqiang, MA Tao, LU Yikai, et al. Land subsidence monitoring in North Henan plain based on SBAS-InSAR technology[J]. Journal of Jilin University (Earth Science Edition), 2019, 49(4): 1182-1191. （in Chinese with English abstract）
[11] FRUNEAU B,ACHACHE J,DELACOURT C. Observation and modelling of the Saint-Étienne-de-Tinée landslide using SAR interferometry[J]. Tectonophysics,1996,265(3/4):181 − 190.
[12] MADSEN S N,ZEBKER H A,MARTIN J. Topographic mapping using radar interferometry:processing techniques[J]. IEEE Transactions on Geoscience and Remote Sensing,1993,31(1):246 − 256. DOI: 10.1109/36.210464
[13] 邓云凯,禹卫东,张衡,等. 未来星载SAR技术发展趋势[J]. 雷达学报,2020,9(1):1 − 33. [DENG Yunkai,YU Weidong,ZHANG Heng,et al. Forthcoming spaceborne SAR development[J]. Journal of Radars,2020,9(1):1 − 33. (in Chinese with English abstract) DENG Yunkai, YU Weidong, ZHANG Heng, et al. Forthcoming spaceborne SAR development[J]. Journal of Radars, 2020, 9（1）: 1-33. （in Chinese with English abstract）
[14] FERRETTI A, PRATI C, ROCCA F. Permanent scatterers in SAR interferometry[C]//IEEE Transactions on Geoscience and Remote Sensing. January 2001, IEEE, 2002: 8 − 20.
[15] BERARDINO P,FORNARO G,LANARI R,et al. A new algorithm for surface deformation monitoring based on small baseline differential SAR interferograms[J]. IEEE Transactions on Geoscience and Remote Sensing,2002,40(11):2375 − 2383. DOI: 10.1109/TGRS.2002.803792
[16] 姚佳明,姚鑫,陈剑,等. 基于InSAR技术的缓倾煤层开采诱发顺层岩体地表变形模式研究[J]. 水文地质工程地质,2020,47(3):135 − 146. [YAO Jiaming,YAO Xin,CHEN Jian,et al. A study of deformation mode and formation mechanism of a bedding landslide induced by mining of gently inclined coal seam based on InSAR technology[J]. Hydrogeology & Engineering Geology,2020,47(3):135 − 146. (in Chinese with English abstract) DOI: 10.16030/j.cnki.issn.1000-3665.201903072 YAO Jiaming, YAO Xin, CHEN Jian, et al. A study of deformation mode and formation mechanism of a bedding landslide induced by mining of gently inclined coal seam based on InSAR technology[J]. Hydrogeology & Engineering Geology, 2020, 47（3）: 135-146. （in Chinese with English abstract） DOI: 10.16030/j.cnki.issn.1000-3665.201903072
[17] 赵富萌,张毅,孟兴民,等. 基于小基线集雷达干涉测量的中巴公路盖孜河谷地质灾害早期识别[J]. 水文地质工程地质,2020,47(1):142 − 152. [ZHAO Fumeng,ZHANG Yi,MENG Xingmin,et al. Early identification of geological hazards in the Gaizi valley near the Karakoran highway based on SBAS-InSAR technology[J]. Hydrogeology & Engineering Geology,2020,47(1):142 − 152. (in Chinese with English abstract) ZHAO Fumeng, ZHANG Yi, MENG Xingmin, et al. Early identification of geological hazards in the Gaizi valley near the Karakoran Highway based on SBAS-InSAR technology[J]. Hydrogeology & Engineering Geology, 2020, 47（1）: 142-152. （in Chinese with English abstract）
[18] ZHANG Y,MENG X M,DIJKSTRA T A,et al. Forecasting the magnitude of potential landslides based on InSAR techniques[J]. Remote Sensing of Environment,2020,241:111738. DOI: 10.1016/j.rse.2020.111738
[19] 季灵运,王庆良,崔笃信,等. 利用SBAS-DInSAR技术提取腾冲火山区形变时间序列[J]. 大地测量与地球动力学,2011,31(4):149 − 153. [JI Lingyun,WANG Qingliang,CUI Duxin,et al. Time series of deformation in Tengchong volcanic area extracted by SBAS-DInSAR[J]. Journal of Geodesy and Geodynamics,2011,31(4):149 − 153. (in Chinese with English abstract) DOI: 10.3969/j.issn.1671-5942.2011.04.034 JI Lingyun, WANG Qingliang, CUI Duxin, et al. Time series of deformation in Tengchong volcanic area extracted by sbas-dinsar[J]. Journal of Geodesy and Geodynamics, 2011, 31（4）: 149-153. （in Chinese with English abstract） DOI: 10.3969/j.issn.1671-5942.2011.04.034
[20] 葛伟丽,李元杰,张春明,等. 基于InSAR技术的内蒙古巴彦淖尔市地面沉降演化特征及成因分析[J]. 水文地质工程地质,2022,49(4):198 − 206. [GE Weili,LI Yuanjie,ZHANG Chunming,et al. An attribution analysis of land subsidence features in the city of Bayannur in Inner Mongolia based on InSAR[J]. Hydrogeology & Engineering Geology,2022,49(4):198 − 206. (in Chinese with English abstract) GE Weili, LI Yuanjie, ZHANG Chunming, et al. An attribution analysis of land subsidence features in the city of Bayannur in Inner Mongolia based on InSAR[J]. Hydrogeology & Engineering Geology, 2022, 49（4）: 198-206. （in Chinese with English abstract）
-
期刊类型引用(3)
1. 李丽华，刘文，白玉霞，王翠英，李双琴. 椰壳纤维-石灰协同作用改良黏土性能试验研究. 水文地质工程地质. 2025(01): 130-140 . 
百度学术
2. 王子杰，郭欢，唐丽霞，吴文丽，廖拉拉. 火棘与紫穗槐根表皮特征对根系拉拔摩擦力的影响. 南京林业大学学报(自然科学版). 2025(03): 181-189 . 百度学术
3. LIANG Shen，WANG Shu，LIU Yabin，PANG Jinghao，ZHU Haili，LI Guorong，HU Xiasong. Mechanism underlying the uprooting of taproot-type shrub species in the loess area of northeastern Qinghai-Xizang Plateau, China. Journal of Arid Land. 2024(10): 1426-1443 . 必应学术
其他类型引用(0)
计量
- 文章访问数: 1408
- HTML全文浏览量: 1317
- PDF下载量: 354
- 被引次数: 3
邮件订阅
RSS