The Impact of Biased Sampling on the Estimation of the Semivariogram Within Fractured Media Containing Multiple Fracture Sets

Monte Carlo simulation was used to examine the error (statistical bias) introduced in estimating a sample semivariogram through application of oriented sampling patterns to variables which are correlated with fracture orientation. Sample semivariograms of the directional components of the water velocity were used to illustrate that oriented sampling schemes can provide biased data sets which result in error in the estimation of the semivariogram, particularly in the estimation of the sill (or variance). Three sampling patterns were used to analyze directional semivariograms of the components of the fluid velocity: sampling along lines parallel to the mean regional hydraulic gradient, sampling among lines perpendicular to the mean regional hydraulic gradient, and sampling along fracture segments. The first two sampling patterns were shown to introduce substantial error in the sills of the velocity variograms. It is argued that this error is due to the combination of unequal sampling of fractures with different orientations (i.e., sampling bias) and the systematic variation in the magnitude of the velocity components with orientation of the fracture. As a consequence, it is suggested that correction factors developed to correct fracture frequency statistics need to be extended to improve estimation of spatial moments of variables which are correlated with fracture orientation.

[1]  D. T. Snow.  ROCK FRACTURE SPACINGS, OPENINGS, AND POROSITIES , 1968 .

[2]  R. D. Terzaghi.  Sources of Error in Joint Surveys , 1965 .

[3]  Design of a Piezometer Network for estimation of the sample variogram of the hydraulic gradient: The role of the instrument , 1997 .

[4]  J. M. Shafer,et al.  Approximation of confidence limits on sample semivariograms from single realizations of spatially correlated random fields. , 1990 .

[5]  Peter K. Kitanidis,et al.  Subsurface Flow and Transport: Groundwater flow in heterogeneous formations , 1997 .

[6]  T. Clemo.  Dual permeability modeling of fractured media , 1994 .

[7]  Yoram Rubin,et al.  Stochastic analysis of boundaries effects on head spatial variability in heterogeneous aquifers: 2. Impervious boundary , 1989 .

[8]  M. Lipschutz,et al.  Meteoroid streams: Evidence from meteorites recovered on Earth , 1995 .

[9]  G. Dagan,et al.  Stochastic analysis of boundaries effects on head spatial variability in heterogeneous aquifers: 1. Constant head boundary , 1988 .

[10]  Marc G. Genton,et al.  Highly Robust Variogram Estimation , 1998 .

[11]  Mei-Cheng Wang,et al.  Hazards regression analysis for length-biased data , 1996 .

[12]  W. Blanckenhorn,et al.  CORRECTING FOR SAMPLING BIAS IN QUANTITATIVE MEASURES OF SELECTION WHEN FITNESS IS DISCRETE , 1999, Evolution; international journal of organic evolution.

[13]  M. Mauldon,et al.  Fracture sampling on a cylinder: From scanlines to boreholes and tunnels , 1997 .

[14]  D. Karlen,et al.  Aerial dry matter and nutrient accumulation comparisons among five soybean experiments , 1995 .

[15]  Brian Berkowitz,et al.  Theory of anomalous chemical transport in random fracture networks , 1998 .

[16]  Estimating the theoretical semivariogram from finite numbers of measurements , 2000 .

[17]  A. Warrick,et al.  Optimization of Sampling Locations for Variogram Calculations , 1987 .

[18]  N. Cressie,et al.  Robust estimation of the variogram: I , 1980 .

[19]  A Y Yakovlev,et al.  A distribution of tumor size at detection and its limiting form. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[20]  P. Bickel,et al.  NONPARAMETRIC INFERENCE UNDER BIASED SAMPLING FROM A FINITE POPULATION , 1992 .

[21]  David Russo,et al.  Design of an Optimal Sampling Network for Estimating the Variogram , 1984 .

[22]  E. Ariga.  Availability and role of multipurpose trees and shrubs in sustainable agriculture in Kenya , 1997 .

[23]  C R Weinberg,et al.  The design and analysis of case-control studies with biased sampling. , 1990, Biometrics.

[24]  R. Webster,et al.  Sample adequately to estimate variograms of soil properties , 1992 .

[25]  Y. Rubin.  Subsurface Flow and Transport: Transport of inert solutes by groundwater: recent developments and current issues , 1997 .