Two-dimensional Vertical Cross Section Research Articles

Rational determination of cone penetration test quantity in a two-dimensional vertical cross-section for site investigation

Site investigation plays a critical role in underground space development and geotechnical engineering practice, and it aims at collecting data for characterization of subsurface soil properties and stratification, which are often simplified as a two-dimensional (2D) vertical cross-section for subsequent geotechnical design and analysis in practice. In-situ tests, such as cone penetration tests (CPT), are often performed during site investigation. Existing tunneling or geotechnical guidelines or design codes only provide conceptual principle for determination of CPT quantity (e.g., accuracy of site investigation improves as the CPT quantity increases) and recommend some rule-of-thumb minimal quantities (e.g., at least one CPT every 100 m). No rational method is available to determine the minimal CPT quantity for achieving a target level of accuracy in an objective and quantitative manner. This paper aims to develop a rational method for addressing this issue. The method developed is based on an extensive parametric study that uses random field simulation to cover typical CPT data encountered in geotechnical engineering practice and Bayesian compressive sampling to interpolate the spatially varying, but correlated, CPT data along both depth and horizontal directions simultaneously. The proposed method is illustrated and validated using real CPT data, and it is shown to perform well.

Nitrification of the lowermost stratosphere during the exceptionally cold Arctic winter 2015–2016

Abstract. The Arctic winter 2015–2016 was characterized by exceptionally low stratospheric temperatures, favouring the formation of polar stratospheric clouds (PSCs) from mid-December until the end of February down to low stratospheric altitudes. Observations by GLORIA (Gimballed Limb Observer for Radiance Imaging of the Atmosphere) on HALO (High Altitude and LOng range research aircraft) during the PGS (POLSTRACC–GW-LCYCLE II–SALSA) campaign from December 2015 to March 2016 allow the investigation of the influence of denitrification on the lowermost stratosphere (LMS) with a high spatial resolution. Two-dimensional vertical cross sections of nitric acid (HNO3) along the flight track and tracer–tracer correlations derived from the GLORIA observations document detailed pictures of wide-spread nitrification of the Arctic LMS during the course of an entire winter. GLORIA observations show large-scale structures and local fine structures with enhanced absolute HNO3 volume mixing ratios reaching up to 11 ppbv at altitudes of 13 km in January and nitrified filaments persisting until the middle of March. Narrow coherent structures tilted with altitude of enhanced HNO3, observed in mid-January, are interpreted as regions recently nitrified by sublimating HNO3-containing particles. Overall, extensive nitrification of the LMS between 5.0 and 7.0 ppbv at potential temperature levels between 350 and 380 K is estimated. The GLORIA observations are compared with CLaMS (Chemical Lagrangian Model of the Stratosphere) simulations. The fundamental structures observed by GLORIA are well reproduced, but differences in the fine structures are diagnosed. Further, CLaMS predominantly underestimates the spatial extent of HNO3 maxima derived from the GLORIA observations as well as the overall nitrification of the LMS. Sensitivity simulations with CLaMS including (i) enhanced sedimentation rates in case of ice supersaturation (to resemble ice nucleation on nitric acid trihydrate (NAT)), (ii) a global temperature offset, (iii) modified growth rates (to resemble aspherical particles with larger surfaces) and (iv) temperature fluctuations (to resemble the impact of small-scale mountain waves) slightly improved the agreement with the GLORIA observations of individual flights. However, no parameter could be isolated which resulted in a general improvement for all flights. Still, the sensitivity simulations suggest that details of particle microphysics play a significant role for simulated LMS nitrification in January, while air subsidence, transport and mixing become increasingly important for the simulated HNO3 distributions towards the end of the winter.

Cone penetration test (CPT)-based subsurface soil classification and zonation in two-dimensional vertical cross section using Bayesian compressive sampling

A novel method is developed in this study for soil classification and zonation in a two-dimensional (2D) vertical cross section using cone penetration tests (CPTs). A CPT is usually performed vertically and the number of CPT soundings in a site is often limited in geotechnical engineering practice. It is, therefore, difficult to properly interpret CPT results along the horizontal direction or accurately estimate the horizontal correlation length of CPT data. The method proposed in this study bypasses the difficulty in estimating horizontal correlation length and provides proper identification of subsurface soil stratification (i.e., soil layer number is constant along horizontal direction) and zonation (i.e., soil layer number varies along horizontal direction) in a 2D vertical cross section directly from a limited number of CPT soundings. The proposed method consists of three key elements: 2D interpolation of CPT data using 2D Bayesian compressive sampling; determination of soil behavior type (SBT) using a SBT chart at every location in the 2D section, including locations with measurements and unsampled locations; and soil layer or zone delineation using an edge detection method. Both simulated and real data examples are used to illustrate the proposed method. Results show that the method performs well even when only five sets of CPT soundings are available.

Unusual chlorine partitioning in the 2015/16 Arctic winter lowermost stratosphere: observations and simulations

Abstract. The Arctic winter 2015/16 was characterized by cold stratospheric temperatures. Here we present a comprehensive view of the temporal evolution of chlorine in the lowermost stratosphere over the course of the studied winter. We utilize two-dimensional vertical cross sections of ozone (O3) and chlorine nitrate (ClONO2), measured by the airborne limb imager GLORIA (Gimballed Limb Observer for Radiance Imaging of the Atmosphere) during the POLSTRACC/GW-LCYCLE II/GWEX/SALSA campaigns, to investigate the tropopause region in detail. Observations from three long-distance flights in January, February, and March 2016 are discussed. ClONO2 volume mixing ratios up to 1100 pptv were measured at 380 K potential temperature in mesoscale structures. Similar mesoscale structures are also visible in O3 measurements. Both trace gas measurements are applied to evaluate simulation results from the chemistry transport model CLaMS (Chemical Lagrangian Model of the Stratosphere) and the chemistry–climate model EMAC (ECHAM5/MESSy Atmospheric Chemistry). These comparisons show agreement within the expected performance of these models. Satellite measurements from Aura/MLS (Microwave Limb Sounder) and SCISAT/ACE-FTS (Atmospheric Chemistry Experiment – Fourier Transform Spectrometer) provide an overview over the whole winter and information about the stratospheric situation above the flight altitude. Time series of these satellite measurements reveal unusually low hydrochloric acid (HCl) and ClONO2 at 380 K from the beginning of January to the end of February 2016, while chlorine monoxide (ClO) is strongly enhanced. In March 2016, unusually rapid chlorine deactivation into HCl is observed instead of deactivation into ClONO2, the more typical pathway for deactivation in the Arctic. Chlorine deactivation observed in the satellite time series is well reproduced by CLaMS. Sensitivity simulations with CLaMS demonstrate the influence of low abundances of O3 and reactive nitrogen (NOy) due to ozone depletion and sedimentation of NOy-containing particles, respectively. On the basis of the different altitude and time ranges of these effects, we conclude that the substantial chlorine deactivation into HCl at 380 K arose as a result of very low ozone abundances together with low temperatures. Additionally, CLaMS estimates ozone depletion of at least 0.4 ppmv at 380 K and 1.75 ppmv at 490 K, which is comparable to other extremely cold Arctic winters. We have used CLaMS trajectories to analyze the history of enhanced ClONO2 measured by GLORIA. In February, most of the enhanced ClONO2 is traced back to chlorine deactivation that had occurred within the past few days prior to the GLORIA measurement. In March, after the final warming, air masses in which chlorine has previously been deactivated into ClONO2 have been transported in the remnants of the polar vortex towards the location of measurement for at least 11 d.

Estimating Travel Time in Bank Filtration Systems from a Numerical Model Based on DTS Measurements.

An approach is presented to determine the seasonal variations in travel time in a bank filtration system using a passive heat tracer test. The temperature in the aquifer varies seasonally because of temperature variations of the infiltrating surface water and at the soil surface. Temperature was measured with distributed temperature sensing along fiber optic cables that were inserted vertically into the aquifer with direct push equipment. The approach was applied to a bank filtration system consisting of a sequence of alternating, elongated recharge basins and rows of recovery wells. A SEAWAT model was developed to simulate coupled flow and heat transport. The model of a two-dimensional vertical cross section is able to simulate the temperature of the water at the well and the measured vertical temperature profiles reasonably well. MODPATH was used to compute flowpaths and the travel time distribution. At the study site, temporal variation of the pumping discharge was the dominant factor influencing the travel time distribution. For an equivalent system with a constant pumping rate, variations in the travel time distribution are caused by variations in the temperature-dependent viscosity. As a result, travel times increase in the winter, when a larger fraction of the water travels through the warmer, lower part of the aquifer, and decrease in the summer, when the upper part of the aquifer is warmer.

Search for low-temperature solvents using nonplanar tie-lines

The technique of searching for the point of quaternary eutectics is imitated by constructing a series of two-dimensional vertical cross sections (method of tie-lines) by the example of a model of eutectic-type T-x-y-z diagram without solid-phase solubility. Errors that arise when mapping cross sections of experimentally studied diagrams are analyzed.

Closed‐Form Approximations for Two‐Dimensional Groundwater Age Patterns in a Fresh Water Lens

Simple closed-form approximations are presented for calculating the steady-state groundwater age distribution in two-dimensional vertical cross sections of idealized fresh water lenses overlying salt water, for aquifers that are vertically semi-infinite and of finite thickness. The approximations are developed on the basis of existing one-dimensional analytical solutions for travel-time calculation in fresh water lenses and approximate streamline formulations. The two-dimensional age distributions based on the closed-form solutions match convincingly with numerical simulations. As expected, notable deviations from the numerical solution are encountered at the groundwater flow divide and when submarine groundwater discharge occurs. Ratios of recharge over hydraulic conductivities are varied to explore how the magnitude of the deviations changes, and it is found that the approximate closed-form solutions perform well over a range of conditions found in natural systems.

Sources of oxygen flux in groundwater during induced bank filtration at a site in Berlin, Germany

The microbial degradation of pharmaceuticals found in surface water used for artificial recharge is strongly dependent on redox conditions of the subsurface. Furthermore the durability of production wells may decrease considerably with the presence of oxygen and ferrous iron due to the precipitation of trivalent iron oxides and subsequent clogging. Field measurements are presented for oxygen at a bank filtration site in Berlin, Germany, along with simplified calculations of different oxygen pathways into the groundwater. For a two-dimensional vertical cross-section, oxygen input has been calculated for six scenarios related to different water management strategies. Calculations were carried out in order to assess the amount of oxygen input due to (1) the infiltration of oxic lake water, (2) air entrapment as a result of water table oscillations, (3) diffusive oxygen flux from soil air and (4) infiltrating rainwater. The results show that air entrapment and infiltrating lake water during winter constitute by far the most important mechanism of oxygen input. Oxygen input by percolating rainwater and by diffusive delivery of oxygen in the gas phase is negligible. The results exemplify the importance of well management as a determining factor for water oscillations and redox conditions during artificial recharge.

Fringe-Controlled Natural Attenuation of Phenoxy Acids in a Landfill Plume: Integration of Field-Scale Processes by Reactive Transport Modeling

Data obtained from a field study of an aquifer contaminated by landfill leachate and related laboratory experiments suggest that natural attenuation of phenoxy acid herbicides such as mecoprop (MCPP) occurs in the transition zone between the anaerobic plume core and the overlying aerobic water body. The location of this transition zone is assumed to be primarily controlled by vertical transverse dispersion processes occurring downstream of the pollution source. A reactive transport modeling study was carried out to evaluate this conceptual model. The transport was simulated for a two-dimensional vertical cross section to quantify the combined physical, geochemical, and microbial processes that affect the fate of the phenoxy acid herbicides. The simulations, showing removal of phenoxy acids, an increase of phenoxy acid degraders in the fringe zone, and a dependency of the results on vertical transverse dispersivity, are compatible with the hypothesis of fringe-controlled aerobic biodegradation of the phenoxy acids.

The physical scale modelling of braided alluvial architecture and estimation of subsurface permeability

ABSTRACTThe quantitative modelling of fluvial reservoirs, especially in the stages of enhanced oil recovery, requires detailed three‐dimensional data at both the scale of the channel belt and within‐channel. Although studies from core, analogue outcrop and modern environments may partially meet these needs, they often cannot provide detail on the smaller‐scale (i.e. channel‐scale) heterogeneity, frequently suffer from limited three‐dimensional exposure and cannot be used to examine the influence of different variables on the process–deposit relationship. Physical modelling offers a complementary technique that can address many of these quantitative requirements and holds great future potential for integration with reservoir modelling. Physical modelling provides the potential to upscale results and derive reservoir information on three‐dimensional facies geometry, connectivity and permeability.This paper describes the development and use of physical modelling, which employs generic Froude‐scaling principles, in an experimental basin that permits aggradation in order to model the morphology and subsurface depositional stratigraphy of coarse‐grained braided rivers. An example is presented of a 1:50 scale model based on the braided Ashburton River, Canterbury Plains, New Zealand and the adjacent late Quaternary braided alluvium exposed in the coastal cliffs. Critically, a full, bimodal grain size distribution (20% sand and 80% gravel) was used to replicate the prototype, which allows the realistic reproduction of the surface morphology and importantly permits grain size sorting during deposition. Uncertainties associated with the compression of time, sediment mass balance and the hydrodynamics of the finest particle sizes do not appear to affect the reproducibility of stratigraphy between experimental and natural environments.Sectioning of the preserved sedimentary sequence in the physical model allows quantification of the geometry, shape, spatial distribution and internal sedimentary structure of the coarse‐ and fine‐grained facies. A six‐fold facies scheme is proposed for the model braided alluvium and a direct link is established between the grain size distribution and facies type: this allows permeability to be estimated for each facies, which can be mapped onto two‐dimensional vertical cross‐sections of the preserved stratigraphy. Results demonstrate the dominance of four facies based on permeability that range over three orders of magnitude in hydraulic conductivity. Quantification of such variability, and linkage to both vertical proportion curves for facies distribution and connectivity presents significant advantages over other methodologies and offers great potential for the modelling of heterogeneous braided river sediments at the within channel‐belt scale. This paper outlines how physical models may be used to develop high‐resolution, geologically‐accurate, object‐based reservoir simulation models.

Using Karhunen–Loéve decomposition and artificial neural network to model miscible fluid displacement in porous media

In this paper, we describe an approach to model fluid displacements in porous media that combines two powerful techniques, namely Karhunen–Loéve (KL) decomposition and artificial neural networks (ANNs). The KL decomposition, for data compression and feature identification, is used to extract coherent structures or eigenfunctions using fluid concentration maps obtained from fine-mesh numerical simulations of miscible fluid displacements of oil by solvent in a two-dimensional vertical cross-section. Twenty KL eigenfunctions that capture 98.8% of the total energy are extracted. Corresponding data coefficients are constructed by projecting the fluid concentration maps of the numerical simulations onto the KL eigenfunctions. Processing these data coefficients through an ANN is found to be a powerful tool in predicting the fluid displacements of the fine-mesh numerical simulations without actually performing these simulations.

Observing gravity wave activity in the mesopause region by means of airglow tomography

Abstract The possibility of reconstructing a two-dimensional vertical cross section of structures in mesospheric airglow by means of tomographic inversion is investigated. A chain of cameras in a line on the ground level is assumed and the data to be inverted consist of luminosity curves from scans of airglow images. Results of inverting simulated data indicates that reconstruction of tilted structures due to gravity waves should be possible. This method should therefore allow the determination of the vertical wavelength as well as the horizontal wave parameters. The desirable number of cameras in a chain is discussed. An example of results obtained from a two-camera experiment using all-sky optics is also presented.

Two-dimensional simulation of clastic and carbonate sedimentation, consolidation, subsidence, fluid flow, heat flow and solute transport during the formation of sedimentary basins

Geological processes that create sedimentary basins or act during their formation can be simulated using the public domain computer code `BASIN'. For a given set of geological initial and boundary conditions the sedimentary basin evolution is calculated in a forward modeling approach. The basin is represented in a two-dimensional vertical cross section with individual layers. The stratigraphic, tectonic, hydrodynamic and thermal evolution is calculated beginning at an initial state, and subsequent changes of basin geometry are calculated from sedimentation rates, compaction and pore fluid mobilization, isostatic compensation, fault movement and subsidence. The sedimentologic, hydraulic and thermal parameters are stored at discrete time steps allowing the temporal evolution of the basin to be analyzed. A maximum flexibility in terms of geological conditions is achieved by using individual program modules representing geological processes which can be switched on and off depending on the data available for a specific simulation experiment. The code incorporates a module for clastic and carbonate sedimentation, taking into account the impact of clastic sediment supply on carbonate production. A maximum of four different sediment types, which may be mixed during sedimentation, can be defined. Compaction and fluid flow are coupled through the consolidation equation and the nonlinear form of the equation of state for porosity, allowing nonequilibrium compaction and overpressuring to be calculated. Instead of empirical porosity-effective stress equations, a physically consistent consolidation model is applied which incorporates a porosity dependent sediment compressibility. Transient solute transport and heat flow are calculated as well, applying calculated fluid flow rates from the hydraulic model. As a measure for hydrocarbon generation, the Time–Temperature Index (TTI) is calculated. Three postprocessing programs are available to provide graphic output in PostScript format: BASINVIEW is used to display the distribution of parameters in the simulated cross-section of the basin for defined time steps. It is used in conjunction with the Ghostview software, which is freeware and available on most computer systems. AIBASIN provides PostScript output for Adobe Illustrator®, taking advantage of the layer-concept which facilitates further graphic manipulation. BASELINE is used to display parameter distribution at a defined well or to visualize the temporal evolution of individual elements located in the simulated sedimentary basin. The modular structure of the BASIN code allows additional processes to be included. A module to simulate reactive transport and diagenetic reactions is planned for future versions. The program has been applied to existing sedimentary basins, and it has also shown a high potential for classroom instruction, giving the possibility to create hypothetical basins and to interpret basin evolution in terms of sequence stratigraphy or petroleum potential.

A neural network prediction model of fluid displacements in porous media

This paper presents the development and design of an artificial neural network that is able to predict the breakthrough oil recovery of immiscible displacement of oil by water in a two-dimensional vertical cross section. The data used in training the neural network was obtained from the results of fine-mesh numerical simulations. Several network architectures were investigated and trained using the back propagation with momentum algorithm. The neural network that gave the best predictive performance was a two-hidden layer network with 8 neurons in the first hidden layer and 8 neurons in the second hidden layer. This network also performed well against a cross validation test. The reservoir simulation data used so far in the training process was for a homogeneous reservoir, the more general case is still under investigation.

An artificial neural network for the prediction of immiscible flood performance

Despite several decades of artificial neural network reserach in other engineering disciplines, only recently work has been reported on its use as a prediction tool in petroleum engineering applications. Existing methods for the prediction of fluid flow in porous medium include numerical simulation techniques and laboratory core flood experiments. Both of these methods are generally expensive and time consuming. However, neural networks, once successfully trained, can be used to predict reservoir performance in a short time with a personal computer. An artificial neural network was developed using data obtained from fine-mesh numerical simulation to predict the breakthrough oil recovery of immiscible displacement of oil by water in a two-dimensional vertical cross section. The network is able to predict the results of the fine-mesh numerical simulations without actually performing these simulation runs. Various neural network connections were investigated using the back-propagation with momentum algorithm for error minimization. This paper describes the design, development, and testing of the neural network.

A shallow attenuating anomaly inside the ring fracture of the Valles Caldera, New Mexico

Spectral ratios of teleseismic direct and scattered P waves observed in the Valles Caldera, New Mexico, show a systematic pattern of low amplitudes at sites inside the caldera relative to sites on or outside the ring fracture. Waveforms recorded at caldera stations are considerably more complex than those recorded outside the caldera. The data used in this study were collected during a passive seismic monitoring experiment conducted in 1987. Twenty-four teleseismic events were recorded on two linear arrays spanning the caldera. To first order, the pattern of low amplitudes did not vary with source incidence angle or azimuth of approach, and could not be explained by anomalous amplification at the ring fracture. This observation suggests the presence of a shallow, attenuating zone associated with the caldera fill material inside of the ring fracture. We estimated the general features of the caldera's near-surface structure for the two-dimensional vertical cross section beneath the array, using a modification of the Aki-Larner discrete-wavenumber method to forward model the observed amplitude variations. Our results indicate that the caldera fill material must be subdivided into at least two distinct zones: a strongly attenuating lower zone, extending to depths in excess of 4 km, and a mildly attenuating surface layer. To fit the data we had to assign an unrealistically low value to seismic Q in the deeper attenuating anomaly. We attribute this to the inability of the Aki-Larner method to account for strong re-direction of energy away from the caldera due to local heterogeneity that we could not include within the low-Q anomaly. This interpretation is consistent with the pervasive, fractured hydrothermal system that is known to exist in the caldera fill material.

Pseudorelative Permeabilities from Fractal Distributions

ABSTRACT Fractal distributions of porosity and permeability in two-dimensional vertical cross sections can be used to represent reservoir heterogeneity. We compare a hybrid finite-difference/streamtube method and pseudo-relative permeabilities as two methods of using these distributions in reservoir simulations. We simulate waterflood displacements in a field case study and some additional sensitivity cases. In all the cases studied, predicted cumulative oil production values agree within 5 percent between the two methods.

Cloud Top Liquid Water from Lidar Observations of Marine Stratocumulus

Marine stratus clouds were simultaneously observed by nadir Nd:YAG lidar measurements and in situ cloud physics measurements. A procedure was applied to derive the two-dimensional vertical cross section of the liquid water from within the cloud top lidar observations. A comparison to direct in-cloud liquid water observations gave good results. The liquid water retrieval was limited to an effective optical depth of 1.5. The true cloud optical thickness was also obtained from the retrieval procedure to a corresponding limit of 3.8. The optical thickness of the observed marine stratus clouds was predominantly below 3.0.