摘要
岩溶基岩裂隙水文研究考虑裂隙非均质和各向异性建模较难。裂隙-连续介质模型(FC)从原理上能够比较真实地刻画实际的渗流形态,是岩溶地区较理想的水文模型。基于裂隙空间形态及分布的统计特性,应用Monte Carlo模拟技术和图论,采用MATLAB程序生成二维裂隙网络原图及三级连通图,将第三级连通图映射到有限差分网格,建立离散裂隙网络(DFN)与MODFLOW相结合的FC渗流模型,编制相应的渗流模拟程序。复现交叉裂隙试验工况,及通过2个多裂隙算例分析有限差分网格分辨率以及裂隙死端对稳定渗流模拟精度的影响,以DFN渗流为准则,验证FC模型代码编写的正确性和有效性。结果表明:FC模型模拟的节点水头、总网格流、裂隙长度在粗网格的高估或低估,通过细化网格可大大消除这种影响,裂隙死端对FC模型的MODFLOW网格流模拟具有不可忽略的影响。本研究实现了裂隙渗透性非均质和各向异性在数值计算单元的表达,将促进对基于连通性的裂隙流的MODFLOW模拟理解。
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 thirdlevel 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 thirdlevel 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.
作者
王晋丽
陈喜
张志才
康建荣
胡晋山
WANG Jinli;CHEN Xi;ZHANG Zhicai;KANG Jianrong;HU Jinshan(School of Geography,Geomatics and Planning,Jiangsu Normal University,Xuzhou,Jiangsu 221116,China;School of Earth System Science,Tianjin University,Tianjin 300072,China;College of Hydrology and Water Resources,Hohai University,Nanjing,Jiangsu 210098,China)
出处
《中国岩溶》
北大核心
2025年第1期1-14,共14页
Carsologica Sinica
基金
国家自然科学基金重点项目(42030506)
国家自然科学基金面上项目(52074133)。
作者简介
第一作者:王晋丽(1970-),女,高级实验师,博士,主要从事地下水数值模拟研究。E-mail:wang_jinli70@163.com。
