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.     

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.     

Would different methodologies for assessing carbon leakage exposure lead to different risk levels? A case study of the Brazilian industry

Energy-intensive industries play an important role in low-carbon development, being particularly exposed to climate policies. Concern over possible carbon leakage in this sector poses a major challenge for designing effective carbon pricing instruments (CPI). Different methodologies for assessing carbon leakage exposure are currently used by different jurisdictions, each of them based on different approaches and indicators. This paper aims to analyse the extent to which the use of different methodologies leads to different results in terms of exposure to the risk of carbon leakage, using the Brazilian industry sector as a case study. Results indicate that carbon leakage exposure is an expected outcome of eventual CPI implementation in Brazilian industry. However, results vary according to the chosen methodology, so the definition of the criteria is paramount for assessing sectoral exposure to the risk of carbon leakage.

Key policy insights

• Despite increasing discussion about the implementation of carbon pricing on the Brazilian industrial sector, the evaluation of carbon leakage risks is still neglected.

• Assessments of the risk of carbon leakage are directly related to the indicators and criteria used by each methodology. Thus, a given subsector may present different levels of exposure to carbon leakage depending on the methodological choice.

• More than a purely technical discussion, the methodological definition of carbon leakage risk is a political discussion – it can be well-conducted, leading to the success of a CPI, or even sabotaged, by implicitly subsidizing energy-intensive industries.

     

Long‐term temporal and spatial patterns in bioeroding sponge distribution at the Abrolhos Bank,Brazil, Southwestern Atlantic

Bioeroding sponges belong to the most dominant bioeroders, significantly contributing to the erosion of coral reefs. Some species are tolerant or even benefit from environmental conditions such as ocean warming, acidification, and eutrophication. In consequence, increases in sponge bioerosion have been observed on some coral reefs over the last decades. The Abrolhos Bank is the largest coral reef system in the South Atlantic. It has been affected by sedimentation, eutrophication, overfishing, and climate change, mainly affecting coastal reefs, and at lesser intensity outer ones as well. This study aimed to describe spatial and temporal patterns in bioeroding sponge distribution in carbonate substrates in the Abrolhos Bank. Photo‐quadrats were used to compare bioeroding sponge abundance between two shallow reefs: a coastal, Pedra de Leste (PL), and an outer reef, Parcel dos Abrolhos (PAB). Each individual was delimitated over the substrate by determining the sponge surface through a line connecting the outermost papillae. The study was conducted over 6 years in 2008–2009 and 2013–2016. Four species of bioeroding sponges were identified: Cliona carteri Ridley, 1881, C. delitrix Pang, 1973, C. cf. schmidtii Ridley, 1881, and Siphonodictyon coralliphagum Rützler, 1971. The distribution and abundance of species varied between the inner and outer reefs and across the years, and displayed certain selectivity for the calcareous substrates recorded. Crustose coralline algae (CCA) were the main substrate excavated by the most abundant bioeroding species, C. carteri, and represented 70% of the substrate types occupied by this sponge (CCA, coral overgrown by CCA and plain coral). The highest abundance of bioeroding sponges observed in photo‐quadrats was 21.3 individuals/m² at the outer reefs (PAB) in 2014. The abundances or areal extents of bioeroding sponges were up to 10 times greater on the outer reefs than on the coastal ones, where sedimentation is higher and more strongly influenced by siliciclastic material. Moreover, a higher herbivorous fish biomass has been reported on outer reefs which could also influence the higher abundance of bioeroding sponges in outer reefs. During the study period of 6 years, an increase in bioeroding sponge abundance was observed at the outer reefs (PAB), with the sea surface temperature increase. As CCA have an important role in reefal cementation and carbonate production in the Abrolhos reefs, a bioerosion impact might be expected, in particular, on the outer reefs.     

Flexural deformation controls on Late Quaternary sediment dispersal in the Garo-Rajmahal Gap,NW Bengal Basin

Subsurface deformation is a driver for river path selection when deformation rates become comparable to the autogenic mobility rate of rivers. Here we combine geomorphology, soil and sediment facies analyses, and geophysical data of the Late Quaternary sediments of the central Garo-Rajmahal Gap in Northwest Bengal to link subsurface deformation with surface processes. We show variable sedimentation characteristics, from slow rates (<0.8 mm/year) in the Tista megafan at the foot of the Himalaya to nondeposition at the exposed surface of the Barind Tract to the south, enabling the development of mature soils. Combined subsidence in the Tista fan and uplift of the Barind Tract are consistent with a N-S flexural response of the Indian plate to loading of the Himalaya Mountains given a low value of elastic thickness (15–25 km). Provenance analysis based on bulk strontium concentration suggests a dispersal of sediment consistent with this flexural deformation—in particular the abandonment of the Barind Tract by a Pleistocene Brahmaputra River and the current extents of the Tista megafan lobes. Overall, these results highlight the control by deeply rooted deformation patterns on the routing of sediment by large rivers in foreland settings.     

Using Clustered Heat Maps in Mineral Exploration to Visualize Volcanic-Hosted Massive Sulfide Alteration and Mineralization

Natural Resources Research - This study proposes an extension of a visualization approach common in biochemistry (the clustered heat maps—CHMs) to geochemical data with the main objective of...     

Population connectivity of an overexploited coastal fish,Argyrosomus coronus (Sciaenidae), in an ocean-warming hotspot

The West Coast dusky kob Argyrosomus coronus is a commercially exploited fish with a distribution confined to the Angola–Benguela Frontal Zone (ABFZ) of the southeastern Atlantic Ocean. A previous study revealed that during a recent period of local warming the species extended its distribution into Namibian waters, where it hybridised with the resident and congeneric Argyrosomus inodorus. Environmental changes are a major threat to marine biodiversity and when combined with overfishing have the potential to accelerate the decline of species. However, little is known regarding the evolutionary history and population structure of A. coronus across the ABFZ. We investigated genetic diversity, population structure and historical demographic changes using mtDNA control region sequences and genotypes at six nuclear microsatellite loci, from 180 individuals. A single, genetically homogeneous population was indicated across the distributional range of A. coronus (φ_ST = 0.041, F_ST = 0.000, D = 0.000; p > 0.05). These findings imply that the oceanographic features within the ABFZ do not appear to significantly influence population connectivity in A. coronus, which simplifies management of the species. However, reconstruction of the demographic history points to a close link between the evolutionary history of A. coronus and the environmental characteristics of the ABFZ. This outcome suggests the species’ vulnerability to the rapid environmental changes being observed across this region, and highlights a pressing need for transboundary management to mitigate the impacts of climate change in this global hotspot of seawater temperature changes.     

Analytical study of the swing-by maneuver in an elliptical system

The objective of the present paper is to derive a set of analytical equations that describe a swing-by maneuver realized in a system of primaries that are in elliptical orbits. The goal is to calculate the variations of energy, velocity and angular momentum as a function of the usual basic parameters that describe the swing-by maneuver, as done before for the case of circular orbits. In elliptical orbits the velocity of the secondary body is no longer constant, as in the circular case, but it varies with the position of the secondary body in its orbit. As a consequence, the variations of energy, velocity and angular momentum become functions of the magnitude and the angle between the velocity vector of the secondary body and the line connecting the primaries. The “patched-conics” approach is used to obtain these equations. The configurations that result in maximum gains and losses of energy for the spacecraft are shown next, and a comparison between the results obtained using the analytical equations and numerical simulations are made to validate the method developed here.