Comparison between cell-centered and nodal-based discretization schemes for linear elasticity

In this paper, we study newly developed methods for linear elasticity on polyhedral meshes. Our emphasis is on applications of the methods to geological models. Models of subsurface, and in particular sedimentary rocks, naturally lead to general polyhedral meshes. Numerical methods which can directly handle such representation are highly desirable. Many of the numerical challenges in simulation of subsurface applications come from the lack of robustness and accuracy of numerical methods in the case of highly distorted grids. In this paper, we investigate and compare the Multi-Point Stress Approximation (MPSA) and the Virtual Element Method (VEM) with regard to grid features that are frequently seen in geological models and likely to lead to a lack of accuracy of the methods. In particular, we look at how the methods perform near the incompressible limit. This work shows that both methods are promising for flexible modeling of subsurface mechanics.     

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.     

Lithospheric mantle beneath NE part of Bohemian Massif and its relation to overlying crust: new insights from Pilchowice xenolith suite,Sudetes, SW Poland

International Journal of Earth Sciences - The Oligocene/Miocene basanite from Pilchowice (Sudetes Mts., SW Poland) carries numerous small xenoliths of mantle peridotite, mostly harzburgite. The...     

Contrasting sources of Late Paleozoic rhyolite magma in the Polish Lowlands: evidence from U–Pb ages and Hf and O isotope composition in zircon

International Journal of Earth Sciences - The Polish Lowlands, located southwest of the Teisseyre–Tornquist Zone, within Trans-European Suture Zone, were affected by bimodal, but dominantly...     

3rd Regional Symposium on Landslides in the Adriatic-Balkan Region (3rd ReSyLAB)—a final report

The 3rd Regional Symposium on Landslides in the Adriatic-Balkan Region (3rd ReSyLAB) was held in Ljubljana, Slovenia, from June 11 to 13, 2017, with 70 participants from nine countries (Austria, Bosnia and Hercegovina, Croatia, Czech Republic, Italy, Republic of Macedonia, Serbia, Slovenia, Spain)—scientists, engineers, researchers, students, experts, politicians, and other decision-makers working in the area of landslide risk reduction in the region. The ReSyLAB is a biannual event organized by the Adriatic-Balkan Network of the International Consortium on Landslides (ICL ABN). Being an important form of activities of this ICL regional network comprising of six ICL members from four countries, it was also a contribution of the International Consortium on Landslides (ICL) to the implementation of the Sendai Framework for Disaster Risk Reduction 2015–2030. This article reports on the main outcomes of the 3rd ReSyLAB Symposium. Altogether, 41 abstracts were published in the symposium book of abstracts, and the symposium proceedings with over 20 reviewed full papers are under preparation to be printed early in 2018. During the 3rd ReSyLAB, a five invited keynote lectures have been presented, and 28 oral presentations are given to the audience. An important part of the symposium was a Round Table entitled “Enhancing cooperation between landslide research community and end users.” On the last day of the symposium, over 30 experts participated in two post-symposium study tours in Slovenia.     

Evidence for the occurrence of hydrogeochemical processes in the groundwater of Khoy plain,northwestern Iran,using ionic ratios and geochemical modeling

Rapid population growth, industrialization, and agricultural expansion in the Khoy area (northwestern Iran) have led to its dependence on groundwater and degradation of groundwater quality. This study attempts to decipher the major processes and factors that degrade the groundwater quality of the Khoy plain. For this purpose, 54 groundwater samples from unconfined and confined aquifers of the plain were collected in July 2017 and analyzed for major cations and anions (Na, K, Ca, Mg, HCO₃, SO_4, and Cl), minor ions (NO₃ and F), and Al. Magnesium and bicarbonate were identified as the dominant cation and anion, respectively. Several ionic ratios and geochemical modeling using PHREEQC indicated that the most important hydrogeochemical processes to affect groundwater quality in the plain were weathering and dissolution of evaporitic and silicate minerals, mixing, and ion exchange. There were smaller effects from evaporation and anthropogenic factors (e.g., industries). Results showed that the high salinity of the groundwater in the northeast area of the plain was due to the high solubility of the evaporitic minerals, e.g., halite and gypsum. Reverse ion exchange and the contribution of mineral dissolution were more significant than ion exchange in the northeastern part of the plain. Elevated salinity of the groundwater in the southeast was attributed mostly to reverse ion exchange and somewhat to evaporation.