Finite difference modeling of Biot's poroelastic equations atseismic frequencies
Across the seismic band of frequencies (loosely defined as<10 kHz), a seismic wave propagating through a porous material willcreate flow in the pore space that is laminar; that is, in thislow-frequency "seismic limit," the development of viscous boundary layersin the pores need not be modeled. An explicit time steppingstaggered-grid finite difference scheme is presented for solving Biot'sequations of poroelasticity in this low-frequency limit. A key part ofthis work is the establishment of rigorous stability conditions. It isdemonstrated that over a wide range of porous material properties typicalof sedimentary rock and despite the presenceof fluid pressure diffusion(Biot slow waves), the usual Courant condition governs the stability asif the problem involved purely elastic waves. The accuracy of the methodis demonstrated by comparing to exact analytical solutions for both fastcompressional waves and slow waves. Additional numerical modelingexamples are also presented.
- Research Organization:
- Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
- Sponsoring Organization:
- USDOE. Assistant Secretary for Fossil Energy.Petroleum
- DOE Contract Number:
- DE-AC02-05CH11231
- OSTI ID:
- 927822
- Report Number(s):
- LBNL-61829; JGREA2; R&D Project: G32801; BnR: AC1005000; TRN: US200816%%1074
- Journal Information:
- Journal of Geophysical Research, Vol. 111, Issue B10305; Related Information: Journal Publication Date: 10/14/2006; ISSN 0148-0227
- Country of Publication:
- United States
- Language:
- English
Similar Records
Finite-difference modeling of Biot's poroelastic equations across all frequencies
A hydrologic view on Biot's theory of poroelasticity