Kuster-Toksoz
Learning objectives
- Read Kuster-Toksoz as the first quantitative inclusion model: long-wavelength scattering off ellipsoidal pores, assembled through Berryman's polarization factors P and Q
- Use the sphere special case as the sanity anchor: spherical inclusions make Kuster-Toksoz reproduce the Hashin-Shtrikman bound exactly, K about 59.79 GPa for 10 percent water in calcite
- State the honesty of the model: inclusions are treated as non-interacting, so it holds only at modest concentrations
- See the dilute failure quantitatively: at 30 percent dry porosity with aspect ratio 0.15, Kuster-Toksoz gives K about 6.75 GPa where the self-consistent differential model gives 10.87
Scattering Off a Pore
The last section argued that pore shape controls compliance. Kuster and Toksoz turned that argument into a formula. Their idea is to treat a long seismic wave passing through the rock as a wave scattering off each pore. A pore is a soft (or empty) ellipsoid embedded in stiff mineral, and it perturbs the passing wavefield in a way that depends on its shape. For a wavelength much longer than the pore, that perturbation can be captured by two numbers per pore shape: a bulk polarization factor and a shear polarization factor . Berryman worked out these factors for an oblate spheroid of any aspect ratio, and they are the machinery this book uses. A round pore and a flat crack scatter the wave differently, so they carry different and , and that is exactly where the aspect ratio enters the effective modulus.
The Sphere Sanity Check
Any inclusion model has one answer it must get right. Make the pores perfect spheres and the model should return the stiffest possible arrangement of that much soft material in that much mineral, which is the Hashin-Shtrikman bound of Part 4. Kuster-Toksoz passes this exactly. Put 10 percent water-filled spherical porosity in calcite and the model returns a bulk modulus of 59.79 GPa, and the Hashin-Shtrikman bound for the same mixture is 59.79 GPa to the last digit. That identity is the anchor: whatever else the model does, at the round-pore limit it collapses onto a bound we already trust. From there, flattening the pores only softens the rock, tracing the ladder of the previous section.
The Price of Being Dilute
Kuster-Toksoz buys its simplicity with one strong assumption: every pore scatters the wave as if it were alone in the mineral, feeling nothing from its neighbors. That is fine when pores are sparse, but real reservoir porosity is not sparse, and when pores crowd together they interact, each sitting in mineral already softened by the others. The model does not know this, so it over-softens. Watch it fail: for dry pores of aspect ratio 0.15 in calcite at 30 percent porosity, Kuster-Toksoz gives a dry bulk modulus of 6.75 GPa, while the differential effective medium model of the next section, which does account for the crowding, gives 10.87. The two agree closely at low porosity (both near 54 GPa at 5 percent) and diverge as porosity climbs, precisely because the non-interacting assumption breaks down. The honest rule is to keep Kuster-Toksoz to modest concentrations and reach for a self-consistent method past them. The next section builds that method by adding the pores a pinch at a time.
References
- Kuster, G. T., & Toksoz, M. N. (1974). Velocity and attenuation of seismic waves in two-phase media: Part I. Theoretical formulations. Geophysics, 39(5), 587-606.
- Berryman, J. G. (1980). Long-wavelength propagation in composite elastic media II. Ellipsoidal inclusions. Journal of the Acoustical Society of America, 68(6), 1820-1831.