共查询到20条相似文献，搜索用时 31 毫秒
1.
微地震事件初至拾取是井下微地震监测数据处理的关键步骤之一.初至误差的存在会使微地震震源定位结果产生较大偏差,进而影响后续的压裂裂缝解释.通常初至拾取过程对所有的微地震事件选择相同的特征函数并采用一致的拾取参数进行统一处理,然而当事件的能量、震源机制、传播路径以及背景噪声等存在明显差异时,所得初至拾取结果差别显著.为了提高微地震事件初至拾取标准一致性,本文提出基于波形相似特征的初至拾取及全局校正方法.该方法首先利用互相关函数对每个事件内的各道记录进行时差校正,得到初始初至信息并形成叠加道,再对所有事件的叠加道进行全局互相关得到事件间初至相对校正量,最终初至结果可以通过各个事件的初始初至信息与其相对校正量相加得到.方法将所有微地震事件初至结果作为一个整体处理,从而能够克服常规方法初至拾取标准一致性差的缺陷.实际资料处理结果表明,相比于常规方法,该方法可以有效提高事件初至拾取和定位结果的一致性. 相似文献
2.
P. L. Swanson L. H. Estey F. M. Boler S. Billington 《Pure and Applied Geophysics》1992,139(3-4):375-404
The accuracy and precision of microseismic event locations were measured, analyzed, and compared for two types of location systems: anolog and digital. In the first system, relative times of first arrival were estimated from analog signals using automated hardware circuitry; station positions were estimated from mine map coordinates; and event locations were determined using the BLD (Blake, Leighton, and Duvall) direct solution method. In the second system, arrival times were manually measured during interactive displays of digital waveforms; station coordinates were surveyed; and the SW-GBM (Salamon and Wiebols; Godson, Bridges, and McKavanagh) direct basis function was used to solve for locations. Both systems assume constant isotropic seismic velocity of slightly different signals data sets, calibration blast signals with known source site and origin time, and microseismic event signals, were recorded by each location system employing the same array of high-frequency (5 kHz) accelerometers with 150 m maximum dimension. The calibration blast tests indicated a location precision of ±2 m and accuracy of ±10 m for the analog system. Location precision and accuracy for the digital system measured ±1 m and ±8 m, respectively. Numerical experiments were used to assess the contributions of errors in velocity, arrival times, and station positions on the location accuracy and precision for each system. Measured and estimated errors appropriate to each system for microseismic events were simulated in computing source locations for comparison with exact synthetic event locations. Discrepancy vectors between exact locations and locations calculated with known data errors averaged 7.7 and 1.4 m for the analog and digital systems, respectively. These averages are probably more representative of the location precision of microseismic events, since the calibration blast tests produce impulsive seismic arrivals resulting in smaller arrival-time pick errors in the analog system. For both systems, location accuracy is limited by inadequate modeling of the velocity structure. Consequently, when isotropic velocity models are used in the travel-time inversions, the increased effort expended with the digital location system does not, for the particular systems studied, result in increased accuracy. 相似文献
3.
A challenge in microseismic monitoring is quantification of survey acquisition and processing errors, and how these errors jointly affect estimated locations. Quantifying acquisition and processing errors and uncertainty has multiple benefits, such as more accurate and precise estimation of locations, anisotropy, moment tensor inversion and, potentially, allowing for detection of 4D reservoir changes. Here, we quantify uncertainty due to acquisition, receiver orientation error, and hodogram analysis. Additionally, we illustrate the effects of signal to noise ratio variances upon event detection. We apply processing steps to a downhole microseismic dataset from Pouce Coupe, Alberta, Canada. We use a probabilistic location approach to identify the optimal bottom well location based upon known source locations. Probability density functions are utilized to quantify uncertainty and propagate it through processing, including in source location inversion to describe the three-dimensional event location likelihood. Event locations are calculated and an amplitude stacking approach is used to reduce the error associated with first break picking and the minimization with modelled travel times. Changes in the early processing steps have allowed for understanding of location uncertainty of the mapped microseismic events. 相似文献
4.
Locating microseismic events using borehole data 总被引：1，自引：0，他引：1
Constraining microseismic hypocentres in and around hydrocarbon reservoirs and their overburdens is essential for the monitoring of deformation related to hydraulic fracturing, production and injection and the assessment of reservoir security for CO2 and wastewater storage. Microseismic monitoring in hydrocarbon reservoirs can be achieved via a variety of surface and subsurface acquisition geometries. In this study we use data from a single, subsurface, vertical array of sensors. We test an existing technique that uses a 1D velocity model to constrain locations by minimizing differential S‐to‐P arrival times for individual sensors. We show that small errors in either arrival time picks or the velocity model can lead to large errors in depth, especially near velocity model discontinuities where events tend to cluster. To address this issue we develop two methods that use all available arrival times simultaneously in the inversion, thus maximizing the number of potential constraints from to N, where N is the number of phase picks. The first approach minimizes all available arrival time pairs whilst the second approach, the equal distance time (EDT) method defines the hypocentre as the point where the maximum number of arrival time surfaces intersect. We test and compare the new location procedures with locations using differential S‐to‐P times at each individual sensor on a microseismic data set recorded by a vertical array of sensors at the Ekofisk reservoir in the North Sea. Specifically, we test each procedure's sensitivity to perturbations in measured arrival times and the velocity model using Monte Carlo analysis. In general, location uncertainties increase with increasing raypath length. We show that errors in velocity model estimates are the most significant source of uncertainty in source location with these experiments. Our tests show that hypocentres determined by the new procedures are less sensitive to erroneous measurements and velocity model uncertainties thus reducing the potential for misinterpretation of the results. 相似文献
5.
《地震研究进展（英文）》2021,1(3):100008
In the presence of background noise, arrival times picked from a surface microseismic data set usually include a number of false picks that can lead to uncertainty in location estimation. To eliminate false picks and improve the accuracy of location estimates, we develop an association algorithm termed RANSAC-based Arrival Time Event Clustering (RATEC) that clusters picked arrival times into event groups based on random sampling and fitting moveout curves that approximate hyperbolas. Arrival times far from the fitted hyperbolas are classified as false picks and removed from the data set prior to location estimation. Simulations of synthetic data for a 1-D linear array show that RATEC is robust under different noise conditions and generally applicable to various types of subsurface structures. By generalizing the underlying moveout model, RATEC is extended to the case of a 2-D surface monitoring array. The effectiveness of event location for the 2-D case is demonstrated using a data set collected by the 5200-element dense Long Beach array. The obtained results suggest that RATEC is effective in removing false picks and hence can be used for phase association before location estimates. 相似文献
6.
震源定位是微地震监测技术要解决的主要问题.目前,井下微地震监测多采用走时拟合法计算震源位置.常规方法受到环境噪声、初至拾取误差、速度模型误差等因素的影响,定位结果存在一定误差.为了提高定位精度,本文提出了一种改进的基于网格搜索的微地震震源定位方法.本文方法根据P波的偏振特征参数计算概率密度函数求取震源方位角,并采用改进的目标函数和搜索算法计算震源的径向距离和深度.模型数据和实际资料的处理结果表明,本文方法具有较强的抗噪性,计算得到的震源方位角更加接近真实值;与常规目标函数相比,本文方法采用的目标函数具有更好的收敛性,其定位结果受初至拾取误差和速度模型误差的影响更小;本文提出的搜索算法能够消除由于错误拾取造成的观测到时中的异常值对定位结果的影响. 相似文献
7.
We present a robust method for the automatic detection and picking of microseismic events that consists of two steps. The first step provides accurate single-trace picks using three automatic phase pickers adapted from earthquake seismology. In the second step, a multi-channel strategy is implemented to associate (or not) the previous picks with actual microseismic signals by taking into account their expected alignment in all the available channels, thus reducing the false positive rate. As a result, the method provides the number of declared microseismic events, a confidence indicator associated with each of them, and the corresponding traveltime picks. Results using two field noisy data records demonstrate that the automatic detection and picking of microseismic events can be carried out with a relatively high confidence level and accuracy. 相似文献
8.
R. P. Young S. Talebi D. A. Hutchins T. I. Urbancic 《Pure and Applied Geophysics》1989,129(3-4):455-474
Rockbursts and mining-induced seismic events have serious socio-economic consequences for the Canadian mining industry, as their mines are extended to greater depths. Automatic multichannel monitoring systems (Electro-Lab MP250s) are routinely, used to detect the arrival times of seismic waves radiated by mining-induced events and sensed on an array of single component transducers installed throughout a mine. These arrival times are then used to locate the events and produce maps of areas of high activity for use in mine planning and design. This approach has limitations in that, it does not allow a detailed analysis of source mechanisms, which could be extracted if whole waveform signals are recorded and analyzed.A major research project, sponsored by the Natural Sciences and Engineering Research Council of Canada (NSERC) with the collaboration of the Canadian mining industry, is aimed at enhancing existing mine seismic monitoring technology in Canada, in order to carry out more advanced processing of data to obtain fundamental scientific information on mining-induced seismic events This paper describes preliminary results from seismic monitoring experiments carried out in a hard rock nickel mine in Sudbury, Canada. Existing seismic monitoring instrumentation was enhanced with a low cost microcomputer-based whole waveform seismic acquisition system. Some of the signals recorded during this experiment indicate anisotropic wave propagation through the mine rock masses, as observed by the splitting of shear waves and the relative arrival of two shear waves polarized in directions which may be related to the structural fabric and/or state of stress in the rock mass. Analysis of compressional wave first motion shows the predominance of shear events, as indicated by focal mechanism studies and is confirmed by spectral analysis of the waveforms. The source parameters were estimated fro typical low magnitude localized microseismic events during the initial monitoring experiments. The seismic moment of these events varied between 106 N.m and 2.108 N.m. with a circular source radius of between 1 m and 2 m with an estimated stress drop of the order of 1 MPa. 相似文献
9.
Anisotropy Effects on Microseismic Event Location 总被引：2，自引：0，他引：2
Seismic anisotropy in sedimentary environments is significant—microseismic waveforms often show strong shear-wave splitting,
with differences reaching 40% between horizontally and vertically-polarized shear-wave velocities. Failure to account for
this anisotropy is shown to result in large microseismic event location errors. A method is presented here for determining
the five elastic parameters of a homogeneous, vertical transverse-isotropic (VTI) model from calibration shot data. The method
can also use data from mining-induced seismic events, which are then simultaneously located. This simple model provides a
good fit to arrival times from coal-environment data, and results in dramatic shifts in interpreted event locations. 相似文献
10.
分析矿震破裂机制及微震的时空分布能够为矿区灾害评估提供更多的有效信息.本研究基于密集台阵观测对2019年11月12日辽宁抚顺2.4级矿震开展震源参数研究,震源机制解显示地震破裂包含明显的非双力偶分量,表现为体积压缩的塌陷机制,且震源深度较浅,最佳拟合矩心深度为0.6 km.同时,对11月3日—25日记录的连续地震波形开展微震扫描,新检测出324个微震事件(-0.5~2.0级),定位结果显示在M 2.4矿震发生前M>1.0级事件显著增多,且在矿震位置存在近南北向的微震条带分布,微震序列随时间向深部迁移(约1.5 km),暗示存在断层活化迹象.结合震源破裂机制,我们认为此次事件与矿区塌陷破裂密切相关,同时伴随先存断裂的剪切滑动.本研究表明,基于密集台阵观测的地震矩张量反演和微震检测研究,对判定矿震类型及防范矿区灾害具有重要的研究意义. 相似文献
11.
Automatic pickings in earthquake real-time monitoring systems often contain noise bursts and/or phases of different event(s)
occurring almost simultaneously. Typically, a locator uses these picks as P and S waves arrival times coming from a single
event and, therefore, should be complemented by a distinctive phase association logic. The method we propose manages to automatically
associate data related to different events and eliminates the influence of spoiled data from single events. The method is
based on “network beamforming”, a robust and stable algorithm, which utilizes a hypocenter grid search for the stack maximum
of a set of complex exponents applied to the P phase readings. The algorithm separates the residual outliers and then uses
them for location. If successful, a hypocenter is established for the interfering event. The solutions obtained are overall
robust and independent from the estimate of origin times. The preliminary epicenter for the grid search is provided by the
intersection of perpendicular bisectors in the modified “arrival order algorithm” or by the modified “Tnow” algorithm, which uses non-arrival information. We applied this method to automatic first arrival phase readings of 915 events
registered by the Hi-net Japan seismic network and our results are statistically promising. Here, we present two interesting
and complicated examples. 相似文献
12.
煤矿冲击矿压现象日益严峻,作为矿山动力灾害的主要监测手段,微震监测系统已在许多矿井广泛使用,为保证矿震定位和能量计算的准确,提高预测预报煤矿冲击矿压的可行性,应建立一套台网布设优化及评价系统.应用微震定位和D值优化设计理论,结合煤矿实际条件研究了影响矿震定位精度的主要因素和不利条件,并提出了采用综合指数法确定煤矿高微震活动区域和区域内矿震发生的概率,制定了台站候选点和监测区域确定的一般原则.通过理论分析震中和震源标准差反映台网定位能力的不足,建立基于数值仿真实验方法的震中与震源误差期望值模型,最终形成台网布设优化及评价系统.实验和现场应用结果表明,该系统能够快速确定台网最优布设方案,准确评价台网定位能力,满足煤矿微震监测的需要. 相似文献
13.
Yang Zhang Leo Eisner William Barker Michael C. Mueller Kevin L. Smith 《Geophysical Prospecting》2013,61(5):919-930
We develop a methodology to obtain a consistent velocity model from calibration shots or microseismicity observed on a buried array. Using a layered 1D isotropic model derived from checkshots as an initial velocity model, we invert P‐wave arrival times to obtain effective anisotropic parameters with a vertical axis of symmetry (VTI). The nonlinear inversion uses iteration between linearized inversion for anisotropic parameters and origin times or depths, which is specific to microseismic monitoring. We apply this technique to multiple microseismic events from several treatments within a buried array. The joint inversion of selected events shows a largely reduced RMS error indicating that we can obtain robust estimates of anisotropic parameters, however we do not show improved source locations. For joint inversion of multiple microseismic events we obtained Thomsen anisotropic parameters ε of 0.15 and δ of 0.05, which are consistent with values observed in active seismic surveys. These values allow us to locate microseismic events from multiple hydraulic fracture treatments separated across thousands of metres with a single velocity model. As a result, we invert the effective anisotropy for the buried array region and are able to provide a more consistent microseismicity mapping for past and future hydraulic fracture stimulations. 相似文献
14.
Accurate Location of Synthetic Acoustic Emissions and Location Sensitivity to Relocation Methods, Velocity Perturbations, and Seismic Anisotropy 总被引：1，自引：0，他引：1
G. A. Jones S. E. J. Nippress A. Rietbrock J. M. Reyes-Montes 《Pure and Applied Geophysics》2008,165(2):235-254
Acoustic emission (AE) monitoring is a non-invasive method of monitoring fracturing both in situ, and in experimental rock deformation studies. Until recently, the major impediment for imaging brittle failure within a
rock mass is the accuracy at which the hypocenters may be located. However, recent advances in the location of regional scale
earthquakes have successfully reduced hypocentral uncertainties by an order of magnitude. The least-squares Geiger, master
event relocation, and double difference methods have been considered in a series of synthetic experiments which investigate
their ability to resolve AE hypocentral locations. The effect of AE hypocenter location accuracy due to seismic velocity perturbations,
uncertainty in the first arrival pick, array geometry and the inversion of a seismically anisotropic structure with an isotropic
velocity model were tested. Hypocenters determined using the Geiger procedure for a homogeneous, isotropic sample with a known
velocity model gave a RMS error for the hypocenter locations of 2.6 mm; in contrast the double difference method is capable
of reducing the location error of these hypocenters by an order of magnitude. We test uncertainties in velocity model of up
to ±10% and show that the double difference method can attain the same RMS error as using the standard Geiger procedure with
a known velocity model. The double difference method is also capable of precise locations even in a 40% anisotropic velocity
structure using an isotropic model for location and attains a RMS mislocation error of 2.6 mm that is comparable to a RMS
mislocation error produced with an isotropic known velocity model using the Geiger approach. We test the effect of sensor
geometry on location accuracy and find that, even when sensors are missing, the double difference method is capable of a 1.43
mm total RMS mislocation compared to 4.58 mm for the Geiger method. The accuracy of automatic picking algorithms used for
AE studies is ±0.5 μs (1 time sample when the sampling rate is 0.2 μs). We investigate how AE locations are effected by the
accuracy of first arrival picking by randomly delaying the actual first arrival by up to 5 time samples. We find that even
when noise levels are set to 5 time samples the double difference method successfully relocates the synthetic AE. 相似文献
15.
An application of waveform denoising for microseismic data using polarization–linearity and time–frequency thresholding 
下载免费PDF全文
下载免费PDF全文 Jubran Akram 《Geophysical Prospecting》2018,66(5):872-893
Noise suppression or signal‐to‐noise ratio enhancement is often desired for better processing results from a microseismic dataset. In this paper, a polarization–linearity and time–frequency‐thresholding‐based approach is used for denoising waveforms. A polarization–linearity filter is initially applied to preserve the signal intervals and suppress the noise amplitudes. This is followed by time–frequency thresholding for further signal‐to‐noise ratio enhancement in the S transform domain. The parameterisation for both polarization filter and time–frequency thresholding is also discussed. Finally, real microseismic data examples are shown to demonstrate the improvements in processing results when denoised waveforms are considered in the workflow. The results indicate that current denoising approach effectively suppresses the background noise and preserves the vector fidelity of signal waveform. Consequently, the quality of event detection, arrival‐time picking, and hypocenter location improves. 相似文献
16.
In downhole microseismic monitoring, accurate event location relies on the accuracy of the velocity model. The model can be estimated along with event locations. Anisotropic models are important to get accurate event locations. Taking anisotropy into account makes it possible to use additional data – two S-wave arrivals generated due to shear-wave splitting. However, anisotropic ray tracing requires iterative procedures for computing group velocities, which may become unstable around caustics. As a result, anisotropic kinematic inversion may become time consuming. In this paper, we explore the idea of using simplified ray tracing to locate events and estimate medium parameters. In the simplified ray-tracing algorithm, the group velocity is assumed to be equal to phase velocity in both magnitude and direction. This assumption makes the ray-tracing algorithm five times faster compared to ray tracing based on exact equations. We present a set of tests showing that given perforation-shot data, one can use inversion based on simplified ray-tracing even for moderate-to-strong anisotropic models. When there are no perforation shots, event-location errors may become too large for moderately anisotropic media. 相似文献
17.
球坐标系下多震相走时三参数(速度、震源位置和反射界面)同时反演需要解决两个关键问题:(1)球坐标系下3D速度模型中多次透射、反射(折射)及转换波精确、快速的射线追踪;(2)同时反演时三种不同参数间的强耦合问题.为此,我们将直角坐标系下分区多步不规则最短路径算法推广至球坐标系中,进行区域或者全球尺度的多震相射线追踪.然后将其与适合多参数同时反演的子空间算法相结合,形成一种球坐标系下联合多震相走时三参数同时反演的方法技术.与双参数(速度和反射界面或速度和震源位置)同时反演的数值模拟对比分析显示:三参数与双参数的同时反演结果大体接近,并且它们对到时数据中可容许的随机噪声不太敏感.结果说明本文中的同时反演成像为一种提高成像分辨率,同时反演速度、震源位置和反射界面的有效方法. 相似文献
18.
基于偏移成像的微地震定位方法由于可避免走时拾取误差以及可实现自动化定位等优点被广泛应用.绕射叠加方法将能量沿走时曲线聚焦到空间网格点上,进行成像时需要搜索发震时刻.干涉成像方法利用互相关提取的走时差信息可避免搜索发震时刻,但定位结果受数据信噪比的影响较大.为了进一步提高干涉成像法的定位精度,本文提出一种同时使用互相关和自相关道集的全干涉成像方法,增加的自相关道集提取的S-P走时差可降低震源-检波器方向的定位误差,提高定位精度.单井监测的理论测试和实际数据的研究结果表明,全干涉成像的定位精度高于仅使用互相关道集的干涉成像方法,同时计算效率高于绕射叠加方法. 相似文献
19.
本文提出了分析区域地震台网地震定位能力的一种通用方法.在分析过程中考虑了如下事实:当改变地震震级和震源在台网区域内的位置时,能读出 P 震相和 S 震相到时的台站组合也随之改变.作为一个实例,具体分析了北京遥测地震台网的定位能力.北京台网现有19个台站,覆盖面积约为300公里400公里.根据观测资料,建立了确定每个台站检测能力(作为台站仪器放大倍数和地震震级的函数)的经验公式.然后,对每个给定震级和震源位置的可能事件,确定其子台网的构成.在这些子台网的基础上,利用奇异值分解技术,对事先设定的到时数据标准误差计算出震源座标的误差分布和线性化条件方程组的条件数分布.对给定震级(例如,ML=1.0,2.0或3.0),将计算结果绘成等值线图.为了进行比较,还分析了一个包括61个台站的北京台网扩建方案,从而可以预先估计未来台网的地震定位能力. 相似文献
20.
The stimulation of a geothermal well in Basel, Switzerland produced a distribution of microseismic event locations with an overall alignment in the direction of the maximum horizontal stress. Fault plane solutions of individual larger events indicated movements on fracture planes at an angle to the maximum horizontal stress that could not be reliably interpreted from the event locations. To obtain higher resolution images of the microseismic event locations, events with similar waveforms have been identified by multiplet analysis. A number of receivers were used in the multiplet processing to ensure each multiplet is represented by a unique group of waveforms. The location accuracy within each multiplet has been significantly improved using cross‐correlation to refine the shear‐wave traveltime picks. The distribution of events within each multiplet can be interpreted as being due to movements on a single fracture or a number of near parallel fractures. It is shown that whilst the overall distribution of events is around the direction of the maximum horizontal stress, the individual multiplets representing fracture planes have a variety of azimuths and dips. 相似文献