中国耕地土壤相对湿度时空分异

以全国653个农业气象站1993-2013年耕地的土壤相对湿度数据为基础,运用地统计方法,分析中国耕地土壤相对湿度时空分异特征与规律。结果表明：自1993年以来全国耕地的土壤相对湿度呈现波动上升趋势。全国耕地的土壤相对湿度普遍大于60%,分布区域自4月中旬开始随夏季风推移不断向北向西扩大,自10月下旬开始向东、南方向缩小。耕地土壤相对湿度值随土壤深度的增加而增加。年际尺度上,耕地的土壤相对湿度在夏秋季上升速度最快,变化幅度随土层深度增加而变小。土壤相对湿度与降水量存在较强的正相关关系,与潜在蒸发量、气温普遍存在较强的负相关关系。土壤相对湿度与各气象要素的相关性随深度加深而减弱。春、夏、秋季气象因素对旱地土壤相对湿度影响较大,冬季气象因素对水田影响更大。     

基于分区间泛化Kriging近似模型的船舶局部结构形状优化

为提高Kriging近似模型在船舶结构性能多维度响应预测方面的适用性,对常规Kriging近似模型进行分区间泛化改进:一是对设计样本点的各个维度(分量)进行划区,并在每一个分区间内采用最优拉丁超立方(OLhd)取样;二是引入比例系数w1组合高斯与指数型相关函数提高模型对数据的泛化能力,从而在每个划分的区间内建立泛化的Kriging近似模型。通过SCH测试函数,验证了构造的分区间泛化Kriging模型有效性。结合参数化建模和改进Kriging近似模型对某过渡肘板和舱口角隅边界进行形状优化,优化第一阶段由改进的近似模型通过多岛遗传算法得到全局初步的最优解,第二阶段在初步解的基础上缩小优化变量范围,由少量的FEM计算即可在小范围内搜寻到精确的最优设计变量。结果表明:分区间泛化Kriging近似模型在预测多维度响应时较常规Kriging模型预测精度更高;分阶段的形状优化流程在保证极小误差和缩小计算成本的情况下能够得到理想的应力分布和重量优化结果,有助于船舶结构的轻量化研究。     

面向稀疏散布数据集的时空Kriging优化

时空Kriging法通过将变异函数向时空域进行扩展得到时空变异函数,有效地利用时空邻近的采样点综合进行插值,由于时空稀疏散布数据集具有单一时刻下样本点数量少以及时空分布不规律的特点,难以满足使用时空Kriging插值法的基本条件,导致插值精度不高,据此本文提出了优化方法：通过多时段叠置拟合空间变异函数的方法,综合利用时空邻域内的采样点以解决单一时刻下空间邻域内数量不足情况;控制时间变异对空间变异函数拟合的误差影响;采用积合式模型构建时空变异函数进行插值。最后使用Argo海温数据进行插值实验,在相同条件下与时空Kriging法以及时空权重法的交叉验证结果对比得出,该方法在保证拟合所需采样点数量要求的同时,有效削减了一般时空Kriging法中时间变异对空间变异函数拟合结果的干扰,插值结果的绝对误差均值从0.5降低至0.2以内,稳定性进一步增强,改善了时空Kriging法在稀疏散布数据条件下精度上的不足。     

序贯高斯模拟在矿石品位估计中的应用研究

随机模拟是地质统计方法的重要内容。在矿石品位估计方法中克里格方法作为一种无偏估计方法,常被用于矿石品位的估计。但克里格法估值存在平滑效应。作者在分析序贯高斯模拟和普通克里格法基本原理的基础上,运用序贯高斯模拟方法和普通克里格方法对某铁矿体内全铁(TFe)品位进行估计,给出了品位估计结果模型。研究从勘探线方向、垂直勘探线方向和竖直方向分别计算变差函数,对球状模型、指数模型、高斯模型的变差函数拟合效果进行了优选,结果表明球型模型拟合效果最好。针对序贯模拟和克里格品位估值效果进行了分析,结果显示:序贯高斯模拟结果在品位分布形态上更接近样品品位分布形态,其平滑效应更小;克里格方法估计与序贯高斯模拟方法相比仅在品位均值方面更接近样本品位均值。因此,认为序贯高斯模拟方法可以更好地刻画矿体内品位分布状态。     

云南个旧多金属矿石门坎—保和地区Sn、Sb元素的成矿信息识别与提取

石门坎—保和地区位于云南个旧多金属成矿区西南部,成矿背景复杂,成矿作用具有多期性和多样性。利用研究区1∶1万土壤化探数据,采用泛克里格法和因子分析相结合的方法,提取不同层次的成矿异常信息。结果表明:主因子F_3包含Mo-Au-Sb等元素,反映成矿过程中的岩体和热液相关信息;主因子F_4包含Sn-Bi-W元素,指示中低温热液矿化环境; Sb沿杨家田断裂呈串珠状分布,反映了由北东向构造控制的空间变化特征,与已知矿点对应良好; Sn元素异常在研究区北部呈双峰状环形分布,推测受隐伏岩体或者断层控制。根据元素组合异常F_3和F_4的分布规律,结合地质资料圈定出2个找矿靶区,为进一步找矿工作提供了方向。     

自适应顺序采样Kriging方法及其在增量动力分析中的应用

复杂结构的增量动力分析（IDA）对于结构的抗震设计和分析有着重要意义,但需对结构进行大量的非线性时程分析,计算量成本高。本文结合Kriging元模型和自适应顺序采样并用于结构增量动力分析以提高其计算效率和精度,其中:Kriging元模型用于预测结构的地震响应,顺序采样根据候选点的熵值补充非线性时程分析逐步增加Kriging模型的预测精度。借助本文方法,IDA曲线可通过少量的时程分析实现较高的精度。为了校验本文方法的可行性与有效性,对二层和九层钢框架结构模型应用直接IDA、hunt&fill方法和本文方法分析并比较了上述三种方法的计算误差、计算效率和IDA曲线差别。在此基础上本文将自适应顺序采样Kriging方法用于考虑结构不确定性参数的IDA分析,并和传统的蒙特卡洛方法进行比较。结果表明:该方法具有较高的计算效率,可以保证IDA曲线的精度。     

福建省晋江市城市热岛强度时空变化特征分析

利用传统的气象站法, 结合空间统计学方法(普通克里金插值法), 对福建省晋江市2010—2014年40个自动气象站逐小时温度资料加以计算处理, 分析了晋江市年、季、昼夜热岛强度时空变化规律。(1)晋江市年、季、昼夜热岛强度都呈带状分布, 等值线呈西南-东北走向, 年、季、昼夜变化趋势显著, 北部热岛强度高于南部。五年间热岛强度持续增强, 但增幅不大, 增速放缓。(2)城市化水平的提高, 会导致热岛强度高值出现季节提前, 故旅游区秋冬季热岛强度高于春夏季, 中心城区和产业经济区夏秋季热岛强度高于冬春季。(3)晋江市热岛效应昼夜空间分布格局差异性大, 夜间热岛强度显著高于白天, 最低值出现在14—16时, 中心城区和产业经济区最低值出现时间较旅游区略推迟, 三个功能区的最高值均出现在凌晨。      

基于MODIS雪盖数据的北疆雪深多元非线性回归克里金插值

为了提高北疆地区雪深时空分布监测的准确性,以该区域48个气象站点2006年12月—2007年1月的月平均雪深观测数据为基础,通过分析月均雪深空间自相关性及其与经纬度、高程的相关性,结合MODIS雪盖数据构建了多元非线性回归克里金插值方法,插值获得了北疆地区较高精度的雪深空间分布数据。将插值雪深数据与普通克里金插值法、考虑高程为辅助变量的协同克里金插值法的预测结果进行比较,结果表明:1相对普通克里金和协同克里金方法,多元非线性回归克里金法的12月份雪深预测精度分别提高了15.14%和9.54%,1月份的提高了4.8%和6.7%;2由于充分利用了经纬度和地形信息,多元非线性回归克里金法的雪深预测结果可提供更多细节信息;3预测结果客观地表达了雪深随经纬度和地形变化的趋势,反映了积雪深度的空间变异性;4基于不显著相关的协变量高程的协同克里金插值法预测的雪深数据精度劣于普通克里金插值法的预测结果。     

三维激光扫描技术在复杂地质隧道岩体信息获取中的优化

为了更高效、高精度地获取岩体结构面信息,借助三维激光扫描技术,提出并建立一种基于Kriging算子和最小二乘法组合的点云插值方法。首先对获取的隧道点云数据利用单一的最小二乘法进行插值,并以此为基础,结合Kriging算子的适应度函数建立组合插值模型。在利用三维激光扫描仪已获取的某岩溶隧道的原始点云数据的基础上,运用组合模型进行插值处理,以此获取岩体结构面信息,并将组合方法得到的结果与单一的插值方法进行对比,发现组合方法插值效果较好,误差也较小。基于Kriging算法和最小二乘法结合的插值方法,在利用三维激光扫描技术获取岩溶隧道岩体几何信息—结构面产状的过程中能更加精确、全面,并可以用于岩体稳定性评价。     

基于Kriging方法的负偏态分布数据的空间插值研究

Kriging插值方法的前提是数据符合或基本符合正态分布的假设,或通过对非正态数据进行正态变换,使其满足要求。但对于拖尾的负偏态分布数据,正态变换方法无法奏效,若直接对数据进行Kriging插值,精度较低。文中在对样本数据深入分析的基础上,借助GIS平台,从不同透视角度对数据进行趋势分析,利用二次多项式模型对数据进行全局趋势拟合,并对随机的短程变异进行Kriging插值,选取不同的变异函数模型进行拟合比较,最终确定最优的模型参数,达到较为理想的插值结果。