克立格方法在区域化探数据处理中的应用

针对克立格方法能表征变量空间结构特征 ,以团结沟金矿研究区的化探数据处理为例 ,讨论了一般多元统计方法、协同克立格方法和普通克立格方法在化探元素间的相关分析结果。协同克立格方法对原始数据进行估值得到的绝对误差的均值要比用其它方法得到的绝对误差的均值低 ,标准差也明显降低。相比之下 ,协同克立格方法在一定程度上较为优越 ,它提高了地质统计学的预测精度 ,缺点是计算工作量较大。     

时空域中的指示克立格理论研究

对时空域内的指示克立格进行了较全面的研究,主要研究内容有：（１） 时空域内的指示克立格（ＳＴＩＫ） 变换;（２） 指示二阶距;（３）单元指示克立格（ＳＴＩＫ）：①简单指示克立格（ＳＴＳＩＫ） ,②普通指示克立格（ＳＴＯＩＫ） ,③中值指示克立格（ＳＴＭＩＫ）;（４）多元指示克立格：①指示协同克立格（ＳＴＣＯＩＫ）,②主分量指示克立格（ＳＴＰＣＩＫ）     

协同克立格法在广西林旺金矿化探数据处理中的应用

协同克立格法,同时兼具化探数据多元性及克立格法表征空间属性的特点。考虑到地质变量两个以上的空间属性,在化探数据处理中运用协同克立格法,可以进一步提高估计精度。运用协同克立格法对广西林旺矿区中金元素成矿进行预测,经相关性分析显示,在矿区中元素的基本组合是Au、Ag、As、Hg,协同克立格插值以Au为主要变量,伴生元素Ag、As、Hg为次要变量。把插值计算后的结果与传统多元统计方法和普通克立格法的计算结果进行比较,结果协同克立格法得到的估值误差较小,预测精度较高,在成矿元素的预测中具有一定程度的优越性。     

纯时间域多元信息地质统计学

随着人们动态地认识客观世界的要求日益增加,时空域中的多元信息地质统计学技术研究受到重视。本文在时间域内对几种纯时间域克立格技术进行了研究,重点对纯时间域中的普通克立格（ＴＯＫ）、简单克立格（ＴＳＫ）及泛克立格（ＴＵＫ）进行了论述,也对纯时间域中的协同克立格（ＴＣＯＫ）和指示克立格（ＴＩＫ）作了简要的概括。     

时空多元协同克立格的理论研究

本文结合地质统计学的最新成果,在空间协同区域化理论的基础上,对时空域中多元信息的协同克立格（ＳＴＣＯＫ）理论进行了较为详细的研究。主要研究内容有：（１）ＳＴＣＯＫ中的互变异函数与互协方差函数;（２）ＳＴＣＯＫ方程组及求解估值权因子的三种方法：①传统普通协同克立格法（ＳＴＴＯＣＯＫ）,②标准普通协同克立格法（ＳＴＳＯＣＯＫ）,③简单协同克立格法（ＳＴＳＣＯＫ）;（３）排列协同克立格（ＳＴＣＯＬＣＯＫ）;（４）指示协同克立格（ＳＴＩＫＣＯＫ）。     

农田土壤水分和盐分空间分布的指示克立格分析评价

运用非参数地质统计学的多元指示克立格法,结合单元指示克立格法,对黄河河套灌区55hm2的盐渍化土地上两个比较关键时期的土壤水分和盐分进行空间分布的分析,给出了同一时期土壤水盐和不同时期土壤盐分满足一定标准的综合概率图,并借助该案例详细介绍了多元指示克立格的原理、步骤及应用范围。指示克立格方法为不同尺度(从农田到区域)上水土资源质量的评价提供了新的方法,能为水土资源的现代决策管理提供指导。     

月尺度的协克立格模型在水文数据插补延长中的应用

月尺度的协同克立格模型能够对水文变量进行线性、无偏和最佳估计[2]。针对水文数据在时间上富足而空间上缺乏以及单变量具有相依性、多变量之间又有相关性的特点,以滦河流域内蒙古段的大河口、多伦和花塘沟三个雨量站的实测月降水量资料为基础,用月尺度的协克立格模型对缺测数据进行了插补延长。经验证明估值精度较高,结果可信。同时还对影响估值精度的数据结构进行了讨论。结果表明,数据结构对估值精度的影响主要取决于观测数据间存在的相关性。     

析取克立格方法分析

非线性空间信息统计学主要包括对数正态克立格、指示克立格、析取克立格、条件期望等理论方法。通过源代码和运行结果对各种克立格方法在数据特点、时间复杂度等方面进行了比较。DK估计较普通克立格更有效。首先，DK考虑到了Z(x)的二维概率分布，充分应用有效信息，克服了线性克立格的前提条件太少，导致数值离散性大时，估计不够准确的缺点；其次，DK利用区域化变量的非线性组合估计，所以可以方便地得到变量的任意函数的估计结果，在局部可回采储量等方面比线性方法更容易、结果更准确。     

协同泛克立格法—数据处理的新方法

协同克立格法是多元地质统计学的一个重要方法，它以协同区域化变量为研究对象，充分考虑了变量的空间相关性和变量间的统计相关性，因而有着单变量克立格法无可比拟的优点。非平稳条件下的协同泛克立格法不仅能对变量进行线性、无偏、最优估计，而且能对漂移和剩余进行估计。协同泛克立格法在数据处理中有着十分乐观的发展前景。     

归来庄金矿床化探信息金异常的协同克立格法研究

根据收集到的归来庄金矿床的化探资料,选取与Au相关系数较大的W、As、Sb等区域化探元素作为变量,对该区金异常空间分布进行协同克立格法研究.在对4个区域化变量进行空间结构分析的基础上,建立空间结构模型并进行估计.综合化探金异常的协同克立格法研究结果表明,该区有多处金异常存在,具有较好的找矿前景.     

储层物性的空间结构特征分析与预测—以青海尕斯库勒油田E3^1油藏为例

储层物性受诸多地质因素的影响，经常表现出强烈的非均质性。这种非均质性在一定的空间尺度上往往具有明显的结构特点。本文以青海尕斯库勒油田E3^1油藏为例，讨论了地质统计学方法在油田储层特性空间结构特征分析与预测中的应用。采用变异函数定量描述了孔隙度和渗透率的空间分布结构特征，并在此基础上利用克立格法进行了最成插值预测。结果表明，研究区储层物性（孔隙度和渗透率）具有显著的空间结构性特点，变程一般在800－2000m之间；不同小层的储层物性具有不同的空间结构方向性。这种特征主要沉积相带空间展布的影响，各小层孔隙度和渗透率的实验半变异函数均可用具有块金效应的球状模型来拟合并进行预测。作为验证，本文还采用“多重趋势面”预测模型对储层的孔隙度和渗透率进行了预测分析。     

储层物性的空间结构特征分析与预测——以青海尕斯库勒油田E_3~1油藏为例

储层物性受诸多地质因素的影响,经常表现出强烈的非均质性。这种非均质性在一定的空间尺度上往往具有明显的结构特点。本文以青海尕斯库勒油田E_3~1油藏为例,讨论了地质统计学方法在油田储层物性空间结构特征分析与预测中的应用。采用变异函数定量描述了孔隙度和渗透率的空间分布结构特征,并在此基础上利用克立格法进行了最优插值预测。结果表明,研究区储层物性(孔隙度和渗透率)具有显著的空间结构性特点,变程一般在800～2000m之间;不同小层的储层物性具有不同的空间结构方向性。这种特征主要受沉积相带空间展布的影响,各小层孔隙度和渗透率的实验半变异函数均可用具有块金效应的球状模型来拟合并进行预测。作为验证,本文还采用“多重趋势面”预测模型对储层的孔隙度和渗透率进行了预测分析。     

基于指示和普通克里金的不连续地层厚度估计方法

因为克里金插值方法考虑了样本的空间统计特征而被广泛应用于地层建模。但对于地层尖灭、缺失等不连续地层情况,常用的克里金方法（如普通克里金方法）通常给出较为平滑的地层厚度估计结果,在地层分布边缘的局部范围内与实际情况相差较大。针对不连续地层,提出一种采用指示克里金和普通克里金相结合的地层厚度估计方法。该方法首先采用指示克里金方法估计地层分布范围,然后采用普通克里金方法估计分布范围内的地层厚度,并根据估计标准差来评价地层厚度的不确定性。将提出的方法应用于上海长江隧道的地层建模,通过交叉验证,证明了该方法在地层分布边缘的局部范围内能明显改善地厚度估计结果,标准差降低了15%～18%,得到了比普通克里金和线性插值方法更接近真实情况的地层厚度估计值。     

克里金参数估值法及其在参数估计分析中的应用

为考虑岩土介质参数的空间分布的结构性和随机性等不确定因素,引入了克里金参数估值法。采用变异函数描述参数在空间结构上的变化,建立其空间变异规律的数学模型,从而实现对岩土参数的估值。通过实例分析揭示了克里金估值法具有反映“过滤效应”和“集团效应”的优点,对不同位置的数据赋予不同的权重系数,能够有效地反映参数空间变异结构,有利于对参数的合理化分析     

克立格估计邻域大小的确定方法

分析了以往地质统计学确定估计领域域时的不足,提出了克立格估计方差、均值权,真值与估计值之间线性表达式的斜率、真值与估计值的相关系数以及克立格估值中的负权样品数等可作为确定克立格估计领域域大小的关键参数,对于脉状脉体而言,除了考虑以上参数外,还应限制矿体厚度方向上的大小。     

结合软数据的同位置协同克里金估值新方法

协同克里金法结合软数据可以提高地质体建模中估值的准确性。通过分析传统的结合软数据同位置协同克里金估值过程中存在的问题及原因,提出并推导了一种新的结合软数据的同位置协同克里金估值方法。该方法采用软数据的期望值以保证外推的可靠性,利用主变量观测值滤去协同变量的期望值以估计变异函数参数。该方法的特点是:当待估点位于主变量相关变程之内时,估值的结果主要反映主变量变化特征;当待估点落在相关变程之外时,估值结果的差别主要来源于各自点上协同变量(软数据)的期望值。最后,用实例验证了新估值方法的优点。      

储集层非均质性分析的计算机方法:以西弗吉尼亚州下密西西比系油田为例

本文使用几种计算法和统计法研究美国西弗吉尼亚洲中部格兰尼克里克油田大Ｉｎｊｕｎ砂岩储集层的结构及其与采油量的关系。计算机程序根据已出版的资料编写，以适合（ｌ）计算渗透率对孔隙度的回归；（２）标绘三维孔隙度；（３）确定和标绘根据地球物理测井记录推断的相；（４）估算原始储量、累积产量和采油率。因为回归分析显示出测井记录和岩心孔隙度及岩心渗透率间微弱但重要的相关，所以，孔隙度可用干地层渗透率中以构成储集层的模式。也使用定量地层对比和多维定算法来估算在缺乏可用数据的情况下构造的影响。于井间使用克里格法内插绘制的剖面，突出了孔隙度较高的地带。使用地球物理测井资料和岩心描述，用聚类分析确定导电相。原始储量的估算结果与累计产量数据相结合，得出采油率的估算值。这些采油量变量图通常呈现在岩相图上观察到的相同的南北走向。     

多孔介质非均质性空间变异结构分析

非均质多孔介质的空间分布具有随机性与结构性的双重性质.地统计学方法在研究具有结构性和随机性双重性变量的空间变异性方面具有极强的科学性.本文在地理信息系统平台ArcGIS上,采用Geostatistical Analyst地统计学分析模块中的kriging插值方法研究非均质性多孔介质的空间变异性.以河北省廊坊市地下水污染物的现场数据为研究对象,利用普通克立格法对氟化物、硫酸盐、氯化物和汞浓度的空间变异性进行分析.     

环渤海低平原水土盐分与水位埋深的空间变异及协同克立格估值

针对环渤海低平原水土盐分和水位埋深,通过选取0～20cm深度内的127个代表性土样进行土壤全盐量测定,并对130个水井的水位埋深及128个地下水样的矿化度进行测定.综合运用普通克立格(OrdinaryKriging)与协同克立格(CoKriging)方法与GIS技术研究水土盐分及水位埋深的空间分布状况.结果表明,土壤盐...     

克立格法在地下水数值模型中的应用

回顾了克立格法在水文地质学中的应用现状及进展，介绍了克立格法的基本原理，论述了克立格法在地下水模拟数值模型中应用的思路与实现方法，举例说明了利用克里格法插值处理含水层结构与初始水头分布。