The normal-to-tangential fracture compliance ratio is usually used as a fracture fluid indicator (FFI) for fluid identification in fractured reservoirs. With a new parameterization for fracture weaknesses, we have defined a new FFI based on azimuthally anisotropic elastic impedance (EI) inversion and fractured anisotropic rock-physics models. First, we derived a new azimuthally anisotropic EI equation with a similar expression for the isotropic and anisotropic EI parts to remove the exponential correction of EI that is attributable to weak anisotropy. Then, we built a fractured anisotropic rock-physics model used for the estimation of well-log parameters for the normal and tangential fracture weaknesses, which built the initial background low-frequency trend of fracture weaknesses. Finally, based on the azimuthally anisotropic EI inversion method with the Cauchy-sparse and low-frequency information regularization, we estimated an FFI applied to fluid identification in fractured reservoirs. Tests on the synthetic and real data demonstrate that the anisotropic parameters related to fracture weaknesses can be estimated reasonably and stably and that our method appears to provide an alternative available for fluid identification in fractured reservoirs.

Horizontally and vertically aligned fractures (joints) permeated in an isotropic background can give rise to a long-wavelength equivalent orthorhombic medium, which is common for naturally fractured reservoirs. We have developed a feasible approach to characterize the horizontal and vertical fractures using observed azimuthal seismic reflection data. For weak anisotropy, we assume that the horizontal and vertical fractures are decoupled. Using a linear-slip orthorhombic model, we first obtain analytic expressions for the effective elastic stiffness matrix and the corresponding perturbed matrix of an effective orthorhombic anisotropic elastic medium. Combining the scattering function and the perturbed matrix of the orthorhombic medium, we then derive a linearized PP-wave reflection coefficient of effective orthorhombic medium in terms of compressional-wave (P-wave) and shear-wave (S-wave) moduli, density, and horizontal- and vertical-fracture-induced normal and tangential weaknesses. To handle the inverse problem for multiple model parameters in an orthorhombic medium, we reexpress the linearized PP-wave reflection coefficient as a Fourier series. According to the sensitivity analysis of Fourier coefficients (FCs) to isotropic background elastic parameters and four fracture weaknesses, we finally establish a three-step inversion approach to describe the orthorhombic model, which involves (1) estimation of the FCs at different incidence angles from observed azimuthal seismic data and determination of the symmetry axis azimuth, (2) an iterative Bayesian inversion for isotropic background elastic parameters and fracture weaknesses of horizontal fractures from the zeroth-order FC, and (3) an iterative Bayesian inversion for fracture weaknesses of vertical fractures from the second-order FC. The proposed approach is validated by tests on synthetic and field data sets, which demonstrates that the inversion results of P- and S-wave moduli, density, and four fracture weaknesses are robust and reasonable for gas-bearing fractured reservoir characterization with orthorhombic symmetry.

The characterization of fracture-induced tilted transverse isotropy (TTI) seems to be more suitable to actual scenarios of geophysical exploration for fractured reservoirs. Fracture weaknesses enable us to describe fracture-induced anisotropy. With the incident and reflected PP-wave in TTI media, we have adopted a robust method of azimuthal amplitude variation with offset (AVO) parameterization and inversion for fracture weaknesses in a fracture-induced reservoir with TTI symmetry. Combining the linear-slip model with the Bond transformation, we have derived the stiffness matrix of a dipping-fracture-induced TTI medium characterized by normal and tangential fracture weaknesses and a tilt angle. Integrating the first-order perturbations in the stiffness matrix of a TTI medium and scattering theory, we adopt a method of azimuthal AVO parameterization for PP-wave reflection coefficient for the case of a weak-contrast interface separating two homogeneous weakly anisotropic TTI layers. We then adopt an iterative inversion method by using the partially incidence-angle-stacked seismic data with different azimuths to estimate the fracture weaknesses of a TTI medium when the tilt angle is estimated based on the image well logs prior to the seismic inversion. Synthetic examples confirm that the fracture weaknesses of a TTI medium are stably estimated from the azimuthal seismic reflected amplitudes for the case of moderate noise. A field data example demonstrates that geologically reasonable results of fracture weaknesses can be determined when the tilt angle is treated as the prior information. We determine that the azimuthal AVO inversion approach provides an available tool for fracture prediction in a dipping-fracture-induced TTI reservoir.

Fluid identification and fracture discrimination play an important role in the exploration and development of an oil-bearing fractured reservoir. The most common fluid indicator in fractured reservoirs, the normal-to-shear fracture compliance ratio, is influenced by the fluid content and the fracture intensity. To reduce the ambiguities in the discrimination of the fluid and fracture parameters, we have aimed to extend the scattering theory to implement the fluid identification and fracture detection by incorporating the azimuthal data in an oil-bearing fractured reservoir via the proposed Bayesian amplitude variation with offset and azimuth (AVOAz) inversion approach. The background medium is, as far as the scattering theory is concerned, an isotropic medium without fractures, and the fractured medium is corresponding to a perturbed medium. The elastic parameters of a saturated anisotropic medium can be parameterized as a perturbation over a homogeneous isotropic background medium. We used the scattering theory to derive a generalized AVOAz approximation that provided the iterative estimates of hydrocarbon fluid indicator, shear modulus, density, and fracture weaknesses in a Bayesian scheme. The inversion algorithm is based on a convolutional model and a weak-contrast and small-weakness PP-wave reflection coefficient. The approach is applied to an oil-bearing field data set from a fractured marlstone reservoir. We observe that reasonable estimates of fluid indicator and fracture weaknesses are inverted, which can be used to perform the discrimination of fluid and fracture parameters. We conclude that the proposed approach provides us a potentially powerful tool to estimate the reservoir fluid and fracture properties in a more straightforward and efficient manner than those previous methods.