High-dimensional geostatistical history matching

Hydrocarbon reservoir modelling and characterisation is a challenging subject within the oil and gas industry due to the lack of well data and the natural heterogeneities of the Earth’s subsurface. Integrating historical production data into the geo-modelling workflow, commonly designated by history matching, allows better reservoir characterisation and the possibility of predicting the reservoir behaviour. We present herein a geostatistical-based multi-objective history matching methodology. It starts with the generation of an initial ensemble of the subsurface petrophysical property of interest through stochastic sequential simulation. Each model is then ranked according the match between its dynamic response, after fluid flow simulation, and the observed available historical production data. This enables building regionalised Pareto fronts and the definition of a large ensemble of optimal subsurface Earth models that fit all the observed production data without compromising the exploration of the uncertainty space. The proposed geostatistical multi-objective history matching technique is successfully implemented in a benchmark synthetic reservoir dataset, the PUNQ-S3, where 12 objectives are targeted.     

Comparison between cell-centered and nodal-based discretization schemes for linear elasticity

In this paper, we study newly developed methods for linear elasticity on polyhedral meshes. Our emphasis is on applications of the methods to geological models. Models of subsurface, and in particular sedimentary rocks, naturally lead to general polyhedral meshes. Numerical methods which can directly handle such representation are highly desirable. Many of the numerical challenges in simulation of subsurface applications come from the lack of robustness and accuracy of numerical methods in the case of highly distorted grids. In this paper, we investigate and compare the Multi-Point Stress Approximation (MPSA) and the Virtual Element Method (VEM) with regard to grid features that are frequently seen in geological models and likely to lead to a lack of accuracy of the methods. In particular, we look at how the methods perform near the incompressible limit. This work shows that both methods are promising for flexible modeling of subsurface mechanics.     

An orthotropic interface damage model for simulating drying processes in soils

The study of drying process in soils has received an increased attention in the last few years. This is very complex phenomenon that generally leads to the formation and propagation of desiccation cracks in the soil mass. In recent engineering applications, high aspect ratio elements have proved to be well suited to tackle this type of problem using finite elements. However, the modeling of interfaces between materials with orthotropic properties that generally exist in this type of problem using standard (isotropic) constitutive model is very complex and challenging in terms of the mesh generation, leading to very fine meshes that are intensive CPU demanding. A novel orthotropic interface mechanical model based on damage mechanics and capable of dealing with interfaces between materials in which the strength depends on the direction of analysis is proposed in this paper. The complete mathematical formulation is presented together with the algorithm suggested for its numerical implementation. Some simple yet challenging synthetic benchmarks are analyzed to explore the model capabilities. Laboratory tests using different textures at the contact surface between materials were conducted to evaluate the strengths of the interface in different directions. These experiments were then used to validate the proposed model. Finally, the approach is applied to simulate an actual desiccation test involving an orthotropic contact surface. In all the application cases the performance of the model was very satisfactory.     

Groundwater characteristics within loessic deposits: the coastal springs of Los Acantilados,Mar del Plata,Argentina

The urban growth of the southern neighborhoods of Mar del Plata City provoked significant changes in the groundwater balance of the loessic sequences. These regional loessic levels with a significant portion of volcanic ash layers were reported subject to fluoride and nitrate concentrations. Residential houses pump from sands located 70 m depth and withdraw the sewages to depths less than 5 m. These effects cause significant local and seasonal (summer) increments of the water table outcropping via springs at certain unconformities of the coastal cliffs. A mathematical model was applied to analyze the water level lowering at the productive levels, while there is a decrease in the quality of the upper levels subject to waste discharges. Much of this groundwater flow is concentrated in unconformities between different types of sediments. Human activities have affected the aquifer dynamics increasing the groundwater pumping rates and the return velocity of the sewages. This should be considered in the management strategies of coastal hydric resources.     

Soil wettability in ground engineering: fundamentals,methods, and applications

Wettability is a fundamental property controlling the extent of wetting in flat and granular solids. In natural soils, wettability affects a wide variety of processes including infiltration, preferential flow and surface runoff. In mineral processing, wettability is paramount in enhancing the efficiency of separation of minerals from gangue. The manipulation of surface wettability is equally crucial in many industrial applications. For instance, superhydrophobic surfaces are those on which water drops roll off easily and as such are used for self-cleaning applications. Therefore, while wettability is strongly cross-disciplinary, its evolution has been discipline-specific with a direct extrapolation or transfer of concepts, approaches, and methods to ground engineering unlikely to remain valid. This paper synthesizes relevant aspects from surface chemistry, materials science, mining engineering, and soil science, and discusses their implications within the context of new granular materials that resist wetting, for use in barriers or ground improvement and, in unsaturated soils, where the effects of wettability have been documented.     

Bioaccumulation of heavy metals in <Emphasis Type="Italic">Limnobium laevigatum</Emphasis> and <Emphasis Type="Italic">Ludwigia peploides</Emphasis>: their phytoremediation potential in water contaminated with heavy metals

Contamination with heavy metals in surface and groundwater is a threat to human health and ecosystems. Due to this, the need arises to remediate water polluted through ecological and profitable technologies, such as phytoremediation. The objective of the work was to evaluate the concentration of lead (Pb) and zinc (Zn) in the floating macrophytes Limnobium laevigatum and Ludwigia peploides, after being exposed to contaminated water experimentally. In this way to be able to determine if these plants have mechanisms that allow them to accumulate the metals in the roots and to perform the translocation of these to different vegetative organs, L. laevigatum and L. peploides were placed in solutions contaminated with Pb ([Pb]?=?5 mg/l) and Zn ([Zn]?=?20 mg/l). The concentrations of metals in water, root and leaf samples were evaluated as a function of time (0, 1, 2 and 4 days). The determination of the metals was performed by the atomic absorption spectrophotometry technique. After 4 days of exposure to Pb and Zn, the plants showed high metal removal efficiencies of water, more to 70% in all cases. Pb was accumulated fundamentally by roots, while Zn was accumulated more in the leaves. In addition, the bioconcentration and translocation factors for each metal were calculated.     

Assessing the performance of WRF model in predicting high-impact weather conditions over Central and Western Africa: an ensemble-based approach

Natural Hazards - For numerical weather prediction over a particular region, it is important to know the best combination of physical parameterizations available in the considered modelling frame...     

Image-aided platform orientation determination with a GNSS/low-cost IMU system using robust-adaptive Kalman filter

The availability of miniaturized sensors with enhanced capabilities, new methods for image processing, and easy access to small and low-weight airborne platforms for data acquisition, including unmanned vehicles, opens new possibilities for geodetic navigation applications and developing new developments in sensor fusion. In this context, the development of efficient methods, based on low-cost sensors, to extract precise georeferenced information from digital cameras is of utmost interest. We present a method to improve the performance of the integration of GNSS/low-cost IMU by exploiting the orientation changes retrieved from digital images. In this work, a robust-adaptive Kalman filter is also introduced to further improve the performance of the method deployed. The adaptive factor and the robust factor accomplished are determined by innovation information and the threshold value of orientation changes between consecutive images. Results from airborne tests used to assess the performance of the method are presented. The results show that using a non-metric camera, the Euler angle estimation accuracy of the GNSS/low-cost IMU integration can be improved to be close to 0.5 degree and an additional improvement, which can reach 59%, can be achieved after using the robust-adaptive Kalman filter.     

The Colider and Roosevelt Volcanic rocks (sw amazonian craton): geochemistry and sm-nd isotope characteristics of a silicic large igneous province

The volcanic rocks of the Colíder and Roosevelt formations are extensively exposed in the south-central portion of the Amazonian Craton where effusive and pyroclastic rocks have been mapped. Both units, topped by chemical sediments and oceanic facies as rhyolite and andesite lavas, rhyodacite, and porphyritic dacite, with frequent intercalations of pyroclastic and epiclastic deposits. Whole-rock geochemistry for 55 samples of rhyolitic to andesitic composition suggests the involvement of fertile mantle-derived components with E-MORB to OIB compositions. The analyzed rocks display calc-alkaline to shoshonitic affinity consistent with generation related to an active continental margin. The whole-rock Sm-Nd isotope data from selected felsic volcanic rocks of the Colíder and Roosevelt formations yield negative initial ε_Nd values between –3 and –9, indicating the predominantly crustal nature of the parental magmas with early Archean to late Paleoproterozoic (ca. 2.5–2.0 Ga) depleted mantle model ages.

     

Single-particle crushing strength under different relative humidity conditions

This paper presents a fundamental study on the effect of the relative humidity on the rockfill crushing strength. This aspect plays an important role in the mechanical behaviour of rockfill, and it is known that certain characteristics of the granular materials, such as compressibility and shear strength, depend on the confining stress, which is a function of the particles crushing. An increased interest has been observed regarding the effect of the relative humidity in the mechanical behaviour of rockfill. Unfortunately, limited research has been conducted until now regarding the study of individual particle crushing. Therefore, this paper thoroughly investigated particle crushing, by performing single-particle crushing tests on rockfill particles divided into four size ranges, under different relative humidity conditions. The experimental results reveal a considerable influence of the relative humidity in the studied rockfill particles, whose strength of the particles with the greatest dimensions in saturated conditions was reduced by half. Consistent macro-mechanical evidence demonstrates that particle’s size and relative humidity conditions depict the most important factors that influence particle crushing strength.