小波去噪的广西加权平均温度插值研究

针对广西地区探空站稀少,难以获得精确的T_m问题,GGOS atmosphere提供了利用ECMWF的相关资料计算而得到的时间分辨率为6h(UTC 00:00:00,06:00:00,12:00:00,18:00:00)、空间分辨率为2.5°×2°的全球T_m格网数据可以在没有气象数据的情况下获得较高时空分辨率的T_m,该文利用GGOS atmosphere T_m格网数据对广西地区4个探空站插值T_m,并用无线电探空数据计算的T_m检验其精度;对误差进行分析后,选取最优小波基与尺度对其残差去噪,利用去噪后得到的曲线建立T_m的改正模型。实验结果表明,插值T_m经基于小波去噪的模型改正后,其RMSE为1.29K;Bevis模型的RMSE为10.71K;GPT2_1W模型的RMSE为3.56K;改正模型精度优于传统模型,可以达到地基GPS反演GNSS-PWV的精度要求。     

复杂地形地区月平均气温(混合)地理加权回归克里格插值

以地形地貌特征复杂、观测站点分布稀疏不均匀的四川省为研究区,引入地形因子（坡度和坡向）和植被指数,采用顾及空间关系非平稳性的（混合）地理加权回归克里格模型（（mixed）geographically weighted regression Kriging,（m）GWRK）进行月尺度平均气温插值方法及精度分析研究。针对不同季节和不同地区,将（m）GWRK插值结果与基于全局回归的回归克里格（regression Kriging,RK）插值结果进行对比。结果表明,RK、GWRK、mGWRK回归关系的决定系数R2分别为0.795、0.922、0.911,均方根误差分别为0.83℃、0.64℃、0.55℃,表明GWRK、mGWRK对目标变量的解释能力以及插值精度都优于RK;GWRK、mGWRK相对于RK对月平均气温插值的改进具有季节与地区差异,冬半年的改进大于夏半年,在地形地貌变化大的地区改进大于地形地貌变化小的地区。     

干热河谷冲沟DEM插值误差的空间分布研究

为了分析云南元谋干热河谷典型冲沟插值误差的空间分布特征,采用反距离加权（inverse distance weighting,IDW）、局部多项式（local polynomial interpolation,LPI）、张力样条（spline with tension,ST）、析取克里格（disjunctive Kriging,DK）以及不规则三角网（triangulated irregular network,TIN）模型方法对高程采样点进行插值,用交叉验证法、相对差系数及沟谷线差异衡量其插值精度。遴选高程误差大于1 m的误差点,用变异系数（coefficient of variation,CV）、全局Moran指数和Getis-Ord G_i*指数分析其空间格局特征。结果表明:TIN和DK精度较高,IDW精度最低;高程误差均呈聚集分布,聚集程度TIN > LPI > DK > ST > IDW;高程误差均呈空间正自相关,TIN模型插值误差的自相关程度最高;误差热点位于坡度大的区域。     

基于厄米特插值和时间序列模型在沉降监测中的应用

针对沉降监测中用时间序列分析数据样本较少、处理精度不高的问题,本文引入厄米特插值,提出用厄米特插值法对数据进行插值,增加数据样本,然后用时间序列对数据进行建模分析。为了验证用厄米特插值后的数据建模精度,将上述理论运用到工程实例中,用时间序列分析分别对厄米特插值处理后的数据和原始数据进行建模分析,得到用厄米特插值处理后的数据进行时间序列分析建模的平均绝对误差比原始数据建模的平均绝对误差高0.083 8 mm。说明厄米特插值和时间序列组合模型精度比较高。     

渤海湾盆地南堡凹陷南部深层沙一段优质碎屑岩储层特征及成因机制

利用岩芯、薄片、扫描电镜和分析化验等资料,对南堡凹陷南部深层古近系沙一段（沙河街组一段）优质碎屑岩储层特征进行了分析并揭示其成因机制。研究表明：①南堡凹陷南部沙一段优质碎屑岩储层属于低孔中渗储层,发育原生孔隙、次生孔隙和微裂缝;②沙一段沉积期处于强水动力沉积环境,辫状河三角洲前缘水下分流河道砂体普遍发育,岩性主要为中粗砂岩、含砾砂岩和砂砾岩等粗岩性,岩石类型主要为岩屑长石砂岩,石英、岩浆岩岩屑等刚性组分含量高,储层在强压实作用下仍保留了局部原生残余粒间孔,部分刚性组分强压条件下形成微裂缝;③沙一段与东三段为区域不整合接触,大气水的酸性流体在成岩初期以及烃源岩成熟之后产生的有机酸在成岩后期通过油源断裂和不整合面进入储层,造成长石、易溶碎屑等组分溶蚀形成次生孔隙,溶蚀作用不仅增加孔隙空间,而且扩大喉道宽度。结论认为：粗岩性和高刚性颗粒含量为原生孔隙的保存提供了物质基础;油源断层、不整合面为酸性流体的进入提供路径,促进了次生孔隙发育。     

基于颗粒识别分析系统的碎屑流堆积物颗粒识别和统计方法研究

碎屑流堆积物颗粒识别和统计是碎屑流灾害的研究重点。文章基于图像处理孔隙（颗粒）及裂隙图像识别与分析系统（PCAS）,以贵州纳雍普洒村崩塌-碎屑流为例,结合纳雍崩塌堆积物粒径实测结果,通过阐释识别过程中阈值、孔喉封闭半径、最小孔隙面积的参数意义,研究PCAS软件在碎屑流颗粒识别与统计中的应用,并提出了颗粒识别时这些参数的选取方法。分析结果表明:（1）PCAS能自动准确地识别碎屑流堆积物颗粒与孔隙,相比人工计数更精细,所识别堆积物各区小颗粒比重较大,0～2 m颗粒粒径各区占比都在50%以上;（2）当阈值为170（像素）时能获得精细的二值图像,颗粒与孔隙得到了准确地区分;（3）不同参数取值下获得堆积物颗粒粒径分布结果不同,碎屑流堆积物颗粒识别宜采用较大的孔喉封闭半径和较小孔隙面积,当二者比值为3/30（像素）时能更好地反应颗粒粒径分布情况;（4）PCAS具有较高的可行性,统计结果显示,各粒径含量变化趋势与人工统计相近,两种统计方法各粒径占比、分布规律基本吻合,说明利用PCAS可以实现对崩塌碎屑流颗粒粒径分布的高效便捷分析。     

基于非均质地层参数表征的颗粒流模型构建

非均质地层参数的获取与表征方法是浅层地质体改造的重要依据。基于非均质地层的地质建模,表征具有随机不确定性和模糊不确定性的土体表观参数和物理力学参数。利用 Weibull 分布统计描述局部地层信息,并应用Diamond-Square分形插值方法,结合差分盒子维计算非均质地层分形维数,进一步合理演算整体地层参数。采用PFC3D数值模拟软件,结合马尔科夫链蒙特卡洛方法采样数据,对非均质地层块石颗粒及土颗粒进行建模。通过非均质地层模拟结果,发现其破坏过程与应力应变曲线要比均质土体复杂得多,由于块石大颗粒的含量、分布特性对土样力学性质影响显著。本文建模方法有助于非均质地层参数表征研究。     

非饱和土中水平入渗方程显式求解

在求解非饱和态土中水分入渗问题时,水力函数是体积含水率或者吸力的函数,致使其控制方程呈现出强非线性的特征,进而使得其求解变得十分困难。基于水分在土体介质中流动耗时取极值路径的选择这一假定,引入时间泛函,基于变分法原理将水平入渗问题转化为泛函极值问题。通过求解Euler–Lagrange方程,结合边界条件,得到非线性瞬态水平入渗问题的显式解析解。结合Brooks-Corey型水力函数,显式地求解出该类型非饱和态土的体积含水率发展分布规律。通过计算4种不同类型土体的水平入渗规律,将求解结果与已有结果以及数值结果进行对比,验证了该方法的有效性。结果表明:体积含水率分布与位置距离和湿润峰距离比值呈幂函数关系,指数取决于土-水特征曲线的形状参数;初始条件与边界条件会对体积含水率分布造成不同程度的影响。     

内蒙古东北部地区一次极端降雪过程的水汽输送特征

利用内蒙古呼伦贝尔市常规观测资料和GDAS、NCEP/NCAR再分析资料,采用欧拉方法分析了2016年春季内蒙古东北部地区一次极端暴雪过程的水汽输送及收支特征,利用HYSPLIT模式和聚类分析模拟计算了此次暴雪天气过程的水汽源地、主要水汽输送通道及其对水汽输送的贡献,并与传统的欧拉方法结果进行对比。结果表明:(1)有3支不同源地的水汽流在内蒙古东北部地区交汇,对呼伦贝尔地区暴雪的发生与维持有重要影响;(2)经向和纬向输送为此次暴雪天气的发生提供了充足的水汽,暴雪区水汽主要源于中高层的南边界和随西风气流的西边界;(3)利用HYSPLIT模式模拟发现,在此次暴雪天气过程中水汽主要来源于新地岛以西洋面、日本海以及巴尔喀什湖,且三者贡献率大致相当。     

平缓地区地形湿度指数的计算方法

地形湿度指数( topographic wetness index) 可定量模拟流域内土壤水分的干湿状况, 在流域 的土壤及分布式水文模型等研究中具有重要的意义。但现有的地形湿度指数计算方法在应用于 地形平缓地区时会得到明显不合理的结果, 即在河谷地区内, 地形湿度指数仅在狭窄的汇水线上 数值较高, 而在汇水线以外的位置则阶跃式地变为异常低的地形湿度指数值。本文针对此问题对 地形湿度指数的计算方法提出改进: 以多流向算法MFD- fg 计算汇水面积, 相应地以最大下坡计 算地形湿度指数, 再基于一个正态分布函数对河谷平原地区内的地形湿度指数进行插值处理。应 用结果表明, 所得地形湿度指数的空间分布不但能合理地反映平缓地区坡面上的水分分布状况, 并且在河谷地区内地形湿度指数值也都比较高, 其空间分布呈平滑过渡, 因而整个研究区域的水 分分布状况得到了比较合理的反映。