Size distribution of geological faults: Model choice and parameter estimation

Geological faults are important in reservoir characterization, since they influence fluid flow in the reservoir. Both the number of faults, or the fault intensity, and the fault sizes are of importance. Fault sizes are often represented by maximum displacements, which can be interpreted from seismic data. Owing to limitations in seismic resolution only faults of relatively large size can be observed, and the observations are biased. In order to make inference about the overall fault population, a proper model must be chosen for the fault size distribution. A fractal (Pareto) distribution is commonly used in geophysics literature, but the exponential distribution has also been suggested. In this work we compare the two models statistically. A Bayesian model is defined for the fault size distributions under the two competing models, where the prior distributions are given as the Pareto and the exponential pdfs, respectively, and the likelihood function describes the sampling errors associated with seismic fault observations. The Bayes factor is used as criterion for the model choice, and is estimated using MCMC sampling. The MCMC algorithm is constructed using pseudopriors to sample jointly the two models. The statistical procedure is applied to a fault size data set from the Gullfaks Field in the North Sea. For this data set we find that the fault sizes are best described by the exponential distribution.

[1]  Hilde Grude Borgos,et al.  Stochastic Modeling and Statistical Inference of Geological Fault Populations and Patterns , 2000 .

[2]  B. E. Shaw,et al.  Experimental evidence for different strain regimes of crack populations in a clay model , 1999 .

[3]  H. Fossen,et al.  Research article: The influence of seismic noise in structural interpretation of seismic attribute maps , 1997 .

[4]  T. Mattfeldt.  Stochastic Geometry and Its Applications , 1996 .

[5]  H. Fossen,et al.  Properties of fault populations in the Gullfaks Field, northern North Sea , 1996 .

[6]  I. Main,et al.  Special issue: Scaling laws for fault and fracture populations - Analyses and applications - Introduction , 1996 .

[7]  S. Chib.  Marginal Likelihood from the Gibbs Output , 1995 .

[8]  P. Green.  Reversible jump Markov chain Monte Carlo computation and Bayesian model determination , 1995 .

[9]  B. Carlin,et al.  Bayesian Model Choice Via Markov Chain Monte Carlo Methods , 1995 .

[10]  D. Sanderson,et al.  Sampling power-law distributions , 1995 .

[11]  W. Ryan,et al.  Quantitative fault studies on the East Pacific Rise : a comparison of sonar imaging techniques , 1994 .

[12]  C. Scholz,et al.  Growth of normal faults: Displacement-length scaling , 1993 .

[13]  Randall Marrett,et al.  Estimates of strain due to brittle faulting : sampling of fault populations , 1991 .

[14]  R. Hatcher.  Structural Geology: Principles, Concepts, and Problems , 1990 .

[15]  K. Heffer,et al.  Scaling Relationships in Natural Fractures: Data, Theory, and Application , 1990 .

[16]  D. Stoyan,et al.  Stochastic Geometry and Its Applications , 1989 .

[17]  John J. Walsh,et al.  Displacement gradients on fault surfaces , 1989 .

[18]  J. Walsh,et al.  Displacement Geometry in the Volume Containing a Single Normal Fault , 1987 .

[19]  J. J. Walsh,et al.  Distributions of cumulative displacement and seismic slip on a single normal fault surface , 1987 .

[20]  P. J. Green,et al.  Density Estimation for Statistics and Data Analysis , 1987 .

[21]  Benoit B. Mandelbrot,et al.  Fractal Geometry of Nature , 1984 .

[22]  W. K. Hastings,et al.  Monte Carlo Sampling Methods Using Markov Chains and Their Applications , 1970 .