岩溶泉域降雨径流水文过程的模拟——以重庆金佛山水房泉为例

通过地貌、水文地质、土壤、植被、地下洞穴管道等实地调查,以及示踪试验确定水房泉泉域范围和地下水文系统特征,并通过代表性点的土壤入渗试验、降雨和流量监测以及DEM数据等获得泉域水文模型所需的面积、结点高程、入渗率、糙率、管道长度、含水层孔隙度、出流系数等参数。选择SCS径流曲线模型估算地表产流,利用SWMM模型模拟泉域对场降雨的径流响应过程。通过运行模型,与实测流量比较,结果显示模拟曲线与观测流量曲线吻合较好,用于校正和验证的两场降雨产流的模拟误差分别为9.5%和12%,表明SWMM可应用于岩溶区以管道流为主要排水系统的含水介质的模拟。     

Mise-a-la-masse method,tracer tests and cave detection analysis of underground-river karst-conduit distribution and structure: a case study of the Dafengdong conduits in Guizhou Province,China

In karst areas, groundwater movement is dominated by conduit flow; thus, understanding the distribution and structure of karst conduits has great significance for water research, groundwater protection and engineering construction. With the Dafengdong (Guizhou Province, China) underground river karst as an example, a mise-a-la-masse method, tracer tests and cave detection were conducted to study the distribution and structure of karst conduits. Combined with information on the geological and hydrogeological conditions, the geological factors that form karst conduits were determined. It can be concluded that: (1) Under the influence of faults, karst conduits usually develop towards or along faults; (2) multiple karst conduits form easily on both sides of tectonic fracture zones; (3) both lithology and geological structure affect the formation of blue holes, and when the directions of karst conduits change, blue holes easily form in weak parts of the structure; (4) at springs where two aquifers intersect, with strong chemical dissolution of the lower karst aquifer and mechanical collapse of the upper aquifer, blue holes also form easily; and (5) integrated mise-a-la-masse method, tracer tests and cave detection can accurately discern the distribution and structure of karst conduits. Geological factors can be used to preliminarily delineate the distribution and structure of karst conduits in similar areas based on hydrogeological conditions. Such methods hold great significance for groundwater extraction and protection and engineering construction in karst areas.

     

岩溶管道结构影响泉流量变化的数值模拟研究

采用MODFLOW-CFP建立管道流数值模型,模拟岩溶含水系统在暴雨期的响应过程,以泉口流量峰值作为因变量,分别对控制管道结构的4个参数以及落水洞集中补给比例进行调整,研究管道结构如何控制泉流量变化过程。结果表明:管道直径、管壁渗透系数和落水洞补给比例对泉口流量峰值均表现为正相关,其中管壁渗透系数影响最大;管道弯曲度值较小时其增大促进出口流量峰值变大,当径流途径变长引起的流量减小幅度超过管壁面积变大引起的流量增大幅度时,流量峰值逐渐变小;管壁粗糙度在达到两个相对应的雷诺数临界值时水流状态发生突变,导致流量峰值也表现为两次突变。      

Modeling hydrological responses of karst spring to storm events: example of the Shuifang spring (Jinfo Mt., Chongqing,China)

Reproduction of hydrographs at karst springs has been an approach of understanding the karst aquifer, which normally acts as drains for the groundwater flow. However, its numerical modeling is difficult since factors for the internal geometry and connectedness are unknown and hard to quantify. Hydrographs of the karst aquifer with well-developed conduits in Shuifang spring catchment were obtained from the automatic gauging station at the spring orifice. Data as to the conduit system were also obtained based on results and analyses of tracer tests. With these data, the hydrological responses of Shuifang spring to storm events were simulated by storm water management model (SWMM) developed by USA EPA (Environmental Protection Agency). Nash–Sutcliffe efficiencies are used to compare the computed flow to the observed, which are 0.95 and 0.92 for calibration and validation. SWMM was verified and applicable in karst conduit drainage system. The model illustrated correctly quick recharge through conduits and slow and low inflow from the fissured aquifer matrix. The SCS-CN (soil conservation service-curve number) infiltration method was used for computation of losses and runoff. Field tests indicated that permeability was extremely high but different in karst area, which was less sensitive to the computed runoff when exceeded the common value provided by SWMM. Therefore, an improved quantitative infiltration model for karst area will make SWMM possible to be a useful tool for assessing and reproducing spring hydrographs.     

The comparison of hydrodynamic characteristics of karst aquifers: application on two karst formations in Zagros (Asmari and Ilam-Sarvak), southwest Iran

In order to study the function, hydrodynamic behavior, and hydraulic properties of the karst aquifers in Izeh, southwest Iran, time series analysis was applied to the precipitation, spring discharge, and piezometric head data of two representative karst systems of Zagros (Ilam-Sarvak and Asmari Formations). Time series analysis was applied to two karst aquifers, those of Asmari and Ilam-Sarvak Formations. The daily precipitations of anticlines were estimated based on the precipitation–elevation function which was applied on digital elevation model (DEM) of the area. The mean estimated daily precipitations were considered in bivariate time series analysis as input data of each karst system. The total length of time series was about 2.7 years, extending from May 2007 to December 2009. During the research, several one-parameter probe data loggers were installed, which daily measure the water surfaces in karst aquifers. Time series analysis was applied for describing Izeh karst aquifers with a focus on both univariate (autocorrelation and spectral analysis) and bivariate (cross-correlation, gain function, and coherency function) methods. The results show that Asmari karst aquifer in Kamarderaz Anticline has a large storage capacity. Because of lacking a well-organized karst network, in the Asmari karst aquifer, baseflow dominates with low contribution of quick-flow. In the Ilam-Sarvak karst aquifer (Shavish and Tanush Anticlines), the karstification occurred in fractures and small diameter conduits, which caused to quick-flow between dense limestone. The Ilam-Sarvak karst aquifer could be regarded as a transition between two extreme types of karst, e.g., highly karstified system and in the opposite, extremely diffused one. The analysis of well hydrograph in Ilam-Sarvak karst aquifer shows that the karst aquifer has a low storage capacity. Unlike Asmari karst aquifer, the fractures and small diameter conduits in Ilam-Sarvak karst aquifer are more enhanced, producing a better developed karst network. Contrary to the typical karst systems, however, diffuse flow and conduit flow coexist in the Asmari Formation.     

Application of a solute transport model under variable velocity conditions in a conduit flow aquifer: Olalde karst system, Basque Country, Spain

A model based on numerical solutions, which allows for solving the dispersion equation under variable recharge and velocity conditions, is developed to simulate solute transport in conduit flow aquifers during flow recession periods. As an example, the evolution of a tracer in the little known karst conduit that links the sinking stream of Oma valley to the Olalde spring is investigated in the karstic region of Santa Eufemia-Ereñozar (Basque Country, Spain). The model, with different hypothetical structures, allows for obtaining series of tracer breakthrough curves, which are fitted to experimental data using an optimization algorithm. These results, although they can be used to simulate the tracer evolution between the two points considered, do not allow for determining the internal structure and spatial disposition of contributions in the aquifer.     

Comparison of a karst groundwater model with and without discrete conduit flow

Karst aquifers exhibit a dual flow system characterized by interacting conduit and matrix domains. This study evaluated the coupled continuum pipe-flow framework for modeling karst groundwater flow in the Madison aquifer of western South Dakota (USA). Coupled conduit and matrix flow was simulated within a regional finite-difference model over a 10-year transient period. An existing equivalent porous medium (EPM) model was modified to include major conduit networks whose locations were constrained by dye-tracing data and environmental tracer analysis. Model calibration data included measured hydraulic heads at observation wells and estimates of discharge at four karst springs. Relative to the EPM model, the match to observation well hydraulic heads was substantially improved with the addition of conduits. The inclusion of conduit flow allowed for a simpler hydraulic conductivity distribution in the matrix continuum. Two of the high-conductivity zones in the EPM model, which were required to indirectly simulate the effects of conduits, were eliminated from the new model. This work demonstrates the utility of the coupled continuum pipe-flow method and illustrates how karst aquifer model parameterization is dependent on the physical processes that are simulated.     

西南某岩溶区地下水系统示踪试验与解析

以西南某典型岩溶区为例,解析示踪试验在岩溶管道连通性以及获取水文地质参数中的应用。选择落水洞为投放点,分别从落水洞西侧和东侧寻找地下水出露点作为接收点,判别落水洞地下径流的实际去向以及落水洞与接收点之间的水力联系。结合Qtracer2软件对示踪试验成果进行定量解析,确定示踪剂回收率、地下水平均流速、最快流速,估算出岩溶管道结构特征和水文地质参数。结果表明:落水洞与接收点JS01、JS03之间不存在直接水力联系;落水洞与接收点JS02存在水力联系且岩溶管道极为发育,含水介质不均匀,地下水运移路径较为通畅,为典型的紊流流态;落水洞地下径流的主要方向是由西向东,但在丰水期雨量较大期间,接收点JS04能够接收到荧光素钠,说明丰水期水位上涨后两者间会有水力联系,导致部分水量向落水洞西侧排泄。     

Assessment of key transport parameters in a karst system under different dynamic conditions based on tracer experiments: the Jeita karst system,Lebanon

Artificial tracer experiments were conducted in the mature karst system of Jeita (Lebanon) under various flow conditions using surface and subsurface tracer injection points, to determine the variation of transport parameters (attenuation of peak concentration, velocity, transit times, dispersivity, and proportion of immobile and mobile regions) along fast and slow flow pathways. Tracer breakthrough curves (TBCs) observed at the karst spring were interpreted using a two-region nonequilibrium approach (2RNEM) to account for the skewness in the TBCs’ long tailings. The conduit test results revealed a discharge threshold in the system dynamics, beyond which the transport parameters vary significantly. The polynomial relationship between transport velocity and discharge can be related to the variation of the conduit’s cross-sectional area. Longitudinal dispersivity in the conduit system is not a constant value (α?=?7–10 m) and decreases linearly with increasing flow rate because of dilution effects. Additionally, the proportion of immobile regions (arising from conduit irregularities) increases with decreasing water level in the conduit system. From tracer tests with injection at the surface, longitudinal dispersivity values are found to be large (8–27 m). The tailing observed in some TBCs is generated in the unsaturated zone before the tracer actually arrives at the major subsurface conduit draining the system. This work allows the estimation and prediction of the key transport parameters in karst aquifers. It shows that these parameters vary with time and flow dynamics, and they reflect the geometry of the flow pathway and the origin of infiltrating (potentially contaminated) recharge.     

Simulation of the development of karst aquifers: role of the epikarst

 The evolution of a karst aquifer is modelled taking into account the karst groundwater flow as well as the dissolution kinetics of calcite. In particular, infiltration of water from the epikarst is simulated which controls the temporal and spatial distribution of recharge to the phreatic zone. The results show that the evolution of karst conduits is initiated in the spring. The existence of preferential flow paths leads to the evolution of highly conductive so-called dendritic cave systems, i.e., single passages which concentrate the flow and drain the catchment. With time, the amount of undersaturated water flowing directly into the conduit system is increased leading to an acceleration of the conduits enlargement. Three phases are identified for the evolution of karst aquifers: (a) an initiation stage; (b) an enlargement stage; and (c) a stagnation phase. Received: 4 August 1997 / Accepted: 19 January 1998     

A Leaky-Conduit Model of Transient Flow in Karstic Aquifers

Karst Flow Model (KFM) simulates transient flow in an unconfined karstic aquifer having a well-developed conduit system. KFM treats the springshed as a two-dimensional porous matrix containing a triangulated irregular network of leaky conduits. The number and location of conduits can be specified arbitrarily, perhaps using field information as a guide, or generated automatically. Conduit networks can be tree-like or braided. Rainwater that has infiltrated down from the surface leaks into the conduits from the adjacent porous matrix at a rate dictated by Darcy’s law, then flows turbulently to the spring via the conduits. KFM is calibrated using the known steady state; geometry and recharge determine the steady fluxes in the conduits, and the head distribution determines conduit gradients and sizes. Spring flow can vary with time due to spatially and temporally variable recharge and due to prescribed variations in the elevation of the spring. KFM is illustrated by four examples run on a test aquifer consisting of 27 nodes, 42 elements, and 26 conduits. Three examples (drought, uniform rainstorm, storm-water input to one element) are simulations, while the fourth uses data from a spring-basin flooding event. The qualitative fit between the predicted and observed spring discharge in the fourth example provides support of the hypothesis that the dynamic behavior of a karst conduit system is an emergent property of a self-organized system, largely independent of the locations and properties of individual conduits.     

Regional-scale analysis of karst underground flow deduced from tracing experiments: examples from carbonate aquifers in Malaga province,southern Spain

Tracer concentration data from field experiments conducted in several carbonate aquifers (Malaga province, southern Spain) were analyzed following a dual approach based on the graphical evaluation method (GEM) and solute transport modeling to decipher flow mechanisms in karst systems at regional scale. The results show that conduit system geometry and flow conditions are the principal factors influencing tracer migration through the examined karst flow routes. Solute transport is mainly controlled by longitudinal advection and dispersion throughout the conduit length, but also by flow partitioning between mobile and immobile fluid phases, while the matrix diffusion process appears to be less relevant. The simulation of tracer breakthrough curves (BTCs) suggests that diffuse and concentrated flow through the unsaturated zone can have equivalent transport properties under extreme recharge, with high flow velocities and efficient mixing due to the high hydraulic gradients generated. Tracer mobilization within the saturated zone under low flow conditions mainly depends on the hydrodynamics (rather than on the karst conduit development), which promote a lower longitudinal advection and retardation in the tracer migration, resulting in a marked tailing effect of BTCs. The analytical advection-dispersion equation better approximates the effective flow velocity and longitudinal dispersion estimations provided by the GEM, while the non-equilibrium transport model achieves a better adjustment of most asymmetric and long-tailed BTCs. The assessment of karst underground flow properties from tracing tests at regional scale can aid design of groundwater management and protection strategies, particularly in large hydrogeological systems (i.e. transboundary carbonate aquifers) and/or in poorly investigated ones.     

Relative importance and chemical effects of diffuse and focused recharge in an eogenetic karst aquifer: an example from the unconfined upper Floridan aquifer, USA

Karst aquifer studies often focus on allogenic water inputs and large conduit flow. However, diffuse recharge can be significant, particularly in unconfined eogenetic karst aquifers that retain high matrix permeability. This study examines an unconfined region of the upper Floridan aquifer (USA) that hosts a sinking stream, its resurgence, and a large conduit system. Daily diffuse recharge was approximated using a water-budget method and ranged from 17% of precipitation during a low precipitation year to >53% during the highest precipitation year, illustrating the highly variable nature of diffuse recharge in this region. The total allogenic input via the sinking stream over the 5 years of the study was significantly larger than the volume of diffuse recharge. However, only about 2% of the allogenic recharge flows from the conduit into the surrounding aquifer. That flow is restricted to storm events when hydraulic heads in the conduits exceed those in the surrounding aquifer. The estimated volume of dissolution is similar for allogenic recharge and diffuse recharge to the unconfined region surrounding the conduits, but dissolution from the diffuse recharge is distributed over a larger area than dissolution from allogenic recharge. These results exemplify how recharge type impacts flow and water–rock interactions in eogenetic karst aquifers.     

Applying a modified conduit flow process to understand conduit-matrix exchange of a karst aquifer

Due to the high heterogeneity and complexity of water flow movement for multiple karst water-bearing mediums,the evaluation,effective development,and utilization of karst water resources are significantly limited.Matrix flow is usually laminar,whereas conduit flow is usually turbulent.The driving mechanisms of water exchange that occur between the karst conduit and its adjacent matrix are not well understood.This paper investigates the hydrodynamic characteristics and the mechanism of flow exchange in dual water-bearing mediums(conduit and matrix)of karst aquifers through laboratory experimentation and numerical simulation.A karst aquifer consisting of a matrix network and a conduit was proposed,and the relationship between the water exchange flux and hydraulic head differences generated from the laboratory experiments was analyzed.Two modes of experimental tests were performed with different fixed water level boundaries in the laboratory karst aquifer.The results indicate that the water exchange capacity was proportional to the square root of hydraulic head differences.The linear exchange term in the conduit flow process(CFP)source program was modified according to experimental results.The modified CFP and the original CFP model experimental data results were compared,and it was found that the modified CFP model had better fitting effects.These results showed that the water exchange mechanism between conduit and matrix is very important for solid-liquid interface reaction,water resource evaluation,and understanding of karst hydrodynamic behavior.     

基于MODFLOW的岩溶管道水流模拟方法探讨与应用

针对岩溶管道水流特征,选取广西寨底岩溶流域为研究对象,分别采用三维有限差分地下水流动模型(MODFLOW)中的Drain和River模块概化模拟岩溶管道特征,结合管道内上游、中游及下游三个观测孔对比分析两种方法的适用性,结果表明:末期等水位线图中位于管道两侧处流线均有形状突变,符合实际管道水流动态特征,且能模拟出管道内水位变化趋势。最后分析两种方法模拟管道水流的原理,Drain模块中岩溶管道仅起排水作用,不允许管道内水流流向周围含水层,而River模块允许岩溶管道水与周围含水层水流交换,因而相对Drain模块而言,应用River模块概化模拟岩溶管道更加精确。       

Simulating flow in karst aquifers at laboratory and sub-regional scales using MODFLOW-CFP

Groundwater flow in a well-developed karst aquifer dominantly occurs through bedding planes, fractures, conduits, and caves created by and/or enlarged by dissolution. Conventional groundwater modeling methods assume that groundwater flow is described by Darcian principles where primary porosity (i.e. matrix porosity) and laminar flow are dominant. However, in well-developed karst aquifers, the assumption of Darcian flow can be questionable. While Darcian flow generally occurs in the matrix portion of the karst aquifer, flow through conduits can be non-laminar where the relation between specific discharge and hydraulic gradient is non-linear. MODFLOW-CFP is a relatively new modeling program that accounts for non-laminar and laminar flow in pipes, like karst caves, within an aquifer. In this study, results from MODFLOW-CFP are compared to those from MODFLOW-2000/2005, a numerical code based on Darcy’s law, to evaluate the accuracy that CFP can achieve when modeling flows in karst aquifers at laboratory and sub-regional (Woodville Karst Plain, Florida, USA) scales. In comparison with laboratory experiments, simulation results by MODFLOW-CFP are more accurate than MODFLOW 2005. At the sub-regional scale, MODFLOW-CFP was more accurate than MODFLOW-2000 for simulating field measurements of peak flow at one spring and total discharges at two springs for an observed storm event.     

Inference of creep mechanism in underground soil loss of karst conduits I. Conceptual model

Karst is widely distributed in the southwest of China, especially in Guizhou Province. The phenomenon of desertification in these areas is very serious. And soil erosion is the key link in the process of desertification. Through field monitoring, underground soil leakage is derived to the main mode of soil loss in this area. Shear strength tests and creep experiments were carried out with the aim of analyzing the creep mechanism in underground soil loss. It is shown that the water content can lead to the great influence on the shear strength of the brown clay. This variation has been combined with creep characteristics besides the structural geology, hydrology condition and brown clay distribution circumstance (field observation). A conceptual creep model of the brown clay sliding along the karst conduits has been unveiled to show the detailed inference of the creep mechanism in the underground soil loss: geology and hydrology control the development of the karst conduit system; and penetration of water induces the weakening of the shear strength of the surface soil and accelerates the creeping and sliding of the brown clay along the karst conduit system. This understanding of the creep mechanism has significant implications for the future management of the soil erosion in the karst area.     

Development of a Darcy-Brinkman model to simulate water flow and tracer transport in a heterogeneous karstic aquifer (Val d’Orléans, France)

Darcy’s law is the equation of reference widely used to model aquifer flows. However, its use to model karstic aquifers functioning with large pores is problematic. The physics occurring within the karstic conduits requires the use of a more representative macroscopic equation. A hydrodynamic model is presented which is adapted to the karstic aquifer of the Val d’Orléans (France) using two flow equations: (1) Darcy’s law, used to describe water flow within the massive limestone, and (2) the Brinkman equation, used to model water flow within the conduits. The flow equations coupled with the transport equation allow the prediction of the karst transfer properties. The model was tested by using six dye tracer tests and compared to a model that uses Darcy’s law to describe the flow in karstic conduits. The simulations show that the conduit permeability ranges from 5?×?10^?6 to 5.5?×?10^?5?m² and the limestone permeability ranges from 8?×?10^?11 to 6?×?10^?10?m². The dispersivity coefficient ranges from 23 to 53 m in the conduits and from 1 to 5 m in the limestone. The results of the simulations carried out using Darcy’s law in the conduits show that the dispersion towards the fractures is underestimated.     

Characterization of flow and transport dynamics in karst aquifers by analyzing tracer test results in conduits and recharge areas (the Egino Massif,Basque Country,Spain): environmental and management implications

Karst aquifers contribute to supplying drinking water to almost a quarter of the world´s population. Their complex dynamics requires specific approaches aimed at recognizing their singularities, analyzing its vulnerability, and ensuring water resources quality. In this paper, the results of processing and modeling five breakthrough tracer curves obtained under different hydrodynamic conditions in the main conduit of Egino karst aquifer (Basque Country, Spain) are analyzed together with those involving pressure injections of the tracer in the saturated zone of the karst massif recharge area. In the conduit, transport is immediate and highly efficient (recovery rates above 84% and dispersion coefficients from 15.04 to 84.35 m²/min); tracer retentions increase as flow rates decrease and no significant contributions to its surroundings are observed. In contrast, tracer transport from the massif recharge area is more complex: after injection at a pressure of 1 MPa, most tracer remains in the surrounding of the injection borehole, retained in a saturated medium of low effective fracture porosity (? _f ?=?1.02?×?10^?4, assuming a radial divergent flow model); subsequently, the main tracer mobilization to the spring was registered with the first rains, with 0.088 m/min mean velocity and high concentrations per unit mass being injected (C _p /M₀?=?0.03 mg/L/kg), which is evidence that the tracer reaches soon the karst conduit network. In any case, a decreasing tracer presence is registered at the injection zone during a hydrological cycle. In both cases, the observed non-linearity of transport processes should be considered in the development of vulnerability approaches, modeling efforts, and mapping. Furthermore, in the case of karst massif recharge areas, as the presence of pollutants may have a significant impact on the springs and persist over time, their management and protection needs must be revised in each specific site. Simultaneously, quality-monitoring programs at the springs must be adapted to the aquifers recognized dynamics.