地震应急灾情服务进展

利用Python语言编写了地震应急专题图快速产出的脚本程序,当地震发生后,由地震应急值班人员启动该脚本程序,该系统能够自动获取地震灾情评估系统产出的地震三要素、地震影响场、震中距等,并调用预制的地震应急指挥专题图模板,能够在震后5分钟内实现地震专题地图的快速产出,有效地解决了地震制图难、效率慢、质量不高的问题,能够为地震应急救援提供有力的技术支持。     

三维粘弹介质地震波场有限差分并行模拟

有限差分法在三维粘弹性复杂介质正演模拟地震波的传播中对计算机内存和计算速度要求比较高,单个PC机或工作站只能计算较少网格内短时间的波场。本文介绍一种基于MPI的并行有限差分法,可在PCCluster上模拟较大规模三维粘弹性复杂介质中地震波传播时的波场;可预测地震波在此类条件下传播时的运动学和动力学性质。对于更好地理解波动传播现象,解释实际地震资料及反问题的解决等均具有重要的理论与实际意义。     

基于Python的地震时大数据拥堵实时传输系统设计

针对当前地震时数据传输系统存在时延高、安全性较差及数据传输拥堵等问题,提出并设计基于Python的地震时大数据拥堵实时传输系统。系统硬件主要由视频转码模块、Python交互平台、数据流传输模块和系统维护模块几部分组成。利用解协议、解封装、解码、视频数据存储、编码、封装等实现视频转码;利用Python交互平台实现测点与地震数据监测中心之间的交互;结合CDMA传输方式与VPDN组网方式作为无线传输模块;由区域中心运维管理和通知下发等单元组成系统维护模块。系统软件分为收包与发包两种数据传输模式,以收包为主完成系统设计。实验结果表明,该系统响应速度快、延迟低,且安全性较为优越,解决了地震数据拥堵的问题。     

大地电磁三维快速松弛反演并行算法研究(英文)

为了进一步提高大地电磁三维快速松弛反演的计算效率,在深入分析大地电磁三维快速松弛反演算法的基础上,结合MPI自身的优越性,确定了并行计算的思路,实现了三维快速松弛反演的并行计算。通过理论模型合成数据和实测数据对实现的三维快速松弛反演并行程序进行了试算,分析对比了在多种情况下程序的执行效率。测试结果表明,所实现的三维快速松弛反演并行程序运行结果正确,效率提高明显。此成果在普通微机上实现,推动了大地电磁三维反演技术的实用化,可为其它地球物理三维正反演研究所借鉴。     

Parallel numerical simulation of two-phase flow model in porous media using distributed and shared memory architectures

A two-phase (water and oil) flow model in a homogeneous porous media is studied, considering immiscible and incompressible displacement. This model is numerically solved using the Finite Volume Method (FVM) and we compare four numerical schemes for the approximation of fluxes on the faces of the discrete volumes. We describe briefly how to obtain the mathematical and computational models applying axiomatic formulations and generic programming. Two strategies of parallelization are implemented in order to reduce the execution time. We study distributed (cluster of CPUs) and shared (Graphics Processing Units) memory architectures. A performance comparison of these two architectures is done along with an analysis of the four numerical schemes, for a water-flooding five-spot pattern model.     

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

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

大地电磁三维数据空间反演并行算法研究

目前大地电磁三维反演实际应用的主要问题是计算效率低.在对大地电磁三维数据空间反演算法进行深入分析的基础上,本文提出了基于频点和矩阵划分的大粒度并行反演方案和具体实现步骤,并在曙光TC5000A高性能计算平台上实现了基于MPI的大地电磁三维数据空间反演并行算法.该算法实现了包括三维正演、灵敏度矩阵、叉积矩阵以及模型改正量的并行执行,不仅计算效率高,而且每个节点机上灵敏度矩阵的存储空间只需原来微机上的2/N(N是参加并行计算的节点机个数),大大地减少了内存开销.通过两个理论模型合成的数据对实现的三维数据空间反演并行算法进行试算,对比分析了多个节点机下程序的执行效率.测试结果表明,所实现的三维数据空间反演并行算法是可行的、高效的,与单机相比,不仅可以提高运行速度,缩短计算时间,而且还可以扩大计算规模,极大地推动了大地电磁三维反演的实用化.     

Evapotranspiration modelling from climatic data using a neural computing technique

Evapotranspiration is one of the basic components of the hydrologic cycle and essential for estimating irrigation water requirements. This paper investigates the modelling of evapotranspiration using the feed‐forward artificial neural network (ANN) technique with the Levenberg–Marquardt (LM) training algorithm. The LM algorithm has never been used in evapotranspiration estimation before. The LM is used for the optimization of network weights, since this algorithm is more powerful and faster than the conventional gradient descent. Various combinations of daily climatic data, i.e. wind speed, air temperature, relative humidity and solar radiation, from three stations in Los Angeles, USA, are used as inputs to the ANN so as to evaluate the degree of effect of each of these variables on evapotranspiration. A comparison is made between the estimates provided by the ANN and those of the following empirical models: Penman, Hargreaves, Turc. Mean square error, mean absolute error and determination coefficient statistics are used as comparing criteria for the evaluation of the models' performances. Based on the comparisons, it was found that the neural computing technique could be employed successfully in modelling evapotranspiration process from the available climatic data. The results also indicate that the Hargreaves method provides better performance than the Penman and Turc methods in estimation of the evapotranspiration. The accuracy of the ANN technique in evapotranspiration estimation using nearby station data was also investigated. Copyright © 2007 John Wiley & Sons, Ltd.     

利用PETSc和FEPG编制海啸数值模拟程序的研究

海啸的数值模拟是海啸研究的一个重要领域, 它对于帮助理解海啸的基本物理特性和预防减灾具有重要意义。 海啸数值模拟程序的编制是一项繁杂的工作, 该文介绍了利用两种功能强大的通用软件(PETSc和FEPG)来进行海啸数值模拟程序编写的研究。 PETSc和FEPG采用有限差分、 有限元、 有限体积等多种离散方法, 可以对数值问题给出稳定的求解。 该文采用了直角坐标系下的非线性浅水波方程作为海啸波传播的控制方程进行离散求解, 并将其结果与TUNAMI N1模型进行了比较, 表明这两种方法方便而且有效。     

Modelling mean annual sediment yield using a distributed approach

In this paper a spatially distributed model for the calculation of sediment delivery to river channels is presented (SEDEM: SEdiment DElivery Model). The model consists of two components: (1) the calculation of a spatial pattern of mean annual soil erosion rates in the catchment using a RUSLE (Revised Soil Erosion Equation) approach; and (2) the routing of the eroded sediment to the river channel network taking into account the transport capacity of each spatial unit. If the amount of routed sediment exceeds the local transport capacity, sediment deposition occurs. An existing dataset on sediment yield for 24 catchments in central Belgium was used to calibrate the transport capacity parameters of the model. A validation of the model results shows that the sediment yield for small and medium sized catchments (10–5000 ha) can be predicted with an average accuracy of 41 per cent. The predicted sediment yield values with SEDEM are significantly more accurate than the predictions using a lumped regression model. Moreover a spatially distributed approach allows simulation of the effect of different land use scenarios and soil conservation techniques. Copyright © 2001 John Wiley & Sons, Ltd.     

Glacier mass balance simulation using SWAT distributed snow algorithm

Application of a temperature-index melt model incorporated into the Soil and Water Assessment Tool (SWAT) is presented to simulate mass balance (MB) and equilibrium line altitude (ELA) of three glaciers. The snow accumulation/melt parameters were adjusted to glacierized and free glacier areas, respectively. The SWAT snow algorithm enabled us to consider spatial variation of snow parameters by elevation bands across the sub-basins, while in the previous studies using SWAT, the related parameters were constant for an entire basin. The results show slight improvement in runoff simulation and significant improvement in simulated MB when considering ELA in model calibration. The results showed that SWAT can be applied to simulate MB, vertical MB distribution and annual ELA, with light calibration efforts for data-scarce catchments. The accuracy of the results depends on the modelled area of ablation zone from which most of the meltwater is released.     

Analysis of subsurface contaminant migration and remediation using high performance computing

Highly resolved simulations of groundwater flow, chemical migration and contaminant recovery processes are used to test the applicability of stochastic models of flow and transport in a typical field setting. A simulation domain encompassing a portion of the upper saturated aquifer materials beneath the Lawrence Livermore National Laboratory was developed to hierarchically represent known hydrostratigraphic units and more detailed stochastic representations of geologic heterogeneity within them. Within each unit, Gaussian random field models were used to represent hydraulic conductivity variation, as parameterized from well test data and geologic interpretation of spatial variability. Groundwater flow, transport and remedial extraction of two hypothetical contaminants were made in six different statistical realizations of the system. The effective flow and transport behavior observed in the simulations compared reasonably with the predictions of stochastic theories based upon the Gaussian models, even though more exacting comparisons were prevented by inherent nonidealities of the geologic model and flow system. More importantly, however, biases and limitations in the hydraulic data appear to have reduced the applicability of the Gaussian representations and clouded the utility of the simulations and effective behavior based upon them. This suggests a need for better and unbiased methods for delineating the spatial distribution and structure of geologic materials and hydraulic properties in field systems. High performance computing can be of critical importance in these endeavors, especially with respect to resolving transport processes within highly variable media.©1998 Elsevier Science Limited. All rights reserved     

大地电磁三维交错网格有限差分数值模拟的并行计算研究

为了更有效的提高大地电磁三维正演的计算速度,引入了并行处理技术.大地电磁三维交错网格有限差分数值模拟是按照不同频率来计算的,各频率之间求取电磁场值的过程是相互独立的.根据这一特点,可以将多个频率的计算任务平均划分为一个或者几个频率的计算子任务,分配到各个计算节点去并行执行,计算完成后将结果汇总.本文通过采用主从并行模式、分频并行计算的方案,在曙光TC5000A高性能并行平台上实现了基于MPI的大地电磁三维正演的并行计算.通过两个理论模型对实现的大地电磁三维正演并行算法进行试算,对比分析了多个节点机下程序的执行效率.测试结果表明,所实现的三维正演并行算法是正确的、高效的,为进一步的大地电磁三维反演并行算法研究奠定了重要基础.     

A new approach for computing a flood vulnerability index using cluster analysis

A Flood Vulnerability Index (FloodVI) was developed using Principal Component Analysis (PCA) and a new aggregation method based on Cluster Analysis (CA). PCA simplifies a large number of variables into a few uncorrelated factors representing the social, economic, physical and environmental dimensions of vulnerability. CA groups areas that have the same characteristics in terms of vulnerability into vulnerability classes. The grouping of the areas determines their classification contrary to other aggregation methods in which the areas' classification determines their grouping. While other aggregation methods distribute the areas into classes, in an artificial manner, by imposing a certain probability for an area to belong to a certain class, as determined by the assumption that the aggregation measure used is normally distributed, CA does not constrain the distribution of the areas by the classes.FloodVI was designed at the neighbourhood level and was applied to the Portuguese municipality of Vila Nova de Gaia where several flood events have taken place in the recent past. The FloodVI sensitivity was assessed using three different aggregation methods: the sum of component scores, the first component score and the weighted sum of component scores.The results highlight the sensitivity of the FloodVI to different aggregation methods. Both sum of component scores and weighted sum of component scores have shown similar results. The first component score aggregation method classifies almost all areas as having medium vulnerability and finally the results obtained using the CA show a distinct differentiation of the vulnerability where hot spots can be clearly identified.The information provided by records of previous flood events corroborate the results obtained with CA, because the inundated areas with greater damages are those that are identified as high and very high vulnerability areas by CA. This supports the fact that CA provides a reliable FloodVI.     

Multilayer shallow water flow using lattice Boltzmann method with high performance computing

A multilayer lattice Boltzmann (LB) model is introduced to solve three-dimensional wind-driven shallow water flow problems. The multilayer LB model avoids the expensive Navier–Stokes equations and obtains stratified horizontal flow velocities as vertical velocities are relatively small and the flow is still within the shallow water regime. A single relaxation time BGK method is used to solve each layer coupled by the vertical viscosity forcing term. To increase solution stability, an implicit step is suggested to obtain flow velocities. The main advantage of using the LBM is that after selecting appropriate equilibrium distribution functions, the LB algorithm is only slightly modified for each layer and retains all the simplicities of the LBM within the high performance computing (HPC) environment. The performance of the parallel LB model for the multilayer shallow water equations is investigated on CPU-based HPC environments using OpenMP. We found that the explicit loop control with cache optimization in LBM gives better performance on execution time, speedup and efficiency than the implicit loop control as the number of processors increases. Numerical examples are presented to verify the multilayer LB model against analytical solutions. We demonstrate the model’s capability of calculating lateral and vertical distributions of velocities for wind-driven circulation over non-uniform bathymetry.     

Upscaling of field-scale soil moisture measurements using distributed land surface modeling

Accurate coarse-scale soil moisture information is required for robust validation of current- and next-generation soil moisture products derived from spaceborne radiometers. Due to large amounts of land surface and rainfall heterogeneity, such information is difficult to obtain from existing ground-based networks of soil moisture sensors. Using ground-based field data collected during the Soil Moisture Experiment in 2002 (SMEX02), the potential for using distributed modeling predictions of the land surface as an upscaling tool for field-scale soil moisture observations is examined. Results demonstrate that distributed models are capable of accurately capturing a significant level of field-scale soil moisture heterogeneity observed during SMEX02. A simple soil moisture upscaling strategy based on the merger of ground-based observations with modeling predictions is developed and shown to be more robust during SMEX02 than upscaling approaches that utilize either field-scale ground observations or model predictions in isolation.     

高层建筑利用TLD减震的实例计算

本文进行了高层建筑利用调频液体阻尼器（TLD）减震的实例计算。计算结果表明,对高层建筑而言,同时控制结构多振型反应比只控制结构第一振型反应效果好。     

Modelling Canadian prairie wetland hydrology using a semi‐distributed streamflow model

A hydrological model (SLURP) that was designed for simulating hydrological processes taking place in large river basins was, with minimal modification, used successfully to simulate water level variations over a 28‐year period (1969–1996) for a 3‐ha prairie wetland in Saskatchewan. The model calculates a water balance based on precipitation, snowmelt, evaporation, surface runoff and subsurface flow on a daily time‐step. The model was first calibrated for two periods (1969–1973 for cropland and 1987–1990 for grassland), then it was applied to records outside the calibration periods. The model reproduced the wetland water level variations during a 28‐year period with good accuracy. The wetland water levels were most sensitive to the infiltration coefficient of surface soil under frozen conditions and to maximum soil moisture storage. The applicability of the model and the calibrated parameters to a smaller wetland, with an area of 0·24 ha, was examined. This simulation indicated that scale effects are important, probably largely in relation to snow redistribution by wind. Copyright © 2000 John Wiley & Sons, Ltd.     

Improved spatial delineation of streambed properties and water fluxes using distributed temperature sensing

A new method was developed for analysing and delineating streambed water fluxes, flow conditions and hydraulic properties using coiled fibre‐optic distributed temperature sensing or closely spaced discrete temperature sensors. This method allows for a thorough treatment of the spatial information embedded in temperature data by creating a matrix visualization of all possible sensor pairs. Application of the method to a 5‐day field dataset reveals the complexity of shallow streambed thermal regimes. To understand how velocity estimates are affected by violations of assumptions of one‐dimensional, saturated, homogeneous flow and to aid in the interpretation of field observations, the method was also applied to temperature data generated by numerical models of common field conditions: horizontal layering, presence of lateral flow and variable streambed saturation. The results show that each condition creates a distinct signature visible in the triangular matrices. The matrices are used to perform a comparison of the behaviour of one‐dimensional analytical heat‐tracing models. The results show that the amplitude ratio‐based method of velocity calculation leads to the most reliable estimates. The minimum sensor spacing required to obtain reliable velocity estimates with discrete sensors is also investigated using field data. The developed method will aid future heat‐tracing studies by providing a technique for visualizing and comparing results from fibre‐optic distributed temperature sensing installations and testing the robustness of analytical heat‐tracing models. Copyright © 2016 John Wiley & Sons, Ltd.     

Evolutionary assimilation of streamflow in distributed hydrologic modeling using in-situ soil moisture data

This study has applied evolutionary algorithm to address the data assimilation problem in a distributed hydrological model. The evolutionary data assimilation (EDA) method uses multi-objective evolutionary strategy to continuously evolve ensemble of model states and parameter sets where it adaptively determines the model error and the penalty function for different assimilation time steps. The assimilation was determined by applying the penalty function to merge background information (i.e., model forecast) with perturbed observation data. The assimilation was based on updated estimates of the model state and its parameterizations, and was complemented by a continuous evolution of competitive solutions.The EDA was illustrated in an integrated assimilation approach to estimate model state using soil moisture, which in turn was incorporated into the soil and water assessment tool (SWAT) to assimilate streamflow. Soil moisture was independently assimilated to allow estimation of its model error, where the estimated model state was integrated into SWAT to determine background streamflow information before they are merged with perturbed observation data. Application of the EDA in Spencer Creek watershed in southern Ontario, Canada generates a time series of soil moisture and streamflow. Evaluation of soil moisture and streamflow assimilation results demonstrates the capability of the EDA to simultaneously estimate model state and parameterizations for real-time forecasting operations. The results show improvement in both streamflow and soil moisture estimates when compared to open-loop simulation, and a close matching between the background and the assimilation illustrates the forecasting performance of the EDA approach.