1.岩溶塌陷发育特征和成灾机制
前人研究者将我国岩溶塌陷形成机制归纳为重力致塌模式、潜蚀致塌模式、真空吸蚀致塌模式、冲爆致塌模式、振动致塌模式、荷载致塌模式、溶蚀致塌模式和根蚀致塌模式等8种模式(徐卫国等,1981;康彦仁,1992;李前银,2009),形成了我国岩溶塌陷研究的基本理论。国外对岩溶塌陷的成因机理研究较为关注的是水的潜蚀和侵蚀作用,影响因素包括暴雨、抽水、地表径流、下水管道渗漏等(杨立中,1997;Benson et al.,2016;Santolo et al.,2018;Gkkaya et al.,2021)。后人在此基础上,结合各自研究区特点,借助数学、力学和数值模拟等手段在机理定量化表达和特殊致塌因素分析等方面进行了深入研究,特别是极限平衡拱理论和盖层薄顶板稳定性评价的引入使得土洞型岩溶塌陷临界条件的量化研究有了明显提升(程星等,2002a,2002b;王滨等,2011;贺可强等,2017;李新华等,2021;Jia et al.;2020;Jia et al.;2021)。
岩溶塌陷影响因素包括暴雨、抽水、地表径流、下水管道渗漏等(杨立中,1997;Benson et al.,2016;Santolo et al.,2018;Gkkaya et al.,2021),其他因素还包括钻探(Tihansky,1999)、爆破和地震(Giampaolo et al.,2016;Kuniansky et al.,2016;Kuniansky et al.,2017)等。
从大量岩溶塌陷的实例可发现,地下水的动态变化是引发塌陷或加速塌陷变形的主要原因。地下水位的变动具体包括两种影响:1)地下水位下降,人类活动区开采地下水、矿井抽水、建筑基坑降水等是引起地下水位下降的主要原因;由于水浮托力的降低和覆盖层内瞬时孔隙水压力梯度提高,增加覆盖层顶板压力的同时,也使得顶板松散物质迁移性増大,大部分的岩溶地面塌陷形成于水位下降环境(雷明堂等,1993;程星等,2002),武汉、广州和唐山等城市抽排地下水,引发了大量岩溶塌陷的案例(刘道涵,2022;贺可强等,2017)。刘之葵等(2004)统计桂林市岩溶区人为塌陷中90%以上为人工抽水所致,水力梯度是决定潜蚀作用的关键,而水位变幅速率是覆盖层崩解速率的关键。死海地区,地下水位下降引起的海水侵入,加速岩溶溶蚀,沿岸以色列和约旦等国的岩溶塌陷明显増加(Yechieli et al.,2006)。2)地下水位抬升;如水利水电工程库区,地下水位抬升普遍存在的,库水可能沿岩溶管道反灌,亦会加大岩溶溶蚀速率,使裂隙加宽(Gutiérrez,2010);降雨入渗会形成垂向渗透力导致盖层土体破坏,但引起的地下水位抬升是否会触发岩溶塌陷的讨论仍存在争议,传统的饱和土渗流分析无法准确描述降雨过程中塌陷体孔隙水压力的动态变化,目前缺乏更直观的观测手段予以支持(贺可强等,2002;冯永,2009;刘长礼等,2013,2021)。
2.岩溶塌陷研究方法
目前,相关学者对覆盖型岩溶进行了物理模型试验。唐万春(2011)研究粉砂土、上层粉质黏土+下层卵石层、上层卵石层+下层粉质砂土3种覆盖层构成,并建立地下水下降与土洞顶部垂直方向位移的关系;张少波等(2019)以福建省龙岩市永定区某岩溶塌陷为研究对象,主要围绕水位波动对覆盖层孔隙水压力和盖层沉降变化的影响,展开室内物理模型试验。
岩溶塌陷具有突发性、隐蔽性、复杂性等特点,往往给工程带来巨大危害。因此,对岩溶塌陷监测预警、对其发育特征和变化规律的掌握具有重要意义。在地表,一般运用激光雷达、合成孔径雷达、热红外(TI)、全球定位系统(GPS)和地面地质雷达等航空遥感及物探方法进行定期扫描监测(Benson et al.,2011;Gutiérrez et al.,2011;Jiang et al.,2016);在地下,监测预警方法主要是对上覆盖层土体的物理力学性质进行监测,包括应力应变监测、土壤含水率监测等,地下水动力条件监测也是常用的地下预警方法(Meng et al.,2014),具有易推广、成本低等优点,但其难点在于预警判据的获取,可以借助地下水数值模拟、力学分析模拟或通过室内物理模型模拟方法确定预警阈值,GMS、FLAC 3D等模拟软件在临界条件研究方面应用较广泛(Van Dyk et al.,2016;蒙彦等,2018)。
3.粒子图像测速技术和试验系统
粒子图像测速技术于1984年首次报道。该技术起源于流体力学,起初是为了更好地理解牛顿流体中的湍流行为,后来引入岩土工程中用来分析颗粒材料变形和位移等相关问题,现如今该技术已经在能源燃烧、流场测量、实时监控、矿业加工、海洋学、地质学等各个领域广泛应用。White等(2003)首次在物理模拟实验中使用粒子图像测速技术,该技术运用在地质灾害物理模拟实验中比传统照相机拍摄图片进行定性或半定量的处理方法在空间和时间分辨率上提升了一个数量级。与粒子图像测速技术前身LDV(激光多普勒测速仪)只能测量单点流速相比,粒子图像测速技术可分析二维(甚至三维)速度矢量场,更有利于研究应变分布和变化,滑坡、泥石流和地震等地质灾害发育、演化。
近年来,越来越多的学者将粒子图像测速技术与构造物理模拟进行结合,与传统分析方法相比,该方法是一种快速且精确的定量分析技术,并且这种技术是一种非入侵性的,不会对被测量物体造成破坏。如在海底水下滑坡物理模拟实验中,该技术能将启滑全过程更好的可视化。
将粒子图像测速技术和图像捕捉与处理分析联合对实验过程中变形前后的图像进行变形及渗流分析(杜建明等,2021)。实验所需装置和材料主要包括光学平台、计算机、激光器、高速CCD相机、滤镜、模型槽、透明颗粒、孔隙流体、荧光染色剂等,详见图2。当物体被激光照射时,由于激光的相干性会在其表面呈现颗粒状结构或光暗斑点,这种结构称为“激光散斑”。由于物体在位移、振动和变形前后的散斑图样不同,因此可以利用激光散斑技术实现土体变形的观测。在物理模型试验中使用激光光源时,避免自然光线对试验结果的影响,应在黑暗或者深色背景下进行(严少华等,2020)。激光散斑图案的变化能够反映土体位移、变形情况。因而有时也可不使用荧光染色剂而直接使用激光照射透明土体模型,也可获取土体内部变化图像,研究土体变形相关问题时,多使用此种方法。相机拍摄方向与激光射人方向垂直,将试验所获图像输入计算机,借助图像分析技术,相关函数等方法构建土体内部三维图像,进而获得土体变形、渗流等信息。
图2 基于粒子测速技术地质灾害物理模型试验图像采集系统
综上,目前岩溶塌陷机理研究可以从以下方面切入:开展历时的、动态的研究,岩溶塌陷的演化是基于水文过程和岩溶过程的量变到质变的发展,有效的监测手段、直接可视化的模型试验对岩溶塌陷的成灾过程重现具有重要的提升效果。
伴随桂林入围全国第二批系统化全域推进海绵城市建设示范城市,城市建设正在提速,重大工程全面布局,岩溶塌陷造成的经济损失愈加严重,严重制约国土规划和社会经济发展。而以往关于桂林市岩溶塌陷技术手段限制还有些问题研究不够深入,对岩溶塌陷灾害的成因机理和演化机制还有待深化。现阶段相关研究仍以定性分析为主,不能满足桂林城市发展的实际需求,因此,基于定量化的岩溶塌陷灾害问题研究十分必要和迫切。
本研究以桂林市覆盖型岩溶为研究对象,以阻-透型盖层结构为原型建立二元结构地质概化模型开展模型试验,布设位移计和孔隙水压力传感器,获取相关渗流和沉降数据,分析降雨和水位波动条件下覆盖型岩溶变化规律。在此基础上,使用反滤层回填法对该模型进行治理,并借助粒子测速技术和透明土试验技术对回填和后续地下水位波动影响土体稳定性进行可视化模拟,评价反滤层法治理效果,揭示桂林典型覆盖型岩溶塌陷的形成模式和成因机制,为岩溶塌陷预测与评估提供科学依据。