A derivative-free methodology with local and global search for the constrained joint optimization of well locations and controls

In oil field development, the optimal location for a new well depends on how it is to be operated. Thus, it is generally suboptimal to treat the well location and well control optimization problems separately. Rather, they should be considered simultaneously as a joint problem. In this work, we present noninvasive, derivative-free, easily parallelizable procedures to solve this joint optimization problem. Specifically, we consider Particle Swarm Optimization (PSO), a global stochastic search algorithm; Mesh Adaptive Direct Search (MADS), a local search procedure; and a hybrid PSO–MADS technique that combines the advantages of both methods. Nonlinear constraints are handled through use of filter-based treatments that seek to minimize both the objective function and constraint violation. We also introduce a formulation to determine the optimal number of wells, in addition to their locations and controls, by associating a binary variable (drill/do not drill) with each well. Example cases of varying complexity, which include bound constraints, nonlinear constraints, and the determination of the number of wells, are presented. The PSO–MADS hybrid procedure is shown to consistently outperform both stand-alone PSO and MADS when solving the joint problem. The joint approach is also observed to provide superior performance relative to a sequential procedure.     

Water effects on rock strength and stiffness degradation

Reduction in strength and stiffness in rocks attributed to an increase in water content has been extensively researched on a large variety of rock types over the past decades. Due to the considerable variations of texture and lithology, the extent of water-weakening effect is highly varied among different rock types, spanning from nearly negligible in quartzite to 90 % of uniaxial compressive strength reduction in shale. Readers, however, often face difficulties in comparing the data published in different sources due to the discrepancy of experimental procedures of obtaining the water saturation state and how the raw laboratory data is interpreted. In view of this, the present paper first reviews the terminologies commonly used to quantify the amount of water stored in rocks. The second part of the paper reviews the water-weakening effects on rock strengths, particularly focusing on uniaxial compressive strength and modulus, as well as tensile strength, under quasi-static loading and dynamic loading. The correlation relationships established among various parameters, including porosity, density and fabric of rocks, and external factors such as strain rate, surface tension and dielectric constant of the saturating liquid, absorption percentage and suction pressure, are reviewed and presented toward the end of the paper.     

Phytoplankton photosynthetic response to light in an internal tide dominated environment

The influence of internal tides on phytoplankton photosynthetic response to light was studied in the Lower St. Lawrence Estuary. Photosynthesis at saturating light intensity responded to variations in the vertical density gradient, which were linked to the internal tides. The photosynthetic response was lag-correlated to the vertical water stratification. This suggests that the link between photosynthesis and the internal tides may have resulted from phytoplankton light adaptation.     

Textural Characterization Using Gy’s Heterogeneity Functions

A solid can be regarded as a set of contiguous elementary units. The distribution within the solid of any properties, measurable within each elementary unit, can be characterized using two parameters. These parameters are built using the constitution and distribution heterogeneities of Gy. The former account for the granularity of the elementary units, whereas the latter assess the spatial distribution of the property. A texture which definition involves several properties can be described using a diagram where both parameters work as variables. Potential applications encompass (i) the textural classification of soils, ore, breccia, and concrete and (ii) the monitoring of textural transformation during process like dolomitization, metamorphism, weathering, deformation, or annealing.     

The Anfeg post-collisional Pan-African high-K calc-alkaline batholith (Central Hoggar, Algeria), result of the LATEA microcontinent metacratonization

The Anfeg batholith (or composite laccolith) occupies a large surface (2000 km²) at the northern tip of the Laouni terrane, just south of Tamanrasset in Hoggar. It is granodioritic to granitic in composition and comprises abundant enclaves that are either mafic microgranular enclaves (MME) or gneissic xenoliths. It intruded an Eburnian (≈2 Ga) high-grade basement belonging to the LATEA metacraton at approximately 608 Ma (recalculated from the U–Pb dating of [Tectonics 5 (1986) 955]) and cooled at approximately 4 kbar, with a temperature of about 750 °C. This emplacement occurred mainly along subhorizontal thrust planes related to Pan-African subvertical mega-shear zones close to the attachment zone of a strike-slip partitioned transpression system. Although affected by some LILE mobility, the Anfeg batholith can be ascribed to a high-K calc-alkaline suite but characterized by low heavy REE contents and high LREE/HREE ratios. The MME belong to the Anfeg magmatic trend while some xenoliths belong to Neoproterozoic island arc rocks.The Anfeg batholith defines a Nd–Sr isotopic initial ratios trend (_Nd/(⁸⁷Sr/⁸⁶Sr)_i from −2.8/0.7068 to −11.8/0.7111) pointing to a mixing between a depleted mantle and an old Rb-depleted granulitic lower crust. Both sources have been identified within LATEA and elsewhere in the Tuareg shield (_Nd/⁸⁷Sr/⁸⁶Sr)_i of +6.2/0.7028 for the depleted mantle, −22/0.708 for the old lower crust.The model proposed relates the above geochemical features to a lithospheric delamination along the subvertical mega-shear zones that dissected the rigid LATEA former passive margin without major crustal thickening (metacratonization) during the general northward tectonic escape of the Tuareg terranes, a consequence of the collision with the West African craton. This delamination allowed the uprise of the asthenosphere. In turn, this induced the melting of the asthenosphere by adiabatic pressure release and of the old felsic and mafic lower crust due to the high heat flow. A gradient in the mantle/crust ratio within the source of the Pan-African magmatism is observed in LATEA from the northeast (Egéré-Aleksod terrane) where rare plutons are rooted within the Archaean/Eburnian basement to the southwest (Laouni terrane) where abundant batholiths, including Anfeg, have a mixed signature. Some mantle melts with only slight crustal contamination (Laouni troctolitic layered intrusions) are even present. This suggests that the southern boundary of LATEA microcontinent is not far south of the Tuareg shield.     

Acoustic monitoring of gas emissions from the seafloor. Part II: a case study from the Sea of Marmara

A rotating, acoustic gas bubble detector, BOB (Bubble OBservatory) module was deployed during two surveys, conducted in 2009 and 2011 respectively, to study the temporal variations of gas emissions from the Marmara seafloor, along the North Anatolian Fault zone. The echosounder mounted on the instrument insonifies an angular sector of 7° during a given duration (of about 1 h). Then it rotates to the next, near-by angular sector and so forth. When the full angular domain is insonified, the “pan and tilt system” rotates back to its initial position, in order to start a new cycle (of about 1 day). The acoustic data reveal that gas emission is not a steady process, with observed temporal variations ranging between a few minutes and 24 h (from one cycle to the other). Echo-integration and inversion performed on the acoustic data as described in the companion paper of Leblond et al. (Mar Geophys Res, 2014), also indicate important variations in, respectively, the target strength and the volumetric flow rates of individual sources. However, the observed temporal variations may not be related to the properties of the gas source only, but reflect possible variations in sea-bottom currents, which could deviate the bubble train towards the neighboring sector. During the 2011 survey, a 4-component ocean bottom seismometer (OBS) was co-located at the seafloor, 59 m away from the BOB module. The acoustic data from our rotating, monitoring system support, but do not provide undisputable evidence to confirm, the hypothesis formulated by Tary et al. (2012), that the short-duration, non-seismic micro-events recorded by the OBS are likely produced by gas-related processes within the near seabed sediments. Hence, the use of a multibeam echosounder, or of several split beam echosounders should be preferred to rotating systems, for future experiments.     

A porosity mapping survey in Hecate Strait using a seafloor electro-magnetic profiling system

A newly developed marine electromagnetic (E-M) system was used to create apparent porosity maps of the uppermost 20 m of sediments on the seafloor. Measurements of electrical conductivity were interpreted to give the porosity of the bottom sediments and underlying units. The data were collected continuously in a surveying mode using two receivers at different spacings to improve the resolution of porosity as a function of depth. The variations in apparent porosity over the area correlate well with information obtained from cores and acoustic profiles.

Acoustic profiles indicated the presence of biogenic gas in the central region of the survey area. The absence of distinct changes in apparent porosity over regions of gas allows us to conclude that the concentration of gas trapped within the Queen Charlotte muds does not exceed about 4% of the sediment by volume if we assume that the gas displaces fluid in the sediment.

A buried, lower porosity layer was detected, its position reflecting changes in the depth to the glaciomarine layer. Mapped variations in apparent porosity are well correlated with features on the acoustic record. This illustrates the system's ability to obtain continuous profiles of apparent porosity over seafloor features. The marine E-M survey provided rapid areal coverage, and, combined with the acoustic profiles, information on the porosity of deeper units.     

Distinctive C isotopic signature distinguishes a novel sea ice biomarker in Arctic sediments and sediment traps

A C₂₅ highly branched isoprenoid (HBI) monoene hydrocarbon, designated IP₂₅, has been proposed previously to originate from diatoms living in Arctic sea ice, while the presence of IP₂₅ in sediments has been suggested to be a proxy for the occurrence of former Arctic sea ice. Here, we show that the ¹³C isotopic composition of IP₂₅ in sea ice, in sediment trap material collected under sea ice, and in high latitude northern sediments, is distinctive (isotopically ‘heavy’) and distinguishable from that of organic matter of planktonic or terrigenous origin. Mean δ¹³C values for IP₂₅ were − 22.3 ± 0.4‰ (sea ice), − 19.6 ± 1.1‰ (sediment traps) and − 19.3 ± 2.3‰ (sediments). These measurements, therefore, support further the proposed use of IP₂₅ as an Arctic sea ice proxy.     

Origin and characteristics of the zooplankton phosphatase activity in a coastal ecosystem of the Mediterranean sea (Toulon Bay)

In Toulon Bay (France), very high phosphatase activities have been found in the zooplankton fraction>90 microm. This work was intended to specify their origin. For that purpose, larvae, juvenile and adult Crustacea (Copepods: Calanoids, Cyclopoids, Branchiopods: Cladocera, and Cirripeds) were isolated. Their activities were measured using paranitrophenyl phosphate dissolved in sea water in order to calculate Km (the enzyme half saturation concentration) and Vmax (the reaction rate when the enzyme is saturated with substrate). Vmax were referred to protein contents of the isolated organisms to calculate specific activities. For all zooplankton groups high and low affinity phosphatase activities were found. The low affinity enzyme was responsible for at least 70% of the total phosphatase activity. Its specific activity was higher for larvae than for copepodites and adults. In Cirriped nauplii this activity was particularly high with values which were several hundred times higher than that in other Crustacea. These enzymes had optimum pH close to 8.4, magnesium requirement and were competitively inhibited by orthophosphate. Experiments with intact and lysed Cirriped nauplii confirmed that living organisms had only a weak external activity and showed that most of the activity of these larvae was primarily intracellular.