Geostatistical analysis of field hydraulic conductivity in compacted clay

Hydraulic conductivity(K) of fractured or porous materials is associated intimately with water flow and chemical transport. Basic concepts imply uniform flux through a homogeneous cross-sectional area. If flow were to occur only through part of the area, actual rates could be considerably different. Because laboratory values ofK in compacted clays seldom agree with field estimates, questions arise as to what the true values ofK are and how they should be estimated. Hydraulic conductivity values were measured on a 10×25 m elevated bridge-like platform. A constant water level was maintained for 1 yr over a 0.3-m thick layer of compacted clay, and inflow and outflow rates were monitored using 10×25 grids of 0.3-m diameter infiltration rings and outflow drains subtending approximately 1×1 m blocks of compacted clay. Variography of inflow and outflow data established relationships between cores and blocks of clay, respectively. Because distributions of outflow rates were much less and bore little resemblance to the distributions of break-through rates based on tracer studies, presence of macropores and preferential flow through the macropores was suspected. Subsequently, probability kriging was applied to reevaluate distribution of flux rates and possible location of macropores. Sites exceeding a threshold outflow of 100×10^–9 m/s were classified as outliers and were assumed to probably contain a significant population of macropores. Different sampling schemes were examined. Variogram analysis of outflows with and without outliers suggested adequacy of sampling the site at 50 randomly chosen locations. Because of the potential contribution of macropores to pollutant transport and the practical necessity of extrapolating small plot values to larger areas, conditional simulations with and without outliers were carried out. Simulated scenarios based on all available data compared well with conditional simulations based on randomly chosen locations.This paper was presented at MGUS 87 Conference, Redwood City, California, 14 April 1987.     

Kepler Equation solver

Kepler's Equation is solved over the entire range of elliptic motion by a fifth-order refinement of the solution of a cubic equation. This method is not iterative, and requires only four transcendental function evaluations: a square root, a cube root, and two trigonometric functions. The maximum relative error of the algorithm is less than one part in 10¹⁸, exceeding the capability of double-precision computer arithmetic. Roundoff errors in double-precision implementation of the algorithm are addressed, and procedures to avoid them are developed.     

SOFTWARE REVIEWS

software reviews are in this article BMDP/PC CONCEPTS COMPUTERIZED ATLAS EVALUATION OF ALTERNATE PROPOSALS (EAP) EXTRA: EXPRESS TRANSIT ANALYSIS MATC CAD MOTORS TRANSPORTATION PACKAGE PCMAP: Thematic Mapping Software for the IBM Personal Computer, Version 2.3 PLACE TRANPLAN . US-ATLAS     

Inverse problem in hydrogeology

The state of the groundwater inverse problem is synthesized. Emphasis is placed on aquifer characterization, where modelers have to deal with conceptual model uncertainty (notably spatial and temporal variability), scale dependence, many types of unknown parameters (transmissivity, recharge, boundary conditions, etc.), nonlinearity, and often low sensitivity of state variables (typically heads and concentrations) to aquifer properties. Because of these difficulties, calibration cannot be separated from the modeling process, as it is sometimes done in other fields. Instead, it should be viewed as one step in the process of understanding aquifer behavior. In fact, it is shown that actual parameter estimation methods do not differ from each other in the essence, though they may differ in the computational details. It is argued that there is ample room for improvement in groundwater inversion: development of user-friendly codes, accommodation of variability through geostatistics, incorporation of geological information and different types of data (temperature, occurrence and concentration of isotopes, age, etc.), proper accounting of uncertainty, etc. Despite this, even with existing codes, automatic calibration facilitates enormously the task of modeling. Therefore, it is contended that its use should become standard practice.

---

Resumen Se sintetiza el estado del problema inverso en aguas subterráneas. El énfasis se ubica en la caracterización de acuíferos, donde los modeladores tienen que enfrentar la incertidumbre del modelo conceptual (principalmente variabilidad temporal y espacial), dependencia de escala, muchos tipos de parámetros desconocidos (transmisividad, recarga, condiciones limitantes, etc), no linealidad, y frecuentemente baja sensibilidad de variables de estado (típicamente presiones y concentraciones) a las propiedades del acuífero. Debido a estas dificultades, no puede separarse la calibración de los procesos de modelado, como frecuentemente se hace en otros campos. En su lugar, debe de visualizarse como un paso en el proceso de entendimiento del comportamiento del acuífero. En realidad, se muestra que los métodos reales de estimación de parámetros no difieren uno del otro en lo esencial, aunque sí pueden diferir en los detalles computacionales. Se discute que existe amplio espacio para la mejora del problema inverso en aguas subterráneas: desarrollo de códigos amigables al usuario, acomodamiento de variabilidad a través de geoestadística, incorporación de información geológica y diferentes tipos de datos (temperatura, presencia y concentración de isótopos, edad, etc), explicación apropiada de incertidumbre, etc. A pesar de esto, aún con los códigos existentes, la calibración automática facilita enormemente la tarea de modelado. Por lo tanto, se sostiene que su uso debería de convertirse en práctica standard.

---

Résumé Létat du problème inverse des eaux souterraines est synthétisé. Laccent est placé sur la caractérisation de laquifère, où les modélisateurs doivent jouer avec lincertitude des modèles conceptuels (notamment la variabilité spatiale et temporelle), les facteurs déchelle, plusieurs inconnues sur différents paramètres (transmissivité, recharge, conditions aux limites, etc.), la non linéarité, et souvent la sensibilité de plusieurs variables détat (charges hydrauliques, concentrations) des propriétés de laquifère. A cause de ces difficultés, le calibrage ne peut être séparé du processus de modélisation, comme cest le cas dans dautres cas de figure. Par ailleurs, il peut être vu comme une des étapes dans le processus de détermination du comportement de laquifère. Il est montré que les méthodes dévaluation des paramètres actuels ne diffèrent pas si ce nest dans les détails des calculs informatiques. Il est montré quil existe une large panoplie de techniques d ‹inversion : codes de calcul utilisables par tout-un-chacun, accommodation de la variabilité via la géostatistique, incorporation dinformations géologiques et de différents types de données (température, occurrence, concentration en isotopes, âge, etc.), détermination de lincertitude. Vu ces développements, la calibration automatique facilite énormément la modélisation. Par ailleurs, il est souhaitable que son utilisation devienne une pratique standardisée.

     

Numerical Exploration of Chermnykh’s Problem

We study the equilibrium points and the zero-velocity curves of Chermnykh’s problem when the angular velocity ω varies continuously and the value of the mass parameter is fixed. The planar symmetric simple-periodic orbits are determined numerically and they are presented for three values of the parameter ω. The stability of the periodic orbits of all the families is computed. Particularly, we explore the network of the families when the angular velocity has the critical value ω = 2√2 at which the triangular equilibria disappear by coalescing with the collinear equilibrium point L₁. The analytic determination of the initial conditions of the family which emanate from the Lagrangian libration point L₁ in this case, is given. Non-periodic orbits, as points on a surface of section, providing an outlook of the stability regions, chaotic and escape motions as well as multiple-periodic orbits, are also computed. Non-linear stability zones of the triangular Lagrangian points are computed numerically for the Earth–Moon and Sun–Jupiter mass distribution when the angular velocity varies.     

Impact of the Mass Parameter on Particle Dynamics in a Ring Configuration of N Bodies

Parameters play a very important and determinative role in the dynamics of a dynamical system as well as in the formation of its particular characteristics. In this paper we investigate the way in which a large scale variation of the mass parameter, influences the behavior of a mass-less particle which moves in the vicinity of a ring arrangement of N-bodies. More precisely, we study the impact of this parameter on periodic motions and their characteristics.     

高煤级煤储层煤层气产能“瓶颈”问题研究

基于山西沁水盆地高煤级煤储层宏观裂隙、显微裂隙的连续观测,孔隙的系统测量,结合应力渗透率、气-水相对渗透率、吸附膨胀等实验成果,分析了高煤级煤储层三级渗流特征,探讨了有效应力和煤基质收缩对高煤级煤储层渗透率的耦合作用,系统揭示了在地面排水降压开发煤层气的过程中,高煤级煤储层初期产气量高,数月后急剧衰减之“瓶颈”现象,找出了造成高煤级煤储层产气缺陷的根本原因。鉴于高煤级煤储层物性的特殊性,指出了高煤级煤储层煤层气开发的技术和措施。     

Global Bifurcations of Periodic Solutions of the Restricted Three Body Problem

We describe global bifurcations from the libration points of non-stationary periodic solutions of the restricted three body problem. We show that the only admissible continua of non-stationary periodic solutions of the planar restricted three body problem, bifurcating from the libration points, can be the short-period families bifurcating from the Lagrange equilibria L ₄, L ₅. We classify admissible continua and show that there are possible exactly six admissible continua of non-stationary periodic solutions of the planar restricted three body problem. We also characterize admissible continua of non-stationary periodic solutions of the spatial restricted three body problem. Moreover, we combine our results with the Déprit and Henrard conjectures (see [8]), concerning families of periodic solutions of the planar restricted three body problem, and show that they can be formulated in a stronger way. As the main tool we use degree theory for SO(2)-equivariant gradient maps defined by the second author in [25].This revised version was published online in October 2005 with corrections to the Cover Date.     

Non-Continuation of Integrals of the Rotating Two-body Problem in Hill's Lunar Problem

We consider Hill's lunar problem as a perturbation of the integrable two-body problem. For this we avoid the usual normalization in which the angular velocity of the rotating frame of reference is put equal to unity and consider as the perturbation parameter. We first express the Hamiltonian H of Hill's lunar problem in the Delaunay variables. More precisely we deduce the expressions of H along the orbits of the two-body problem. Afterwards with the help of the conserved quantities of the planar two-body problem (energy, angular momentum and Laplace–Runge–Lenz vector) we prove that Hill's lunar problem does not possess a second integral of motion, independent of H, in the sense that there exist no analytic continuation of integrals, which are linear functions of in the rotating two-body problem. In connection with the proof of this main result we give a further restrictive statement to the nonintegrability of Hill's lunar problem.     

The restricted two-body problem in constant curvature spaces

We perform the bifurcation analysis of the Kepler problem on and . An analog of the Delaunay variables is introduced. We investigate the motion of a point mass in the field of a Newtonian center moving along a geodesic on and (the restricted two-body problem). For the case of a small curvature, the pericenter shift is computed using the perturbation theory. We also present the results of numerical analysis based on an analogy with the motion of a rigid body.