Computing derivative information of sequentially coupled subsurface models

A generic framework for the computation of derivative information required for gradient-based optimization using sequentially coupled subsurface simulation models is presented. The proposed approach allows for the computation of any derivative information with no modification of the mathematical framework. It only requires the forward model Jacobians and the objective function to be appropriately defined. The flexibility of the framework is demonstrated by its application in different reservoir management studies. The performance of the gradient computation strategy is demonstrated in a synthetic water-flooding model, where the forward model is constructed based on a sequentially coupled flow-transport system. The methodology is illustrated for a synthetic model, with different types of applications of data assimilation and life-cycle optimization. Results are compared with the classical fully coupled (FIM) forward simulation. Based on the presented numerical examples, it is demonstrated how, without any modifications of the basic framework, the solution of gradient-based optimization models can be obtained for any given set of coupled equations. The sequential derivative computation methods deliver similar results compared to FIM methods, while being computationally more efficient.     

Precursory slope distress prior to the 2010 Mount Meager landslide,British Columbia

In 2010, the south flank of Mount Meager failed catastrophically, generating the largest (53 ± 3.8 × 10⁶ m³) landslide in Canadian history. We document the slow deformation of the edifice prior to failure using archival historic aerial photographs spanning the period 1948–2006. All photos were processed using Structure from Motion (SfM) photogrammetry. We used the SfM products to produce pre-and post-failure geomorphic maps that document changes in the volcanic edifice and Capricorn Glacier at its base. The photographic dataset shows that the Capricorn Glacier re-advanced from a retracted position in the 1980s then rapidly retreated in the lead-up to the 2010 failure. The dataset also documents 60 years of progressive development of faults, toe bulging, and precursory failures in 1998 and 2009. The 2010 collapse was conditioned by glacial retreat and triggered by hot summer weather that caused ice and snow to melt. Meltwater increased pore water pressures in colluvium and fractured rocks at the base of the slope, causing those materials to mobilize, which in turn triggered several secondary failures structurally controlled by lithology and faults. The landslide retrogressed from the base of the slope to near the peak of Mount Meager involving basement rock and the overlying volcanic sequence. Elsewhere on the flanks of Mount Meager, large fractures have developed in recently deglaciated areas, conditioning these slopes for future collapse. Potential failures in these areas have larger volumes than the 2010 landslide. Anticipated atmospheric warming over the next several decades will cause further loss of snow and glacier ice, likely producing additional slope instability. Satellite- and ground-based monitoring of these slopes can provide advanced warning of future landslides to help reduce risk in populated regions downstream.     

Ecotoxicological evaluation of dredged sediments from Rodrigo de Freitas Lagoon (Rio de Janeiro State,Brazil) using bioassays with earthworms and collembolans

The toxicity and metal bioavailability were studied in dredged sediments from Rodrigo de Freitas Lagoon (Rio de Janeiro, RJ, Brazil) using acute and avoidance tests with Eisenia andrei, and reproduction tests with Folsomia candida. The sediment was mixed with an artificial soil, and two natural soils (ferralsol and chernosol—representative Brazilian tropical soils) to obtain the following doses: 1, 3, 6, 12, 18, 24 and 30%. Total metal concentrations were determined in the sediment to support the interpretation of ecotoxicological data. Metal concentrations in the mixtures were in agreement with the threshold limits established by Brazilian law. However, significant avoidance responses were found on doses ≥?3% and were the most sensitive endpoint. Earthworm mortality found in artificial soil mixtures (LC50?=?3.9) suggests higher toxicity levels than those obtained in ferralsol (LC50?=?7.6%) and chernosol (11.0%) treatments. Earthworm mortality, avoidance responses and collembolan reproduction levels found in ferralsol mixtures (LC50?=?9.2; avoidance EC50?=?2.3%; reproduction EC50?=?2.8%) were higher compared to chernosol treatments (LC50?=?11.0%; avoidance EC50?=?4.3%; reproduction EC50?=?4.9%). The reduction of toxicity levels in chernosol mixtures is probably due to the abundance of expansive clay minerals in chernosols with capacity of adsorbing metals and other xenobiotic substances from soil pore water, decreasing metal bioavailability. Finally, threshold limits defined by Brazilian legislation for soil quality and land disposal of dredged sediments are not sufficient to prevent noxious effects on soil fauna and should be complemented with a preliminary ecotoxicological evaluation.     

Construction and Characterization of Transportable Calibration Pads for Calibration of Portable and Airborne Gamma Ray Spectrometers in India

Characterization of the Panandhro lignite deposits from western Indian state of Gujarat, based on the geochemical and palynological evidences, has been performed to assess the floral composition, maturity and hydrocarbon potential of the sequence. Elementally, the lignites consist of moderate carbon, low hydrogen and moderate sulfur contents. The samples are characterized by high TOC contents (lignite: av. 46.43 wt.%, resin: 62.47 wt.%). The average HI values for the lignite is 136 mg HC/g TOC, and that of the associated resin is 671 mg HC/g TOC. The highest T_max is recoded in lignite (422°C) and lowest in the resin (39°C) samples. The FTIR spectrum of lignite is characterized by highly intense OH stretching peak ~3350 cm^-1, aliphatic CH_x stretching peaks between 3000-2800 cm^-1, aromatic C=O stretching and an aromatic C=C stretching. The spectrum of resin shows strongest absorption due to aliphatic CHx stretching between 2940-2915 cm^-1 and 2870-2850 cm^-1, and deformation by the medium peak between 1450 and 1650 cm^-1. The recovered palynofloral assemblage indicates the dominance of angiosperm pollen grains with maximum abundance of Arecaceae family, and subdominant pteridophytic spores. Marine influence is indicated by the presence of abundant dinoflagellate cysts. The occurrence of flora from a variety of ecological niches suggests a luxuriant diverse vegetation pattern existed in the vicinity of depositional site under humid tropical conditions. The overall characteristics of the lignite deposits point towards their ability to generate (upon maturation) hydrocarbons as they have types III–II admixed kerogen (organic matters).     

Zoning and exsolution in alkali feldspars from Laacher See volcano (Western Germany): constraints on temperature history prior to eruption

Compositional zoning and exsolution patterns of alkali feldspars in carbonatite-bearing cognate syenites from the 6.3 km³ (D.R.E) phonolitic Laacher See Tephra (LST) deposit in western Germany (12.9 ka) are reported. These rocks represent the cooler outer portion and crystal-rich products of a cooling magma reservoir at upper crustal levels. Major and trace-element difference between cores and rims in sanidine crystals represent two generations of crystal growth separated by unmixing of a carbonate melt. Trace-element differences measured by LA–ICP–MS are in accordance with silicate–carbonate unmixing. Across the core–rim boundary, we extracted gray-scale profiles from multiple accumulations of back-scattered electron images. Gray scales directly represent K/Na ratios owing to low concentrations of Ba and Sr (<?30 ppm). Diffusion gradients are modeled to solve for temperature using known pre-eruptive U–Th zircon ages (0–20 ky) of each sample (Schmitt et al., J Petrol 51:1053–1085, 2010). Estimated temperatures range from 630 °C to 670 °C. For the exsolution boundaries, a diffusive homogenization model is constrained by the solvus temperature of ~ 712^_725 °C and gives short time scales of only 15–50 days. Based on our results, we present a model for the temperature–time history of these rocks. The model also constrains the thermal variation across the cooling crystal-rich carapace of the magma reservoir over 20 ka and suggests a thermal reactivation of cumulates, the cooling carapace, and probably the entire system only a few years prior to the explosive eruption of the remaining molten core of the phonolitic magma reservoir.