In geo-engineering, hydraulic fracturing (HF) is used to improve the oil/gas production and optimize well stimulation in low permeability reservoirs. To design HF treatments, it is necessary to predict the pattern of fracture geometry as a function of treatment parameters. Present day HF simulators used in the industry are based on rather old approaches, relying on questionable assumptions about fracture geometry and disregarding hydro-mechanical coupling when dealing with fluid pressure evolution. As a consequence, most of the commercial codes fail in the reproduction of the complex intricate fracture patterns shown by the field acoustic measurements [1].
In this contribution we present the application of a recently developed model of brittle damage of confined rock masses to simulate HF treatments. The model is based on the explicit micromechanical construction of connected patterns of parallel equi-spaced cracks, by extending to saturated porous media the dry multi-scale brittle damage model firstly introduced in [2]. A relevant feature of the model is that the fracture patterns are not arbitrary, but their inception, orientation and spacing follow from energetic consideration. The capability of the model to predict the mechanical response of brittle rocks and the related porosity and permeability evolution has been presented in [3] and [4].
The model, based on the Terzaghi effective stress concepts, has been implemented into a coupled hydro-mechanical finite element code, where the linear momentum and the fluid mass balance equations are numerically solved via a staggered approach. The coupled code is used to simulate laboratory and field scale fracturing processes induced by an increase in pore pressure, as in hydraulic fracturing jobs for reservoir stimulation. A reference case, based on typical data of a HF process, has been used to ascertain the influence of different operational parameters on the outcomes of the process. In particular, with reference to a simple geometry, a parametric study is performed by varying the number of wellbore perforations, the number of fracking processes, the distance between the different perforation slots, the wellbore deviation from the minimum principal stress direction axis, fluid pressure, fluid density, and proppant size.
The examples show the capability of the model in reproducing three-dimensional multiscale complex fracture patterns and permeability enhancement in the damaged porous medium.
Non-Extensive Statistical Physics (NESP) is based on a generalization of the standard Boltzmann-Gibbs expression of the entropy. It was introduced a few decades ago by Tsallis [1] in an attempt to describe and enlighten phenomena with “anomalous” behaviour for which the statistical mechanical concepts of the Boltzmann-Gibbs approach have been proven inadequate. A typical class of such phenomena are the ones involving long-range interactions and memory effects. Given that fracture is characterized by both long-range interactions and memory effects [2], it is reasonable to examine whether the mechanical response of brittle building materials at load levels approaching these causing fractures could be described in terms of concepts based on NESP. The challenge is to detect proper and clear indicators that could be considered as early warning signals of the upcoming entrance of the system (loaded specimen or structure) into its “critical stage”, i.e., that of impending fracture. In this direction, advantage is taken of experimental data gathered from a long series of experimental protocols with specimens made of brittle building materials submitted to a variety of loading schemes (direct tension, uniaxial compression, three-point bending, shear) either monotonic or stepwise. The protocols include both elementary and structural tests. The study is carried out in terms of characteristic parameters of the Acoustic Emissions detected during loading. It is concluded that the time evolution of the entropic parameters of NESP provide reliable pre-failure indicators, in good agreement to the respective ones obtained from completely different analyses of the acoustic activity (as it is, for example, the F-function [3] or the Ib-value) or even from different monitoring tools (as it is, for example, the Pressure Stimulated Currents technique [4]).
Keywords:This paper is focused on numerical analysis of the coupled hydro-mechanical response of geomaterials that contain pre-existing or newly developing zones of localized deformation, such as faults, macrocracks, shear bands, etc. In terms of mechanical response, the existence of discontinuities results in the reduction of material strength that is triggered by sliding/separation along the defects. In addition, the pre-existing cracks act as stress concentrators prompting the formation of new macrocracks that may propagate through the domain. In terms of flow properties, the hydraulic conductivity is also strongly affected by the fracture pattern and displays anisotropy at the macroscale.
The numerical analysis of flow through fractured porous media is usually conducted by employing the Extended Finite Element Method (XFEM). The approach typically involves the assumption that the fluid pressure is continuous across the discontinuity, while the pressure gradient is discontinuous. Such formulation allows for the transport and storage of fluid inside the crack. The jump in the pressure gradient is achieved by partitioning the pressures at both sides of the discontinuity by a signed distance function. Although the approach is accurate, it is computationally very inefficient. The latter stems from incorporation of enriched DOSs, i.e. additional degrees of freedom that account for the presence of discontinuities, as well the need for partitioning of the domain with triangular sub-elements for the Gaussian integration scheme.
Given the problems inherent to XFEM approach, a new formulation is developed here that employs the averaging of the field operators within the referential volume adjacent to macrocrack. This leads to an enriched form of Darcy’s law, which incorporates the notion of equivalent conductivity. The latter is defined as a symmetric second-order tensor whose components are function of hydraulic properties of constituents (viz. intact material and fractured region) as well as the internal length parameter. Such an approach does not require any additional degrees of freedom to account for the presence of discontinuities, which significantly improves the computational efficiency as compared to XFEM.
The mechanical analysis incorporates an enhanced embedded discontinuity approach, which is conceptually similar to that employed for specification of equivalent hydraulic conductivity. It employs the same ‘characteristic dimension’ related to geometry of fractures and enables a discrete tracing of the propagation of new cracks.
The formulation is illustrated by a number of examples. In particular, a series of compression tests on pre-fractured rock-like samples is simulated. Various geometric configurations of pre-existing as well as new propagating cracks are considered and the results are compared with the experimental data. In addition, a steady-state flow through the fractured domain under a prescribed hydraulic gradient is examined for different geometries of fractures. Finally, a coupled problem is considered involving a transient flow under constant traction boundary conditions.
The gradient approach is applied to discuss disintegration and size effects in cavities and deep level tunnels in geomechanics. Previous results of classical theories are revisited and further extended by adjusting the newly introduced gradient coefficient which multiplies the gradient term and accounts for nonlocality and underlying microstructural heterogeneity. Qualitative comparison with related observations is discussed.
Keywords:It is known that the classical elastic beam model fails in the indentation analysis of an elastic beam by a rigid indenter. A simple higher-order beam model developed based on the Kerr-type differential relation between the indentation pressure and the deflection of pressured surface of the beam is applied to study the indentation of an elastic beam by a rigid circular cylinder. The proposed method is validated by comparing its predicted results with known data, and the merit of the proposed model is demonstrated by some new easy-to-use explicit formulas and numerical results. In particular, the present model confirms that the contact zone becomes two separate strips when its width increases and exceeds a certain critical value. It is expected that this simple higher-order beam model could be found useful in the mechanical indentation analysis of some elastic beam problems for which the classical beam model fails.
Keywords:Granular materials are ubiquitous in nature and industry. From landslides to extrusion of animal feed, particulate materials are relevant for all aspects of human life. Despite our frequent exposure to these materials, their handling still confronts us with unique engineering challenges. Notably the flow behavior of these materials is spatiotemporally heterogeneous and non-local in nature, both of which present both experimental and theoretical complexities. Granular materials however also present us with an opportunity to systematically investigate this behavior, as we are nowadays uniquely capable of tuning particle properties to a great extent, and can study their mechanical behavior in ever greater detail. We foresee that studying granular materials thus offers a window on the behavior of a broad range of classes of materials that display heterogeneous behavior. In my talk I will discuss the latest experimental developments on granular flow measurements and related modeling attempts.
We present recent advances in the modelling and simulation of fracture in heterogeneous materials, either natural, such as geomaterials, or manmade, such as aerospace and energy materials.
We discuss how imaged microstructures can be used to model the macroscopic behaviour of engineering structures used in various areas of engineering.
The main directions discussed include:
- acceleration of multi-scale methods
- machine-learning-enhanced model order reduction
- applications in chemistry, engineering, materials and medicine
- stochastic approaches for uncertainty quantification
We show how recent developments are used in practice and transferred to companies and start-ups and discuss innovation and knowledge transfer in the field of computational methods for engineering and sciences.
We present recent advances in the modelling and simulation of fracture in heterogeneous materials, either natural, such as geomaterials, or manmade, such as aerospace and energy materials. We discuss how imaged microstructures can be used to model the macroscopic behaviour of engineering structures used in various areas of engineering. The main directions discussed include: acceleration of multi-scale methods, machine-learning-enhanced model order reduction, applications in chemistry, engineering, materials and medicine, and stochastic approaches for uncertainty quantification. We show how recent developments are used in practice and transferred to companies and start-ups and discuss innovation and knowledge transfer in the field of computational methods for engineering and sciences.
A statistical study of precursor activity in rain- as well as earthquake-induced landslides by means of spring block models enhanced with displacement gradients and stochasticity is performed. This way, a robust 2D model can be formulated for studying triggered landslides. A cellular automaton is utilized, in order to examine the dynamic behavior and the stability of rock/soil slopes due to neighboring heavy rainfall or earthquake activity. The type and nature of the failure plane, as well as the triggering mechanism is studied. Moreover, the different dynamic evolution modes of the slope can be mapped to specific shape parameters of the corresponding distributions of the incremental displacements. These parameters are calibrated through comparison with statistical data on landslides events available in the literature. The calibrated model can then be used as a means for understanding, predicting and mitigating the impact of catastrophic landslides and its theoretical predictions are compared with the respective predictions of alternative techniques for studying slope stability.
Keywords:The broader area of the Aegean, i.e. the southern Balkans, mainland Greece, the Aegean sea with more than three thousand islands and islets and the Western Asia Minor (Anatolia), is particularly geodynamically active today. It has a long geological history of more than 250 million years. It raised from Pangaea and the deep sediments of the Tethys sea, to today's mountain ranges, with the seas and the many bays, the peninsulas and the impressive underwater relief. This place, where Western Civilization has its deep roots, in the Cycladic culture and classical antiquity, the concept of Europe was shaped by mythology and philosophy. Ancient Greek mythology created gods and deities as manifestations of natural phenomena associated with geological processes and as geomythology they are presented and interpreted today. This long-suffering, hard and beautiful place, crossroads of gods and people, myths and history, culture and wars, scientific geological laboratory, place of art and grandeur, prosperity and decline, “amaranth sea”, where “life harmonizes ceaselessly and the pulse of the earth beats calmly”, it was praised by poets, from Homer to the modern Nobel laureate Odysseus Elytis in his great poem "AXION ESTI” (Worthy Being), in which all its geological history is transcendentally described up to the struggles of its inhabitants.
