Numerical simulation of critical flow velocity in the covered karst soil cave formation
-
摘要:
覆盖型岩溶塌陷灾害具有隐蔽性、突发性、不确定性、机理复杂性等特点,致使其预防、治理工作十分困难。覆盖型岩溶塌陷过程中,覆盖层土体形成土洞,土洞发育、扩展直至地表塌陷。现阶段土洞塌陷研究方法多为有限单元法等数值模拟,但有限单元法无法解决颗粒迁移、土体崩落等岩体大变形问题。在浙江省江山市岩溶塌陷地质条件分析的基础上,利用颗粒流数值模拟方法(PFC)对覆盖型岩溶土洞的形成过程进行了模拟,从微观角度揭示土洞形成过程中特殊土体颗粒的位移、裂纹的发展趋势和系统不平衡力的变化情况等,再现外力作用下岩溶土洞形成的多场变化特征。结果表明:(1)岩溶土洞发育过程可概化为裂纹产生、连通面形成、颗粒剥落3个阶段;(2)岩溶土洞开始发育时,土体颗粒位移曲线、裂纹的贯通程度、不平衡力三者之间具有良好的一致性,且垂直方向扩展速度大于水平速度;(3)覆盖层土体为黏土时初始土洞形成的临界水流速度为0.104 cm/s。研究结果可为岩溶塌陷治理提供参考数据。
Abstract:The covered karst collapse has the characteristics of concealment, suddenness, uncertainty, and mechanism complexity, which leads to its prevention and treatment more difficult. In the process of covered karst collapse, soil caves formed in the overlying soil, develop and expand until the surface collapse. Numerical simulation such as the finite element method is an effective method of soil cave collapse; however, the finite element method can not solve the problem of large deformation of rock mass such as particle migration and soil collapse. Based on the analysis of the geological conditions of karst collapse in Jiangshan city, Zhejiang province, this study used the particle flow code (PFC) to simulate the formation process of covered karst soil cave and then to reveal the displacement of special soil particles, the development trend of cracks, and the change of system imbalance force on a micro level. The multi-field variation characteristics of the formation of covered karst soil cave under external force conditions were reproduced. The results show that the development process of karst soil cave could be summarized into three stages including crack generation, connected surface formation, and particle peeling. As the karst soil cave began to develop, there was a good consistency among the displacement curve of soil particles, the penetration degree of cracks, and the unbalance force, and the vertical expansion speed was larger than the horizontal speed. The critical flow velocity of the initial soil cavity is 0.104 cm/s when the overburden soil is clay. This study can provide basic information for the treatment of karst collapse.
-
Key words:
- particle flow code /
- covered /
- karst collapse /
- soil caves /
- critical flow velocity
-
-
表 1 土体力学参数标定结果
Table 1. calibration results of soil mechanical parameters
参数名称 黏聚力/ kN 内摩擦角 初始弹性模量 泊松比 抗拉强度/ kPa 渗透系数/(cm·s−1) 实际值 28 24 10.4 0.22 5 5.6×10−7 模拟值 30 26 11.2 0.2 4.8 5.7×10−7 误差 7.1% 8.3% 7.7% 9.1% 4% 1.7% 表 2 土体微观参数取值表
Table 2. soil microscopic parameter table
参数 法向刚度/(N·m−1) 切向刚度/(N·m−1) 摩擦系数 法向黏结强度/Pa 切向黏结强度/Pa 传导系数/(cm·s−1) 取值 1 $ \text{×} $ 1075 $ \text{×} $ 1060.3 8 $ \text{×} $ 1042 $ \text{×} $ 1045 $ \text{×} $ 10−11 -
[1] 金晓文,陈植华,曾斌,等. 岩溶塌陷机理定量研究的初步思考[J]. 中国岩溶,2013,32(4):437 − 446. [JIN Xiaowen,CHEN Zhihua,ZENG Bin,et al. Preliminary thinking of quantitative research on the mechanism of Karst collapse[J]. Carsologica Sinica,2013,32(4):437 − 446. (in Chinese with English abstract)]
JIN Xiaowen, CHEN Zhihua, ZENG Bin, et al. Preliminary thinking of quantitative research on the mechanism of Karst collapse[J]. Carsologica Sinica, 2013, 32(4): 437 − 446. (in Chinese with English abstract)
[2] 简文彬,吴振祥. 地质灾害及其防治[M]. 北京:人民交通出版社,2015. [JIAN Wenbin,WU Zhenxiang. Geological hazard and control[M]. Beijing:China Communications Press,2015. (in Chinese)]
JIAN Wenbin, WU Zhenxiang. Geological hazard and control[M]. Beijing: China Communications Press, 2015. (in Chinese)
[3] 周长松,邹胜章,朱丹尼,等. 广昆铁路复线秀宁隧道大皮坡—中村段岩溶塌陷成因[J]. 水文地质工程地质,2019,46(3):146 − 152. [ZHOU Changsong,ZOU Shengzhang,ZHU Danni,et al. An analysis of the cause of Karst collapses near the Dapipo-Zhongcun section of the Xiuning tunnel of the Guangzhou-Kunming railway[J]. Hydrogeology & Engineering Geology,2019,46(3):146 − 152. (in Chinese with English abstract)]
ZHOU Changsong, ZOU Shengzhang, ZHU Danni, et al. An analysis of the cause of Karst collapses near the Dapipo-Zhongcun section of the Xiuning tunnel of the Guangzhou-Kunming railway[J]. Hydrogeology & Engineering Geology, 2019, 46(3): 146 − 152. (in Chinese with English abstract)
[4] 蒋小珍,冯涛,郑志文,等. 岩溶塌陷机理研究进展[J]. 中国岩溶,2023,42(3):517 − 527. [JIANG Xiaozhen,FENG Tao,ZHENG Zhiwen,et al. A review of Karst collapse mechanisms[J]. Carsologica Sinica,2023,42(3):517 − 527. (in Chinese with English abstract)] doi: 10.11932/karst20230304
JIANG Xiaozhen, FENG Tao, ZHENG Zhiwen, et al. A review of Karst collapse mechanisms[J]. Carsologica Sinica, 2023, 42(3): 517 − 527. (in Chinese with English abstract) doi: 10.11932/karst20230304
[5] 罗小杰,罗程. 岩溶地面塌陷三机理理论及其应用[J]. 中国岩溶,2021,40(2):171 − 188. [LUO Xiaojie,LUO Cheng. Three-Mechanism Theory(TMT) of Karst ground collapse and its application[J]. Carsologica Sinica,2021,40(2):171 − 188. (in Chinese with English abstract)] doi: 10.11932/karst2021y001
LUO Xiaojie, LUO Cheng. Three-Mechanism Theory(TMT) of Karst ground collapse and its application[J]. Carsologica Sinica, 2021, 40(2): 171 − 188. (in Chinese with English abstract) doi: 10.11932/karst2021y001
[6] 罗小杰. 也论覆盖型岩溶地面塌陷机理[J]. 工程地质学报,2015,23(5):886 − 895. [LUO Xiaojie. Further discussion on mechanism of covered Karst ground collapse[J]. Journal of Engineering Geology,2015,23(5):886 − 895. (in Chinese with English abstract)]
LUO Xiaojie. Further discussion on mechanism of covered Karst ground collapse[J]. Journal of Engineering Geology, 2015, 23(5): 886 − 895. (in Chinese with English abstract)
[7] 任新红,郭永春,王清海,等. 覆盖型岩溶潜蚀塌陷临界裂隙开度模型试验研究[J]. 水文地质工程地质,2012,39(5):84 − 87. [REN Xinhong,GUO Yongchun,WANG Qinghai,et al. Experiment research of the critical fissure width of collapse of covered Karst during underground erosion[J]. Hydrogeology & Engineering Geology,2012,39(5):84 − 87. (in Chinese with English abstract)]
REN Xinhong, GUO Yongchun, WANG Qinghai, et al. Experiment research of the critical fissure width of collapse of covered Karst during underground erosion[J]. Hydrogeology & Engineering Geology, 2012, 39(5): 84 − 87. (in Chinese with English abstract)
[8] 唐灵明. 强降雨作用下覆盖型岩溶塌陷模型试验研究[D]. 桂林:桂林理工大学,2021. [TANG Lingming. Experimental study on covering karst collapse model under heavy rainfall[D]. Guilin:Guilin University of Technology,2021. (in Chinese with English abstract)]
TANG Lingming. Experimental study on covering karst collapse model under heavy rainfall[D]. Guilin: Guilin University of Technology, 2021. (in Chinese with English abstract)
[9] 陈余道, 朱学愚, 蒋亚萍. 粘性土土洞形成的水化学侵蚀实验[J]. 水文地质工程地质,1997,24(1):29 − 32. [CHEN Yudao, ZHU Xueyu, JIANG Yaping. Hydrochemical erosion test of formation of clayey soil hole[J]. Hydrogeology & Engineering Geology,1997,24(1):29 − 32. (in Chinese with English abstract)]
CHEN Yudao, ZHU Xueyu, JIANG Yaping. Hydrochemical erosion test of formation of clayey soil hole[J]. Hydrogeology & Engineering Geology, 1997, 24(1): 29 − 32. (in Chinese with English abstract)
[10] 李慎奎. 武汉地区沙漏型岩溶塌陷数值分析与模型试验研究[J]. 隧道建设(中英文),2020,40(7):981 − 987. [LI Shenkui. Numerical analysis and model test of hourglass-shaped Karst collapse in Wuhan Area[J]. Tunnel Construction,2020,40(7):981 − 987. (in Chinese with English abstract)]
LI Shenkui. Numerical analysis and model test of hourglass-shaped Karst collapse in Wuhan Area[J]. Tunnel Construction, 2020, 40(7): 981 − 987. (in Chinese with English abstract)
[11] 陶小虎,叶明,龚建师,等. 基于LBM-DEM细观数值模拟的水力诱导覆盖型岩溶地面塌陷发育过程分析[J]. 中国地质灾害与防治学报,2024,35(1):124 − 131. [TAO Xiaohu,YE Ming,GONG Jianshi,et al. Analysis of the formation process of the covered Karst ground collapse induced by groundwater changes based on the coupled LBM-DEM numerical simulation at micro scale[J]. The Chinese Journal of Geological Hazard and Control,2024,35(1):124 − 131. (in Chinese with English abstract)]
TAO Xiaohu, YE Ming, GONG Jianshi, et al. Analysis of the formation process of the covered Karst ground collapse induced by groundwater changes based on the coupled LBM-DEM numerical simulation at micro scale[J]. The Chinese Journal of Geological Hazard and Control, 2024, 35(1): 124 − 131. (in Chinese with English abstract)
[12] 张鑫. 覆盖型岩溶塌陷模型试验与数值模拟研究[D]. 合肥:合肥工业大学,2017. [ZHANG Xin. Model test and numerical simulation of overlying Karst collapse[D]. Hefei:Hefei University of Technology,2017. (in Chinese with English abstract)]
ZHANG Xin. Model test and numerical simulation of overlying Karst collapse[D]. Hefei: Hefei University of Technology, 2017. (in Chinese with English abstract)
[13] 李万有, 陈立龙, 李爽. 覆盖层的抗剪强度与土洞型岩溶塌陷高度关系的数值分析研究[J]. 林业科技情报,2010,42(1):102 − 104. [LI Wanyou, CHEN Lilong, LI Shuang. Study on numerical analysis of the relationship between shear strength of cover layer and the soil cave Karst collapse height[J]. Forestry Science and Technology Information,2010,42(1):102 − 104. (in Chinese with English abstract)]
LI Wanyou, CHEN Lilong, LI Shuang. Study on numerical analysis of the relationship between shear strength of cover layer and the soil cave Karst collapse height[J]. Forestry Science and Technology Information, 2010, 42(1): 102 − 104. (in Chinese with English abstract)
[14] 程星,彭世寿. 岩溶区地下水位下降致塌的数值模拟研究[J]. 地球与环境,2005,33(增刊1):119 − 123. [CHENG Xing,PENG Shishou. The numerical simulating study on the Karst collapse of the fall of groundwater table[J]. Earth and Environment,2005,33(Sup 1):119 − 123. (in Chinese with English abstract)]
CHENG Xing, PENG Shishou. The numerical simulating study on the Karst collapse of the fall of groundwater table[J]. Earth and Environment, 2005, 33(Sup 1): 119 − 123. (in Chinese with English abstract)
[15] 王文轩,夏源. 真空吸蚀致岩溶塌陷的稳定性分析及其数值模拟[J]. 地下水,2023,45(5):18 − 20. [WANG Wenxuan,XIA Yuan. Stability analysis and numerical simulation of Karst collapse caused by vacuum suction[J]. Ground Water,2023,45(5):18 − 20. (in Chinese with English abstract)]
WANG Wenxuan, XIA Yuan. Stability analysis and numerical simulation of Karst collapse caused by vacuum suction[J]. Ground Water, 2023, 45(5): 18 − 20. (in Chinese with English abstract)
[16] 于林弘,颜嘉良,于晓静,等. 基于FLAC3D的岩溶土洞演化及数值模拟分析[J]. 地下水,2020,42(4):55 − 57. [YU Linhong,YAN Jialiang,YU Xiaojing,et al. Evolution and numerical simulation of Karst caves based on FLAC3D[J]. Ground Water,2020,42(4):55 − 57. (in Chinese with English abstract)]
YU Linhong, YAN Jialiang, YU Xiaojing, et al. Evolution and numerical simulation of Karst caves based on FLAC3D[J]. Ground Water, 2020, 42(4): 55 − 57. (in Chinese with English abstract)
[17] 孙金辉. 覆盖型岩溶塌陷临界参数模型试验与数值模拟研究[D]. 成都:西南交通大学,2011. [SUN Jinhui. Study on Critical Parameter in Cover Karst Collapse by Model Experiment and Numerical Simulation[D]. Chengdu:Southwest Jiaotong University,2011. (in Chinese with English abstract)]
SUN Jinhui. Study on Critical Parameter in Cover Karst Collapse by Model Experiment and Numerical Simulation[D]. Chengdu: Southwest Jiaotong University, 2011. (in Chinese with English abstract)
[18] 薛明明,陈学军,宋宇,等. 基于FLAC3D不同降雨速率下土洞致塌规律研究[J]. 中国岩溶,2022,41(6):905 − 914. [XUE Mingming,CHEN Xuejun,SONG Yu,et al. A study on collapse law of soil cave with different rainfall rates based on FLAC3D[J]. Carsologica Sinica,2022,41(6):905 − 914. (in Chinese with English abstract)] doi: 10.11932/karst20220605
XUE Mingming, CHEN Xuejun, SONG Yu, et al. A study on collapse law of soil cave with different rainfall rates based on FLAC3D[J]. Carsologica Sinica, 2022, 41(6): 905 − 914. (in Chinese with English abstract) doi: 10.11932/karst20220605
[19] 苏添金,洪儒宝,简文彬. 覆盖型岩溶土洞致灾过程的数值模拟与预测[J]. 自然灾害学报,2018,27(5):179 − 187. [SU Tianjin,HONG Rubao,JIAN Wenbin. Numerical simulation and prediction of covered Karst collapse[J]. Journal of Natural Disasters,2018,27(5):179 − 187. (in Chinese with English abstract)]
SU Tianjin, HONG Rubao, JIAN Wenbin. Numerical simulation and prediction of covered Karst collapse[J]. Journal of Natural Disasters, 2018, 27(5): 179 − 187. (in Chinese with English abstract)
[20] 苟颉龙,楼谦谦,杨小兰,等. 浙江省江山市地质灾害成因分析与防治对策研究[J]. 科技通报,2018,34(10):253 − 258. [GOU Jielong,LOU Qianqian,YANG Xiaolan,et al. Cause analysis and corresponding countermeasures of geological hazards in the Jiangshan city,Zhejiang Province[J]. Bulletin of Science and Technology,2018,34(10):253 − 258. (in Chinese with English abstract)]
GOU Jielong, LOU Qianqian, YANG Xiaolan, et al. Cause analysis and corresponding countermeasures of geological hazards in the Jiangshan city, Zhejiang Province[J]. Bulletin of Science and Technology, 2018, 34(10): 253 − 258. (in Chinese with English abstract)
[21] 曾斌,杨木易,邵长杰,等. 基于层次分析法的杭长高速岩溶塌陷易发性评价[J]. 安全与环境工程,2018,25(1):29 − 38. [ZENG Bin,YANG Muyi,SHAO Changjie,et al. Susceptibility assessment of Karst collapse of Hangchang expressway projects based on analytic hierarchy process[J]. Safety and Environmental Engineering,2018,25(1):29 − 38. (in Chinese with English abstract)]
ZENG Bin, YANG Muyi, SHAO Changjie, et al. Susceptibility assessment of Karst collapse of Hangchang expressway projects based on analytic hierarchy process[J]. Safety and Environmental Engineering, 2018, 25(1): 29 − 38. (in Chinese with English abstract)
[22] 黄立勇. 江山市岩溶塌陷机理及易发性模拟预测研究[D]. 武汉:中国地质大学,2015. [HUANG Liyong. Research on the Simulation Prediction and Mechanism of Karst Collapse in Jiangshan[D]. Wuhan:China University of Geosciences,2015. (in Chinese with English abstract)]
HUANG Liyong. Research on the Simulation Prediction and Mechanism of Karst Collapse in Jiangshan[D]. Wuhan: China University of Geosciences, 2015. (in Chinese with English abstract)
[23] 周健,贾敏才. 土工细观模型试验与数值模拟[M]. 北京:科学出版社,2008. [ZHOU Jian,JIA Mincai. Meso-model test and numerical simulation of geotechnical engineering[M]. Beijing:Science Press,2008. (in Chinese with English abstract)]
ZHOU Jian, JIA Mincai. Meso-model test and numerical simulation of geotechnical engineering[M]. Beijing: Science Press, 2008. (in Chinese with English abstract)
[24] Itasca Consulting Group Inc. PFC2D particle flow code in 2 dimensions-optional features[M]. Minneapolies :Itasca Consulting Group Inc,2005.
[25] Itasca Consulting Group Inc. Verification Problems and Example Applications[M] Minneapolies:Itasca Consulting Group Inc,2005.
[26] CUNDALL P A. PFC2D User’s Manual (Version2. 0)[R]. Minnesota:Itasca Consulting Group Inc,1999.
[27] 陈小茜,曾斌,王春晖,等. 基于颗粒流理论的层状非均质介质等效渗透系数研究[J]. 水文地质工程地质,2018,45(2):7 − 12. [CHEN Xiaoxi,ZENG Bin,WANG Chunhui,et al. A study of the equivalent hydraulic conductivity of layered heterogeneous media with the particle flow theory[J]. Hydrogeology & Engineering Geology,2018,45(2):7 − 12. (in Chinese with English abstract)]
CHEN Xiaoxi, ZENG Bin, WANG Chunhui, et al. A study of the equivalent hydraulic conductivity of layered heterogeneous media with the particle flow theory[J]. Hydrogeology & Engineering Geology, 2018, 45(2): 7 − 12. (in Chinese with English abstract)
[28] 雷明堂,蒋小珍,李瑜. 岩溶塌陷模型试验——以武昌为例[J]. 地质灾害与环境保护,1993,4(2):39 − 44. [LEI Mingtang,JIANG Xiaozhen,LI Yu. Model test of Karst collapse:taking Wuchang as an example[J]. Journal of Geological Hazards and Environment Preservtion,1993,4(2):39 − 44. (in Chinese)]
LEI Mingtang, JIANG Xiaozhen, LI Yu. Model test of Karst collapse: taking Wuchang as an example[J]. Journal of Geological Hazards and Environment Preservtion, 1993, 4(2): 39 − 44. (in Chinese)
[29] 万志清,秦四清,李志刚,等. 土洞形成的机理及起始条件[J]. 岩石力学与工程学报,2003,22(8):1377 − 1382. [WAN Zhiqing,QIN Siqing,LI Zhigang,et al. Formation mechanism and initial condition of soil cavity[J]. Chinese Journal of Rock Mechanics and Engineering,2003,22(8):1377 − 1382. (in Chinese with English abstract)] doi: 10.3321/j.issn:1000-6915.2003.08.028
WAN Zhiqing, QIN Siqing, LI Zhigang, et al. Formation mechanism and initial condition of soil cavity[J]. Chinese Journal of Rock Mechanics and Engineering, 2003, 22(8): 1377 − 1382. (in Chinese with English abstract) doi: 10.3321/j.issn:1000-6915.2003.08.028
-