甘肃省地震局高性能计算系统

高性能计算是现代科学研究、工程技术开发和大规模数据处理的关键支撑技术,在地震系统发挥着日益重要的作用。甘肃省地震局高性能计算集群系统由60个节点、采用40 GB InfiniBand网络互连组成,峰值计算速度达9.2万亿次/秒。作为地震系统重要的科技服务平台之一,为地震反演、数据分析、数值模拟等领域以及科研创新平台提供高性能计算环境。     

Three-dimensional forward modeling of DC resistivity using the aggregation-based algebraic multigrid method

To speed up three-dimensional (3D) DC resistivity modeling, we present a new multigrid method, the aggregation-based algebraic multigrid method (AGMG). We first discretize the differential equation of the secondary potential field with mixed boundary conditions by using a seven-point finite-difference method to obtain a large sparse system of linear equations. Then, we introduce the theory behind the pairwise aggregation algorithms for AGMG and use the conjugate-gradient method with the V-cycle AGMG preconditioner (AGMG-CG) to solve the linear equations. We use typical geoelectrical models to test the proposed AGMG-CG method and compare the results with analytical solutions and the 3DDCXH algorithm for 3D DC modeling (3DDCXH). In addition, we apply the AGMG-CG method to different grid sizes and geoelectrical models and compare it to different iterative methods, such as ILU-BICGSTAB, ILU-GCR, and SSOR-CG. The AGMG-CG method yields nearly linearly decreasing errors, whereas the number of iterations increases slowly with increasing grid size. The AGMG-CG method is precise and converges fast, and thus can improve the computational efficiency in forward modeling of three-dimensional DC resistivity.     

A method of reconstructing complex stratigraphic surfaces with multitype fault constraints

The construction of complex stratigraphic surfaces is widely employed in many fields, such as petroleum exploration, geological modeling, and geological structure analysis. It also serves as an important foundation for data visualization and visual analysis in these fields. The existing surface construction methods have several deficiencies and face various difficulties, such as the presence of multitype faults and roughness of resulting surfaces. In this paper, a surface modeling method that uses geometric partial differential equations (PDEs) is introduced for the construction of stratigraphic surfaces. It effectively solves the problem of surface roughness caused by the irregularity of stratigraphic data distribution. To cope with the presence of multitype complex faults, a two-way projection algorithm between threedimensional space and a two-dimensional plane is proposed. Using this algorithm, a unified method based on geometric PDEs is developed for dealing with multitype faults. Moreover, the corresponding geometric PDE is derived, and an algorithm based on an evolutionary solution is developed. The algorithm proposed for constructing spatial surfaces with real data verifies its computational efficiency and its ability to handle irregular data distribution. In particular, it can reconstruct faulty surfaces, especially those with overthrust faults.     

Variation in reach‐scale thalweg‐migration intensity in a braided reach of the lower Yellow River in 1986–2015

Thalweg migration of an alluvial river plays a key role in channel evolution, which may influence the effect of existing river training works and biodiversity on floodplains, and cause losses in riparian land and property. The braided reach of the Lower Yellow River underwent continuous channel aggradation during the period from 1986 to 1999, and then remarkable channel degradation in 1999–2015 owing to the state of operation of the Xiaolangdi Reservoir in 1999. Here we quantify associated thalweg migration changes and identify the key influencing factor in the braided reach. Thalweg‐migration distances and intensities at section‐ and reach‐scales were calculated during the past 30 years from 1986 to 2015, in order to investigate the characteristics of thalweg migration in the reach. There was a 47% reduction in the reach‐scale thalweg‐migration distance and a 35% reduction in the corresponding migration intensity after the reservoir operation. It is also revealed that fluvial erosion intensity is a dominant factor in controlling the thalweg migration, based on the investigation into various influencing factors in the study reach. The thalweg‐migration intensity of the braided reach can be expressed as a power function of the previous four‐year average fluvial erosion intensity. The calculated thalweg‐migration intensities in 1986–2015 using the proposed relation generally agree with the observed data. Copyright © 2017 John Wiley & Sons, Ltd.     

宾川地区的地壳厚度和泊松比分布特征研究

This paper presents the changes of crust thickness and Poissons ratios distribution in the Binchuan region,where the first air-gun transmitting station and its a small dense array were deployed.From September 2011 to January 2014,more than 239 teleseismic events of M≥6.0 were recorded in 16 stations in the Binchuan region.Their P-wave receiver functions were analyzed respectively.The high spatial resolution result shows that the average crust thickness of Binchuan region is 45.3km,it follows the rule of"deeper in the north and east part,shallower in the south and west part."The deepest region is in Xiaoyindian Station;the crust thickness is 47.9km;the shallowest region is in Paiying Station,it has the thickness of 42.1km.It shows that the deeper Moho surface nearby the Chenghai fault and shallower nearby the Honghe fault;the isoline distribution of thickness changes greatly nearby the Chenghai fault and slowly nearby the Honghe fault.From the distribution of Poissons ratios,it is unevenly in the study area with a great difference from the north part to the south part,which shows a characteristic of"lower in the north,higher in the south".The Poissons ratio nearby the Honghe fault is medium too high(0.26≤σ≤0.29);lower nearby the Chenghai fault(σ≤0.26).This paper concludes the possible reason of different characteristic between Poissons ratio and crust thickness is thicker in the upper crust in the Binchuan region.     

Variation in soil saturated hydraulic conductivity along the hillslope of collapsing granite gullies

Saturated hydraulic conductivity (K_s) affects the soil hydrological process and is influenced by many factors that exhibit strong spatial variations. To accurately measure K_s and its scale, spatial variability and relationship with collapsing gullies, we analysed four double-ring infiltrometer diameters in three soil layers during in situ experiments designed to measure K_s in two typical collapsing gullies (three slope sites) in Tongcheng County of China. The results showed that K_s increased with increasing inner ring diameter, but no significant difference existed between inner diameters of 30 and 40 cm. The K_s in red soil layers was higher than that in sandy soil layers, the transition layers had the lowest values. K_s also varied with slope position, gradually decreasing with distance from the gully head. The suggestion is that the spatial variation in K_s is affected not only by the intrinsic soil properties but also by the interaction with the collapsing gully.     

蓟县地震台小辛庄山洞定点形变观测干扰

整理蓟县地震台小辛庄山洞定点形变观测资料,分析强降雨、气压变化、机井抽水、爆破、仪器故障和标定、人为干扰等因素对观测数据的影响,强化认识干扰特征和一般规律,为地震前地球物理异常判定提供准确可靠依据,为同类型地震观测台站提供借鉴。     

永清MS 4.3地震前永清井水温异常分析

2017年12月永清井观测温度急剧下降,采用该井水温、水位多年观测数据,从仪器测量原理及精细温度梯度测量出发,结合井孔资料,认为该变化为永清M_S 4.3地震（井震距<30 km）发生前的异常信息。利用正弦累加模型计算温度异常持续时间和幅度,发现该异常变化存在2个周期：①周期42天,温度变化幅度为0.042 12℃;②周期16.77天,温度变化幅度为0.018 62℃。分析认为,区域应力状态发生变化,使含水层渗透性受到影响,从而间接影响到井下温度传感器安放处水温。地震发生后近距离观测到水温异常变化尚属首次,利用模型对变化时间与幅度进行量化提取,可为经验预报、统计预报提供经典震例。     

Development of a predictive model for <Emphasis Type="Italic">Clostridium difficile</Emphasis> infection incidence in hospitals using Gaussian mixture model and Dempster–Shafer theory

Clostridium difficile infection is one of the major patient safety concerns in hospitals worldwide. Clostridium difficile infection can have high economic burden to patients, hospitals, and government. Limited work has been done in the area of predictive modeling. In this article, A new predictive model based on Gaussian mixture model and Dempster–Shafter theory is proposed to predict Clostridium difficile infection incidence in hospitals. First, the Gaussian mixture model and expectation–maximization algorithms are used to generate explicit probability criteria of risk factors based on the given data. Second, Dempster–Shafter theory is used to predict the Clostridium difficile infection incidence based on the generated probability criteria that have different beliefs attributing to their different credits. The main procedure includes (1) generate the probability criteria model using Gaussian mixture model and expectation–maximization algorithm; (2) determine the credit of the probability criteria; (3) generate the basic probability assignment; (4) discount the evidences; (5) aggregate the evidences using Dempster combining rule; (6) predict Clostridium difficile infection incidence using pignistic probability transformation. Results show that the model has a higher accuracy than an existing model. The proposed model can generate the criteria ratings of risk factors automatically, which would potentially prevent the imprecision caused by the subjective judgement of experts. The proposed model can assist risk managers and hospital administrators in the prediction and control of Clostridium difficile infection incidence with optimizing their resources.     

An enhanced bit-wise parallel algorithm for real-time GPS software receiver

Software-based global positioning system (GPS) receivers perform all the baseband signal processing and the high level functions on a general purpose processor. The heavy computational loads of the signal correlation in baseband processing make it difficult for software receivers to operate in real time. In order to improve the real-time performance, an enhanced bit-wise parallel algorithm has been developed in this study. The enhanced algorithm has been implemented and tested in a 12 channels real-time GPS software receiver. The system consists of a radio frequency front end, a data acquisition board and software that runs on a laptop with a Pentium-M 1.5 GHz processor running the Window^® XP operating system. The data acquisition board packs the 2-bit intermediate frequency samples with a 2-bit in/8-bit out shift register and transfers the packed samples to laptop through a USB port. The software running on the laptop performs all the baseband and navigation processing in real time. The test results show that the enhanced algorithm significantly improves the real-time performance of the software receiver by reducing the computational operations for signal correlation by 50% compared with the existing bit-wise parallel algorithm. Furthermore, the enhanced algorithm also reduces the amount of required memory for storing data for signal correlation.