2025 Vol. 44, No. 1
Article Contents

WANG Jinli, CHEN Xi, ZHANG Zhicai, KANG Jianrong, HU Jinshan. MODFLOW-based analysis on seepage in discrete fissure networks[J]. Carsologica Sinica, 2025, 44(1): 1-14. doi: 10.11932/karst2024y043
Citation: WANG Jinli, CHEN Xi, ZHANG Zhicai, KANG Jianrong, HU Jinshan. MODFLOW-based analysis on seepage in discrete fissure networks[J]. Carsologica Sinica, 2025, 44(1): 1-14. doi: 10.11932/karst2024y043

MODFLOW-based analysis on seepage in discrete fissure networks

  • Studies on karst bedrock fissures face challenge in modeling due to the heterogeneity and anisotropy of fissures. Groundwater primarily flows through fissure and pipe networks within karst aquifer systems, while the bedrock pores and micro-fissures mainly serve as water storage. Therefore, the ideal hydrological model for karst areas is the continuum–fissure–pipe model. Fissures, as the main pathway for groundwater flow, significantly influence seepage in fissure media due to their connectivity. The fissure–continuum model (FC), which integrates the advantages of continuum models and discrete fissure network (DFN) models, can realistically characterize the dual texture of karst aquifers. The FC model can also account for matrix diffusion and water exchange between the matrix and fissures. However, this model remains challenging to accurately determine water exchange between rock matrix and fissure media.

    This study was based on the statistical distribution of fissure geometry and hydraulic parameters. It utilized the Monte Carlo stochastic simulation technique and MATLAB program to generate a two-dimensional fissure network diagram that mirrored the distribution of actual rock fissures. The adjacency matrix of an undirected graph of graph theory was employed to represent the intersection relationships between fissures in the fissure network. In addition, this study utilized percolation theory to eliminate isolated fissures and clusters of isolated fissures that cannot conduct fluid. This approach facilitated the generation of a first-level connectivity diagram. Fissures with single nodes were then eliminated to create a second-level connectivity diagram, and dead-ends in fissures were eliminated to generate a third-level connectivity diagram. Based on the connectivity diagrams of fissure networks, the cubic law and the continuity equation for seepage were applied to establish a single-phase, saturated, and stable DFN model. Each fissure in the connectivity diagram was mapped onto a finite difference grid to develop an FC model that combines DFN with MODFLOW. When eliminating the fissure dead-ends, computational rounding may result in the loss of some nodes. Therefore, the second-level and third-level connectivity diagrams were separately used for DFN seepage simulation and FC seepage simulation, respectively. Simulation programs were developed to reproduce cross-fissure test scenarios, and two multi-fissure cases were analyzed to investigate the effects of finite difference grid resolution and fissure dead-ends on the accuracy of stable seepage simulations. The DFN was used as a metric to evaluate the correctness and validity of FC model code writing. This study was based on the statistical distribution of fissure geometry and hydraulic parameters. It utilized the Monte Carlo stochastic simulation technique and MATLAB program to generate a two-dimensional fissure network diagram that mirrored the distribution of actual rock fissures. The adjacency matrix of an undirected graph of graph theory was employed to represent the intersection relationships between fissures in the fissure network. In addition, this study utilized percolation theory to eliminate isolated fissures and clusters of isolated fissures that cannot conduct fluid. This approach facilitated the generation of a first-level connectivity diagram. Fissures with single nodes were then eliminated to create a second-level connectivity diagram, and dead-ends in fissures were eliminated to generate a third-level connectivity diagram. Based on the connectivity diagrams of fissure networks, the cubic law and the continuity equation for seepage were applied to establish a single-phase, saturated, and stable DFN model. Each fissure in the connectivity diagram was mapped onto a finite difference grid to develop an FC model that combines DFN with MODFLOW. When eliminating the fissure dead-ends, computational rounding may result in the loss of some nodes. Therefore, the second-level and third-level connectivity diagrams were separately used for DFN seepage simulation and FC seepage simulation, respectively. Simulation programs were developed to reproduce cross-fissure test scenarios, and two multi-fissure cases were analyzed to investigate the effects of finite difference grid resolution and fissure dead-ends on the accuracy of stable seepage simulations. The DFN was used as a metric to evaluate the correctness and validity of FC model code writing. Results show that the total flow simulation values and head simulation of the fine grid under the two test conditions are in good agreement with the theoretical values of DFN. However, the node heads and total grid flow simulated by the FC model in two cases are either overestimated or underestimated in the coarse grid. Although elongated paths of grid flow Lg between cells are corrected by increasing the permeability coefficients of the cells, the coarse grid based on MODFLOW cannot accurately analyze the length of inclined fissures Lf on the grid. If the size of selected coarse grid is larger than the distance between two adjacent fissures, part of the two adjacent fissures will be merged in the grid, thereby shortening the actual flow path. Secondly, the length of inclined fissure Lf. is either overestimated or underestimated. The inaccurate analysis of flow paths will lead to an increase in the estimation error of flow rates and heads, especially in dense fissure zones near the center of the model domain. This effect can be greatly alleviated by refining the grid. The fissure dead-ends have non-negligible effects on MODFLOW grid flow. In this study, the MODFLOW-based FC model can be used to solve the pressure distributions between interconnected network of fissures and the rock matrix, overcoming the difficulty of determining water exchange between rock matrix and fissure media. This study can realize the expression of heterogeneity and anisotropy of fissure permeability in numerical computational units, which will promote the understanding of MODFLOW simulation of fissure flow based on connectivity.

    Furthermore, The size of the MODFLOW grid should be determined in combination with the size of the study area, which should be small enough to capture the detailed features in the flow, and large enough to reduce the computational time and cost. Generally, the grid size is determined by a trial algorithm. When FC approach is applied to the actual watershed, the spatial variability of topography and landform factors should also be considered, such as coupling surface elevation DEM data.

  • 加载中
  • [1] WILLIMAS P W. The role of the epikarst in karst and cave hydrogeology: a review[J]. International Journal of Speleology, 2008, 37(1): 1-10. doi: 10.5038/1827-806X.37.1.1

    CrossRef Google Scholar

    [2] 陈喜, 张志才. 喀斯特地区地球关键带科学与生态水文学发展综述[J]. 中国岩溶, 2022, 41(3):356-364.

    Google Scholar

    CHEN Xi, ZHANG Zhicai. An overview on the development of science and ecological hydrology of the earth critical zones in karst area[J]. Carsologica Sinica, 2022, 41(3): 356-364.

    Google Scholar

    [3] 张君, 陈洪松, 聂云鹏, 付智勇, 连晋姣, 王发, 罗紫东, 王克林. 西南喀斯特关键带结构及其水文过程研究进展[J]. 应用生态学报, 2024, 35(4):985-996.

    Google Scholar

    ZHANG Jun, CHEN Hongsong, NIE Yunpeng, FU Zhiyong, LIAN Jinjiao, WANG Fa, LUO Zidong, Wang Kelin. Research progress on structure and hydrological processes in the karst critical zone, southwest China[J]. Chinese Journal of Applied Ecology, 2024, 35(4): 985-996.

    Google Scholar

    [4] 蒋忠诚. 中国南方表层岩溶带的特征及形成机理[J]. 热带地理, 1998, 18(4):322-326. doi: 10.3969/j.issn.1001-5221.1998.04.007

    CrossRef Google Scholar

    JIANG Zhongcheng. Features of epikarst zone south China and formation mechanism[J]. Tropical Geography, 1998, 18(4): 322-326. doi: 10.3969/j.issn.1001-5221.1998.04.007

    CrossRef Google Scholar

    [5] 刘晓丽, 王恩志, 王思敬, 樊赟赟. 裂隙岩体表征方法及岩体水力学特性研究[J]. 岩石力学与工程学报, 2008, 27(9):1814-1821. doi: 10.3321/j.issn:1000-6915.2008.09.008

    CrossRef Google Scholar

    LIU Xiaoli, WANG Enzhi, WANG Sijing, FAN Yunyun. Representation method of fractured rock mass and its hydraulic properties study[J]. Chinese Journal of Rock Mechanics and Engineering, 2008, 27(9): 1814-1821. doi: 10.3321/j.issn:1000-6915.2008.09.008

    CrossRef Google Scholar

    [6] YE Zuyang, FAN Xincheng, ZHANG Jun, SHENG Jianlong, CHEN Yuting, FAN Qingli, QIN Huikai. Evaluation of connectivity characteristics on the permeability of two-dimensional fracture networks using geological entropy[J]. Water Resources Research, 2021, 57(10): 1-16.

    Google Scholar

    [7] HU Mengsu, RUTQVIST J, WANG Yuan. A practical model for fluid flow in discrete-fracture porous media by using the numerical manifold method[J]. Advances in Water Resources, 2016, 97: 38-51. doi: 10.1016/j.advwatres.2016.09.001

    CrossRef Google Scholar

    [8] 罗明明, 尹德超, 张亮, 陈植华, 周宏, 韩兆丰, 史婷婷. 南方岩溶含水系统结构识别方法初探[J]. 中国岩溶, 2015, 34(6):543-550.

    Google Scholar

    LUO Mingming, YIN Dechao, ZHANG Liang, CHEN Zhihua, ZHOU Hong, HAN Zhaofeng, SHI Tingting. Identifying methods of karst aquifer system structure in South China[J]. Carsologica Sinica, 2015, 34(6): 543-550.

    Google Scholar

    [9] 徐中平, 周训, 崔相飞, 拓明明, 王昕昀, 张颖. 岩溶区地下水数值模拟研究进展[J]. 中国岩溶, 2018, 37(4):475-483.

    Google Scholar

    XU Zhongping, ZHOU Xun, CUI Xiangfei, TA Mingming, WANG Xinyun, ZHANG Ying. Research advances of numerical simulation of groundwater in karst areas[J]. Carsologica Sinica, 2018, 37(4): 475-483.

    Google Scholar

    [10] 常勇. 裂隙-管道二元结构的岩溶泉水文过程分析与模拟[D]. 南京: 南京大学, 2015.

    Google Scholar

    CHANG Yong. Analysis and simulation of the hydrological process of the karst aquifer with fracture-conduit dual structure [D]. Nanjing: Nanjing University, 2015.

    Google Scholar

    [11] 党志文, 邵景力, 崔亚莉, 李军, 宫志强, 赵良杰, 梁永升. 基于MODFLOW-CFP的贵州大井流域岩溶地下水数值模拟[J]. 中国岩溶, 2023, 42(2):266-276.

    Google Scholar

    DANG Zhiwen, SHAO Jingli, CUI Yali, LI Jun, GONG Zhiqiang, ZHAO Liangjie, LIANG Yongsheng. Numerical simulation of karst groundwater in Dajing basin of Guizhou Province based on MODFLOW-CFP[J]. Carsologica Sinica, 2023, 42(2): 266-276.

    Google Scholar

    [12] 赵良杰, 夏日元, 杨杨, 邵景力, 易连兴, 王喆. 基于MODFLOW的岩溶管道水流模拟方法探讨与应用[J]. 中国岩溶, 2017, 36(3):346-351.

    Google Scholar

    ZHAO Liangjie, XIA Riyuan, YANG Yang, SHAO Jingli, YI Lianxing, WANG Zhe. Discussion and application of simulation methods for karst conduit flow based on MODFLOW[J]. Carsologica Sinica, 2017, 36(3): 346-351.

    Google Scholar

    [13] 王晋丽, 陈喜, 黄远洋, 张志才. 岩体裂隙网络随机生成及连通性研究[J]. 水文地质工程地质, 2013, 40(2):30-35.

    Google Scholar

    WANG Jinli, CHEN Xi, HUANG Yuanyang, ZHANG Zhicai. A study of stochastic generation and connectivity of fracture network in rock mass[J]. Hydrogeology & Engineering Geology, 2013, 40(2): 30-35.

    Google Scholar

    [14] BOTROS F E, HASSAN A E, REEVES D M, POHLL G. On mapping fracture networks onto continuum[J]. Water Resources Research, 2008, 44,W08435: 1-17.

    Google Scholar

    [15] 王晋丽, 陈喜, 常彦伟, 张志才, 康建荣, 程勤波. 基于裂隙网络连通性的裂隙流的MODFLOW模拟方法, CN108399290B[P]. 2021-07-06.

    Google Scholar

    WANG Jinli, CHEN Xi, CHANG Yanwei, ZHANG Zhicai, KANG Jianrong, CHENG Qinbo. MODFLOW simulation method for fracture flow based on fracture network connectivity, CN108399290B[P]. 2021-07-06.

    Google Scholar

    [16] 方涛, 柴军瑞, 胡海浪, 徐文彬. Monte Carlo方法在岩体裂隙结构面模拟中的应用[J]. 露天采矿技术, 2007(1):7-9+13. doi: 10.3969/j.issn.1671-9816.2007.01.003

    CrossRef Google Scholar

    FANG Tao, CHAI Junrui, HU Hailang, XU Wenbin. Application of Monte Carlo method to simulating the fracture interconnectivity of rock mass[J]. Surface Mining Technology, 2007(1): 7-9+13. doi: 10.3969/j.issn.1671-9816.2007.01.003

    CrossRef Google Scholar

    [17] LONG J C S, BILLAUX D M. From field data to fracture network modeling: An example incorporating spatial structure[J]. Water Resources Research, 1987, 23(7): 1201-1216. doi: 10.1029/WR023i007p01201

    CrossRef Google Scholar

    [18] 王晋丽, 陈喜, 常彦伟, 张志才, 康建荣. 一种基于图论的裂隙网络连通性及渗流计算的方法[P], 公告号:CN106991244B.申请(专利)号:CN201710240011.0, 2019-06-28.

    Google Scholar

    WANG Jinli, CHEN Xi, CHANG Yanwei, ZHANG Zhicai, KANG Jianrong. A method based on graph theory for fracture network connectivity and seepage calculating[P], CN106991244B. 2019-06-28.

    Google Scholar

    [19] BERKOWITZ B. Analysis of fracture network connectivity using percolation theory[J]. Mathematical Geology, 1995, 27(4): 467-483. doi: 10.1007/BF02084422

    CrossRef Google Scholar

    [20] SNOW D T. Anisotropie permeability of fractured media[J]. Water Resources Research, 1969, 5(6): 1273-1289. doi: 10.1029/WR005i006p01273

    CrossRef Google Scholar

    [21] PRIEST S D. Discontinuity analysis for rock engineering [M]. London: Chapman and Hall, 1993.

    Google Scholar

    [22] 王晋丽, 陈喜, 张志才, 康建荣, 高满. 基于离散裂隙网络模型的裂隙水渗流计算[J]. 中国岩溶, 2016, 35(4):363-371.

    Google Scholar

    WANG Jinli, CHEN Xi, ZHANG Zhicai, KANG Jianrong, GAO Man. Computation of fracture water flow based on discrete fracture network model[J]. Carsologica Sinica, 2016, 35(4): 363-371.

    Google Scholar

    [23] REEVES D M, BENSON D A, MEERSCHAERT M M. Transport of conservative solutes in simulated fracture networks: 1. Synthetic data generation[J]. Water Resources Research, 2008, 44,W05404: 1-10.

    Google Scholar

    [24] 詹美礼, 速宝玉. 交叉裂隙水流N-S方程有限元分析[J]. 水科学进展, 1997(1):4-11. doi: 10.3321/j.issn:1001-6791.1997.01.001

    CrossRef Google Scholar

    ZHAN Meili, SU Baoyu. Finite element analysis of N-S equation of intersection fissure flow[J]. Advances in Water Science, 1997(1): 4-11. doi: 10.3321/j.issn:1001-6791.1997.01.001

    CrossRef Google Scholar

    [25] 速宝玉, 詹美礼, 郭笑娥. 交叉裂隙水流的模型实验研究[J]. 水利学报, 1997(5):2-7.

    Google Scholar

    SU Baoyu, ZHAN Meili, GUO Xiaoe. Experiment research of cross fracture flow[J]. Journal of Hydraulic Engineering, 1997(5): 2-7.

    Google Scholar

  • 加载中
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Figures(9)

Tables(7)

Article Metrics

Article views(52) PDF downloads(17) Cited by(0)

Access History

Catalog

    /

    DownLoad:  Full-Size Img  PowerPoint