共享内存环境下的大块干涉相位图解缠算法

针对大块干涉相位图解缠效率低的问题,提出了一种共享内存环境下的大块干涉相位图解缠方法。利用质量引导和最小不连续优化策略设计了复合相位解缠算法,并在共享内存环境下对算法中的质量图计算、高低质量区域分割和低质量区域优化3个步骤进行了并行化。为充分利用多核计算资源,针对低质量区域的并行优化,设计了动态计算任务分配方法。对仿真和干涉合成孔径雷达相位图的解缠试验结果表明,所提算法在提高相位解缠精度的同时,也极大地优化了相位解缠效率。     

基于超算的全球电离层模型快速并行解算

为应对日益丰富的观测数据以及数据再处理对高性能计算的需求,开发了基于OpenMP以及MPI（Message Passing Interface）并行计算的全球电离层快速建模算法。采用武汉大学超级计算机对全球电离层建模效率进行了不同并行计算方案的实验。结果表明,采用多节点MPI并行计算能够极大地提高数据处理效率,相比传统单节点串行计算提高了近30倍,相比单节点OpenMP并行计算提高了近3~4倍。MPI并行计算方案充分利用了丰富的计算机资源来提高全球电离层建模效率,对电离层建模算法的快速测试、产品的重新再处理具有重要作用,对多系统全球导航卫星系统（Global Navigation Satellite System,GNSS）快速精密定轨、大规模GNSS网解也有较好的参考价值。     

非差模式的GNSS数据并行解算设计及实现

基于精密单点定位技术的非差模式是当前GNSS数据处理的主要策略之一。随着测站规模的增大,非差模式的处理时间也线性递增,传统的串行处理方法需消耗大量的计算时间。采用工厂模式和责任链模式实现了非差精密单点定位;利用轻量级的并行编程技术从底层设计并实现了基于任务的非差多核并行解算;进一步在网络多节点环境中建立并发布非差计算服务,实现了网络多节点协同并行解算GNSS数据。通过大量数据的测试与试验,验证了多核多节点的非差并行解算方案的高效性。试验结果表明,单节点多核并行、双节点网络并行、四节点网络并行、六节点网络并行的计算效率分别比单节点串行方案平均提高了2.74,5.30,9.38和14.69倍。     

iGMAS全球电离层延迟模型及并行计算策略

针对全球电离层延迟建模中传统串行处理方法效率低等问题,研究了基于全球分布的IGS跟踪站和iGMAS跟踪站观测数据实现全球电离层建模并行解算的基本方法、流程及策略。在Bernese软件基础上研制了一套iGMAS全球电离层延迟建模软件。为了验证并行解算方法的正确性和计算效率,利用全球200个左右IGS跟踪站和6个iGMAS跟踪站2014-08-20-2014-09-06共7周的观测数据,解算了快速电离层TEC格网。与IGS,CODE以及ESA最终电离层格网比较,结果表明:基于该方法解算的快速电离层TEC格网,与CODE,ESA以及IGS最终电离层TEC格网的互差,统计不同纬度带内偏差的均方根误差,全球范围内偏差的均方根误差均在1.5~2.5 TECu之间,南北半球高纬度地区在0.5~1.5 TECu之间,所有地区均优于5 TECu,整体精度与IGS,CODE以及ESA最终电离层TEC格网精度产品相当。     

利用并行计算方法实现地震活动速率变化参数Z值的空间扫描处理及其计算效能评价

近30多年的研究表明,地震平静可能是中长期的地震前兆现象,而地震平静的内在特征体现了地震活动速率的变化. Z 检验法是检验地震活动速率变化应用最为广泛的方法之一,其度量值即 Z 值. 由于 Z 值空间扫描中采用等地震数的空间扫描窗口,每个空间扫描窗口的大小通常都是变化的,这样带来了大计算量的问题,因此,本文在自行组建的4节点PC机群上实现了 Z 值空间扫描的并行化处理. 初步测试结果说明,并行化处理中获取了很高的加速比和并行计算效率,可以将此种方法应用于其他地震活动性参数的空间扫描处理.     

桥梁-桩-土相互作用三维地震反应分析方法

利用大型有限元通用软件ANSYS/LS-DYNA及其并行运算功能,提出了一种大跨度斜拉桥在地震作用下考虑SSI效应的三维整体时程反应分析方法,比较了考虑桥梁结构线性和非线性两种情形的计算结果,并对其计算效率进行了探讨。     

重叠区域伪谱法计算非均匀地球介质地震波传播

伪谱法是一种高效、高精度计算非均匀介质地震波传播的数值算法,由于它的微分算子的全局性,一般认为该方法不适于并行计算. 本文介绍了并行计算非均匀介质中地震波传播的重叠区域分解算法,给出了一种基于傅里叶伪谱法的并行算法. 文中给出的算法将介质划分为相互重叠的若干区域,在各个子域上单独求解,利用重叠部分的解的传递,将各个子域连接起来,实现了伪谱法在分布式并行处理机上的计算. 文中给出了一个将二维区域分解的算例,比较了并行算法和整体算法的结果,分析了并行算法的计算精度. 结果表明,并行算法会有效降低计算时间,并且保证计算精度. 该方法在大规模三维非均匀介质的地震波场模拟方面有应用价值.      

采用动态负载均衡的LiDAR数据生成DEM并行算法

随着高性能计算的发展,并行技术已经广泛应用于LiDAR数据的分析处理。本文针对现有LiDAR数据生成DEM并行算法所存在的负载不均衡问题,设计并实现了动态负载均衡的LiDAR数据生成DEM并行算法。该算法采用主从式并行策略,管理进程负责LiDAR点云的高效自适应条带划分,计算进程负责LiDAR点云生成DEM的计算。本文设计了任务量的动态调度策略:首先,由所有进程并行创建任务量由大到小排列的待处理任务队列;然后,管理进程根据计算进程的反馈对待处理任务进行动态分配,以达到负载均衡。在24 核集群环境下,用30 GB(约12 亿点)LiDAR数据对本文算法进行测试,生成分辨率为1 m的格网DEM,算法加速比峰值达到15.16;同时,与静态调度策略进行对比实验,结果显示本文的动态负载均衡策略可更好地保证进程间的负载均衡,有效地提高了LiDAR生成DEM并行算法的整体效率。     

二维矢量线符号在三维地形表面的贴合渲染方法

二维矢量符号是二维地图的重要组成部分。如何将二维矢量符号准确地贴合渲染在三维地形模型表面是地理信息科学的研究热点之一。针对线状矢量要素（简单线和周期线）中现有方法存在的周期线周期分布不均、渲染精度低、性能受地形模型复杂度影响大、矢量要素与地表贴合不紧密等问题,本文提出了一种基于屏幕像素进行投影反算的渲染方法,以及基于地表延伸长度的周期判断方法。该方法通过预处理计算线要素贴合在地表的真实长度,并将线要素按其实际宽度进行面域化处理;实时计算每个屏幕像素在二维矢量平面内的覆盖范围,通过该范围与矢量面域的位置关系筛选位于其内部的像素;进一步通过线的实际长度进行周期判断以确定像素色彩值。整个过程通过CPU-GPU并行异构计算模型实现,提高了场景交互性能。实验表明,本方法具有周期分布均匀、渲染精度高、贴合结果真实准确和性能不受地形模型影响等优点。     

基于MPICH环境下的复杂重力异常体并行正演模拟

任意多面体重力异常正演公式常用于解决复杂几何形态地质体的正演问题。 本文以均匀物性多面体重力异常正演公式为基础, 应用有限元技术中的网格离散化思想, 以任意四面体为基本单元, 通过并行计算技术在MPICH环境下实现了任意连续空间物性分布复杂异常体网格模型的重力异常正演模拟, 通过并行处理可以有效加速正演计算速度。 本文研究结果对于联合重力异常场正演建模和开展复杂模型网格的重力场计算有一定参考意义。