Determination of lower and upper bounds of predicted production from history-matched models

We present a method to determine lower and upper bounds to the predicted production or any other economic objective from history-matched reservoir models. The method consists of two steps: 1) performing a traditional computer-assisted history match of a reservoir model with the objective to minimize the mismatch between predicted and observed production data through adjusting the grid block permeability values of the model. 2) performing two optimization exercises to minimize and maximize an economic objective over the remaining field life, for a fixed production strategy, by manipulating the same grid block permeabilities, however without significantly changing the mismatch obtained under step 1. This is accomplished through a hierarchical optimization procedure that limits the solution space of a secondary optimization problem to the (approximate) null space of the primary optimization problem. We applied this procedure to two different reservoir models. We performed a history match based on synthetic data, starting from a uniform prior and using a gradient-based minimization procedure. After history matching, minimization and maximization of the net present value (NPV), using a fixed control strategy, were executed as secondary optimization problems by changing the model parameters while staying close to the null space of the primary optimization problem. In other words, we optimized the secondary objective functions, while requiring that optimality of the primary objective (a good history match) was preserved. This method therefore provides a way to quantify the economic consequences of the well-known problem that history matching is a strongly ill-posed problem. We also investigated how this method can be used as a means to assess the cost-effectiveness of acquiring different data types to reduce the uncertainty in the expected NPV.     

The Influence of Bioturbation on Iron and Sulphur Cycling in Marine Sediments: A Model Analysis

The geochemical cycles of iron and sulphur in marine sediments are strongly intertwined and give rise to a complex network of redox and precipitation reactions. Bioturbation refers to all modes of transport of particles and solutes induced by larger organisms, and in the present-day seafloor, bioturbation is one of the most important factors controlling the biogeochemical cycling of iron and sulphur. To better understand how bioturbation controls Fe and S cycling, we developed reactive transport model of a coastal sediment impacted by faunal activity. Subsequently, we performed a model sensitivity analysis, separately investigating the two different transport modes of bioturbation, i.e. bio-mixing (solid particle transport) and bio-irrigation (enhanced solute transport). This analysis reveals that bio-mixing and bio-irrigation have distinct—and largely opposing effects on both the iron and sulphur cycles. Bio-mixing enhances transport between the oxic and suboxic zones, thus promoting the reduction of oxidised species (e.g. iron oxyhydroxides) and the oxidation of reduced species (e.g. iron sulphides). Through the re-oxidation of iron sulphides, bio-mixing strongly enhances the recycling of Fe and S between their reduced and oxidised states. Bio-irrigation on the other hand removes reduced solutes, i.e. ferrous iron and free sulphide, from the sediment pore water. These reduced species are then reoxidised in the overlying water and not recycled within the sediment column, which leads to a decrease in Fe and S recycling. Overall, our results demonstrate that the ecology of the macrofauna (inducing bio-mixing or bio-irrigation, or both) matters when assessing their impact on sediment geochemistry. This finding seems particularly relevant for sedimentary cycling across Cambrian transition, when benthic fauna started colonizing and reworking the seafloor.     

The sedimentary and tectonic evolution of the Amur River and North Sakhalin Basin: new evidence from seismic stratigraphy and Neogene–Recent sediment budgets

The North Sakhalin Basin in the western Sea of Okhotsk has been the main site of sedimentation from the Amur River since the Early Miocene. In this article, we present regional seismic reflection data and a Neogene–Recent sediment budget to constrain the evolution of the basin and its sedimentary fill, and consider the implications for sediment flux from the Amur River, in particular testing models of continental‐scale Neogene drainage capture. The Amur‐derived basin‐fill history can be divided into five distinct stages: the first Amur‐derived sediments (>21–16.5 Ma) were deposited during a period of transtension along the Sakhalin‐Hokkaido Shear Zone, with moderately high sediment flux to the basin (71 Mt year^?1). The second stage sequence (16.5–10.4 Ma) was deposited following the cessation of transtension, and was characterised by a significant reduction in sediment flux (24 Mt year^?1) and widespread retrogradation of deltaic sediments. The third (10.4–5.3 Ma) and fourth (5.3–2.5 Ma) stages were characterised by progradation of deltaic sediments and an associated increase in sediment flux (48–60 Mt year^?1) to the basin. Significant uplift associated with regional transpression started during this time in southeastern Sakhalin, but the north‐eastward propagating strain did not reach the NE shelf of Sakhalin until the Pleistocene (<2.5 Ma). This uplift event, still ongoing today, resulted in recycling of older deltaic sediments from the island of Sakhalin, and contributed to a substantially increased total sediment flux to the adjacent basinal areas (165 Mt year^?1). Adjusted rates to discount these local erosional products (117 Mt year^?1) imply an Amur catchment‐wide increase in denudation rates during the Late Pliocene–Pleistocene; however, this was likely a result of global climatic and eustatic effects, combined with tectonic processes within the Amur catchment and possibly a smaller drainage capture event by the Sungari tributary, rather than continental‐scale drainage capture involving the entire upper Amur catchment.     

Seasonally and spatially referenced estimates of recreational shore-angling effort,catch composition,catch rates and total catch in the Goukamma Marine Protected Area,South Africa

A roving creel survey of the recreational shore fishery along the 16.4-km coastline in the Goukamma Marine Protected Area on the south coast of South Africa was conducted from 2009 to 2011. Some 838 patrols were stratified equally among months, areas and years, but intentionally biased towards weekends. Angler densities at Buffalo Bay and Groenvlei were 0.59 and 0.28 anglers km?1, respectively. Weekend densities were double to quadruple weekday densities and fishing during winter was more popular than during summer. Area, habitat and distance to access points explained variation in angler densities. Shannon–Wiener diversity in catches declined from 2.18 in an earlier (1993–2002) survey to 1.79. Although the order of species abundance in the catches remained largely unchanged, blacktail Diplodus capensis dominance increased to 57.3% by number, at the expense of galjoen Dichistius capensis. Habitat explained 27% of the variance in catch composition. The catch per unit effort (CPUE) for the top nine species ranged from 0.19 to 6.35 fish 100-h?1. The CPUE of all species, except spotted grunter Pomadasys commersonnii, declined. Blacktail and galjoen CPUE declined by 17% and 77%, respectively. The total catch estimate was 2 986 fish y?1. Transgressions of size limits were common. The results suggest that the fishery is overexploited and that catch rates are declining.     

Carbon Isotope Characteristics as Evidence of an External Source of High-Temperature Granitoids in Granulite Complexes

Doklady Earth Sciences - Data on the carbon isotope composition of graphite and CO2 from inclusions in quartz of granitoids of the Southern Marginal Zone (SMZ) of the Limpopo granulite belt, South...     

A distributed parallel multiple-relaxation-time lattice Boltzmann method on general-purpose graphics processing units for the rapid and scalable computation of absolute permeability from high-resolution 3D micro-CT images

Digital rock physics (DRP) is a rapidly evolving technology targeting fast turnaround times for repeatable core analysis and multi-physics simulation of rock properties. We develop and validate a rapid and scalable distributed-parallel single-phase pore-scale flow simulator for permeability estimation on real 3D pore-scale micro-CT images using a novel variant of the lattice Boltzmann method (LBM). The LBM code implementation is designed to take maximum advantage of distributed computing on multiple general-purpose graphics processing units (GPGPUs). We describe and extensively test the distributed parallel implementation of an innovative LBM algorithm for simulating flow in pore-scale media based on the multiple-relaxation-time (MRT) model that utilizes a precise treatment of body force. While the individual components of the resulting simulator can be separately found in various references, our novel contributions are (1) the integration of all of the mathematical and high-performance computing components together with a highly optimized code implementation and (2) the delivery of quantitative results with the simulator in terms of robustness, accuracy, and computational efficiency for a variety of flow geometries including various types of real rock images. We report on extensive validations of the simulator in terms of accuracy and provide near-ideal distributed parallel scalability results on large pore-scale image volumes that were largely computationally inaccessible prior to our implementation. We validate the accuracy of the MRT-LBM simulator on model geometries with analytical solutions. Permeability estimation results are then provided on large 3D binary microstructures including a sphere pack and rocks from various sandstone and carbonate formations. We quantify the scalability behavior of the distributed parallel implementation of MRT-LBM as a function of model type/size and the number of utilized GPGPUs for a panoply of permeability estimation problems.     

Precursory slope distress prior to the 2010 Mount Meager landslide,British Columbia

In 2010, the south flank of Mount Meager failed catastrophically, generating the largest (53 ± 3.8 × 10⁶ m³) landslide in Canadian history. We document the slow deformation of the edifice prior to failure using archival historic aerial photographs spanning the period 1948–2006. All photos were processed using Structure from Motion (SfM) photogrammetry. We used the SfM products to produce pre-and post-failure geomorphic maps that document changes in the volcanic edifice and Capricorn Glacier at its base. The photographic dataset shows that the Capricorn Glacier re-advanced from a retracted position in the 1980s then rapidly retreated in the lead-up to the 2010 failure. The dataset also documents 60 years of progressive development of faults, toe bulging, and precursory failures in 1998 and 2009. The 2010 collapse was conditioned by glacial retreat and triggered by hot summer weather that caused ice and snow to melt. Meltwater increased pore water pressures in colluvium and fractured rocks at the base of the slope, causing those materials to mobilize, which in turn triggered several secondary failures structurally controlled by lithology and faults. The landslide retrogressed from the base of the slope to near the peak of Mount Meager involving basement rock and the overlying volcanic sequence. Elsewhere on the flanks of Mount Meager, large fractures have developed in recently deglaciated areas, conditioning these slopes for future collapse. Potential failures in these areas have larger volumes than the 2010 landslide. Anticipated atmospheric warming over the next several decades will cause further loss of snow and glacier ice, likely producing additional slope instability. Satellite- and ground-based monitoring of these slopes can provide advanced warning of future landslides to help reduce risk in populated regions downstream.     

Small shards and long distances — three cryptotephra layers from the Nahe palaeolake including the first discovery of Laacher See Tephra in Schleswig-Holstein (Germany)

Investigations of Lateglacial to Early Holocene lake sediments from the Nahe palaeolake (northern Germany) provided a high-resolution palynological record. To increase the temporal resolution of the record a targeted search for cryptotephra was carried out on the basis of pollen stratigraphy. Three cryptotephra horizons were detected and geochemically identified as G10ka series tephra (a Saksunarvatn Ash), Vedde Ash and Laacher See Tephra. Here we present the first geochemically confirmed finding of the ash from the Laacher See Eruption in Schleswig-Holstein—extending the so far detected fallout fan of the eruption further to the north-west. These finds enable direct stratigraphical correlations and underline the potential of the site for further investigations.     

Correction: Hydraulic tomography analysis of municipal-well operation data with geology-based groundwater models

Hydrogeology Journal - A Correction to this paper has been published: https://doi.org/10.1007/s10040-021-02351-x