MATLAB 平台下基于克立格理论的固体资源储量估算 系统的开发与实践

文章以MATLAB技术为平台,以西藏自治区雄村铜矿床为实例,详细讨论了如何开发以克立格理论为基础的固体资源储量估算系统,其中包括模块:数据预处理、变异函数分析、克立格估值、储量估算、地质可靠度评价等.论文使用该系统对西藏自治区雄村铜矿床进行储量预测,得出结果数据表格;并对原生矿体、次生矿体、氧化矿体、总矿体、矿体的自然形态以及在工业品位规范下的矿体形态进行拟合.通过与专业储量估算软件Minesjght的估算结果对比,证明了本系统的可靠性和实用性.     

基于地统计学的江西省年降水量插值研究

本文选取江西省86个气象站点2008年的年降水量数据,基于地统计学理论对数据进行探索性空间数据分析,得出该组气象数据插值的理论最优模型。采用不同全局趋势及变差函数模型进行插值,通过交叉性检验证实了该理论模型的有效性。最终结果显示,采用3次函数剔除全局趋势并以球状模型作为变差函数模型的考虑高程的Co-Kriging方法插值结果最优。     

丘陵地貌区高程内插适宜方法选择——以长春净月潭为例

对不同的地貌区进行空间内插时选择的适宜方法不同,以典型丘陵地貌区长春市净月潭地区的高程为研究对象,用ArcGIS 9.2分别作反距离权重法、最近邻点法、趋势面分析和克里金法做高程内插。应用交叉检验方法对其内插的精度作出分析,检验其科学性和可靠性,得出反距离权重法、最近邻点法、趋势面法、样条函数法和克里金法高程内插的不同精度,并对各种方法的适用性进行讨论。结果表明最近邻点法精度最高,以下依次为普通克里金插值法、样条函数、反距离权重法、趋势面法。为类似于净月潭地区的丘陵地貌区高程内插提供选择参考。     

基于谱烈度改进算法的研究与应用

本文基于Housner谱烈度公式中存在的不足和问题提出了一种改进算法,首先将相对速度反应谱积分公式上下限范围扩大为0.1～10.1 s,通过统计分析集集地震、汶川地震和芦山地震反应谱积分值和宏观烈度之间的关系,给出8组分档线性公式来计算改进谱烈度值I_s,并采用普通克里金插值法绘制谱烈度分布图。本文提出的改进算法在地震动加速度时程积分与谱烈度之间建立起一组长周期、连续性和精细化的函数关系,经鲁甸地震对比验证,该改进算法得出的谱烈度分布图与宏观烈度图整体变化趋势一致,谱烈度值与宏观烈度吻合率(偏差±1度以内)约为90%,具有一定的科学性和准确性,可为灾后应急救援以及决策部署等提供科学的数据支撑。     

贵州百花湖表层水中不同形态汞的分布规律

均匀布点采集百花湖表层水样44个,对水样不同形态汞含量进行了分析,并现场测定水质参数。采用经典统计学与半变异函数拟合相结合的方法,研究了百花湖表层水中不同形态汞含量的空间变异性。根据半变异函数拟合的结果,在ArcGIS软件的支持下进行克立格插值,以此来揭示百花湖表层水中不同形态汞含量的空间分布规律。研究表明:百花湖表层水中的汞污染主要来自上游,随水体流动汞污染减少;但也有一些区域可能因富营养化而使颗粒态汞和总汞含量增加。     

中国大陆速度场模型建立方法研究

将中国大陆划分成东北、华北、华南、青藏、西北5个块体,基于克里金插值法、最小二乘配置法、最小二乘配置移动拟合推估法等3种方法及速度场数据研究中国大陆速度场模型建立方法。利用外检点分析不同建模方法的精度情况,最小二乘配置法在大部分块体（华北、华南、东北块体）的拟合精度都较高,表现稳定|克里金插值法在青藏、华北块体的拟合精度较高。利用克里金插值法与最小二乘配置法获取中国大陆1°×1°格网中心点速度场成果,两种方法成果差异较小。     

云南地区近期地壳活动特征

基于云南省的GNSS连续观测资料,通过最小二乘、克里金插值等方法进行速度场、应变参数的求解,获得了自2010年起至今云南地区的地壳活动特征。结果表明,云南地区最大剪应变变化速率存在西强东弱,北强南弱的特征,目前滇西北地区最大剪应变增强最快,滇南至滇西南等区域挤压应变积累最明显,且2011年缅甸7.2级地震的发生未改变这一地区的应变积累趋势。     

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

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

基于地统计方法的气候要素空间插值研究

在回顾了地统计学的产生、发展及其基本原理的基础上,对于目前众多可以提供计算网格的气候要素的空间插值方法中,具体探讨了普通克里格法和协同克里格法。将甘肃省1961—1990年30年平均降水量和蒸发量作为区域化变量,根据不同的半变异函数理论模型,采用普通克里格法和双变量协同克里格法,通过对比分析得到：（1）不论是多年均降水量还是多年平均蒸发量在空间上都呈现明显的梯度变化,二者的空间变程都很大,而降水量变化幅度更大。降水量从东南部向西北部逐渐减少,蒸发量则相反,从东南向西北逐渐增加。（2）基于地统计的插值方法,根据半变异函数云图和试验方差最小的原理,选择合适的半变异函数理论模型进行变量的空间插值,能够较好地模拟区域化变量的空间连续分布格局,并取得较好的效果。对比普通克里格法和协同克里格法,后者增加了高度对降水量和蒸发量的影响,在空间分布上更为合理,插值的精度也要明显好于普通克里格法。（3）采用地统计方法虽然在总体上能够较好地反映气候要素的空间分布格局,但检验显示,两种方法空间插值的精度都还不是很高,插值的精度还有待进一步提高。     

基于块克里金的土壤水分点观测向像元尺度的尺度上推研究

土壤水分是陆面水文过程的一个重要分量,土壤水分的地面点观测与卫星观测的尺度是不匹配的,用点观测数据进行遥感反演结果的验证或者融合这两种观测都需要开展点观测数据向卫星像元尺度的尺度上推研究.土壤水分时空异质性的分析和描述是进行尺度上推的基础,地统计方法是描述连续随机变量空间结构的经典方法.应用块克里金法分别将黑河寒区遥感试验阿柔试验区2008年4月1日与L波段微波辐射计同步的地面液态含水量和含冰量点观测数据经过尺度转换,得到与遥感像元相匹配的像元平均估计值和标准差,可用于该天L波段微波辐射计土壤液态含水量和含冰量反演结果的真实性检验,估计结果充分利用了像元临近位置的观测,得到了比直接的采样平均更合理的块估计结果,对块克里金估计值与采样平均值进行比较,发现两种结果趋势是一致的,块克里金法提供了更为合理的块估计结果,土壤水分空间结构的时间变异、小尺度的土壤特性的变化和测量误差都会对估计结果带来一定的不确定性.