Discretizing the fracture-matrix interface to simulate solute transport   总被引：1，自引：0，他引：1  

Weatherill D  Graf T  Simmons CT  Cook PG  Therrien R  Reynolds DA 《Ground water》2008,46(4):606-615

This article examines the required spatial discretization perpendicular to the fracture-matrix interface (FMI) for numerical simulation of solute transport in discretely fractured porous media. The discrete-fracture, finite-element model HydroGeoSphere ( Therrien et al. 2005 ) and a discrete-fracture implementation of MT3DMS ( Zheng 1990 ) were used to model solute transport in a single fracture, and the results were compared to the analytical solution of Tang et al. (1981) . To match analytical results on the relatively short timescales simulated in this study, very fine grid spacing perpendicular to the FMI of the scale of the fracture aperture is necessary if advection and/or dispersion in the fracture is high compared to diffusion in the matrix. The requirement of such extremely fine spatial discretization has not been previously reported in the literature. In cases of high matrix diffusion, matching the analytical results is achieved with larger grid spacing at the FMI. Cases where matrix diffusion is lower can employ a larger grid multiplier moving away from the FMI. The very fine spatial discretization identified in this study for cases of low matrix diffusion may limit the applicability of numerical discrete-fracture models in such cases.  相似文献   

The value of subsidence data in ground water model calibration   总被引：2，自引：0，他引：2  

Yan T  Burbey TJ 《Ground water》2008,46(4):538-550

The accurate estimation of aquifer parameters such as transmissivity and specific storage is often an important objective during a ground water modeling investigation or aquifer resource evaluation. Parameter estimation is often accomplished with changes in hydraulic head data as the key and most abundant type of observation. The availability and accessibility of global positioning system and interferometric synthetic aperture radar data in heavily pumped alluvial basins can provide important subsidence observations that can greatly aid parameter estimation. The aim of this investigation is to evaluate the value of spatial and temporal subsidence data for automatically estimating parameters with and without observation error using UCODE-2005 and MODFLOW-2000. A synthetic conceptual model (24 separate cases) containing seven transmissivity zones and three zones each for elastic and inelastic skeletal specific storage was used to simulate subsidence and drawdown in an aquifer with variably thick interbeds with delayed drainage. Five pumping wells of variable rates were used to stress the system for up to 15 years. Calibration results indicate that (1) the inverse of the square of the observation values is a reasonable way to weight the observations, (2) spatially abundant subsidence data typically produce superior parameter estimates under constant pumping even with observation error, (3) only a small number of subsidence observations are required to achieve accurate parameter estimates, and (4) for seasonal pumping, accurate parameter estimates for elastic skeletal specific storage values are largely dependent on the quantity of temporal observational data and less on the quantity of available spatial data.  相似文献   

Assessment of chemical dispersant effectiveness in a wave tank under regular non-breaking and breaking wave conditions     

Li Z  Lee K  King T  Boufadel MC  Venosa AD 《Marine pollution bulletin》2008,56(5):903-912

Current chemical dispersant effectiveness tests for product selection are commonly performed with bench-scale testing apparatus. However, for the assessment of oil dispersant effectiveness under real sea state conditions, test protocols are required to have hydrodynamic conditions closer to the natural environment, including transport and dilution effects. To achieve this goal, Fisheries and Oceans Canada and the US Environmental Protection Agency (EPA) designed and constructed a wave tank system to study chemical dispersant effectiveness under controlled mixing energy conditions (regular non-breaking, spilling breaking, and plunging breaking waves). Quantification of oil dispersant effectiveness was based on observed changes in dispersed oil concentrations and oil-droplet size distribution. The study results quantitatively demonstrated that total dispersed oil concentration and breakup kinetics of oil droplets in the water column were strongly dependent on the presence of chemical dispersants and the influence of breaking waves. These data on the effectiveness of dispersants as a function of sea state will have significant implications in the drafting of future operational guidelines for dispersant use at sea.  相似文献   

Characterization of Damage in Sandstones along the Mojave Section of the San Andreas Fault: Implications for the Shallow Extent of Damage Generation     

Ory Dor  Judith S. Chester  Yehuda Ben-Zion  James N. Brune  Thomas K. Rockwell 《Pure and Applied Geophysics》2009,166(10-11):1747-1773

Following theoretical calculations that suggest shallow generation of rock damage during an earthquake rupture, we measure the degree of fracture damage in young sedimentary rocks from the Juniper Hills Formation (JHF) that were displaced 21 km along the Mojave section of San Andreas Fault (SAF) and were not exhumed significantly during their displacement. In exposures adjacent to the fault, the JHF typically displays original sedimentary fabrics and little evidence of bulk shear strain at the mesoscopic scale. The formation is, however, pervasively fractured at the microscopic scale over a zone that is about a 100 m wide on the southwest side of the SAF near Little Rock. The abundance of open fractures, the poor consolidation, and the shallow inferred burial depth imply that the damage was generated close to the surface of the Earth. The spatial correlation of this damage with a seismically active trace of the SAF suggests that it was generated by SAF slip events that by assumption were of a seismic nature throughout the displacement history of the JHF. Thus the JHF provides a very shallow upper bound for the generation of brittle damage in a seismic fault zone. The fracture fabric is characterized by preferred orientations of fractures that split grains between contact points and is consistent with overall deformation under directed compression. However, the available results cannot be used to distinguish between proposed off-fault damage mechanisms. Fracture orientations are compatible with a maximum compressive stress oriented at a high angle to the fault at about 10 m, and at a lower, more variable angle farther away from the fault. The fracture distribution and fabric are consistent with observations made of the microscale damage characteristics of the Hungry Valley Formation in the northwestern section of the SAF in the Mojave, and with previous observations of exhumed, ancestral strands of the SAF.  相似文献   

Kinematic transformations for planar multi‐directional pseudodynamic testing     

Oya Mercan  James M. Ricles  Richard Sause  Thomas Marullo 《地震工程与结构动力学》2009,38(9):1093-1119

The pseudodynamic (PSD) test method imposes command displacements to a test structure for a given time step. The measured restoring forces and displaced position achieved in the test structure are then used to integrate the equations of motion to determine the command displacements for the next time step. Multi‐directional displacements of the test structure can introduce error in the measured restoring forces and displaced position. The subsequently determined command displacements will not be correct unless the effects of the multi‐directional displacements are considered. This paper presents two approaches for correcting kinematic errors in planar multi‐directional PSD testing, where the test structure is loaded through a rigid loading block. The first approach, referred to as the incremental kinematic transformation method, employs linear displacement transformations within each time step. The second method, referred to as the total kinematic transformation method, is based on accurate nonlinear displacement transformations. Using three displacement sensors and the trigonometric law of cosines, this second method enables the simultaneous nonlinear equations that express the motion of the loading block to be solved without using iteration. The formulation and example applications for each method are given. Results from numerical simulations and laboratory experiments show that the total transformation method maintains accuracy, while the incremental transformation method may accumulate error if the incremental rotation of the loading block is not small over the time step. A procedure for estimating the incremental error in the incremental kinematic transformation method is presented as a means to predict and possibly control the error. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

MODFLOW-CFP: A New Conduit Flow Process for MODFLOW–2005     

reviewed by Thomas Reimann  and Melissa E. Hill 《Ground water》2009,47(3):321-325

  相似文献   

Time–space continuity of daily maps of fractional snow cover and albedo from MODIS     

Jeff Dozier  Thomas H. Painter  Karl Rittger  James E. Frew 《Advances in water resources》2008

Using reflectance values from the seven MODIS “land” bands with 250 or 500 m resolution, along with a corresponding cloud product, we estimate the fraction of each 500 m pixel that snow covers, along with the albedo of that snow. The daily products have data gaps and errors because of cloud cover and sensor viewing geometry. Rather than make users interpolate and filter these patchy daily maps without completely understanding the retrieval algorithm and instrument properties, we use the daily time series to improve the estimate of the measured snow properties for a particular day. We use a combination of noise filtering, snow/cloud discrimination, and interpolation and smoothing to produce our best estimate of the daily snow cover and albedo. We consider two modes: one is the “predictive” mode, whereby we estimate the snow-covered area and albedo on that day using only the data up to that day; the other is the “retrospective” mode, whereby we reconstruct the history of the snow properties for a previous period.  相似文献   

Hydrodynamics and sediment dynamics in the German Bight. A focus on observations and numerical modelling in the East Frisian Wadden Sea     

Joanna Staneva  Emil V. Stanev  Jörg-Olaf Wolff  Thomas H. Badewien  Rainer Reuter  Burghard Flemming  Alexander Bartholomä  Karsten Bolding 《Continental Shelf Research》2009

This work deals with analysis of hydrographic observations and results of numerical simulations. The data base includes acoustic Doppler current profilers (ADCP) observations, continuous measurements on data stations and satellite data originating from the medium resolution imaging spectrometer (MERIS) onboard the European Space Agency (ESA) satellite ENVISAT with a spatial resolution of 300 m. Numerical simulations use nested models with horizontal resolutions ranging from 1 km in the German Bight to 200 m in the East Frisian Wadden Sea coupled with a suspended matter transport model. Modern satellite observations have now a comparable horizontal resolution with high-resolution numerical model of the entire area of the East Frisian Wadden Sea allowing to describe and validate new and so far unknown patterns of sediment distribution. The two data sets are consistent and reveal an oscillatory behaviour of sediment pools to the north of the back-barrier basins and clear propagation patterns of tidally driven suspended particulate matter outflow into the North Sea. The good agreement between observations and simulations is convincing evidence that the model simulates the basic dynamics and sediment transport processes, which motivates its further use in hindcasting, as well as in the initial steps towards forecasting circulation and sediment dynamics in the coastal zone.  相似文献   

The November 2002 eruption of Piton de la Fournaise,Réunion: tracking the pre-eruptive thermal evolution of magma using melt inclusions     

Nathalie Vigouroux  A. E. Williams-Jones  Paul Wallace  Thomas Staudacher 《Bulletin of Volcanology》2009,71(9):1077-1089

The November 2002 eruption of Piton de la Fournaise in the Indian Ocean was typical of the activity of the volcano from 1999 to 2006 in terms of duration and volume of magma ejected. The first magma erupted was a basaltic liquid with a small proportion of olivine phenocrysts (Fo₈₁) that contain small numbers of melt inclusions. In subsequent flows, olivine crystals were more abundant and richer in Mg (Fo_83–84). These crystals contain numerous melt and fluid inclusions, healed fractures, and dislocation features such as kink bands. The major element composition of melt inclusions in this later olivine (Fo_83–84) is out of equilibrium with that of its host as a result of extensive post-entrapment crystallization and Fe²⁺ loss by diffusion during cooling. Melt inclusions in Fo₈₁ olivine are also chemically out of equilibrium with their hosts but to a lesser degree. Using olivine–melt geothermometry, we determined that melt inclusions in Fo₈₁ olivine were trapped at lower temperature (1,182 ± 1°C) than inclusions in Fo_83–84 olivine (1,199–1,227°C). This methodology was also used to estimate eruption temperatures. The November 2002 melt inclusion compositions suggest that they were at temperatures between 1,070°C and 1,133°C immediately before eruption and quenching. This relatively wide temperature range may reflect the fact that most of the melt inclusions were from olivine in lava samples and therefore likely underwent minor but variable amounts of post-eruptive crystallization and Fe²⁺ loss by diffusion due to their relatively slow cooling on the surface. In contrast, melt inclusions in tephra samples from past major eruptions yielded a narrower range of higher eruption temperatures (1,163–1,181°C). The melt inclusion data presented here and in earlier publications are consistent with a model of magma recharge from depth during major eruptions, followed by storage, cooling, and crystallization at shallow levels prior to expulsion during events similar in magnitude to the relatively small November 2002 eruption.  相似文献   

Structural features of Panarea volcano in the frame of the Aeolian Arc (Italy): Implications for the 2002–2003 unrest     

Valerio Acocella  Marco Neri  Thomas R. Walter   《Journal of Geodynamics》2009,47(5):288-292

Panarea, characterized by gas unrest in 2002–2003, is the volcanic island with the least constrained structure in the eastern-central Aeolian Arc (Italy). Based on structural measurements, we define here its deformation pattern relative to the Arc. The main deformations are subvertical extension fractures (63% of data), normal faults (25%) and dikes (12%). The mean orientation of the extension fractures and faults is N38°E, with a mean opening direction of N135° ± 8°, implying extension with a moderate component of dextral shear. These data, matched with those available for Stromboli volcano (pure opening) and Vulcano, Lipari and Salina volcanoes (predominant dextral motions) along the eastern-central Arc, suggest a progressive westward rotation of the extension direction and an increase in the dextral shear. The dextral shear turns into compression in the western arc. The recent unrest at Panarea, coeval to that of nearby Stromboli, may also be explained by the structural context, as both volcanoes lie along the portion of the Arc subject to extension.  相似文献