Effects of grain morphology on critical state: a computational analysis

We introduce a new DEM scheme (LS-DEM) that takes advantage of level sets to enable the inclusion of real grain shapes into a classical discrete element method. Then, LS-DEM is validated and calibrated with respect to real experimental results. Finally, we exploit part of LS-DEM potentiality by using it to study the dependency of critical state (CS) parameters such as critical state line (CSL) slope \(\lambda \), CSL intercept \(\varGamma \), and CS friction angle \(\varPhi _{\mathrm{CS}}\) on the grain’s morphology, i.e., sphericity, roundness, and regularity. This study is carried out in three steps. First, LS-DEM is used to capture and simulate the shape of five different two-dimensional cross sections of real grains, which have been previously classified according to the aforementioned morphological features. Second, the same LS-DEM simulations are carried out for idealized/simplified grains, which are morphologically equivalent to their real counterparts. Third, the results of real and idealized grains are compared, so the effect of “imperfections” on real particles is isolated. Finally, trends for the CS parameters (CSP) dependency on sphericity, roundness, and regularity are obtained as well as analyzed. The main observations and remarks connecting particle’s morphology, particle’s idealization, and CSP are summarized in a table that is attempted to help in keeping a general picture of the analysis, results, and corresponding implications.     

Strength criterion for cross-anisotropic sand under general stress conditions

By incorporating the fabric effect and Lode’s angle dependence into the Mohr–Coulomb failure criterion, a strength criterion for cross-anisotropic sand under general stress conditions was proposed. The obtained criterion has only three material parameters which can be specified by conventional triaxial tests. The formula to calculate the friction angle under any loading direction and intermediate principal stress ratio condition was deduced, and the influence of the degree of the cross-anisotropy was quantified. The friction angles of sand in triaxial, true triaxial, and hollow cylinder torsional shear tests were obtained, and a parametric analysis was used to detect the varying characteristics. The friction angle becomes smaller when the major principal stress changes from perpendicular to parallel to the bedding plane. The loading direction and intermediate principal stress ratio are unrelated in true triaxial tests, and their influences on the friction angle can be well captured by the proposed criterion. In hollow cylinder torsional shear tests with the same internal and external pressures, the loading direction and intermediate principal stress ratio are related. This property results in a lower friction angle in the hollow cylinder torsional shear test than that in the true triaxial test under the same intermediate principal stress ratio condition. By comparing the calculated friction angle with the experimental results under various loading conditions (e.g., triaxial, true triaxial, and hollow cylinder torsional shear test), the proposed criterion was verified to be able to characterize the shear strength of cross-anisotropic sand under general stress conditions.     

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.     

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.     

Crustal melting and recycling: geochronology and sources of Variscan syn-kinematic anatectic granitoids of the Tormes Dome (Central Iberian Zone). A U–Pb LA-ICP-MS study

International Journal of Earth Sciences - In this study, we report U–Pb Laser Ablation ICP-MS zircon and ID-TIMS monazite ages for peraluminous granitoid plutons...     

A Cambrian mixed carbonate–siliciclastic platform in SW Gondwana: evidence from the Western Sierras Pampeanas (Argentina) and implications for the early Paleozoic paleogeography of the proto-Andean margin

International Journal of Earth Sciences - The Western Sierras Pampeanas (WSP) of Argentina record a protracted geological history from the Mesoproterozoic assembly of the Rodinia supercontinent to...     

Pico- and Nanophytoplankton Dynamics in Two Coupled but Contrasting Coastal Bays in the NW Mediterranean Sea (France)

Due to its ecological context, the Toulon bay represents a site of scientific interest to study temporal plankton distribution, particularly pico- and nanophytoplankton dynamics. A monthly monitoring was performed during a two-year cycle (October 2013–December 2015) at two coupled sampling sites, referred to as Little and Large bays, which had different morphometric characteristics and human pressures. Flow cytometry analyses highlighted the fact that pico- and nanophytoplankton were more abundant in the eutrophic Little bay. Furthermore, it evidenced two community structures across the Toulon bays: at times, a co-dominance of picoeukaryotes, nanoeukaryotes, Synechococcus 1-like cells and Prochlorococcus-like cells was found, and at other times, a Synechococcus 1-like dominated community existed. The alternation of one structure or the other can be explained by a combined action of temperature regime, nutrient conditions and degree of contamination. This study showed that pico- and nanophytoplankton dynamics were mainly driven by temperature in both sites, as in other temperate Mediterranean regions. Thus, the community was mainly composed of picoeukaryotes and Prochlorococcus-like cells in the winter (<?15 °C), while it was dominated by Synechococcus 1-like cells in the summer (>?20 °C). Additionally, the multiple human stressors in the Little bay seemed to affect the increase in abundance of Synechococcus 1-like cells as they were preferentially observed in the Large bay.     

Rainfall thresholds for the triggering of landslides considering previous soil moisture conditions (Asturias,NW Spain)

Given its geological and climatic conditions and its rugged orography, Asturias is one of the most landslide prone areas in the North of Spain. Most of the landslides occur during intense rainfall episodes. Thus, precipitation is considered the main triggering factor in the study area, reaching average annual values of 960 mm. Two main precipitation patterns are frequent: (i) long-lasting periods of moderate rainfall during autumn and winter and (ii) heavy short rainfall episodes during spring and early summer. In the present work, soil moisture conditions in the locations of 84 landslides are analysed during two rainfall episodes, which represent the most common precipitation patterns: October–November 2008 and June 2010. Empirical data allowed the definition of available water capacity percentages of 99–100% as critical soil moisture conditions for the landslide triggering. Intensity-duration rainfall thresholds were calculated for each episode, considering the periods with sustained high soil moisture levels before the occurrence of each analysed landslide event. For this purpose, data from daily water balance models and weather stations were used. An inverse relationship between the duration of the precipitation and its intensity, consistent with published intensity-duration thresholds, was observed, showing relevant seasonal differences.