地统计方法学研究进展

doi:10.11821/yj2008050022

摘要/Abstract

摘要：

地统计方法学已经成为空间预测与不确定性分析的关键性工具。本研究从文献计量学和方法学演变过程两个角度开展了1967~2005年的地统计方法学研究综述研究。首先从宏观角度分析其发展趋势、应用情况与模式,然后总结其演化规律、适宜性和选择原则。研究表明:地统计方法学的演化规律表现为稳态向非稳态演变、单变量向多变量(含二次信息)演变、参数与非参数方法相互补充、线性向非线性方法演变和空间静态向时空动态演变;此外,未来研究发展方向集中在半变异函数估计新方法、不确定性地统计学、时空地统计学与多点地统计学、机理模型与地统计学耦合研究和基于地统计学模拟的不确定性决策等。

关键词: 地统计学, 文献计量学, 演变过程, 半变异函数, 不确定性

Abstract:

Geostatistical methodology becomes an important tool for spatial prediction and uncertainty analysis. In terms of scientometric analysis and methodology development, advanced research for geostatistical methodology from 1967 to 2005 has been proposed in this study. First, its developing trend, application and patterns were identified; its development mode, suitability and choosing principles were then summarized. Development mode could be summarized as five transformations from stationary to non-stationary, single-variable to multivariate, parametric to non-parametric, linear to non-linear and spatial static to spatiotemporal dynamic geostatistics. Finally, the future researches were preliminarily discussed in this field, which is of great significance to geostatistical methodology and applications in the future. It includes: 1) developing new methods for variogram estimation to reduce the analytical complexity, 2) refining uncertain geostatistics to reflect hybrid uncertainties (input/output dataset, model structure, model parameters), 3) exploiting spatiotemporal/multi-point geostatistics and its algorithm/software to simulate complex realizations, 4) developing hybrid approach of process-mechanism and geostatistics to comprehensively uncover potential process-mechanism, in spite of the difficulties within the process simulations, and 5) expanding geostatistical-based decision-making under uncertainty, which supports the risk-based decision-making in social work, environmental pollution, agricultural production, public health, etc.

Key words: geostatisitcs, scientometrics, development mode, variogram, uncertainty

郭怀成, 周丰, 刀谞. 地统计方法学研究进展[J]. 地理研究, 2008, 27(5): 1191-1202.

GUO Huai-cheng, ZHOU Feng, DAO Xu. State-of-art on geostatistical methodology[J]. GEOGRAPHICAL RESEARCH, 2008, 27(5): 1191-1202.

参考文献

[1] 侯景儒, 黄竞先. 地质统计学及其在矿产储量计算中的应用.北京: 地质出版社,1982.

[2] Krige D G. A statistical analysis of some of the borehole values in the Orange Free State gold field. Journal of the Chemical and Metallurgical Society of South Africa,1952,53:47~64.

[3] Matheron G. Principles of geostatistics. Economic Geology, 1963,58:1250~1266.

[4] Journel A G. Geostatistics for conditional simulation of ore bodies. Economic Geology,1974,69(5):673~687.

[5] Johnson G A, Mortensen D A, Gotway C A. Spatial and temporal analysis of weed seedling populations using geostatistics. Weed Science,1996,44(3):704~710.

[6] Schlesinger W H, Raikes J A, Hartley A E, et al. On the spatial pattern of soil nutrients in desert ecosystems. Ecology,1996, 77 (2): 364~374.

[7] Goovaerts P. Geostatistics in soil science: state-of-the-art and perspectives.Geoderma,1999,89(1-2):1~45.

[8] Goovaerts P. Geostatistical approaches for incorporating elevation into the spatial interpolation of rainfall.Journal of Hydrology,2000,228(1-2):113~129.

[9] 王政权.地统计学及在生态学中的应用. 北京:科学技术出版社,1999.

[10] Kunzli N, Jerrett M, Mack W J, et al.Ambient air pollution and atherosclerosis in Los Angeles.Environmental Health Perspectives, 2005,113(2):201~206.

[11] Isaaks E H, Srivastava R M. An Introduction to Applied Geostatistics. New York: Oxford University Press,1989.

[12] Cressie N A C. Statistics for Spatial Data. New York: Wiley-Interscience,1993.

[13] Goovaerts P. Geostatistical modelling of uncertainty in soil science. Geoderma,2001,103(1-2):3~26.

[14] 王仁铎, 胡光道. 线性地质统计学. 北京: 地质出版社,1988.

[15] 孙洪泉. 地质统计学及其应用.徐州:中国矿业大学出版社,1990.

[16] Tao S. Factor score mapping of soil trace element contents for the Shenzhen area. Water Air and Soil Pollution,1998,02 (3-4): 415~425.

[17] 柏延臣,李新,冯学智.空间数据分析与空间模型.地理研究,1999,18(2):185~190.

[18] Li B G, Cao J, Liu W X, et al. Geostatistical analysis and kriging of Hexachlorocyclohexane residues in topsoil from Tianjin, China. Environmental Pollution,2006, 142 (3): 567~575.

[19] Gustafson E J.Quantifying landscape spatial pattern: What is the state of the art? Ecosystems,1998,1(2):143~156.

[20] McLaughlin D, Townley L R. A reassessment of the groundwater inverse problem. Water Resources Research,1996,32 (5):1131~1161.

[21] McBratney A B, Santos M L M, Minasny B. On digital soil mapping.Geoderma,2003,117 (1-2): 3~52.

[22] Atkinson P M, Lewis P. Geostatistical classification for remote sensing: an introduction. Computers & Geosciences,2000,26(4):361~371.

[23] 刘瑞民, 王学军,陶澍,等. 天津地区表土PAHs含量的多尺度空间相关性.地理研究,2005, 24(4):565~571.

[24] Wang Y, Ma T, Luo Z. Geostatistical and geochemical analysis of surface water leakage into groundwater on a regional scale: A case study in the Liulin karst system, northwestern China. Journal of Hydrology,2001,246 (1-4): 223~234.

[25] Sun B, Zhou S L, Zhao Q G. Evaluation of spatial and temporal changes of soil quality based on geostatistical analysis in the hill region of subtropical China.Geoderma,2003,115(1-2): 85~99.

[26] Zhou F, Guo H C, Ho Y S, et al.Scientometric analysis of geostatistics using multivariate methods. Scientometrics,2007,73(3): 265~279.

[27] Goovaerts P. Geostatistical analysis of disease data: Visualization and propagation of spatial uncertainty in cancer mortality risk using Poisson kriging and p-field simulation. International Journal of Health Geographics,2006, 5: 7~32.

[28] Krieger N, Chen J T, Waterman P D, et al. Geocoding and monitoring of US socioeconomic inequalities in mortality and cancer incidence: Does the choice of area-based measure and geographic level matter? The Public Health Disparities Geocoding Project. American Journal of Epidemiology,2002,156(5):471~482.

[29] Krige D G.Some basic considerations in the application of geostatistics to the valuation of ore in South African gold mines: South African Institute of Mining and Metallurgy,1976,76(9):383~391.

[30] Hoeksema R J, Kitanidis P K.An application of the geostatistical approach to the inverse problem in two-dimensional groundwater modeling.Water Resources Research,1984,20(7):1003~1020.

[31] Hoeksema R J, Kitanidis, P K.Comparison of Gaussian conditional mean and Kriging estimation in the geostatistical approach to the inverse problem. Water Resources Research, 1985, 21(6), 825~836.

[32] Journel A G. Geostatistics: Models and tools for the earth sciences. Mathematical Geology, 1986, 18: 119~141.

[33] Rouhani S, Wackemagel H. Multivariate geostatistical approach to space-time data analysis. Water Resources, 1990,26:585~591.

[34] Dowd P A.A review of recent developments in geostatistics. Computer & Geosciences, 1992, 17(10): 1481~1500.

[35] Christakos G, Li X Y.Bayesian maximum entropy analysis and mapping: A farewell to kriging estimators? Mathmetical Geology, 1998, 30:435~462.

[36] Wingle W L, Poeter E P, Mckenna S A.UNCERT: Geostatistics, uncertainty analysis and visualization software applied to flow groundwater and contaminant transport modeling. Computers and Geosciences,1999,25:365~376.

[37] Lark RM. Designing sampling grids from imprecise information on soil variability: An approach based on the fuzzy kriging variance. Geoderma, 2000,98 (1-2):35~59.

[38] Saito H, Goovaerts P. Accounting for source location and transport direction into geostatistical prediction of contaminants. Environmental Science and Technology, 2001,35(24): 4823~4829.

[39] Pardo-Iguzquiza E, Dowd P A. Empirical maximum likelihood kriging: The general case. Mathematical Geology,2005,37(5):477~492.

[40] Abbaspour K C, vanGenuchten M T, Schulin R, et al. A sequential uncertainty domain inverse procedure for estimating subsurface flow and transport parameters.Water Resources Research, 1997,33(8):1879~1892.

[41] Russo D, Russo I, Laufer A. On the spatial variability of parameters of the unsaturated hydraulic conductivity. Water Resources Research,1997,33(5):947~956.

[42] Yeh T C J, Jin M H, Hanna S. An iterative stochastic inverse method: Conditional effective transmissivity and hydraulic head fields. Water Resources Research, 1996,32 (1): 85~92.

[43] Yu W H, Harvey C M, Harvey C F. Arsenic in groundwater in Bangladesh: A geostatistical and epidemiological framework for evaluating health effects and potential remedies. Water Resources Research, 2003,39 (6): 1146~1162.

[44] de Marsily G, Delay F, Goncalves J, et al. Dealing with spatial heterogeneity. Hydrogeology Journal,2005,13 (1): 161~183.

[45] Journel A, Zhang T F. The necessity of a multiple-point prior model. Mathematical Geology,2006,38 (5): 591~610.

[46] Kyriakidis P C, Journel A G. Stochastic modeling of atmospheric pollution: A spatial time-series framework. Part I: Methodology. Atmospheric Environment,2001,35 (13): 2331~2337.

[47] Liu S Y, Yeh T C J. An integrative approach for monitoring water movement in the vadose zone. Vadose Zone Journal,2004,3(2):681~692.

[48] Brus D J, Heuvelink G B M. Optimization of sample patterns for universal kriging of environmental variables. Geoderma,2007,138 (1-2): 86~95.

[49] Pardo-Iguzquiza E, Dowd P A. The second-order stationary universal kriging model revisited. Mathematical Geology,1998,30 (4): 347~378.

[50] Wackernagel H. Multivariate Geostatistics: An Introduction with Applications.Berlin:Springer-Verlag,2003.

[51] Juang K W, Lee DY. Simple indicator kriging for estimating the probability of incorrectly delineating hazardous areas in a contaminated site. Environmental Science & Technology,1998,32 (17): 2487~2493.

[52] Saito H, Goovaerts P. Accounting for measurement error in uncertainty modeling and decision-making using indicator kriging and p-field simulation: Application to a dioxin contaminated site. Environmetrics,2002,13(5-6): 555~567.

[53] Gomez-Hernanez J J, Sahuquillo A, Capilla J E. Stochastic simulation of transmissivity fields conditional to both transmissivity and piezometric data:I. Theory. Journal of Hydrology,1997,203 (1-4): 162~174.

[54] Armstrong M. Basic Linear Geostatistics. Berlin: Springer-Verlag, 1998.

[55] Deutsch C V. A sequential indicator simulation program for categorical variables with point and block data: BlockSIS. Computers & Geosciences,2006,32 (10):1669~1681.

[56] Emery X. Properties and limitations of sequential indicator simulation. Stochastic Environmental Research and Risk Assessment,2004,18 (6): 414~424.

[57] 孙英君,王劲峰,柏延臣.地统计学方法进展研究.地球科学进展,2004,19(2):268~274.

[58] 侯景儒,王志民,潘汉军,等. 时间空间域中多元信息的地质统计学. 北京科技大学学报,1995 ,17(2) :101~106.

[59] Christakos G, Serre M L. Spatiotemporal analysis of environmental exposure-health effect associations.Journal of Exposure Analysis and Environmental Epidemiology,2000,10 (2): 168~187.

[60] Akita Y, Carter G, Serre M L. Spatiotemporal nonattainment assessment of surface water tetrachloroethylene in New Jersey. Journal of Environmental Quality, 2007,36 (2):508~520.

[61] Douaik A, Van Meirvenne M, Toth T. Soil salinity mapping using spatio-temporal kriging and Bayesian maximum entropy with interval soft data.Geoderma,2005,128 (3-4): 234~248.

[62] Olea R A. A six-step practical approach to semivariogram modeling.Stochastic Environmental Research and Risk Assessment,2006,20(5):307~318

[63] Lark R M.Estimating variograms of soil properties by the method-of-moments and maximum likelihood. European Journal of Soil Science,2000,51 (4): 717~728 DEC 2000.

[64] Minasny B, McBratney A B. The Matern function as a general model for soil variograms.Geoderma,2005,128(3-4): 192~207.

[65] Todini E, Ferraresi M.Influence of parameter estimation uncertainty in Kriging.Journal of Hydrology,1996,175(1-4):555~566.

[66] Pardo-Iguzquiza E, Dowd P.Variance-covariance matrix of the experimental variogram: Assessing variogram uncertainty.Mathematical Geology. 2001,33(4): 397~419.

[67] Ortiz C J, Deutsch C V.Calculation of uncertainty in the variogram.Mathematical Geology,2002,34(2):169~183.

[68] Marchant B P, Lark R M. Estimating variogram uncertainty.Mathematical Geology,2004,36(8):867~898.