首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Pressure pulsing technology is an innovative method that has been developed with the aim of overcoming preferred flow paths associated with remediation techniques that rely on the injection of reagents. Numerical and field experiments were conducted to assess how pressure pulsing affects groundwater flow and solute transport during reagent injection. A series of field experiments were performed at two field sites where a monitoring network designed to capture the breakthrough of solutes delivered from an injection well was installed. Pressure pulsing and conventional injection methods were used at each site. One site was comprised of fine sand with low heterogeneity, while the other was moderately heterogeneous with discrete layers varying from fine sand to silt. The data suggest that breakthrough was more uniform for the pressure pulsing injections; however, this difference was minor and complicated by sorption of some of the tracers employed. The groundwater flow and solute transport modeling exercise simulated the rapid boundary pressure modulation that occurs in association with pressure pulsing. Two‐dimensional (2D) simulations revealed that repeated sudden onset of injection cessation produces brief periods of gradient reversal and the development of a mixing zone near the injection well. The spatial extents of this mixing zone were found to be highly dependent upon the hydraulic diffusivity of the medium, with medium heterogeneity and pulsing frequency playing secondary roles. Three‐dimensional (3D) numerical simulations were used to benchmark the observations from one of the field sites. The results from the modeling effort showed that solute breakthrough from a pressure pulsing injection is more dispersed relative to a conventional injection as a result of the mixing zone phenomenon; however, we were unable to directly observe this mixing zone using the instrumentation deployed at the two field sites.  相似文献   

2.
The effects of the scale of measurement, i.e., the field of view, on the interpretation of fracture properties from seismic wave propagation was investigated using an acoustic lens system to produce a pseudo-collimated wavefront. The incident wavefront had a controllable beam diameter that set the field of view at 15 mm, 30 mm and 60 mm. On a smaller scale, traditional acoustic scans were used to probe the fracture in 2 mm increments. This laboratory approach was applied to two limestone samples, each containing a single induced fracture and compared to an acrylic control sample. From the analysis of the average coherent sum of the signals measured on each scale, we observed that the scale of the field of view affected the interpretation of the fracture specific stiffness. Many small-scale measurements of the seismic response of a fracture, when summed, did not predict the large-scale response of the fracture. The change from a frequency-independent to frequency-dependent fracture stiffness occurs when the scale of the field of view exceeds the spatial correlation length associated with fracture geometry. A frequency-independent fracture specific stiffness is not sufficient to classify a fracture as homogeneous. A nonuniform spatial distribution of fracture specific stiffness and overlapping geometric scales in a fracture cause a scale-dependent seismic response, which requires measurements at different field of views to fully characterize the fracture.  相似文献   

3.
While hydraulic tomography (HT) is a mature aquifer characterization technology, its applications to characterize hydrogeology of kilometer‐scale fault and fracture zones are rare. This paper sequentially analyzes datasets from two new pumping tests as well as those from two previous pumping tests analyzed by Illman et al. (2009) at a fractured granite site in Mizunami, Japan. Results of this analysis show that datasets from two previous pumping tests at one side of a fault zone as used in the previous study led to inaccurate mapping of fracture and fault zones. Inclusion of the datasets from the two new pumping tests (one of which was conducted on the other side of the fault) yields locations of the fault zone consistent with those based on geological mapping. The new datasets also produce a detailed image of the irregular fault zone, which is not available from geological investigation alone and the previous study. As a result, we conclude that if prior knowledge about geological structures at a field site is considered during the design of HT surveys, valuable non‐redundant datasets about the fracture and fault zones can be collected. Only with these non‐redundant data sets, can HT then be a viable and robust tool for delineating fracture and fault distributions over kilometer scales, even when only a limited number of boreholes are available. In essence, this paper proves that HT is a new tool for geologists, geophysicists, and engineers for mapping large‐scale fracture and fault zone distributions.  相似文献   

4.
In this work, we address the problem of characterizing the heterogeneity and uncertainty of hydraulic properties for complex geological settings. Hereby, we distinguish between two scales of heterogeneity, namely the hydrofacies structure and the intrafacies variability of the hydraulic properties. We employ multiple-point geostatistics to characterize the hydrofacies architecture. The multiple-point statistics are borrowed from a training image that is designed to reflect the prior geological conceptualization. The intrafacies variability of the hydraulic properties is represented using conventional two-point correlation methods, more precisely, spatial covariance models under a multi-Gaussian spatial law. We address the different levels and sources of uncertainty in characterizing the subsurface heterogeneity, and explore their effect on groundwater flow and transport predictions. Typically, uncertainty is assessed by way of many images, termed realizations, of a fixed statistical model. However, in many cases, sampling from a fixed stochastic model does not adequately represent the space of uncertainty. It neglects the uncertainty related to the selection of the stochastic model and the estimation of its input parameters. We acknowledge the uncertainty inherent in the definition of the prior conceptual model of aquifer architecture and in the estimation of global statistics, anisotropy, and correlation scales. Spatial bootstrap is used to assess the uncertainty of the unknown statistical parameters. As an illustrative example, we employ a synthetic field that represents a fluvial setting consisting of an interconnected network of channel sands embedded within finer-grained floodplain material. For this highly non-stationary setting we quantify the groundwater flow and transport model prediction uncertainty for various levels of hydrogeological uncertainty. Results indicate the importance of accurately describing the facies geometry, especially for transport predictions.  相似文献   

5.
Major challenges exist in delineating bedrock fracture zones because these cause abrupt changes in geological and hydrogeological properties over small distances. Borehole observations cannot sufficiently capture heterogeneity in these systems. Geophysical techniques offer the potential to image properties and processes in between boreholes. We used three‐dimensional cross borehole electrical resistivity tomography (ERT) in a 9 m (diameter) × 15 m well field to capture high‐resolution flow and transport processes in a fractured mudstone contaminated by chlorinated solvents, primarily trichloroethylene. Conductive (sodium bromide) and resistive (deionized water) injections were monitored in seven boreholes. Electrode arrays with isolation packers and fluid sampling ports were designed to enable acquisition of ERT measurements during pulsed tracer injections. Fracture zone locations and hydraulic pathways inferred from hydraulic head drawdown data were compared with electrical conductivity distributions from ERT measurements. Static ERT imaging has limited resolution to decipher individual fractures; however, these images showed alternating conductive and resistive zones, consistent with alternating laminated and massive mudstone units at the site. Tracer evolution and migration was clearly revealed in time‐lapse ERT images and supported by in situ borehole vertical apparent conductivity profiles collected during the pulsed tracer test. While water samples provided important local information at the extraction borehole, ERT delineated tracer migration over spatial scales capturing the primary hydrogeological heterogeneity controlling flow and transport. The fate of these tracer injections at this scale could not have been quantified using borehole logging and/or borehole sampling methods alone.  相似文献   

6.
Kurtzman D  Nativ R  Adar EM 《Ground water》2005,43(2):250-258
Fracture trends (defined as kilometer-scale linear features interpolated between field observations of fractures along their strikes) often have a dominant orientation. Finding a correlation between this orientation and hydraulic data could shed light on their hydraulic influence. A significant correlation between head residuals from first-order regional drift and the orientation of 2- to 4-km-long fracture trends was found in a study site in the Negev, Israel, using the semivariogram cloud analysis. Correlation of head residuals rather than the head itself implies that the orientation of the fracture trends controls the anisotropy and heterogeneity at this scale, mainly because the fracture trends define the orientation of blocks, which differ in their hydraulic properties. Preferential transmissive pathways are probably shorter than the full extent of the fracture trends, causing a relatively high head difference along the trends on the 2- to 4-km scale. Fracture trend density and additional data from short-range hydraulic tests helped characterize two blocks separated by a fault zone. The identification of hydraulic features on a kilometer scale is necessary for better modeling of regional ground water flow and transport. Hydraulic tests at this scale are not feasible, thereby rendering combined analyses of head and structural data, such as the one presented here, essential.  相似文献   

7.
Traditionally a streambed is treated as a layer of uniform thickness and low saturated hydraulic conductivity (K) in surface‐ and ground‐water studies. Recent findings have shown a high level of spatial heterogeneity within a streambed and such heterogeneity directly affects surface‐ and ground‐water exchange and can have ecological implications for biogeochemical transformations, nutrient cycling, organic matter decomposition, and reproduction of gravel spawning fish. In this study a detailed field investigation of K was conducted in two selected sites in Touchet River, a typical salmon spawning stream in arid south eastern Washington, USA. In‐stream slug tests were conducted to determine K following the Bouwer and Rice method. For the upper and lower sites, each 50 m long and 9 m wide and roughly 20 m apart, a sampling grid of 5 m longitudinally and 3 m transversely was used. The slug tests were performed for each horizontal coordinate at 0·3–0·45, 0·6–0·75, 0·9–1·05 and 1·2–1·35 m depth intervals unless a shallower impenetrable obstruction was encountered. Additionally, water levels were measured to obtain vertical hydraulic gradient (VHG) between each two adjacent depth intervals. Results indicated that K ranged over three orders of magnitude at both the upper and lower sites and differed between the two sites. At the upper site, K did not differ significantly among different depth intervals based on nonparametric statistical tests for mean, median, and empirical cumulative distribution, but the spatial pattern of K varied among different depth intervals. At the lower site, K for the 0·3–0·45 m depth interval differed statistically from those at other depth intervals, and no similar spatial pattern was found among different depth intervals. Zones of upward and downward water flow based on VHG also varied among different depth intervals, reflecting the complexities of the water flow regime. Detailed characterization of the streambed as attempted in this study should be helpful in providing information on spatial variations of streambed hydraulic properties as well as surface‐ and ground‐water interaction. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

8.
Fluid conductivity and elastic properties in fractures depend on the aperture geometry – in particular, the roughness of fracture surfaces. In this study, we have characterized the surface roughness with a log-normal distribution and investigated the transport and flow behaviour of the fractures with varying roughness characteristics. Numerical flow and transport simulations have been performed on a single two-dimensional fracture surface, whose aperture geometry changes with different variances and correlation lengths in each realization. We have found that conventional measurement of hydraulic conductivity alone is insufficient to determine these two parameters. Transient transport measurements, such as the particle breakthrough time, provide additional constraints to the aperture distribution. Nonetheless, a unique solution to the fracture aperture distribution is still under-determined with both hydraulic conductivity and transport measurements. From numerical simulations at different compression states, we have found that the flow and transport measurements exhibit different rates of changes with respect to changes in compression. Therefore, the fracture aperture distribution could be further constrained by considering the flow and transport properties under various compression states.  相似文献   

9.
Cross-borehole flowmeter tests have been proposed as an efficient method to investigate preferential flowpaths in heterogeneous aquifers, which is a major task in the characterization of fractured aquifers. Cross-borehole flowmeter tests are based on the idea that changing the pumping conditions in a given aquifer will modify the hydraulic head distribution in large-scale flowpaths, producing measurable changes in the vertical flow profiles in observation boreholes. However, inversion of flow measurements to derive flowpath geometry and connectivity and to characterize their hydraulic properties is still a subject of research. In this study, we propose a framework for cross-borehole flowmeter test interpretation that is based on a two-scale conceptual model: discrete fractures at the borehole scale and zones of interconnected fractures at the aquifer scale. We propose that the two problems may be solved independently. The first inverse problem consists of estimating the hydraulic head variations that drive the transient borehole flow observed in the cross-borehole flowmeter experiments. The second inverse problem is related to estimating the geometry and hydraulic properties of large-scale flowpaths in the region between pumping and observation wells that are compatible with the head variations deduced from the first problem. To solve the borehole-scale problem, we treat the transient flow data as a series of quasi-steady flow conditions and solve for the hydraulic head changes in individual fractures required to produce these data. The consistency of the method is verified using field experiments performed in a fractured-rock aquifer.  相似文献   

10.
— A numerical fracture flow simulation based on the lubrication approximation is used to investigate the influence of roughness on the flow inside a rough fracture, at low Reynolds number. Facing surfaces are described as self-affine topographies with identical roughness magnitude. Resolution of the Reynolds equation is achieved using two distinct numerical schemes, with consistency. Fracture closure is studied assuming perfect plastic contact between facing surfaces. Long-range correlations are shown to exist in the local aperture field due to the fracture geometry and subsequently in the local fluxes inside the fracture. Flow channeling is the result of these correlations in terms of spatial distribution of the flow, and is responsible for either flow-enhancing or flow-inhibiting behavior of the fracture. Matching between the two surfaces at scales larger than a mismatch scale is studied. The mismatch scale is the upper limit scale for the local apertures scale invariance. It appears to control flow channeling and the related dispersion of the possible behaviors over a large statistics of fractures with identical statistical features. Hydraulic anisotropy of a given fracture is investigated: the dependence of the fracture transmittivity on the pressure drop orientation is proved to be sinusoidal, with an amplitude that is controlled by the mismatch scale.  相似文献   

11.
Indirect bedload transport measurements have been made with the Swiss plate geophone system in five gravel‐bed mountain streams. These geophone sensors record the motion of bedload particles transported over a steel plate mounted flush with the channel bed. To calibrate the geophone system, direct bedload transport measurements were undertaken simultaneously. At the Erlenbach in Switzerland, a moving‐basket sampler was used. At the Fischbach and Ruetz streams in Austria, a Helley–Smith type bedload sampler provided the calibration measurements. A Bunte‐type bedload trap was used at the Rofenache stream in Austria. At the Nahal Eshtemoa in Israel, Reid‐type slot bedload samplers were used. To characterize the response of the geophone signal to bedload particles impacting on the plate, geophone summary values were calculated from the raw signal and stored at one second intervals. The number of impulses, i.e. the number of peaks above a pre‐defined threshold value of the geophone output signal, correlated well with field measured gravel transport loads and was found to be a robust parameter. The relations of impulses to gravel transport loads were generally near‐linear, but the steepness of the calibration relations differed from site to site. By comparing the calibration measurements from the different field sites and utilizing insights gained during preliminary flume experiments, it has been possible to identify the main factors that are responsible for site specific differences in the calibration coefficient. The analysis of these calibration measurements indicates that the geophone signal also contains some information about the grain size distribution of bedload. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

12.
Outcrop studies reveal a common occurrence of tabular zones of significantly‐increased fracture intensity affecting otherwise well‐lithified rocks. These zones, called fracture corridors, can have a profound effect on multi‐phase fluid flow in the subsurface. Using standard geo‐modelling tools, it is possible to generate 3D realizations of reservoirs that contain distributions of such fracture corridors that are consistent with observations, including the vertical frequency in pseudo‐wells inserted into the model at random locations. These models can generate the inputs to flow simulation. The approach adopted here is to run the flow simulations in a single‐porosity representation where the flow effects of fractures are upscaled into equivalent cell‐based properties, preserving a clear spatial relationship between the input geology and the resulting cellular model. The simulated reservoir performance outcomes are very similar to those seen in real oilfields: extreme variability between wells, early water breakthrough, disappointing recoveries and patchy saturation distributions. Thus, a model based on fracture corridors can provide an explanation for the observed flow performance of a suitable field. However, the use of seismics to identify fracture corridors is not an easy task. New work is needed to predict the seismic responses of fracture corridor systems to be able to judge whether it is likely that we can robustly detect and characterize these flow‐significant features adequately.  相似文献   

13.
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.  相似文献   

14.
Accurate quantification of in situ heterogeneity and flow processes through fractured geologic media remains elusive for hydrogeologists due to the complexity in fracture characterization and its multiscale behavior. In this research, we demonstrated the efficacy of tracer-electrical resistivity tomography (ERT) experiments combined with numerical simulations to characterize heterogeneity and delineate preferential flow paths in a fractured granite aquifer. A series of natural gradient saline tracer experiments were conducted from a depth window of 18 to 22 m in an injection well (IW) located inside the Indian Institute of Technology Hyderabad campus. Tracer migration was monitored in a time-lapse mode using two cross-sectional surface ERT profiles placed in the direction of flow gradient. ERT data quality was improved by considering stacking, reciprocal measurements, resolution indicators, and geophysical logs. Dynamic changes in subsurface electrical properties inferred via resistivity anomalies were used to highlight preferential flow paths of the study area. Temporal changes in electrical resistivity and tracer concentration were monitored along the vertical in an observation well located at 48 m to the east of the IW. ERT-derived tracer breakthrough curves were in agreement with geochemical sample measurements. Fracture geometry and hydraulic properties derived from ERT and pumping tests were further used to evaluate two mathematical conceptualizations that are relevant to fractured aquifers. Results of numerical analysis conclude that dual continuum model that combines matrix and fracture systems through a flow exchange term has outperformed equivalent continuum model in reproducing tracer concentrations at the monitoring wells (evident by a decrease in RMSE from 199 to 65 mg/L). A sensitivity analysis on model simulations conclude that spatial variability in hydraulic conductivity, local-scale dispersion, and flow exchange at fracture-matrix interface have a profound effect on model simulations.  相似文献   

15.
In general, the accuracy of numerical simulations is determined by spatial and temporal discretization levels. In fractured porous media, the time step size is a key factor in controlling the solution accuracy for a given spatial discretization. If the time step size is restricted by the relatively rapid responses in the fracture domain to maintain an acceptable level of accuracy in the entire simulation domain, the matrix tends to be temporally over-discretized. Implicit sub-time stepping applies smaller sub-time steps only to the sub-domain where the accuracy requirements are less tolerant and is most suitable for problems where the response is high in only a small portion of the domain, such as within and near the fractures in fractured porous media. It is demonstrated with illustrative examples that implicit sub-time stepping can significantly improve the simulation efficiency with minimal loss in accuracy when simulating flow and transport in fractured porous media. The methodology is successfully applied to density-dependent flow and transport simulations in a Canadian Shield environment, where the flow and transport is dominated by discrete, highly conductive fracture zones.  相似文献   

16.
  相似文献   

17.
A quasi-three-dimensional particle tracking model is developed to characterize the spatial and temporal effects of advection, molecular diffusion, Taylor dispersion, fracture wall deposition, matrix diffusion, and co-transport processes on two discrete plumes (suspended monodisperse or polydisperse colloids and dissolved contaminants) flowing through a variable aperture fracture situated in a porous medium. Contaminants travel by advection and diffusion and may sorb onto fracture walls and colloid particles, as well as diffuse into and sorb onto the surrounding porous rock matrix. A kinetic isotherm describes contaminant sorption onto colloids and sorbed contaminants assume the unique transport properties of colloids. Sorption of the contaminants that have diffused into the matrix is governed by a first-order kinetic reaction. Colloids travel by advection and diffusion and may attach onto fracture walls; however, they do not penetrate the rock matrix. A probabilistic form of the Boltzmann law describes filtration of both colloids and contaminants on fracture walls. Ensemble-averaged breakthrough curves of many fracture realizations are used to compare arrival times of colloid and contaminant plumes at the fracture outlet. Results show that the presence of colloids enhances contaminant transport (decreased residence times) while matrix diffusion and sorption onto fracture walls retard the transport of contaminants. Model simulations with the polydisperse colloids show increased effects of co-transport processes.  相似文献   

18.
Amplitudes of seismic waves increase significantly as they pass through soft soil layers near the earth's surface. This phenomenon, commonly known as site amplification, is a major factor influencing the extent of damage on structures. It is crucial that site amplification is accounted for when designing structures on soft soils. The characteristics of site amplification at a given site can be estimated by analytical models, as well as field tests. Analytical models require as inputs the geometry of all soil layers from surface to bedrock, their dynamic properties (e.g. density, wave velocity, damping), and the incident bedrock motions. Field tests involve recording and analyzing the dynamic response of sites to artificial excitations, ambient forces, and actual earthquakes. The most reliable estimates of site amplification are obtained by analyzing the recorded motions of the site during strong earthquakes. This paper presents a review of the types and the generating mechanisms of site amplification, and the models and methods that are used to characterize them from earthquake records.  相似文献   

19.
I. Haltas 《水文研究》2012,26(22):3448-3458
Recognizing the spatial heterogeneity of hydraulic parameters, many researchers have studied the solute transport by both groundwater and channel flow in a stochastic framework. One of the methodologies used to up‐scale the stochastic solute transport equation, from a point‐location scale to a grid scale, is the cumulant expansion method combined with the calculus for the time‐ordered exponential and the calculus for the Lie operator. When the point‐location scale transport equation is scaled up to the grid scale, using the cumulant expansion method, a new dispersion coefficient emerges in the dispersive term of the solute transport equation in addition to the molecular dispersion coefficient. This velocity driven dispersion is called ‘macrodispersion’. The macrodispersion coefficient is the integral function of the time‐ordered covariance of the random velocity field. The integral is calculated over a Lagrangian trajectory of the flow. The Lagrangian trajectory depends on the following: (i) the spatial origin of the particle; (ii) the time when the macrodispersion is calculated; and (iii) the mean velocity field along the trajectory itself. The Lagrangian trajectory is a recursive function of time because the location of the particle along the trajectory at a particular time depends on the location of the particle at the previous time. This recursive functional form of the Lagrangian trajectory makes the calculation of the macrodispersion coefficient difficult. Especially for the unsteady, spatially non‐stationary, non‐uniform flow field, the macrodispersion coefficient is a highly complex expression and, so far, calculated using numerical methods in the discrete domains. Here, an analytical method was introduced to calculate the macrodispersion coefficient in the discrete domain for the unsteady and steady, spatially non‐stationary flow cases accurately and efficiently. This study can fill the gap between the theory of the ensemble averaged solute transport model and its numerical implementations. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

20.
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.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号