中国西部主要台站积雪深度的空间插值研究

应用空间插值方法,对中国西部地区(27.17～48.05°N,79.05～103.57°E)110个气象台站的观测数据进行空间内插,研究积雪分布特征。结果表明,普通Kriging法和Cokriging法都能够反映出研究区积雪深度分布的空间结构特征,与实际情况比较吻合。但相比之下,Cokriging法的精度更高,局部特征的反映更佳。分析发现,影响插值结果精度的主要因素是研究区内气象台站稀少,且空间分布很不均匀。通过合理的采样设计,选择合适的插值方法,并考虑地形、气候等影响积雪分布的因素将有利于改善空间插值精度。     

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

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

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

对不同的地貌区进行空间内插时选择的适宜方法不同,以典型丘陵地貌区长春市净月潭地区的高程为研究对象,用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%,得到了比普通克里金和线性插值方法更接近真实情况的地层厚度估计值。     

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

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

基于三维地质模型的排土场边坡整体稳定性探究

目前所广泛采用的平面应变分析手段,针对排土场边坡稳定性评价时有时会出现某剖面安全系数高于设定允许值,实际中却遭破坏的现象,除了常被提及的力学参数、本构模型等因素,排土场凹形夹持效应和凸形发散作用突出以及滑体端部效应的影响,是分析结果不能准确反映实际的根本原因之一。鉴于此,从平面高程数据、地质剖面、钻孔资料等条件,配合插值效果较为稳定的Kriging法对空间数据进行插值,建立排土场三维地质整体模型。然后,采用自编FISH语言接口程序,将地质模型导入FLAC3D中离散化并计算,较为简单、有效地实现三维地质建模、模型可视化、岩土工程计算分析三位一体功能,有效地克服了平面分析无法解决的上述问题。以攀钢兰尖铁矿尖山第7排土场为实,对其整体稳定性进行应力、应变分析,采用三维强度折减法计算排土场三维整体安全系数,与现场调查及平面计算结果较为一致。基于该模型的计算结果可以有效避免排土场空间效应和滑体端部效应,得到的排土场全域应力-应变规律,对露天矿排土场的生产管理,具有一定的指导意义。