历史地震图数据库及共享平台

历史地震图纸是记录历史地震的重要资料,在进行电子化之前,图纸的唯一性限制了其大范围共享的可能。电子化处理是抢救和保护历史地震图纸的重要手段,可以有效发挥其使用及研究价值,具有重要意义。为此研制历史地震图数据库及共享平台,对北京国家地球观象台现存的60万张历史地震图纸进行高精度(600 dpi)扫描和规则化命名,利用数据库和FastDFS分布式文件系统进行统一存储,并采用多级缩略图和网络服务方式提供数据共享服务。     

Data acquisition and pre-processing required for simultaneous P-SV inversion

The goal of seismic reflection surveys is the derivation of petrophysical subsurface parameters from surface measurements. Today's well established technique in data acquisition, as well as processing terms, is based on the acoustic approximation to the real world's wave propagation. In recent years a lot of work has been done to extend the technique to the elastic approximation. There was especially an important trend towards elastic inversion techniques operating on plane-wave seismograms, called simultaneous P-SV inversion (or short P-SV inversion) within this paper. Being still under investigation, some important aspects of P-SV inversion concerning data acquisition as well as pre-processing, should be pointed out. To fit the assumptions of P-SV inversion schemes, at least a two-dimensional picture of the reflected wavefield with vertical and in-line horizontal receivers has to be recorded. Moreover, the theoretical work done suggests that in addition to a survey with a compressional wave source, a second survey should be done using sources radiating vertically polarized shear waves, is needed. Finally, proper slant stacking must be performed to get plane-wave seismograms. The P/S separated plane-wave seismograms are then well prepared for feeding into the inversion algorithms. P/S separated planewave seismograms are then well prepared for feeding into the inversion algorithm.s In this paper, a tutorial overview of the data acquisition and pre-processing in accordance with the P-SV inversion philosophy is given and illustrated using synthetic seismograms. A judgement on the feasibility of the P-SV inversion philosophy must be left to ongoing research.     

模拟地震图纸数字化存储信息录入的设计与实现

本文针对中国地震局正在实施的模拟地震资料抢救项目,对我国现存各类模拟地震图纸进行分类,设计实现了适用于模拟地震图纸扫描基础信息录入软件。首先,通过对图纸图章的分析设计了软件的需求功能,并确定了界面风格;其次,根据图纸扫描工作的保密属性,设计运用xml文件格式将扫描图纸的信息保存在本地文件中;最后,设计了图纸信息录入校验规则以及图纸的命名规则,并实现了异常自动提取和一键生成报表等功能,而且易于使用。本文开发设计的软件不仅降低了大量信息人工录入的错误率,而且显著提升了模拟地震图纸扫描录入与校验的工作效率。通过这项工作所获得的成果数据,便于地震科研工作效率的提高。      

Digitized Database of Old Seismograms Recorder in Romania

The aim of this paper is to describe a managing system for a unique Romanian database of historical seismograms and complementary documentation (metadata) and its dissemination and analysis procedure. For this study, 5188 historical seismograms recorded between 1903 and 1957 by the Romanian seismological observatories (Bucharest-Filaret, Foc?ani, Bac?u, Vrincioaia, Câmpulung-Muscel, Ia?i) were used. In order to reconsider the historical instrumental data, the analog seismograms are converted to digital images and digital waveforms (digitization/ vectorialisation). First, we applied a careful scanning procedure of the seismograms and related material (seismic bulletins, station books, etc.). In a next step, the high resolution scanned seismograms will be processed to obtain the digital/numeric waveforms. We used a Colortrac Smartlf Cx40 scanner which provides images in TIFF or JPG format. For digitization the algorithm Teseo2 developed by the National Institute of Geophysics and Volcanology in Rome (Italy), within the framework of the SISMOS Project, will be used.     

应用人工神经元网络方法识别近震与远震

利用人工神经元网络方法,提出了一种在实时系统中识别近震与远震的算法,并用广东省台网记录的地震波形数据进行训练和测试。结果表明,该算法能准确可靠地识别远、近震,可以应用于台网的实时处理系统。     

GDS－1000宽频带数字地震仪野外数据回收处理系统

在笔记本计算机上运行的ＧＤＳ－１０００宽频带数字地震仪地震数据回收处理系统软件，用于在野外条件下执行现场仪器参数设置、数据回收、数据解编及地震图的开窗显示等功能。近一年的野外实际应用结果表明，该系统软件具有数据通讯快速可靠、使用灵活方便等优点。     

A novel technique for automatic seismic data processing using both integral and local feature of seismograms

A novel technique for automatic seismic data processing using both integral and local feature of seismograms was presented in this paper. Here, the term integral feature of seismograms refers to feature which may depict the shape of the whole seismograms. However, unlike some previous efforts which completely abandon the DIAL approach, i.e., signal detection, phase identification, association, and event localization, and seek to use envelope cross-correlation to detect seismic events directly, our technique keeps following the DIAL approach, but in addition to detect signals corresponding to individual seismic phases, it also detects continuous wave-trains and explores their feature for phase-type identification and signal association. More concrete ideas about how to define wave-trains and combine them with various detections, as well as how to measure and utilize their feature in the seismic data processing were expatiated in the paper. This approach has been applied to the routine data processing by us for years, and test results for a 16 days’ period using data from the Xinjiang seismic station network were presented. The automatic processing results have fairly low false and missed event rate simultaneously, showing that the new technique has good application prospects for improvement of the automatic seismic data processing.     

频率域中的合成地震记录

波动理论的合成地震记录是地震资料处理和解释中的重要研究课题。 目前,制作波动理论的合成记录主要是用Trorey的克希霍夫方法,正在研究的有有限差分法和有限元素法。本文提出一种新的方法,在频率域中制作波动理论的合成记录,给出了与克希霍夫绕射波方程等价的频率域中的绕射波方程。在国产DJS-11型计算机上,应用两种绕射波方程实现了合成地震记录,并且进行了比较,认为频率域方法是有效的,可以用于波动方程偏移和地震资料解释,进一步工作有可能获得多层变速模型的波动理论记录。     

单台识别核爆的方法初探

根据乌鲁木齐地震台记录的国外６个主要核试验基地的核爆地震记录资料，通过对振动波形、首波振幅Ｐ１、最大纵波振幅Ｐｍ、最大横波振幅Ｓｍ及体波震级ｍｂ、面波震级Ｍｓ之间一些关系的研究，找出了核爆记录与天然地震记录的差别，为单台正确识别核爆提供了基本方法。     

短周期仪器记录在远震分析中的应用

地震泼由多种频率成份组成,在其传播过程中。由于地球介质的滤波作用,随着传播距离的增大。高频成份逐渐被衰减。周期越大的成份传播越远。因此,在地震观测上,不同频带的地震仪用于观测不同震中距范围的地震,使仪器较好地响应,就我国台网而言,短周期仪用于地方震、近震的监测,中长周期仪和长周期仪用于远震,极远震的监测,其中长周期仪偏重于记录极远震。但作者在实际分析工作中注意到,短周期仪的记录应用于远震,极远震的分析,显示出独特的作用,充分挖掘和利用短仪资料,将会更加丰富远震、极远震的震相资料。本文以高台地震台的资料为依据,从震中距和震级角度,与中长仪和长仪作比较,对短仪记录远震、极远震的震相及特征作了统计分析与初步探讨,并对新疆地区“影区”地震S波的记录、日本地区地震_PPcP等震相的出现提出了现象的存在,给出了相应的观测结果,以待步一步研究。     

Teleseismic receiver function using stacking and smoothing of multi seismic-records at a single station

Receiver functions (RFs) obtained using teleseismic wave records at a seismic station and synthetic seismograms indicate that RF with a single teleseismic wave record is related to the selection of record section and to the calculating parameters of the RF. The scatter noise contained in the seismogram also affects the quality of RF. A new method for calculating receiver function, stacking and smoothing multi-seismic records in a single station, is presented in this paper. The RF results using some records and some synthetic seismograms with different noises indicate prominent mantle discontinuity and thus prove that the method is effective and satisfied.     

A cloud-based synthetic seismogram generator implemented using Windows Azure

Synthetic seismograms generated by solving the seismic wave equation using numerical methods are being widely used in seismology. For fully three-dimensional seismic structure models, the generation of these synthetic seismograms may require large amount of computing resources. Conventional high-performance computer clusters may not provide a cost-effective solution to this type of applications. The newly emerging cloud-computing platform provides an alternative solution. In this paper, we describe our implementation of a synthetic seismogram generator based on the reciprocity principle using the Windows Azure cloud application framework. Our preliminary experiment shows that our cloud-based synthetic seismogram generator provides a cost-effective and numerically efficient approach for computing synthetic seismograms based on the reciprocity principle.     

密度和压缩系数的散射层析成像法

本文在速度成像的基础上研究了同时对密度和压缩系数成像的散射波层析成像法.对不同散射角度的计算可以得到一系列反演图像,拟合这些图像,从而可以有效地达到对密度和压缩系数（或速度）成像的目的.与单纯的速度成像相比,增加了反演的难度.首先是对资料的方位性要求增加;其次是对资料的利用率下降.即便如此,从对较少量的炮点和检波点资料的数值计算来看,仍取得了满意的成像结果.我们对组成字母“A”的散射体结构进行了成像计算,结果能够同时再现密度和压缩系数,成像清晰,表明了方法的可行性,并能应用于复杂结构的成像问题.     

Focus-of-attention techniques in the automatic interpretation of seismograms

The focus-of-attention techniques implemented in SNA2, a knowledge-based system for seismogram interpretation, are presented. They consist of data compression of the input digital records, scanning of the compressed traces to detect candidate seismograms and extraction of seismogram features. A criterion is given to rate the clarity of seismograms; the clarity defines the order in which the system will consider them to build up the interpretation. The proposed techniques are simple and fast; they allow quick rejection of noise and focussing the attention of the system on the portions of traces containing relevant information.     

弱各向异性介质弹性波的准各向同性近似正演模拟

准各向同性(QI)近似可用于弱各向异性介质的正演模拟.本文通过运用QI方法的零阶和一阶近似,计算了VTI介质模型的地震记录.得出的地震记录与标准各向同性射线理论(IRT)和基于伪谱法的三维地震正演模拟得出的地震记录作了比较,可以认为是精确的合成地震记录.     

地震转换波对叠前同时反演的影响研究

叠前同时反演是油气探测的一种有效工具.其理论基础是平面P波Zoeppritz方程计算的反射系数的近似,是入射角的函数.叠前同时反演可以利用三项或两项Fatti方程进行反演分析.本文针对实际油田的测井数据,利用反射率法模拟了仅包含P波一次反射记录,包含P波一次反射和P波层间多次波记录以及全波场地震记录,再利用叠前同时反演对合成地震记录进行反演研究.研究结果表明,在大偏移距处P波主要反射受到其它模式波的污染,从而影响了叠前同时反演结果的精度.对于薄互层介质当转换波影响严重时,使用小角度数据的两项AVO反演比使用大角度数据的三项AVO反演更合理可靠.     

Determining the focal mechanisms and depths of relatively small earthquakes using a few stations by full‐waveform modelling

Determining the focal mechanism of earthquakes helps us to better define faults and understand the stress regime. This technique can be helpful in the oil and gas industry where it can be applied to microseismic events. The objective of this paper is to find double couple focal mechanisms, excluding scalar seismic moments, and the depths of small earthquakes using data from relatively few local stations. This objective is met by generating three‐component synthetic seismograms to match the observed normalized velocity seismograms. We first calculate Green's functions given an initial estimate of the earthquake's hypocentre, the locations of the seismic recording stations and a 1D velocity model of the region for a series of depths. Then, we calculate the moment tensor for different combinations of strikes, dips and rakes for each depth. These moment tensors are combined with the Green's functions and then convolved with a source time function to produce synthetic seismograms. We use a grid search to find the synthetic seismogram with the largest objective function that best fits all three components of the observed velocity seismogram. These parameters define the focal mechanism solution of an earthquake. We tested the method using three earthquakes in Southern California with moment magnitudes of 5.0, 5.1 and 4.4 using the frequency range 0.1–2.0 Hz. The source mechanisms of the events were determined independently using data from a multitude of stations. Our results obtained, from as few as three stations, generally match those obtained by the Southern California Earthquake Data Center. The main advantage of this method is that we use relatively high‐frequency full‐waveforms, including those from short‐period instruments, which makes it possible to find the focal mechanism and depth of earthquakes using as few as three stations when the velocity structure is known.     

关于“反射法”计算综合地震图的实施方案

本文针对反射法计算综合地震图中的运算速度和存储安排等问题,提出一种有效的计算实施方案。文中附有三种模型的体波综合地震图计算实例。     

Application of the Maslov Seismogram Method in Three Dimensions

Asymptotic methods provide an efficient way to compute seismograms in heterogeneous media. However, zeroth-order ray theory, the simplest of the asymptotic methods, often fails because of the presence of caustics. Maslov theory is an extension of zeroth-order ray theory, which gives a uniformly valid expression of the wavefield everywhere, including the caustics. This result is given in terms of an integral of ray data over one or two ray parameters. It is shown in this paper how geometrical arrivals are constructed in the one and two-parameter Maslov integrals.In practice Maslov seismograms have been computed using only one ray parameter. However, in three-dimensional media two parameters are needed to uniquely define a ray. In this paper we present an efficient algorithm to compute two-parameter Maslov integrals. The Maslov integral is evaluated by computing the frequency-to-time Fourier transform prior to integration over the ray parameters. The wavefield is then discretized by smoothing with a boxcar function. The resulting expression, which only requires the results of ordinary kinematic and dynamic ray tracing, cen be computed efficiently and robustly. A numerical example is given that illustrates the use of this algorithm.     

终止相的理论地震图-兼论唐山地震震源破裂过程

本文首先较为系统地回顾了终止相研究的历史、进展,对震源破裂过程产生终止震相的物理机制进行了讨论,并详细讨论了主事件分析法以及合成包含终止相的理论图的方法。在此基础上应用主事件分析法处理了唐山地震的实际资料,给出唐山地震的震源破裂过程的大致时空图象,同时计算了包括终止事件的合成图。结果表明,考虑破裂过程的终止事件可以较好地拟合P波观测图,从而进一步确认唐山地震震源破裂过程包括两次终止事件的可能性。此外本文还对唐山地震的震源环境以及调整单元的让位情况进行了讨论。