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.     

Tectono-sedimentary evolution of Jurassic–Cretaceous diapiric structures: Miravete anticline,Maestrat Basin,Spain

Integration of extensive fieldwork, remote sensing mapping and 3D models from high-quality drone photographs relates tectonics and sedimentation to define the Jurassic–early Albian diapiric evolution of the N–S Miravete anticline, the NW-SE Castel de Cabra anticline and the NW-SE Cañada Vellida ridge in the Maestrat Basin (Iberian Ranges, Spain). The pre shortening diapiric structures are defined by well-exposed and unambiguous halokinetic geometries such as hooks and flaps, salt walls and collapse normal faults. These were developed on Triassic salt-bearing deposits, previously misinterpreted because they were hidden and overprinted by the Alpine shortening. The Miravete anticline grew during the Jurassic and Early Cretaceous and was rejuvenated during Cenozoic shortening. Its evolution is separated into four halokinetic stages, including the latest Alpine compression. Regionally, the well-exposed Castel de Cabra salt anticline and Cañada Vellida salt wall confirm the widespread Jurassic and Early Cretaceous diapiric evolution of the Maestrat Basin. The NE flank of the Cañada Vellida salt wall is characterized by hook patterns and by a 500-m-long thin Upper Jurassic carbonates defining an upturned flap, inferred as the roof of the salt wall before NE-directed salt extrusion. A regional E-W cross section through the Ababuj, Miravete and Cañada-Benatanduz anticlines shows typical geometries of salt-related rift basins, partly decoupled from basement faults. These structures could form a broader diapiric region still to be investigated. In this section, the Camarillas and Fortanete minibasins displayed well-developed bowl geometries at the onset of shortening. The most active period of diapiric growth in the Maestrat Basin occurred during the Early Cretaceous, which is also recorded in the Eastern Betics, Asturias and Basque-Cantabrian basins. This period coincides with the peak of eastward drift of the Iberian microplate, with speeds of 20 mm/year. The transtensional regime is interpreted to have played a role in diapiric development.     

Geological and geophysical studies of the Agoudal impact structure (Central High Atlas,Morocco): New evidence for crater size and age

Since the discovery of shatter cones (SCs) near the village of Agoudal (Morocco, Central High Atlas Mountains) in 2013, the absence of one or several associated circular structures led to speculation about the age of the impact event, the number, and the size of the impact crater or craters. Additional constraints on the crater size, age, and erosion rates are obtained here from geological, structural, and geophysical mapping and from cosmogenic nuclide data. Our geological maps of the Agoudal impact site at the scales of 1:30,000 (6 km²) and 1:15,000 (2.25 km²) include all known occurrences of SCs in target rocks, breccias, and vertical to overturned strata. Considering that strata surrounding the impact site are subhorizontal, we argue that disturbed strata are related to the impact event. Three types of breccias have been observed. Two of them (br1‐2 and br2) could be produced by erosion–sedimentation–consolidation processes, with no evidence for impact breccias, while breccia (br1) might be impact related. The most probable center of the structure is estimated at 31°59′13.73?N, 5°30′55.14?W using the concentric deviation method applied to the orientation of strata over the disturbed area. Despite the absence of a morphological expression, the ground magnetic and electromagnetic surveys reveal anomalies spatially associated with disturbed strata and SC occurrences. The geophysical data, the structural observations, and the area of occurrence of SCs in target rocks are all consistent with an original size of 1.4–4.2 km in diameter. Cosmogenic nuclide data (³⁶Cl) constrain the local erosion rates between 220 ± 22 m Ma^?1 and 430 ± 43 m Ma^?1. These erosion rates may remove the topographic expression of such a crater and its ejecta in a time period of about 0.3–1.9 Ma. This age is older than the Agoudal iron meteorite age (105 ± 40 kyr). This new age constraint excludes the possibility of a genetic relationship between the Agoudal iron meteorite fall and the formation of the Agoudal impact site. A chronolgy chart including the Atlas orogeny, the alternation of sedimentation and erosion periods, and the meteoritic impacts is presented based on all obtained and combined data.     

Barred Spirals,Double Galaxies and massive Halos1

It is shown that there is no significant dependence of the frequency with which barred spirals occur on environment. Therefore, the difference between normal and barred spirals must be largely due to initial conditions. The excess of early type galaxie (in particular ellipticals and SO's) among the statistically selected binary samples is probably caused by a contamination by cluster members and will therefore result in a systematic overestimate of the masses of binary galaxies. The large fraction of multiple first-ranked galaxies found in rich clusters of BAUTZ -MORGAN type I is surprising in view of tidal effects and dynamical friction, which should lead to the destruction of such massive binary systems on a time scale much shorter than the age of the Universe.     

不同光生物反应器培养坛紫菜 丝状体的比较研究

Vegetative cells of Porphyra haitanensis sporophytes were grown in a bubble-column and a draft-tube internal-loop airlift bioreactor at 20℃ with an illuminance of 40μmol photons m-2 s-1 for 20days. The draft-tube internal-loop airlift bioreactor yielded maximum cell density of 4850mg DCW (Dried Cell Weight)/L and higher growth rate of 0.094/d, compared to the bubble-column, 3950mg DCW/L and 0.086/d, respectively. The advantages were related to the more homogeneous fluid dynamic characteristics in the draft-tube internal-loop airlift bioreactor than the bubble-column bioreactor.     

Evaluation of the performance of HPLC–CHEMTAX analysis for determining phytoplankton biomass and composition in a turbid estuary (Schelde, Belgium)

In the upper Schelde estuary in 2002, phytoplankton biomass and community composition were studied using microscopic and pigment analyses. Chlorophyll a concentration was a good predictor of phytoplankton biomass estimated from cell counts and biovolume measurements. The phytoplankton carbon to chlorophyll a ratio, however, was often unrealistically low (<10). CHEMTAX was used to estimate the contribution of the major algal groups to total chlorophyll a. The dominant algal groups were diatoms and chlorophytes. While diatom equivalents in chlorophyll a predicted diatom biomass relatively well, chlorophyte equivalents in chlorophyll a were only weakly related to chlorophyte biomass. The pigment-based approach to study phytoplankton overestimated phytoplankton biomass in general and chlorophyte biomass in particular in late autumn and winter, when phytoplankton biomass was low. A possible explanation for this overestimation may be the presence of large amounts of vascular plant detritus in the upper Schelde estuary. Residual chlorophyll a, chlorophyll b and lutein in this detritus may result in an overestimation of total phytoplankton and chlorophyte biomass when the contribution of phytoplankton to total particulate organic matter is low.     

Inter-annual variability in phytoplankton summer blooms in the freshwater tidal reaches of the Schelde estuary (Belgium)

The inter-annual variability in phytoplankton summer blooms in the upper reaches of the Schelde estuary was investigated between 1996 and 2005 by monthly sampling at 10 stations. The large inter-annual variations of the chlorophyll a concentration in the freshwater tidal reaches were independent from variations in chlorophyll a in the tributary river Schelde. Summer mean chlorophyll a concentrations were significantly negatively correlated with flushing rate (Spearman correlation: r = −0.67, p = 0.05, n = 9) but not with temperature, irradiance and suspended particulate matter or dissolved silica (DSi) concentrations. During dry summers, low flushing rates permitted the development of dense phytoplankton populations in the upper part of the estuary, while during wet summers high flushing rates prevented the development of dense phytoplankton blooms. Flushing rate was also found to be important for the phytoplankton community composition. At low flushing rates, the community was dominated by diatoms that developed within the upper estuary. At high flushing rates, chlorophytes imported from the tributary river Schelde became more important in the phytoplankton community. The position of the chlorophyll a maximum shifted from the head of the estuary when flushing rates were low, to more downstream when flushing rates were high. Although DSi concentrations tended to be lower during years of high phytoplankton (mainly diatom) biomass, the relation with flushing rate was not significant.     

Kinetic and equilibrium studies of copper-dissolved organic matter complexation in water column of the stratified Krka River estuary (Croatia)

An interaction of dissolved natural organic matter (DNOM) with copper ions in the water column of the stratified Krka River estuary (Croatia) was studied. The experimental methodology was based on the differential pulse anodic stripping voltammetric (DPASV) determination of labile copper species by titrating the sample using increments of copper additions uniformly distributed on the logarithmic scale. A classical at-equilibrium approach (determination of copper complexing capacity, CuCC) and a kinetic approach (tracing of equilibrium reconstitution) of copper complexation were considered and compared. A model of discrete distribution of organic ligands forming inert copper complexes was applied. For both approaches, a home-written fitting program was used for the determination of apparent stability constants (K_i^equ), total ligands concentration (L_iT) and association/dissociation rate constants (k_i¹,k_i^- 1).A non-conservative behaviour of dissolved organic matter (DOC) and total copper concentration in a water column was registered. An enhanced biological activity at the freshwater–seawater interface (FSI) triggered an increase of total copper concentration and total ligand concentration in this water layer. The copper complexation in fresh water of Krka River was characterised by one type of binding ligands, while in most of the estuarine and marine samples two classes of ligands were identified. The distribution of apparent stability constants (log K₁^equ: 11.2–13.0, log K₂^equ:8.8–10.0) showed increasing trend towards higher salinities, indicating stronger copper complexation by autochthonous seawater organic matter.Copper complexation parameters (ligand concentrations and apparent stability constants) obtained by at-equilibrium model are in very good accordance with those of kinetic model. Calculated association rate constants (k₁¹:6.1–20 × 10³ (M s)^− 1, k₂¹: 1.3–6.3 × 10³ (M s)^− 1) indicate that copper complexation by DNOM takes place relatively slowly. The time needed to achieve a new pseudo-equilibrium induced by an increase of copper concentration (which is common for Krka River estuary during summer period due to the nautical traffic), is estimated to be from 2 to 4 h.It is found that in such oligotrophic environment (dissolved organic carbon content under 83 µM_C, i.e. 1 mg_CL^− 1) an increase of the total copper concentration above 12 nM could enhance a free copper concentration exceeding the level considered as potentially toxic for microorganisms (10 pM).     

High-resolution record of tectonic and sedimentary processes in growth strata

Growth strata are used to determine the kinematics of synsedimentary structures such as faults. Classical methods of analysis such as thickness versus throw plot consider that the available space created by fault slip in the hanging wall of faults is instantaneously filled up by sediments. This has lead many previous works to identify a cyclic activity for growth faults. Here we perform a careful analysis of the variation of strata thicknesses on each side of a very well documented normal growth fault in the Niger delta. We show that these thickness variations are induced by the alternation of sedimentary processes during continuous fault slip. Suspended-load processes induce either uniform or slightly variable thickness of a large majority of mudstone layers. Bedload processes result in a preferential thickening of sand layers in the hanging wall. These high quality data thus provide strong grounds for doubting the polycyclic growth diagnosed for some faults at the scale of sedimentary cycles and supports the notion that fault displacement rates can be very well behaved. Our study emphasizes the important conclusion that stable fault growth, and related displacement rates, can appear to be punctuated when viewed at the scale of sedimentary cycles. It follows that care should be taken when attempting to derive displacement rates on temporal scales equivalent to those of alternating sedimentological cycles.     

A recent history of science cases for optical interferometry

Optical long-baseline interferometry is a unique and powerful technique for astronomical research. Since the 1980’s (with I2T, GI2T, Mark I to III, SUSI, ...), optical interferometers have produced an increasing number of scientific papers covering various fields of astrophysics. As current interferometric facilities are reaching their maturity, we take the opportunity in this paper to summarize the conclusions of a few key meetings, workshops, and conferences dedicated to interferometry. We present the most persistent recommendations related to science cases and discuss some key technological developments required to address them. In the era of extremely large telescopes, optical long-baseline interferometers will remain crucial to probe the smallest spatial scales and make breakthrough discoveries.