某区地表水质量空间信息统计分析研究.docVIP

某市市区地表水质量空间信息统计分析研究 高加庆1， 2 （1. 苏州科技学院土木工程学院，苏州 215011；2. 苏州科技学院 土木工程学院，苏州 215011） 摘 要： 介绍了空间信息统计分析原理及计算方法；运用空间信息统计的原理，以变差函数为工具，对苏州市地表水中总氮（NH4+-N）含量迸行了分折，得到两个不同方向上的实验变差函数及理论变差函数模型，并用克立格估值法，对其含量的空间分布进行了模拟，直观显示了其变化的空间分布特征，为总氮（NH4+-N）的含量提供直观的空间变化形态。 关键词： 空间信息统计；变差函数；克立格估值 The Geostatistics Analytical Investigation： 分析研究 of Surface Water Quality in City GAO Jia-qing1, SUN Hong-quan2 (1.School of Civil Engineering , Suzhou University of Science and Technology , Suzhou 215011; China 2. School of Civil Engineering , Suzhou University of Science and Technology , Suzhou 215011 , China) Abstract： The principles and methods of geostatistcs are introduced. Using the principles of geostatistcs and with the help of variogram, the NH4+-N content in rivers is analysed. The models of variogram in the different directions are obtained. By using the Kringing, the special distribution of the NH4+-N content is simulated. The three dimensional shape of the water content is shown intuitively. Keywords: geostatistics; variogram; kringing 空间信息统计是年创立并发展起来的一门新兴边缘科学。1.1 变差函数 为一组观测值，其中为第i个观测点上的坐标，i个观测点上的观测值。称下列公式 （1） 为实验变差函数。 其中：——为实验变差函数值；——观测点的坐标，在二维空间中它代表；——为i点上的观测值；——为距i点的距离为h的那一点的观测值；h ——为两观测点之间距离，称为滞后距；——为相距为h的点对的数目。 最小的滞后距称为基本滞后距。在实际计算中不同的滞后距一般是成倍增加的。对不同的滞后距h，公式（3）可以算出相应的实验变差函数值。对每一滞后距hi，把诸点[hi, ]在h—图上标出，所提到的图形称为实验变差函数图，或实验变差图。实验变差函数具有方向性，按照某一方向取滞后距h，可得该方向上的实验变差函数。 1.2 理论变差函数[hi, ]拟合成理论变差函数模型。常用的理论变差函数模型是球状模型。球状模型变差函数的数学表达式为 （2） 其中：a为变程；C0为块金值；C为拱高；C0+C为基台值。球状模型变差函数的图形可以表示为 从球状模型理论变差函数图上我们可以看出：变差函数在（变程）范围内，随的增大而增加，a变程反映了变量相关性的范围；当时，变差函数就不再单调增加了，而是稳定在C0+C附近C0+C反映了变量变化性的大小；块金值反映了变量连续性优劣。 1.3 克立格估值与克立格方差 该估值方法是南非采矿工程师克立格（D.G.Krige）首次提出来的，。用已知点上的数值来预测（估计）未知点上的值称为估值。通常采用的线性估计式为 （） 为未知点上的估计值；为权系数的确定。克立格估值的权系数是在“无偏”和“方差最小”的前提下导出的，所以它是一种最优线性无偏估计量（the Best Linear Unbiased Estimator）。 ① 无偏条件为为该点上的真实值，估计总会有偏差的，但我们总可以使这偏差值的数学期望满足无偏性条件： （4） ② 方差最小