GPS精密单点定位中精密星历和钟差的内插算法选取

利用了IGS中心给出的精密星历,计算了几种插值以及拟合算法,分析和比较了插值与拟合的质量与运算效率,基于精密单点定位软件验证与对比了其精度。     

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

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

电导率各向异性的海洋电磁三维有限单元法正演

本文提出了一种基于非结构化网格的海洋电磁有限单元正演算法.为了回避场源奇异性, 文中选用二次场算法, 将背景电阻率设置为水平层状且各向异性, 场源在水平层状各向异性介质中所激发的一次场通过汉克尔积分得到.基于Coulomb规范得到二次矢量位和标量位所满足的Maxwell方程组, 通过Galerkin加权余量法形成大型稀疏有限元方程, 采用不完全LU分解(ILU)预条件因子的quasi-minimum residual(QMR)迭代解法对有限元方程进行求解得到二次矢量位和标量位; 进而, 利用滑动平均方法得到二次矢量位和标量位在空间的导数, 由此得到二次电磁场; 通过一维模型对算法的可靠性进行验证, 与此同时, 针对实际复杂海洋电磁模型, 比较有限元模拟结果与积分方程模拟结果, 进一步验证算法精度.若干计算结果均表明, 文中算法具有良好的通用性, 适用于井中电磁、航空电磁, 环境地球物理等非均匀且各向异性介质中的电磁感应基础研究.     

基于Kriging算法的三棱柱地质体三维可视化实现

采用了一种基于Kriging算法的三棱柱三维可视化方法对地质体进行了可视化,其基本思想是首先用Kriging算法将空间任意的离散的数据点生成规则的点数据,将点数据划分成多个三棱柱,然后以三棱柱为基本单元进行地质体可视化。本系统是在Visual C++/OpenGL平台上开发而成,具有精度高、交互好的特点,经过实践运行,切实可行。     

气候模式模拟中海量数据的并行抽取与可视化分析诊断

在地球动力学和气候模拟等领域, 数值模拟产生的数据规模达到Tb至Pb量级。 实现这些海量数据的实时可视化和实时诊断分析面临很大的困难, 时空多尺度数据抽取可以解决这一瓶颈。 高精度数据的可视化结果展示需要高分辨率的显示设备, 并行的大屏幕显示技术是解决这一问题的有效手段。 地球科学数值模拟、 并行数据抽取和高分辨率显示都需要搭建高性能计算机集群。 本文在基于Lagrange插值的多维度、 多尺度、 多分辨率并行数据抽取算法的同时, 利用并行计算节点及LCD显示器, 基于Rocks cluster系统搭建起一个176核, 4×10×1024×1280分辨率的高性能计算模拟、 数据抽取和并行显示输出的集成平台, 并将该平台成功应用到气候模式模拟产生的海量数据的并行抽取和并行显示。     

局部最优形态参数的RBF分块地形插值方法与实验

径向基函数（Radial Basis Function,RBF）是一种不需对数据做任何假设,能准确逼近任意维度数据的空间插值方法。其特别适合于复杂地形的数字高程模型（Digital Elevation Model,DEM）插值重建,但随着已知点数量的增加会导致插值模型求解困难或缓慢。针对这个问题,本文基于二叉树自适应递归分块原理,采用局部最优形态参数的RBF分块插值方法进行DEM插值重建。首先,设定子区域最小点数阈值和相邻子区域的重叠率,自顶向下,对研究区域进行递归分块,构建区域分块二叉树,对二叉树叶子节点区域,采用逐点交叉验证（Leave One Out Cross Validation,LOOCV）方法求解其最优形态参数,建立局部RBF最优插值模型;然后,根据单元分解原理,采用加权平均方法对相应叶子节点区域内的待插值点高程进行加权求和,自底向上递归求解,得到待插值点最终高程值。以云南某地区DEM进行插值实验表明,采用本文方法进行DEM插值重建,稳定性较好,插值精度高。     

DEM内插算法对坡度的影响

本文以干旱、半干旱地区的渭干河-库车河三角洲绿洲为研究区域,利用样区的1:50000的矢量地形图为数据源,通过对不同插值算法生成的DEM计算坡度,坡度计算像元值分别采用25m、50m、100m、1.50m、200m,5个分级.通过对比研究发现:在坡度小于27°时,不同的算法对坡度具有显著的影响,而大干27°时,算法不会...     

弹道扰动引力的有限元逼近算法

为了克服多项式逼近弹道扰动引力的缺点，根据有限元插值的原理，采用了对弹道周围空间区域进行有限元剖分的方法，利用剖分单元各顶点的扰动引力分量内插出弹道点对应的扰动引力分量值。结果表明，文中提出的逼近算法能够快速精确可靠地逼近弹道扰动引力，是一种具有应用价值的方法。     

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

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

重庆风暴尺度集合预报系统降水预报性能分析

基于在重庆市气象局业务运行的对流可分辨尺度（3 km）集合预报系统,在已有初值扰动、模式物理过程扰动和侧边界扰动的基础上,对不同集合成员采用不同地形插值方案和地形平滑方案实现对模式静态地形高度的扰动,体现数值模式中地形转换过程的不确定性,开展集合预报批量平行试验。结果表明：（1）实现对模式静态地形高度的扰动后,各集合成员地形高度的离散度与实际地形的起伏程度对应关系较好,两者空间分布特征非常相似,地形较平坦的平原地区离散度较小,而地形较复杂的高原地区或山区离散度较大;（2）加入模式地形扰动方案后,集合扰动能量总体上有所增大,且低层比中、高层更明显,能量增幅在较短预报时效（12 h）内最显著,随着预报时效延长呈逐渐减小趋势,且能量增幅大值中心主要出现在地形较复杂、集合成员地形高度离散度较大地区;（3）模式地形扰动方案一定程度上能提高降水概率预报技巧及改进集合平均降水预报,在对高空要素和2 m温度、10 m风场等近地面要素的集合平均均方根误差和集合离散度无负面影响的前提下,能一定程度上优化集合分布。