并行计算技术及其在勘探地球物理学中的现状与展望

本文详细阐述了并行计算技术及其在地球物理勘探数据处理中的发展现状和发展趋势,分析了几个代表性的并行算法实例.这些结果表明,在拥有强大的并行机的基础上,基于并行计算开发环境(MPI和PVM等)设计高效的并行算法,通过分配合理的并行粒度、通信开销、负载平衡等执行高效的并行计算,可以有效加快处理速度、降低成本.目前,并行算法在地震数据处理中应用已较为成熟,近年来向更实用的基于PC机群的并行技术发展.然而,在非地震方法中,并行算法应用较少见文献报道,研究尚处于初级研究阶段.在大地电磁的二维和三维正、反演问题上,并行计算技术逐渐得到越来越多关注和重视.随着资源和能源需求的增长,地球物理勘探向深度和广度快速发展,大幅增长的数据量使得高性能并行计算机和高效的并行算法在勘探地球物理学中的发展和应用将占据愈来愈重要的地位.     

基于数字信号处理器的Linux网络远程实时监控系统的实现

针对现有市场上传统视频监控系统的不足,以及大部分数字信号处理器在控制和内存管理方面的局限性,设计了一个面向网络设备的嵌入式监控模型。该模型能最大限度的发挥数字信号处理器的性能,在单片数字信号处理器下,实现网络远程实时监控。首先对该监控系统的设计方案进行阐述;接着针对该系统模型的重点部分,数字信号处理器下的μCLinux操作系统移植,视频压缩模块和实时传输模块实现及优化进行了具体设计;最后对该模型的整体性能进行测试,证明该模型能满足网络远程实时监控系统的需要。     

随机介质理论预测近距离平行盾构引起的地表沉降

提出近距离双线水平平行盾构施工引起的总的地表沉降曲线符合正态分布规律。基于单线随机介质理论简化公式,针对现有不足,建立修正的随机介质理论简化计算公式,可以计算近距离双线平行盾构施工引起的地表沉降,同时能够考虑沉降曲线的不对称性,对计算参数取值进行了探讨。算例分析结果表明,文中方法预测值与实测值比较吻合,且计算比较简便。实测数据统计结果表明,先行隧道开挖后双线隧道的半径收敛值比先行隧道半径收敛值要大,而主要影响角则变小;提出可直接用于近距离界定的式（6）,认为双线隧道的地表沉降曲线形状与P值大小有关,近距离的界定范围为 。     

藏北羌塘盆地晚三叠世地层特征与对比

羌塘盆地上三叠统日干配错组、土门各拉组和肖茶卡组、那底岗日组分别分布于南、北羌塘凹陷。在以往研究和笔者在中国地质调查局青藏重点沉积盆地油气资源潜力分析项目工作基础上,通过对各地层单元的详细研究和区域对比,从古生物、岩石组合、沉积系列、地层单元的接触关系、年代地层进行了综合分析。结果显示北羌塘凹陷那底岗日组与南羌塘凹陷日干配错组和土门各拉组的时代主要为晚三叠世诺利期—瑞替期,三者为羌塘晚三叠世诺利期—瑞替期裂陷盆地形成过程中的同期异相的火山-沉积系列。北羌塘凹陷肖茶卡组的是残留陆表海沉积产物,其时代主要为晚三叠世卡尼期-诺利期早期,与那底岗日组、日干配错组和土门各拉组均为上下对比关系。     

High-performance three-horizon composition algorithm for large-scale terrains

This work presents a high-performance algorithm to compute the horizon in very large high-resolution DEMs. We used Stewart's algorithm as the core of our implementation and considered that the horizon has three components: the ground, near, and far horizons. To eliminate the edge-effect, we introduced a multi-resolution halo method. Moreover, we used a new data partition approach, to substantially increase the parallelism in the algorithm. In addition, several optimizations have been applied to considerably reduce the number of arithmetical operations in the core of the algorithm. The experimental results have demonstrated that by applying the above-described contributions, the proposed algorithm is more than twice faster than Stewart's algorithm while maintaining the same accuracy.     

A general-purpose parallel raster processing programming library test application using a geographic cellular automata model

A general-purpose parallel raster processing programming library (pRPL) was developed and applied to speed up a commonly used cellular automaton model with known tractability limitations. The library is suitable for use by geographic information scientists with basic programming skills, but who lack knowledge and experience of parallel computing and programming. pRPL is a general-purpose programming library that provides generic support for raster processing, including local-scope, neighborhood-scope, regional-scope, and global-scope algorithms as long as they are parallelizable. The library also supports multilayer algorithms. Besides the standard data domain decomposition methods, pRPL provides a spatially adaptive quad-tree-based decomposition to produce more evenly distributed workloads among processors. Data parallelism and task parallelism are supported, with both static and dynamic load-balancing. By grouping processors, pRPL also supports data–task hybrid parallelism, i.e., data parallelism within a processor group and task parallelism among processor groups. pSLEUTH, a parallel version of a well-known cellular automata model for simulating urban land-use change (SLEUTH), was developed to demonstrate full utilization of the advanced features of pRPL. Experiments with real-world data sets were conducted and the performance of pSLEUTH measured. We conclude not only that pRPL greatly reduces the development complexity of implementing a parallel raster-processing algorithm, it also greatly reduces the computing time of computationally intensive raster-processing algorithms, as demonstrated with pSLEUTH.     

Massively parallel spatial point pattern analysis: Ripley’s K function accelerated using graphics processing units

This study presents a massively parallel spatial computing approach that uses general-purpose graphics processing units (GPUs) to accelerate Ripley’s K function for univariate spatial point pattern analysis. Ripley’s K function is a representative spatial point pattern analysis approach that allows for quantitatively evaluating the spatial dispersion characteristics of point patterns. However, considerable computation is often required when analyzing large spatial data using Ripley’s K function. In this study, we developed a massively parallel approach of Ripley’s K function for accelerating spatial point pattern analysis. GPUs serve as a massively parallel platform that is built on many-core architecture for speeding up Ripley’s K function. Variable-grained domain decomposition and thread-level synchronization based on shared memory are parallel strategies designed to exploit concurrency in the spatial algorithm of Ripley’s K function for efficient parallelization. Experimental results demonstrate that substantial acceleration is obtained for Ripley’s K function parallelized within GPU environments.     

全国前兆台网“九五”系统台站接入的设计与实施

介绍了“全国前兆台网‘九五’系统台站接入改造”项目的技术方案,即采用接入地震前兆通信协议转换器或更换“九五”前兆公用数据采集器方式将“九五”前兆仪器接入“十五”前兆系统,并详述了该方案的技术特点;介绍了方案的部署实施情况,测试平台运行情况以及通过对测试平台的监控而发现并解决的问题.并行运行结果显示,接入改造系统运行稳定,产出数据可靠.     

施甸盆地三维地震动研究

本文利用三维施甸盆地模型以及川滇地区的地壳速度结构,采用谱元法和并行计算技术,考虑地形的影响,模拟了盆地外围中等设定地震作用下施甸盆地地震波的传播,分析了面波的产生和传播,研究了该盆地不同位置剖面地震动的差异和原因.研究表明,小型盆地在中等地震的作用下盆地内部地震动空间差异很大,盆地中较深的小型凹陷更容易汇聚面波,形成更强的振动,位于覆盖层厚度剧烈变化区域的地震动强度也会很大,甚至大于覆盖层厚的区域.较基岩地震动,盆地内地震动持续时间明显增加,地震动的放大作用明显,这也是盆地震害严重的主要原因之一.     

大区域地下水模拟的预优并行GMRES(m)算法研究

大区域研究区由于涉及范围大、水文地质参数复杂多变,一直是进行地下水数值模拟的热点和难点。针对大区域地下水模拟的特点,在MPI环境中对Krylov子空间GMRES(m)算法的并行性进行分析,提出基于区域分解法的并行实现策略,并对不同的预条件子的加速效果进行比较。数值实验结果表明:并行GMRES(m)算法在求解大区域三维地下水模型时可以显著的加快求解速度,且具有较好的可扩展性。另外,Jacobi预条件子与GMRES算法的组合具有更优的加速比和执行效率,是一种求解大型化、复杂化地下水水流问题的可行方案。