共查询到20条相似文献，搜索用时 31 毫秒
1.
Abstract: Recent rapid progress in cyberinfrastructure in geosciences is providing seismologists an enormous boost for addressing multi-physical phenomena of regional seismic activities. The inherent nature of their multi-scale properties, from temporal to spatial spaces, makes it inevitably to be solved using large-scale computations and distributed parallel data processing schemes. Under such circumstance, using the advanced numerical algorithms and unstructured mesh generation technologies become the obstacles for modern seismologists. The main objective of this paper is to present a framework, which includes a parallel finite element simulation and distributed data infrastructure, to address the novel algorithms, state-of-the-art modeling and their implementation in regional seismicgenic systems. We also discuss and implement this framework to analyze the strong earthquake evolution processes in the Sichuan-Yunnan region. This study is the key to long-term seismic risk by estimates, providing a platform for predictive large-scale numerical simulation modeling of regional earthquake activities. 相似文献
2.
Huai Zhang Shaolin Chen Shi Chen Si Chen Huimin Jing David A. Yuen Yaolin Shi 《Visual Geosciences》2008,13(1):117-124
Today parallel visualization of massive datasets from observation and numerical simulation of seismic waves is one of the
major goals of geoscience community. A majority of these datasets are time-varying volume data (TVVD), also known as 4D field
data. The difficulty of visualizing them on distributed parallel system mainly lies in the algorithm designing for distributed
preprocessing of raw datasets, hierarchical point-to-point or collective communication implementation based on distributed
data allocation, synchronous volume rendering techniques. In this work we present viable solutions for preprocessing of raw
data sets, novel algorithms of parallel rendering and display matrix. Our main objective is focused on the parallel visualization
of results coming from full 4D seismic wave propagation simulations. 相似文献
3.
Yifeng Cui Reagan Moore Kim Olsen Amit Chourasia Philip Maechling Bernard Minster Steven Day Yuanfang Hu Jing Zhu Thomas Jordan 《Acta Geotechnica》2009,4(2):79-93
Earthquakes are among the most complex terrestrial phenomena, and modeling of earthquake dynamics is one of the most challenging
computational problems in science. Computational capabilities have advanced to a state where we can perform wavefield simulations
for realistic three-dimensional earth models, and gain more insights into the earthquakes that threaten California and many
areas of the world. The Southern California Earthquake Center initiated a major earthquake research program called TeraShake
to perform physics-based numerical simulations of earthquake processes for large geographical regions, at high resolution,
and for high frequencies. For a large scale simulation such as TeraShake, optimization problems tend to emerge that are not
significant in smaller scale simulations. This involves both large parallel computation and also massive data management and
visualization coordination. In this paper, we describe how we performed single-processor optimization of the TeraShake AWM
application, optimization of the I/O handling, and optimization of initialization. We also look at the challenges presented
by run-time data archive management and visualization. The improvements made to the TeraShake AWM code enabled execution on
the 40k IBM Blue Gene processors and have created a community code that can be used by seismologists to perform petascale
earthquake simulations. 相似文献
4.
Reza Tavakoli Gergina Pencheva Mary F. Wheeler Benjamin Ganis 《Computational Geosciences》2013,17(1):83-97
We present a parallel framework for history matching and uncertainty characterization based on the Kalman filter update equation for the application of reservoir simulation. The main advantages of ensemble-based data assimilation methods are that they can handle large-scale numerical models with a high degree of nonlinearity and large amount of data, making them perfectly suited for coupling with a reservoir simulator. However, the sequential implementation is computationally expensive as the methods require relatively high number of reservoir simulation runs. Therefore, the main focus of this work is to develop a parallel data assimilation framework with minimum changes into the reservoir simulator source code. In this framework, multiple concurrent realizations are computed on several partitions of a parallel machine. These realizations are further subdivided among different processors, and communication is performed at data assimilation times. Although this parallel framework is general and can be used for different ensemble techniques, we discuss the methodology and compare results of two algorithms, the ensemble Kalman filter (EnKF) and the ensemble smoother (ES). Computational results show that the absolute runtime is greatly reduced using a parallel implementation versus a serial one. In particular, a parallel efficiency of about 35 % is obtained for the EnKF, and an efficiency of more than 50 % is obtained for the ES. 相似文献
5.
为了解决长时间正演模拟的不稳定性问题, 实现天然地震波场长时间数值模拟.在高阶有限差分数值模拟的基础上, 给出了多轴向完全匹配吸收边界(M-PML)二维划分方式, 讨论了M-PML吸收边界转换系数P的取值对其吸收能力与截断误差对数值模拟有效信号的影响.实现了二维弹性波高阶有限差分的长时间数值模拟.采用2013年12月16日湖北省巴东县地震及余震信息, 对该地震数据进行了数值模拟, 并与武汉、秭归两台站数据进行了分析和讨论, 验证该正演模拟方法的正确性和有效性.为天然地震波场传播规律、成像和震相识别及震源定位等研究提供了基础. 相似文献
6.
在高精度三维地震资料基础上，通过信号强化、大比例显示及蚂蚁体追踪技术，在剖面上建立了小断层的识别标志：1）同相轴的轻微错断；2）同相轴的分叉、合并；3）同相轴的挠曲、强相位转换；4）反射零乱或出现空白带；5）同相轴突然增减或消失，波组间隔突然变化，表现为不连续；6）借助蚂蚁体剖面识别小断层.通过分析临北地区不同规模断层的空间组合样式，利用平衡剖面原理，恢复了不同规模断层的演化过程.构造物理模拟显示，临北地区小断层主要为主干断层活动的派生或伴生构造，大断层附近的区域应力场是其主要形成原因，形成的小断层多与大断裂平行，呈平行式分布. 相似文献
7.
在高精度三维地震资料基础上,通过信号强化、大比例显示及蚂蚁体追踪技术,在剖面上建立了小断层的识别标志：1）同相轴的轻微错断;2）同相轴的分叉、合并;3）同相轴的挠曲、强相位转换;4）反射零乱或出现空白带;5）同相轴突然增减或消失,波组间隔突然变化,表现为不连续;6）借助蚂蚁体剖面识别小断层.通过分析临北地区不同规模断层的空间组合样式,利用平衡剖面原理,恢复了不同规模断层的演化过程.构造物理模拟显示,临北地区小断层主要为主干断层活动的派生或伴生构造,大断层附近的区域应力场是其主要形成原因,形成的小断层多与大断裂平行,呈平行式分布. 相似文献
8.
Quantification of seismic activity is one of the most challenging problems faced by earthquake engineers in probabilistic seismic hazard analysis. Currently, this problem has been attempted using empirical approaches which are based on the regional earthquake recurrence relations from the available earthquake catalogue. However, at a specified site of engineering interest, these empirical models are associated with large number of uncertainties due to lack of sufficient data. Due to these uncertainties, engineers need to develop mechanistic models to quantify seismic activity. A wide range of techniques for modeling continental plates provides useful insights on the mechanics of plates and their seismic activity. Among the different continental plates, the Indian plate experiences diffused seismicity. In India, although Himalaya is regarded as a plate boundary and active region, the seismicity database indicates that there are other regions in the Indian shield reporting sporadic seismic activity. It is expected that mechanistic models of Indian plate, based on finite element method, simulate stress fields that quantify the seismic potential of active regions in India. This article explores the development of a finite element model for Indian plate by observing the simulated stress field for various boundary conditions, geological and rheological conditions. The study observes that the magnitude and direction of stresses in the plate is sensitive to these conditions. The numerical analysis of the models shows that the simulated stress field represents the active seismic zones in India. 相似文献
9.
Nicole Bethoux Monica Segovia Viviana Alvarez Jean-Yves Collot Philippe Charvis Audrey Gailler Tony Monfret 《Journal of South American Earth Sciences》2011,31(1):139-152
A seismic study of a segment of the convergent margin of Ecuador is presented. During the SISTEUR campaign a network of 24 Ocean Bottom Seismometers (OBS) was deployed on the Carnegie Ridge, one line along the main axes of the ridge and two lines across the strike of the edge of the ridge, during one month. This marine network was complemented with a land network of 20 stations distributed in two lines: one parallel to the margin and the other perpendicular to it.The seismic event recorded by these networks, were located using different crustal models defined from the wide-angle seismic data modeling. Relative location techniques were used to improve earthquake locations. Seismogram waveform modeling allowed us to constrain hypocentral location for events farther than ～50 km from the network. This modeling also provided additional information to constrain the focal mechanisms of these events. The upper limit of the Interplate Seismogenic Zone (ISZ) is estimated to be at a 10 km depth in the region. The background seismic activity of the upper plate provided new insights:1) A seismic cluster that reaches the base of the overriding plate is linked to the Jipijapa-Portoviejo fault. The reactivation of this Quaternary fault is confirmed by focal mechanisms that provide rupture planes parallel to its superficial projection (N10°–N25°).2) The focal mechanisms presented in this study are compatible with a homogeneous regional stress field corresponding to an E–W to ESE–WNW compression and an NNE–SSW extension. The presence of strike-slip deformation, with a reverse component, corresponds to the NNE escape of the North Andean Block. Normal faulting accommodating this movement suggests that this part of the North Andean Block cannot be considered as a rigid block. 相似文献
10.
Dave A.Yuen Benjamin J. Kadlec Evan F. Bollig Witold Dzwinel Zachary A. Garbow Cesar R. S. da Silva 《Visual Geosciences》2005,10(1):1-12
We present a web client-server service WEB-IS, which we have developed for remote analysis and visualization of seismic data consisting of both small magnitude events and large earthquakes. We show that the problem-solving environment (PSE) intended for prediction of large magnitude earthquakes can be based on this WEB-IS idea. The clustering schemes, feature generation, feature extraction techniques and rendering algorithms form a computational framework of this environment. On the other hand, easy and fast access both to the seismic data distributed among distant computing resources and to computational and visualization resources can be realized in a GRID framework. We discuss the usefulness of NaradaBrokering (iNtegrated Asynchronous Real-time Adaptive Distributed Architecture) as a middleware, allowing for flexibility and high throughput for remote visualization of geophysical data. The WEB-IS functionality was tested both on synthetic and the actual earthquake catalogs. We consider the application of similar methodology for tsunami alerts. 相似文献
11.
A Numerical Framework for Modeling Folds in Structural Geology 总被引：1，自引：1，他引：0
A numerical framework for modeling folds in structural geology is presented. This framework is based on a novel and recently published Hamilton–Jacobi formulation by which a continuum of layer boundaries of a fold is modeled as a propagating front. All the fold classes from the classical literature (parallel folds, similar folds, and other fold types with convergent and divergent dip isogons) are modeled in two and three dimensions as continua defined on a finite difference grid. The propagating front describing the fold geometry is governed by a static Hamilton–Jacobi equation, which is discretized by upwind finite differences and a dynamic stencil construction. This forms the basis of numerical solution by finite difference solvers such as fast marching and fast sweeping methods. A new robust and accurate scheme for initialization of finite difference solvers for the static Hamilton–Jacobi equation is also derived. The framework has been integrated in simulation software, and a numerical example is presented based on seismic data collected from the Karama Block in the North Makassar Strait outside Sulawesi. 相似文献
12.
13.
历史强震对渭河中游群发大型滑坡的诱发效应反演 总被引：1，自引：0，他引：1
以渭河中游地区为例,探索提出了开展历史地震对区域群发滑坡诱发效应反演研究的思路和方法。首先,基于汶川地震在渭河中游地区形成的高烈度异常和震害启示,通过区域活动构造和斜坡带断裂控滑分析,指出历史强震对区内群发大型滑坡的诱发效应不容忽视。然后,利用强震诱发滑坡的最远致灾震中距分析法,筛选出研究区周边300 km范围内需要重点考察其诱发效应的4次关键历史强震:公元前780年岐山MS7.0级地震、1654年天水南MS8.0级地震、1556年华县MS8.25级地震及1920年海原MS8.5级地震。随后,以岐山地震为例,具体阐述了基于Newmark位移模型的地震诱发滑坡位移及危险性反演评估方法;同时反演了其他3次历史强震诱发区内滑坡位移及危险性。最后,定量比较了反演历史强震诱发滑坡的位移与实际大型滑坡分布的空间匹配程度,结果显示天水南MS8.0级地震对渭河中游现存群发大型滑坡的诱发效应最强。 相似文献
14.
地震滑坡危险性概念和基于力学模型的评估方法探讨 总被引：1，自引：0，他引：1
在我国大陆地区运用基于力学模型的Newmark位移分析方法开展地震滑坡危险性定量评估,尚处在起步阶段.为了进一步明确地震滑坡危险性概念和改进推广基于力学模型的评估方法,首先阐明了狭义的地震滑坡危险性预测评估与震后反演评估的关系; 同时为了应对地震应急、震后重建及潜在地震条件下的不同评估需求,初步提出了广义的地震滑坡危险性评估框架.随后申述了基于Newmark位移分析的地震滑坡危险性评估方法的理论基础、方法分类及最新进展,并以汶川地震滑坡危险性快速评估为例,剖析了目前影响评估有效性的不确定性及空间数据质量等问题,指出了基于力学模型的地震滑坡危险性评估方法的改进方向.建议开展潜在地震及其诱发滑坡危险性的耦合评估,建立适用于我国大陆地区地震滑坡位移分析的经验模型,以便为国家层面的地震滑坡危险性区划服务. 相似文献
15.
Assessment of seismic vulnerability of urban areas provides fundamental information for activities of planning and management of emergencies. The main difficulty encountered when extending vulnerability evaluations to urban contexts is the definition of a framework of assessment appropriate for the specific characteristics of the site and providing reliable results with a reasonable duration of surveys and post-processing of data. The paper proposes a new procedure merging different typologies of information recognized on the territories investigated and for this reason called “hybrid.” Knowledge of historical events influencing urban evolution and analysis of recurrent building technologies are used to evaluate the vulnerability indexes of buildings and building stocks. On the other hand, a vulnerability model is calibrated by means of experimental and numerical investigations on prototype buildings representative of the most recurrent typologies. In the final framework, the vulnerability index, calculated through simplified assessment forms, is linked to the seismic intensity expressed by the peak ground acceleration and associated with an index of damage expressing the economical loss. The procedure has been tested on the urban center of Lampedusa island (Italy) providing as the output vulnerability index maps, vulnerability curves, critical PGA maps, and estimation of the economical damage associated with different earthquake scenarios. The application of the procedure can be suitably repeated for medium-to-small urban areas, typically recurring in the Mediterranean by carrying out each time a recalibration of the vulnerability model.
相似文献16.
台湾大地构造环境与地震 总被引：1，自引：0，他引：1
台湾、吕宋和琉球等岛弧区是中新世以来地震多发区.1999年9月21日台湾南投大地震以其震级大、震期长和损失惨重引起国内外关注.本文以台湾所处大地构造环境特征探讨发震机制和9.21大地震对地学界的警示. 相似文献
17.
多波地震勘探技术能够同时获得纵波和转换波资料,与常规的纵波勘探相比,能够提供更多的地下介质信息。且转换波对于埋深较浅的小断层、小幅度构造有更高的分辨率,充分利用多分量地震资料可以有效的提高地震勘探的精度。本文将多波地震技术引入到陆域冻土区天然气水合物勘探中,对陆域冻土区水合物进行多波地震数值模拟研究。采用弹性波有限差分方法进行多分量正演模拟,利用仿射坐标转换法将原始的多分量地震数据进行纵横波叠前分离,并通过多波高斯束偏移方法分别对波场分离后的PP波和转换PS波进行叠前深度偏移成像。根据我国祁连山冻土区水合物实际地震地质资料,建立逼近真实地质情况的天然气水合物数值模型。对上述水合物模型进行多波地震数值试算表明,多波地震勘探技术是一种有效的陆域冻土区水合物探测方法,充分利用多波地震资料有利于查明陆域冻土区水合物精细的地质构造,取得更好的勘探效果。 相似文献
18.
基于近场地震记录频谱分析,结合地震危险性分析得到的基岩场地反应谱,进行了高山峡谷地区大型地下洞室群的人工地震动加速度时程的合成。该方法使用三角级数迭加法,利用攀枝花地震动记录进行傅里叶变换,得到离散点的傅里叶幅值谱,选择两个卓越频率作为控制傅里叶谱的峰值控制点,通过反复迭代生成了初始人工地震波,根据高山峡谷地形进行了深度修正、基线校正和高频滤波,最终生成了可以用于FLAC和ABAQUS等数值分析软件的人工地震动加速度-时程。数值分析表明,此方法拟合的人工地震波可以较好地反映不同频谱对地下洞室群强震响应的影响。 相似文献
19.
20.
Assessment of seismic vulnerability of urban areas provides fundamental information for activities of planning and management of emergencies. The main difficulty encountered when extending vulnerability evaluations to urban contexts is the definition of a framework of assessment appropriate for the specific characteristics of the site and providing reliable results with a reasonable duration of surveys and post-processing of data. The paper proposes a new procedure merging different typologies of information recognized on the territories investigated and for this reason called “hybrid.” Knowledge of historical events influencing urban evolution and analysis of recurrent building technologies are used to evaluate the vulnerability indexes of buildings and building stocks. On the other hand, a vulnerability model is calibrated by means of experimental and numerical investigations on prototype buildings representative of the most recurrent typologies. In the final framework, the vulnerability index, calculated through simplified assessment forms, is linked to the seismic intensity expressed by the peak ground acceleration and associated with an index of damage expressing the economical loss. The procedure has been tested on the urban center of Lampedusa island (Italy) providing as the output vulnerability index maps, vulnerability curves, critical PGA maps, and estimation of the economical damage associated with different earthquake scenarios. The application of the procedure can be suitably repeated for medium-to-small urban areas, typically recurring in the Mediterranean by carrying out each time a recalibration of the vulnerability model. 相似文献