考虑权值非负约束的泛克立格算法

泛克立格法是地质统计学的一种重要估值方法,当区域化变量在空间变异几何域内非平稳时常用泛克立格法来估值。然而,用泛克立格方程组计算估计权值时由于未对权的符号作任何限制,从而使得计算出的权值经常出现负权现象,而负权的存在有许多弊端,所以在许多应用中有必要对权值作非负要求。本文基于线性规划方法提出了一种考虑权值非负约束的泛克立格算法,该算法既考虑到了权值的非负约束条件,又利用了线性规划方法求解简便快捷的优点。     

一种考虑权值非负约束的克立格算法

基于线性规划方法提出了一种考虑权值非负约束的克立格算法,该算法既利用了线性规划求解简便、快捷的优点,又克服了其它正克立格法计算工作量大以及受主观因素影响的缺点。     

次生晕数据的对数正态泛克立格法研究及异常评价

首先介绍了对数正态泛克立格法的基本理论和方法,包括对数正态分布,三参数对数正态分布,求解估值(Z_v)~*所需的权系数λ_α.及求解漂移值(m_v)~*所需的权系数ρ_α的对数正态泛克立格方程组.其次,笔者应用对数正态泛克立格法探讨了华北某测区的次生晕数据、元素统计分布特征以及变异函数及结构分析.最后,给出测区估计结果,元素综合异常图,划分出5个异常区,并结合地质条件对异常进行了综合评价.     

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

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

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

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

基于下限原理有限元的强度折减法

对于岩土工程中常用的强度折减系数,其规划问题是非线性的,不能直接利用线性规划进行求解。基于四边形单元的下限原理有限元法,根据强度折减系数与超载系数近似符合双曲函数的特点,通过调整强度参数使得超载系数逼近于1[1],可将边坡稳定性分析中常用的强度折减系数的非线性规划求解转化为线性规划求解问题。分析表明,采用拟合双曲线插值法求解强度折减系数的计算效率高于常规的二分法及割线法,且具有较好的收敛性;该方法能够充分利用当前高效的线性规划算法,便于工程应用。     

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

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

化探制图中校准分析偏倚的克里格法

高水平上的分析偏倚,是地球化学制图中遇到过的严重问题。以往认为这种偏倚的校准是非常困难的。在实践基础上,本文提出应用“克立格法”——“地质统计学”的原理,可以有效的校准这种偏倚。本方法实施的要点是模拟实验半变异函数模型,建立克立格方程组求解权系数,并利用斯米尔诺夫法求出校正系数,根据校正数据进行地球化学制图。     

限制性克立格法在阿舍勒铜矿储量计算中的研究与应用

为提高矿产储量计算的精确性，本文运用了限制性克立格法。给出了限制性克立格法的计算理论，并设计出限制性克立格法在阿舍勒铜矿中的储量计算流程，实现了对阿舍勒铜矿体的限制性克立格储量计算。结果表明，限制性克立格法的储量计算结果（783654．39t）要比普通克立格法的储量计算结果（706541．53t）更接近于实际勘探值（919454t），限制性克立格储量计算方法能够取得较为理想的预测结果。     

克立格法在确定地下水位监测网点最优增补位置中的应用

本文运用克立格法（Kriging)对新监测点增补前后的已有监测点水位估值的克立格方差进行了计算,在此基础上利用精度选择法对监测点最优增补位置进行了选择。     

Kriging without negative weights

Under a constant drift, the linear kriging estimator is considered as a weighted average ofn available sample values. Kriging weights are determined such that the estimator is unbiased and optimal. To meet these requirements, negative kriging weights are sometimes found. Use of negative weights can produce negative block grades, which makes no practical sense. In some applications, all kriging weights may be required to be nonnegative. In this paper, a derivation of a set of nonlinear equations with the nonnegative constraint is presented. A numerical algorithm also is developed for the solution of the new set of kriging equations.     

理论变异函数球状模型的加权线性规划法似合

对理论变异函数球状模型及其套合结构拟合这一问题作了探讨，提出了加权线性规划拟合法。同于该法在目标函数中既可考虑到对不同滞后h下所得实验变异函数值进行加权，又可保证拟合成功，同时还可进行人工干预，因此，方法综合了现有加权多项式拟合法及线性规划拟合法的各自优点，且在计算上较目标规划拟合法更为简单。     

克立格算法的若干思考

克立格法是一种最优、线性(或非线性)、无偏内插估计量的方法.随着信息社会的发展,克立格法的应用范围越来越广.文中介绍了克立格法相关的计算方法,并基于Windows编写克立格算法流程图,提出改进克立格算法的若干方法,包括:①强化数据预处理;②将基本滞后距和实验变差函数计算结合起来一并运算;③防止数据溢出;④选择快速有效算法求解克立格方程组.     

A note on programs performing kriging with nonnegative weights

This note deals with the problem of solving kriging systems with nonnegative weights. Mathematically the question is considered as a case of quadratic programming. Hints concerning computational approaches and software are given.     

基于Windows克立格算法的可视化编程与应用

随着信息社会的发展,克立格法已被广泛应用于矿产储量计算,物探、化探、遥感地质数据处理,环境科学、水文地质、农林科学和农田水利等诸多方面.文章综合分析克立格各种计算方法,对克立格复杂计算过程进行模块分解,基于Windows编写了克立格算法流程图并在Visual C 6.0平台上实现可视化编程,并对二维普通克立格算法的优化进行了有益的探索.     

用协同—指示克立格预测储层物性参数：以胜利油田罗家砂四上砾岩为例

指示克立格是一种非参数统计方法,它能在不舍弃特异值的条件下进行有效的空间估计,协同克立格以协同区域化变量的研究对象,充分考虑了变量空间相关性和变量间的统计相关性,有着单变量克立格法可可比拟的优点,因而,本文提出协同－指示克立格法这一方法,在指示克立格中考虑多变量的空间相关性和变量间的统计相关性,使指示克拉立格与协同克立格相互补充,模拟现实,提高预测精度,并以胜利油田罗家砂四上砾岩为例进行参数预测,     

地质三维属性建模及其可视化

从运畴学的角度出发,提出了基于带权线性规划的方法拟合变差函数,有效地保证了求解结果的非负性。将样本点建立空间三维网格索引,大大简化了待插点邻域样本的搜索过程,提高了Kriging属性建模的速度。最后结合OpenGL图形库,实现了属性模型的三维可视化,从而有效地实现了地质信息的定量化。     

变差函数的参数模拟

本文给出利用线性方程组非负解理论，进行地质统计学变差函数理论模型参数最优拟合，并给出了计算实例，从而为实现地质统计学计算过程的自动化提供了重要方法。