Frequency-dependent anisotropy due to non-orthogonal sets of mesoscale fractures in porous media

Da Shuai,Alexey Stovas,Yang Zhao,Xingxing Huang,Lingyi He

DOI: https://doi.org/10.1093/gji/ggab341

IF: 3.352

2021-08-20

Geophysical Journal International

Abstract:SUMMARY The monoclinic medium with a horizontal symmetry plane is gradually being studied for seismic anisotropy characterization. The principle goal of this paper is to investigate the effect of fracture parameters (azimuth angle, density, aspect ratio, scale) on the exact and approximate monoclinic anisotropy parameters. We derive the monoclinic porous media based on the Chapman model which accounts for the wave-induced fluid flow and give the expressions of the Thomsen-style anisotropy parameters (nine orthorhombic anisotropy parameters: VP0, VS0, ε1, ε2, γ1, γ2, δ1, δ2, δ3, three exact monoclinic parameters: ζ1, ζ2, ζ3 and three approximate monoclinic parameters: $\widetilde{\zeta _{1}}, \widetilde{\zeta _{2}}, \widetilde{\zeta _{3}}$). The dependence of Thomsen-style anisotropy parameters associated with azimuth angle between two fracture sets is analysed. The orthorhombic anisotropy parameters and monoclinic anisotropy parameters have the same period (π) on the azimuth angle between two fracture sets. The exact and approximate monoclinic anisotropy parameters responsible for the rotation of the P-wave NMO ellipse have a similar trend versus the azimuth angle, while those responsible for the rotation of the S1- and S2-wave NMO ellipses have significant discriminations. The influence of fracture density, aspect ratio, and scale on the monoclinic parameters are also analysed. The monoclinic anisotropy parameters responsible for the rotation of the P-wave NMO ellipse decrease with fracture density and aspect ratio increasing from 0 to 0.1, while those responsible for the rotation of S1- and S2-wave NMO ellipses increase with the fracture parameters. The fracture density has a bigger influence on the monoclinic anisotropy parameters than the fracture aspect ratio. When saturated with different fluids (water and CO2), the monoclinic parameters have a similar behaviour versus the azimuth angle between two fracture sets.

geochemistry & geophysics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to explore the influence of fracture parameters (azimuth, density, aspect ratio, scale) on the exact and approximate expressions of monoclinic anisotropic parameters in porous media. Specifically, the author aims to study the anisotropic characteristics of porous media containing two sets of non - orthogonal meso - scale fractures at different frequencies, and takes into account the wave - induced fluid flow (WIFF) effect through the Chapman model. In addition, the author also compares the performance of the exact and approximate monoclinic anisotropic parameters under different fracture parameters to evaluate the impact of these parameters on seismic data processing and reservoir characterization.

### Main problem analysis

1. **Influence of fracture parameters on anisotropy**:
- The paper focuses on analyzing how the azimuth, density, aspect ratio, and scale of fractures affect monoclinic anisotropic parameters.
- The author derives the porous monoclinic medium based on the Chapman model and gives the exact and approximate expressions of monoclinic anisotropic parameters.

2. **Wave - induced fluid flow (WIFF) effect**:
- The wave - induced fluid flow effect is introduced, which is a factor that has been rarely considered in previous studies. The author believes that WIFF is the main cause of attenuation and dispersion in the seismic frequency band.

3. **Application of Thomsen - style anisotropic parameters**:
- Thomsen - style anisotropic parameters are used to describe transversely isotropic (VTI), orthotropic (ORT), and monoclinic anisotropic media.
- The differences between the approximate and exact monoclinic anisotropic parameters under different fracture parameters are compared to evaluate their impact on practical applications.

4. **Numerical simulation and verification**:
- Numerical examples are used to show the specific influence of fracture azimuth, density, aspect ratio, and scale on anisotropic parameters.
- The results show that fracture parameters significantly affect the changing trend of anisotropic parameters. In particular, under different azimuths, the approximate parameters underestimate the influence of azimuth on the NMO ellipse rotation of S1 and S2 waves.

### Summary of mathematical formulas

- **Fluid pressure change rate equation**:
\[
\frac{\partial m_f^{(1)}}{\partial t}=\rho_f^{(0)} k a_1^2 \frac{p^*-f_1}{\eta \ell}
\]
\[
\frac{\partial m_f^{(2)}}{\partial t}=\rho_f^{(0)} k a_2^2 \frac{p^*-f_2}{\eta \ell}
\]

- **Fluid pressure change rate equation in pores**:
\[
\frac{\partial m_p}{\partial t}=\frac{3\rho_f^{(0)} p_v}{4\mu}\left[(1 + K_p)\frac{\partial p^*}{\partial t}-\frac{1-\nu}{1 + \nu}\frac{\partial \sigma_{ll}}{\partial t}\right]
\]

- **Monoclinic anisotropic parameters**:
\[
\zeta_1=\frac{C_{16}(C_{33}-C_{55})-C_{36}(C_{13}+C_{55})}{2C_{55}(C_{33}-C_{55})}
\]
\[
\zeta_2=\frac{C_{26}(C_{33}-C_{44})-C_{36}(C_{23}+C_{44})}{2C_{44}(C_{33}-C_{44})}
\]
\[
\zeta_3=\frac{C_{36}}{2C_{33}\left(\frac{C_{13}+C_{55}}{C_{33}-C_{55}}+\frac{C_{23}+C_{44}}{C_{33}-C_{44}}\right)}