A rapid and accurate two-point ray tracing method in horizontally layered velocity model

Introduction The determination of the ray path between a source and a receiver is fundamental to many seismic problems such as earthquake location, travel time inversion, computation of synthetic seismograms, depth migration, and seismic tomography. The calculation of travel times and their derivatives depend on the determination of ray paths, which is the most time-consuming work in seismic tomography. Therefore, a fast and accurate ray tracing method is particularly important. Many ray-traci…     

Domain decomposition based on the spectral element method for frequency-domain computational elastodynamics

We propose a domain decomposition method based on the spectral element method(DDM-SEM)for elastic wave computation in frequency domain.It combines the high accuracy of the spectral element method and the high degree of parallelism of a domain decomposition technique,which makes this method suitable for accurate and efficient simulations of large scale problems in elastodynamics.In the DDM-SEM,the original large-scale problem is divided into a number of well designed subdomains.We use the spectral element method independently for each subdomain,and the neighboring subdomains are connected by a frequency-domain version of Riemann transmission condition(RTC)for elastic waves.For the proposed method,we can employ the non-conforming meshes and different interpolation orders in different subdomains to maximize the efficiency.By separating the internal and boundary unknowns of each subdomain,an efficient and naturally parallelizable block LDU direct solver is developed to solve the final system matrix.Numerical experiments verify its accuracy and efficiency,and show that the proposed DDM-SEM can be a promising numerical tool for accurately and effectively solving large and multi-scale problems of elastic waves.It is potentially valuable for the frequency domain seismic inversion where multiple source illuminations are required.     

The Program Application on the Relative Method of Fuzzy Mathematics

The method of fuzzy mathematics for simultaneous assessment of time and intensity of earthquake hazards has been studied.This method is based on fundamental statistical indices of regional seismicity.Applying the retrieval method of fuzzy information,we can classify the time and intensity into several intervals and classes of seismic activity,then the possible time interval of large earthquakes with magnitude of M≥Ms can be estimated in a given region.Based on the preceding idea,an FRPP program is constructed.For the automatic data processing when this method is used,it is very important to design the statistical process of each index decomposition so that the program could be fit to a different sample discussed.There are some functions in the FRPP program.The man-made impact on results is reduced to the minimum as far as possible.Computation time is saved.There is a menu on which time interval,index,intensity class,and output data all can be selected.The catalog input that can be displayed on the scre     

Analysis and application of pseudo-offset method in the converted-wave prestack time migration

Pseudo-offset migration （POM） is a new method for prestack time migration of converted waves that improves on equivalent-offset migration （EOM）. The mapping of POM is different than EOM but the purpose of the two methods is to map the input samples to the common conversion scatter point （CCSP） gathers. This paper introduces the principles of the two migration methods and model parameter sensitivity tests for both POM and EOM. At large offset-to-depth ratios the hyperbolic approximation, the three-term approximation and the double square root （DSR） equation are used to NMO-correct the mapped POM gathers to obtain more accurate migration velocities. These equations were derived and calculated by small pseudo offset. POM is then used to image complex structure and prestack time migration.     

Prospects for the Development and Application of the Earthquake Monitoring Network

With the rapid development of the economy in China, the seismic network has been changing rapidly, in that the capability of instruments, technological systems and network density are approaching those of developed countries and a large quantity of observation data has been accumulated. How to apply these resources to economic construction and public safety has become an important issue worth studying. In order to improve earthquake prediction and earthquake emergency response, it is suggested in this paper that extracting valuable precursor information, improving earthquake rapid reporting ability and extending rapid intensity reporting function are key issues. Integrating network resources, building unified standards and a multifunction seismic monitoring network are preconditions of establishing a public safety service platform and earthquake observation resources will contribute significantly to the fields of engineering, ocean, meteorology, and environmental protection. Thus, the future directions of the development of the seismic network are exploring monitoring resources, enhancing independent innovation, constructing a technological platform and enlarging the service field.     

Sediment loads and processes in large rivers

<正>The pressure and impact on large rivers have increased greatly in recent years,as a consequence of their exploitation to meet various human needs.Large rivers are particularly exposed to problems of multiple uses,often with conflicting aims.At the global scale river systems are increasingly altered by dams,leading to interruption of continuity and a disrupted sediment regime,which results in far-reaching problems.These challenges point also to the need for more integrated management approaches and consideration of catchment-river interactions emphasizing land use management and effects of these uses.At the global scale,there is currently no overview assessment of the current status of the world's large rivers,the conflicting demands on such rivers,and likely future anthropogenic impacts,as well as the potential for restoration,improvements in integrated management and the associated problems caused by their multiple uses.The special issue in the International Journal of Sediment Research on"Sediment loads and processes in large rivers"aims to provide a global forum for a wide-ranging discussion of key issues related to research on large rivers and to their effective and sustainable management,involving both scientists and decision makers.The sequence of papers in the special issue highlights the current situation in different river systems with regard to morphodynamic reactions of     

The Present Status and Prospect of Earthquake Focal Depth Locating

Locating an earthquakes focal depth is always a key project in seismology. Precise focal depth is of critical importance for evaluating seismic hazards, deciphering dynamic mechanisms of earthquake generating,estimating aftershock evolutions and risk,as well as monitoring nuclear tests. However,how we determine an accurate focal depth is always a challenge in seismological studies. Aiming to solve these problems, we analyzed and summarized the present status and the future development of earthquake focal depth locating. In this paper we first reviewed the present status of focal depth locating in the world,and summarized the frequently-used relocating methods and ideas at present,and introduced two types of focal depth relocating ideas: arrival time relocating and waveform modeling methods. For these ideas,we systematically described the S-P and the Pn-Pg methods that belong to arrival time method,and polarization focal depth locating and amplitude focal depth locating that belongs to waveform modeling,and further analyzed the advantages and limitations of these methods. Since the depth phase methods are highly sensitive to focal depth,and are relatively free from the uncertainties of crustal models,we mainly reviewed the depth phases of s Pm P,s PL,s Pn,and s Sn,and quantitatively evaluated their availabilities and characteristics. Second,we also discussed the effects of crustal velocity models on the reliability of focal depth locating,and reviewed the advancements of seismic tomography techniques over recent years. Finally,based on the present status of the progress on the focal depth locating,and studies of seismic velocity structures,we proposed an idea of combining multiple datasets and relocating methods,jointly utilizing seismologic and geodetic techniques to relocate focal depth,which should be the major research field in investigating focal depth and source parameters in the near future.     

Watershed science:Bridging new advances in hydrological science with good management of river basins

<正>Watershed science traditionally refers to the themes of hydrology and water resource management.Watershed science has been experiencing a rapid evolution that thrives on a forceful superimposition of multi-discipline and innovative earth observing and information techniques.The water and its interactions with other systems in a watershed is increasingly becoming a focus in scientific communities,and several new disciplines such as ecohydrology,ecoeconomics,environmental hy-     

Variation of specific yield with depth in an alluvial aquifer of the Platte River valley, USA

Estimation of specific yield (Sy) of an aquifer is of great importance in water resource management. In this study, an experimental drainage method was developed to determine specific yield of an alluvial aquifer of the Platte River valley, Nebraska, USA. Sets of sediment cores with continuous interval depth were collected to plastic tubes using the direct push technique and then taken to the laboratory. During the Sy experiment, those sediment cores were re-saturated by placing them vertically in a large water tank. Sy was determined by the water drained from the sediments by force of gravity in a bracket. Our results show that the values of Sy varied largely with depth at each site and the variability for Sy with interval depth between the test sites is also observed. This spatial heterogeneity in Sy might result from the variation of grain size, grain shape, sorting and compaction of sediments in different cores with interval depth. The Sy for all sediment cores ranged from 0.01 to 0.18 and the mean value was 0.08±0.04. Our drainage method can functionally preserve the sedimentary structures in their original state and it is easier to experiment at a lower expense.     

中国地震网格试验系统中资源管理与作业调度研究-资源管理和作业动态分配模型ProRMJS的建立me

通过分析与研究建立资源管理和作业动态分配ProRMJS模型,解决了中国地震灾害仿真网格试验系统（CEDAGrid）中网格科学计算平台在初期建设中存在的一些资源管理方面的问题. 针对网格科学计算平台默认各个计算节点均能提供计算服务的问题,ProRMJS通过ldquo;计算池rdquo;这一虚拟组织为调度器分发作业提供支撑;而调度器则根据计算池中各个节点的计算能力和状态,动态分配任务,这样就保证了计算平台的稳定性. 针对默认各个计算节点均能完成任务的问题, ProRMJS通过监控各个节点所负担作业的运行进度和设定作业时间阈值来管理作业运行. 对于各个节点计算能力大小的区别, ProRMJS通过对各个节点计算性能的权值量化,实行按ldquo;需rdquo;分配作业. 对各个节点计算性能的权值量化, ProRMJS是通过综合考虑各个因素对机器计算能力大小的影响后通过加权计算的方式完成的,从而提高了整个计算平台的工作效率. 最后通过实现地震应力触发科学计算实例验证了方案的有效性,为网格技术在地震领域内的进一步应用作了有意义的探索.      

基于GT3的网格计算平台的构建及应用

基于CEDAGrid项目的需求, 通过应用Web服务技术和网格计算技术, 构建了基于GT3的网格计算应用平台, 在有4个节点的网格计算平台上(其中只有3个节点参与了科学计算任务)进行了地震应力触发计算实例(以川滇地区为例, 共有17个断层), 系统耗时为13min25s, 如果只在一个节点上进行计算, 则需要时间为35min10s。 该实例的成功运行证明, 基于网格技术构建的计算平台可以解决传统的地震领域内的科学计算问题, 在时间效率上可以提高许多, 如果应力触发计算本身越复杂, 网格的计算资源越丰富, 网格计算的执行效率会体现的越优越。 通过这一科学计算试验为网格技术在地震领域的应用做了有意的探索。     

关于网格资源调度的探讨

随着计算机网络技术的飞速发展, 出现了第三代Internet技术网格计算(Grid Computing)。 网格是一种异构、 动态、 分布式的计算环境。 网格资源管理是网格的核心组件, 文中介绍了网格资源管理的特点, 分析了Globus工具集网格资源管理的核心组件。 网格资源调度是网格资源管理的核心部分, 也是网格资源管理中极其繁琐复杂的问题。 基于分析网格计算环境中资源管理的结构特点, 分析和总结了网格调度策略, 提出了多级资源调度模式, 给出了多级资源调度的模式。     

基于网格的地震模拟工作流调度

网格为地震模拟等大规模工作流应用提供了强大的计算能力,如何保障此类应用的软实时性要求是工作流调度问题的一个挑战.利用排队模型来描述网格资源的动态负载压力,提出了评估工作流健壮性量化指标及其计算方法.基于DAG图转换得出的任务执行优先级,根据最大健壮性优先的思想,确定了执行工作流子任务的候选资源;将工作流全局截止时间划分问题描述为一个约束下的非线性规划问题并通过已有方法求解该问题,提出了工作流全局截止时间动态划分方法;最后,提出了一种健壮性增强的地震模拟工作流调度算法RESAESW.仿真实验采用实际地震模拟工作流应用和实际系统数据来验证提出算法的性能表现,实验结果表明本文算法在网格环境的自适应性和地震模拟工作流应用的截止时间要求方面优于其他两个实际网格系统中的调度算法.     

地震预测重点实验室高性能计算平台计费服务系统的设计与实现

本文简要回顾了高性能计算在地震研究中的应用,介绍了中国地震局地震预测重点实验室高性能计算平台及其应用情况,并指出作业计费的必要性。其次,详细介绍了1种高性能计算平台计费服务系统的设计与实现。该系统基于TORQUE和LSF集群作业管理系统日志信息的分析处理,具备作业查询、作业统计、作业计费、作业预约、信息管理、留言管理等功能,有益于解决高性能计算平台计算收费和运维支出平衡的问题,提高了平台的透明度和利用率,促进信息共享和技术交流,并具有较好的扩展性和应用前景。最后,对高性能计算过渡到高性能云计算进行了展望。     

An alternative direct method towards mean dynamic topography computations

An alternative “direct method” to “mean dynamic topography” (MDT) computations using satellite altimetry-derived “mean sea surface” (MSS) and “global geopotential model” (GGM), without direct application of the geoid, is devised. The developed approach, which is based on derivation of an equipotential surface of the gravity field of the Earth that fits to global MSS in least squares sense, is formulated via a constrained optimization problem. The validity of our method is numerically tested by computing a global MDT model based on DNSC08 MSS model and EGM2008 GGM as input data.     

网格技术研究进展

介绍了网格的概念,以及网格体系结构、网格资源管理和网格应用等网格技术的研究现状, 最后对中国地震减灾仿真网格试验系统的研究概况进行了介绍。     

高性能集群计算系统的构建

高性能集群计算系统因其具有强大的运算能力、 较高的I/O性能、 高性能管理和较强的系统扩展能力而广受关注。 本文简要介绍了高性能计算的发展现状和高性能集群计算系统构成, 比较详细地描述了系统环境、 系统性能及集群管理与应用情况, 对高性能计算的新技术进行展望。 并以一个实例阐述了集群系统构建等内容。 该系统由192个计算节点、 4个I/O节点、 1个登录节点、 1个管理节点、 1个Windows集群服务节点、 1组盘阵、 1组计算与I/O网络和1组管理网络组成。 系统计算能力在2009年和2010年中国高性能计算机性能TOP100排行榜分别为第36位和第63位, 具有相当的竞争力, 能够在近几年内为地震科研提供有力支持。     

国家地震应急救援物资综合管理平台建设

针对目前中国地震应急救援物资管理中办法互有差异、系统各有不同,信息化和标准化程度较低,难以统计物资储备的现状,且不易在震后应急救援中对物资进行快速合理地调度等问题,提出数据资源化、业务标准化和服务智能化的解决方案,以此研究并建设面向全国各地震应急救援物资管理单位,统一化的服务于地震应急救援保障与决策的物资信息综合管理平台。该平台已在中国地震应急搜救中心上线,运行结果表明,有效提升了地震应急救援物资综合管理、使用和分析的水平,可解决目前存在的问题,对地震应急救援保障与决策具有重要意义。关键词地震应急;应急救援物资;管理平台     

Ground-motion prediction equations for induced seismicity in the main anticline and main syncline,Upper Silesian Coal Basin,Poland

The purpose of the research was to determine parameters of ground-motion models for two areas characterized by considerable induced seismicity and different geology. Fifty-nine events collected from surface seismological stations of coal mine “Bielszowice” (at the Main Anticline, South Poland) and 144 events from coal mine “Ziemowit” (at the Main Syncline, South Poland) were used for computation. For both areas, simple ground-motion prediction equations (GMPEs) without site effects were derived, but the model was acceptable only for “Bielszowice” area. The GMPE was calculated once again for “Ziemowit”, but this time we took into consideration the amplification coefficient, which significantly improved the model solution. Finally, the theoretical value of amplification was calculated. Knowing that the amplification is associated with subsurface layers, we used three different models of overburden: (i) with Quaternary sediments only, (ii) with a complex of Quaternary-Tertiary sediments, and (iii) with a complex of Quaternary-Tertiary-Triassic sediments and Carboniferous as a basement. Usually, the amplification of vibrations appears in the Quaternary sediments. However, theoretical calculations of amplification were consistent with the results obtained from GMPE when a rigid Carboniferous substratum was applied.