Ensemble clustering for efficient robust optimization of naturally fractured reservoirs

In the development of naturally fractured reservoirs (NFRs), the existence of natural fractures induces severe fingering and breakthrough. To manage the flooding process and improve the ultimate recovery, we propose a numerical workflow to generate optimal production schedules for smart wells, in which the inflow control valve (ICV) settings can be controlled individually. To properly consider the uncertainty introduced by randomly distributed natural fractures, the robust optimization would require a large ensemble size and it would be computationally demanding. In this work, a hierarchical clustering method is proposed to select representative models for the robust optimization in order to avoid redundant simulation runs and improve the efficiency of the robust optimization. By reducing the full ensemble of models into a small subset ensemble, the efficiency of the robust optimization algorithm is significantly improved. The robust optimization is performed using the StoSAG scheme to find the optimal well controls that maximize the net-present-value (NPV) of the NFR’s development. Due to the discrete property of a natural fracture field, traditional feature extraction methods such as model-parameter-based clustering may not be directly applicable. Therefore, two different kinds of clustering-based optimization methods, a state-based (e.g., s _w profiles) clustering and a response-based (e.g., production rates) clustering, are proposed and compared. The computational results show that the robust clustering optimization could increase the computational efficiency significantly without sacrificing much expected NPV of the robust optimization. Moreover, the performance of different clustering algorithms varies widely in correspondence to different selections of clustering features. By properly extracting model features, the clustered subset could adequately represent the uncertainty of the full ensemble.     

The Middle Triassic evolution of the Bangong–Nujiang Tethyan Ocean: evidence from analyses of OIB-type basalts and OIB-derived phonolites in northern Tibet

International Journal of Earth Sciences - In this paper, we present new major and trace element chemical data for the basalts and phonolites of the Nare ocean island fragment (NaOI), as well as...     

Tectonic significance of Triassic mafic rocks in the June Complex,Sanandaj–Sirjan zone,Iran

This study concentrates on the petrological and geochemical investigation of mafic rocks embedded within the voluminous Triassic June Complex of the central Sanandaj–Sirjan zone (Iran), which are crucial to reconstruct the geodynamics of the Neotethyan passive margin. The Triassic mafic rocks are alkaline to sub-alkaline basalts, containing 43.36–49.09 wt% SiO₂, 5.19–20.61 wt% MgO and 0.66–4.59 wt% total alkalis. Based on MgO concentrations, the mafic rocks fall into two groups: cumulates (Mg# = 51.61–58.94) and isotropic basaltic liquids (Mg# = 24.54–42.66). In all samples, the chondrite-normalized REE patterns show enrichment of light REEs with variable (La/Yb)_N ratios ranging from 2.48 to 9.00, which confirm their amalgamated OIB-like and E-MORB-like signatures. Enrichment in large-ion lithophile elements and depletion in high field strength elements (HFSE) relative to the primitive mantle further support this interpretation. No samples point to crustal contamination, all having undergone fractionation of olivine + clinopyroxene + plagioclase. Nevertheless, elemental data suggest that the substantial variations in (La/Sm)_PM and Zr/Nb ratios can be explained by variable degrees of partial melting rather than fractional crystallization from a common parental magma. The high (Nb/Yb)_PM ratio in the alkaline mafic rocks points to the mixing of magmas from enriched and depleted mantle sources. Abundant OIB alkaline basalts and rare E-MORB appear to be linked to the drifting stage on the northern passive margin of the Neotethys Ocean.     

Comparative study for environmental performances of traditional manufacturing and directed energy deposition processes

Additive manufacturing is considered more sustainable than traditional manufacturing due to its efficient energy and materials usage. However, previous literature indicates that this suggestion is applicable only for the polymer materials, and the environmental issues of additive manufacturing with metallic materials are still not clear. With the method of life cycle assessment, this paper analyzes and compares the energy consumptions and environmental impacts of direct energy deposition and traditional machining processes for a typical metal part. Further, the article attempts to identify the significant issues in the two manufacturing options that contribute most to the environmental impacts. Six environmental impacts were assessed in this study: global warming potential (GWP); acidification potential (AP); eutrophication potential; ozone depletion potential (ODP); photochemical ozone creation potential (POCP); and abiotic depletion potential (ADP). The results show that the gear laser fabrication process consumes more energy and releases more negative emissions compared with traditional gear manufacturing processes. The results of GWP, AP, ODP, ADP and POCP of the traditional gear manufacturing are only 30.33, 43.42, 17, 65.05 and 54.68% of the gear laser fabrication.     

Photodegradation of dimethyl phthalate (DMP) by UV–TiO<Subscript>2</Subscript> in aqueous solution: operational parameters and kinetic analysis

In this paper, the adsorption and degradation phenomenon involved in the photocatalytic degradation of dimethyl phthalates (DMPs) by titanium dioxide (TiO₂) was studied. A variety of operating variables were selected firstly. Then, it was proved that even for such weak adsorption properties molecules as DMP, adsorption was still an important prerequisite for photolysis. A surface-mediated reaction process was proposed that the photodegradation of DMP assisted by TiO₂ particles occurred primarily at the surface of the photocatalyst rather than in the homogeneous phase. According to Langmuir–Hinshelwood model, the adsorption constant determined from the dark adsorption was far less than that obtained in the light condition. Enhanced DMP adsorption on the surface of TiO₂ under irradiation was the possible reason for the improvement of photodegradation efficiency. Under the irradiation of light, a synergistic mechanism of adsorption and photocatalysis was responsible for DMP degradation. The quantitative analysis by adding scavengers indicated that ^·OH radical was primarily responsible for the photodegradation of DMP. It was further verified that ^·OH was produced much more from conduction band electrons rather than valance band holes toward photodegradation of DMP by adding foreign Cu²⁺.     

Up-to-date urban rainstorm intensity formulas considering spatial diversity in China

With accelerating urbanization in China, urban waterlogging has had a serious impact on urban sustainable development and citizen welfare. Simple urban rainstorm intensity formulas with a monotonous frequency distribution type cannot meet the practical needs of urban drainage planning and design. This study focuses on the development of urban rainstorm intensity formulas based on spatial diversity in China. Using the annual maximum sampling method, rainstorm data of 607 cities throughout China were collected into a database, with a total of 24,933 rainfall samples (annual observations) under various specified precipitation durations. The database was used to verify that integrating the Pearson III and Gumbel distributions would constitute an optimal theoretical distribution type, owing to its small error and increased fitting precision. Modification and coordination of four important parameters in the rainstorm intensity formula were done using a digital elevation model, which improved the accuracy of the formula. In addition, precipitation distributions in China were treated from the perspective of topographic features to validate the calculations from up-to-date formulas. Accuracy assessment was accomplished using a national code (GB-50014-2006), GIS-based isograms, and authoritative results from the Hydrological Bureau of the Ministry of Water Resources. This work provides a comprehensive foundation for the establishment of an up-to-date rainstorm intensity formula for China, which can be used widely in different cities.     

Satellite altimetry for measuring river stages in remote regions

The advent in satellite altimetry with the most accurate satellite radar altimeter since 1992 and its successive missions have enabled the routine global monitoring of water-level (or stage) for surface waters and changes in the quantities of dammed water reservoirs. However, satellite altimeter measurements typically have spatial resolution capable of observing only large water bodies, such as major lakes and rivers. This paper addresses the challenges of how to investigate water levels in medium (~?1 km in width) to small (~?100 m and narrower) rivers. Comparisons between the ENVISAT altimetry ICE-1 waveform retracking height and standard water-level measurements for multiple sections of Ohio River, Columbia River, and Red River of the North in the United States (US) reveals that the satellite altimetry measured water levels agree well with those observed at nearby US Geological Survey gaging stations over the 10-year period starting from 2002. The significant results include those obtained at Thompson, North Dakota (ND, correlation coefficient or R value of 0.76 between satellite and in situ water-level measurements) and Fargo, ND (R?=?0.74), where the stream channels of Red River are merely?~?50 m and ~?40 m wide, respectively, under normal climatic conditions. In addition, demonstrations of the approach over largely inaccessible portions of Tigris–Euphrates Rivers and Helmand River in the Middle East aided in understanding hydrology in these systems. This study demonstrates the ability of satellite radar altimetry to characterize rivers in these study regions which are much narrower than 100 m in width.     

Preparation of integrative cubes as a novel biological permeable reactive barrier medium for the enhancement of in situ aerobic bioremediation of nitrobenzene-contaminated groundwater

Recent studies have indicated that in situ aerobic bioremediation is one of the most effective methods for remediating groundwater contaminated with nitrobenzene (NB). Aerobic bioremediation is largely dependent on the maintenance of adequate dissolved oxygen (DO) levels in the groundwater. Traditional in situ aeration has various disadvantages, such as high operational costs and nonuniform aeration. In this study, integrative cubes were prepared and utilized as a novel biological permeable reactive barrier (bio-PRB) medium to enhance the aerobic bioremediation of NB-contaminated groundwater. The results revealed an NB removal rate greater than 98.68% after 15–20 days of continuous oxygen release from the bio-PRB medium. DO concentrations reached 8.0 mg/L during treatment, and NB-degrading bacteria were able to tolerate a range of pH conditions. This multifaceted bio-PRB medium can simultaneously adsorb and biodegrade NB, release oxygen, and neutralize the pH with phosphate buffer. The results of this study suggest that this bio-PRB medium represents a highly effective in situ bioremediation method for NB-contaminated groundwater.     

An insight into the surface velocity of Inylchek Glacier and its effect on Lake Merzbacher during 2006–2016 with Landsat time-series imagery

Mountain glacier is one of the extremely sensitive indicators for climate change, and its surface motion distribution and corresponding variation are valuable information for understanding ice mass exchange and glacier dynamics. This paper presents the long-term ice velocity distributions of Inylchek Glacier in the Tianshan region by pixel-tracking algorithm with time-series Landsat imagery acquired during 2006–2016. Then the monitored ice motion fields of Inylchek Glacier were carefully analyzed and revealed a generally similar spatial distribution characteristic. Most of the ice of the North Inylchek Glacier remains in a stagnant state except for the upstream part, but a relatively high velocity of 20–40 cm/day with an RMSE of 3 cm/day was observed on most part of the South Inylchek Glacier, except for the slow-moving glacier terminus. We also state the glacier dynamics around Lake Merzbacher and their possible effect on its glacier lake outburst flood (GLOF) risk. Besides, the surface velocity distribution on South Inylchek Glacier surface during the ablation period from 2014 to 2016 was also established and also compared with annual velocity. The corresponding difference yields that there is a positive relation between ice motion and temperature variation. Therefore, the time-series ice surface motion yielded by the Landsat imagery thus could provide us an efficient and low-cost way to analyze the current state and changes in glaciers, thanks to the continuous and regular spaceborne observations provided by the Landsat satellites.     

Development and verification of the comprehensive model for physical properties of hydrate sediment

Natural gas hydrate is widely distributed all over the world and may be a potential resource in the near future, whereas hydrate dissociation during the development affects wellbore stability and drilling safety. However, the present modeling of hydrate reservoir parameters ignored the influence of effective stress and only considered the hydrate saturation. In this paper, a series of stress sensitivity experiments for the unconsolidated sandstone were carried out, and the influence of mean effective stress on physical parameters was obtained; a comprehensive model for the physical parameters of hydrate reservoir was developed subsequently. With the help of ABAQUS finite element software, the established comprehensive model was verified by the use of the wellbore stability numerical model of hydrate reservoir. The verification results show that ignoring the effect of mean effective stress on the parameters of hydrate formation aggravates the invasion of drilling fluid into the hydrate formation. Besides, ignoring the stress sensitivity of reservoir physical parameters will underestimate the wellbore instability during hydrate drilling, which will be a threat to the safety of gas hydrate drilling. At the end of the drilling operation, the maximum plastic strain of the model for considering and not considering stress sensitivity was 0.0145 and 0.0138, respectively. Therefore, the established comprehensive model will provide a theoretical support for accurately predicting the engineering geological disasters in hydrate development process.