CPL - Chalmers Publication Library
| Utbildning | Forskning | Styrkeområden | Om Chalmers | In English In English Ej inloggad.

Numerical identification of a viscoelastic substitute model for heterogeneous poroelastic media by a reduced order homogenization approach

R. Janicke ; Fredrik Larsson (Institutionen för tillämpad mekanik, Material- och beräkningsmekanik) ; Kenneth Runesson (Institutionen för tillämpad mekanik, Material- och beräkningsmekanik) ; H. Steeb
Computer Methods in Applied Mechanics and Engineering (0045-7825). Vol. 298 (2016), p. 108-120.
[Artikel, refereegranskad vetenskaplig]

The paper deals with the computational homogenization of pressure diffusion processes in a poroelastic medium. The underlying physical phenomena are of interest for the interpretation of seismic data with applications in hydrocarbon production and geothermal energy. Pressure diffusion is assumed to take place on a length scale much smaller than the observer scale. Thus, the macroscopic observer is not able to measure the properties of the poroelastic medium directly but notices an intrinsic viscous attenuation. Under these circumstances, the macro-scale can be interpreted as a single-phase solid with (apparent) viscoelastic properties. In this paper, we establish a numerical upscaling procedure based on a volume averaging concept. This enables us to identify the material properties of the viscoelastic substitute model in a numerically efficient manner. For this purpose, the poroelastic medium on the small scale is modelled in terms of the momentum balance of the biphasic mixture and a coupled diffusion equation. We approximate the poroelastic pressure field on the small scale by a linear combination of pressure modes forming a reduced orthogonal basis and being identified by a Proper Orthogonal Decomposition (POD) technique. From the superposition principle, the evaluation of the poroelastic continuity equation results in a proper identification of the evolution equations defining the apparent viscoelastic model. In comparison to the nested FE2 solution schemes, the reduced order approach only requires a small amount of "off-line" precomputations and, therefore, causes very low numerical costs. The proposed method is validated for the simple setup of a layered porous rock with alternating water-and gas-saturated zones.

Nyckelord: Computational homogenization, Poroelasticity, Order reduction,

Denna post skapades 2016-01-20.
CPL Pubid: 231009


Läs direkt!

Länk till annan sajt (kan kräva inloggning)

Institutioner (Chalmers)

Institutionen för tillämpad mekanik, Material- och beräkningsmekanik (2005-2017)


Teknisk mekanik

Chalmers infrastruktur