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

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

导纳函数的中国南海海底地形模型

针对目前大面积海域依然存在海深测量数据空白的问题,该文选取中国南海4°×4°(12~16°N,115~119°E)海域范围为研究区域,通过频率域上的相关性分析,得到重力异常在30~130km波长范围与研究海域海底地形的相关程度较高。以ETOPO1海深模型作为先验模型,应用导纳函数实现了在无船测数据情况下对海深的反演,最终得到1′×1′的海深模型。将反演结果与检核点进行比较发现,模型1在2 500m以上海域相对误差较小,反演精度较高;2 500m以下海域相对误差变化大,反演精度较低。无船测数据环境下,适当加入一定数量船测海深值作为控制点得到的海深模型2,相较于未加控制点的模型1,在1 000m以下水深处的标准差明显优于模型1,与检核点的差值精度最大提高了45%左右。     

附不等式约束平差模型的一种快速搜索算法

大地测量中常存在一些先验不等式约束信息,充分利用它们可以保证参数解的唯一性和稳定性。然而,现有的不等式约束平差算法主要是基于优化理论,算法通常比较复杂,需要选取有效约束或建立罚函数。在最小二乘平差准则基础上,把不等式约束看成是一个可行域,借助Fisher函数在可行域中快速搜索使误差平方和达到最小的最优解,推导出了可行解为最优解的充分必要条件。建立了基于Wolfe-Powell算法的非精确快速搜索算法,从而减小了搜索算法的计算量,得到了一种新的不等式约束平差计算方法。该算法的平差准则与最小二乘平差准则一致,不需要矩阵求逆运算,可适用于维数较大的平差问题解算。     

利用模拟退火算法反演多波束测量声速剖面

针对多波束测量声速剖面站点布设密度不够而引起的声速剖面代表性误差问题,提出一种利用模拟退火算法反演声速剖面的方法。首先,对测区已有声速剖面序列进行经验正交函数（empirical orthogonal function,EOF）分析,利用声速扰动矩阵和前几阶EOF求得EOF重构系数及其变化范围。其次,采用模拟退火算法对EOF重构系数进行迭代优化,以多波束测得的海底地形畸变量大小为依据构建目标函数,并设定合理的退火控制参数。最后,得到待反演区域的声速剖面数据。实例分析表明,该方法反演的声速剖面较时间就近原则选取的替代声速剖面更接近真实声速剖面,且利用反演声速剖面改正后的海底地形更接近真实地形,有效削弱了声速剖面代表性误差的影响,显著提高了多波束测量精度及数据处理效率。     

基于样条函数的海岸带变化监测影像局部配准优势分析

海岸带变化监测以遥感影像为矢量更新依据。在实际应用中,由于数据源复杂,全色影像与底图套合、多光谱影像与全色配准容易出现局部超限,为项目生产带来一定的难度。本文主要介绍利用样条函数进行影像的局部配准以满足项目精度要求。     

确定砂土动剪切模量和阻尼比的方法对比

确定土体动剪切模量的常用方法有规范法、Kumar法和自相关函数法,确定相应阻尼比的方法有规范法、Das and Luo法、Kokusho法、Kumar法和互相关函数法,为了分析不同方法所产生差异,实现定量化对比分析,笔者以福建标准砂（粒径为0.5~1.0 mm）为研究对象,采用不排水的应力控制动三轴试验,探讨不同的确定土体动剪切模量和阻尼比方法的差异性,并给出了不同土体条件建议选用的方法。结果表明：1）3种方法确定动剪切模量的结果有一定的差异,随剪应变的增大结果的差异逐渐增大,有效围压对结果的差异性有所影响,当剪应变为4×10^-3,有效围压为100 kPa时,3种方法差异显著,相对误差最大接近20%;2）而5种方法确定阻尼比的结果差异显著,随着剪应变的增大,5种方法确定的阻尼比相对误差大体上均在迅速减小,只有规范法在有效围压为100 kPa时,其相对误差有较小的增大趋势;5种方法中,Kumar法确定的阻尼比最接近平均阻尼比,互相关函数法远高于平均阻尼比,Das and Luo法和Kokusho法确定的阻尼比基本一致但低于平均阻尼比。建议以后的工程应用中,加载方式为应力控制时,可采用自相关函数法确定动剪切模量,采用Kumar法确定阻尼比,二者确定的动剪切模量和阻尼比均最接近平均值。     

寒区渠基黏土热参数最优概率分布

为分析寒区渠基黏土热参数的随机分布特征及概率分布模型,以寒区渠基黏土的导热系数为样本,结合经典分布拟合法、多项式逼近法、最大熵法和正态信息扩散法,分别对寒区渠基黏土热参数的概率分布规律进行了研究。首先通过分析热参数的离散性,并比较概率分布曲线、拟合检验值和累计概率分布值,对不同方法描述热参数随机性的优劣进行了评价;然后,基于寒区渠基黏土热学参数对温度的敏感性,提出了一个可以达到理想拟合精度的寒区渠基黏土热参数概率推断的区间取值标准。研究结果表明：寒区渠基黏土的热参数具有随机变量的特征;正态信息扩散法可以描述热参数样本的随机波动性;在4种方法中,正态信息扩散法的拟合精度最高。使用3.5σ法,将[μ-3.5σ,μ+3.5σ]（μ为随机变量的均值,σ为标准差）作为概率函数推断时的取值区间,同时考虑偏度的影响,可使得累计概率值达到1.000 0的精度,能够较准确地推断热参数的概率分布函数。     

基于省区协调发展的GDP预测模型

基于柯布-道格拉斯生产函数与自回归移动平均模型(ARIMA)构建出一个GDP综合预测模型,并且考虑十九大全面建成小康社会与实现共同富裕的精神与国家关于技术、资本、劳动力等方面的区域平衡发展战略调整模型的参数,计算了2016—2050年中国分省的GDP总量与人均GDP,进一步通过计算省区间人均GDP的基尼系数来分析省区协调发展的水平。研究结果表明,在考虑省区协调发展时,各省区在2016—2050年间的GDP总量与人均GDP的差距逐渐缩小,省区间人均GDP的基尼系数将从2015年的0.219下降到2030年的0.176和2050年的0.137,未来区域间发展不均衡的态势在实现经济稳步增长同时可以得到缓解。     

温室气体排放的环境库兹涅茨曲线检验与减排路径的冲击动态——基于黑龙江省统计核算数据

基于2001-2015年黑龙江省温室气体排放统计核算数据,对地区GDP与温室气体排放的环境库兹涅茨曲线关系检验呈现倒U型,预期2019年达到理论拐点;通过偏最小二乘回归模型得到4个减排路径的年平均减排效果顺序依次为单位GDP化石能源消费量减少、经济结构调整、人均GDP增长、贸易结构变化;减排路径对应脉冲响应函数的动态冲击效果分别为波动性增排、收敛性减排、发散性减排、转变的排放作用;推动黑龙江省温室气体减排的路径顺序为控制化石能源消费量、优化经济结构、发展低碳经济、调整贸易结构。