首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
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.  相似文献   

2.
The objective of this study was to characterize the behavior of a groundwater contaminant (trichloroethene, TCE) plume after implementation of a source‐containment operation at a site in Arizona. The plume resides in a quasi‐three‐layer system comprising a sand/gravel unit bounded on the top and bottom by relatively thick silty clayey layers. The system was monitored for 60 months beginning at start‐up in 2007 to measure the change in contaminant concentrations within the plume, the change in plume area, the mass of the contaminant removed, and the integrated contaminant mass discharge (CMD). The concentrations of TCE in groundwater pumped from the plume extraction wells have declined significantly over the course of operation, as have concentrations for groundwater sampled from 40 monitoring wells located within the plume. The total CMD associated with operation of the plume extraction wells peaked at 0.23 kg/d, decreased significantly within 1 year, and thereafter began an asymptotic decline to a current value of approximately 0.03 kg/d. Despite an 87% reduction in contaminant mass and a comparable 87% reduction in CMD for the plume, the spatial area encompassed by the plume has decreased by only approximately 50%. This is much less than would be anticipated based on ideal flushing and mass‐removal behavior. Simulations produced with a simplified three‐dimensional (3D) numerical model matched reasonably well to the measured data. The results of the study suggest that permeability heterogeneity, back diffusion, hydraulic factors associated with the specific well field system, and residual discharge from the source zone are all contributing to the observed persistence of the plume, as well as the asymptotic behavior currently observed for mass removal and for the reduction in CMD.  相似文献   

3.
Numerical modeling was employed to study the performance of thermal conductive heating (TCH) in fractured shale under a variety of hydrogeological conditions. Model results show that groundwater flow in fractures does not significantly affect the minimum treatment zone temperature, except near the beginning of heating or when groundwater influx is high. However, fracture and rock matrix properties can significantly influence the time necessary to remove all liquid water (i.e., reach superheated steam conditions) in the treatment area. Low matrix permeability, high matrix porosity, and wide fracture spacing can contribute to boiling point elevation in the rock matrix. Consequently, knowledge of these properties is important for the estimation of treatment times. Because of the variability in boiling point throughout a fractured rock treatment zone and the absence of a well-defined constant temperature boiling plateau in the rock matrix, it may be difficult to monitor the progress of thermal treatment using temperature measurements alone.  相似文献   

4.
Hydraulic tomography for detecting fracture zone connectivity   总被引：1，自引：0，他引：1  
Hao Y  Yeh TC  Xiang J  Illman WA  Ando K  Hsu KC  Lee CH 《Ground water》2008,46(2):183-192
Fracture zones and their connectivity in geologic media are of great importance to ground water resources management as well as ground water contamination prevention and remediation. In this paper, we applied a recently developed hydraulic tomography (HT) technique and an analysis algorithm (sequential successive linear estimator) to synthetic fractured media. The application aims to explore the potential utility of the technique and the algorithm for characterizing fracture zone distribution and their connectivity. Results of this investigation showed that using HT with a limited number of wells, the fracture zone distribution and its connectivity (general pattern) can be mapped satisfactorily although estimated hydraulic property fields are smooth. As the number of wells and monitoring ports increases, the fracture zone distribution and connectivity become vivid and the estimated hydraulic properties approach true values. We hope that the success of this application may promote the development and application of the new generations of technology (i.e., hydraulic, tracer, pneumatic tomographic surveys) for mapping fractures and other features in geologic media.  相似文献   

5.
The spontaneous imbibition of water and other liquids into gas-filled fractures in variably-saturated porous media is important in a variety of engineering and geological contexts. However, surprisingly few studies have investigated this phenomenon. We present a theoretical framework for predicting the 1-dimensional movement of water into air-filled fractures within a porous medium based on early-time capillary dynamics and spreading over the rough surfaces of fracture faces. The theory permits estimation of sorptivity values for the matrix and fracture zone, as well as a dispersion parameter which quantifies the extent of spreading of the wetting front. Quantitative data on spontaneous imbibition of water in unsaturated Berea sandstone cores were acquired to evaluate the proposed model. The cores with different permeability classes ranging from 50 to 500 mD and were fractured using the Brazilian method. Spontaneous imbibition in the fractured cores was measured by dynamic neutron radiography at the Neutron Imaging Prototype Facility (beam line CG-1D, HFIR), Oak Ridge National Laboratory. Water uptake into both the matrix and the fracture zone exhibited square-root-of-time behavior. The matrix sorptivities ranged from 2.9 to 4.6 mm s−0.5, and increased linearly as the permeability class increased. The sorptivities of the fracture zones ranged from 17.9 to 27.1 mm s−0.5, and increased linearly with increasing fracture aperture width. The dispersion coefficients ranged from 23.7 to 66.7 mm2 s−1 and increased linearly with increasing fracture aperture width and damage zone width. Both theory and observations indicate that fractures can significantly increase spontaneous imbibition in unsaturated sedimentary rock by capillary action and surface spreading on rough fracture faces. Fractures also increase the dispersion of the wetting front. Further research is needed to investigate this phenomenon in other natural and engineered porous media.  相似文献   

6.
Matrix diffusion can attenuate the rate of plume migration in fractured bedrock relative to the rate of ground water flow for both conservative and nonconservative solutes of interest. In a system of parallel, equally spaced constant aperture fractures subject to steady-state ground water flow and an infinite source width, the degree of plume attenuation increases with time and travel distance, eventually reaching an asymptotic level. The asymptotic degree of plume attenuation in the absence of degradation can be predicted by a plume attenuation factor, beta, which is readily estimated as R' (phi(m)/phi(f)), where R' is the retardation factor in the matrix, phi(m) is the matrix porosity, and phi(f) is the fracture porosity. This dual-porosity relationship can also be thought of as the ratio of primary to secondary porosity. Beta represents the rate of ground water flow in fractures relative to the rate of plume advance. For the conditions examined in this study, beta increases with greater matrix porosity, greater matrix fraction organic carbon, larger fracture spacing, and smaller fracture aperture. These concepts are illustrated using a case study where dense nonaqueous phase liquid in fractured sandstone produced a dissolved-phase trichloroethylene (TCE) plume approximately 300 m in length. Transport parameters such as matrix porosity, fracture porosity, hydraulic gradient, and the matrix retardation factor were characterized at the site through field investigations. In the fractured sandstone bedrock examined in this study, the asymptotic plume attenuation factors (beta values) for conservative and nonconservative solutes (i.e., chloride and TCE) were predicted to be approximately 800 and 12,210, respectively. Quantitative analyses demonstrate that a porous media (single-porosity) solute transport model is not appropriate for simulating contaminant transport in fractured sandstone where matrix diffusion occurs. Rather, simulations need to be conducted with either a discrete fracture model that explicitly incorporates matrix diffusion, or a dual-continuum model that accounts for mass transfer between mobile and immobile zones. Simulations also demonstrate that back diffusion from the matrix to fractures will likely be the time-limiting factor in reaching ground water cleanup goals in some fractured bedrock environments.  相似文献   

7.
Deep-well injection into fractured sandstone is an option for the disposal of contaminated mine dewatering discharge from an open pit uranium mine. As part of the assessment of potential contaminant migration from deep-well injection, the effect of matrix diffusion was evaluated. An analytical mathematical model was developed for the simulation of the radial movement of a contaminant front away from an injection point under steady flow conditions in a planar fracture with uniform properties. The model includes the effects of advection in the fracture, diffusion of contaminants from the fracture into the rock matrix, and equilibrium adsorption on the fracture surface as well as in the rock matrix. Effective diffusion coefficients obtained from laboratory experiments on 11 intact core samples varied from 3.4 × 10−8 to 3.2 × 10−7 cm2/s. Model simulations were made with diffusion coefficient values in this range and with single-fracture injection rates estimated from fracture frequencies in boreholes, and from bulk hydraulic conductivity values obtained from field tests. Because of matrix diffusion, the rate of outward movement of the front of the nonreactive contaminants from the injection well is much slower than the rate of water flow in the fractures. Simulations of the movement of contaminants that undergo adsorption indicate that even a small distribution coefficient for the rock matrix causes the contaminants to remain very close to the injection well during the one-year period. The results of the simplified model demonstrate that matrix diffusion is an important process that cannot be neglected in the assessment of a waste disposal scheme located in fractured porous rock. However, in order to make a definitive assessment of the capability of matrix diffusion and associated matrix adsorption to significantly limit the extent of contaminant migration around injection wells, it would be necessary to conduct field tests such as a preliminary or experimental injection.  相似文献   

8.
Boreholes drilled through contaminated zones in fractured rock create the potential for vertical movement of contaminated ground water between fractures. The usual assumption is that purging eliminates cross contamination; however, the results of a field study conducted in a trichloroethylene (TCE) plume in fractured sandstone with a mean matrix porosity of 13% demonstrates that matrix-diffusion effects can be strong and persistent. A deep borehole was drilled to 110 m below ground surface (mbgs) near a shallow bedrock well containing high TCE concentrations. The borehole was cored continuously to collect closely spaced samples of rock for analysis of TCE concentrations. Geophysical logging and flowmetering were conducted in the open borehole, and a removable multilevel monitoring system was installed to provide hydraulic-head and ground water samples from discrete fracture zones. The borehole was later reamed to complete a well screened from 89 to 100 mbgs; persistent TCE concentrations at this depth ranged from 2100 to 33,000 microg/L. Rock-core analyses, combined with the other types of borehole information, show that nearly all of this deep contamination was due to the lingering effects of the downward flow of dissolved TCE from shallower depths during the few days of open-hole conditions that existed prior to installation of the multilevel system. This study demonstrates that transfer of contaminant mass to the matrix by diffusion can cause severe cross contamination effects in sedimentary rocks, but these effects generally are not identified from information normally obtained in fractured-rock investigations, resulting in potential misinterpretation of site conditions.  相似文献   

9.
We use particle tracking to determine contributing areas (CAs) to wells for transient flow models that simulate cyclic domestic pumping and extreme recharge events in a small synthetic watershed underlain by dipping sedimentary rocks. The CAs consist of strike-oriented bands at locations where the water table intersects high-hydraulic conductivity beds, and from which groundwater flows to the pumping well. Factors that affect the size and location of the CAs include topographic flow directions, rock dip direction, cross-bed fracture density, and position of the well relative to streams. For an effective fracture porosity (ne) of 10−4, the fastest advective travel times from CAs to wells are only a few hours. These results indicate that wells in this type of geologic setting can be highly vulnerable to contaminants or pathogens flushed into the subsurface during extreme recharge events. Increasing ne to 10−3 results in modestly smaller CAs and delayed well vulnerability due to slower travel times. CAs determined for steady-state models of the same setting, but with long-term average recharge and pumping rates, are smaller than CAs in the models with extreme recharge. Also, the earliest-arriving particles arrive at the wells later in the steady-state models than in the extreme-recharge models. The results highlight the importance of characterizing geologic structure, simulating plausible effective porosities, and simulating pumping and recharge transience when determining CAs in fractured rock aquifers to assess well vulnerability under extreme precipitation events.  相似文献   

10.
This study presents the impact of fractures on CO2 transport, capillary pressure and storage capacity by conducting both experimental and numerical studies. A series of laboratory experiment tests was designed with both a homogeneous and a fractured core under CO2 storage conditions. The experimental results reveal a piston-like brine displacement with gravity override effects in the homogeneous core regardless of CO2 injection rates. In the fractured core, however, two distinctive types of brine displacements were observed; one showing brine displacement only in the fracture whereas the other shows brine displacement both in the fracture and matrix with different rates, which were dependent on the magnitude of the pressure build-up in the matrix. The injectivity in the fractured core was twice of the homogeneous core, while the amount of calculated CO2 in the homogeneous core was over 1.5 times greater than the fractured core. Salt precipitation, which is likely to occur near injection wells, was observed in the experiments; X-ray images enabled the observation of salt-precipitation during CO2-flooding tests. Finally, numerical simulations predict free-phase CO2 transfer between fracture and matrix in a fracture-matrix system. Pressure gradients between the fracture and matrix enforced CO2 to transfer from the fracture into matrix at the front of the CO2 plume, whereas, the reversal of pressure gradients at the rear zone of the CO2 plume reversed the transfer process. The variation of CO2 saturation within the fracture was caused by fracture aperture variations, and local variations of fracture permeability control the free-phase CO2 transfer between the fracture and matrix.  相似文献   

11.
The hydrogeologic framework of fractured sedimentary bedrock at the former Naval Air Warfare Center (NAWC), 1 1 For more information on the USGS Toxics Substances Hydrology Program at the Naval Air Warfare Center visit the NAWC website at http://nj.usgs.gov/nawc/
Trenton, New Jersey, a trichloroethylene (TCE)-contaminated site in the Newark Basin, is developed using an understanding of the geologic history of the strata, gamma-ray logs, and rock cores. NAWC is the newest field research site established as part of the U.S. Geological Survey Toxic Substances Hydrology Program, Department of Defense (DoD) Strategic Environmental Research and Development Program, and DoD Environmental Security Technology Certification Program to investigate contaminant remediation in fractured rock. Sedimentary bedrock at the NAWC research site comprises the Skunk Hollow, Byram, and Ewing Creek Members of the Lockatong Formation and Raven Rock Member of the Stockton Formation. Muds of the Lockatong Formation that were deposited in Van Houten cycles during the Triassic have lithified to form the bedrock that is typical of much of the Newark Basin. Four lithotypes formed from the sediments include black, carbon-rich laminated mudstone, dark-gray laminated mudstone, light-gray massive mudstone, and red massive mudstone. Diagenesis, tectonic compression, off-loading, and weathering have altered the rocks to give some strata greater hydraulic conductivity than other strata. Each stratum in the Lockatong Formation is 0.3 to 8 m thick, strikes N65°E, and dips 25° to 70°NW. The black, carbon-rich laminated mudstone tends to fracture easily, has a relatively high hydraulic conductivity and is associated with high natural gamma-ray count rates. The dark-gray laminated mudstone is less fractured and has a lower hydraulic conductivity than the black carbon-rich laminated mudstone. The light-gray and the red massive mudstones are highly indurated and tend to have the least fractures and a low hydraulic conductivity. The differences in gamma-ray count rates for different mudstones allow gamma-ray logs to be used to correlate and delineate the lithostratigraphy from multiple wells. Gamma-ray logs and rock cores were correlated to develop a 13-layer gamma-ray stratigraphy and 41-layer lithostratigraphy throughout the fractured sedimentary rock research site. Detailed hydrogeologic framework shows that black carbon-rich laminated mudstones are the most hydraulically conductive. Water-quality and aquifer-test data indicate that groundwater flow is greatest and TCE contamination is highest in the black, carbon- and clay-rich laminated mudstones. Large-scale groundwater flow at the NAWC research site can be modeled as highly anisotropic with the highest component of permeability occurring along bedding planes.  相似文献   

12.
Rock fractures are of great practical importance to petroleum reservoir engineering because they provide pathways for fluid flow, especially in reservoirs with low matrix permeability, where they constitute the primary flow conduits. Understanding the spatial distribution of natural fracture networks is thus key to optimising production. The impact of fracture systems on fluid flow patterns can be predicted using discrete fracture network models, which allow not only the 6 independent components of the second‐rank permeability tensor to be estimated, but also the 21 independent components of the fully anisotropic fourth‐rank elastic stiffness tensor, from which the elastic and seismic properties of the fractured rock medium can be predicted. As they are stochastically generated, discrete fracture network realisations are inherently non‐unique. It is thus important to constrain their construction, so as to reduce their range of variability and, hence, the uncertainty of fractured rock properties derived from them. This paper presents the underlying theory and implementation of a method for constructing a geologically realistic discrete fracture network, constrained by seismic amplitude variation with offset and azimuth data. Several different formulations are described, depending on the type of seismic data and prior geologic information available, and the relative strengths and weaknesses of each approach are compared. Potential applications of the method are numerous, including the prediction of fluid flow, elastic and seismic properties of fractured reservoirs, model‐based inversion of seismic amplitude variation with offset and azimuth data, and the optimal placement and orientation of infill wells to maximise production.  相似文献   

13.
The potential for trichloroethene (TCE) biodegradation in a fractured dolomite aquifer at a former chemical disposal site in Smithville, Ontario, Canada, is assessed using chemical analysis and TCE and cis‐DCE compound‐specific isotope analysis of carbon and chlorine collected over a 16‐month period. Groundwater redox conditions change from suboxic to much more reducing environments within and around the plume, indicating that oxidation of organic contaminants and degradation products is occurring at the study site. TCE and cis‐DCE were observed in 13 of 14 wells sampled. VC, ethene, and/or ethane were also observed in ten wells, indicating that partial/full dechlorination has occurred. Chlorine isotopic values (δ37Cl) range between 1.39 to 4.69‰ SMOC for TCE, and 3.57 to 13.86‰ SMOC for cis‐DCE. Carbon isotopic values range between ?28.9 and ?20.7‰ VPDB for TCE, and ?26.5 and ?11.8‰ VPDB for cis‐DCE. In most wells, isotopic values remained steady over the 15‐month study. Isotopic enrichment from TCE to cis‐DCE varied between 0 and 13‰ for carbon and 1 and 4‰ for chlorine. Calculated chlorine‐carbon isotopic enrichment ratios (?Cl/?C) were 0.18 for TCE and 0.69 for cis‐DCE. Combined, isotopic and chemical data indicate very little dechlorination is occurring near the source zone, but suggest bacterially mediated degradation is occurring closer to the edges of the plume.  相似文献   

14.
The horizontal reactive media treatment well (HRX Well®) uses directionally drilled horizontal wells filled with a treatment media to induce flow-focusing behavior created by the well-to-aquifer permeability contrast to passively capture proportionally large volumes of groundwater. Groundwater is treated in situ as it flows through the HRX Well and downgradient portions of the aquifer are cleaned via elution as these zones are flushed with clean water discharging from the HRX Well. The HRX Well concept is particularly well suited for sites where long-term mass discharge control is a primary performance objective. This concept is appropriate for recalcitrant and difficult-to-treat constituents, including chlorinated solvents, per- and polyfluoroalkyl substances (PFAS), 1,4-dioxane, and metals. A full-scale HRX Well was installed and operated to treat trichloroethene (TCE) with zero valent iron (ZVI). The model-predicted enhanced flow through the HRX Well (compared to the flow in and equivalent cross-sectional area orthogonal to flow in the natural formation before HRX Well installation) and treatment zone width was consistent with flows and widths estimated independently by point velocity probe (PVP) testing, HRX Well tracer testing, and observed treatment in downgradient monitoring wells. The actual average capture zone width was estimated to be between 45 and 69 feet. Total TCE mass discharge reduction was maintained through the duration of the performance monitoring period and exceeded 99.99% (%). Decreases in TCE concentrations were observed at all four downgradient monitoring wells within the treatment zone (ranging from 50 to 74% at day 436), and the first arrival of treated water was consistent with model predictions. The field demonstration confirmed the HRX Well technology is best suited for long-term mass discharge control, can be installed under active infrastructure, requires limited ongoing operation and maintenance, and has low life cycle energy and water requirements.  相似文献   

15.
This paper presents the largest In Situ Thermal Desorption (ISTD) project completed to date. The redevelopment of a former aerospace manufacturing facility adjacent to a commercial airport was the main driver, requiring relatively rapid reduction of several chlorinated volatile organic compounds (CVOC) in a 3.2‐acre source zone. The source zone was divided into four quadrants with differing treatment depths, heated simultaneously using a total of 907 thermal conduction heater wells. Five different depths were selected across the area, according to the depth of contaminant impact. Prior to implementation, a risk and optimization study led to placement of a vertical sheet‐pile wall around the treatment zone to minimize groundwater flow, and a pilot test of a novel direct‐drive method for installation of the heater casings. Because of a shallow water table, a layer of clean fill was placed over the treatment zone, and partial dewatering was necessary prior to heating. A network of vertical multiphase extraction wells and horizontal vapor extraction wells was used to establish hydraulic and pneumatic control and to capture the contaminants. The site was split into four decision units, each with a rigorous soil sampling program which included collecting a total of 270 confirmatory soil samples from locations with the highest pretreatment CVOC concentrations requiring reduction to below 1 mg/kg for each contaminant. Temperature monitoring and mass removal trends were used to trigger the sampling events. Eventually, a small area near the center of the site required the installation of four additional heaters before the soil goals were reached after 238 days of heating. The total energy usage for heating and treating the source area was 23 million kWh—slightly lower than the estimated 26.5 million kWh. Actual energy losses and the energy removal associated with the extracted steam were lower than anticipated. An estimated 13,400 kg (29,800 lbs) of CVOC mass was removed, and all soil goals were met. This paper presents the challenges associated with a project of this scale and describes the solutions to successfully complete the ISTD remedy.  相似文献   

16.
A new lumped-parameter approach to simulating unsaturated flow processes in dual-porosity media such as fractured rocks or aggregated soils is presented. Fluid flow between the fracture network and the matrix blocks is described by a non-linear equation that relates the imbibition rate of the local difference in liquid-phase pressure between the fractures and the matrix blocks. Unlike a Warren-Root-type equation, this equation is accurate in both the early and late time regimes. The fracture/matrix interflow equation has been incorporated into an existing unsaturated flow simulator, to serve as a source/sink term for fracture gridblocks. Flow processes are then simulated using only fracture gridblocks in the computational grid. This new lumped-parameter approach has been tested on two problems involving transient flow in fractured/porous media, and compared with simulations performed using explicit discretisation of the matrix blocks. The new procedure seems to accurately simulate flow processes in unsaturated fractured rocks, and typically requires an order of magnitude less computational time than do simulations using fully-discretised matrix blocks.  相似文献   

17.
Remediation of groundwater contaminated by chlorinated hydrocarbons via in situ technologies such as direct injection of nanoscale zero valent iron (ZVI, Fe(O)) particles is increasingly common. However, assessing target compound degradation by abiotic processes is difficult because (1) the injection may displace the contaminant plume so that concentration measurements alone are often inconclusive and (2) biodegradation may also occur, making it challenging to identify and evaluate the abiotic degradation component. In this study, trichloroethylene (TCE) and 1,1,1-trichloroethane (1,1,1-TCA) were treated in a highly heterogeneous hydrogeologic setting. The purpose of this study was to evaluate the potential for compound-specific stable isotope analysis (CSIA) to monitor the effectiveness of ZVI injection by assessing TCE and 1,1,1-TCA degradation. Prior to ZVI injection, carbon isotope measurements demonstrated biodegradation of TCE by native microorganisms. This in situ biodegradation was quantified by measuring the enrichment of 13C in TCE samples downstream of the suspected source. When ZVI was injected through only two injection wells, no changes in TCE and 1,1,1-TCA isotope signatures were detected compared to preinjection values. In contrast, when ZVI was injected through 11 wells covering a greater portion of the contaminated area, 5 out of 10 monitoring wells showed further enrichment of 13C in either TCE or 1,1,1-TCA, indicating additional target compound transformation. The abiotic nature of this TCE transformation was confirmed through temporal trends in carbon isotope values of the putative transformation products cis-dichloroethylene (cis-DCE), ethene and ethane. This demonstrates the usefulness of CSIA in distinguishing abiotic vs. biotic transformation in the field.  相似文献   

18.
An equivalent medium model for wave simulation in fractured porous rocks   总被引：3，自引：0，他引：3  
Seismic wave propagation in reservoir rocks is often strongly affected by fractures and micropores. Elastic properties of fractured reservoirs are studied using a fractured porous rock model, in which fractures are considered to be embedded in a homogeneous porous background. The paper presents an equivalent media model for fractured porous rocks. Fractures are described in a stress‐strain relationship in terms of fracture‐induced anisotropy. The equations of poroelasticity are used to describe the background porous matrix and the contents of the fractures are inserted into a matrix. Based on the fractured equivalent‐medium theory and Biot's equations of poroelasticity, two sets of porosity are considered in a constitutive equation. The porous matrix permeability and fracture permeability are analysed by using the continuum media seepage theory in equations of motion. We then design a fractured porous equivalent medium and derive the modified effective constants for low‐frequency elastic constants due to the presence of fractures. The expressions of elastic constants are concise and are directly related to the properties of the main porous matrix, the inserted fractures and the pore fluid. The phase velocity and attenuation of the fractured porous equivalent media are investigated based on this model. Numerical simulations are performed. We show that the fractures and pores strongly influence wave propagation, induce anisotropy and cause poroelastic behaviour in the wavefields. We observe that the presence of fractures gives rise to changes in phase velocity and attenuation, especially for the slow P‐wave in the direction parallel to the fracture plane.  相似文献   

19.
We present a method to determine equivalent permeability of fractured porous media. Inspired by the previous flow-based upscaling methods, we use a multi-boundary integration approach to compute flow rates within fractures. We apply a recently developed multi-point flux approximation Finite Volume method for discrete fracture model simulation. The method is verified by upscaling an arbitrarily oriented fracture which is crossing a Cartesian grid. We demonstrate the method by applying it to a long fracture, a fracture network and the fracture network with different matrix permeabilities. The equivalent permeability tensors of a long fracture crossing Cartesian grids are symmetric, and have identical values. The application to the fracture network case with increasing matrix permeabilities shows that the matrix permeability influences more the diagonal terms of the equivalent permeability tensor than the off-diagonal terms, but the off-diagonal terms remain important to correctly assess the flow field.  相似文献   

20.
A three-dimensional two-phase flow model is coupled to a non-linear reactive transport model to study the efficacy of potassium permanganate treatment on dense, non-aqueous phase liquid (DNAPL) source removal in porous media. A linear relationship between the soil permeability (k) and concentration of manganese dioxide precipitate ([MnO2(s)]), k = ko + Srind [MnO2(s)], is utilized to simulate nodal permeability reductions due to precipitate formation. Using published experimental column studies, an Srind = −5.5 × 10−16 m2 L/mg was determined for trichloroethylene (TCE) DNAPL. This Srind was then applied to treatment simulations on three-dimensional TCE DNAPL source zones comprising either DNAPL at residual saturations, or DNAPL at pooled saturations.  相似文献   

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

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