一种利用Nyström离散与FFT快速褶积的散射地震波并行计算方法

利用数值方法解Lippermann-Schwinger (L-S)方程的主要困难在于系数矩阵存储和线性方程组求解.这主要是因为L-S方程的积分部分是一个空间褶积,在离散后将导致一个满秩矩阵,进而形成一个大型或超大型代数方程组.因此,在利用L-S解决地震波散射问题时,一般是利用散射级数法而非数值方法.然而,散射级数法的计算精度和收敛性强烈地依赖于速度扰动的强度,而克服这种依赖性的一个可能的途径就是对现有的数值方法进行改进或是建立新的数值求解方案.在这种思想指导下,首先对L-S方程进行改写,得到一个与原L-S方程等价的积分方程(等价L-S方程).然后,对等价L-S方程进行逐点归一化处理,并利用Nyström法对经归一化处理的等价L-S方程(归一化等价L-S方程)进行离散,并用FFT计算空间褶积.之所以这样选择是由于归一化等价L-S方程经Nyström法离散生成的系数阵为一个Toeplitz阵,可利用其Toeplitz性质降低存储空间;而FFT可以将矩矢空间褶积转化为乘积,且积分核部分只要计算一次即可.进一步,为节约正演计算时间,设计了进程级和线程级相结合的MPI+OpenMP并行模式.数值试验表明,与传统的积分方程数值算法相比,利用等价L-S方程、Nyström离散和FFT快速褶积的计算方案可极大地降低存储需求,进而在保证精度的同时提高计算效率.

     

分区统计地图上符号的视觉搜寻和认知工效研究

根据某地区的分区统计数据和有关图集,总结设计出4种方案的分区统计地图,在自行编制的认知测试软件上进行了目标搜寻测试。结果表明,在分区统计地图上搜寻目标时,视觉系统运用了平行加工策略和序列加工策略,4个方案的目标搜寻平均反应时间和正确率受到这两种策略的影响。     

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

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

二郎山公路隧道轴线地应力测试研究

介绍了二郎山公路隧道轴线上3个钻孔的地应力测试结果,应用水压致裂法在钻孔中测试,所测的两个不同大小的水平主应力都大于垂直主应力（除了地表以下很浅处外）,这个地应力与区域大地构造方向一致,在此基础上评价了隧道碉室围岩的稳定性。     

Analytical stress and displacement due to twin tunneling in an elastic semi‐infinite ground subjected to surcharge loads

The construction of twin tunnels at shallow depth has become increasingly common in urban areas. In general, twin tunnels are usually near each other, in which the interaction between tunnels is too significant to be ignored on their stability. The equivalent arbitrarily distributed loads imposed on ground surface were considered in this study, and a new analytical approach was provided to efficiently predict the elastic stresses and displacements around the twin tunnels. The interaction between 2 tunnels of different radii with various arrangements was taken into account in the analysis. We used the Schwartz alternating method in this study to reduce the twin‐tunnel problem to a series of problems where only 1 tunnel was contained in half‐plane. The convergent and highly accurate analytical solutions were achieved by superposing the solutions of the reduced single‐tunnel problems. The analytical solutions were then verified by the good agreement between analytical and numerical results. Furthermore, by the comparison on initial plastic zone and surface settlement between analytical solution and numerical/measured results of elastoplastic cases, it was proven that the analytical solution can accurately predict the initial plastic zone and its propagation direction and can qualitatively provide the reliable ground settlements. A parametric study was finally performed to investigate the influence of locations of surcharge load and the tunnel arrangement on the ground stresses and displacements. The new solution proposed in this study provides an insight into the interaction of shallow twin tunnels under surcharge loads, and it can be used as an alternative approach for the preliminary design of future shallow tunnels excavated in rock or medium/stiff clay.     

Analytical study of ground responses induced by the excavation of quasirectangular tunnels at shallow depths

The construction of quasirectangular tunnels at shallow depths is becoming increasingly common in urban areas to efficiently utilize underground space and reduce the need for backfilling. To clarify the mechanical mechanism of the stresses and displacements around the tunnels, this study proposes analytical solutions that precisely account for quasirectangular tunnel shapes, the ground surface, the tunnel depth, and the ground's elastic/viscoelastic properties. The Schwarz alternating method combined with complex variable theory is employed to derive the elastic solution, and convergent and highly accurate solutions are obtained by superposing the solutions in the alternating iterations. Based on the solution and the extended corresponding principle for the viscoelastic problem, the time-dependent analytical solutions for the displacement are obtained for the ground assuming any viscoelastic model. The analytical solutions agree well with the finite element method (FEM) numerical results for models that are completely consistent, and qualitatively agree with field data. Furthermore, based on the stress solution combined with the Mohr-Coulomb failure criterion, the predicted initial plastic zone and propagation directions around the tunnels are qualitatively consistent with those determined by the limit analysis. A parametric study is performed to investigate the influences of the rectangular/quasirectangular tunnel shape, burial depth, and supporting pressure on the ground stresses and displacements.     

Nonlinear wave-wave interactions and wedge waves

A tetrad mechanism for exciting long waves, for example edge waves, is described based on nonlinear resonant wave-wave interactions. In this mechanism, resonant interactions pass energy to an edge wave,from the three participating gravity waves. The estimated action flux into the edge wave can be orders of magnitude greater than the transfer fluxes derived from other competing mechanisms, such as triad interactions. Moreover, the numerical results show that the actual transfer rates into the edge wave from the three participating gravity waves are two- to three- orders of magnitude greater than bottom friction.     

时空大数据背景下并行数据处理分析挖掘的进展及趋势

随着互联网、物联网和云计算的高速发展,与时间、空间相关的数据呈现出“爆炸式”增长的趋势,时空大数据时代已经来临。时空大数据除具备大数据典型的“4V”特性外,还具备丰富的语义特征和时空动态关联特性,已经成为地理学者分析自然地理环境、感知人类社会活动规律的重要资源。然而在具体研究应用中,传统数据处理和分析方法已无法满足时空大数据高效存取、实时处理、智能挖掘的性能需求。因此,时空大数据与高性能计算/云计算融合是必然的发展趋势。在此背景下,本文首先从大数据的起源出发,回顾了大数据概念的发展历程,以及时空大数据的特有特征;然后分析了时空大数据研究应用产生的性能需求,总结了底层平台软硬件的发展现状;进而重点从时空大数据的存储管理、时空分析和领域挖掘3个角度对并行化现状进行了总结,阐述了其中存在的问题;最后指出了时空大数据研究发展趋势。     

A parallel method to accelerate spatial operations involving polygon intersections

Polygon intersection is an important spatial data-handling process, on which many spatial operations are based. However, this process is computationally intensive because it involves the detection and calculation of polygon intersections. We addressed this computation issue based on two perspectives. First, we improved a method called boundary algebra filling to efficiently rasterize the input polygons. Polygon intersections were subsequently detected in the cells of the raster. Owing to the use of a raster data structure, this method offers advantages of reduced task dependence and improved performance. Based on this method, we developed parallel strategies for different procedures in terms of workload decomposition and task scheduling. Thus, the workload across different parallel processes can be balanced. The results suggest that our method can effectively accelerate the process of polygon intersection. When addressing datasets with 1,409,020 groups of overlapping polygons, our method could reduce the total execution time from 987.82 to 53.66 s, thereby obtaining an optimal speedup ratio of 18.41 while consistently balancing the workloads. We also tested the effect of task scheduling on the parallel efficiency, showing that reducing the total runtime is effective, especially for a lower number of processes. Finally, the good scalability of the method is demonstrated.     

数字地形分析方法研究的维度——精准、高效、易用

数字地形分析基于数字高程模型计算各种地形信息,在诸多地理现象和地理过程的定量刻画和模拟中都发挥着不可或缺的重要作用。作为地理信息空间分析的一个基本组成内容、GIS的必备功能之一,数字地形分析的方法研究不断深化、扩展。本文从一个新的视角,即方法研究的维度,对该领域已有工作及发展趋势进行梳理,将数字地形分析方法研究工作归纳为依次显现出的3个研究维度：① 精准,即如何设计新的数字地形分析方法以更准确、精确地计算出符合地理学认知、满足应用领域问题求解所需的地形信息;② 高效,即如何使数字地形分析的计算更高效;③ 易用,即如何使数字地形分析对于应用者（尤其是相关应用领域中非数字地形分析专家的用户）更易用。特别是近年新显现出的“易用”研究维度,通过研究设计数字地形分析领域应用建模知识的形式化表达及相应的推理应用方法,有望实现数字地形分析应用工作流的智能化构建,满足实际应用中对数字地形分析的易用性需求;同时,其研究思路和方法也可为地理信息空间分析中诸多其他分支领域解决所面临的类似问题提供借鉴。