A mixed finite-element method for solving the poroelastic Biot equations with electrokinetic coupling

Abstract

SUMMARY In this paper, a numerical method for solving the Biot poroelastic equations is developed. These equations comprise acoustic (typically water) and elastic (porous medium frame) equations, which are coupled mainly through fluid/solid drag terms. This wave solution is coupled to a simplified form of Maxwell’s equations, which is solved for the streaming potential resulting from electrokinesis. The ultimate aim is to use the generated electrical signals to provide porosity, permeability and other information about the formation surrounding a borehole. The electrical signals are generated through electrokinesis by seismic waves causing movement of the fluid through pores or fractures of a porous medium. The focus of this paper is the numerical solution of the Biot equations in displacement form, which is achieved using a mixed finite-element formulation with a different finite-element representation for displacements and stresses. The mixed formulation is used in order to reduce spurious displacement modes and fluid shear waves in the numerical solutions. These equations are solved in the time domain using an implicit unconditionally stable time-stepping method using iterative solution methods amenable to solving large systems of equations. The resulting model is embodied in the MODELLING OF ACOUSTICS, POROELASTICS AND ELECTROKINETICS (MAPEK) computer model for electroseismic analysis.