Abstract:We describe a methodology for quantitatively characterizing the fractured nature of a hydrocarbon or geothermal reservoir from surface seismic data under a Bayesian inference framework. The method combines different kinds of measurements of fracture properties to find a best-fitting model while providing estimates of the uncertainty of model parameters. Fractures provide pathways for fluid flow in a reservoir, and hence, knowledge about a reservoir's fractured nature can be used to enhance production from the reservoir. The fracture properties of interest in this study (to be inferred) are fracture orientation and excess compliance, where each of these properties are assumed to vary spatially over a 2-D horizontal grid which is assumed to represent the top of a reservoir. The Bayesian framework in which the inference problem is cast has the key benefits of (1) utilization of a prior model that allows geological information to be incorporated, (2) providing a straightforward means of incorporating all measurements (across the 2-D spatial grid) into the estimates at each gridpoint, (3) allowing different types of measurements to be combined under a single inference procedure and (4) providing a measure of uncertainty in the estimates. The observed data are taken from a 2-D array of surface seismic receivers responding to an array of surface sources. Well understood features from the seismic traces are extracted and treated as the observed data, namely the P-wave reflection amplitude variation with acquisition azimuth and offset (amplitude versus azimuth data) and fracture transfer function (FTF) data. Amplitude versus azimuth data are known to be more sensitive to fracture properties when the fracture spacing is significantly smaller than the seismic wavelength, whereas FTF data are more sensitive to fracture properties when the fracture spacing is on the order of the seismic wavelength. Combining these two measurements has the benefit of allowing inferences to be made about fracture properties over a larger range of fracture spacing than otherwise attainable. Geophysical forward models for the measurements are used to arrive at likelihood models for the data. The prior distribution for the fracture variables is obtained by defining a Markov random field over the lateral 2-D grid where we wish to obtain fracture properties, where this method for defining the prior has the added benefit of allowing for non-stationarity in the resulting model covariance. The fracture variables are then inferred by application of loopy belief propagation to yield approximations for the posterior marginal distributions of the fracture properties, as well as the maximum a posteriori and Bayes least-squares (posterior mean) estimates of these properties. Verification of the inference procedure is performed using a synthetic data set, where the estimates obtained are shown to be at or near ground truth for the full range of fracture spacings for fracture orientation and at low fracture spacings for excess compliance estimates.

Improved spatial understanding of induced seismicity hazard from the discretization of a curved fault surface

Induced seismicity potential based on probabilistic geomechanics for the San Juan Basin CarbonSAFE project

The Role of Fault-Zone Architectural Elements on Pore Pressure Propagation and Induced Seismicity

Probabilistic Assessment of Seismic Stability of a Rock Slope by Combining the Simulation of Stochastic Ground Motion with Permanent Displacement Analysis

Integration of rupture directivity models for the US National Seismic Hazard Model

Geomechanical Modeling of Fault Systems Using the Material Point Method – Application to the Estimation of Induced Seismicity Potential to Bolster Hydraulic Fracturing Social License

A rupture directivity adjustment model and its application in seismic hazard

A Bayesian Framework for Fracture Characterization from Surface Seismic Data

Separating Hydraulic Fracturing Microseismicity From Induced Seismicity by Bayesian Inference of Non‐Linear Pressure Diffusivity

Uncertainty quantification for fault slip inversion

New probabilistic fault displacement hazard models for reverse faulting

Fault rupture propagation through stratified sand–clay deposits and engineered earth structures: a meshfree and critical-state modeling approach

Bayesian Updating of Seismic Ground Failure Estimates via Causal Graphical Models and Satellite Imagery

Validating induced seismicity forecast models - Induced Seismicity Test Bench

Evaluation of the fault reactivation dynamics using fully coupled hydro-mechanical models

Mitigation and optimization of induced seismicity using physics-based forecasting

A simplified seismicity model of the bookshelf fault system of the Southwest Iceland transform zone

Tectonic Landform and Lithologic Age Impact Uncertainties in Fault Displacement Hazard Models

Coupling 3D geodynamics and dynamic earthquake rupture: fault geometry, rheology and stresses across timescales

Hazard implications of synchronous fault rupture: Case study of the Wasatch and West Valley faults, Utah

Coupled Boundary Element and Finite Volume Methods for Modeling Fluid-Induced Seismicity in Fault Networks within Low-Permeability Rocks