Citation: | Lu Haifeng, Ning Mingcheng, Zhang Yuan, Liang Xiuyu. Stochastic generation and numerical gridding of dissolved pore-fracture composite networks[J]. Bulletin of Geological Science and Technology, 2023, 42(5): 264-272. doi: 10.19509/j.cnki.dzkq.tb20220461 |
The problem of water inrush caused by karst fissure aquifers in coal seam floor has been widely concerned. The key to reduce the occurrence of such accidents is mastering the seepage characteristics of aquifer. Before that, an accurate and reasonable aquifer model should be established first.
The karst fractured aquifer has the characteristics of composite network. At present, the numerical simulation method is mainly used to resolve this kind of problem. However, most of the existing numerical models only simulate the fracture network of aquifer, ignoring the karst phenomena such as dissolved pores with stochastic distribution, which is different from the practice. Therefore, this paper focuses on the simulation of randomly distributed dissolution pores in limestone aquifers, and combines them with the existing fracture network to obtain a more accurate karst featured aquifer network model, then convert it into a numerical ones. For this reason, firstly, some drilling and outcrop data are selected, and the contour of dissolved pores is depicted by AutoCAD software to establish corresponding morphological library, which provides sufficient samples for the stochastic generation of dissolved pore. When there are enough dissolved pore samples, the authenticity of the morphological library of dissolved pores are maintained. Secondly, the multi-factor directional block cutting technology is used to establish a storage space for the subsequent random dissolved pores. The dissolved pores are randomly extracted from corresponding morphological library, after scaling and rotating these samples into the storage space according to a certain scaling ratio, and finally the stochastic generation of dissolved pores can be realized. Subsequently, superimposed it with the fracture network generation based on the Monte Carlo method to obtain a randomly generated limestone aquifer composite grid model with a certain dissolved porosity, fracture density and spacing. Finally, the network model is transformed into FLAC3D numerical gird based on AutoCAD software, and the numerical model of limestone can be established.
According to the above method, the composite network model of limestone aquifer with different solution porosity, fracture density and spacing can be generated. The density and ductility of the fractures in the fracture network are labelled as low, moderately low, medium, moderately high and high respectively. The dissolution rate of dissolved pore is classified into 0.05, 0.10, 0.15, 0.20, 0.25 and 0.35 respectively. Taking high density-low ductility fracture network as an example. Two models with dissolved pore rate of 0.05 and 0.25 are superimposed, and statistial analysis was carried out on the trace length of fracture, long axis ratio and dip angle of dissolved pore. The obtained results are in line with expectations, which proves the effectiveness of the generated model.
The research results have certain significance for the establishment of three-dimensional random dissolution pore model, and lay a foundation for the subsequent analysis of seepage characteristics of karst fracture aquifer based on numerical simulation.
[1] |
Wu Q, Guo X M, Shen J J, et al. Relevamiento del riesgo de irrupción de agua desde los acuíferos subyacentes de la mina de carbón Gushuyuan, China[J]. Mine Water and the Environment, 2017, 36(1): 96-103. doi: 10.1007/s10230-016-0410-8
|
[2] |
袁亮. 煤及共伴生资源精准开采科学问题与对策[J]. 煤炭学报, 2019, 44(1): 1-9. https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201901001.htm
Yuan L. Scientific problems and countermeasures of precise mining of coal and associated resources[J]. Chinese Journal of Coal, 2019, 44(1): 1-9(in Chinese with English abstract). https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201901001.htm
|
[3] |
鲁海峰, 孟祥帅, 张元, 等. 采场底板层状结构关键层隔水性能力学分析[J]. 中国矿业大学学报, 2020, 49(6): 1057-1066. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGKD202006004.htm
Lu H F, Meng X S, Zhang Y, et al. Mechanical analysis of waterproof performance of key layers of layered structure of stope floor[J]. Journal of China University of Mining and Technology, 2020, 49(6): 1057-1066(in Chinese with English abstract). https://www.cnki.com.cn/Article/CJFDTOTAL-ZGKD202006004.htm
|
[4] |
梁永平, 申豪勇, 高旭波. 中国北方岩溶地下水的研究进展[J]. 地质科技通报, 2022, 41(5): 199-219. doi: 10.19509/j.cnki.dzkq.2022.0199
Liang Y P, Shen H Y, Gao X B. Research progress of karst groundwater in northern China[J]. Bulletin of Geological Science and Technology, 2022, 41(5): 199-219(in Chinese with English abstract). doi: 10.19509/j.cnki.dzkq.2022.0199
|
[5] |
宋晓晨, 徐卫亚. 裂隙岩体渗流模拟的三维离散裂隙网络数值模型(Ⅰ): 裂隙网络的随机生成[J]. 岩石力学与工程学报, 2004, 23(12): 2015-2020. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200412014.htm
Song X C, Xu W Y. Three-dimensional discrete fracture network numerical model for seepage simulation of fractured rock mass (Ⅰ): Random generation of fracture network[J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(12): 2015-2020(in Chinese with English abstract). https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200412014.htm
|
[6] |
王家臣, 常来山, 陈亚军, 等. 露天矿节理岩体三维网络模拟与概率损伤分析[J]. 北京科技大学学报, 2005, 27(1): 1-4. https://www.cnki.com.cn/Article/CJFDTOTAL-BJKD200501001.htm
Wang J C, Chang L S, Chen Y J, et, al. Three-dimensional network simulation and probabilistic damage analysis of jointed rock mass in open pit mines[J]. Journal of University of Science and Technology Beijing, 2005, 27(1): 1-4(in Chinese with English abstract). https://www.cnki.com.cn/Article/CJFDTOTAL-BJKD200501001.htm
|
[7] |
石祥超, 张茹, 高明忠, 等. 顶板随机裂隙对采动煤岩体支承压力的影响[J]. 中国矿业大学学报, 2013, 42(6): 948-953. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGKD201306010.htm
Shi X C, Zhang R, Gao M Z, et al. Influence of random cracks in the roof on the bearing pressure of mining coal and rock mass[J]. Journal of China University of Mining and Technology, 2013, 42(6): 948-953(in Chinese with English abstract). https://www.cnki.com.cn/Article/CJFDTOTAL-ZGKD201306010.htm
|
[8] |
李爱华, 朱江. 基于二维裂隙网络模拟的岩块搜索与岩层追踪方法[J]. 水利水运工程学报, 2014, 36(6): 65-70. https://www.cnki.com.cn/Article/CJFDTOTAL-SLSY201406010.htm
Li A H, Zhu J. Rock block search and rock formation tracking method based on two-dimensional fracture network simulation[J]. Chinese Journal of Water Conservancy and Water Transport Engineering, 2014, 36(6): 65-70(in Chinese with English abstract). https://www.cnki.com.cn/Article/CJFDTOTAL-SLSY201406010.htm
|
[9] |
Li L H, Huang B X, Li Y Y, et al. Multi-scale modeling of shale laminas and fracture networks in the Yanchang Formation, southern Ordos Basin, China[J]. Engineering Geology, 2018, 243: 231-240.
|
[10] |
Xu W T, Zhang Y S, Li X Z, et al. Study on three-dimensional fracture network connectivity path of rock mass and seepage characteristics based on equivalent pipe network[J]. Environmental Earth Sciences, 2019, 78(16): 1-21.
|
[11] |
张学文. 断层带岩体裂隙网络蒙特卡罗模拟及应用[J]. 陕西煤炭, 2012, 31(1): 59-61. https://www.cnki.com.cn/Article/CJFDTOTAL-SXMJ201201026.htm
Zhang X W. Monte Carlo simulation and application of rock fracture network in fault zone[J]. Shaanxi Coal, 2012, 31(1): 59-61(in Chinese with English abstract). https://www.cnki.com.cn/Article/CJFDTOTAL-SXMJ201201026.htm
|
[12] |
于本昌. 基于蒙特卡罗随机裂隙的岩溶水工隧洞稳定性及渗流特性研究[D]. 成都: 西南交通大学, 2015.
Yu B C. Study on the stability and seepage characteristics of karst hydraulic tunnel based on Monte Carlo random fractures[D]. Chengdu: Southwest Jiaotong University, 2015(in Chinese with English abstract).
|
[13] |
Ni X M, Chen W X, Li Z Y, et al. Reconstruction of different scales of pore-fractures network of coal reservoir and its permeability prediction with Monte Carlo method[J]. International Journal of Mining Science and Technology, 2017, 27(4): 693-699.
|
[14] |
刘耀儒, 杨强, 覃振朝. 基于统计模型的裂隙岩体渗流场的并行数值模拟[J]. 岩石力学与工程学报, 2008, 27(4): 736-742. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200804014.htm
Liu Y R, Yang Q, Qin Z C. Parallel numerical simulation of seepage field in fractured rock mass based on statistical model[J]. Chinese Journal of Rock Mechanics and Engineering. 2008, 27(4): 736-742(in Chinese with English abstract). https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200804014.htm
|
[15] |
Han C N, Zhou D H, Liao W Y, et al. Probability statistical model and network simulation technology of random fractured rock mass[J]. Applied Mechanics and Materials, 2014, 580/583: 3307.
|
[16] |
敖雪菲, 王晓玲, 赵梦琦, 等. 坝基裂隙岩体三维灌浆数值模拟[J]. 水利学报, 2017, 48(8): 945-954. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201708007.htm
Ao X F, Wang X L, Zhao M Q, et al. Three-dimensional grouting numerical simulation of fissure rock mass in dam foundation[J]. Journal of Hydraulic Engineering, 2017, 48(8): 945-954(in Chinese with English abstract). https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201708007.htm
|
[17] |
Snow D T. A parallel plate model of fractured permeable media[D]. Berkeley: University of California, 1965.
|
[18] |
方艺翔, 李卓, 范光亚, 等. 监测资料、压水试验与综合物探法在某心墙坝渗漏识别中的应用研究[J]. 水利水电技术: 中英文, 2022, 53(2): 87-97. https://www.cnki.com.cn/Article/CJFDTOTAL-SJWJ202202009.htm
Fang Y X, Li Z, Fan G Y, et al. Application of monitoring data, pressurized water test and comprehensive geophysical method in leakage identification of a core wall dam[J]. Water Conservancy and Hydropower Technology: Chinese and English, 2022, 53(2): 87-97(in Chinese with English abstract). https://www.cnki.com.cn/Article/CJFDTOTAL-SJWJ202202009.htm
|
[19] |
胡成, 陈刚, 曹孟雄, 等. 基于离散裂隙网络法和水流数值模拟技术的地下水封洞库水封性研究[J]. 地质科技通报, 2022, 41(1): 119-126, 136. doi: 10.19509/j.cnki.dzkq.2022.0029
Hu C, Chen G, Cao M X, et al. Study on water sealing of underground water-sealed caverns based on discrete fracture network method and flow numerical simulation technology[J]. Bulletin of Geological Science and Technology, 2022, 41(1): 119-126, 136(in Chinese with English abstract). doi: 10.19509/j.cnki.dzkq.2022.0029
|
[20] |
刘新荣, 杜立兵, 邓志云, 等. 基于多因素椭圆堆叠的岩土细观建模方法[J]. 岩土力学, 2020, 41(11): 3797-3809. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202011032.htm
Liu X R, Du L B, Deng Z Y, et al. Geotechnical mesoscopic modeling method based on multi-factor elliptical stacking[J]. Rock and Soil Mechanics, 2020, 41(11): 3797-3809(in Chinese with English abstract). https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202011032.htm
|
[21] |
刘新荣, 杜立兵, 邓志云, 等. 基于闵科夫斯基差和优化波前法的二维天然堆石料生成方法及应用[J]. 岩石力学与工程学报, 2020, 39(9): 1832-1846. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202009010.htm
Liu X R, Du L B, Deng Z Y, et al. Two-dimensional natural rockfill generation method and application based on Minkowski difference and optimized wavefront method[J]. Chinese Journal of Rock Mechanics and Engineering, 2020, 39(9): 1832-1846(in Chinese with English abstract). https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202009010.htm
|
[22] |
Wang C, Liang V. A packing generation scheme for the granular assemblies with planar elliptical particles[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 1997, 21(5): 347-358.
|
[23] |
Feng Y T, Han K D, Owen R J. Filling domains with disks: An advancing front approach[J]. International Journal for Numerical Methods in Engineering, 2003, 56(5): 699-713.
|
[24] |
王涛, 韩煊, 苏凯, 等. FLAC3D数值模拟方法及工程应用: 深入剖析FLAC3D 5.0[M]. 第2版. 北京: 中国建筑工业出版社, 2019.
Wang T, Han X, Su K, et al. FLAC3D numerical simulation method and engineering application: In-depth analysis of FLAC3D 5.0[M]. Second Edition. Beijing: China Construction Industry Press, 2019(in Chinese).
|