均质边坡离心试验锚杆的应力响应及布设探究_晏长根.pdf

第 45 卷 第 6 期 煤田地质与勘探 Vol. 45 No.6 2017 年 12 月 COAL GEOLOGY P642.22 文献标识码A DOI 10.3969/j.issn.1001-1986.2017.06.016 Exploration of stress responses and arrangement of fully grouted bolts based on centrifugal model test of homogeneous rock slope YAN Changgen1, SUN Weifeng1, YUAN Tong2, GU Liangjun3, GUO Zhiqi4, SHI Yuling5 1. School of Highway, Chang’an University, Xi’an 710064, China; 2. Guangxi Communications Planning Surveying and Designing Institute, Nanning 530029, China; 3. Ningbo Zhengxin Detection Technology Co. Ltd., Ningbo 315105, China; 4. Ningbo Communication Engineering Construction Group Co. Ltd., Ningbo 315105, China; 5. School of Geo- logical Engineering and Geometrics, Chang’an University, Xi’an 710054, China Abstract Using low-carbon steel sheets as fully grouted bolts, centrifugal model test of homogeneous rock slope was carried out and stress responses of fully grouted bolts were studied. Based on the stress responses of fully grouted bolts from the test and practical deation and failure law of homogeneous rock slope, the preferred po- sition for single bolt stress monitoring of one homogeneous rock slope was also explored by using finite-difference numerical simulation technology and a suggested local strength reduction on relaxation region. The results are as follows deation and failure of the homogeneous rock slope started at the shear failure at the toe of the slope, earlier and obvious stress response of the bolt close to the toe of the slope will be triggered and the stress growth rate of the bolt will also continue to increase with the deation and failure development of the slope; the closer the distance between stress monitoring point installed in one bolt and potential sliding surface plastic de- ation zone of the slope, the earlier the time of the bolt stress response and the larger the bolt stress is. It is suggested that the monitoring bolt be put on the slope with 0.250.35 H H is slope height from the toe of homo- geneous rock slope, and concrete analysis must be done for a concrete problem. Results from this paper can provide some references to some extent for the bolts monitoring work of homogeneous rock slope. Keywords stress responses; centrifugal model test; homogeneous rock slope; relaxation region; strength reduction; fully grouted bolts; monitoring; arrangement ChaoXing 第 6 期 晏长根等 均质边坡离心试验锚杆的应力响应及布设探究 97 在复杂的地质条件作用下，边坡会演变为滑坡等 地质灾害，并因此而或多或少地导致一定的损失[1-5]。 为了边坡的防灾、减灾及治灾，避免不必要的损失， 开展必要的边坡监测预警研究一直是岩土工程领域 的热点。 传统的边坡监测手段为位移变形监测方法， 科学技术的不断进步衍生了多种多类的位移监测仪 器[5-8]和预警模型[9-10]，长期工程应用也积累了丰富 的位移监测经验。但实际情况表明，边坡的变形是 边坡失稳的必要条件而非充分条件，边坡的失稳是 其内部应力调整失败的产物。对于数量多、单体规 模小且高度较低的线性工程边坡，全长黏结锚杆作 为常见的边坡加固结构，将其作为探测边坡稳定性 的“媒介”，同时构建无线远程监测预警系统，不失 为一种经济、及时且有效的边坡监测方法。有研究 者[11-12]意识到位移监测方法不但花费较大且反馈较 为滞后，为此开展边坡应力监测方法研究探索。刘 祖强等[13]对工程边坡锚杆监测数据进行多元回归分 析，评价岩体趋于稳定。王涛等[14]研究全长黏结锚 杆应力变化具有边坡开挖效应、时效变形效应和温 度效应等特性，坡体的反复变形会一定程度弱化锚 杆的锚固效果及其对边坡时效变形的控制。梁桂兰 等[15]采用 PSO-RBFNN 模型对某实际边坡一测点处 的锚杆应力增长曲线进行拟合分析。何满潮等[7]建 立滑坡摄动力学监测模型，并给出了相应的预警准 则，在矿山边坡等较为重大的工程项目中开展了应 用研究。以上研究重在现场实际应用，本文以低碳 钢片模拟全长黏结锚杆，开展均质岩质边坡离心模 型试验，研究了锚杆的应力响应规律。并结合试验 模型及实际工程，采用有限差分数值模拟技术和局 部强度折减法，研究了均质岩质边坡单锚应力监测 的较佳位置。研究成果可为均质岩质边坡锚杆应力 监测的推广应用提供参考。 1 试验方法 1.1 试验原理 在保持岩性物理力学参数和变形参数与实际边 坡基本近似的前提下，离心模型试验通过增加边坡 模型的加速度，即可达到与实际边坡极为相似的应 力状态，反映边坡变形破坏的全过程[16]。图 1 为一 单位宽度厚的均质岩质边坡模型，试验过程中，潜 在滑体重力为 WWρaV，ρ 为密度，a 为离心加速 度，V 为体积；坡面从上自下布设 k 根全长黏结锚 杆，第 j 根锚杆入射角为 βj，设第 i 根锚杆在潜在滑 面处所能提供的极限轴力为 Fj。基于瑞典圆弧条分 法，则试验过程中边坡安全系数 K 为 111 r s 1 [cossin]tancos sin nkk ijjjjjj ijj n i i ρVaαFβγφcLFβγ F K F ρVaα        式中 αi为潜在滑体第 i 个分条的滑动面法线与竖直 线的夹角； γj为与第 j 根锚杆相交点处的滑动面法向 与竖直线的夹角；φ 为潜在滑动面的内摩擦角；c 为潜在滑动面的黏聚力；L  为潜在滑动面的弧长；n 为潜在滑体的分条数。 图 1 均质岩质边坡离心模型示意图 Fig.1 Schematic diagram of centrifugal model of homogene- ous rock slope 试验不断加载时，式中的唯一变量离心加速度 a 不断增加，滑动力 Fs以比抗滑力 Fr更大的速率增 长，直至边坡模型趋近临界状态而产生变形破坏。 1.2 试验内容 以低碳钢片模拟全长黏结锚杆，开展均质岩质 边坡离心模型试验，研究边坡直至出现变形破坏迹 象过程中的锚杆应力响应规律。 采用长安大学 TLJ-3 型土工离心机，最大容载 60g t，有效半径 2.0 m， 模型箱尺寸 50 cm36 cm40 cm。图 2 沿坡面由上 自下布置 1 号、2 号和 3 号监测锚杆，各锚杆长度 为 8 cm，1 号锚杆上布设应变计监测点 1、2 和 3， 2 号锚杆上布设应变计监测点 4、5、6，3 号锚杆上 布设应变计监测点 7、8 和 9。 1.3 模型材料及制作 a. 锚杆 为了便于布设应变计，采用激光精密切割横截 面积为 10 mm1 mm、长为 8 cm 的低碳钢片模拟全 长黏结锚杆，钢片弹性模量为 200 GPa，其表面按 ChaoXing 98 煤田地质与勘探 第 45 卷 图 2 均质边坡模型及监测锚杆布置图单位cm Fig.2 Layout of monitoring bolts and homogeneous slope model 图 2 布设相应的应变计。钢片先采用细纱布打磨除 锈，再用无水乙醇清洗干燥，接着将应变计用薄薄 的 AB 胶黏接固定， 后再于钢片表层涂抹薄薄的 AB 胶并均匀撒上细沙，用以防水和增加锚杆与岩石的 黏结效果。边坡模型制作时，将锚头端钢片插入设 计的模板坡面板方孔横截面积为 10 mm1 mm中， 且其孔轴线方向与坡面板法向一致。 b. 边坡模型 边坡模型材料为重晶石粉、石英砂、石膏和水， 配比为 110.270.25， 用以模拟均质岩质边坡。 边坡 模型物理力学参数见表 1。在模型箱内侧均匀涂抹 凡士林并粘贴保鲜膜以尽量消除边界摩阻效应。将 模型箱透明玻璃板侧外法向竖直向上， 卸除玻璃板， 插入设计的边坡模型制作模板，安装锚杆定位后， 少量多份地拌制模型材料进行连续浇筑，必要时采 用细木棒捣密。 表 1 均质边坡模型物理力学参数 Table 1 Mechanical properties of homogeneous slope model 密度/kgm-3 泊松比 弹性模量/GPa 黏聚力/kPa 内摩擦角/ 2 500 0.26 0.5 43 32 1.4 测量及观测方法 a. 锚杆应变监测 锚杆上布设的应变计通过线 缆连接至离心机传感器面板，室内电脑配置的软件 可实时接收应变监测数据。 整理锚杆应变监测数据， 根据胡克定律可计算锚杆的应力值。 b. 摄像观测 试验时，高速摄像机透过模型箱 侧方透明玻璃板可进行摄像，用以观察边坡剖面有 无显著变形破坏迹象。 1.5 试验加载 模型制作完成， 调试应变计及摄像机工作正常， 离心加速度 20g→40g→60g→80g→100g→120g 进 行分级试验加载，加载过程见表 2。 表 2 试验分级加载过程 Table 2 Graded loading process of test 离心加速度 a/ms-2 开始加载时刻/s 到达加载时刻/s 20g 0 270 40g 380 540 60g 970 1 130 80g 1 360 1 480 100g 1 800 1 910 120g 2 140 2 310 注120g 维持至试验运行 2 790 s 时开始减速停止。 2 试验结果与分析 2.1 试验现象 试验加载至约 120g 时，坡脚出现小范围的剪切 破坏现象，但边坡的整体稳定性尚好，说明坡脚应力 较为集中， 边坡变形破坏起始于坡脚的局部剪切破坏。 2.2 锚杆应变监测结果 由胡克定律计算 3 根锚杆各监测点的应力历时 曲线图 3，图中横坐标为试验加载的累积时间。由 图 3a 可知1 号锚杆的监测点 3 约 400 s 时产生了 明显的应力响应，而监测点 2 和 3 则于约 750 s 时 产生明显的应力响应；3 个监测点产生明显的应力 响应后，监测点 3 处的锚杆应力值始终是 3 个监测 点中最大值；监测点 1、2 和 3 处的锚杆应力值最终 都趋于收敛值， 分别约为 9 MPa、 13 MPa 和 23 MPa。 由图 3b 可知2 号锚杆的 3 个监测点约于 685 s 同 时产生明显的应力响应，之后监测点 6 处的锚杆应 力值始终是 3 个监测点中的最大值；随着试验的不 断加载，监测点 4、5 和 6 表现不断增长的趋势。由 图 3c 可知3 号锚杆的监测点 8 约 400 s 时产生了 明显的应力响应，而监测点 7 和 9 则于约 685 s 产 生了明显的应力响应；约 6851 685 s 范围内，监测 点 9 处的锚杆应力值略大于监测点 7， 之后监测点 8 处的锚杆应力值始终是 3 个监测点中的最大值；随 着试验的不断加载，监测点 7 最终表现趋于收敛的 趋势，而监测点 8 和 9 则表现出不断增长的趋势。 对比图 3a3c 可知3 号锚杆上的监测点 8 开 始产生明显应力响应的时刻不但最早，并且其监测 锚杆应力增长速率也基本是 3 根锚杆各监测点中最 大的。这表明对于均质岩质边坡，为尽早捕捉边坡 变形破坏的发展趋势，在靠近坡脚区域布设监测锚 杆是比较适宜的。 ChaoXing 第 6 期 晏长根等 均质边坡离心试验锚杆的应力响应及布设探究 99 图 3 各锚杆各监测点应力历时曲线 Fig.3 Duration curves of stress of every bolt at different points 2.3 数值模拟分析 为什么监测点 3、监测点 6 和监测点 8 处的应 力在相应的锚杆上应力响应时间较早且量值最大 采用 FLAC3D进行数值模拟分析。不考虑锚杆的作 用，建立图 2 所示且走向为 0.1 cm 的边坡模型，按 平面应变问题考虑，模型的边界条件与图 1 相同， 模型参数见表 1。先采用“弹性求解法”[17]进行模型 初始地应力生成，再采用摩尔库伦本构模型计算， 划分网格边长约 0.25 cm， 计算终止条件最大不平衡 力比率 ratio 小于 110-5。 不考虑锚杆的作用时，分别模拟模型在离心试 验 30g 和 80g 工况，查看相应工况下边坡的剪切应 变增量云图图 4， 并在图中标示了 9 个监测点的位 置。由图可知试验加载的过程，边坡模型的变形 破坏始于坡脚的剪切破坏；随着边坡潜在滑动面的 开展，锚杆监测点 3 和监测点 6 靠近潜在滑动面， 而监测点 7 正好位于潜在滑动面上。边坡潜在滑动 面相当于塑性变形区，研究表明[18-19]塑性变形区是 锚杆应力峰值所在区域，从而解释了监测点 3 和监 测点 6 处的应力分别是 1 号锚杆和 2 号锚杆上应力 响应较早且量值最大的原因。而监测点 7 虽正好位 于潜在滑动面上，但因其靠近坡面，且伴随着坡脚 剪切破坏的开展，3 号锚杆锚头部分出现解耦，故 监测点 8 处的应力是 3 号锚杆上应力响应时间较早 且量值最大的。因此，实际边坡监测时，宜在靠近 边坡潜在滑动面塑性变形区处的锚杆上布设应力 监测计。 3 锚杆较佳监测布设位置探讨 根据试验结果得知在靠近坡脚区域布设监测锚 杆是比较适宜的，但具体布设在什么位置有待于进 一步深入研究。实际工程边坡变形破坏演变与离心 模型试验的原理有所差异，采用下面的方法对锚杆 较佳监测布设位置进行探讨。将图 2 的边坡断面线 性比尺扩大 40 倍得图 5 所示的分析模型， 考虑锚杆 图 4 两种离心加速度下模型剪切应变增量云图 Fig.4 Cloud chart of shear strain increment of model under two centrifugal accelerations 的影响范围在边坡走向取 3 m，模型的边界条 件与图 1 相同。模型参数见表 1，另假设边坡岩石 为软岩，饱和单轴抗压强度取 15 MPa。先采用“弹 性求解法”[17]进行模型初始地应力生成， 再采用摩尔 库伦本构模型计算，划分网格边长约 0.50 m，计算 终止条件 ratio 小于 110-5。锚杆入射角为 10，长 度取 9 m，钻孔直径取 75 mm，采用 FLAC3D中的 cable 单元[17]模拟锚杆，cable 单元参数见表 3。 研究表明[20-21]工程边坡开挖的过程中存在一定 的开挖松弛区。而整体强度折减法存在计算的塑性 ChaoXing 100 煤田地质与勘探 第 45 卷 图 5 锚杆监测布设位置探讨分析模型单位m Fig.5 Analysis model of bolt position 表 3 Cable 单元参数 Table 3 Parameters of Cable element 参数 Emod/ GPa gr_coh/ MPam-1 gr_fric / gr_per /mm gr_k/ MPa Xrea/ mm2 Ytens /kN 数值 210 1.213 50 235.5 8.24 615.44246 区偏大，且与实际的边坡失稳破坏并不相符[22]。确 定开挖松弛区下界比较困难，此处根据没有锚杆时 边坡整体强度折减法计算后的 x 向位移量级一致 性，较保守地推测图 5 中的边坡开挖松弛区的下界 符合三次多项式曲线 32 0.00550.1151.30350.4777zxxx， 6.0715.74x≤≤。其中，图 5 的水平向和竖直向分 别为 x 轴和 z 轴，坐标原点在边坡模型的左下角处。 对于开挖完本来稳定的均质岩质边坡，其开挖 松弛区最容易受外界温度变化、降雨及风化等作用 的影响，为松弛区岩体质量不断劣化而变形破坏提 供了条件。故探讨将单个全长黏结锚杆分别布设在 距离坡脚 ζ HH 为坡高，ζ 取值 0.1、0.2、0.3、0.4、 0.5 和 0.6的坡面上，编制 fish 程序仅对开挖松弛区 下界以上岩体的抗剪强度参数进行折减，采用 FLAC3D中的 history 命令监测锚杆最大拉力值， 得锚 杆最大拉力随边坡局部折减系数的变化曲线图 6。 由图 6 可知随着边坡局部折减系数的不断增 加，布设在不同位置的锚杆最大拉力值由线性变化 阶段急剧转变为加速变化阶段，锚杆应力响应比较 显著；锚杆布设在距离坡脚 0.3H 坡面上时，最大拉 力值线性变化阶段的速率和加速变化阶段起始段的 速率都是最大的。这表明当外界不利作用引起边坡 开挖松弛区岩体抗剪强度参数不断劣化时，在距离 上述边坡坡脚 0.3H 的坡面处布设监测锚杆， 其最大 拉力值及其变化速率均最显著，更容易引起监测人 员的注意与警惕，以便提前预警边坡潜在风险。 图 6 锚杆最大拉力随边坡局部折减系数的变化曲线 Fig.6 Variation curves of bolt’s maximum tension vs local reduction factor of slope 笔者基于离心试验监测锚杆应力响应规律的启 发，并考虑到一般实际边坡监测对象尺寸一定，锚 杆灌浆质量和入射角规范[23]都有相应规定，均质岩 质边坡不同坡高岩体的应力应变存在差异性，而外 界影响可能不断劣化边坡松弛区岩体质量而造成变 形破坏。故采用上述控制变量法，数值模拟分析了 监测锚杆布设在不同坡高坡面位置时的应力响应规 律， 探讨监测锚杆布设的较佳位置位于 0.3H 的坡面 处。但应当注意监测锚杆最大拉力及其变化速率 与锚杆入射角、锚杆布设位置、灌浆质量、外界 影响、边坡尺寸、岩体性质及其演变等诸多因素 相互关联[19,24-26]。对于特定的监测边坡对象，要想 确定更为准确的锚杆监测布设位置，有待具体问题 具体分析， 上述的研究探讨方法可作参考。 故根据上 述锚杆较佳监测布设位置探讨结果， 又考虑到施工的 方便性，建议在距离坡脚约 0.250.35H 高度的坡面 处布设监测锚杆， 为进一步准确把握均质岩质边坡锚 杆监测的合理监测位置提供了一定参考借鉴。 4 结 论 a. 均质岩质边坡变形破坏起始于坡脚的剪切 破坏，且随着边坡变形破坏的开展，靠近坡脚的全 长黏结锚杆会较早产生明显的应力响应，锚杆应力 增长速率会不断增加。 b. 锚杆应力监测点距离均质岩质边坡潜在滑 动面塑性变形区越近，锚杆应力响应时间越早且 量值越大。实际边坡监测时，宜在靠近边坡潜在滑 动面塑性变形区处的锚杆上布设应力监测计。 c. 基于试验监测锚杆的响应规律及实际均质 岩质边坡演化的客观规律，提出仅针对开挖松弛 区进行局部折减，数值模拟分析确定监测锚杆相 对于边坡高度的较佳布设位置。兼顾锚杆受力的 影 响 因 素 及 施 工 方 便 性 ， 建 议 在 距 离 坡 脚 约 0.250.35H 高度的坡面处布设监测锚杆，为进一 ChaoXing 第 6 期 晏长根等 均质边坡离心试验锚杆的应力响应及布设探究 101 步准确把握均质岩质边坡锚杆监测布设位置提供 一定的参考借鉴。 d. 监测锚杆最大拉力及其变化速率与锚杆入 射角、锚杆布设位置、灌浆质量、外界影响、边坡 尺寸、岩体性质及其演变等诸多因素相互关联，对 于特定的监测边坡对象，要想确定更为准确的锚杆 监测布设位置，提倡具体问题具体分析。 参考文献 [1] 廖安杰，李天斌，薛德敏，等. 张家沟滑坡治理工程及数值模 拟[J]. 煤田地质与勘探，2016，445126–130. LIAO Anjie，LI Tianbin，XUE Demin，et al. Engineering treatment and numerical simulation of Zhangjiagou landslide[J]. Coal Geology Exploration，2016，445126–130. [2] 韩琨. 青藏高原山区某公路滑坡滑动机理及其稳定性分析[J]. 煤田地质与勘探，2013，41163–66. HAN Kun. Landslide sliding mechanism and stability analysis of highway in Qinghai-Tibet mountainous area[J]. Coal Geology Exploration，2013，41163–66. [3] 李元雄. 随机地震作用下滑坡稳定性可靠度分析[J]. 煤田地 质与勘探，2010，38351–54. LI Yuanxiong. Study on reliability of slope under the random earthquake effect[J]. Coal Geology Exploration， 2010， 383 51–54. [4] 雍睿， 朱志明， 邹宗兴， 等. 万塘滑坡治理效果预测评价[J]. 煤 田地质与勘探，2012，40155–59. YONG Rui， ZHU Zhiming， ZOU Zongxing， et al. uation of Wantang landslide treatment effect[J]. Coal Geology Explora- tion，2012，40155–59. [5] 李胜，齐嘉义，胡海永，等. 露天矿边坡自动化监测关键技术 研究[J]. 煤田地质与勘探，2016，446124–128. LI Sheng，QI Jiayi，HU Haiyong，et al. Research on key technology of automatic monitoring of surface mine slope[J]. Coal Geology Exploration，2016，446124–128. [6] DUNNING S A， ROSSER N J， MASSEY C I. The integration of terrestrial laser scanning and numerical modeling in landslide investigations[J]. Quarterly Journal of Engineering Geology and Hydrogeology，2010，432233–247. [7] 董文文，朱鸿鹄，孙义杰，等. 边坡变形监测技术现状及新进 展[J]. 工程地质学报，2016，2461088–1095. DONG Wenwen， ZHU Honghu， SUN Yijie， et al. Current status and new progress on slope deation monitoring technolo- gies[J]. Journal of Engineering Geology ， 2016 ， 246 1088–1095. [8] XU J B， YAN C G， ZHAO X， et al. Monitoring of train-induced vibrations on rock slopes[J]. International Journal of Distributed Sensor Networks，2017，1311–6. [9] 郭云开，姚瑶，苟叶培，等. 山区高速公路边坡变形监测模型[J]. 长沙理工大学学报自然科学版，2014，11142–47. GUO Yunkai，YAO Yao，GOU Yepei，et al. High slope de- ation monitoring model in mountainous area[J]. Journal of Changsha University of Science and TechnologyNatural Sci- ence，2014，11142–47. [10] 张德成. 滑坡预测预报研究以丁家坟滑坡为例[D]. 昆 明昆明理工大学，2014. [11] 何满潮. 滑坡地质灾害远程监测预报系统及其工程应用[J]. 岩石力学与工程学报，2009，2861081–1089. HE Manchao. Real-time remote monitoring and forecasting sys- tem for geological disasters of landslides and its engineering ap- plication[J]. Chinese Journal of Rock Mechanics and Engineer- ing，2009，2861081–1089. [12] 郭永建， 王少飞， 李文杰. 应力监测在公路岩质边坡中的应用 研究[J]. 岩土力学，2013，3451397–1342. GUO Yongjian， WANG Shaofei， LI Wenjie. Application of stress monitoring to highway rock slope[J]. Rock and Soil Mechanics， 2013，3451397–1342. [13] 刘祖强， 廖勇龙， 朱全平. 三峡永久船闸一期高边坡监测锚杆 应力分析[J]. 长江科学院院报，2004，21140–44. LIU Zuqiang， LIAO Yonglong， ZHU Quanping. Analysis of first stage monitoring anchor bar stress on high slope of TGP perma- nent lock during operation[J]. Journal of Yangtze River Scientific Research Institute，2004，21140–44. [14] 王涛， 朱焕春， 邰淑彩. 三峡临时船闸高边坡系统锚杆应力回 归分析[J]. 人民长江，2001，3257–9. WANG Tao， ZHU Huanchun， TAI Shucai. Bolt stress regression on high slope of TGP Ship lock[J]. Yangtze River，2001，325 7–9. [15] 梁桂兰，徐卫亚，何育智，等. PSO–RBFNN 模型及其在岩土 工程非线性时间序列预测中的应用[J]. 岩土力学，2008， 292995–1000. LIANG Guilan，XU Weiya，HE Yuzhi，et al. Study and appli- cation of PSO-RBFNN model to nonlinear time series forecast- ing for geotechnical engineering[J]. Rock and Soil Mechanics， 2008，292995–1000. [16] 程展林，龚壁卫，李波，等. 砂土边坡稳定性离心模型试验研 究[J]. 水利学报，2015，46增刊 112–17. CHENG Zhanlin，GONG Biwei，LI Bo，et al. Centrifuge modelling of the stability of the sandy soil slope[J]. Journal of Hydraulic Engineering，2015，46S112–17. [17] 陈育民，徐鼎平. FLAC/FLAC3D 基础与工程实例[M]. 北京 中国水利水电出版社，2009. [18] 胡玉银. 利用锚杆轴力量测定围岩破坏区的方法[J]. 工程地 质学报，1995，3367–72. HU Yuyin. for determining failure zone of surrounding rock by means of axial stress measurement of bolt[J]. Journal of Engineering Geology，1995，3367–72. [19] 伍国军， 陈卫忠， 王永刚. 基于岩体蠕变效应的锚杆应力分布 及其变化规律研究[J]. 岩土力学， 2010， 31增刊 1 150–155. 下转第 110 页 ChaoXing 110 煤田地质与勘探 第45卷 133519–524. [10] 周文生， 吴振宇， 刘海燕. 无人机遥感在矿山地质环境调查中 的应用[J]. 地下水，2014，362128–129. ZHOU Wensheng，WU Zhenyu，LIU Haiyan. Application of UAV remote sensing in the investigation of mine geological en- vironment[J]. Groundwater，2014，362128–129. [11] 魏长婧，汪云甲，王坚，等. 无人机影像提取矿区地裂缝信息 技术研究[J]. 金属矿山，2012，411090–92. WEI Changjing， WANG Yunjia， WANG Jian， et al. Research on the ination technology of mining area in the mining area of UAV image extraction[J]. Metal Mine，2012，411090–92. [12] 赵星涛，胡奎，卢晓攀，等. 无人机低空航摄的矿山地质灾害 精细探测方法[J]. 测绘科学，2014，39649–52. ZHAO Xingtao，HU Kui，LU Xiaopan，et al. The fine detection of mine geological hazards in low altitude aerial photog- raphy[J]. Science of Surveying and Mapping，2014，396 49–52. [13] 张启元. 无人机航测技术在青藏高原地质灾害调查中的应 用[J]. 青海大学学报自然科学版，2015，33267–72. ZHANG Qiyuan. The UAV aerial technology in Qinghai Tibet Pla- teau geological disaster investigation application[J]. Journal of Qing- hai UniversityNatural Science Edition，2015，33267–72. [14] 镕李冰，刘源，刘素红，等. 基于低空无人机遥感的冬小麦覆 盖度变化监测[J]. 农业工程学报，2012，2813160–165. LI Bing，LIU Rongyuan，LIU Suhong，et al. Monitoring vegetation coverage variation of winter wheat by low-altitude UAV remote sensing system[J]. Transactions of the Chinese Society of Agricultural Engineering， 2012， 2813 160–165. [15] 毕凯，李英成，丁晓波，等. 轻小型无人机航摄技术现状及发 展趋势[J]. 测绘通报，2015327–31. BI Kai，LI Yingcheng，DING Xiaobo，et al. Aerial photo- grammetric technology of light small UAVStatus and trend of development[J]. Bulletin of S