Core-log integration: Optimal geostatistical signal reconstruction from secondary information

Abstract

This paper deals with the problem of reconstructing a significantly under-sampled signal of interest (primary signal) using the few data available for that signal and an exhaustively sampled secondary signal. This is a typical situation in various geophysical settings, as for example in core-log integration. Core measurements are made in the laboratory and are very precise but, because core recovery may be fragmentary along the borehole, there are large gaps in the information these measurements provide. On the other hand, wireline logging offers a complete set of measurements along the borehole. The variable of interest (e.g. porosity) can be measured directly in the cores (i.e. primary variable) while the logging variables are proxy variables (e.g. neutron porosity, electrical resistivity, sonic velocity). The optimal combination of both types of measurements can improve the assessment of the variable of interest. Geostatistics offers a methodology for merging such different types of information. Assuming a linear relation between the primary and secondary variables, various geostatistical techniques are applicable: simple kriging with an estimated local mean, kriging with an external drift, cokriging and Bayesian integration. These methods are compared in a simulation experiment in which the signal is considered to be discrete with 1024 data values and reconstruction is assessed for the case when 10% of the data are known (i.e. 102 data values) and for the case when 50% of the signal is known (512 data values). The various methods are assessed using statistics of the real error involved in the reconstruction and the ratio between the real error and the estimated error. The methods are ranked according to the values of these statistics and their inference and computational requirements. The methodology is demonstrated on a case study in which a few tens of core samples of effective porosity are integrated with a much more exhaustive set of wireline neutron porosity.

Keywords

GEOSTATISTICS GEOPHYSICS DATA FUSION OPTIMAL INTERPOLATION SIGNAL RECONSTRUCTION SIMPLE KRIGING KRIGING WITH AN EXTERNAL DRIFT COKRIGING BAYESIAN INTEGRATION

Get full access to this article

View all access options for this article.

References

Anonymous: ‘The geology and hydrogeology of Sellafield', Nirex report 263, UK Nirex Ltd, 1992.

Anonymous: ‘The geology of the region around Sellafield', Nirex report 264, UK Nirex Ltd, Oxfordshire, 1992.

Anonymous: ‘Structural geology of the Sellafield region', Nirex report 265, UK Nirex Ltd, Oxfordshire, 1992.

Anonymous: ‘Sellafield: Site-district Geology', Nirex Report 266. UK Nirex Ltd, Oxfordshire, 1992.

Anonymous: ‘The geological structure of the Sellafield site', Report S/97/007, UK Nirex Ltd, 1997.

Brereton

N. R.

, Rogers

S. F.

and Evans

C. J.

: ‘Spatial heterogeneity of rock mass properties', BGS Technical report WK/95/10C, Geol. Cat. No. 01/1464, UK Nirex Ltd, Oxfordshire, 1996.

Chen

, Hubbard

and Rubin

: ‘Estimating the hydraulic conductivity at the South Oyster site from geophysical tomo-graphic data using Bayesian techniques based on the normal linear regression model', Water Resourc. Res., 2001, 37, (6), 1603–1613.

Deutsch

C. V.

and Journel

A. G.

: ‘Geostatistical software library and user's guide'; 1992, New York, Oxford University Press.

Dowd

P. A.

: ‘Geostatistical characterization of three-dimensional spatial heterogeneity of rock properties at Sellafield', Trans. Inst. Min. Metall A, 1997, 10A, A133—A147.

10.

Domenico

P. A.

and Schwartz

F. W.

: ‘Physical and chemical hydrogeology'; 1990, New York, John Wiley and Sons.

11.

Doyen

P. M.

: ‘Porosity from seismic data: a geostatistical approach', Geophyiscs, 1988, 53, (10), 1263–1275.

12.

Draper

and Smith

: ‘Applied regression analysis', 2nd edn; 1981, New York, John Wiley & Sons.

13.

Harvey

P. K.

and Lovell

M. A.

(eds.): ‘Core-log integration', Special Publication 136; 1998, London, Geological Society.

14.

Journel

A. G.

and Huijbregts

: ‘Mining geostatistics', 600; 1978, New York, Academic Press.

15.

Kingdon

, Rogers

S. F.

, Evans

C. J.

and Breton

N. R.

: ‘The comparison of core and geophysical log measurements obtained in the Nirex investigation of the Sellafield region', in ‘Core-log integration', (ed. Harvey

P. K.

and Lovell

M. A.

), Special Publication 136, 97–113; 1998, London, Geological Society.

16.

Lee

P. M.

: ‘Bayesian statistics: an introduction', 294; 1989, New York, Oxford University Press.

17.

Myers

R. H.

: ‘Classical and modern regression analysis', 2nd edn, 488; 1990, Belmont, CA, Duxbury Press.

18.

Proceedings of the Yorkshire Geological Society, Special issue on ‘Advances in West Cumbrian geology resulting from the Nirex investigations around Sellafield', Vol. 50, Part 1; 1994, London, Geol. Soc. Publishing House.

19.

Robinson

G. K.

: ‘A role for variograms', Aust. J. Statist., 1990, 32, (3), 327–335.

20.

Wackernagel

: ‘Multivariate geostatistics', 256; 1995, Berlin, Springer-Verlag.

21.

, Tran

T. T.

, Srivastava

R. M.

and Journel

A. G.

: ‘Integrating seismic data in reservoir modelling: the collocated cokriging alternative', Proc. 67th Ann. Tech. Conf. Exhb., Washington, DC, USA, October 1992, SPE24742, 833–842.