Preliminary analysis on rainfall thresholds for early warning of the rainfall induced landslides based on “step” deformation characteristics
-
摘要: 由于滑坡灾害的发生具有随机性和复杂性，经常给人们生活造成严重的经济损失，甚至威胁生命安全，故建立滑坡预警判据及模型是减小损失的重要手段。文章首先确定了有效降雨量-变形速率-改进切线角为降雨型阶跃滑坡的预警关键判据。然后运用斜率单变点法确定滑坡“阶跃”变形曲线上的“破坏拐点”和“稳定拐点”并提取变形加速区间的详细信息，以此求解有效降雨量-变形速率-改进切线角综合预警判据阈值。最后，结合麻池村1号滑坡各发展阶段宏观变形特征，将时空演化规律有机结合，建立了4级递进式综合判据分级预警模型。结果表明：麻池村1号滑坡的运移形式为牵引式，在降雨因素的激发下，累计位移曲线呈典型的“阶跃”变形特征；麻池村1号滑坡变形受前期降雨影响，前期降雨与当期降雨组成一个降雨过程的时间间隔阈值为10 d，“前期+当期降雨”模式下对应的降雨阈值分别为24 mm、32 mm；基于监测数据精细化分析以及改进切线角法获取的麻池村1号滑坡的位移速率阈值为 12 mm/d。研究可为此类型滑坡灾害的预警提供理论基础和技术指导。Abstract: Since the occurrence of landslide disasters is random and complex, which often causes serious economic losses to people's lives and even threatens life safety, establishing landslide early warning criteria and models is an important means to reduce losses. In this paper, the effective rainfall-deformation rate-improved tangent angle is determined as the key criterion for the early warning of rainfall-type step landslide. Then, the slope single-variable point method is applied to determine the “damage inflection point” and “stability inflection point” on the landslide “step” deformation curve to extract detailed information of the deformation acceleration interval, so as to solve the effective The threshold value of rainfall-deformation rate-improved tangent angle integrated warning criterion. Finally, a 4-stage progressive integrated warning model was established by combining the macroscopic deformation characteristics of the landslide at each stage of development and the spatial and temporal evolution laws. The results show that the transport form of landslide No. 1 in Machi Village is traction type, and the cumulative displacement curve shows typical “step” deformation characteristics under the excitation of rainfall factors. The deformation of landslide No. 1 in Machi village is affected by the early rainfall, and the time interval threshold between the early rainfall and the current rainfall is 10 d, and the corresponding rainfall thresholds in the “early + current rainfall” mode are 24 mm and 32 mm respectively. The study can provide theoretical basis and technical guidance for early warning of this type of landslide hazard.
-
Keywords:
- rainfall landslide /
- step deformation /
- pre-rainfall /
- threshold /
- early warning model
-
0. 引言
随着我国经济社会的高速发展，尤其是大规模的城市化进程导致各大城市人地关系紧张。位于广州市西北部的广花盆地分布有大面积的隐伏岩溶，目前区内有大量工程建设项目，岩溶塌陷等问题成为制约城市建设的重大地质因素。查明岩溶分布的范围，摸清岩溶发育的规律，分析岩溶塌陷的致灾因素，探索行之有效的解决办法，已成为岩溶发育区工程建设领域的重大技术难点之一。
近年来，国际上有大量的相关岩溶塌陷实例研究文章发表。DOĞAN等[1]对土耳其Obruk地区的19处岩溶塌陷点展开了相关调查，结果表明，引发当地岩溶塌陷的主要因素是由农业活动抽排地下水，导致水位下降。BENAC等[2]对克罗地亚两个典型岩溶地区的形成条件进行了对比研究，一个是上覆型岩溶地貌，一个是在古河流作用下发育的裸露型岩溶地貌。 HEIDARI等[3]研究了伊朗西部地区的水文地质资料与岩溶发育的关系，结果表明地下水循环强烈影响岩溶的发育，尤其是断裂发育的位置影响更为显著。RONGIERA等[4]模拟了岩溶管道的空间三维分布。DURINGER等[5]研究了东南亚地区的岩溶发育及影响因素。SONG等[6]研究了韩国的岩溶发育，EPTING等[7]研究了瑞士的岩溶与地下水，FIDELIBUS等[8]、COOPER等[9]和NOVEL等[10]也分别对意大利、英国和希腊等国的岩溶地下水及岩溶塌陷实例进行了研究，取得了较多的认识和成果。
国内除广州以外，武汉市和桂林市也是岩溶塌陷高发的城市，胡亚波等[11]以武汉市烽火村塌陷为例分析了岩溶地面塌陷的形成机理，探讨了塌陷的变化过程；罗小杰[12]对武汉市的浅层岩溶发育特征进行了研究并对岩溶塌陷灾害提出了防治建议；陈学军等[13]以广西桂林西城区为例进行了岩溶地区破坏性抽水致塌试验研究；康春景等[14]对岩溶发育特征及其对铁路线性工程的不利影响进行了研究。由此可见针对岩溶高发育地区展开深入研究对增强城市规划、工程设计及施工建设的前瞻性、高效性和安全性具有重大的现实意义。
随着城市用地紧张，近年广花盆地内因岩溶问题引发的工程事故屡见不鲜，已引起了相关专家和学者的高度重视，相关研究也相继展开，如广花盆地南部的金沙洲地区因岩溶发育产生了大量的岩溶塌陷事故[15]，黄健民等[16-17]、郑小战等[18-19]和郭宇等[20-21]等对其进行了大量调查研究，对引发岩溶塌陷地质灾害的成因机理做了详细的研究和分析，认为广花盆地南部金沙洲地区广泛发育的岩溶是地质灾害发生的物质基础，武广高铁金沙洲隧道以及当地一些工程建设抽排地下水导致水位波动是地质灾害发生的诱发因素，相关研究在长期监测地下水位的基础上结合塌陷的发生时间，分析了塌陷形成的内在机理和外部动力，为类似研究积累了宝贵的数据和资料。除金沙洲地区以外，广花盆地南部另一重要地区大坦沙岛也饱受岩溶塌陷的困扰，苏扣林等[22]对广州市荔湾区大坦沙岛的岩溶稳定性进行了评价，提供了该岛的未来规划建设性的建议。广花盆地除南部多发岩溶塌陷外，西北部也是重点多发区，骆荣等[23]对广花盆地西北部的赤坭地区进行了实地调查，并对岩溶塌陷地质灾害的形成机制进行了深入的分析，为当地的地质灾害防治提出了指导性的意见。
广州市地质调查院在2009年对夏茅区进行了应急调查和普查工作，2013年开展了地质灾害详细调查，部署了钻探实物工作量，后期该地块进行了地质灾害治理和原地复建结合并于2019年完工，目前正在开展地质环境监测工作，本文整理了夏茅村岩溶地面塌陷的现场资料和有关数据，从地质灾害的基本情况、岩溶发育特征以及致灾主控因素等方面，深入研究夏茅地区岩溶发育、演化及致塌的机理原因，为其他城市和地区在城市规划建设、运营管理和地质灾害防治等方面提供研究实例以供借鉴。
1. 研究区地质灾害特征
研究区夏茅村位于广州市白云区，广花盆地中南部。东面与均禾、嘉禾、鹤龙街，南与黄石、新市、石井相邻，西面与石门街接壤，北面为白云区江高镇。交通十分便利，东有机场高速公路及广花一级公路，南有黄石西路可达市区，西有广清高速公路，北有华南城市快速路。
2008年12月19日16时，夏茅村向西街5巷18号开展房基勘探，钻孔将地下土洞打穿，受此影响向西街及沙园坊华富街发生大规模岩溶地面塌陷。据现场调查，共发育5个地面塌陷坑，分别位于向西街8巷16号、6巷16号、3巷13号、5巷北端及沙园坊华富街7巷6号（图1），其中向西街6巷16号房低层建筑全部沉入塌陷坑中，沙园坊华富街7巷6号房底部大半沉入塌陷坑中（基础为10～13 m桩基），5巷北端塌陷坑面积最大，达163 m2，其它单个陷坑面积为70.4～86.0 m2，坑内均被水充填，水面距地面约0.5 m，地下水没有流动迹象。
![]()
图 1 夏茅村岩溶地面塌陷坑及地面变形裂缝位置分布图Figure 1. Distribution of karst ground collapse pits and ground deformation cracks in Xiamao Village本次地面塌陷致使8幢房屋直接损毁，数十处房屋墙壁及地面出现变形和开裂，事故未造成人员伤亡，但影响范围较大约5000 m2。根据有关方面的鉴定，受本次岩溶地面塌陷影响损伤的村民自建房有81栋，建筑面积23308 m2。其中38栋房屋受到轻微损坏，其主体承重构件保持完好；15栋住宅楼存在一般性损伤，主体承重构件少量损坏，围护墙体稍有裂缝和地基基础有轻微不均匀沉降现象；2栋房屋承重构件严重损坏，显著影响整体承载功能和使用功能，影响正常使用；18栋住宅房屋属于危房，应立即停止使用；8栋直接损毁。
2. 研究区地质背景
研究区位于白云区石井街夏茅村向西北街与沙园坊华富街交界一带，区内地势较为平坦，地面标高6.5～8.3 m，地貌类型属广花盆地河谷冲积平原地貌，岩性为可溶性石灰岩，地层时代为石炭系壶天群，灰岩为隐晶质结构或微晶结构，厚层状构造，溶洞是区内岩溶发育的主要形态，溶隙发育带和溶沟也较为常见。区内发育北北东向F10和近东西向F12两条断裂构造带，断层带物质为构造角砾岩，岩体较为破碎，断层破碎带及影响带为地下水运移提供了良好的通道，同时沿断层破碎带灰岩的溶蚀作用强烈，溶洞、溶沟及溶槽等岩溶形态较为发育。
3. 岩溶发育特征
研究区属于广花断陷盆地，形成于白垩纪以前，受广州−从化活动深断裂控制。广花复式向斜的可溶性碳酸盐岩受褶皱和断裂构造的影响，在特定的地质历史时期，由于地表水和地下水的侵蚀溶蚀作用，形成了溶洞、溶蚀沟槽、岩溶洼地、岩溶漏斗甚至岩溶地下河。后期受构造沉降作用，原来位于地表的岩溶地貌被上覆土体覆盖而隐埋于地下，形成隐伏岩溶，隐伏岩溶在地下水动力条件的作用下，溶洞进一步发育，第四系松散覆盖层中伴生的土洞也极易出现顶板坍塌，从而造成岩溶地面塌陷。
3.1 基岩面发育特征
区内岩溶基岩面起伏剧烈，高程变化较大，基岩面发育的形态各异，有溶沟、溶槽、石芽、漏斗和溶脊等（表1，图2）。据钻孔资料显示，基岩面的溶蚀纹沟宽度6～8 cm，深度3～5 cm。基岩面的高程变化较大，相对高差一般4.3～10 m，最大高差达22.3 m，最大坡度高达70°，局部地段存在突变性溶沟陡坎，如BZK2与BZK2−1两孔之间水平距离不足1 m，岩面高差可达3 m，BZK8与BZK8−1两孔之间水平距离约0.6 m，而高差竟高达7.7 m。此外，在XMZK4～XMZK1−1（四巷与五巷交界部位）两孔之间还有一条溶槽，宽5.0～12.1 m，深7.8～11.7 m，长约13.0 m，向北北东方向延伸，其东侧存在一条溶脊，位于BZK5～BZK12（三巷与四巷交界部位），溶脊走向北北东，脊背宽约5～8 m，两侧基岩面坡度陡峭，约为50°～70°。
表 1 夏茅村岩溶地面塌陷区钻孔揭露基岩面统计表Table 1. Statistics of bedrock surface exposed by boreholes in karst ground collapse area of Xiamao Village序号 钻孔编号 X Y 基岩标高/m 序号 钻孔编号 X Y 基岩标高/m 1 BZK1 40232.48 36294.55 −8.8 15 BZK13 40230.63 36326.76 −17.56 2 BZK2 40221 36303.48 −17.41 16 BZK14 40220.25 36318.33 −12.6 3 BZK2-1 40221 36304.1 −14.41 17 BZK15 40233.53 36302.07 −10.58 4 BZK3 40226.03 36306.16 −12 18 XM3 40237.16 36341.17 −22.83 5 BZK4 40238.71 36331.12 −24.38 19 XMZK1-1 40236.63 36342.35 −11.13 6 BZK5 40232.83 36343.35 −6.8 20 XM6 40222.02 36322.91 −14.03 7 BZK6 40221.68 36353.88 −15.2 21 XMZK4 40242 36328 −16.58 8 BZK7 40214.63 36355.35 −15.3 22 XMZK6 40226.51 36329.41 −14.55 9 BZK8 40242.65 36405.57 −18.9 23 XMZK7 40216.91 36315.78 −11.6 10 BZK8-1 40242.72 36406.85 −26.6 24 ZK3 40274 36316 −23.3 11 BZK9 40239.29 36404.27 −18.9 25 ZK5 40266 36337 −48 12 BZK10 40236.45 36411.06 −41.2 26 ZK8 40256 36357 −38 13 BZK11 40225.66 36291.71 −13.7 27 BZK16 40258 36301 −16.5 14 BZK12 40226.89 36339.67 −7.4 ![]()
图 2 BZK15孔岩面溶沟形态和BZK13孔岩面溶沟形态Figure 2. Rock surface karst ditch morphology of BZK15 and BZK13 hole根据浅层地震反射数据，塌陷区基岩面起伏变化较大，埋深在6.65～88.74 m ，区内存在一个基岩凹陷区，深度范围在30～88.74 m，面积约4066.91 m2，走向北西，中部有向北东方向延伸的趋势，断裂F10恰好从基岩凹陷区中部穿过，使其成为最复杂的局部地带。基岩凹陷区的西南部和东部相对平缓，深度分别在16～30 m和20～35 m之间变化。
3.2 溶洞发育特征
2013年钻探施工的22个钻孔中，19个揭露到溶洞，见洞率高达86%，总计发现溶洞58个，单个钻孔普遍揭露到2～5层溶洞，最多的一个钻孔揭露到8个溶洞，多层溶洞呈串珠状相连通，单孔线岩溶率10.22%～91.52%，岩洞埋深普遍在18.0～30.0 m。溶洞的顶板厚度一般0.2～2.6 m，较薄，溶洞一般都漏水、大多数为半充填状，少量全充填状，很少见未充填的溶洞。溶洞充填物一般为砾质粉质黏土，砾石的成分较杂，有石英砂岩、灰岩和铝土质岩等，砾径0.2～10 cm，土质以软塑为主，少数为可塑（图3）。
据钻探施工所揭露的58个溶洞情况显示，溶洞高度一般为0.5～9.5 m，最大高度可达15.4m（编号14），最小的高度为0.2 m，高度0.2～4 m的最多，占总量的50%， 其中0.2～2 m的占38%，2～4 m的占12%。溶洞规模分布情况见图4。
钻探施工过程中，当实施BZK5孔钻进时，附近岩溶裂隙水监测孔水位突然下降了1.68 m，这说明地下的溶洞存在连通性。当施工BZK12钻孔至36.5～37.3 m溶洞时，发生了严重的漏水导致无法施工，而相邻的XM6号CT孔口冒水（图5），同样说明了两孔之间溶洞是连通的。在施工BZK13钻孔至28.0～29.2 m段溶洞时，孔内严重漏水，难以施钻，而相邻的XMZK4号CT孔口冒水（图5），同样说明两孔之间的溶洞是连通的。在施工BZK14钻孔至深度23.8～24.60 m溶洞段时，孔内发生了严重漏水，而相邻的XMZK7号CT孔口先冒清水，后冒泥浆水（图5），说明两孔之间溶洞的连通或两孔打到的溶洞为同一溶洞。此外，在BZK2−1孔钻至深度31.40 m时，可见有人工混凝土块与灰岩胶结紧密，说明该孔溶洞与外围岩溶系统相连通（图5）。
3.3 岩溶垂向发育特征
根据钻孔揭露溶洞在垂向上的分布统计，溶洞洞顶埋深基本分布在15.50～35.8 m，占到溶洞总数的91.4%，洞顶埋深在35.0～45.0 m的有5个，占8.6%，洞顶埋深小于15 m和大于45 m的溶洞个数为0，溶洞在垂向上的分布特征有单层、双层和多层发育，多层发育呈串珠状、蜂窝状分布，溶洞之间的距离一般为0.2～8.7 m不等，首层溶洞顶板厚度（距基岩面）一般为0.40～3.30 m，占总数的81%，少数首层溶洞顶板厚度大于5 m，占比为14%。岩面以下即基岩埋深15 m以内的岩溶最发育，其线岩溶率达26.43%～51.43%（图6），多层溶洞发育的特征也反映了广花断陷盆地地壳的间歇性沉降。
3.4 物探剖面上岩溶异常特征
从研究区跨孔弹性波CT反演波速影像（图7）可看出，岩溶极为发育，溶洞规模大小不一，空间分布形态极不规则。据物探专项调查，跨孔弹性波CT法完成了29对剖面，按实际位置粘贴成17条物探测线，发现了70处低速异常，通过测线相对位置推断了43个溶洞，2个土洞，溶洞顶板高程在−6.86～−37.44 m，洞底高程范围在−10.14～−36.61 m，其溶洞规模大小不等，溶洞宽度一般都在10 m以内，占总数的80%，溶洞高度一般小于10 m，占比达76.3%，最大溶洞达31 m（宽）×22.1 m（高），溶洞为半填充或全填充砂及黏性土，半充填溶洞占总数的56%，在重点勘查区推断的土洞和溶洞中，编号T002的土洞最大（洞高和水平半径都在10 m以上）；编号R001、R005和R009的溶洞规模较大，其中R009溶洞最大（垂向和横向都延伸了20 m以上）。
![]()
图 7 研究区跨孔弹性波CT反演波速影像Figure 7. Cross hole elastic wave CT inversion wave velocity image in the study area跨孔弹性波CT法推测的溶洞和土洞基本接近实际地层情况，部分溶洞相互连通，部分溶洞应该与土层相连，形成陡坎或“鹰嘴”。该方法揭示的岩溶发育情况更为真实，是本区岩溶探测较好的方法之一。
4. 致灾因素分析
4.1 规模巨大土洞是产生岩溶地面塌陷的物质基础
根据钻探资料，位于塌陷区中部地段的基岩受F10和F12两条断裂构造的影响曾经受到了严重的风化剥蚀，形成了局部低洼地带，该处后期沉积了较厚的土层，一般为20.6～31 m，最厚处达到了50.2 m。该地段岩溶发育又受两条断裂构造的叠加交汇影响，局部地下水活跃，容易形成土洞，BZK4及XMZK3两个钻孔揭露的土洞都属于该地段的范围。由于该处土层的厚度较大，一般规模的土洞在上覆土体的自重荷载作用下不易发生冲剪失稳破坏而产生塌陷。但是当土洞逐渐发育并向上发展扩大、上部土层被掏空而变薄，再加之后期其他动力的作用，才可能引起盖层发生冲剪破坏而失稳。2008年12月19日发生的塌陷，塌陷产生的坑深度大、平面直径大，足以说明原土洞空间的规模较大。
4.2 大量修筑房屋加大土体盖层的荷载
塌陷区和其他农村地区一样，以低矮疏散的建筑布局为主，但随着经济社会的发展，村民新建的房屋越来越高，逐渐由二三层加高至六七层，甚至八九层，房屋密集使地面容积率逐渐加大，增加了地面的外荷载，使盖层土体的支撑逐渐趋于临界状态，存在着一触即发的严重隐患。
4.3 地下水的抽取改变了土洞边界的应力状态
研究区外围存在排污管道的检查施工及石井河泵站施工开挖建设，这些施工建设都曾对第四系松散岩类孔隙水和岩溶水进行过不同程度的抽排水。虽然第四系水抽取量不大，但随着地下水位的波动，已成形土洞内的充水状况将被改变，尤其是土洞内地下水位突然降低，会产生瞬时真空状态，从而改变土洞周边的应力状态，导致局部应力集中而失稳掉块，最终使土洞逐渐向上发展扩大。此外，研究区内村民在挖桩、开挖竖井和开挖建基面等过程中抽排水，造成了区内地下水位的波动变化，从而加速了水土临界状态的失稳。
4.4 区内楼房基础施工扰动是地面塌陷的诱发因素
经调查了解，研究区近年来新建的房屋多采用钻孔桩，频繁的桩基础施工极大扰动了区内地下水和土体，区内的岩、土、水和气等平衡状态受到严重破坏。本次发生的大规模岩溶塌陷就是由于钻探施工从砂性土层钻穿下部的黏性土层，黏性土层形成的土洞盖板被钻穿，加上钻孔的流水循环，使砂性土层很快变形破坏，导致土洞迅速扩大，在上部土体和外荷载的作用下，失稳致塌，即5巷18号钻孔桩施工破坏了区内的地下水、气平衡状态，诱发了本次严重的岩溶塌陷发生。
5. 地质灾害风险初步定性评估
根据国土资源部《地质灾害危险性评估规范》（DZ/T 0286-2015）以及广东省地质灾害防治协会《广东省地质灾害危险性评估实施细则》（2019年修订版）的相关规定，结合本区域岩溶地面塌陷地质灾害的形成和发育特征，对建设用地——夏茅村集中居住区的岩溶地面塌陷地质灾害风险进行初步的定性评估。
根据地质灾害风险评估一般原则，采用如下概念模型进行定性评估：
风险度 = 易发度 × 易损度
对于本地区来说，岩溶地面塌陷地质灾害易发度主要考虑7项地质环境条件指标，即：岩溶发育程度、岩溶地下水位及动态变化、覆盖土层岩性与结构、覆盖土层厚度、地面地貌特征、地质构造、人类工程活动（钻探、抽排水、基坑开挖、采矿等），其分级及具体内容见表2。地质灾害社会经济易损性指标主要分5项，即：土地利用类型、潜在经济损失、受威胁人数、房屋地基基础、房屋结构，其分级及具体内容见表3。易发度与易损度叠加即为风险度。本次评估暂不考虑各个指标的分级数值、权重以及具体计算公式，只考虑各个指标的级别进行地质灾害风险的初步定性评估。
表 2 岩溶地面塌陷易发度指标Table 2. Karst ground collapse susceptibility index稳定性
分级
评价要素不稳定
（高易发）较不稳定
（中易发）基本稳定
（低易发）稳定
（不易发）岩溶发育程度 强发育 中等发育 弱发育 不发育 岩溶地下水位及动态变化 水位埋深<5 m，
水位变化大水位埋深5～10 m，
水位变化较大水位埋深10～15 m，
水位变化较小水位埋深≥15 m ,
水位变化小覆盖土层岩性、结构 均一松散砂层或软土 均一稍密-中密砂土；双层或
多层砂土，底为砂砾双层或多层，粘性土
与砂土互层均一可塑-硬塑粘性土 覆盖土层厚度/m <10 10～20 20～30 ≥30 地面地貌特征 低洼地带临近
地表水体平原、谷地、低阶地 山前缓坡、中高阶地 台地、坡地、高阶地 地质构造 两组及以上交叉断裂 两组平行断裂 一组断裂 无断裂 人类工程活动（钻探、抽排水、
基坑开挖、采矿等）100 m以内有活动 200 m以外有活动 300 m以外有活动 无扰动地下水活动 表 3 岩溶地面塌陷易损度指标Table 3. Vulnerability index of karst ground collaps损失分级
评价要素损失大
（高易损）损失中等
（中易损）损失小
（低易损）无损失
（不易损）土地利用类型 建成区 耕地（广义） 林草地 空地、裸地 潜在经济损失/万元 >500 100～500 <小于100 无 受威胁人数/人 >30 3～30 <3以下 无 房屋地基基础 天然或浅基础 入土桩基础 入岩桩基础 无 房屋结构 砖混结构 预制板结构 框架结构 无 根据前述的地质环境条件以及地质灾害发育的具体资料来看，夏茅村岩溶发育强度为强发育，本次钻探施工的27个钻孔有19个揭露到溶洞，见洞率高达86%，单个钻孔普遍揭露到2～5层溶洞，最多的一个钻孔揭露到8个溶洞，多层溶洞呈串珠状相连通，单孔线岩溶率可达10.22%～91.52%；此外，土洞也十分发育。覆盖层岩性为双层或多层砂质土与淤泥互层。覆盖土层厚度15～20 m。地下水埋深小于1 m，在基岩面附近波动。地貌为岩溶盆地中河流低阶地，临近小河流主河道。塌陷区发现近东西向和北东向两组断裂破碎带，岩溶地面塌陷触发点正位于两条断裂带交汇处。人类活动强烈，全部塌陷坑及威胁损害房屋均位于钻探工程造成的塌陷触发点中心四周100 m以内。这些基本数据多数位于稳定性分级的“不稳定”级别，少数属于“较不稳定”级别。因此，将夏茅村岩溶地面塌陷地质灾害易发度定性为不稳定——“高易发度”。
受本次岩溶地面塌陷影响损伤的村民自建房有81栋，如果以每间房一户3人计算，则潜在受威胁人数约243人。如果以每间房屋基本价钱平均为15万元计算，则潜在经济损失约1230万元。对照表3，这两项对应“损失大”的级别。同时，夏茅村多数村民建房的地基基础为天然浅基础，少数为桩基础；房屋结构多为预制件或者砖混结构，少数为框架结构。因此，总体来说，本地的地质灾害易损性定性为“高易损性”。
根据上述地质灾害风险度概念模型，将易发度（高）与易损度（高）相叠加，则可以定性地评估夏茅村岩溶地面塌陷的潜在地质灾害风险度为大。
6. 结论与建议
6.1 结论
（1）研究区地质环境条件极其复杂，隐伏基岩为石炭系壶天群灰岩，岩溶强烈发育，基岩面起伏剧烈，高程变化较大，溶洞普遍为多层且联通性良好，上覆第四系覆盖土层厚度大，底部粉质黏土层土洞发育，形成产生岩溶地面塌陷的基础条件。
（2）水动力条件和外部荷载的变化，加剧了区内岩、土、水和气等平衡状态的打破，尤其是村民建房钻探施工从砂性土层钻穿下部的黏性土层，黏性土层形成的土洞盖板被钻穿，加上钻孔的流水循环，使砂性土层很快变形破坏，导致土洞迅速扩大，在上部土体和外荷载的作用下，发生了此次严重的岩溶地面塌陷地质灾害。
（3）运用7项地质环境条件指标确定了研究区岩溶地面塌陷的易发度高，运用5项经济人口指标定性确定其易损度也为高，由此，将二者叠加则确定本区岩溶地面塌陷地质灾害风险大。
6.2 建议
（1）建议当地后期建筑物的布局在设计和施工实施时统一对溶、土洞进行处理，损毁房屋复建时尽量考虑将单栋结构体连成片，从而增大房屋结构的整体性，以防万一出现塌陷，对上部房屋造成大的损失。
（2）塌陷区的地质环境条件复杂且脆弱，溶洞和土洞发育，当地必须严格管理控制地下水的抽排，限制地面荷载的增加，规范私人建房的技术处理，最大限度降低引发岩溶地面塌陷地质灾害发生的外部风险。
（3）塌陷区所处区域人口及建筑密度大，人类工程活动极易破坏地质环境，引发地质灾害，建议下一步加紧开展周边区域地质灾害风险评估，为精准开展地质灾害防控提供依据。
-
![]()
图 1 “阶跃型”滑坡位移曲线变形阶段划分
Figure 1. Classification of deformation stages for “step-like” landslide displacement curve
![]()
图 2 麻池村1号滑坡全貌图（镜向235°）
Figure 2. Full view of landslide No.1 in Machi Village (lens direction at 235°)
![]()
图 3 麻池村1号滑坡工程地质平面图
Figure 3. Geological plan of the landslide project No.1 in Machi Village
![]()
图 4 麻池村1号滑坡1-1’剖面布置图
Figure 4. Landslide section layout profile 1-1’ of landslide No.1 in Machi Village
![]()
图 5 麻池村1号滑坡累计位移-降雨量-时间关系图
Figure 5. Cumulative displacement-rainfall-time relationship of landslide No.1 in Machi Village
![]()
图 6 麻池村1号滑坡7次阶跃变形区间图
Figure 6. Interval graph of 7-step deformation of landslide No. 1 in Machi Village
![]()
图 7 有效降雨—改进切线角—位移速率判据曲线变化特征图
Figure 7. Characteristics variation of the criterion curve for effective rainfall-improved tangent angle-displacement rate
表 1 G6监测点破坏拐点与稳定拐点详细信息表
Table 1 Detailed information table on the damage inflection point and stable inflection point of monitoring point G6
变形区间 破坏拐点 稳定拐点 变形期次 区间天数/d 起始位移/mm 结束位移/mm 位移增量/mm 位移速率/（mm·d−1） T1 2020-07-20 2020-07-30 [416,426] 11 31.28 42.26 10.98 1.00 T2 2020-08-21 2020-08-30 [448,456] 10 53.11 70.30 17.19 1.72 T3 2020-10-03 2020-10-11 [490,498] 9 79.20 92.55 13.36 1.48 T4 2021-06-25 2021-06-30 [756,761] 6 124.59 133.23 8.64 1.44 T5 2021-08-13 2021-08-20 [805,812] 8 172.65 211.82 39.17 4.90 T6 2021-08-26 2021-09-06 [818,829] 12 213.47 538.06 324.59 27.05 T7 2021-09-19 2021-09-25 [842,848] 7 561.99 889.01 327.02 46.72 
下载: 导出CSV
表 2 一次独立降雨位移回落时间
Table 2 Summary for single independent rainfall-induced displacement pullback time
降雨日期 降雨量/mm 回落时间/d 2021-08-13 45.8 9 2021-09-28 18.1 10 2021-10-15 26.7 11
下载: 导出CSV
表 3 G6监测点“前期降雨+当期降雨”模式详细信息表
Table 3 Detailed information table on the “pre-rainfall + current rainfall” pattern details table for monitoring point G6
变形期次 降雨日期 当期降雨/mm 前期降雨/mm T1=[416,426] 2020-07-21 29.2 23.41 2020-07-22 35.1 41.44 T2=[448,456] 2020-08-21 19.8 37.42 T3=[490,498] 2020-10-03 22.9 21.92 T4=[756,761] 2021-06-25 9.8 7.80 2021-06-26 45.0 11.15 2021-06-27 10.1 44.59 T5=[805,812] 2021-08-13 45.8 13.50 2021-08-20 14.4 11.82 T6=[818,829] 2021-08-26 34.7 16.46 2021-08-28 17.4 32.67 2021-08-29 78.7 37.55 T7=[842,848] 2021-09-19 49.5 5.82
下载: 导出CSV
表 4 改进切线角—变形速率（G6）阈值及变形过程综合预警判据统计表
Table 4 Statistical table of the threshold for improved tangent angle-deformation rate (G6) and comprehensive early warning criteria for deformation process
变形阶段 改进切线
角阈值/（°）位移速率
/（mm·d−1）预警级别 初始变形阶段 α<45 v<12 等速变形阶段 α≈45 v≈12 注意级别（Ⅰ） 加速
变形
阶段初加速阶段 45<α<80 12<v<68.06 警示级别（Ⅱ） 中加速阶段 80≤α<85 68.06≤v<137.16 警戒级别（Ⅲ） 加加速阶段
（临滑）85≤α≤90 v≥137.16 警报级别（Ⅳ）
下载: 导出CSV
-
[1] 卫童瑶,殷跃平,高杨,等. 三峡库区巫山县塔坪H1滑坡变形机制[J]. 水文地质工程地质,2020,47(4):73 − 81. [WEI Tongyao,YIN Yueping,GAO Yang,et al. Deformation mechanism of the Taping H1 landslide in Wushan County in the Three Gorges Reservoir area[J]. Hydrogeology & Engineering Geology,2020,47(4):73 − 81. (in Chinese with English abstract) WEI Tongyao, YIN Yueping, GAO Yang, et al. Deformation mechanism of the Taping H1 landslide in Wushan County in the Three Gorges Reservoir area[J]. Hydrogeology & Engineering Geology, 2020, 47（4）: 73-81. （in Chinese with English abstract）
[2] 胡爱国,周伟. 地震与强降雨作用下堆积体滑坡变形破坏机理及防治方案分析[J]. 中国地质灾害与防治学报,2022,33(1):27 − 34. [HU Aiguo,ZHOU Wei. Deformation and failure mechanism and analysis on prevention measures of colluction landslide under earthquake and heavy rainfall[J]. The Chinese Journal of Geological Hazard and Control,2022,33(1):27 − 34. (in Chinese with English abstract) HU Aiguo, ZHOU Wei. Deformation and failure mechanism and analysis on prevention measures of colluction landslide under earthquake and heavy rainfall[J]. The Chinese Journal of Geological Hazard and Control, 2022, 33（1）: 27-34. （in Chinese with English abstract）
[3] 黄晓虎,雷德鑫,夏俊宝,等. 降雨诱发滑坡阶跃型变形的预测分析及应用[J]. 岩土力学,2019,40(9):3585 − 3592. [HUANG Xiaohu,LEI Dexin,XIA Junbao,et al. Forecast analysis and application of stepwise deformation of landslide induced by rainfall[J]. Rock and Soil Mechanics,2019,40(9):3585 − 3592. (in Chinese with English abstract) HUANG Xiaohu, LEI Dexin, XIA Junbao, et al. Forecast analysis and application of stepwise deformation of landslide induced by rainfall[J]. Rock and Soil Mechanics, 2019, 40（9）: 3585-3592. （in Chinese with English abstract）
[4] 陈春利,方志伟. 福建省地质灾害气象预警有效降雨模型研究[J]. 地质力学学报,2023,29(1):99 − 110. [CHEN Chunli,FANG Zhiwei. Research on an effective rainfall model for geological disaster early warning in Fujian Province, China[J]. Journal of Geomechanics,2023,29(1):99 − 110. (in Chinese with English abstract) [CHEN Chunli, FANG Zhiwei. Research on an effective rainfall model for geological disaster early warning in Fujian Province, China[J]. Journal of Geomechanics, 2023, 29(1): 99-110.(in Chinese with English abstract)
[5] 胡华,吴轩,张越. 基于模拟试验的强降雨条件下花岗岩残积土斜坡滑塌破坏机理分析[J]. 中国地质灾害与防治学报,2021,32(5):92 − 97. [HU Hua,WU Xuan,ZHANG Yue. Experimental study on slope collapse characteristics of granite residual soil slope under heavy rainfall[J]. The Chinese Journal of Geological Hazard and Control,2021,32(5):92 − 97. (in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.2021.05-11 HU Hua, WU Xuan, ZHANG Yue. Experimental study on slope collapse characteristics of granite residual soil slope under heavy rainfall[J]. The Chinese Journal of Geological Hazard and Control, 2021, 32（5）: 92-97. （in Chinese with English abstract） DOI: 10.16031/j.cnki.issn.1003-8035.2021.05-11
[6] 邓夕胜,张元,唐煜. 基于失效概率的边坡降雨阈值曲面探讨[J]. 中国地质灾害与防治学报,2021,32(3):70 − 75. [DENG Xisheng,ZHANG Yuan,TANG Yu. Investigation on slope rainfall threshold surface based on failure probablolity[J]. The Chinese Journal of Geological Hazard and Control,2021,32(3):70 − 75. (in Chinese with English abstract) DENG Xisheng, ZHANG Yuan, TANG Yu. Investigation on slope rainfall threshold surface based on failure probablolity[J]. The Chinese Journal of Geological Hazard and Control, 2021, 32（3）: 70-75. （in Chinese with English abstract）
[7] 方然可,刘艳辉,苏永超,等. 基于逻辑回归的四川青川县区域滑坡灾害预警模型[J]. 水文地质工程地质,2021,48(1):181 − 187. [FANG Ranke,LIU Yanhui,SU Yongchao,et al. A early warning model of regional landslide in Qingchuan County,Sichuan Province based on logistic regression[J]. Hydrogeology & Engineering Geology,2021,48(1):181 − 187. (in Chinese with English abstract) FANG Ranke, LIU Yanhui, SU Yongchao, et al. A early warning model of regional landslide in Qingchuan County, Sichuan Province based on logistic regression[J]. Hydrogeology & Engineering Geology, 2021, 48（1）: 181-187. （in Chinese with English abstract）
[8] 李秀珍,许强,黄润秋,等. 滑坡预报判据研究[J]. 中国地质灾害与防治学报,2003,14(4):5 − 11. [LI Xiuzhen,XU Qiang,HUANG Runqiu,et al. Research on prediction criterion for temporary prediction of landslide[J]. The Chinese Journal of Geological Hazard and Control,2003,14(4):5 − 11. (in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.2003.04.002 LI Xiuzhen, XU Qiang, HUANG Runqiu, et al. Research on prediction criterion for temporary prediction of landslide[J]. The Chinese Journal of Geological Hazard and Control, 2003, 14（4）: 5-11. （in Chinese with English abstract） DOI: 10.16031/j.cnki.issn.1003-8035.2003.04.002
[9] 陈志强,王亮清,刘顺昌. 基于斜率阈值法的滑坡预报分析[J]. 人民长江,2010,41(10):38 − 41. [CHEN Zhiqiang,WANG Liangqing,LIU Shunchang. Research on landslide prediction based on slope threshold of rainfall[J]. Yangtze River,2010,41(10):38 − 41. (in Chinese with English abstract) DOI: 10.3969/j.issn.1001-4179.2010.10.011 CHEN Zhiqiang, WANG Liangqing, LIU Shunchang. Research on landslide prediction based on slope threshold of rainfall[J]. Yangtze River, 2010, 41（10）: 38-41. （in Chinese with English abstract） DOI: 10.3969/j.issn.1001-4179.2010.10.011
[10] 王涛,刘甲美,栗泽桐,等. 中国地震滑坡危险性评估及其对国土空间规划的影响研究[J]. 中国地质,2021,48(1):21 − 39. [WANG Tao,LIU Jiamei,LI Zetong,et al. Seismic landslide hazard assessment of China and its impact on national territory spatial planning[J]. Geology in China,2021,48(1):21 − 39. (in Chinese with English abstract) [WANG Tao, LIU Jiamei, LI Zetong, et al. Seismic landslide hazard assessment of China and its impact on national territory spatial planning[J]. Geology in China, 2021, 48(1): 21-39.(in Chinese with English abstract)
[11] 许强,曾裕平,钱江澎,等. 一种改进的切线角及对应的滑坡预警判据[J]. 地质通报,2009,28(4):501 − 505. [XU Qiang,ZENG Yuping,QIAN Jiangpeng,et al. Study on a improved tangential angle and the corresponding landslide pre-warning criteria[J]. Geological Bulletin of China,2009,28(4):501 − 505. (in Chinese with English abstract) DOI: 10.3969/j.issn.1671-2552.2009.04.011 XU Qiang, ZENG Yuping, QIAN Jiangpeng, et al. Study on a improved tangential angle and the corresponding landslide pre-warning criteria[J]. Geological Bulletin of China, 2009, 28（4）: 501-505. （in Chinese with English abstract） DOI: 10.3969/j.issn.1671-2552.2009.04.011
[12] 许强,汤明高,徐开祥,等. 滑坡时空演化规律及预警预报研究[J]. 岩石力学与工程学报,2008,27(6):1104 − 1112. [XU Qiang,TANG Minggao,XU Kaixiang,et al. Research on space-time evolution laws and early warning-prediction of landslides[J]. Chinese Journal of Rock Mechanics and Engineering,2008,27(6):1104 − 1112. (in Chinese with English abstract) XU Qiang, TANG Minggao, XU Kaixiang, et al. Research on space-time evolution laws and early warning-prediction of landslides[J]. Chinese Journal of Rock Mechanics and Engineering, 2008, 27（6）: 1104-1112. （in Chinese with English abstract）
-
期刊类型引用(1)
1. 杨强国, 罗刚, 吴茂林, 张倩. 宜昌至郑万高铁联络线施工阶段隧道进出口落石风险再评估. 中国地质灾害与防治学报. 2025(04) 
本站查看
其他类型引用(0)
计量
- 文章访问数: 352
- HTML全文浏览量: 814
- PDF下载量: 177
- 被引次数: 1
邮件订阅
RSS