Computational Techniques for Efficient Solution of Discretized Biot's Theory for Fluid Flow in Deformable Porous Media

In soil and rock mechanics, coupling effects between geomechanics field and fluid-flow field are important to understand many physical phenomena. Coupling effects in fluid-saturated porous media comes from the interaction between the geomechanics field and the fluid flow. Stresses subjected on the porous material result volumetric strains and fluid diffusion in the pores. In turn, pore pressure change cause effective stresses change that leads to the deformation of the geomechanics field. Coupling effects have been neglected in traditional geotechnical engineering and petroleum engineering however, it should not be ignored or simplified to increases reliability of the results. The coupling effect in porous media was theoretically established in the poroelasticity theory developed by Biot, and it has become a powerful theory for modeling three-dimensional consolidation type of problem. The analysis of the porous media with fully-coupled simulations based on the Biot's theory requires intensive computational effort due to the large number of interacting fields. Therefore, advanced computational techniques need to be exploited to reduce computational time. In order to solve the coupled problem, several techniques are currently available such as one-way coupling, partial-coupling, and full-coupling. The fully-coupled approach is the most rigorous approach and produces the most correct results. However, it needs large computational efforts because it solves the geomechanics and the fluid-flow unknowns simultaneously and monolithically. In order to overcome this limitation, staggered solution based on the Biot's theory is proposed and implemented using a modular approach. In this thesis, Biot's equations are implemented using a Finite Element method and/or Finite Difference method with expansion of nonlinear stress-strain constitutive relation and multi-phase fluid flow. Fully-coupled effects are achieved by updating the compressibility matrix and by using an additional source term in the conventional fluid flow equation. The proposed method is tested in multi-phase FE and FD fluid flow codes coupled with a FE geomechanical code and numerical results are compared with analytical solutions and published results.

Name thesis:degree_name

Ph. D.

Level thesis:degree_level

doctoral

Discipline thesis:degree_discipline

Civil Engineering

Department dc:contributor.department

Civil Engineering

Grantor dc:publisher

Virginia Tech

Year dc:date.issued

2008

Author dc:creator

• Lee, Im Soo

Chair dc:contributor.committeechair

• Gutierrez, Marte S.

Committee members dc:contributor.committeemember

• Filz, George M.
• Duncan, James Michael
• Chin, Lee

• Biot's theory
• fully-coupled solution
• finite elements

Statement dc:rights

• In Copyright

Licence dc:rights.uri

• http://rightsstatements.org/vocab/InC/1.0/

Dc Identifier Other

etd-07012008-115136

OAI identifier oai:identifier

oai:vtechworks.lib.vt.edu:10919/28179