From discrete to continuum modelling of boundary value problems in geomechanics: An integrated FEM-DEM approach

Double-scale numerical methods constitute an effective tool for simultaneously representing the complex nature of geomaterials and treating real-scale engineering problems such as a tunnel excavation or a pressuremetre at a reasonable numerical cost. This paper presents an approach coupling discrete elements (DEM) at the microscale with finite elements (FEM) at the macroscale. In this approach, a DEM-based numerical constitutive law is embedded into a standard FEM formulation. In this regard, an exhaustive discussion is presented on how a 2D/3D granular assembly can be used to generate, step by step along the overall computation process, a consistent Numerically Homogenised Law. The paper also focuses on some recent developments including a comprehensive discussion of the efficiency of Newton-like operators, the introduction of a regularisation technique at the macroscale by means of a second gradient framework, and the development of parallelisation techniques to alleviate the computational cost of the proposed approach. Some real-scale problems taking into account the material spatial variability are illustrated, proving the numerical efficiency of the proposed approach and the benefit of a particle-based strategy.     

History matching time-lapse seismic data using the ensemble Kalman filter with multiple data assimilations

The ensemble Kalman filter (EnKF) has become a popular method for history matching production and seismic data in petroleum reservoir models. However, it is known that EnKF may fail to give acceptable data matches especially for highly nonlinear problems. In this paper, we introduce a procedure to improve EnKF data matches based on assimilating the same data multiple times with the covariance matrix of the measurement errors multiplied by the number of data assimilations. We prove the equivalence between single and multiple data assimilations for the linear-Gaussian case and present computational evidence that multiple data assimilations can improve EnKF estimates for the nonlinear case. The proposed procedure was tested by assimilating time-lapse seismic data in two synthetic reservoir problems, and the results show significant improvements compared to the standard EnKF. In addition, we review the inversion schemes used in the EnKF analysis and present a rescaling procedure to avoid loss of information during the truncation of small singular values.     

The Exploitation of Plant Resources by Early Homo sapiens: The Phytolith Record from Pinnacle Point 13B Cave,South Africa

The capacity to detect remnants of fire usage in archaeological sites, along with a better understanding of its use by prehistoric populations, can help us to shed light on hominin cognition, social organization, and technology. The application of phytolith studies to understand fire and its use can be widely applicable since plants typically contain phytoliths that are identifiable to different taxonomic levels. Pinnacle Point 13B (PP13B) is one of several South African Middle Stone Age (MSA) sites that shows an abundant presence of lenses of burnt material and intact hearths. Phytolith analyses at PP13B have been used to identify the types of plants used as fuel in the hearths. The phytoliths overall show high levels of alteration, and varying alteration is used to identify areas with higher alkaline conditions related to dripping water. In some areas of the site, there is excellent preservation of multicellular structures from the epidermal leaves of dicotyledonous plants, suggesting in situ structures with practically no chemical or post‐depositional alteration. A notable pattern is the abundance of dicotyledonous leaves from the rear of the cave, which might indicate specific leaf‐fuels for the fires, short‐term fire activities, or other actions such as cooking.     

Imaging X-ray Polarimeter for Solar Flares (IXPS)

We describe the design of a balloon-borne Imaging X-ray Polarimeter for Solar flares (IXPS). This novel instrument, a Time Projection Chamber (TPC) for photoelectric polarimetry, will be capable of measuring polarization at the few percent level in the 20?C50 keV energy range during an M- or X-class flare, and will provide imaging information at the ??10 arcsec level. The primary objective of such observations is to determine the directivity of nonthermal high-energy electrons producing solar hard X-rays, and hence to learn about the particle acceleration and energy release processes in solar flares. Secondary objectives include the separation of the thermal and nonthermal components of the flare X-ray emissions and the separation of photospheric albedo fluxes from direct emissions.     

Geochemical and isotopic insights into the assembly,evolution and disruption of a magmatic plumbing system before and after a cataclysmic caldera-collapse eruption at Ischia volcano (Italy)

New geochemical and isotopic data on volcanic rocks spanning the period ~75–50 ka BP on Ischia volcano, Italy, shed light on the evolution of the magmatic system before and after the catastrophic, caldera-forming Monte Epomeo Green Tuff (MEGT) eruption. Volcanic activity during this period was influenced by a large, composite and differentiating magmatic system, replenished several times with isotopically distinct magmas of deep provenance. Chemical and isotopic variations highlight that the pre-MEGT eruptions were fed by trachytic/phonolitic magmas from an isotopically zoned reservoir that were poorly enriched in radiogenic Sr and became progressively less radiogenic with time. Just prior to the MEGT eruption, the magmatic system was recharged by an isotopically distinct magma, relatively more enriched in radiogenic Sr with respect to the previously erupted magmas. This second magma initially fed several SubPlinian explosive eruptions and later supplied the climactic, phonolitic-to-trachytic MEGT eruption(s). Isotopic data, together with erupted volume estimations obtained for MEGT eruption(s), indicate that >5–10 km³ of this relatively enriched magma had accumulated in the Ischia plumbing system. Geochemical modelling indicates that it accumulated at shallow depths (4–6 km), over a period of ca. 20 ka. After the MEGT eruption, volcanic activity was fed by a new batch of less differentiated (trachyte-latite) magma that was slightly less enriched in radiogenic Sr. The geochemical and Sr–Nd-isotopic variations through time reflect the upward flux of isotopically distinct magma batches, variably contaminated by Hercynian crust at 8–12 km depth. The deep-sourced latitic to trachytic magmas stalled at shallow depths (4–6 km depth), differentiated to phonolite through crystal fractionation and assimilation of a feldspar-rich mush, or ascended directly to the surface and erupted.     

Re-interpreting the Devonian ophiolites involved in the Variscan suture: U–Pb and Lu–Hf zircon data of the Moeche Ophiolite (Cabo Ortegal Complex,NW Iberia)

New U–Pb zircon data of a mylonitic greenschist from the Moeche Ophiolite, one of the mafic units involved in the Variscan suture in the Cabo Ortegal Complex (NW of the Iberian Massif), yielded an age of 400 ± 3 Ma. Consequently, this unit can be considered one of the Devonian ophiolites, the most extended group of oceanic units in the Variscan belt. The mafic rocks show transitional compositions between N-MORB and island-arc tholeiites, although Lu–Hf isotope signatures of its zircons clearly indicate contribution from an old continental source. εHf values in the analysed zircons are negative (generally below εHf = ?5), and thence, they are not compatible with their generation from a juvenile mantle source. Accordingly, the igneous protoliths were generated in a setting where juvenile mafic magmas interacted with an old continental crust. The Devonian ophiolites from the Variscan suture have been repeatedly interpreted as remnants of the Rheic Ocean. However, the presence of a continental source in the origin of the mafic rocks of the Moeche Ophiolite allows discarding an intraoceanic setting for their generation, at least for the NW Iberian counterparts. The tectonic setting for the Devonian ophiolites of NW Iberia is very likely represented by an ephemeral oceanic basin opened within a continental realm. Herein, the real Rheic Ocean suture could only be located west of the terrane represented by the upper units of the allochthonous complexes. Apparently that suture is not represented in NW Iberia.     

A universal nutrient application strategy for the bioremediation of oil-polluted beaches

Biostimulation by nutrient application is a viable technology for restoring oil-contaminated beaches. Maximizing the nutrient residence time is key for achieving a rapid cost-effective cleanup. We considered the nutrient injection strategy through a perforated pipe at the high tide line and we simulated numerically beach hydraulics, which allowed us to estimate the optimal injection flow rate of nutrient solution. Our results indicate that the optimal application is one that starts following the falling high tide and lasts for half tidal cycle. The saturated wet-front of the nutrient solution on the beach surface would move seaward with the same speed of the falling tide keeping a constant distance with the tide line. The numerical results were generalized to beaches of wide ranges of hydraulic and tidal properties using a novel dimensionless formulation for water flow and solute transport in porous media. Nomographs were presented to provide the flow rate based on four parameters: The beach slope and hydraulic conductivity, and tidal amplitude and period.     

Closed-loop reservoir management on the Brugge test case

This paper proposes an augmented Lagrangian method for production optimization in which the cost function to be maximized is defined as an augmented Lagrangian function consisting of the net present value (NPV) and all the equality and inequality constraints except the bound constraints. The bound constraints are dealt with using a trust-region gradient projection method. The paper also presents a way to eliminate the need to convert the inequality constraints to equality constraints with slack variables in the augmented Lagrangian function, which greatly reduces the size of the optimization problem when the number of inequality constraints is large. The proposed method is tested in the context of closed-loop reservoir management benchmark problem based on the Brugge reservoir setup by TNO. In the test, we used the ensemble Kalman filter (EnKF) with covariance localization for data assimilation. Production optimization is done on the updated ensemble mean model from EnKF. The production optimization resulted in a substantial increase in the NPV for the expected reservoir life compared to the base case with reactive control.     

Isobar Attenuation Using Anion–Gas Reactions for Accelerator Mass Spectrometry and Application to 36Cl

Preliminary tests of a prototype radio-frequency quadrupole (RFQ) collision cell system, known as an isobar separator for anions (ISA), for the removal of isobaric interferences for accelerator mass spectrometry (AMS) and for studies of anion–gas interactions are reported. The ISA decelerated a mass-analysed beam of anions from an energy (∼ 20 keV) typically generated by an AMS ion source to < 10 eV. RFQs and electrostatic lenses then guided the ions through the collision cell where ion-gas collisions reduced both the energy and energy spread of the ion beam (cooled the ions) and ion-gas reactions attenuated most of the unwanted isobars. The anions were then re-accelerated to their original energy for injection into the rest of the AMS system. With the ISA installed on a full 3 MV AMS system, attenuations of ³²S^-, ¹²C₃^- and ³⁹K^- by six, seven and greater than ten orders of magnitude, respectively were achieved using 0.7–1 Pa NO₂ gas in the collision cell, while maintaining approximately 10–30% of the chlorine anion transmission. A further measurement of a ³⁶Cl/Cl = 4.1 × 10^-11 RM is also described. The results suggest that the ³⁶Cl/Cl lower detection limit of the current system was 10^-14–10^-15 for samples that could be prepared with S/Cl ratios below 10 μg g^-1.