Mechanics of stratified anisotropic poroelastic media |

where is the symmetric compliance matrix. The numbers 1,2,3 always indicate Cartesian axes (say, , , respectively). The -direction is usually chosen as the layering direction, which could be oriented any direction in the earth. But, in many geological and geophysical applications, the 3-axis (or -axis) is also taken to be the vertical direction, and I conform to this convention here. The principal stresses are , , , in the directions 1,2,3, respectively. Similarly, the principal strains are , , . The stresses , , are the torsional shear stresses, associated with rotation-based strains around the 1, 2, or 3 axes, respectively. The corresponding torsional strains are , , and , where the torsional motion is again a rotational straining motion around the 1, 2, or 3 axes. The compliance matrix is symmetric, so , and this fact could have been used when displaying the matrix. The axis pairs in the subscripts , , , , , and , are often labelled (again following the conventions of Voigt) as 1,2,3,4,5,6, respectively.

The important contribution made by Backus (1962) [also see Postma (1955)] is the observation that, in a horizontally layered system, there are certain strains and stresses that are necessarily continuous across boundaries between layers, while the others are not necessarily continuous. I have been implicitly (and now explicitly by calling this fact out) assuming that the interfaces between layers are in welded contact, which means practically that the in-plane strains are always continuous: so if axis 3 (or ) is the symmetry axis (as is most often chosen for our layering problem), I have , , and are all continuous. Similarly, in welded contact, I must have continuity of the all the stresses involving the 3 (or ) direction: so , , and must all be continuous.

Then, following Backus (1962) and/or Schoenberg and Muir (1989) but -- for present purposes considering instead the compliance (inverse of stiffness) matrix -- I have rearranged the statement of the problem so that:

Note that this equation, although similar to (40) is nevertheless quite different because of the rearrangement of the matrix elements and the reordering of the strains and stresses. The expression in (41) is general for all elastic media. In the main text I restrict the discussion to orthotropic media. Assuming then that I am using the correct set of axes as the symmetry axes in the presentation, all

Expression of (41) can be made more compact by writing it as:

where

and

with (with superscript indicating the matrix transpose). Also I have

and

It is important to distinguish between ``slow'' and ``fast'' variables
in this analysis, since this distinction makes it clear when
and how averaging should be performed.
The ``slow'' variables, *i.e.*, those that are continuous across the (here assumed horizontal)
boundaries and also essentially constant for the present quasi-static application, are
those contained in
and
. So, after averaging
along the layering direction, I have:

where , and all the starred quantities are the

Of the three final results, the two easiest ones to compute are:

where is the layer average of some quantity. These results follow from this equation:

which also followed immediately from the formula

multiplying through first by the inverse of , and then performing the layer average. [Note that and are both normally square and invertible matrices, whereas for most systems the off-diagonal matrix is not invertible. But, this fact does not cause problems in the analysis, because I do not need to invert in order to solve the averaging problem at hand.] These averages are meaningful because, when the matrix equations presented are multiplied out, there never appear any cross products of two quantities that are both unknown. [From this view point, Eq. (51) is an equation for , just as the unaveraged version of (51) is an equation for in each layer.] So simple layer-averaging suffices (thereby providing the main motivation and value of this method). Multiplying (51) through by then gives the results (49) and (50).

The remaining result is more tedious to compute, since it requires several intermediate steps in its derivation. But the final result is given by the formula:

To provide some clues to the derivation, again consider:

which is just a rearrangement of (52). The point is that is then given immediately in terms of the quantities and , which are both ``slow'' variables and therefore essentially constant. An intermediate result that helps to explain the form of this relation (53) is:

Substituting for from (54) into

and then averaging, I find that

an expression that completely determines the remaining coefficients. After some more algebra, the formula giving the final result is:

Equation (58) contains all the information needed to produce the third and final result found in (53).

Another check on these formulas is to compare them directly to those found by Schoenberg and Muir (1989). However, direct comparison is not so easy, since their analysis focuses on the stiffness version of these equations. My treatment makes use of the compliance version instead. Since the symmetries of the two forms of the equations nevertheless are nearly identical, cross-checks and comparisons will be left to the interested reader.

Mechanics of stratified anisotropic poroelastic media |

2010-05-19