月球构造新认识：基于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.     

Investigation of highly efficient algorithms for solving linear equations in the discontinuous deformation analysis method

Large‐scale engineering computing using the discontinuous deformation analysis (DDA) method is time‐consuming, which hinders the application of the DDA method. The simulation result of a typical numerical example indicates that the linear equation solver is a key factor that affects the efficiency of the DDA method. In this paper, highly efficient algorithms for solving linear equations are investigated, and two modifications of the DDA programme are presented. The first modification is a linear equation solver with high efficiency. The block Jacobi (BJ) iterative method and the block conjugate gradient with Jacobi pre‐processing (Jacobi‐PCG) iterative method are introduced, and the key operations are detailed, including the matrix‐vector product and the diagonal matrix inversion. Another modification consists of a parallel linear equation solver, which is separately constructed based on the multi‐thread and CPU‐GPU heterogeneous platforms with OpenMP and CUDA, respectively. The simulation results from several numerical examples using the modified DDA programme demonstrate that the Jacobi‐PCG is a better iterative method for large‐scale engineering computing and that adoptive parallel strategies can greatly enhance computational efficiency. Copyright © 2015 John Wiley & Sons, Ltd.     

聚合物驱组合段塞分注实验研究

萨尔图油田南二区东部140号断层以东地区葡11—4油层，以多段多韵律和正韵律沉积为主，平均孔隙度30．0％，平均空气渗透率1980×10^-3μm2。水驱结束时该油层平均含水92％，采出程度39．5％。计划对该区块实施聚合物驱开采。以往经验认为，组合段塞调整注聚技术能够有效改善流度控制作用。然而对组合注入情况下如何进一步提高驱油效果、注入井分注时机、层段注入强度对开发效果有何影响却很少探讨。通过室内驱油实验，在聚驱的组合段塞注入不同时机（水驱空白阶段、含水下降阶段、含水稳定阶段、含水回升阶段）下调整注入强度实施分注以及在含水回升分注基础上进一步调整注入强度，同时对比不分注实验，评价各方案驱油效果。实验表明：不实施分注时聚驱采收率在水驱基础上提高了19．07个百分点，而组合段塞注入时实施分注均能够在不分注基础上进一步提高驱油效果；空白水驱时实施分注采收率提高幅度最大，为26．55个百分点；随着分注时机的滞后采收率提高幅度增加值依次递减，分注时机越早越好。     

关于精密地形改正的若干注记

引用斜率作为新的变量，首次使地形改正一类奇异积分非奇异，并在此基础上详细讨论了精密地形改正计算的若干技术问题，地形模型和某山区的试算结果均表明，以往线性近似误差较大，应避免使用。     

对可控硅调压整流电路若干问题的分析

本文比较了两种可控硅整流电路,在这两种电路中把反并联可控硅接入阳极变压器。文章介绍了电路的过压和过流的保护方法并给出了电路的计算公式和经验数据。     

利用地层空间的相对稳定性、井资料、地震资料模拟高分辨率波阻抗资料

由于地震资料是带限的,而井资料约束反演过程中的地震反演频带的展宽也是有一定限度的,使得反演地震波阻抗也是带限的。如何进一步拓宽反演地震波阻抗的频带宽度,提高地震资料的分辨率是文中研究的内容。利用井约束下的地震反演波阻抗资料、地震资料及井资料,借鉴地质模拟技术,研究如何利用地层的空间相对稳定性特征,通过相关函数的分析,估计反演更高分辨率的地震波阻抗资料。考虑到地质体空间的相关性及分辨率,相关半径取相关函数从最大点到第1个极值点的范围。在目标函数的计算过程中,既考虑当前计算点的权系数信息,又考虑井外推时的权系数信息,这样使计算结果更加稳定可靠。通过综合模拟计算,能够识别较小砂体     

基于R树的分布式并行空间索引机制研究

为提高分布式并行计算环境下海量空间数据管理与并行化处理的效率,基于并行空间索引机制的研究,设计一种多层并行R树空间索引结构。该索引结构以高效率的并行空间数据划分策略为基础,以经典的并行计算方法论为依据,使其结构设计在保证能够获得较好的负载平衡性能的前提下,更适合于海量空间数据的并行化处理。以空间范围查询并行处理的系统响应时间为性能评估指标,通过实验证明并行空间索引结构具有设计合理、性能高效的特点。     

光学测微器的原理、检定及误差分析

着重讨论了高精度光学测微器的基本工作原理、检定方法及误差分析。