Figure 5.8 Fluid substitution modeling of the impact of fracture compliance using isotropic Gassmann models. Compare with the experimental results in Figure 5.5a and b. For (a) the core‐parallel fracture case, fluid substitution is assumed to happen only within the fracture. For (b) the core‐perpendicular case, it is assumed for the entire core. For better agreement between the models and the experiment, the shear modulus was slightly increased for the both cases, and the results are shown in dotted lines.

      In Figure 5.8, qualitatively consistent with the experiment, large η _D resulted in Young's modulus values that did not depend upon the presence of scCO₂ in the fracture. The model, however, does not show good quantitative agreement with the experiment, mostly underestimating the effect of the fracture and scCO₂ injection. This may be attributed to the additional effect of compliant interfaces between the sample and the metal bars as mentioned earlier, but also to the use of an isotropic and homogeneous Gassmann model for representing possibly anisotropic rock. Particularly, small errors in the isotropic shear modulus can have a large impact on the fluid substitution. In Figure 5.8a, in broken lines, examples are shown for the case with a shear modulus, which is ~2% larger than the experimental value, resulting in better agreement of the Young's modulus values for the intact sample.

      5.4.2. Frequency‐Dependent Fluid Pressure Diffusion Effect for Core‐Perpendicular Fracture Cases

      Compared with the core‐parallel fracture case, the modeling results of the core‐perpendicular case are qualitatively different from the experiment (Fig. 5.8b). Although the impact of the fluid substitution on the Young's modulus (through changes in the drained normal fracture compliance and the storage fracture compliance) is clearly seen, the model predicts much smoother and more monotonic changes than the experiment. (Note that, in Fig. 5.8b, similar to the core‐parallel case, assuming a larger shear modulus by 4.5% results in better predictions for the intact core case.) This is because equation 5.10 assumes that reductions in the fluid bulk modulus by scCO₂ injection occur throughout the sample (i.e., quasi‐static assumption). In reality, because of the finite frequency of the experiment, the effect of bulk fluid modulus reduction is limited by the fluid pressure diffusion length (the relative amplitude decays to e ^–1 over the distance) in the water‐saturated part of the rock. This is given approximately by (e.g., Pride, 2003)

      (5.12)

where k ₀ is the static permeability of the rock (34.5 md) and η _f is the fluid viscosity (0.47 cP). For a frequency of 1.5 kHz and the material properties in Table 5.2, δ _D is computed as 12.5 mm (Carbon Tan #1 core) and 13.6 mm (Carbon Tan #2 core).

      Attenuation of the Young's modulus appears to be more strongly affected by the finite frequency effect than the stiffness (the real part of the moduli), as seen in the increasing attenuation during the experiment (Fig. 5.5 a and c). Particularly, for the core‐perpendicular fracture case (Fig. 5.5 b and d), the frequency effect is quite dramatic. Figure 5.9 shows the observed changes during scCO₂ injection into Frac IIa and Frac IIb samples, extracted from Figure 5.5b and d and correlated to the X‐ray CT images of scCO₂. In contrast to the behavior predicted by the Gassmann model, the largest changes in the Young's modulus occur once scCO₂ reaches the fracture (pore space CO₂ saturation of the core ~2% for Frac IIa and ~5% for Frac IIb), corresponding to a peak and a subsequent rapid drop in the attenuation. Attenuation exhibits secondary increases as scCO₂ invades more of the fracture and the rock matrix.

Table 5.2 Modeling parameters

Parameter	Symbol	Value
Undrained frame Young's modulus	E U	20.5 GPa
Matrix frame shear modulus	G	7.60 GPa
Matrix Biot‐Willis effective stress coefficient	α	0.769
Matrix porosity	φ	16.7%
Mineral bulk modulus	K s	38 GPa
Mineral density	ρ s	2,650 kg/m3
Water density		989 kg/m3
scCO2 density		535 kg/m3
Water bulk modulus		2.46 GPa
scCO2 bulk modulus		0.04 GPa
Sample length	H	10.0 cm
Sample diameter	a	3.81 cm
Fracture aperture (core parallel)	h	0.54 mm