Next: EQUATIONS OF MOTION WITH Up: REVIEW OF AVERAGING THEOREM Previous: Dynamical equations

## Quasistatic constitutive relations for isotropic materials

Taking the trace of equations (solidaveragestress) and (fluidaveragestress) gives the following results. The constitutive relations for dilatations and porosity are

- p_sK_m = B<>u_s - 1-   and

- p_fK_f = u_f + ,   where the partial averages and are related to the full volume averages by and .One assumption implicit in (solids) and (fluidf) is that changes much more slowly in space than the displacement variables and . This assumption allows us to remove the factors involving the porosity from the divergence terms, but we will show in a later publication that this assumption is not crucial to our analysis.

It is important at this point to understand the interpretations of all the symbols appearing in the last two equations. First, the variable is just the change in the average fluid pressure throughout the fluid phase. The change in average solid pressure is related to the macroscopic confining pressure change by the averaging relation .Thus, is just the average change in solid pressure experienced by the solid. Since and may be viewed as the pressures we can control, is the (solid volume) weighted average of the confining pressure after subtracting that part of the confining pressure supported by the fluid pressure. The change in porosity is given by .The porosity change occurs naturally in these expressions because

= 1V_I_f B<>n_fu_fdS = -1V_I_s B<>n_su_sdS.   This result is demonstrated in the next section.

The divergence of the average solid displacement is properly interpreted as the dilatation of the porous solid frame (not the dilitation of the solid alone). This interpretation is not obvious, but it follows from the fact that the term arises from the external surface integral [c.f., Eq.(gradofave)]

<u_s> = 1V _E_s B<>n_s u_sdS,   which is exactly the surface integral needed to define the overall behavior of the porous solid frame. Thus, in terms of the definitions of Brown and Korringa (1975),

B<>u_s = VV = e = - p_dK^* - p_fK_m.   This interpretation is the same one reached by Pride et al. (1992) using a combination of the standard thought experiments (jacketed and unjacketed) of Biot and Willis (1957). To check that this is so, we can easily show that

1- = - [1-K^* - 1K_m] (p_s - p_f),   using either approach when a single constituent is present so that the Brown and Korringa unjacketed constants satisfy .Note that (porositycheck) can also be written as

= - [1-K^*-1K_m]p_d,   emphasizing that porosity is constant if differential pressure is constant -- a general result for microhomogeneous porous frames, but not true otherwise. Thus, the left hand side of (solids) is just the solid dilatation , while the two terms on the right hand side are .

Similarly, it is important to understand that the expression is not just a fluid dilatation, but also includes the effects of fluid motion in and out of the volume. In fact, this is already apparent from (fluidf) since the strict fluid dilatation satisfies

-V_fV_f = p_fK_f,   yet (fluidf) contains an additional term related to changes in porosity. The correct physical interpretation of is provided by its relation to the increment of fluid content

= (B<>u_s - B<>u_f),   where is defined as

V_-V_fV = (VV - V_fV_f) +   and has the interpretation (Biot, 1973; Berryman and Thigpen, 1985) of the relative change in fluid mass per unit volume of initial fluid mass. Note that (fluiddil) and (zetadefined) are in agreement with (zetarelation) if the averaging equation (fluidf) is also satisfied.

The equations (solids), (fluidf), and (BandKe) are sufficient to arrive at the standard form of the equations relating e and to the macroscopic pressures and for a single constituent porous medium given by

e - = 1/K^* & 1/K_m - 1/K^* 1/K_m - 1/K^* & 1/K^* + /K_f - (1+)/K_m -p_c - p_f .   These equations are completely consistent with the results of Pride et al. (1992) as can be demonstrated by substituting the definitions given above into their formulas (48) and (49), and then doing a straightforward (though somewhat tedious) matrix inversion. The relationship between these matrix elements and the coefficients H, C, and M in Biot's equations are given by Berryman (1992a). They are: ,, and , where the Biot-Willis parameter is defined by .

Next: EQUATIONS OF MOTION WITH Up: REVIEW OF AVERAGING THEOREM Previous: Dynamical equations
Stanford Exploration Project
11/12/1997