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.     

Molybdenum Mass Fractions and Isotopic Compositions of International Geological Reference Materials

The double‐spike method with multi‐collector inductively coupled plasma‐mass spectrometry was used to measure the Mo mass fractions and isotopic compositions of a set of geological reference materials including the mineral molybdenite, seawater, coral, as well as igneous and sedimentary rocks. The long‐term reproducibility of the Mo isotopic measurements, based on two‐year analyses of NIST SRM 3134 reference solutions and seawater samples, was ≤ 0.07‰ (two standard deviations, 2s, n = 167) for δ^98/95Mo. Accuracy was evaluated by analyses of Atlantic seawater, which yielded a mean δ^98/95Mo of 2.03 ± 0.06‰ (2s, n = 30, relative to NIST SRM 3134 = 0‰) and mass fraction of 0.0104 ± 0.0006 μg g^?1 (2s, n = 30), which is indistinguishable from seawater samples taken world‐wide and measured in other laboratories. The comprehensive data set presented in this study serves as a reference for quality assurance and interlaboratory comparison of high‐precision Mo mass fractions and isotopic compositions.     

Why GPS makes distances bigger than they are

Global navigation satellite systems such as the Global Positioning System (GPS) is one of the most important sensors for movement analysis. GPS is widely used to record the trajectories of vehicles, animals and human beings. However, all GPS movement data are affected by both measurement and interpolation errors. In this article we show that measurement error causes a systematic bias in distances recorded with a GPS; the distance between two points recorded with a GPS is – on average – bigger than the true distance between these points. This systematic ‘overestimation of distance’ becomes relevant if the influence of interpolation error can be neglected, which in practice is the case for movement sampled at high frequencies. We provide a mathematical explanation of this phenomenon and illustrate that it functionally depends on the autocorrelation of GPS measurement error (C). We argue that C can be interpreted as a quality measure for movement data recorded with a GPS. If there is a strong autocorrelation between any two consecutive position estimates, they have very similar error. This error cancels out when average speed, distance or direction is calculated along the trajectory. Based on our theoretical findings we introduce a novel approach to determine C in real-world GPS movement data sampled at high frequencies. We apply our approach to pedestrian trajectories and car trajectories. We found that the measurement error in the data was strongly spatially and temporally autocorrelated and give a quality estimate of the data. Most importantly, our findings are not limited to GPS alone. The systematic bias and its implications are bound to occur in any movement data collected with absolute positioning if interpolation error can be neglected.     

Unraveling the Peruvian Phase of the Central Andes: stratigraphy,sedimentology and geochronology of the Salar de Atacama Basin (22°30–23°S), northern Chile

The Salar de Atacama Basin holds important information regarding the tectonic activity, sedimentary environments and their variations in northern Chile during Cretaceous times. About 4000 m of high‐resolution stratigraphic columns of the Tonel, Purilactis and Barros Arana Formations reveal braided fluvial and alluvial facies, typical of arid to semi‐arid environments, interrupted by scarce intervals with evaporitic, aeolian and lacustrine sedimentation, displaying an overall coarsening‐upward trend. Clast‐count and point‐count data evidence the progressive erosion from Mesozoic volcanic rocks to Palaeozoic basement granitoids and deposits located around the Cordillera de Domeyko area, which is indicative of an unroofing process. The palaeocurrent data show that the source area was located to the west. The U/Pb detrital zircon geochronological data give maximum depositional ages of 149 Ma for the base of the Tonel Formation (Agua Salada Member), and 107 Ma for its middle member (La Escalera Member); 79 Ma for the lower Purilactis Formation (Limón Verde Member), and 73 Ma for the Barros Arana Formation. The sources of these zircons were located mainly to the west, and comprised from the Coastal Cordillera to the Precordillera. The ages and pulses record the tectonic activity during the Peruvian Phase, which can be split into two large events; an early phase, around 107 Ma, showing uplift of the Coastal Cordillera area, and a late phase around 79 Ma indicating an eastward jump of the deformation front to the Cordillera de Domeyko area. The lack of internal deformation and the thicknesses measured suggest that deposition of the units occurred in the foredeep zone of an eastward‐verging basin. This sedimentation would have ended with the K‐T phase, recognized in most of northern Chile.     

Activated sludge acclimation for toluene and DEHP degradation in a two-phase partitioning bioreactor

The effect of activated sludge acclimation on the biodegradation of toluene in the presence of a biodegradable non-aqueous phase liquid, di (2-ethylhexyl) phthalate (DEHP), in a two-phase partitioning bioreactor was characterized. The influence of the presence of DEHP, at a ratio of 0.1 % (volume ratio), and of the acclimation of activated sludge (AS) on the biodegradation of hydrophobic VOC was studied. AS was acclimated to both toluene and DEHP simultaneously. Using acclimated cells, 73 and 96 % improvement of the mean biodegradation rates was recorded for toluene and the organic solvent (DEHP), respectively, if compared to the values recorded in the absence of acclimation, during tests performed in Erlenmeyer flasks. Degradation rates were further improved by the use of acclimated AS in a reactor with a large head space; degradation yields for toluene and DEHP were above 99 and 89 %, respectively.     

Nonparametric Estimation of Undrained Shear Strength for Normally Consolidated Clays

For decades, various well-known empirical correlations have been established and commonly used to estimate the undrained shear strength of normally consolidated clays in geotechnical practice. However, there are still many contradictory cases regarding the validity of these correlations, so their applicability and reliability are questionable. In this study, an extensive geotechnical database for normally consolidated clays has been compiled. Next, a Bayesian nonparametric method is exploited to identify the most critical soil index properties in inferring the undrained shear strength and a new predictive model is proposed. Finally, the performance of the proposed predictive model is assessed and compared with two well-known correlations using an independent testing database. It confirms that the proposed model not only possesses high predictability, but also provides superior performance over other existing empirical correlations.     

Identifying influence areas with connectivity analysis – application to the local perturbation of heterogeneity distribution for history matching

Numerical representations of a target reservoir can help to assess the potential of different development plans. To be as predictive as possible, these representations or models must reproduce the data (static, dynamic) collected on the field. However, constraining reservoir models to dynamic data – the history-matching process – can be very time consuming. Many uncertain parameters need to be taken into account, such as the spatial distribution of petrophysical properties. This distribution is mostly unknown and usually represented by millions of values populating the reservoir grid. Dedicated parameterization techniques make it possible to investigate many spatial distributions from a small number of parameters. The efficiency of the matching process can be improved from the perturbation of specific regions of the reservoir. Distinct approaches can be considered to define such regions. For instance, one can refer to streamlines. The leading idea is to identify areas that influence the production behavior where the data are poorly reproduced. Here, we propose alternative methods based on connectivity analysis to easily provide approximate influence areas for any fluid-flow simulation. The reservoir is viewed as a set of nodes connected by weighted links that characterize the distance between two nodes. The path between nodes (or grid blocks) with the lowest cumulative weight yields an approximate flow path used to define influence areas. The potential of the approach is demonstrated on the basis of 2D synthetic cases for the joint integration of production and 4D saturation data, considering several formulations for the weights attributed to the links.     

Pico- and Nanophytoplankton Dynamics in Two Coupled but Contrasting Coastal Bays in the NW Mediterranean Sea (France)

Due to its ecological context, the Toulon bay represents a site of scientific interest to study temporal plankton distribution, particularly pico- and nanophytoplankton dynamics. A monthly monitoring was performed during a two-year cycle (October 2013–December 2015) at two coupled sampling sites, referred to as Little and Large bays, which had different morphometric characteristics and human pressures. Flow cytometry analyses highlighted the fact that pico- and nanophytoplankton were more abundant in the eutrophic Little bay. Furthermore, it evidenced two community structures across the Toulon bays: at times, a co-dominance of picoeukaryotes, nanoeukaryotes, Synechococcus 1-like cells and Prochlorococcus-like cells was found, and at other times, a Synechococcus 1-like dominated community existed. The alternation of one structure or the other can be explained by a combined action of temperature regime, nutrient conditions and degree of contamination. This study showed that pico- and nanophytoplankton dynamics were mainly driven by temperature in both sites, as in other temperate Mediterranean regions. Thus, the community was mainly composed of picoeukaryotes and Prochlorococcus-like cells in the winter (<?15 °C), while it was dominated by Synechococcus 1-like cells in the summer (>?20 °C). Additionally, the multiple human stressors in the Little bay seemed to affect the increase in abundance of Synechococcus 1-like cells as they were preferentially observed in the Large bay.     

Hydrogeochemical Features of Thermal Waters of South Trungbo (Central Vietnam)

The results of studying the features of the hydrogeological structure and chemical and isotope composition of thermal waters from the central part of Vietnam that are characterized by intense manifestations of intrusive magmatism are presented. It is established that low–and high–thermal waters with temperature varying within 30–85°C are developed in the area under study. The value of total mineralization of the hydrotherms ranges from 0.05 to 10.05 g/dm³. It is assumed that the circulation of thermal waters that are different in temperature and chemical composition occurs at two levels. The regular change of the hydrotherm composition in the direction from mineralized chloride sodium, including with increased Ca content, to fresh sodium bicarbonate is revealed. The ratio of δ¹⁸O–δ²H isotopes indicates that the water component is based on meteoric water. In the coastal areas, there is an isotope shift towards the ocean waters, which is also confirmed by the hydrogeochemical data. The key factors for forming the chemical composition of the thermal waters in South Trungbo are their genetic type, the interaction processes in the “water–rock–gas–organic substance” system, and their equilibrium–nonequilibrium state.     

Estimation of commercial diamond grades based on microdiamonds: a case study of the Koidu diamond mine,Sierra Leone

This paper documents the application of a microdiamond-based approach to the estimation of diamond grade in the Pipe 1 kimberlite at the Koidu mine in Sierra Leone. A geological model of Pipe 1 was constructed to represent the distribution and volume of the dominant kimberlite units within the pipe. Bulk samples, along with representative microdiamond samples, were collected from these units at surface and were used to define the ratio between microdiamond stone frequency (+212 μm stones per kilogram) and recoverable macrodiamond grade (+1.2 mm carats per tonne; 1 carat = 0.2 g). These ratios were applied to a comprehensive, spatially representative microdiamond sample dataset and were combined with a spatial model of country-rock xenolith dilution within the pipe to estimate +1.2 mm recoverable grades. The resource estimate was reconciled with subsequent production results in the elevation range 160 to 100 m above sea level. Production results for each of the six 10 m benches covering this elevation range were compared to the estimated average grades for these zones in the pipe. For the five cases where most of the kimberlite mass on a given bench is represented in the production data, the results show a maximum discrepancy of 6% between predicted and reported production grade with no indication of any consistent bias. This indicates that, when supported by a sound geological model and suitable microdiamond and macrodiamond data, the microdiamond-based estimation approach can provide reliable constraints on macrodiamond grade, even in the case of geologically complex bodies such as Koidu Pipe 1.