共查询到20条相似文献，搜索用时 46 毫秒
1.
提出一种基于信号整体与局部特征的地震数据自动处理新方法, 该方法不同于以往仅采用包络线互相关来直接检测事件. 新方法依然按照检测、 识别、 关联和定位等4个步骤进行处理, 但在进行单个震相信号检测的同时, 也检测信号波列并利用其包络线特征来识别和关联震相. 文中详细阐述了数据处理过程中如何定义一个波列及抽取和应用其特征. 相关的数据处理技术目前已成功应用于区域台网的日常数据处理分析中. 作为例子, 给出了对新疆区域台网连续16天数据进行测试处理的结果. 实际应用结果表明, 这种新方法可以大幅度降低自动处理结果的误检、 漏检率, 在实际应用中具有很好的前景. 相似文献
2.
Farnoush Forghani‐Arani Jyoti Behura Seth S. Haines Mike Batzle 《Geophysical Prospecting》2013,61(4):778-787
In studies on heavy oil, shale reservoirs, tight gas and enhanced geothermal systems, the use of surface passive seismic data to monitor induced microseismicity due to the fluid flow in the subsurface is becoming more common. However, in most studies passive seismic records contain days and months of data and manually analysing the data can be expensive and inaccurate. Moreover, in the presence of noise, detecting the arrival of weak microseismic events becomes challenging. Hence, the use of an automated, accurate and computationally fast technique for event detection in passive seismic data is essential. The conventional automatic event identification algorithm computes a running‐window energy ratio of the short‐term average to the long‐term average of the passive seismic data for each trace. We show that for the common case of a low signal‐to‐noise ratio in surface passive records, the conventional method is not sufficiently effective at event identification. Here, we extend the conventional algorithm by introducing a technique that is based on the cross‐correlation of the energy ratios computed by the conventional method. With our technique we can measure the similarities amongst the computed energy ratios at different traces. Our approach is successful at improving the detectability of events with a low signal‐to‐noise ratio that are not detectable with the conventional algorithm. Also, our algorithm has the advantage to identify if an event is common to all stations (a regional event) or to a limited number of stations (a local event). We provide examples of applying our technique to synthetic data and a field surface passive data set recorded at a geothermal site. 相似文献
3.
4.
Automatic picking of P and S phases using a neural tree 总被引：2，自引：1，他引：2
The large amount of digital data recorded by permanent and temporary seismic networks makes automatic analysis of seismograms and automatic wave onset time picking schemes of great importance for timely and accurate event locations. We propose a fast and efficient P- and S-wave onset time, automatic detection method based on neural networks. The neural networks adopted here are particular neural trees, called IUANT2, characterized by a high generalization capability. Comparison between neural network automatic onset picking and standard, manual methods, shows that the technique presented here is generally robust and that it is capable to correctly identify phase-types while providing estimates of their accuracies. In addition, the automatic post processing method applied here can remove the ambiguity deriving from the incorrect association of events occurring closely in time. We have tested the methodology against standard STA/LTA phase picks and found that this neural approach performs better especially for low signal-to-noise ratios. We adopt the recall, precision and accuracy estimators to appraise objectively the results and compare them with those obtained with other methodologies.Tests of the proposed method are presented for 342 earthquakes recorded by 23 different stations (about 5000 traces). Our results show that the distribution of the differences between manual and automatic picking has a standard deviation of 0.064 s and 0.11 s for the P and the S waves, respectively. Our results show also that the number of false alarms deriving from incorrect detection is small and, thus, that the method is inherently robust.This paper has not been submitted elsewhere in identical or similar form, nor will it be during the first three months after its submission to Journal of Seismology. 相似文献
5.
6.
利用数字图像处理技术提高地震剖面图像信噪比 总被引：1，自引：2，他引：1
提出了利用数字图像处理技术提高地震剖面信噪比的新方法，首先根据数字图像处理要求的格式，对地震剖面数据进行转换，得到地震剖面图像，分析了地震数据特点和初步地震图像的实验结果后，设计了新的预处理方法——“二维沿层滤波”，在此基础上，利用可以计算帧间运动速度及其变化都较大的改进的光流分析技术，计算出多幅地震剖面对应点的偏移量，然后应用图像积累技术对这多幅地震剖面进行积累，实现对三维地震数据体提高信噪比的处理，该方法充分利用了三维地震信息，不但可以提高整个数据体的信噪比，而且可以减少信号能量的损失，并保持原来的信号能量关系，使地震剖面的质量得到明显提高，为地震解释奠定良好的基础。 相似文献
7.
地震震相拾取是地震数据自动处理的首要环节,包括了信号检测、到时估计和震相识别等过程,震相拾取的准确性直接影响到后续事件关联处理的性能,影响观测报告的质量.为了提高震相拾取的准确性,进而提高观测报告质量,本文采用深度卷积神经网络方法来解决震相拾取问题,构建了多任务卷积神经网络模型,设计了分类和回归的联合损失函数,定义了基于加权的分类损失函数,以三分量地震台站的波形数据作为输入,同时实现对震相的检测识别和到时的精确估计.利用美国南加州地震台网的200万条震相和噪声数据对模型进行训练、验证和测试,对于测试集中直达波P、S震相识别的查全率达到98%以上,到时估计的标准偏差分别为0.067s,0.082s.利用迁移学习和数据增强,将模型用于对我国东北地区台网的6个台站13000条数据的训练、验证和测试中,对该数据集P、S震相查全率分别达到91.21%、85.65%.基于迁移训练后的模型,设计了用于连续数据的震相拾取方法,利用连续的地震数据对该算法进行了实际应用测试,并与国家数据中心和中国地震局的观测报告进行比对,该方法的震相检测识别率平均可达84.5%,验证了该方法在实际应用中的有效性.本文所提出的方法展示了深度神经网络在地震震相拾取中的优异性能,为地震震相和事件的检测识别提供了新的思路. 相似文献
8.
An interesting weak seismic event occurred on August 20, 1983 in the Caspian Sea region in southwestern U.S.S.R., which produced at Gauribidanur array two distinct seismic signals, viz. the leading P signal and a similar signal after an interval of nearly 12 s. Preliminary data corresponding to these two signals led to an initial assumption that the event was a double underground nuclear explosion. However, the focal depth of this event from P arrival-time data of 21 stations was found to be 45 km. Further processing of the Gauribidanur array seismograms confirmed that the second signal, detected exactly 11.85 s after the first arrival, was none other than the surface-reflected P pulse which lent excellent support to the source depth estimate. By using three discriminants, viz. signal complexity, third moment of frequency and relative efficiency of generation of core-reflected P, the event was objectively compared with 10 small presumed underground nuclear explosions in southwestern U.S.S.R. to conclude that the August 20, 1983 event was a shallow-focus small earthquake in the Caspian Sea area. 相似文献
9.
正演模拟是叠前弹性波反演的基础.采用慢度法计算层状介质的叠前地震记录,分别对频率和慢度进行积分变换得到时-空域的地震道集,并对在慢度积分过程中产生的计算噪音提出了解决方案.为得到高精度合成地震记录,需将地层细分,但地层层数很多时,计算量较大;而对地层粗分虽然会大大加快运算速度,但合成记录会丢失很多信息,文中给出了地层的划分原则.该方法能够计算出包括转换波和多次反射在内的全地震响应.但在提高合成记录精度的同时,也导致计算量增大、计算效率降低,因此,本文对基于慢度法全波场模拟进行了并行算法设计,采用计算域分割、工作池并行技术,建立了慢度法全波场正演模拟的并行算法,使得弹性波正演问题求解更加高效,为充分利用叠前地震资料进行叠前反演提供了研究基础. 相似文献
10.
Synthetic seismograms can be very useful in aiding understanding of wave propagation through models of real media, verification of geologic models derived from interpretation of field seismic data, and understanding the nature and complexity of wave phenomena. If meaningful results are to be obtained from synthetic seismograms, the method of their computation must, in general, include three-dimensional geometrical spreading of wavefronts associated with highly concentrated (i.e., point) sources. The method should also adequately represent the seismic response of solid-layered media by including enough primaries, multiples, and converted phases to accurately approximate the total wavefield. In addition to these features, it is also very helpful, although not always essential, if the method of seismogram computation provides for explicit identification of wave type and ray path for each arrival. Various seismograms, computed via asymptotic ray theory and an automatic ray generation scheme, are presented for a highly simplified North Sea velocity structure. This is done to illustrate the importance of the above features and to demonstrate the inadequacy of the plane-wave synthesis method of seismogram computation for point sources and the limitations of acoustic models of solid-layered media. 相似文献
11.
S. V. Baranov 《Izvestiya Physics of the Solid Earth》2007,43(2):177-188
The paper is devoted to such a priority direction of digital seismology as the detection of seismic signals in a noisy medium. Various formulations of the problem of seismic event detection are analyzed. The reasons for the ineffectiveness of the classical formulation of the detection problem are revealed, and an adequate general formulation of the problem in the frequency-time space is proposed. A criterion for the automatic determination of the length of a seismic signal is proposed. The time-scale STA/LTA (Short Time Average to Long Time Average) detector, utilizing the expansion of a seismogram in a packet of wavelets, is proposed. The method eliminates the shortcomings of the standard scheme and enables the detection of seismic events with a low signal-to-noise ratio. The problems of the choice of the mother wavelet and the decomposition level of seismograms are discussed from the standpoint of physical properties of seismic wave fields. 相似文献
12.
The accuracy of automatic procedures for locating earthquakes is influenced by several factors such as errors in picking seismic phases, network geometry, modeling errors and velocity model uncertainties. The main purpose of this work is to improve the performances of the automatic procedure employed for the “quasi-real-time” location of seismic events in North Western Italy by developing a procedure based on a waveform similarity analysis and by using only one seismic station.To detect “earthquake families” a cross-correlation technique was applied to a data set of seismic waveforms recorded in the period 1985-2002, in a small test area (1600 km2) located in the South Western Alps (Italy). Normalized cross-correlation matrices were calculated using about 2700 seismic events, selected on the basis of the signal to noise ratio, manually picked and located by using the Hypoellipse code. The waveform similarity analysis, based on the bridging technique, allowed grouping about 65% of the selected events into 80 earthquake families (multiplets) located inside the area considered. For each earthquake family a master event is selected, manually re-picked and re-located by using Hypoellipse code. Having chosen a reference station (STV) on the basis of the completeness of the available data set, an automatic procedure has been developed with the aim of cross-correlating new seismic recordings (automatically picked) to the waveforms of the events belonging to the detected families. If the new event is proved to belong to a family (on the basis of the cross-correlation values), its hypocenter co-ordinates are defined by the location of the master event of the associated family. The performance of the proposed procedure is tested and demonstrated using a data set of 104 selected earthquakes recorded in the period January 2003-June 2004 and located in the test area. The automatic procedure is able to locate, associating events with the multiplets detected by the waveform similarity analysis, about 50% of the test events, almost independently of the accuracy of the automatic phase picker and without the biasing of the network geometry and of the velocity model uncertainties. 相似文献
13.
基于小波包和峰度赤池信息量准则的P波震相自动识别方法 总被引：2，自引：2，他引：0
基于小波包变换和峰度赤池信息量准则(AIC), 提出了一种新的自动识别P波震相的综合方法, 即小波包-峰度AIC方法． 首先对由加权长短时窗平均比(STA/LTA)法粗略确定的P波到时前后3 s的记录进行小波包三尺度的分解与重构, 分别计算每个尺度重构信号的峰度AIC曲线并将其叠加, 叠加曲线的最小值则为P波震相到时; 然后对原始地震记录进行有限冲激响应自适应滤波以提高信噪比和识别精度; 最后将小波包-峰度AIC方法应用到合成理论地震图及实际地震记录的P波初至自动识别中. 结果表明: 初至清晰度对识别精度的影响比信噪比对其影响更大; 与单独使用加权STA/LTA方法和峰度AIC法相比, 小波包-峰度AIC法具有更强的抗噪能力, 识别精度更高; 当初至清晰时, 小波包-峰度AIC法自动识别与人工识别的P波到时平均绝对差值为(0.077±0.075) s． 相似文献
14.
The automatic detection of geological features such as faults and channels is a challenging problem in today's seismic exploration industry. Edge detection filters are generally applied to locate features. It is desirable to reduce noise in the data before edge detection. The application of smoothing or low‐pass filters results in noise suppression, but this causes edge blurring as well. Edge‐preserving smoothing is a technique that results in simultaneous edge preservation and noise suppression. Until now, edge‐preserving smoothing has been carried out on rectangular sampled seismic data. In this paper, an attempt has been made to detect edges by applying edge‐preserving smoothing as a pre‐processing step in the hexagonally sampled seismic‐data spatial domain. A hexagonal approach is an efficient method of sampling and has greater symmetry than a rectangular approach. Here, spiral architecture has been employed to handle the hexagonally sampled seismic data. A comparison of edge‐preserving smoothing on both rectangular and hexagonally sampled seismic data is carried out. The data used were provided by Saudi Aramco. It is shown that hexagonal processing results in well‐defined edges with fewer computations. 相似文献
15.
16.
We present a new approach to enhancing weak prestack reflection signals without sacrificing higher frequencies. As a first step, we employ known multidimensional local stacking to obtain an approximate ‘model of the signal’. Guided by phase spectra from this model, we can detect very weak signals and make them visible and coherent by ‘repairing’ corrupted phase of original data. Both presented approaches – phase substitution and phase sign corrections – show good performance on complex synthetic and field data suffering from severe near-surface scattering where conventional processing methods are rendered ineffective. The methods are mathematically formulated as a special case of time-frequency masking (common in speech processing) combined with the signal model from local stacking. This powerful combination opens the avenue for a completely new family of approaches for multi-channel seismic processing that can address seismic processing of land data with nodes and single sensors in the desert environment. 相似文献
17.
Spectral classification methods in monitoring small local events by the Israel seismic network 总被引：1，自引：1，他引：0
We use the dense Israel Seismic Network (ISN) to discriminate between low magnitude earthquakes and explosions in the Middle East region. This issue is important for CTBT monitoring, especially when considering small nuclear tests which may be conducted under evasive conditions. We explore the performance of efficient discriminants based on spectral features of seismograms using waveforms of 50 earthquakes and 114 quarry and underwater blasts with magnitudes 1.0–2.8, recorded by ISN short-period stations at distances up to 200 km. The single-station spectral ratio of the low and high-frequency seismic energy shows an overlap between explosions and earthquakes. After averaging over a subnet of stations, the resolving power is enhanced and the two classes of events are separated. Different frequency bands were tested; the (1–3 Hz)/(6–8 Hz) ratio provided the best discriminant performance. We also estimated normalized r.m.s. spectral amplitudes in several sequential equal frequency windows within the 1–12 Hz band and applied multiparametric automatic classification procedures (Linear Discrimination Function and Artificial Neural Network) to the amplitudes averaged over a subnetwork. A leave-one-out test showed a low rate of error for the multiparametric procedures. An innovative multi-station discriminant is proposed, based on spectral modulation associated with ripple-firing in quarry blasts and with the bubbling effect in underwater explosions. It utilizes a distinct azimuth-invariant coherency of spectral shapes for different stations in the frequency range (1–12 Hz). The coherency is measured by semblance statistics commonly used in seismic prospecting for phase correlation in the time domain. After modification, the statistics applied to the network spectra provided event separation. A new feature of all the above mentioned procedures is that they are based on smoothed (0.5 Hz window), instrument-corrected FFT spectra of the whole signal; they are robust to the accuracy of onset time estimation and, thus well suited to automatic event identification. 相似文献
18.
提出一种自适应协方差的时频域极化滤波方法。该方法在广义S变换时频方法的基础上,构造时频域自适应协方差矩阵,通过特征分析计算时频域瞬时极化参数,设计极化滤波器,实现多分量地震极化分析和滤波。其优势在于协方差矩阵的分析时窗的长度由多分量地震数据的瞬时频率确定,可以自适应于有效信号的周期,在每个时频点计算极化参数不需要进行插值处理;结合时间频率信息,解决在时间域或频率域波形或频率重叠的信号具有明显的直观性。模型数据及实际三分量台站地震数据处理结果表明,该极化滤波方法在台站地震资料分析和处理方面具有很好的直观性和较高的分辨率。 相似文献
19.
提出了一种新的基于信号子空间分解技术的3 D地震资料高分辨率处理方法.利用信号子空间分解技术,不仅可以分离3 D地震资料中的信号和噪声,而且可以进一步根据地震同相轴的倾角不同,将混合信号进行分解得到单个同相轴信号.通过丢弃噪声子空间,只对不同信号子空间重构的信号利用谱白化技术进行高分辨率处理,然后累加所有处理结果,从而达到既提高地震资料的分辨率,又提高地震资料信噪比的目的.对合成资料和实际地震资料的处理结果表明,此法具有好的应用前景. 相似文献
20.
The receiver function (RF) technique is a well-established method to investigate the crustal and upper mantle structures based on three-component seismograms of teleseismic events. In the present study, we propose a modified automatic procedure to determine the back azimuth and polarization angles of a teleseismic event based on the RF technique. The method is tested for the recording of 3 permanent and 3 temporary broadband seismic stations located in the vicinity of Poland. Additionally, the analysis of Rayleigh wave polarization is conducted to show that the new procedure is not sensitive to incorrect seismometer orientation. The synthetic modelling of RF by a modified ray-tracing method for 2.5D models beneath each seismic station down to a depth of 60 km is performed to show the effectiveness of the proposed method in the calculation of RF for a complex structure with dipping layers. 相似文献