Dispersive Particle Transport: Identification of Macroscale Behavior in Heterogeneous Stratified Subsurface Flows

Dispersive mass transport processes in naturally heterogeneous geological formations (porous media) are investigated based on a particle approach to mass transport and on its numerical implementation using LPT3D, a Lagrangian Particle Tracking 3D code. We are currently using this approach for studying microscale and macroscale space–time behavior (advection, diffusion, dispersion) of tracer plumes, solutes, or miscible fluids, in 1,2,3-dimensional heterogeneous and anisotropic subsurface formations (aquifers, petroleum reservoirs). Our analyses are based on a general advection-diffusion model and numerical scheme where concentrations and fluxes are discretized in terms of particles. The advection-diffusion theory is presented in a probabilistic framework, and in particular, a numerical analysis is developed for the case of advective transport and rotational flows (numerical stability of the explicit Euler scheme). The remainder of the paper is devoted to the behavior of concentration, mass flux density, and statistical moments of the transported tracer plume in the case of heterogeneous steady flow fields, where macroscale dispersion occurs due to geologic heterogeneity and stratification. We focus on the case of perfectly stratified or multilayered media, obtained by generating many horizontal layers with a purely random transverse distribution of permeability and horizontal velocity. In this case, we calculate explicitly the exact mass concentration field C(x, t), mass flux density field f(x, t), and moments. This includes spatial moments and dispersion variance ² _x (t) on a finite domain L, and temporal moments on a finite time scale T, e.g., the mass variance of arrival times ² _T (x). The moments are related to flux concentrations in a way that takes explicitly into account finite space–time scales of analysis (time-dependent tracer mass; spatially variable flow through mass). The multilayered model problem is then used in numerical experiments for testing different ways of recovering information on tracer plume migration, dispersion, concentration and flux fields. Our analyses rely on a probabilistic interpretation that emerges naturally from the particle approach; it is based on spatial moments (particle positions), temporal moments (mass weighted arrival times), and probability densities (both concentrations and fluxes). Finally, as an alternative to direct estimations of the flux and concentration fields, we formulate and study the Moment Inverse Problem. Solving the MIP yields an indirect method for estimating the space–time distribution of flux concentrations based on observed or estimated moments of the plume. The moments may be estimated from field measurements, or numerically computed by particle tracking as we do here.     

Natural tracers in recent groundwaters from different Alpine aquifers

Groundwater with underground residence times between days and a few years have been investigated over more than 20 years from 487 remote sites located in different aquifer types in the Alpine belt. Analysis of the data reveals that groundwaters evolved in crystalline, evaporite, carbonate, molasse, and flysch aquifers can be clearly distinguished based on their major and trace element composition and degree of mineralisation. A further subdivision can be made even within one aquifer type based on the trace element compositions, which are characteristic for the lithologic environment. Major and trace element concentrations can be quantitatively described by interaction of the groundwater with the aquifer-specific mineralogy along the flow path. Because all investigated sites show minimal anthropogenic influences, the observed concentration ranges represent the natural background concentrations and can thus serve as a geo-reference for recent groundwaters from these five aquifer types. This geo-reference is particularly useful for the identification of groundwater contamination. It further shows that drinking water standards can be grossly exceeded for critical elements by purely natural processes.

---

Resumen Durnate más de 20 años se ha investigado aguas subterráneas con una residencia subterránea con una duración de días a varios años en 487 puntos remotos localizados en diferentes tipos de acuíferos en la cadena alpina. El análisis de los datos revela que las aguas subterráneas que han evolucionado dentro de acuíferos cristalinos, evaporíticos, carbonatos, flysch, y molasse se pueden distignuir claramente en base a la composición de sus elementos mayores y marcadores y al grado de minerlización. Asimismo es posible hacer una subdivisión más específica incluso dentro de un tipo de acuífero en base a las composiciones de los elementos marcadores los cuales soncaracterísticos del ambiente litológico. Las concentraciones de los elementos marcadores se pueden describir cuantitativamente por la interacción de las aguas subterráneas con la mineralogía específica del acuífero a lo largo del trayecto del flujo. Puesto qze todos los puntos investigados muestran mínimas influencias antropogénicas, los rangos deconcentraciones observados representan las concentraciones delescenario natural y, por tanto, pueden servir como georeferencia para aguas subterráneas recientes que forman parte de estos cinco tipos de acuíferos. Esta georeferencia es particularmente útil para la identificcación de contaminación de aguas subterráneas. Asimismo muestra que los estándares de agua potable pueden mostrar excesos en elementos críticos por procesos puramente naturales.

---

Résumé Leau souterraine ayant résidé sous la surface du sol entre quelques jours et quelques années a été étudiée sur une période de plus de 20 ans à partir de 487 sites éloignés situés dans différents types daquifères de la ceinture alpine. Lanalyse des données révèle que leau qui a évoluée dans des aquifères cristallins, évaporitiques, carbonatés, molassiques et composés de flysch, peut être facilement différenciée sur la base de sa composition en éléments majeurs et traces, ainsi que par son degré de minéralisation. Une subdivision supplémentaire peut être apportée à lintérieur même dun type daquifère en se basant sur la composition des éléments traces, lesquels sont caractéristiques de lenvironnement lithologique. La concentration en éléments majeurs et traces peut être expliquée quantitativement par linteraction de leau souterraine avec la composition minéralogique spécifique de laquifère le long des lignes découlement. Puisque la majorité des sites étudiés ne montrent que très peu dinfluences anthropogéniques, le registre de concentrations observées représente la concentration de fond naturelle et peut ainsi servir comme géo-référence pour les eaux souterraines récentes dans ces cinq types daquifères. Ces géo-références sont particulièrement utiles pour lidentification de la contamination des eaux souterraines. Par ailleurs, cela démontre que les standards deau potable peuvent être excédés pour certains éléments en raison de processus purement naturels.

     

温度探测土坝圆柱状集中渗漏模型研究

渗漏是堤坝破坏主要形式之一，只有确定渗漏发生的确切位置和集中渗漏通道的大小和强度，才能对其有效治理。本文基于渗流和渗漏情况下堤坝温度场特征及温度示踪基本原理，对复杂的实际问题进行科学简化合理假设，把集中渗漏通道简化为圆柱状并作为边界条件，结合钻孔探测温度实际，运用稳定热传导理论，建立了圆柱状堤坝渗漏温度示踪模型。利用测井温度曲线的最大异常温度，该模型可以计算渗漏通道的位置，渗漏通道大小，流速等参数。工程实践表明，对于垂向流速较小的地下水渗漏进行温度场示踪是有效的方法。     

The use of laboratory experiments for the study of conservative solute transport in heterogeneous porous media

 Laboratory experiments on heterogeneous porous media (otherwise known as intermediate scale experiments, or ISEs) have been increasingly relied upon by hydrogeologists for the study of saturated and unsaturated groundwater systems. Among the many ongoing applications of ISEs is the study of fluid flow and the transport of conservative solutes in correlated permeability fields. Recent advances in ISE design have provided the capability of creating correlated permeability fields in the laboratory. This capability is important in the application of ISEs for the assessment of recent stochastic theories. In addition, pressure-transducer technology and visualization methods have provided the potential for ISEs to be used in characterizing the spatial distributions of both hydraulic head and local water velocity within correlated permeability fields. Finally, various methods are available for characterizing temporal variations in the spatial distribution (and, thereby, the spatial moments) of solute concentrations within ISEs. It is concluded, therefore, that recent developments in experimental techniques have provided an opportunity to use ISEs as important tools in the continuing study of fluid flow and the transport of conservative solutes in heterogeneous, saturated porous media. Received, December 1996 · Revised, July 1997 · Accepted, August 1997     

地下水含水介质的弥散度测定

以水动力弥散理论为基础,根据改进的水动力弥散装置,以氯化钠作为示踪剂,用电导仪间接地测定溶质的突破曲线,最后用解析法确定地下水多孔介质中溶质迁移的弥散度。     

论铅同位素的地质指示作用

在地球壳、慢物质交换过程中,铅同位素具有从均-趋向不均-或从非均-趋向新的均-的倾向,而使其具备作为指示剂的条件。研究时若将铅同位素的时空演化与成矿作用有机结合,则其指示作用更明显。本文以东秦岭和粤北韶关地区的多种地质体的铅同位素成分资料为基础,探讨了铅同位素在构造环境、物质来源及成矿作用等方面的指示作用。     

Constrains on gravity wave induced diffusion in the middle atmosphere

A review of the important constraints on gravity wave induced diffusion of chemical tracers, heat, and momentum is given. Ground-based microwave spectroscopy measurements of H₂O and CO and rocket-based mass spectrometer measurements of Ar constrain the eddy diffusion coefficient for constituent transport (K _zz ) to be (1–3)×10⁵ cm²s^–1 in the upper mesosphere. Atomic oxygen data also limitsK _zz to a comparable value at the mesopause. From the energy balance of the upper mesosphere the eddy diffusion coefficient for heat transport (D _H ) is, at most 6×10⁵ cm²s^–1 at the mesopause and decreasing substantially with decreasing altitude. The available evidence for mean wind deceleration and the corresponding eddy diffusion coefficient for momentum stresses (D _M ) suggests that it is at least 1×10⁶ cm²s^–1, in the upper mesosphere. Consequently the eddy Prandtl number for macroscopic scale lengths is >3.