Optimization of nonconventional wells under uncertainty using statistical proxies

The determination of the optimal type and placement of a nonconventional well in a heterogeneous reservoir represents a challenging optimization problem. This determination is significantly more complicated if uncertainty in the reservoir geology is included in the optimization. In this study, a genetic algorithm is applied to optimize the deployment of nonconventional wells. Geological uncertainty is accounted for by optimizing over multiple reservoir models (realizations) subject to a prescribed risk attitude. To reduce the excessive computational requirements of the base method, a new statistical proxy (which provides fast estimates of the objective function) based on cluster analysis is introduced into the optimization process. This proxy provides an estimate of the cumulative distribution function (CDF) of the scenario performance, which enables the quantification of proxy uncertainty. Knowledge of the proxy-based performance estimate in conjunction with the proxy CDF enables the systematic selection of the most appropriate scenarios for full simulation. Application of the overall method for the optimization of monobore and dual-lateral well placement demonstrates the performance of the hybrid optimization procedure. Specifically, it is shown that by simulating only 10% or 20% of the scenarios (as determined by application of the proxy), optimization results very close to those achieved by simulating all cases are obtained.     

Timing of deformation in the Norseman-Wiluna Belt, Yilgarn Craton, Western Australia

Establishing relative and absolute time frameworks for the sedimentary, magmatic, tectonic and gold mineralisation events in the Norseman-Wiluna Belt of the Archean Yilgarn Craton of Western Australia, has long been the main aim of research efforts. Recently published constraints on the timing of sedimentation and absolute granite ages have emphasized the shortcomings of the established rationale used for interpreting the timing of deformation events. In this paper the assumptions underlying this rationale are scrutinized, and it is shown that they are the source of significant misinterpretations. A revised time chart for the deformation events of the belt is established. The first shortening phase to affect the belt, D₁, was preceded by an extensional event D_1e and accompanied by a change from volcanic-dominated to plutonic-dominated magmatism at approximately 2685–2675 Ma. Later extension (D_2e) controlled deposition of the ca 2655 Ma Kurrawang Sequence and was followed by D₂, a major shortening event, which folded this sequence. D₂ must therefore have started after 2655 Ma—at least 20 Ma later than previously thought and after the voluminous 2670–2655 Ma high-Ca granite intrusion. Younger transcurrent deformation, D₃–D₄, waned at around 2630 Ma, suggesting that the crustal shortening deformation cycle D₂–D₄ lasted approximately 20–30 Ma, contemporaneous with low-volume 2650–2630 Ma low-Ca granites and alkaline intrusions. Time constraints on gold deposits suggest a late mineralisation event between 2640–2630 Ma. Thus, D₂–D₄ deformation cycle and late felsic magmatism define a 20–30 Ma long tectonothermal event, which culminated with gold mineralisation. The finding that D₂ folding took place after voluminous high-Ca granite intrusion led to research into the role of competent bodies during folding by means of numerical models. Results suggest that buoyancy-driven doming of pre-tectonic competent bodies trigger growth of antiforms, whereas non-buoyant, competent granite bodies trigger growth of synforms. The conspicuous presence of pre-folding granites in the cores of anticlines may be a result from active buoyancy doming during folding.     

Controls on Quaternary coastal evolution of the east-northeastern coast of Brazil: roles of sea-level history, trade winds and climate

East-northeastern Brazil has a wave-dominated, micro- to meso-tidal coast, lying entirely within the southern Atlantic trade wind belt. Integration of geologic mapping, radiocarbon dating and vibracoring data shows that the Quaternary coastal evolution of this area was controlled by three major factors: (1) sea-level history; (2) trade winds; and (3) climate change.

Sea-level history. Along the east-northeastern coast of Brazil, relative sea level has fallen approximately 5 m during the last 5000 y. Correlation of this sea-level history with the evolution of beach-ridge, lagoonal and coastal plain deposits shows that: (1) sea-level rise favours the formation of barrier island—lagoonal systems and the construction of intralagoonal deltas; (2) sea-level lowering is not conductive to barrier island formation. Rather, lagoons and bays become emergent and beach-ridge plains rapidly prograde.

Trade winds. Sediment dispersal systems along the coastal zone of east-northeastern Brazil have been highly persistent since Pleistocene time, as deduced from beach-ridge orientation. This persistence results from the fact that sediment dispersal in wave-dominated settings is ultimately controlled by atmospheric circulation which, for the east-northeastern coast of Brazil is associated with the South Atlantic high-pressure cell. The remarkable stability of this cell through time, has allowed the accumulation of extensive beach-ridge plains at the longshore drift sinks located along the coast.

Climate change. Effects of Quaternary climate changes on coastal sedimentation are twofold. Climate changes may affect rainfall patterns, thus exerting an important control on coastal dune development. Along the coast of northeastern Brazil, active coastal dunes are only present in those areas in which at least four consecutive dry months occur during the year. Mapping of these areas has shown that dune development during the Holocene has been episodic, these episodes being probably controlled by variations in rainfall patterns associated with climate changes. Secondly, despite its overall stability, the position of the high-pressure cell has experienced small shifts in position during the Holocene in response to climate changes. Changes in wind direction associated with these shifts have induced modifications in the coastal dispersion system, which are recorded in the strandplains as small truncations in the beach-ridge alignments.

These results have important implications in understanding accumulation of ancient sandstone shoreline sequences.     

Quantifying strain birefringence halos around inclusions in diamond

The pressure and temperature conditions of formation of natural diamond can be estimated by measuring the residual stress that an inclusion remains under within a diamond. Raman spectroscopy has been the most commonly used technique for determining this stress by utilising pressure-sensitive peak shifts in the Raman spectrum of both the inclusion and the diamond host. Here, we present a new approach to measure the residual stress using quantitative analysis of the birefringence induced in the diamond. As the analysis of stress-induced birefringence is very different from that of normal birefringence, an analytical model is developed that relates the spherical inclusion size, R _i, host diamond thickness, L, and measured value of birefringence at the edge of the inclusion, \Updelta n(R_\texti )_\textav \Updelta n(R_{\text{i}} )_{\text{av}} , to the peak value of birefringence that has been encountered; to first order \Updelta n_\textpk = (3/4)(L/R_\texti )  \Updelta n(R_\texti )_\textav \Updelta n_{\text{pk}} = (3/4)(L/R_{\text{i}} ) \, \Updelta n(R_{\text{i}} )_{\text{av}} . From this birefringence, the remnant pressure (P _i) can be calculated using the photoelastic relationship \Updelta n_\textpk = - (3/4)n³ q_\textiso P_\texti \Updelta n_{\text{pk}} = - (3/4)n^{3} q_{\text{iso}} P_{\text{i}} , where q _iso is a piezo-optical coefficient, which can be assumed to be independent of crystallographic orientation, and n is the refractive index of the diamond. This model has been used in combination with quantitative birefringence analysis with a MetriPol system and compared to the results from both Raman point and 2D mapping analysis for a garnet inclusion in a diamond from the Udachnaya mine (Russia) and coesite inclusions in a diamond from the Finsch mine (South Africa). The birefringence model and analysis gave a remnant pressure of 0.53 ± 0.01 GPa for the garnet inclusion, from which a source pressure was calculated as 5.7 GPa at 1,175°C (temperature obtained from IR analysis of the diamond host). The Raman techniques could not be applied quantitatively to this sample to support the birefringence model; they were, however, applied to the largest coesite inclusion in the Finsch sample. The remnant pressure values obtained were 2.5 ± 0.1 GPa (birefringence), 2.5 ± 0.3 GPa (2D Raman map), and 2.5–2.6 GPa (Raman point analysis from all four inclusions). However, although the remnant pressures from the three methods were self-consistent, they led to anomalously low source pressure of 2.9 GPa at 1,150°C (temperature obtained from IR analysis) raising serious concerns about the use of the coesite-in-diamond geobarometer.     

Three dimensional numerical simulation of residential building on shrink–swell soils in response to climatic conditions

Shrink–swell soils can cause distresses in buildings, and every year, the economic loss associated with this problem is huge. This paper presents a comprehensive system for simulating the soil–foundation–building system and its response to daily weather conditions. Weather data include rainfall, solar radiation, air temperature, relative humidity, and wind speed, all of which are readily available from a local weather station or the Internet. These data are used to determine simulation flux boundary conditions. Different methods are proposed to simulate different boundary conditions: bare soil, trees, and vegetation. A coupled hydro‐mechanical stress analysis is used to simulate the volume change of shrink–swell soils due to both mechanical stress and water content variations. Coupled hydro‐mechanical stress‐jointed elements are used to simulate the interaction between the soil and the slab, and general shell elements are used to simulate structural behavior. All the models are combined into one finite element program to predict the entire system's behavior. This paper first described the theory for the simulations. A site in Arlington, Texas, is then selected to demonstrate the application of the proposed system. Simulation results are shown, and a comparison between measured and predicted movements for four footings in Arlington, Texas, over a 2‐year period is presented. Finally, a three‐dimensional simulation is made for a virtual residential building on shrink–swell soils to identify the influence of various factors. Copyright © 2015 John Wiley & Sons, Ltd.     

Impacts of hydrous manganese oxide on the retention and lability of dissolved organic matter

Minerals constitute a primary ecosystem control on organic C decomposition in soils, and therefore on greenhouse gas fluxes to the atmosphere. Secondary minerals, in particular, Fe and Al (oxyhydr)oxides—collectively referred to as “oxides” hereafter—are prominent protectors of organic C against microbial decomposition through sorption and complexation reactions. However, the impacts of Mn oxides on organic C retention and lability in soils are poorly understood. Here we show that hydrous Mn oxide (HMO), a poorly crystalline δ-MnO₂, has a greater maximum sorption capacity for dissolved organic matter (DOM) derived from a deciduous forest composite O_i, O_e, and O_a horizon leachate (“O horizon leachate” hereafter) than does goethite under acidic (pH 5) conditions. Nonetheless, goethite has a stronger sorption capacity for DOM at low initial C:(Mn or Fe) molar ratios compared to HMO, probably due to ligand exchange with carboxylate groups as revealed by attenuated total reflectance-Fourier transform infrared spectroscopy. X-ray photoelectron spectroscopy and scanning transmission X-ray microscopy–near-edge X-ray absorption fine structure spectroscopy coupled with Mn mass balance calculations reveal that DOM sorption onto HMO induces partial Mn reductive dissolution and Mn reduction of the residual HMO. X-ray photoelectron spectroscopy further shows increasing Mn(II) concentrations are correlated with increasing oxidized C (C=O) content (r = 0.78, P < 0.0006) on the DOM–HMO complexes. We posit that DOM is the more probable reductant of HMO, as Mn(II)-induced HMO dissolution does not alter the Mn speciation of the residual HMO at pH 5. At a lower C loading (2 × 10² μg C m^?2), DOM desorption—assessed by 0.1 M NaH₂PO₄ extraction—is lower for HMO than for goethite, whereas the extent of desorption is the same at a higher C loading (4 × 10² μg C m^?2). No significant differences are observed in the impacts of HMO and goethite on the biodegradability of the DOM remaining in solution after DOM sorption reaches steady state. Overall, HMO shows a relatively strong capacity to sorb DOM and resist phosphate-induced desorption, but DOM–HMO complexes may be more vulnerable to reductive dissolution than DOM–goethite complexes.     

Greatlakean Substage: A replacement for Valderan Substage in the Lake Michigan basin

New evidence from recent field and seismic investigations in the Lake Michigan basin and in the type areas of the Valders, Two Creeks and Two Rivers deposits necessitates revision of late-glacial ice-front positions, rock- and time-stratigraphic nomenclature and climatic interpretations and deglaciation patterns for the period ca. 14,000–7,000 radiocarbon years B.P. The previously reported and long accepted pattern of deglaciation for the Lake Michigan basin started with a regular retreat from the Lake Border Morainic System, with a minor oscillation marked by the Port Huron moraine(s) and then an extensive Twocreekan deglaciation followed by a major (320 km) post-Twocreekan advance (Valders). However, we now record a major retreat between the times of the Lake Border and Port Huron moraines, followed by a gradual retreat from the Port Huron limit and interrupted by a minor standstill (deposition of Manitowoc Till), a retreat (Twocreekan) and a readvance (Two Rivers Till). No Woodfordian or younger readvance was as extensive as had been the preceding one. This sequence argues for a normal, climatically controlled, progressive deglaciation rather than one interrupted by a major post-Twocreekan (formerly Valderan) surge. This revision appears finally to harmonize the geologic evidence and the palynological record for the Great Lakes region. Our investigations show that Valders Till from which the Valderan Substage was named is late-Woodfordian in age. We propose the term “Greatlakean” as a replacement for the now misleading time-stratigraphic term “Valderan”. The type section and the definition of the upper and lower boundaries of the Greatlakean Substage remain the same as those originally proposed for the Valderan Substage but the name is changed.     

Neogene upper‐crustal cooling of the Olympus range (northern Aegean): major role of Hellenic back‐arc extension over propagation of the North Anatolia Fault Zone

The North Anatolian Fault Zone (NAFZ) is one of the most hazardous active faults on Earth, yet its Pliocene space‐time propagation across the north Aegean domain remains poorly constrained. We use low‐temperature multi‐thermochronology and inverse thermal modelling to quantify the cooling history of the upper crust across the Olympus range. This range is located in the footwall of a system of normal faults traditionally interpreted as resulting from superposed Middle–Late Miocene N–S stretching, related to the back‐arc extension of the Hellenic subduction zone, and a Pliocene‐Quaternary transtensional field, attributed to the south‐westward propagation of the NAFZ. We find that accelerated exhumational cooling occurred between 12 and 6 Ma at rates of 15–35 °C Ma^?1 and decreased to <3 °C Ma^?1 by 8–6 Ma. The absence of significant Plio‐Pleistocene cooling across Olympus suggests that crustal exhumation there is driven by late Miocene back‐arc extension, while the impact of the NAFZ remains limited.     

Adjoint formulation and constraint handling for gradient-based optimization of compositional reservoir flow

An adjoint formulation for the gradient-based optimization of oil–gas compositional reservoir simulation problems is presented. The method is implemented within an automatic differentiation-based compositional flow simulator (Stanford’s Automatic Differentiation-based General Purpose Research Simulator, AD-GPRS). The development of adjoint procedures for general compositional problems is much more challenging than for oil–water problems due to the increased complexity of the code and the underlying physics. The treatment of nonlinear constraints, an example of which is a maximum gas rate specification in injection or production wells, when the control variables are well bottom-hole pressures, poses a particular challenge. Two approaches for handling these constraints are presented—a formal treatment within the optimizer and a simpler heuristic treatment in the forward model. The relationship between discrete and continuous adjoint formulations is also elucidated. Results for four example cases of increasing complexity are presented. Improvements in the objective function (cumulative oil produced) relative to reference solutions range from 4.2 to 11.6 %. The heuristic treatment of nonlinear constraints is shown to offer a cost-effective means for obtaining feasible solutions, which are, in some cases, better than those obtained using the formal constraint handling procedure.     

Les Facies de Transgression du Trias sur la Bordure Nord-est du Massif Central Français

Résumé La transgression triasique est suivie d'Est en Ouest sur la bordure du Massif Central français. Les formations sont de plus en plus récentes vers l'Ouest, indépendamment de leur faciès. Les limites du Buntsandstein et du Trias supérieur semblent locales et d'origine tectonique, tandis que la limite du Ladinien supérieur parait plus générale et correspond à la période où le caractère marin de la transgression sur les massifs hercyniens était le plus franc.Le début de la transgression aniso-ladinienne est daté pour l'Anisien par des empreintes de Reptiles et pour les faciès les plus marins du Ladinien par des Foraminifères et des pollens. A cette période, des microfaunes mésogéennes de mer ouverte apparaissent sur les bordures.Dans les faciès de transgression réduits, le Trias est constitué de grès fins silteux à pseudomorphoses de sel gemme et empreintes de Reptiles, étalés par une mer pelliculaire sur une surface plate. Le socle sous-jacent a fourni peu de matériel clastique et les courants sont orientés du Nord au Sud. Les évaporites sont plus concentrées sur la bordure que plus au large, ce qui implique une certaine aridité pendant toute la période de la transgression.Les faciès de bordure les plus avancés sur le socle sont formés d'une dolomie paralique supratidale déposée dans des lagunes évaporitiques. Sur les massifs hercyniens, des fractures et des failles ouvertes ont permis des circulations auxquelles sont liés des remplissages carbonates et des minéralisations.

---

Conditions of the triassic transgression are analysed from East to West across the eastern border of the Central Massif. Formations directly overlying the basement get younger and younger westwards, without relation with the environment. Buntsandstein and upper Triassic boundaries appear to be syntectonic local limits. The upper Ladininian boundary seems to be more general and connected with the main period of the marine environment for the transgression on the variscan basement rocks. The anisoladinian transgression beginns during the anisian age as indicated by fooprints, while ladinian age is supported by mesogean foraminifera and palynology during the climax of the transgression, in an open sea environment.Silty fine sandstones with footprints and salt casts constitute the triassic transgression facies on the flattened border of the Central Massif. The underlying basement supplied only a little arenaceous material and the currents are oriented from Nord to South. Evaporites are more concentrated at the boundary of the transgression than towards the basin, what means aridity.Paralic supratidal dolomite constitutes the most external deposits at the transgression limit, in a syntectonic environment. Dolomite infillings and ore deposits are connected with water circulations along faults and fractures.

---

Zusammenfassung Die Triastransgressionsschichten werden von Osten nach Westen am Nordrand des französischen Zentralmassivs beschrieben. Die auf dem Grundgebirge auflagernden Basisschichten werden, unabhängig von ihrer Faziesausbildung, nach Westen stets jünger. Die Begrenzungen der Ablagerungen des Buntsandsteins und der oberen Trias werden nur lokal, wahrscheinlich in Abhängigkeit von synsedimentären, tektonischen Strukturen angetroffen. Die Beckenausdehnung am Top des oberen Ladiniums ist dagegen zusammenhängend zu erkennen. Sie entspricht den eindeutigeren marinen Transgressionsverhältnissen auf dem variskischen Sockel.Die einsetzende Transgression des Aniso-Ladiniums wird durch Reptilienfährten als anisisch und die stärker marinen Ablagerungen durch Foraminiferen und Pollen als ladinisch datiert. In diesem Zeitabschnitt dringen offenmarine, mesogäische Mikrofossilien bis zum Beckenrand vor.Die geringmächtigen Transgressionsablagerungen bestehen aus siltigen Feinsanden mit Steinsalzpseudomorphosen und Reptilienfährten. Diese Sande wurden von einem sehr seichten Meer auf dem stark eingeebneten Grundgebirge verteilt. Das Grundgebirge hat nur sehr wenige klastische Abtragungsprodukte geliefert. Die Strömungen waren von Norden nach Süden ausgerichtet. Evaporite treten stärker in den randlichen Gebieten als zum offenen Becken hin auf. Dies zeigt auf eine ausgeprägte Aridität während der gesamten Transgressionszeit hin.Die am weitesten zum Festland hin vorgeschobenen Transgressionssedimente bestehen aus paralischen, subpratidalen Dolomiten, die sich in evaporitischen Lagunen bildeten. Die in den offenen Spalten und Verwerfungen des Grundgebirges zirkulierenden Lösungen führten zu Dolomitausfällungen und zu schichtgebundenen Vererzungen.

---

. , , , . , , , . . . - , — . , , . . , . . , . . , , . , , .