首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
A key challenge for predictive modeling of transverse mixing and reaction of solutes in groundwater is to determine values of transverse dispersivity (αT)(αT) in heterogeneous flow fields that accurately describe mixing and reaction at the pore scale. We evaluated the effects of flow focusing in high permeability zones on mixing enhancement using experimental micromodel flow cells and pore-scale lattice-Boltzmann-finite-volume model (LB-FVM) simulations. Micromodel results were directly compared to LB-FVM simulations using two different pore structures, and excellent agreement was obtained. Six different flow focusing pore structures were then systematically tested using LB-FVM, and both analytical solutions and a two-dimensional (2D) continuum-scale model were used to fit αTαT values to pore-scale results. Pore-scale results indicate that the overall rate of mixing-limited reaction increased by up to 40% when flow focusing occurred, and it was greater in pore structures with longer flow focusing regions and greater porosity contrast. For each pore structure, αTαT values from analytical solutions of transverse concentration profiles or total product at a given longitudinal location showed good agreement for nonreactive and reactive solutes, and values determined in flow focusing zones were always smaller than those downgradient after the flow focusing zone. Transverse dispersivity values from the 2D continuum model were between values within and downgradient from the flow focusing zone determined from analytical solutions. Also, total product and transverse concentration profiles along the entire pore structure from the 2D continuum model matched pore scale results. These results indicate that accurate quantification of pore-scale flow focusing with transverse dispersion coefficients is possible only when the entire flow and concentration fields are considered.  相似文献   

2.
Geochemical reaction rate laws are often measured using crushed minerals in well-mixed laboratory systems that are designed to eliminate mass transport limitations. Such rate laws are often used directly in reactive transport models to predict the reaction and transport of chemical species in consolidated porous media found in subsurface environments. Due to the inherent heterogeneities of porous media, such use of lab-measured rate laws may introduce errors, leading to a need to develop methods for upscaling reaction rates. In this work, we present a methodology for using pore-scale network modeling to investigate scaling effects in geochemical reaction rates. The reactive transport processes are simulated at the pore scale, accounting for heterogeneities of both physical and mineral properties. Mass balance principles are then used to calculate reaction rates at the continuum scale. To examine the scaling behavior of reaction kinetics, these continuum-scale rates from the network model are compared to the rates calculated by directly using laboratory-measured reaction rate laws and ignoring pore-scale heterogeneities. In this work, this methodology is demonstrated by upscaling anorthite and kaolinite reaction rates under simulation conditions relevant to geological CO2 sequestration. Simulation results show that under conditions with CO2 present at high concentrations, pore-scale concentrations of reactive species and reaction rates vary spatially by orders of magnitude, and the scaling effect is significant. With a much smaller CO2 concentration, the scaling effect is relatively small. These results indicate that the increased acidity associated with geological sequestration can generate conditions for which proper scaling tools are yet to be developed. This work demonstrates the use of pore-scale network modeling as a valuable research tool for examining upscaling of geochemical kinetics. The pore-scale model allows the effects of pore-scale heterogeneities to be integrated into system behavior at multiple scales, thereby identifying important factors that contribute to the scaling effect.  相似文献   

3.
We consider 3D steady flow of fresh water over a salt water body in a confined aquifer of constant thickness D, with application to a pumping well in a coastal aquifer. With neglect of mixing, a sharp interface separates the two fluid bodies and an existing analytical solution, based on the Dupuit assumption, is adopted. The aim is to solve for the mixing between the fresh and salt waters for αT/D  1 (αT transverse dispersivity), as field studies indicate that αT = O(10−3 − 10−2 m). The mixing zone around the interface is narrow and solutions by existing codes experience numerical difficulties. The problem is solved by the boundary layer (BL) approximation, extending a method, applied previously to two-dimensional flows. The BL equations of variable-density flow are solved by using the Von Karman integral method, to determine the BL thickness and the rate of entrainment of salt water along the interface. Application to the pumping well problem yields the salinity of the pumped water, as function of the parameters of the problem (well discharge, seaward discharge, well distance from the coast and density difference).  相似文献   

4.
Thermal diffusivity (D) was measured using laser-flash analysis on pristine and remelted obsidian samples from Mono Craters, California. These high-silica rhyolites contain between 0.013 and 1.10?wt% H2O and 0 to 2?vol% crystallites. At room temperature, D glass varies from 0.63 to 0.68?mm2?s?1, with more crystalline samples having higher D. As T increases, D glass decreases, approaching a constant value of ??0.55?mm2?s?1 near 700?K. The glass data are fit with a simple model as an exponential function of temperature and a linear function of crystallinity. Dissolved water contents up to 1.1?wt% have no statistically significant effect on the thermal diffusivity of the glass. Upon crossing the glass transition, D decreases rapidly near ??1,000?K for the hydrous melts and ??1,200?K for anhydrous melts. Rhyolitic melts have a D melt of ??0.51?mm2?s?1. Thermal conductivity (k?=?D·??·C P) of rhyolitic glass and melt increases slightly with T because heat capacity (C P) increases with T more strongly than density (??) and D decrease. The thermal conductivity of rhyolitic melts is ??1.5?W?m?1?K?1, and should vary little over the likely range of magmatic temperatures and water contents. These values of D and k are similar to those of major crustal rock types and granitic protoliths at magmatic temperatures, suggesting that changes in thermal properties accompanying partial melting of the crust should be relatively minor. Numerical models of shallow rhyolite intrusions indicate that the key difference in thermal history between bodies that quench to obsidian, and those that crystallize, results from the release of latent heat of crystallization. Latent heat release enables bodies that crystallize to remain at high temperatures for much longer times and cool more slowly than glassy bodies. The time to solidification is similar in both cases, however, because solidification requires cooling through the glass transition in the first case, and cooling only to the solidus in the second.  相似文献   

5.
A number of recent optically stimulated luminescence (OSL) studies have cited post-depositional mixing as a dominant source of equivalent dose (De) scatter across a range of sedimentary environments, including those previously considered ‘best suited’ for OSL dating. The potentially insidious nature of sediment mixing means that this problem may often only be identifiable by careful statistical analysis of De data sets. This study aims to address some of the important issues associated with the characterisation and statistical treatment of mixed De distributions at the multi-grain scale of analysis, using simulated De data sets produced with a simple stochastic model. Using this Monte Carlo approach we were able to generate theoretical distributions of single-grain De values, which were then randomly mixed together to simulate multi-grain aliquot De distributions containing a known number of mixing components and known corresponding burial doses. A range of sensitivity tests were undertaken using sediment mixtures with different aged dose components, different numbers of mixing components, and different types of dose component distributions (fully bleached, heterogeneously bleached and significantly overdispersed De distributions). The results of our modelling simulations reveal the inherent problems encountered when dating mixed sedimentary samples with multi-grain De estimation techniques. ‘Phantom’ dose components (i.e. discrete dose populations that do not correspond to the original single-grain mixing components) are an inevitable consequence of the ‘averaging’ effects of multi-grain De analysis, and prevent the correct number of mixing components being identified with the finite mixture model (FMM) for all of the multi-grain mixtures tested. Our findings caution against use of the FMM for multi-grain aliquot De data sets, even when the aliquots consist of only a few grains.  相似文献   

6.
  相似文献   

7.
Upscaling pore-scale processes into macroscopic quantities such as hydrodynamic dispersion is still not a straightforward matter for porous media with complex pore space geometries. Recently it has become possible to obtain very realistic 3D geometries for the pore system of real rocks using either numerical reconstruction or micro-CT measurements. In this work, we present a finite element–finite volume simulation method for modeling single-phase fluid flow and solute transport in experimentally obtained 3D pore geometries. Algebraic multigrid techniques and parallelization allow us to solve the Stokes and advection–diffusion equations on large meshes with several millions of elements. We apply this method in a proof-of-concept study of a digitized Fontainebleau sandstone sample. We use the calculated velocity to simulate pore-scale solute transport and diffusion. From this, we are able to calculate the a priori emergent macroscopic hydrodynamic dispersion coefficient of the porous medium for a given molecular diffusion Dm of the solute species. By performing this calculation at a range of flow rates, we can correctly predict all of the observed flow regimes from diffusion dominated to convection dominated.  相似文献   

8.
Two end-member types of pyroclastic density current are commonly recognized: pyroclastic surges are dilute currents in which particles are carried in turbulent suspension and pyroclastic flows are highly concentrated flows. We provide scaling relations that unify these end-members and derive a segregation mechanism into basal concentrated flow and overriding dilute cloud based on the Stokes number (ST), the stability factor (ΣT) and the dense-dilute condition (DD). We recognize five types of particle behaviors within a fluid eddy as a function of ST and ΣT: (1) particles sediment from the eddy, (2) particles are preferentially settled out during the downward motion of the eddy, but can be carried during its upward motion, (3) particles concentrate on the periphery of the eddy, (4) particles settling can be delayed or ‘fast-tracked’ as a function of the eddy spatial distribution, and (5) particles remain homogeneously distributed within the eddy. We extend these concepts to a fully turbulent flow by using a prototype of kinetic energy distribution within a full eddy spectrum and demonstrate that the presence of different particle sizes leads to the density stratification of the current. This stratification may favor particle interactions in the basal part of the flow and DD determines whether the flow is dense or dilute. Using only intrinsic characteristics of the current, our model explains the discontinuous features between pyroclastic flows and surges while conserving the concept of a continuous spectrum of density currents.  相似文献   

9.
Field tracer tests indicate that advective porosity, the quantity relating advective velocity to Darcy flux, may exhibit directional dependence. Hydraulic anisotropy explains some but not all of the reported directional results. The present paper shows mathematically that directional variations in advective porosity may arise simply from incomplete mixing of an inert tracer between directional flow channels within a sampling (or support) volume ω of soil or rock that may be hydraulically isotropic or anisotropic. In the traditional fully homogenized case, our theory yields trivially a scalar advective porosity equal to the interconnected porosity ϕ, thus explaining neither the observed directional effects nor the widely reported experimental finding that advective porosity is generally smaller than ϕ. We consider incomplete mixing under conditions in which the characteristic time tD of longitudinal diffusion along channels across ω is much shorter than the characteristic time tH required for homogenization through transverse diffusion between channels. This may happen where flow takes place preferentially through relatively conductive channels and/or fractures of variable orientation separated by material that forms a partial barrier to diffusive transport. Our solution is valid for arbitrary channel Peclet numbers on a correspondingly wide range of time scales tD ⩽ t ≪ tH. It shows that the tracer center of mass is advected at a macroscopic velocity which is generally not collinear with the macroscopic Darcy flux and exceeds it in magnitude. These two vectors are related through a second-rank symmetric advective dispersivity tensor Φ. If the permeability k of ω is a symmetric positive-definite tensor, so is Φ. However, the principal directions and values of these two tensors are generally not the same; whereas those of k are a fixed property of the medium and the length-scale of ω, those of Φ depend additionally on the direction and magnitude of the applied hydraulic gradient. When the latter is large, diffusion has negligible effect on Φ and one may consider tracer mass to be distributed between channels in proportion to the magnitude of their Darcy flux. This is made intuitive through a simple example of an idealized fracture network. Our analytical formalism reveals the properties of Φ but is too schematic to allow predicting the latter accurately on the basis of realistic details about the void structure of ω and tracer mass distribution within it. Yet knowing the tensorial properties of Φ is sufficient to allow determining it indirectly on the basis of ω-scale hydraulic and tracer data, including concentrations that represent homogenized samples extracted from (or sensed externally across) an ω-scale plume.  相似文献   

10.
A spectral analysis of simultaneous diurnal variations in the E z component of the quasi-static electric field in the near-Earth atmosphere, VLF radio noise, and the horizontal component of the geomagnetic field, observed at Kamchatka in September 1999, has been performed. These geophysical parameters are indirectly used to study wave processes in the near-Earth atmosphere and in the ionospheric D and dynamo regions within the band of periods of internal gravity waves (T = 0.5?3.5 h). The correlation method in the frequency region is used to analyze the interrelation between the wave processes in these atmospheric regions. The power cross-spectra of various pairs of geophysical parameters have been studied depending on meteorological, seismic, and geomagnetic activities. It is shown that the oscillations in the power spectra in the T ～ 1–1.5 h band of periods are caused by the sources of internal gravity waves in the near-Earth atmosphere and by the remote sources above the dynamo region of the ionosphere within the T ～ 1.5–3 h band of periods.  相似文献   

11.
In this study, seismic data recorded during the period 01/01/1996 to 09/01/2009 has been used to evaluate the seismic hazard potential along the Alborz region, Northern Iran. The technique of mapping local recurrence time, T L, is used to map major asperities, which are considered as the areas with maximum hazard. We calculated T L from a and b values which are in turn derived from the frequency–magnitude relation constants within a radius of 30 km about every corner point of a 10-km spacing grid. Since b value is inversely related to applied stress, the areas with lowest b values and/or shortest T L are interpreted to locate the asperities or the areas of maximum seismic hazard. To test this method, we computed T L map using seismic catalogues before and after the 2004 Baladeh earthquake of M w 6.2. The local recurrence time map before the earthquake shows anomalously short T L in the epicentral region of the Baladeh earthquake a decade before its occurrence. The T L map after the earthquake indicates that this large event has redistributed the applied stress in the Alborz region. The microseismicity of the region after the Baladeh earthquake, however, suggests that there are two anomalies in T L map positioned in Alborz. The places where these anomalies are observed can be considered as the areas with maximum seismic hazard for future large earthquake in the Alborz region.  相似文献   

12.
The role of water depth and bottom boundary layer turbulence upon lee-wave generation in sill regions is examined. Their effect upon vertical mixing is also considered. Calculations are performed using a non-hydrostatic model in cross-section form with a specified tidal forcing. Initial calculations in deeper water and a sill height such that the sill top is well removed from the surrounding bed region showed that downstream lee-wave generation and associated mixing increased as bottom friction coefficient k increased. This was associated with an increase in current shear across the sill. However, for a given k, increasing vertical eddy viscosity A v reduced vertical shear in the across sill velocity, leading to a reduction in lee-wave amplitude and associated mixing. Subsequent calculations using shallower water showed that for a given k and A v, lee-wave generation was reduced due to the shallower water depth and changes in the bottom boundary layer. However, in this case (unlike in the deepwater case), there is an appreciable bottom current. This gives rise to bottom mixing which in shallow water extends to mid-depth and enhances the mid-water mixing that is found on the lee side of the sill. Final calculations with deeper water but small sill height showed that lee waves could propagate over the sill, thereby reducing their contribution to mixing. In this case, bottom mixing was the major source of mixing which was mainly confined to the near bed region, with little mid-water mixing.  相似文献   

13.
To study the amount of heat generated by radioactive decay in the continental crust, the usual practice in the literature is to fit to the heat flow and radioactivity data a relationship of the form: Q = Qr + D · A where Q and A are the observed heat flow and radiogenic heat production. Qr is the “reduced” heat flow and D is a depth scale. This procedure implicitly assumes that uranium, thorium and potassium have identical distributions in the crust. We suggest that significant information may be lost as the three radioelements may in fact be affected by processes operating over different depths.Data published for four heat flow provinces throughout the world are used to estimate the distributions of uranium, thorium and potassium in the continental crust. These distributions are characterized by a depth scales defined as follows: Di =∫0h Ci(z)Ci(0)dz where h is the thickness of the layer containing the bulk of radioactivity and Ci(z) the concentration of element i at depth z. Three depth scales are computed from a least-squares fit to the following relationship: Q = Qr + DU · AU + DT · AT + DK · AT where Q is the observed heat flow and Qr some constant (a reduced heat flow). Ai is the heat generation rate due to the radioactive decay of element i, and Di is the corresponding depth scale.The analysis suggests that the three distributions are different and that they have the same basic features in the four provinces considered. The depth scale for potassium is large in granitic areas, that for thorium is small and that for uranium lies between the other two.We propose a simple model according to which each radioelement essentially provides a record for one process. Potassium gives a depth scale for the primary differentiation of the crust. Thorium gives the depth scale of magmatic or metamorphic fluid circulation. Finally, the uranium distribution reflects the late effects of alteration due to meteoric water. We show that the heat flow and radioactivity data are compatible with this model.Our analysis and numerical results are supported by data from deep boreholes and by geochemical evidence, such as detailed investigations of plutonic series and studies of U-Th-Pb systematics.  相似文献   

14.
One-dimensional nuclear magnetic resonance (1D NMR) logging technology is limited for fluid typing, while two-dimensional nuclear magnetic resonance (2D NMR) logging can provide more parameters including longitudinal relaxation time (T 1) and transverse relaxation time (T 2) relative to fluid types in porous media. Based on the 2D NMR relaxation mechanism in a gradient magnetic field, echo train simulation and 2D NMR inversion are discussed in detail. For 2D NMR inversion, a hybrid inversion method is proposed based on the damping least squares method (LSQR) and an improved truncated singular value decomposition (TSVD) algorithm. A series of spin echoes are first simulated with multiple waiting times (T W s) in a gradient magnetic field for given fluid models and these synthesized echo trains are inverted by the hybrid method. The inversion results are consistent with given models. Moreover, the numerical simulation of various fluid models such as the gas-water, light oil-water, and vicious oil-water models were carried out with different echo spacings (T E s) and T W s by this hybrid method. Finally, the influences of different signal-to-noise ratios (SNRs) on inversion results in various fluid models are studied. The numerical simulations show that the hybrid method and optimized observation parameters are applicable to fluid typing of gas-water and oil-water models.  相似文献   

15.
  相似文献   

16.
The macroscopic spreading and mixing of solute plumes in saturated porous media is ultimately controlled by processes operating at the pore scale. Whilst the conventional picture of pore-scale mechanical dispersion and molecular diffusion leading to persistent hydrodynamic dispersion is well accepted, this paradigm is inherently two-dimensional (2D) in nature and neglects important three-dimensional (3D) phenomena. We discuss how the kinematics of steady 3D flow at the pore scale generate chaotic advection—involving exponential stretching and folding of fluid elements—the mechanisms by which it arises and implications of microscopic chaos for macroscopic dispersion and mixing. Prohibited in steady 2D flow due to topological constraints, these phenomena are ubiquitous due to the topological complexity inherent to all 3D porous media. Consequently 3D porous media flows generate profoundly different fluid deformation and mixing processes to those of 2D flow. The interplay of chaotic advection and broad transit time distributions can be incorporated into a continuous-time random walk (CTRW) framework to predict macroscopic solute mixing and spreading. We show how these results may be generalised to real porous architectures via a CTRW model of fluid deformation, leading to stochastic models of macroscopic dispersion and mixing which both honour the pore-scale kinematics and are directly conditioned on the pore-scale architecture.  相似文献   

17.
Earthquakes began to occur in Koyna region (India) soon after the filling of Koyna Dam in 1962. In the present study, three datasets 1964–1993, 1993–1995, and 1996–1997 are analyzed to study the b-value and fractal dimension. The b-value is calculated using the Gutenberg–Richter relationship and fractal dimension D corr. using correlation integral method. The estimated b-value and D corr. of this region before 1993 are found to be in good agreement with previously reported studies. In the subsequent years after 1995, the b-value shows an increase. The estimated b-values of this region are found within the limits of global average. Also, the pattern of spatial clustering of earthquakes show increase in clustering and migration along the three zones called North-East Zone, South-East Zone (SEZ), and Warna Seismic Zone. The earthquake events having depth ≤5 km are largely confined to SEZ. After 1993, the D corr. shows decrease, implying that earthquake activity gets clustered. This seismic clustering could be helpful for earthquake forecasting.  相似文献   

18.
Abstract

First the exact numerical solutions of DIA system of equations describing the transportation of magnetic field in an infinite medium are presented. It is assumed that the turbulence is stationary, homogeneous, isotropic and incompressible. The spectra of turbulence of δ-type and Kolmogorov's type were used. The steady state values of magnetic field diffusivity DT and the α-effect coefficient α T were calculated for various values of space-scale and lifetimes of these spectra and the spectra of helicity. Also investigated is the dependence of DT and α T on the degree of helicity. The corrections to the α T -coefficient due to the contribution of four-order velocity correlators are given. The results are compared with those due to the self-consistent technique.  相似文献   

19.
We used a wavelet formulation of the classical spectral isostatic analysis to invert satellite-derived gravity and topography/bathymetry for elastic thickness (Te) over South America and its surrounding plates. To provide a homogeneous representation of the gravity field for this vast region, we corrected free-air anomalies derived from a combination of terrestrial/marine gravity data with data from the GRACE and CHAMP satellite missions (model EIGEN-CG03C) by a simple Bouguer slab using a smoothed representation of surface relief (wavelengths > 125 km). The resulting Bouguer anomaly compares well with terrestrial data acquired in the Central Andes and allows Te to be confidently estimated for values greater than 10 km. The Te map resolves regional-scale features that are well-correlated with known surface structures and shows maximum values of 100 ± 15 km over the Archean–Neoproterozoic core of the continent, decreasing to less than 30 km around continental margins. Several regions of the oceanic plates and continental margins have an elastic thickness less than 10 km. We performed a quantitative analysis by comparing the elastic thickness with the thermal structure predicted from the age of oceanic crust and igneous–metamorphic rocks. This demonstrates that oceanic plates have been weakened by thermal interaction with hotspots and locally by fracturing and hydration near the trench. We observe that only the nucleus of the continent has resisted the thermomechanical weakening induced by the rifting of Africa and South America along the passive margin and the Andean orogeny along the active margin. This latter region shows along-strike variations in Te that correlate with the geotectonic segmentation of the margin and with the pattern of crustal seismicity. Our results reveal that the rigidity structure follows the segmentation of the seismogenic zone along the subduction fault, suggesting a causal relationship that should be investigated in order to improve the understanding and predictability of great earthquakes and tsunamis.  相似文献   

20.
  相似文献   

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

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