Wave Scattering in Porous Media and Shape Reconstruction of Scatterers
Wave Scattering in Porous Media and Shape Reconstruction of Scatterers
Wave Scattering in Porous Media and Shape Reconstruction of Scatterers
Koji Yamamoto
Japan Oil, Gas and Metals National Corporation, Chiba, Japan
Masaki Yamamda
Department of Civil Engineering, Tohoku University, Sendai, Japan
ABSTRACT: The Linearized Inverse Scattering Analysis (LISA) with Born approximation is utilized for a
shape reconstruction of inhomogeneties in fluid filed porous materials using elastic waves. The mechanical
characters of the media are modeled with Biots theory of poroelasticity. Integral expression of the scattered
field is the basis of this method, and inverse Fourier transform of the approximate form of the equations gives the
geometry of the scatterer from the observable values of scattered amplitude. Numerical examples of the shape
reconstruction that are derived from simulated scattered amplitude computed by the boundary integral equation.
In the examples, images of the scatterers such as a cavity in porous media or a domain in which pore fluid is
different from surrounding region are presented.
235
Copyright 2005 Taylor & Francis Group plc, London, UK
D (,,,f /m,1/M,1/m)
x
Dc (,, + , + ,
f /m + (f /m)1/M + (1/M), 1/m + (1/m))
q(s)
q(i)
p(i) where ui , p is the total field that is the sum of the
p(s) incident and the scattered waves.
y By applying the far field approximation of the
fundamental solution (Yamamoto & Kitahara 2003):
Figure 1. Wave scattering by a scatterer Dc in a poroelastic
medium.
236
Copyright 2005 Taylor & Francis Group plc, London, UK
unknown value of the righthand side is solely the shape
of the scatterer (x), and the scattering amplitude is a Transmitter-
receiver
linear function of (x). It is noted that the Born approx- 3a
imation is well established in the low-frequency range 2a
in which wavelengths of incident and scattered waves
is longer than the typical length of the scatterer. Flat ellipsoid Scatterer
Let us multiply both sides of Equation 6 by the par-
ticle displacement vector of the scattering wave dK ,
and the equation is rewritten as 3a
2a
Long ellipsoid
(s) (i)
where, sp = kw(s) p kw(i) p . Figure 2. Schematics of a Pulse-Echo mode measurement.
Since Equation 8 has a form of Fourier transform,
F(k) = f (x)exp(ik x)dx, the inverse Fourier For the scatterer, solely Biots modulus M (Biot &
transform f (x) = 1/(2)3 F(x) exp(ik x)dk gives Willis 1957) is different from the matrix. Among four
the shape function (x). In the case of Biot material, moduli of poroelastic system, only the modulus is
the wavenumber vector s is complex, thus the follow- affected by the compressibility of the pore fluid, and
ing modified form of the inverse Fourier transform the effect is written as follows (Stoll & Kan 1981):
is used:
237
Copyright 2005 Taylor & Francis Group plc, London, UK
(x) (x)
x3
x3
x1
x1
(a) Distribution of (x) in x1-x3
(a) Distribution of (x) in x1-x3
x3
x3 x2
x2 x1
x1
(b) Isopleth surface that (x) is 50% of its peak value
238
Copyright 2005 Taylor & Francis Group plc, London, UK
(x) 0.2- 4.0krad/sec
0.2- 8.0krad/sec
0.2-12.0krad/sec
0.2-16.0krad/sec
x3
0.2-20.0krad/sec
Reference
x1
(x)
(a) Distribution of (x) in x1-x3
x(y = 0, z = 0)
L2-wave
4 CONCLUSION
239
Copyright 2005 Taylor & Francis Group plc, London, UK
REFERENCES Nur, A. M. and Wang, Z. 1987. In situ seismic monitor-
ing EOR: The petrophysical basis. SPE16865. Proc. 62nd
Biot, M. A. 1956. Theory of propagation of elastic waves AnnualTechnical Conference and Exhibition of the Society
in a fluid-saturated porous solid. I. Lowfrequency range. of Petroleum Engineers. Dallas, TX. 307314.
J. Acoust. Soc. Am. 28: 168178. Ohsita, T. 2002. Utilization of depleted gas/oil reservoirs. J.
Biot, M. A. and Willis, D. G. 1957. The elastic coefficients of Japanese Asoc. Petro. Tech. 67: 538546. (in Japanese).
the theory of consolidation. J . Appl. Mech. 24: 594601. Rose, J. H. and Opsal, J. L. 1982. Inversion of ultra-
Biot, M. A. 1962. Mechanics of deformation and acous- sonic scattering data. Review of Progress in Quantitative
tic propagation in porous media. J. Appl. Phys. 33: Nondestructive Evaluation. 1: 187195.
14821498. Rose, J. H. 1984. Exterior reconstruction of a threedimen-
Chadwick, O., Zweigel, P., Gregersen, U., Kirby, G. A., sional scatterer, Wave Motion. 8: 149154.
Holloway, S. and Jhannessen, P. N. 2002. Geological Schmerr, L.W., Sedov,A. and Chiou, C.P. 1989.A unified con-
characterization of CO2 storage site: Lessons from Sleip- strained inversion model for ultrasonic flaw sizing. Res.
ner, northern North sea. Proc. 6th Green House Gas Nondestr. Eval. 1: 7797.
Technology Conference. Kyoto, Japan: B13. Stoll, R. D. and Kan, T.-K. 1981. :flection of acoustic waves
Cheng, A. H.-D., Badmus, T. and Beskos, D. E. 1991. Inte- at a water-sediment interface. J. Acoust. Soc. Am.. 70:
gral equation for dynamic poroelasticity in frequency 149156.
domain with BEM solution. J. Eng. Mech., ASCE. 117: Takahashi, H., Yonezawa, T. and Fercho, E. 2003. Opera-
11361157. tion overview of the 2002 Malik gas hydrate production
Detournay, E. and Cheng, A. H.-D. 1993. Fundamental poroe- research well program at the Machenzie Delta in the Cana-
lasticity. Comprehensive Rock Engineering: Principles, dian Arctic. OTC15124. Proc. 2003 Offshore Technology
Practice, and Projects, Vol. 2 Analysis and Design Meth- Conference. Houston, TX.
ods. C. Fairhurst eds. Pergamon Press, Oxford. 113171. Yamamoto, K. and Kitahara, M. 2003. Characterization of
Dominguez, J. 1992. Boundary element approach for scatterers in poroelastic media using far field integral
dynamic poroelastic problems. Int. J. Numer. Meth. Eng. representation. Proc. 16th Engineering Mechanics Con-
35: 307324. ference, ASCE, July 1618, 2003. Seattle, Washington.
Fukui, T., Funato, K. and Inoue, K. 1996. Frequency domain Yamamoto, K. and Kitahara, M. 2004. A Numerical method
boundary element method for wave propagation in Biot for wave scattering in poroelastic media. Structural
material. J. Boundary Element Methods, JASCOME. 13: Eng./Earthquake Eng., JSCE. Vol.21, No.2: 143s158s.
149152. (in Japanese).
Kitahara, M., Nakahata, K. and Hirose, S. 2002. Elast-
dynamic inversion for shape reconstruction and type
classification of flaws. Wave Motion, 36, 443455.
240
Copyright 2005 Taylor & Francis Group plc, London, UK