On the parallel solution of hydro-mechanical problems with fracture networks and contact conditions

The paper presents a numerical method for simulating flow and mechanics in fractured rock. The governing equations that couple the effects in the rock mass and in the fractures are obtained using the discrete fracture-matrix approach. The fracture flow is driven by the cubic law, and the contact conditions prevent fractures from self-penetration. A stable finite element discretization is proposed for the displacement-pressure-flux formulation. The resulting nonlinear algebraic system of equations and inequalities is decoupled using a robust iterative splitting into the linearized flow subproblem, and the quadratic programming problem for the mechanical part. The non-penetration conditions are solved by means of dualization and an optimal quadratic programming algorithm. The capability of the numerical scheme is demonstrated on a benchmark problem for tunnel excavation with hundreds of fractures in 3D. The paper's novelty consists in a combination of three crucial ingredients: (i) application of discrete fracture-matrix approach to poroelasticity, (ii) robust iterative splitting of resulting nonlinear algebraic system working for real-world 3D problems, and (iii) efficient solution of its mechanical quadratic programming part with a large number of fractures in mutual contact by means of own solvers implemented into an in-house software library

Dynamic development of hydrofracture

Many natural examples of complex joint and vein networks in layered sedimentary rocks are hydrofractures that form by a combination of pore fluid overpressure and tectonic stresses. In this paper, a two-dimensional hybrid hydro-mechanical formulation is proposed to model the dynamic development of natural hydrofractures. The numerical scheme combines a discrete element model (DEM) framework that represents a porous solid medium with a supplementary Darcy based pore-pressure diffusion as continuum description for the fluid. This combination yields a porosity controlled coupling between an evolving fracture network and the associated hydraulic field. The model is tested on some basic cases of hydro-driven fracturing commonly found in nature, e.g., fracturing due to local fluid overpressure in rocks subjected to hydrostatic and nonhydrostatic tectonic loadings. In our models we find that seepage forces created by hydraulic pressure gradients together with poroelastic feedback upon discrete fracturing play a significant role in subsurface rock deformation. These forces manipulate the growth and geometry of hydrofractures in addition to tectonic stresses and the mechanical properties of the porous rocks. Our results show characteristic failure patterns that reflect different tectonic and lithological conditions and are qualitatively consistent with existing analogue and numerical studies as well as field observations. The applied scheme is numerically efficient, can be applied at various scales and is computational cost effective with the least involvement of sophisticated mathematical computation of hydrodynamic flow between the solid grains

Fluid Flow Simulation in Fractured Reservoirs

The purpose of this study is to analyze fluid flow in fractured reservoirs. In most petroleum reservoirs, particularly carbonate reservoirs and some tight sands, natural fractures play a critical role in controlling fluid flow and hence production. Uncertainties involved in the understanding of fracture architecture and properties often propagate in the construction of reservoir flow models. The state-of-the-art reservoir simulation packages used widely in the oil industry often do not take into account the complex random geometry of real fracture systems that can vary from one grid-block to another, and sometimes even within a single grid-block. The reason for this is two-fold: there exists no technology as yet to image the micro-fractures in-situ, and most of the reservoir modeling software does not use micro-scale flow equations to model the change of flow variables. Fractures are highly conductive channels for flow among all types of porous-permeable media, Flow through them can vary widely depending on different fracture properties, such as fracture apertures, densities, asperities, etc. The objective of this study is to understand fluid flow in fractures using a finite difference approach, and to analyze the effects of fracture properties on flow mainly through visualization. Contrary to the conventional macro-scale modeling approach, micro-scale simulation is carried out. The ultimate goal is to incorporate this information into reservoir scale modeling schemes.United States. Dept. of Energy (Grant DE-FG0s-00ER15041)Massachusetts Institute of Technology. Earth Resources LaboratoryEni S.p.A. (Firm

Sequential approach to joint flow-seismic inversion for improved characterization of fractured media

Seismic interpretation of subsurface structures is traditionally performed without any account of flow behavior. Here we present a methodology for characterizing fractured geologic reservoirs by integrating flow and seismic data. The key element of the proposed approach is the identification—within the inversion—of the intimate relation between fracture compliance and fracture transmissivity, which determine the acoustic and flow responses of a fractured reservoir, respectively. Owing to the strong (but highly uncertain) dependence of fracture transmissivity on fracture compliance, the modeled flow response in a fractured reservoir is highly sensitive to the geophysical interpretation. By means of synthetic models, we show that by incorporating flow data (well pressures and tracer breakthrough curves) into the inversion workflow, we can simultaneously reduce the error in the seismic interpretation and improve predictions of the reservoir flow dynamics. While the inversion results are robust with respect to noise in the data for this synthetic example, the applicability of the methodology remains to be tested for more complex synthetic models and field cases.Eni-MIT Energy Initiative Founding Member ProgramKorea (South). Ministry of Land, Transportation and Maritime Affairs (15AWMP-B066761-03

Software for evaluating probability-based integrity of reinforced concrete structures

In recent years, much research work has been carried out in order to obtain a more controlled durability and long-term performance of concrete structures in chloride containing environment. In particular, the development of new procedures for probability-based durability design has proved to give a more realistic basis for the analysis. Although there is still a lack of relevant data, this approach has been successfully applied to several new concrete structures, where requirements to a more controlled durability and service life have been specified. A probability-based durability analysis has also become an important and integral part of condition assessment of existing concrete structures in chloride containing environment. In order to facilitate the probability-based durability analysis, a software named DURACON has been developed, where the probabilistic approach is based on a Monte Carlo simulation. In the present paper, the software for the probability-based durability analysis is briefly described and used in order to demonstrate the importance of the various durability parameters affecting the durability of concrete structures in chloride containing environment

Earthquakes: from chemical alteration to mechanical rupture

In the standard rebound theory of earthquakes, elastic deformation energy is progressively stored in the crust until a threshold is reached at which it is suddenly released in an earthquake. We review three important paradoxes, the strain paradox, the stress paradox and the heat flow paradox, that are difficult to account for in this picture, either individually or when taken together. Resolutions of these paradoxes usually call for additional assumptions on the nature of the rupture process (such as novel modes of deformations and ruptures) prior to and/or during an earthquake, on the nature of the fault and on the effect of trapped fluids within the crust at seismogenic depths. We review the evidence for the essential importance of water and its interaction with the modes of deformations. Water is usually seen to have mainly the mechanical effect of decreasing the normal lithostatic stress in the fault core on one hand and to weaken rock materials via hydrolytic weakening and stress corrosion on the other hand. We also review the evidences that water plays a major role in the alteration of minerals subjected to finite strains into other structures in out-of-equilibrium conditions. This suggests novel exciting routes to understand what is an earthquake, that requires to develop a truly multidisciplinary approach involving mineral chemistry, geology, rupture mechanics and statistical physics.Comment: 44 pages, 1 figures, submitted to Physics Report

Introduction of an Electromagnetism Module in LS-DYNA for Coupled Mechanical Thermal Electromagnetic Simulations

A new electromagnetism module is being developed in LS-DYNA for coupled mechanical/thermal/electromagnetic simulations. One of the main applications of this module is Electromagnetic Metal Forming. The electromagnetic fields are solved using a Finite Element Method for the conductors coupled with a Boundary Element Method for the surrounding air/insulators. Both methods use elements based on discrete differential forms for improved accuracy. The physics, numerical methods and capabilities of this new module are presented in detail as well as its coupling with the mechanical and thermal solvers of LS-DYNA. This module is then illustrated on an Electromagnetic Metal Forming case