多功能水文频率分析计算软件的研制和应用

介绍了多功能水文频率分析计算软件的设计及功能扩充，尤其是增加了采用鼠标拉动曲线直接屏幕适线的功能。本软件在生产实践中，获得了较好的应用效果，不仅功能多，而且实用方便。     

Numerical simulations for coupled rock deformation and gas leak flow in parallel coal seams

Based on the new viewpoint of interaction mechanics for solid and gas, gas leakage in parallel deformable coal seams can be understood. That is, under the action of varied geophysical fields, the methane gas flow in a double deformable coal seam can be essentially considered to be compressible with time-dependent and mixed permeation and diffusion through a pore-cleat deformable, heterogeneous and anisotropic medium. From this new viewpoint, coupled mathematical models for coal seam deformation and gas leak flow in parallel coal seams were formulated and the numerical simulations for slow gas emission from the parallel coal seams are presented. It is found that coupled models might be close to reality. Meanwhile, a coupled model for solid deformation and gas leak flow can be applied to the problems of gas leak flow including mining engineering, gas drainage engineering and mining safety engineering in particular the prediction of the safe range using protective layer mining where coal and gas outbursts can efficiently be prevented. This revised version was published online in July 2006 with corrections to the Cover Date.     

Parallelized combined finite-discrete element (FDEM) procedure using multi-GPU with CUDA

This paper focuses on the efficiency of finite discrete element method (FDEM) algorithmic procedures in massive computers and analyzes the time-consuming part of contact detection and interaction computations in the numerical solution. A detailed operable GPU parallel procedure was designed for the element node force calculation, contact detection, and contact interaction with thread allocation and data access based on the CUDA computing. The emphasis is on the parallel optimization of time-consuming contact detection based on load balance and GPU architecture. A CUDA FDEM parallel program was developed with the overall speedup ratio over 53 times after the fracture from the efficiency and fidelity performance test of models of in situ stress, UCS, and BD simulations in Intel i7-7700K CPU and the NVIDIA TITAN Z GPU. The CUDA FDEM parallel computing improves the computational efficiency significantly compared with the CPU-based ones with the same reliability, providing conditions for achieving larger-scale simulations of fracture.     

并行处理技术在地理信息数据处理中的应用

随着地理信息存储量的飞速增长,传统的单进程、集中式的数据处理方式已不能满足基于网络的地理信息服务的效能要求。分析对比了OpenMP,MPI和MapReduce等主流并行编程模式,将关系型数据库与分布式空间数据管理系统相结合,提出了面向并行处理的地理信息存储模型和数据组织模型,将该模型与传统模型进行了对比分析,并基于MapReduce实现了地理空间数据并行处理框架,选取了矢量数据装载、影像数据装载以及数据切片作为典型数据处理案例开展对比实验,该技术方案的处理效率均数倍于传统技术方案。实验表明,该模型能够很好地支持并行处理框架,可为分布式环境下数据处理中心构建提供一个有效解决方案。     

月球构造新认识：基于60m高精度全月月形图的月球构造形态学研究

一直以来月球上都没有发现像地球上一样的全球板块构造现象,被认为是属于单板块构造的行星体并且是不活动的。随着月球地形探测精度的提高,人们逐渐发现了一些地质现象表明月球并不像之前想象的一样是不活动的。高精度的月形数据显示出了月球上除了存在众所周知的(撞击作用产生的)复杂的撞击构造外,仍然存在许多可能由内生作用产生的构造特征:如裂谷、山脊、火山群等。本文利用最新高精度月球激光高度数据(LOLA数据)建立了全月60m(相当于7秒网格)高精度月形图。由于数据量超大(TB级别),作者编写了基于并行架构的集群计算处理代码,其计算效率提高了近300倍。基于这些图件,作者在全球范围内找出许多重要的构造现象,并根据其构造及地质特征进行分类研究。研究结果表明并不像之前认识一样,月球上的确存在许多地质构造体是由月球的内生地质作用所产生,并建立对于月球结构和构造的新认识。此外,通过对比发现,像地球上全球尺度的线性山脉只是板块构造运动的结果;另一方面,在无板块构造的行星上,就像月球一样,会产生大面积的高地(月陆)和低地(月海地区),这种地质现象很有可能是行星地壳岩浆演化的直接结果。在文章最后,作者还将制作出的60m分辨率的全月月形系列图件附后并公开,以方便其他科研人员进一步开展相关研究;并以此敬献中国地质科学院60周年庆典。     

A parallel algorithm for coverage optimization on multi-core architectures

Location siting is an important part of service provision, with much potential to impact operational efficiency, safety, security, system reliability, etc. A class of location models seeks to optimize coverage of demand for service that is continuously distributed across space. Decision-making and planning contexts include police/fire resource allocation for a community, siting cellular towers to support cell phone signal transmission, locating emergency warning sirens to alert the public of severe weather and other related dangers, and many others as well. When facilities can be sited anywhere in continuous space to provide coverage to an entire region, this is a very computationally challenging problem to solve because potential demand for service is everywhere and there are an infinite number of potential facility sites to consider. This article develops a new parallel solution approach for this location coverage optimization problem through an iterative bounding scheme on multi-core architectures. The developed approach is applied to site emergency warning sirens in Dublin, Ohio, and fire stations in Elk Grove, California. Results demonstrate the effectiveness and efficiency of the proposed approach, enabling real-time analysis and planning. This work illustrates that the integration of cyberinfrastructure can significantly improve computational efficiency in solving challenging spatial optimization problems, fitting the themes of this special issue: cyberinfrastructure, GIS, and spatial optimization.     

A new model of edge-to-edge contact for three dimensional discontinuous deformation analysis

A robust contact theory can be regarded as key to three dimensional discontinuous deformation analyses (3D-DDA). Not only must this theory provide an efficient algorithm to judge the type and location of contacts but also it must be able to present comprehensive formulations for every kind of contact (open, sliding and locked contact). There are six types of contact in three dimensional discontinuous deformation analyses (vertex-to-vertex, vertex-to-edge, vertex-to-face, edge-to-edge, edge-to-face and face-to-face) that can be converted to vertex-to-face and edge-to-edge contacts. This paper presents a new model of edge-to-edge contact to three dimensional discontinuous deformation analyses (3D-DDA). This new model considers both kinds of edge-to-edge contact (cross-over and parallel edge-to-edge contact) and presents a criterion for inter-penetration. Sub matrices of normal and shear spring and friction force are derived by geometrical analysis and penalty method. This new model is implemented in a 3D-DDA computer programme, and the example results demonstrate the validity of the model.     

科学计算程序库WHIGG F90 LIB及其地学应用

WHIGG F90LIB(WFL)是一款用FORTRAN90语言编写的用于地球学科的计算程序库。应用WFL可以减少重复性编程,方便地解决科学计算中数据的预处理问题,结合多种分析方法,对科学问题进行定量分析。利用WFL在地球重力场、气候涛动和数字信号处理3个方面的应用实例,以验证WFL在科学计算中应用的有效性。     

地震层析成像LSQR算法的并行化

讨论了地震层析成像的LSQR算法(最小二乘QR分解). 在建立偏导数矩阵方程组时,对区内地震在方程中保留震源项,引入正交投影算子进行参数分离,对区外远震采用传统的平滑处理方式,用LSQR法求解联立的方程组. 由于区内地震的正交分解处理和区外远震的平滑处理,使得偏导数矩阵中的非零元素成倍增加,对于大型反演问题,这些非零元素常常达到几十GB到几百GB的数量级,巨量的内存占用成为LSQR算法的瓶颈. 针对这一问题,本文研究了偏导数矩阵中非零元素的分布规律,设计出合理的存储结构,采用分布式存储进行矩阵计算,提出了LSQR算法的并行化方案,并在联想深腾6800超级计算机上实现. 导出了LSQR算法的并行效率估算公式. 对两个地区的实际地震层析成像数据进行了效率测试.     

钢筋混凝土柱考虑损伤累积的反复荷载-位移关系分析

为能在反复荷载作用下钢筋混凝土柱的荷载-位移关系分析中考虑柱低周疲劳性能,提出了一损伤模型,对柱中纵向受力钢筋和混凝土的损伤状态作评估与记录。将这一记录材料性能信息的损伤指标带入材料各自的恢复力模型以考虑产生损伤后材料的强度和刚度退化。基于多弹簧模型对不同变幅加载路径下及等幅低周疲劳加载下钢筋混凝土柱的空间反应进行了数值计算模拟。与已有试验结果比较表明,所提材料层次上的损伤累积模型以及考虑损伤累积效应的柱构件空间荷载-位移关系分析方法具有一定的精度,为钢筋混凝土柱的抗震性能分析提供了一个辅助工具。