Numerical Modeling Of Coupled Seismic And Electromagneticwaves In Fluid-Saturated Porous Media
Abstract
This work presents a finite element procedure for the numerical approximation of the electroseismic
problem at seismic frequencies. The electroseismic equations used are the ones derived by S.
Pride, consisting in the coupled Biot’s equations of motion and Maxwell’s equations where the seismoelectric
feedback is being neglected. The modeled domain comprises two half-spaces, one being air, and
the other one a horizontally layered medium. The chosen electromagnetic source is an infinite plane of
time dependent electric current located above the surface of the Earth. Under these two assumptions, the
electric and magnetic fields and the solid and fluid displacements depend only on one coordinate, namely
the one chosen to describe the vertical variations of the Earth; therefore the problem can be considered
to be one dimensional.
The existence of a unique solution for both the continuous and discrete weak problems is analyzed, considering
a finite computational domain by recoursing to absorbing boundary conditions.
The finite element procedure is carried out by using C0 linear functions for the electric field and the
solid displacements, and piecewise constant functions for the magnetic field and fluid displacements,
respectively. A synthetic example showing the capabilities of the numerical method is presented.
problem at seismic frequencies. The electroseismic equations used are the ones derived by S.
Pride, consisting in the coupled Biot’s equations of motion and Maxwell’s equations where the seismoelectric
feedback is being neglected. The modeled domain comprises two half-spaces, one being air, and
the other one a horizontally layered medium. The chosen electromagnetic source is an infinite plane of
time dependent electric current located above the surface of the Earth. Under these two assumptions, the
electric and magnetic fields and the solid and fluid displacements depend only on one coordinate, namely
the one chosen to describe the vertical variations of the Earth; therefore the problem can be considered
to be one dimensional.
The existence of a unique solution for both the continuous and discrete weak problems is analyzed, considering
a finite computational domain by recoursing to absorbing boundary conditions.
The finite element procedure is carried out by using C0 linear functions for the electric field and the
solid displacements, and piecewise constant functions for the magnetic field and fluid displacements,
respectively. A synthetic example showing the capabilities of the numerical method is presented.
Full Text:
PDFAsociación Argentina de Mecánica Computacional
Güemes 3450
S3000GLN Santa Fe, Argentina
Phone: 54-342-4511594 / 4511595 Int. 1006
Fax: 54-342-4511169
E-mail: amca(at)santafe-conicet.gov.ar
ISSN 2591-3522