Evaluation of Combined Direct-Push Methods Used for Aquifer Model Generation

Most established methods to characterize aquifer structure and hydraulic conductivities of hydrostratigraphical units are not capable of delivering sufficient information in the spatial resolution that is desired for sophisticated numerical contaminant transport modeling and adapted remediation design. With hydraulic investigation methods based on the direct-push (DP) technology such as DP slug tests, DP injection logging, and the hydraulic profiling tool, it is possible to rapidly delineate hydrogeological structures and estimate their hydraulic conductivity in shallow unconsolidated aquifers without the need for wells. A combined application of these tools was used for the investigation of a contaminated German refinery site and for the setup of hydraulic aquifer models. The quality of DP investigation and the models was evaluated by comparisons of tracer transport simulations using these models and measured breakthroughs of two natural gradient tracer tests. Model scenarios considering the information of all tools together showed good reproduction of the measured breakthroughs, indicating the suitability of the approach and a minor impact of potential technical limitations. Using the DP slug tests alone yielded significantly higher deviations for the determined hydraulic conductivities compared to considering two or three of the tools. Realistic aquifer models developed on basis of such combined DP investigation approaches can help optimize remediation concepts or identify flow regimes for aquifers with a complex structure.     

Prestack shot-gather depth migration by a rigid flow of Gaussian wave packets

We introduce a new method for prestack depth migration of seismic common-shot gathers. The computational procedure follows standard steps of the reverse-time migration, i.e., downward continuation of the source and the receiver wavefields, followed by application of an imaging condition (e.g. zero-lag cross-correlation of these fields). In our method we first find a sparse data representation with a small number of Gaussian wave packets. We then approximate the downward wavefield propagation (for the source and the receiver fields) by a rigid flow of these wave packets along seismic rays. In this case, the wave packets are simply translated and rotated according to the ray geometry. One advantage of using Gaussian wave packets is that analytic formulas can be used for translation, rotation, and the application of the cross-correlation imaging condition. Moreover, they allow more sparse representations than competing methods. Finally we formulate a computationally and memory efficient migration procedure, as only few rays have to be traced, and since it is cheap to compute the cross-correlation for the intersecting rays.     

A discrete particle representation of hillslope hydrology: hypothesis testing in reproducing a tracer experiment at Gårdsjön,Sweden

Despite the long history of the continuum equation approach in hydrology, it is not a necessary approach to the formulation of a physically based representation of hillslope hydrology. The Multiple Interacting Pathways (MIPs) model is a discrete realization that allows hillslope response and transport to be simultaneously explored in a way that reflects the potential occurrence of preferential flows and lengths of pathways. The MIPs model uses random particle tracking methods to represent the flow of water within the subsurface alongside velocity distributions that acknowledge preferential flows and transition probability matrices, which control flow pathways. An initial realization of this model is presented here in application to a tracer experiment carried out in Gårdsjön, Sweden. The model is used as an exploratory tool, testing several hypotheses in relation to this experiment. Copyright © 2011 John Wiley & Sons, Ltd.     

Distribution surfaces in solvent extraction systems

In the solvent extraction of metals, the distribution of the metal is usually dependent on the concentrations of several different species in the two phases. Despite this fact, distribution isotherms are usually presented as plots of C_o as a function of C_aq and these are used for a McCabe-Thiele type evaluation of a number of stages and flow rates, etc., each individual isotherm being valid only for one set of initial conditions; this treatment limits the possibilities of optimizing the system. In the present paper an extension of three-dimensional plots is suggested with special reference to the Cu/LIX64N extraction system. The distribution surface is defined and the general properties of the plots are outlined.A procedure is outlined whereby the surface can be represented on a two dimensional plot and a theoretical model for the surface can be superimposed.The paper also shows how distribution ratios (D values) are related to the surface.     

Drag force on a circular cylinder in unseparated laminar high Reynolds number nonlinear random oscillatory flow

The stochastic properties of the drag force maxima on a circular cylinder subjected to nonlinear random waves are investigated. Unseparated laminar high Reynolds number flow is considered. A simplified approach based on second order Stokes waves is presented, including the sum-frequency effect only. It is demonstrated how a drag force formula valid for regular linear waves can be used to find the cumulative distribution function of individual drag force maxima for nonlinear irregular waves. Here the [Wang, 1968] drag force coefficient is used.     

Industrial minerals and rocks, aggregates and natural stones in the Nordic countries

The Nordic countries, including Greenland, have a long tradition in mining. The industrial minerals sector is expanding in most Nordic countries and extensive development has taken place during the last few years. The main commodities mined are carbonate rocks, quartz, feldspar, apatite, olivine and talc.
A number of diferent types of dimension stones are quarried in all countries. Rock aggregates are increasingly important, replacing sand and gravel aggregate as construction materials in some countries due to the need to protect ground water supplies.     

Glacial stratigraphy of East and Central Jutland

Glaciostratigraphic investigations at one key locality (Haldum), 9 major and about 160 minor localities in East and Central Jutland, Denmark, together with laboratory work, have led to the establishment of a stratigraphy consisting of 10 till units, usually separated by meltwater deposits. The stratigraphy is in some degree supported by thermoluminescence datings. The complete sequence includes one till unit with associated meltwater deposits of Menapian age, three till units with intercalated meltwater deposits of Elsterian age, marine sediments deposited during the Holsteininan, and three till units with intercalated of Elsterian age, marine sediments deposited during the Holsteinian, and three till units with intercalated glaciofluvial sedimants of Saalian age. Eemian deposits are present above this level, and the whole sequence is capped by till and meltwater deposits related to three glacial advances during the Weichselian.     

Buried pockmarks on the Top Chalk surface of the Danish North Sea and their potential significance for interpreting palaeocirculation patterns

Three-dimensional (3D) seismic data from the southern Danish North Sea were used to analyse the morphology and spatial distribution of depressions in the Danian Chalk deposit. Previously, these depressions were either interpreted as karst structures or pockmarks. The observed depressions occur in an interval from 25 ms below to 12 ms above the Top Chalk surface. Three types of depressions were differentiated based on their plan-view geometry and their degree of symmetry: Type 1, comprising sub-circular and symmetrical depressions, is the dominant group (ca. 70 %). Type 2, elongated and symmetrical depressions, represents only a small fraction (ca. 5 %). The elongated and asymmetrical depressions of Type 3 compose ca. 25 %. In cross section, each depression type can be either characterised by a V- or a U-shape. The maximum size of the depressions ranges from 50 to 580 m, with an average internal depth of 10 m. We interpret the depressions as pockmarks formed by the expulsion of biogenic or thermogenic fluids at the Danian seafloor. Likely, the initial form of the pockmarks has been circular (Type 1) and was subsequently modified for Types 2 and 3 to an elongated form by currents. The long axis of the pockmarks is interpreted to represent the effective current direction. The inferred direction is sub-parallel to the palaeobathymetric contours. The 3D seismic interpretation of pockmarks presented in this paper contributes to the understanding of fluid migration and palaeocirculation patterns during the sedimentation of the terminal Chalk Group in the southern Danish North Sea.     

Reservoir description using dynamic parameterisation selection with a combined stochastic and gradient search

There is a correspondence between flow in a reservoir and large scale permeability trends. This correspondence can be derived by constraining reservoir models using observed production data. One of the challenges in deriving the permeability distribution of a field using production data involves determination of the scale of resolution of the permeability. The Adaptive Multiscale Estimation (AME) seeks to overcome the problems related to choosing the resolution of the permeability field by a dynamic parameterisation selection. The standard AME uses a gradient algorithm in solving several optimisation problems with increasing permeability resolution. This paper presents a hybrid algorithm which combines a gradient search and a stochastic algorithm to improve the robustness of the dynamic parameterisation selection. At low dimension, we use the stochastic algorithm to generate several optimised models. We use information from all these produced models to find new optimal refinements, and start out new optimisations with several unequally suggested parameterisations. At higher dimensions we change to a gradient-type optimiser, where the initial solution is chosen from the ensemble of models suggested by the stochastic algorithm. The selection is based on a predefined criterion. We demonstrate the robustness of the hybrid algorithm on sample synthetic cases, which most of them were considered insolvable using the standard AME algorithm.