共查询到20条相似文献，搜索用时 31 毫秒
1.
Understanding the Geometry of Connected Fracture Flow with Multiperiod Oscillatory Hydraulic Tests 
下载免费PDF全文
下载免费PDF全文 An understanding of the spatial and hydraulic properties of fast preferential flow pathways in the subsurface is necessary in applications ranging from contaminant fate and transport modeling to design of energy extraction systems. One method for the characterization of fracture properties over interwellbore scales is Multiperiod Oscillatory Hydraulic (MOH) testing, in which the aquifer response to oscillatory pressure stimulations is observed. MOH tests were conducted on isolated intervals of wells in siliciclastic and carbonate aquifers in southern Wisconsin. The goal was to characterize the spatial properties of discrete fractures over interwellbore scales. MOH tests were conducted on two discrete fractured intervals intersecting two boreholes at one field site, and a nest of three piezometers at another field site. Fracture diffusivity estimates were obtained using analytical solutions that relate diffusivity to observed phase lag and amplitude decay. In addition, MOH tests were used to investigate the spatial extent of flow using different conceptual models of fracture geometry. Results indicated that fracture geometry at both field sites can be approximated by permeable two‐dimensional fracture planes, oriented near‐horizontally at one site, and near‐vertically at the other. The technique used on MOH field data to characterize fracture geometry shows promise in revealing fracture network characteristics important to groundwater flow and transport. 相似文献
2.
Studying the Flow Dynamics of a Karst Aquifer System with an Equivalent Porous Medium Model 总被引：2，自引：0，他引：2
The modeling of groundwater flow in karst aquifers is a challenge due to the extreme heterogeneity of its hydraulic parameters and the duality in their discharge behavior, that is, rapid response of highly conductive karst conduits and delayed drainage of the low‐permeability fractured matrix after recharge events. There are a number of different modeling approaches for the simulation of the karst groundwater dynamics, applicable to different aquifer as well as modeling problem types, ranging from continuum models to double continuum models to discrete and hybrid models. This study presents the application of an equivalent porous model approach (EPM, single continuum model) to construct a steady‐state numerical flow model for an important karst aquifer, that is, the Western Mountain Aquifer Basin (WMAB), shared by Israel and the West‐Bank, using MODFLOW2000. The WMAB was used as a catchment since it is a well‐constrained catchment with well‐defined recharge and discharge components and therefore allows a control on the modeling approach, a very rare opportunity for karst aquifer modeling. The model demonstrates the applicability of equivalent porous medium models for the simulation of karst systems, despite their large contrast in hydraulic conductivities. As long as the simulated saturated volume is large enough to average out the local influence of karst conduits and as long as transport velocities are not an issue, EPM models excellently simulate the observed head distribution. The model serves as a starting basis that will be used as a reference for developing a long‐term dynamic model for the WMAB, starting from the pre‐development period (i.e., 1940s) up to date. 相似文献
3.
Fractures in porous media have been documented extensively. However, they are often omitted from groundwater flow and mass transport models due to a lack of data on fracture hydraulic properties and the computational burden of simulating fractures explicitly in large model domains. We present a MATLAB toolbox, FracKfinder, that automates HydroGeoSphere (HGS), a variably saturated, control volume finite-element model, to simulate an ensemble of discrete fracture network (DFN) flow experiments on a single cubic model mesh containing a stochastically generated fracture network. Because DFN simulations in HGS can simulate flow in both a porous media and a fracture domain, this toolbox computes tensors for both the matrix and fractures of a porous medium. Each model in the ensemble represents a different orientation of the hydraulic gradient, thus minimizing the likelihood that a single hydraulic gradient orientation will dominate the tensor computation. Linear regression on matrices containing the computed three-dimensional hydraulic conductivity (K) values from each rotation of the hydraulic gradient is used to compute the K tensors. This approach shows that the hydraulic behavior of fracture networks can be simulated where fracture hydraulic data are limited. Simulation of a bromide tracer experiment using K tensors computed with FracKfinder in HGS demonstrates good agreement with a previous large-column, laboratory study. The toolbox provides a potential pathway to upscale groundwater flow and mass transport processes in fractured media to larger scales. 相似文献
4.
Abstract Saltwater intrusion is a naturally occurring phenomenon that is exacerbated significantly by excessive groundwater exploitation in coastal aquifers. In order to determine the extent of saltwater intrusion in a karstified coastal aquifer in Crete, Greece, a three-dimensional, density-dependent groundwater flow and transport model was developed and compared to the more traditional sharp-interface approach. The karstified medium was modelled using a combination of the equivalent porous medium approach (for lower-order fractures) and a discrete fracture approach (for the main fractures/faults). The model takes into consideration the geomorphologic characteristics of the karstic system, such as the depth and orientation of the fault network, and the diffusion phenomena associated with the variable densities of freshwater and saltwater—parameters that create a complex system, inducing uncertainty in the model. The model results showed that the orientation of the fractures, the pumping activity and the fluid density effects drive the seawater intrusion front asymmetrically inland. Editor Z.W. Kundzewicz Citation Dokou, Z. and Karatzas, G.P., 2012. Saltwater intrusion estimation in a karstified coastal system using density-dependent modelling and comparison with the sharp-interface approach. Hydrological Sciences Journal, 57 (5), 985–999. 相似文献
5.
Evaluation of the MODFLOW-2005 Conduit Flow Process 总被引：1，自引：0，他引：1
The recent development of the Conduit Flow Process (CFP) by the U.S. Geological Survey (USGS) provides hydrogeologic modelers with a new tool that incorporates the non-Darcian, multiporosity components of flow characteristic of karst aquifers. CFP introduces new parameters extending beyond those of traditional Darcian groundwater flow codes. We characterize a karst aquifer to collect data useful for evaluating this new tool at a test site in west-central Florida, where the spatial distribution and cross-sectional area of the conduit network are available. Specifically, we characterize: (1) the potential for Darcian/non-Darcian flow using estimates of specific discharge vs. observed hydraulic gradients, and (2) the temporal variation for the direction and magnitude of fluid exchange between the matrix and conduit network during extreme hydrologic events. We evaluate the performance of CFP Mode 1 using a site-scale dual-porosity model and compare its performance with a comparable laminar equivalent continuum model (ECM) using MODFLOW-2005. Based on our preliminary analyses, hydraulic conductivity coupled with conduit wall conductance improved the match between observed and simulated discharges by 12% to 40% over turbulent flow alone (less than 1%). 相似文献
6.
Bioremediation in Fractured Rock: 1. Modeling to Inform Design,Monitoring, and Expectations 下载免费PDF全文
下载免费PDF全文 Claire R. Tiedeman Allen M. Shapiro Paul A. Hsieh Thomas E. Imbrigiotta Daniel J. Goode Pierre J. Lacombe Mary F. DeFlaun Scott R. Drew Carole D. Johnson John H. Williams Gary P. Curtis 《Ground water》2018,56(2):300-316
Field characterization of a trichloroethene (TCE) source area in fractured mudstones produced a detailed understanding of the geology, contaminant distribution in fractures and the rock matrix, and hydraulic and transport properties. Groundwater flow and chemical transport modeling that synthesized the field characterization information proved critical for designing bioremediation of the source area. The planned bioremediation involved injecting emulsified vegetable oil and bacteria to enhance the naturally occurring biodegradation of TCE. The flow and transport modeling showed that injection will spread amendments widely over a zone of lower‐permeability fractures, with long residence times expected because of small velocities after injection and sorption of emulsified vegetable oil onto solids. Amendments transported out of this zone will be diluted by groundwater flux from other areas, limiting bioremediation effectiveness downgradient. At nearby pumping wells, further dilution is expected to make bioremediation effects undetectable in the pumped water. The results emphasize that in fracture‐dominated flow regimes, the extent of injected amendments cannot be conceptualized using simple homogeneous models of groundwater flow commonly adopted to design injections in unconsolidated porous media (e.g., radial diverging or dipole flow regimes). Instead, it is important to synthesize site characterization information using a groundwater flow model that includes discrete features representing high‐ and low‐permeability fractures. This type of model accounts for the highly heterogeneous hydraulic conductivity and groundwater fluxes in fractured‐rock aquifers, and facilitates designing injection strategies that target specific volumes of the aquifer and maximize the distribution of amendments over these volumes. 相似文献
7.
As a result of rock dissolution processes, karst aquifers exhibit highly conductive features such as caves and conduits. Within these structures, groundwater flow can become turbulent and therefore be described by nonlinear gradient functions. Some numerical groundwater flow models explicitly account for pipe hydraulics by coupling the continuum model with a pipe network that represents the conduit system. In contrast, the Conduit Flow Process Mode 2 (CFPM2) for MODFLOW-2005 approximates turbulent flow by reducing the hydraulic conductivity within the existing linear head gradient of the MODFLOW continuum model. This approach reduces the practical as well as numerical efforts for simulating turbulence. The original formulation was for large pore aquifers where the onset of turbulence is at low Reynolds numbers (1 to 100) and not for conduits or pipes. In addition, the existing code requires multiple time steps for convergence due to iterative adjustment of the hydraulic conductivity. Modifications to the existing CFPM2 were made by implementing a generalized power function with a user-defined exponent. This allows for matching turbulence in porous media or pipes and eliminates the time steps required for iterative adjustment of hydraulic conductivity. The modified CFPM2 successfully replicated simple benchmark test problems. 相似文献
8.
A tracer test in a carbonate aquifer is analyzed using the method of moments and two analytical advection-dispersion models (ADMs) as well as a numerical model. The numerical model is a coupled continuum-pipe flow and transport model that accounts for two different flow components in karstified carbonate aquifers, i.e., rapid and often turbulent conduit flow and Darcian flow in the fissured porous rock. All techniques employed provide reasonable fits to the tracer breakthrough curve (TBC) measured at a spring. The resulting parameter estimates are compared to investigate how each conceptual model of flow and transport processes that forms the basis of the analyses affects the interpretation of the tracer test. Numerical modeling results suggest that the method of moments and the analytical ADMs tend to overestimate the conduit volume because part of the water discharged at the spring is wrongly attributed to the conduit system if flow in the fissured porous rock is ignored. In addition, numerical modeling suggests that mixing of the two flow components accounts for part of the dispersion apparent in the measured TBC, while the remaining part can be attributed to Taylor dispersion. These processes, however, cannot reasonably explain the tail of the TBC. Instead, retention in immobile-fluid regions as included in a nonequilibrium ADM provides a possible explanation. 相似文献
9.
The aim of this research was to refine the actual conceptual model related to the activation of high‐altitude temporary springs within the carbonate Apennines in southern Italy. The research was carried out through geophysical, hydrogeological, hydrochemical and isotopic investigations at the Acqua dei Faggi experimental site during five hydrologic years. The research demonstrated that, in carbonate aquifers where low‐permeability faults cause the aquifer system to be compartmentalized, high‐altitude temporary springs may be recharged by groundwater. In such settings, neither surface water infiltration in karst systems nor perched temporary aquifers play a role of utmost importance. The rare (once or a few time a year) activation of such springs is due to the fact that groundwater unusually reach the threshold head that allows the spring to flow. The activation of the studied high‐altitude temporary spring also depended on relationships between a low‐permeability fault core and a karst system that locally interrupts the low‐permeability barrier. In fact, when the hydraulic head did not reach the karst system, the concentrated head loss within the fault core did not allow the spring to flow, because the groundwater entirely flowed through the fault towards the downgradient compartment. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
10.
Alluvial fans are potential sites of potable groundwater in many parts of the world. Characteristics of alluvial fans sediments are changed radially from high energy coarse-grained deposition near the apex to low energy fine-grained deposition downstream so that patchy wedge-shaped aquifers with radial heterogeneity are formed. The hydraulic parameters of the aquifers (e.g. hydraulic conductivity and specific storage) change in the same fashion. Analytical or semi-analytical solutions of the flow in wedge-shaped aquifers are available for homogeneous cases. In this paper we derive semi-analytical solutions of groundwater flow to a well in multi-zone wedge-shaped aquifers. Solutions are provided for three wedge boundary configurations namely: constant head–constant head wedge, constant head–barrier wedge and barrier–barrier wedge. Derivation involves the use of integral transforms methods. The effect of heterogeneity ratios of zones on the response of the aquifer is examined. The results are presented in form of drawdown and drawdown derivative type curves. Heterogeneity has a significant effect on over all response of the pumped aquifer. Solutions help understanding the behavior of heterogeneous multi-zone aquifers for sustainable development of the groundwater resources in alluvial fans. 相似文献
11.
Stephen J. Gale 《Journal of Hydrology》1984,70(1-4):309-327
Analysis of the morphology of dissolution-bedform assemblages and hydraulically-transported sediments found within conduits in carbonate aquifers in North America and the British Isles has allowed the hydraulic conditions under which conduit flows occur to be established. Mean values of flow velocity, boundary-shear stress, conduit Reynolds number, conduit Froude number, boundary friction factor, boundary roughness and flow power have been calculated. The values obtained are in agreement with other evidence in the conduits and are comparable to those obtained from other, similar, hydraulic systems. 相似文献
12.
A systematic investigation of the effect of configurations of stochastically distributed fracture networks on hydraulic behaviour for fractured rock masses could provide either quantitative or qualitative correlation between the structural configuration of the fracture network and its corresponding hydraulic behaviour, and enhance our understanding of appropriate application of groundwater flow and contaminant transport modelling in fractured rock masses. In this study, the effect of block sizes, intersection angles of fracture sets, standard deviations of fracture orientation, and fracture densities on directional block hydraulic conductivity and representative elementary volume is systematically investigated in two dimensions by implementing a numerical discrete fracture fluid flow model and incorporating stochastically distributed fracture configurations. It is shown from this investigation that the configuration of a stochastically distributed fracture network has a significant quantitative or qualitative effect on the hydraulic behaviour of fractured rock masses. Compared with the deterministic fracture configurations that have been extensively dealt with in a previous study, this investigation is expected to be more practical and adequate, since fracture geometry parameters are inherently stochastically distributed in the field. Moreover, the methodology and approach presented in this study may be generally applied to any fracture system in investigating the hydraulic behaviours from configurations of the fracture system while establishing a ‘bridge’ from the discrete fracture network flow modelling to equivalent continuum modelling in fractured rock masses. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
13.
One of the challenges in characterizing fractured-rock aquifers is determining whether the equivalent porous medium approximation is valid at the problem scale. Detailed hydrogeologic characterization completed at a small study site in a densely fractured dolomite has yielded an extensive data set that was used to evaluate the utility of the continuum and discrete-fracture approaches to aquifer characterization. There are two near-vertical sets of fractures at the site; near-horizontal bedding-plane partings constitute a third fracture set. Eighteen boreholes, including five coreholes, were drilled to a depth of approximately 10.6 m. Borehole geophysical logs revealed several laterally extensive horizontal fractures and dissolution zones. Flowmeter and short-interval packer testing identified which of these features were hydraulically important. A monitoring system, consisting of short-interval piezometers and multilevel samplers, was designed to monitor four horizontal fractures and two dissolution zones. The resulting network consisted of >70 sampling points and allowed detailed monitoring of head distributions in three dimensions. Comparison of distributions of hydraulic head and hydraulic conductivity determined by these two approaches suggests that even in a densely fractured-carbonate aquifer, a characterization approach using traditional long-interval monitoring wells is inadequate to characterize ground water movement for the purposes of regulatory monitoring or site remediation. In addition, traditional multiwell pumping tests yield an average or bulk hydraulic conductivity that is not adequate for predicting rapid ground water travel times through the fracture network, and the pumping test response does not appear to be an adequate tool for assessing whether the porous medium approximation is valid. 相似文献
14.
15.
The groundwater flow in a fissured chalky environment at the northern border of the Paris Basin depends on several geological and hydrogeological parameters. Although the studied sector of the basin presents a homogeneous rock type, it is affected by a fracture network. In this type of environment, in which the permeability is low, the groundwater flow displays significant disruption, which is localized in the Fruges region (northern France). The interconnection of the discontinuities (network of fault and/or joints) is reliant on the structural control of groundwater flow through increases in the hydraulic connection between the unsaturated and the saturated zone. The methodology developed herein makes use of microstructural and regional analysis of the fracture patterns, and allowed consideration of the piezometric variations of the chalk aquifer during periods of low and high groundwater levels (April and October 2001) and a diagraphic representation of the estimated physical parameters (electrical resistivity). This enabled us to construct a ‘flow structure’ conceptual model in which we identify two types of faults: tight walls and flow paths that control the piezometric heads and the flow rate. Model validation was carried out on a similar sector. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
16.
Tsuyoshi Nohara Hidemi Tanaka Kunio Watanabe Noboru Furukawa Akira Takami 《Island Arc》2006,15(4):537-545
Abstract The spatial hydrogeological and structural character of the active Mozumi-Sukenobu Fault (MSF) was investigated along a survey tunnel excavated through the MSF in the Kamioka Mine, central Japan. Major groundwater conduits on both sides of the MSF are recognized. One is considered to be a subvertical conduit between the tunnel and the surface, and the other is estimated to be a major reservoir of old meteoric water alongside the MSF. Our studies indicate that part of the MSF is a sub-vertical continuous barrier that obstructs younger meteoric water observed in the south-eastern part of the Active Fault Survey Tunnel (AFST) and water recharge to the rock mass intersected by the north-western part of the AFST. The MSF might be a continuous barrier resulting in the storage of a large quantity of older groundwater to the northwest. The observations and results of in situ hydraulic tests indicate that the major reservoir is not the fault breccia associated with the northeast–southwest trending faults of the MSF, but the zone in which blocks of fractured rocks occur beside high angle faults corresponding to X shears whose tectonic stress field coincides with the present regional stress field and antithetic Riedel shears of the MSF. The results from borehole investigations in the AFST indicate that secondary porosity is developed in the major reservoir due to the destruction of filling minerals and fracture development beside these shears. The increase in hydraulic conductivity is not directly related to increased density of fractures around the MSF. Development of secondary porosity could cause the increase in hydraulic conductivity around the MSF. Our results suggest that minor conduits of the fracture network are sporadically distributed in the sedimentary rocks around the MSF in the AFST. 相似文献
17.
Lucia Mastrorillo Michele Saroli Stefano Viaroli Francesca Banzato Daniela Valigi Marco Petitta 《水文研究》2020,34(5):1167-1181
A sustained increase in spring discharges was monitored after the 2016 Central Italy seismic sequence in the fractured carbonate aquifer of Valnerina–Sibillini Mts. The groundwater surplus recorded between August 2016 and November 2017 was determined to be between 400 and 500 × 106 m3. In fractured aquifers, the post-seismic rise in spring discharges is generally attributed to an increase in bulk permeability caused by the fracture cleaning effect, which is induced by pore pressure propagation. In the studied aquifers, the large amount of additional discharge cannot only be attributed to the enhanced permeability, which was evaluated to be less than 20% after each main seismic event. A detailed analysis of the spring discharge hydrographs and of the water level at five gauging stations was carried out to determine the possible causes of this sudden increase in groundwater outflow. Taking into account the geological and structural framework, a conceptual model of a basin-in-series has been adopted to describe the complex hydrogeological setting, where the thrusts and extensional faults have clearly influenced the groundwater flow directions before and after the seismic sequence. The prevalent portion of the total post-seismic discharge surplus not explained by the increase in permeability has been attributed to changes in the hydraulic gradient that caused seismogenic fault rupture and the disruption in the upgradient sector of the aquifer. The additional flow calculated through the breach of the pre-existing hydrostructural barrier corresponds to approximately 470 × 106 m3. This value is consistent with the total discharge increase measured in the whole study area, validating the proposed conceptual model. Consequently, a shift in the piezometric divide of the hydrogeological system has been induced, causing a potentially permanent change that lowers the discharge amount of the eastern springs. 相似文献
18.
In variably confined carbonate platforms, impermeable confining units collect rainfall over large areas and deliver runoff to rivers or conduits in unconfined portions of platforms. Runoff can increase river stage or conduit heads in unconfined portions of platforms faster than local infiltration of rainfall can increase groundwater heads, causing hydraulic gradients between rivers, conduits and the aquifer to reverse. Gradient reversals cause flood waters to flow from rivers and conduits into the aquifer where they can dissolve limestone. Previous work on impacts of gradient reversals on dissolution has primarily emphasized individual caves and little research has been conducted at basin scales. To address this gap in knowledge, we used legacy data to assess how a gradient of aquifer confinement across the Suwannee River Basin, north‐central Florida affected locations, magnitudes and processes of dissolution during 2005–2007, a period with extreme ranges of discharge. During intense rain events, runoff from the confining unit increased river stage above groundwater heads in unconfined portions of the platform, hydraulically damming inputs of groundwater along a 200 km reach of river. Hydraulic damming allowed allogenic runoff with SICAL < ?4 to fill the entire river channel and flow into the aquifer via reversing springs. Storage of runoff in the aquifer decreased peak river discharges downstream and contributed to dissolution within the aquifer. Temporary storage of allogenic runoff in karst aquifers represents hyporheic exchange at a scale that is larger than found in streams flowing over non‐karst aquifers because conduits in karst aquifers extend the area available for exchange beyond river beds deep into aquifers. Post‐depositional porosity in variably confined carbonate platforms should thus be enhanced along rivers that originate on confining units. This distribution should be considered in models of porosity distribution used to manage water and hydrocarbon resources in carbonate rocks. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
19.
In cases when an equivalent porous medium assumption is suitable for simulating groundwater flow in bedrock aquifers, estimation of seepage into underground mine workings (UMWs) can be achieved by specifying MODFLOW drain nodes at the contact between water bearing rock and dewatered mine openings. However, this approach results in significant numerical problems when applied to simulate seepage into an extensive network of UMWs, which often exist at the mine sites. Numerical simulations conducted for individual UMWs, such as a vertical shaft or a horizontal drift, showed that accurate prediction of seepage rates can be achieved by either applying grid spacing that is much finer than the diameter/width of the simulated openings (explicit modeling) or using coarser grid with cell sizes exceeding the characteristic width of shafts or drifts by a factor of 3. Theoretical insight into this phenomenon is presented, based on the so-called well-index theory. It is demonstrated that applying this theory allows to minimize numerical errors associated with MODFLOW simulation of seepage into UMWs on a relatively coarse Cartesian grid. Presented examples include simulated steady-state groundwater flow from homogeneous, heterogeneous, and/or anisotropic rock into a vertical shaft, a horizontal drift/cross-cut, a ramp, two parallel drifts, and a combined system of a vertical shaft connected to a horizontal drift. 相似文献
20.
Numerical simulations of variable-density flow and solute transport have been conducted to investigate dense plume migration for various configurations of 2D fracture networks. For orthogonal fractures, simulations demonstrate that dispersive mixing in fractures with small aperture does not stabilize vertical plume migration in fractures with large aperture. Simulations in non-orthogonal 2D fracture networks indicate that convection cells form and that they overlap both the porous matrix and fractures. Thus, transport rates in convection cells depend on matrix and fracture flow properties. A series of simulations in statistically equivalent networks of fractures with irregular orientation show that the migration of a dense plume is highly sensitive to the geometry of the network. If fractures in a random network are connected equidistantly to the solute source, few equidistantly distributed fractures favor density-driven transport. On the other hand, numerous fractures have a stabilizing effect, especially if diffusive transport rates are high. A sensitivity analysis for a network with few equidistantly distributed fractures shows that low fracture aperture, low matrix permeability and high matrix porosity impede density-driven transport because these parameters reduce groundwater flow velocities in both the matrix and the fractures. Enhanced molecular diffusion slows down density-driven transport because it favors solute diffusion from the fractures into the low-permeability porous matrix where groundwater velocities are smaller. For the configurations tested, variable-density flow and solute transport are most sensitive to the permeability and porosity of the matrix, which are properties that can be determined more accurately than the geometry and hydraulic properties of the fracture network, which have a smaller impact on density-driven transport. 相似文献