Robust numerical implementation of a 3D rate-dependent model for reservoir geomechanical simulations

A 3D elasto-plastic rate-dependent model for rock mechanics is formulated and implemented into a Finite Element (FE) numerical code. The model is based on the approach proposed by Vermeer and Neher (A soft soil model that accounts for creep. In: Proceedings of the International Symposium “Beyond 2000 in Computational Geotechnics,” pages 249-261, 1999). An original strain-driven algorithm with an Inexact Newton iterative scheme is used to compute the state variables for a given strain increment.The model is validated against laboratory measurements, checked on a simplified test case, and used to simulate land subsidence due to groundwater and hydrocarbon production. The numerical results prove computationally effective and robust, thus allowing for the use of the model on real complex geological settings.     

A seismic source zone model for the seismic hazard assessment of the Italian territory

We designed a new seismic source model for Italy to be used as an input for country-wide probabilistic seismic hazard assessment (PSHA) in the frame of the compilation of a new national reference map.

We started off by reviewing existing models available for Italy and for other European countries, then discussed the main open issues in the current practice of seismogenic zoning.

The new model, termed ZS9, is largely based on data collected in the past 10 years, including historical earthquakes and instrumental seismicity, active faults and their seismogenic potential, and seismotectonic evidence from recent earthquakes. This information allowed us to propose new interpretations for poorly understood areas where the new data are in conflict with assumptions made in designing the previous and widely used model ZS4.

ZS9 is made out of 36 zones where earthquakes with M_w > = 5 are expected. It also assumes that earthquakes with M_w up to 5 may occur anywhere outside the seismogenic zones, although the associated probability is rather low. Special care was taken to ensure that each zone sampled a large enough number of earthquakes so that we could compute reliable earthquake production rates.

Although it was drawn following criteria that are standard practice in PSHA, ZS9 is also innovative in that every zone is characterised also by its mean seismogenic depth (the depth of the crustal volume that will presumably release future earthquakes) and predominant focal mechanism (their most likely rupture mechanism). These properties were determined using instrumental data, and only in a limited number of cases we resorted to geologic constraints and expert judgment to cope with lack of data or conflicting indications. These attributes allow ZS9 to be used with more accurate regionalized depth-dependent attenuation relations, and are ultimately expected to increase significantly the reliability of seismic hazard estimates.     

A Refinement of the McMillen (1988) Recursive Digital Filter for the Analysis of Atmospheric Turbulence

We present a refinement of the recursive digital filter proposed by McMillen (Boundary-Layer Meteorol 43:231–245, 1988), for separating surface-layer turbulence from low-frequency fluctuations affecting the mean flow, especially over complex terrain. In fact, a straightforward application of the filter causes both an amplitude attenuation and a forward phase shift in the filtered signal. As a consequence turbulence fluctuations, evaluated as the difference between the original series and the filtered one, as well as higher-order moments calculated from them, may be affected by serious inaccuracies. The new algorithm (i) produces a rigorous zero-phase filter, (ii) restores the amplitude of the low-frequency signal, and (iii) corrects all filter-induced signal distortions.     

旋转扰动流体中壁面湍流不稳定及其拟序结构研究Ⅰ:实验

重点介绍和讨论了中性条件下旋转扰动流体中边界层强迫不稳定及其相关的一些问题,阐述了旋转体系中切变驱动边界层不稳定的动力学特征.这些不稳定状态的研究在大气物理学、流体动力学、海洋学等多个领域中引起科学家极大的兴趣,近年来在实验和理论研究中都得到了不断的发展.意大利都灵大学基础物理系地球科学实验组通过水槽旋转实验方法,不断改变水槽启动或结束时的旋转运动速度,以及底部壁面粗糙度等要素,所得到的实验结果与SDBL理论非常一致.     

A high‐altitude temporary spring in a compartmentalized carbonate aquifer: the role of low‐permeability faults and karst conduits

The aim of this research was to refine the actual conceptual model related to the activation of high‐altitude temporary springs within the carbonate Apennines in southern Italy. The research was carried out through geophysical, hydrogeological, hydrochemical and isotopic investigations at the Acqua dei Faggi experimental site during five hydrologic years. The research demonstrated that, in carbonate aquifers where low‐permeability faults cause the aquifer system to be compartmentalized, high‐altitude temporary springs may be recharged by groundwater. In such settings, neither surface water infiltration in karst systems nor perched temporary aquifers play a role of utmost importance. The rare (once or a few time a year) activation of such springs is due to the fact that groundwater unusually reach the threshold head that allows the spring to flow. The activation of the studied high‐altitude temporary spring also depended on relationships between a low‐permeability fault core and a karst system that locally interrupts the low‐permeability barrier. In fact, when the hydraulic head did not reach the karst system, the concentrated head loss within the fault core did not allow the spring to flow, because the groundwater entirely flowed through the fault towards the downgradient compartment. Copyright © 2009 John Wiley & Sons, Ltd.     

Measure of the exponential splitting of the homoclinic tangle in four-dimensional symplectic mappings

Using four-dimensional symplectic maps as a model problem, we numerically compute the unstable manifolds of the hyperbolic manifolds of the phase space related to the single resonances. We measure an exponential dependence of the size of the lobes of these manifolds through many orders of magnitude of the perturbing parameter. This is an indirect numerical verification of the exponential decay of the normal form, as predicted by the Nekhoroshev theorem. The variation of the size of the lobes turns out to be correlated to the diffusion coefficient.     

A numerical study of the hyperbolic manifolds in a priori unstable systems. A comparison with Melnikov approximations

Using numerical methods we study the hyperbolic manifolds in a model of a priori unstable dynamical system. We compare the numerically computed manifolds with their analytic expression obtained with the Melnikov approximation. We find that, at small values of the perturbing parameter, the topology of the numerically computed stable and unstable manifolds is the same as in their Melnikov approximation. Increasing the value of the perturbing parameter, we find that the stable and unstable manifolds have a peculiar topological transition. We find that this transition occurs near those values of the perturbing parameter for which the error terms of Melnikov approximations have a sharp increment. The transition value is also correlated with a change in the behaviour of dynamical quantities, such as the largest Lyapunov exponent and the diffusion coefficient.     

Volcanic geomorphology and tectonics of the Aeolian archipelago (Southern Italy) based on integrated DEM data

The geomorphological and morphometric analysis of the sea floor topography surrounding the Aeolian Islands, South Tyrrhenian Sea, Italy, provides insights into the relationships between the volcanological evolution of the islands and their tectonic features. We constructed geomorphological maps of the submarine portions of the seven large edifices constituting the islands on the basis of a DEM with a 5 m resolution step. These maps include constructional and destructional landforms such as submarine volcanic vents located west of Lipari and north of Alicudi, and hummocky surfaces recognised north of Lipari and Salina. The latter landforms, together with the occurrence of large scars affecting the main edifices on land, suggest that sector collapses affected some islands. Geomorphological data indicate that the location of subaerial and submarine vents is strongly controlled by local tectonic structures striking WNW-ESE (Alicudi-Filicudi sector), NNW-SSE (Salina-Lipari-Vulcano sector) and NE-SW (Panarea-Stromboli sector). The islands can be divided into two groups on the basis of some morphometric parameters: a first group with a pancake-like shape, Dp/D (abrasion platform diameter/basal diameter) higher than 0.40 and H/D (total height/basal diameter) lower than 0.13, and a second group with a conical shape, characterised by Dp/D lower than 0.34 and H/D higher than 0.14. These ratios and other morphometric parameters reflect the different volcanological and structural evolution of the Aeolian Islands. The pancake-like shaped complexes have been created, in addition to their submarine stage, by extrusive and highly explosive activity, whereas the cone-shaped edifices have been characterised by effusive or moderate explosive activity.Editorial responsibility: C Kilburn