An improved HASM method for dealing with large spatial data sets

Surface modeling with very large data sets is challenging. An efficient method for modeling massive data sets using the high accuracy surface modeling method (HASM) is proposed, and HASM_Big is developed to handle very large data sets. A large data set is defined here as a large spatial domain with high resolution leading to a linear equation with matrix dimensions of hundreds of thousands. An augmented system approach is employed to solve the equality-constrained least squares problem (LSE) produced in HASM_Big, and a block row action method is applied to solve the corresponding very large matrix equations. A matrix partitioning method is used to avoid information redundancy among each block and thereby accelerate the model. Experiments including numerical tests and real-world applications are used to compare the performances of HASM_Big with its previous version, HASM. Results show that the memory storage and computing speed of HASM_Big are better than those of HASM. It is found that the computational cost of HASM_Big is linearly scalable, even with massive data sets. In conclusion, HASM_Big provides a powerful tool for surface modeling, especially when there are millions or more computing grid cells.     

气候模式模拟中海量数据的并行抽取与可视化分析诊断

在地球动力学和气候模拟等领域, 数值模拟产生的数据规模达到Tb至Pb量级。 实现这些海量数据的实时可视化和实时诊断分析面临很大的困难, 时空多尺度数据抽取可以解决这一瓶颈。 高精度数据的可视化结果展示需要高分辨率的显示设备, 并行的大屏幕显示技术是解决这一问题的有效手段。 地球科学数值模拟、 并行数据抽取和高分辨率显示都需要搭建高性能计算机集群。 本文在基于Lagrange插值的多维度、 多尺度、 多分辨率并行数据抽取算法的同时, 利用并行计算节点及LCD显示器, 基于Rocks cluster系统搭建起一个176核, 4×10×1024×1280分辨率的高性能计算模拟、 数据抽取和并行显示输出的集成平台, 并将该平台成功应用到气候模式模拟产生的海量数据的并行抽取和并行显示。     

基于移动大数据的地震人口热力实时获取系统实现与应用——以泸定MS6.8地震为例

震后第一时间获取震中的人口情况对于开展应急处置决策非常关键。获取海量移动互联网人口分布数据,采用分布式集群进行海量数据的实时处理和网格化,为确保震后系统能够达到秒级产出的效率,设计了多级的多策略抽稀计算并基于HeatMap实现可视化展示,震后实现秒级震中周边指定范围的热力人口数据计算和分布图的自动绘制。以泸定 M_S6.8 地震为例进行应用,体现出基于海量移动互联网用户数据获取震中人口分布的优势。     

基于主成分分析方法的海量地震数据属性降维优化

针对传统的地震数据属性降维优化方法所选取的地震数据属性特征贡献率低导致降维过程计算量大、CPU占用率高等问题,提出一种基于主成分分析的海量地震数据属性降维优化方法。首先根据地震样本特征建立地震数据特征矩阵,把矩阵中的特征进行聚类,运用降序法排列聚类结果,选取前几项数据作为地震数据属性特征选取结果,对其结果评估分类信息量;通过特征积分准则(FSC)修正分类信息量,获取海量地震数据属性特征节点;运用主成分分析方法对地震数据属性特征节点主成分添加标签,确定Fisher判别分析与PCA可变动选择不确定关系,建立半监督降维的全局最优化形式,运用特征值分解计算降维结果,克服海量地震数据属性降维过程中的过拟合问题,融合主成分分析算法与Fisher判别分析算法实现海量地震数据属性降维优化。实验结果证明,所提方法选取的属性特征精度及贡献率较高,降维过程中CPU占用率较低。     

Identification of accelerometric stations in ITACA with distinctive features in their seismic response

In this work, we analyze the records of the Italian strong motion database (ITACA, ) with the aim of identifying stations affected by site effects that are not captured by standard seismic classification schemes. In particular, we consider four different site classifications, two of them based on geological/geophysical characteristics and two driven by data. For each classification we develop a ground motion prediction equation using a random effect approach to isolate the between-station and within-station distribution of errors. The site coefficients obtained for the different classes confirm that site amplification effects are significant for both the horizontal and vertical components. The between-station error normalized to the standard deviation of the between-station error distribution is then used to identify stations characterized by large errors, attributable to site effects not accounted for by the classification schemes. The results show that large errors can affect the predictions when the site effects are not uniquely related to the reduction of the seismic impedance in the uppermost layers. For example, amplifications of ground motion over the long period range are observed for stations installed within alluvial closed-shape basins, as consequence of locally generated surface waves. For these stations, classifications based on the horizontal to vertical response spectra ratio are not reliable, since amplifications are also affecting the vertical component. Another interesting feature which emerges from the analysis is the significant de-amplification of short period spectral ordinates that seems to be related to stations typically set in at the foundation level of massive structures. To increase the usefulness of the data set, the most important distinctive features of the strong motion stations are documented in the ITACA database reports containing the instrument information and the available geological-geotechnical data.     

大型数据平台搭建技术在地震行业中的应用与研究

依据地震行业业务系统的数据处理和存储需求,提出了数据模型和数据平台建设思路。结合地震行业网中数据平台的建设,在国家中心、区域中心和台站节点分别采用了相应的数据库平台和数据存储平台技术,设计了各类节点的数据备份策略,实现了满足海量数据流高效运行的大型数据平台。     

Merging active and passive data sets in traveltime tomography: The case study of Campi Flegrei caldera (Southern Italy)

We propose a strategy for merging both active and passive data sets in linearized tomographic inversion. We illustrate this in the reconstruction of 3D images of a complex volcanic structure, the Campi Flegrei caldera, located in the vicinity of the city of Naples, southern Italy. The caldera is occasionally the site of significant unrests characterized by large ground uplifts and seismicity. The P and S velocity models of the caldera structure are obtained by a tomographic inversion based on travel times recorded during two distinct experiments. The first data set is composed of 606 earthquakes recorded in 1984 and the second set is composed of recordings for 1528 shots produced during the SERAPIS experiment in 2001. The tomographic inversion is performed using an improved method based on an accurate finite‐difference traveltime computation and a simultaneous inversion of both velocity models and earthquake locations. In order to determine the adequate inversion parameters and relative data weighting factors, we perform massive synthetic simulations allowing one to merge the two types of data optimally. The proper merging provides high resolution velocity models, which allow one to reliably retrieve velocity anomalies over a large part of the tomography area. The obtained images confirm the presence of a high P velocity ring in the southern part of the bay of Pozzuoli and extends its trace inland as compared to previous results. This annular anomaly represents the buried trace of the rim of the Campi Flegrei caldera. Its shape at 1.5 km depth is in good agreement with the location of hydrothermalized lava inferred by gravimetric data modelling. The Vp/Vs model confirms the presence of two characteristic features. At about 1 km depth a very high Vp/Vs anomaly is observed below the town of Pozzuoli and is interpreted as due to the presence of rocks that contain fluids in the liquid phase. A low Vp/Vs body extending at about 3–4 km depth below a large part of the caldera is interpreted as the top of formations that are enriched in gas under supercritical conditions.     

Potential of legacy 2D seismic data for deep targeting and structural imaging at the Neves–Corvo massive sulphide-bearing deposit,Portugal

Seismic methods are becoming an established choice for deep mineral exploration after being extensively tested and employed for the past two decades. To investigate whether the early European mineral-exploration datasets had potential for seismic imaging that was overlooked, we recovered a low-fold legacy seismic dataset from the Neves–Corvo mine site in the Iberian Pyrite Belt in southern Portugal. This dataset comprises six 4–6 km long profiles acquired in 1996 for deep targeting. Using today's industry-scale processing algorithms, the world-class, ca. 150 Mt, Lombador massive sulphide and other smaller deposits were better imaged. Additionally, we also reveal a number of shallow but steeply dipping reflections that were absent in the original processing results. This study highlights that legacy seismic data are valuable and should be revisited regularly to take advantage of new processing algorithms and the experiences gained from processing such data in hard-rock environments elsewhere. Remembering that an initial processing job in hard rock should always aim to first obtain an overall image of the subsurface and make reflections visible, and then subsequent goals of the workflow could be set to, for example understanding relative amplitude ratios. The imaging of the known mineralization implies that this survey could likely have been among one of the pioneer studies in the world that demonstrated the capability of directly imaging massive sulphide deposits using the seismic method.     

连续重力观测潮汐与非潮汐信号提取与处理平台设计与实现

介绍了连续重力观测潮汐与非潮汐信号提取与处理平台的平台架构、数据库、主要功能模块、技术特点与功能特性等。该平台面向中国地震重力站网的大网运维业务需求和海量数据处理需求,以连续重力数据处理中的潮汐信号处理、非潮汐信号处理、时频变换等科学计算为核心,基于大数据技术实现了中国地震重力站网海量数据的自动清理、数据质量自动评估、数据产品自动产出,产出重力潮汐因子、时频图等24种数据产品,提升了连续重力观测数据的应用效能和地震重力站网的运行、管理和服务水平。     

佳木斯地块及东缘岩石圈电性结构特征

佳木斯地块及东缘是中国东北地区的重要地质构造单元. 为探测该地区地壳深部结构与构造关系,沿桦南—饶河实施了240 km的大地电磁测深的探测研究. 采用光滑模型二维反演方法对桦南—饶河大地电磁剖面的探测数据进行了二维反演和综合地球物理解释. 研究结果揭示了研究区断面域的电性结构：(1)剖面西段具有稳定的高阻特征且具有稳定的岩石圈厚度（约90km）,在十几公里深度范围存在壳内高导层,该区段对应佳木斯地块. (2)剖面中部具有明显的电性梯度带,该梯度带反映了佳木斯地块的东界位置及其深部的结构形态. (3)剖面东段电性特征揭示了佳木斯地块以东区段浅部为逆冲推覆体,深部为多个高阻块体和低阻条带相间的电性结构,这些高阻块体可能为早期俯冲的岩石圈残片.     

Coral bleaching and mortality on artificial and natural reefs in Maldives in 1998, sea surface temperature anomalies and initial recovery

The bleaching and subsequent mortality of branching and massive corals on artificial and natural reefs in the central atolls of Maldives in 1998 are examined with respect to sea surface temperature (SST) anomalies. SST normally peaks in April-May in Maldives. The UK Meteorological Office's Global sea-Ice and SST data set version 2.3 b shows that in 1998 monthly mean SST was 1.2-4 S.D. above the 1950-1999 average during the warmest months (March-June), with the greatest anomaly in May of +2.1 degrees C. Bleaching was first reported in mid-April and was severe from late April to mid-May with some recovery evident by late-May. At least 98% of branching corals (Acroporidae, Pocilloporidae) on artificial structures deployed on a reef flat in 1990 died whereas the majority of massive corals (Poritidae, Faviidae, Agariciidae) survived the bleaching. The pre-bleaching coral community on the artificial reefs in 1994 was 95% branching corals and 5% massives (n = 1589); the post-bleaching community was 3% branching corals and 97% massives (n = 248). Significant reductions in live coral cover were seen at all natural reefs surveyed in the central atolls, with average live coral cover decreasing from about 42% to 2%, a 20-fold reduction from pre-bleaching levels. A survey of recruitment of juvenile corals to the artificial structures 10 months after the bleaching event showed that 67% of recruits (> or = 0.5 cm diameter) were acroporids and pocilloporids and 33% were from massive families (n = 202) compared to 94% and 6%, respectively, in 1990-1994 (n = 3136). Similar post-bleaching dominance of recruitment by branching corals was seen on nearby natural reef (78% acroporids and pocilloporids; 22% massives). A linear regression of April mean monthly SST against year was highly significant (p < 0.001) and suggests a rise of 0.16 degree C per decade. If this trend continues, by 2030 mean April SST in the central atolls will normally exceed the anomaly level at which corals appear there are susceptible to mass bleaching.     

三维立体可视化震情会商系统平台设计与实现

随着政府和社会公众对建立现代震情会商系统的需求日益增强,本文尝试应用三维GIS技术和体绘制可视化方法,解决地震会商过程中的多源、海量地学数据的协同可视化问题.设计并实现了三维立体可视化会商系统平台,并在其上将汶川余震分布、数值模拟和华北精细地下三维速度结构等研究成果进行了集成测试.     

时序数据库在海量地震波形数据分布式存储与处理中的应用初探

通过分析时序数据库的技术特点和海量地震数据存储和处理的需求,设计和实现了基于时序数据库的地震波形数据存储和处理系统。经过开发、测试和应用比较,新的地震时序数据库系统具有高可用、实时存储吞吐量大、历史数据压缩率高、分布式处理速度快等优势,能够较好地解决地震波形数据体量激增带来的存储和处理困难问题。     

前兆观测异常数据检测方法研究

借鉴数据挖掘中的算法,本文设计了一种前兆时序模式表示方法,利用该方法可以快速检测数据序列中大幅突跳、阶跃等比较明显的异常数据。实际观测数据应用结果表明,该方法对于大量数据的异常检测效率很高,对前兆数据的预处理工作具有积极意义。     

Virtual ChuanDian——A parallel numerical modeling of Sichuan-Yunnan regional strong earthquake activities: Model construction and parallel simulation

Reasonable use of mathematical modeling and numerical simulation of regional earthquake activities has attracted more and more attentions in realistic seismicity studies in recent years. The inherent nature of their attendant multi-scale properties, from either the temporal or the spatial direction, necessitates the use of the cutting-edge massive parallel data processing techniques. Our study here is focused on the model construction and numerical simulation of the regional stress evolution after large ear...     

基于树形结构的地震数据库系统设计方法

地震观测数据的大量增长给数据管理提出了新的要求. 本文将树形结构引入到数据库设计中, 提出了基于树形结构设计关系数据库的方法. 该方法简单、 直观、 易于数据的组织, 可以简化数据库的设计过程. 在ldquo;长江三峡工程诱发地震监测系统地震监测总站数据库rdquo;设计中应用了该方法, 并取得了良好的效果.      

The method of design for seismic data database system based on tree structure

Introduction In recent years, with the development of earthquake observation technology, increasing ofobservation stations and improvement of observation instruments which are digitalized and con-nected by network, a great deal of data is recorded, to some extend, which brings difficulty for usto manage and store them. It leads to traditional methods, which use the file system to managemass data, could not satisfy our requirement, and it is necessary to find a new method in place oftradit…     

滇西区域水化观测台网近源场震兆特征

在本文中，作者用变幅滑动平均、滑动残差、群体突跳率３种方法，对滇西地区水化观测台网建网以来的资料进行了系统的综合处理，并与邻近地区１０多次５级以上地震进行相关分析，结果表明，在台网内发生的地震之前，多数观测点及观测项目都有背景异常显示，异常持续时间长短及幅度大小与地震强度、震中距呈正相关。短期异常表现为群体突跳，当群体突跳率达峰值并发生转折后４０－５０天地震发生。该台网有较强的映震能力，各观测点的     

Deep massive sulphide exploration using 2D and 3D geoelectrical and induced polarization data in Skellefte mining district,northern Sweden

Geoelectrical and induced polarization data from measurements along three profiles and from one 3D survey are acquired and processed in the central Skellefte District, northern Sweden. The data were collected during two field campaigns in 2009 and 2010 in order to delineate the structures related to volcanogenic massive sulphide deposits and to model lithological contacts down to a maximum depth of 1.5 km. The 2009 data were inverted previously, and their joint interpretation with potential field data indicated several anomalous zones. The 2010 data not only provide additional information from greater depths compared with the 2009 data but also cover a larger surface area. Several high‐chargeability low‐resistivity zones, interpreted as possible massive sulphide mineralization and associated hydrothermal alteration, are revealed. The 3D survey data provide a detailed high‐resolution image of the top ～450 m of the upper crust around the Maurliden East, North, and Central deposits. Several anomalies are interpreted as new potential prospects in the Maurliden area, which are mainly concentrated in the central conductive zone. In addition, the contact relationship between the major geological units, e.g., the contact between the Skellefte Group and the Jörn Intrusive Complex, is better understood with the help of 2010 deep‐resistivity/chargeability data. The bottommost part of the Vargfors basin is imaged using the 2010 geoelectrical and induced polarization data down to ～1‐km depth.