首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Microseismic monitoring has proven invaluable for optimizing hydraulic fracturing stimulations and monitoring reservoir changes. The signal to noise ratio of the recorded microseismic data varies enormously from one dataset to another, and it can often be very low, especially for surface monitoring scenarios. Moreover, the data are often contaminated by correlated noises such as borehole waves in the downhole monitoring case. These issues pose a significant challenge for microseismic event detection. In addition, for downhole monitoring, the location of microseismic events relies on the accurate polarization analysis of the often weak P‐wave to determine the event azimuth. Therefore, enhancing the microseismic signal, especially the low signal to noise ratio P‐wave data, has become an important task. In this study, a statistical approach based on the binary hypothesis test is developed to detect the weak events embedded in high noise. The method constructs a vector space, known as the signal subspace, from previously detected events to represent similar, yet significantly variable microseismic signals from specific source regions. Empirical procedures are presented for building the signal subspace from clusters of events. The distribution of the detection statistics is analysed to determine the parameters of the subspace detector including the signal subspace dimension and detection threshold. The effect of correlated noise is corrected in the statistical analysis. The subspace design and detection approach is illustrated on a dual‐array hydrofracture monitoring dataset. The comparison between the subspace approach, array correlation method, and array short‐time average/long‐time average detector is performed on the data from the far monitoring well. It is shown that, at the same expected false alarm rate, the subspace detector gives fewer false alarms than the array short‐time average/long‐time average detector and more event detections than the array correlation detector. The additionally detected events from the subspace detector are further validated using the data from the nearby monitoring well. The comparison demonstrates the potential benefit of using the subspace approach to improve the microseismic viewing distance. Following event detection, a novel method based on subspace projection is proposed to enhance weak microseismic signals. Examples on field data are presented, indicating the effectiveness of this subspace‐projection‐based signal enhancement procedure.  相似文献   

2.
Comparison of surface and borehole locations of induced seismicity   总被引:1,自引:0,他引:1  
Monitoring of induced microseismic events has become an important tool in hydraulic fracture diagnostics and understanding fractured reservoirs in general. We compare microseismic event and their uncertainties using data sets obtained with surface and downhole arrays of receivers. We first model the uncertainties to understand the effect of different acquisition geometries on location accuracy. For a vertical array of receivers in a single monitoring borehole, we find that the largest part of the final location uncertainty is related to estimation of the backazimuth. This is followed by uncertainty in the vertical position and radial distance from the receivers. For surface monitoring, the largest uncertainty lies in the vertical position due to the use of only a single phase (usually P‐wave) in the estimation of the event location. In surface monitoring results, lateral positions are estimated robustly and are not sensitive to the velocity model. In this case study, we compare event location solutions from two catalogues of microseismic events; one from a downhole array and the second from a surface array of 1C geophone. Our results show that origin time can be reliably used to find matching events between the downhole and surface catalogues. The locations of the corresponding events display a systematic shift consistent with a poorly calibrated velocity model for downhole dataset. For this case study, locations derived from surface monitoring have less scatter in both vertical and horizontal directions.  相似文献   

3.
We present a novel denoising scheme via Radon transform-based adaptive vector directional median filters named adaptive directional vector median filter (AD-VMF) to suppress noise for microseismic downhole dataset. AD-VMF contains three major steps for microseismic downhole data processing: (i) applying Radon transform on the microseismic data to obtain the parameters of the waves, (ii) performing S-transform to determine the parameters for filters, and (iii) applying the parameters for vector median filter (VMF) to denoise the data. The steps (i) and (ii) can realize the automatic direction detection. The proposed algorithm is tested with synthetic and field datasets that were recorded with a vertical array of receivers. The P-wave and S-wave direct arrivals are properly denoised for poor signal-to-noise ratio (SNR) records. In the simulation case, we also evaluate the performance with mean square error (MSE) in terms of signal-to-noise ratio (SNR). The result shows that the distortion of the proposed method is very low; the SNR is even less than 0 dB.  相似文献   

4.
微地震事件初至拾取是井下微地震监测数据处理的关键步骤之一.初至误差的存在会使微地震震源定位结果产生较大偏差,进而影响后续的压裂裂缝解释.通常初至拾取过程对所有的微地震事件选择相同的特征函数并采用一致的拾取参数进行统一处理,然而当事件的能量、震源机制、传播路径以及背景噪声等存在明显差异时,所得初至拾取结果差别显著.为了提高微地震事件初至拾取标准一致性,本文提出基于波形相似特征的初至拾取及全局校正方法.该方法首先利用互相关函数对每个事件内的各道记录进行时差校正,得到初始初至信息并形成叠加道,再对所有事件的叠加道进行全局互相关得到事件间初至相对校正量,最终初至结果可以通过各个事件的初始初至信息与其相对校正量相加得到.方法将所有微地震事件初至结果作为一个整体处理,从而能够克服常规方法初至拾取标准一致性差的缺陷.实际资料处理结果表明,相比于常规方法,该方法可以有效提高事件初至拾取和定位结果的一致性.  相似文献   

5.
Distributed acoustic sensing uses an optical fibre together with an interrogator unit to perform strain measurements. The usage of distributed acoustic sensing in geophysics is attractive due to its dense spatial sampling and low operation cost if the optical fibre is freely accessible. In the borehole environment, optical fibres for distributed acoustic sensing are often readily available as a part of other sensing tools, such as for temperature and pressure. Although the distributed acoustic sensing system promises great potential for reservoir monitoring and surface seismic acquisition, the single axial strain measurement of distributed acoustic sensing along the fibre is inadequate to fully characterise the different wave modes, thus making reservoir characterisation challenging. We propose an acquisition system using five equally spaced helical optical fibres and a straight optical fibre to obtain six different strain projections. This system allows us to reconstruct all components of the 3D strain tensor at any location along the fibre. Analysing the condition number associated with the geometry of the optical fibre, we can systematically search for the optimum design parameters for our configuration. Numerical examples demonstrate the effectiveness of our proposed method to successful reconstruction of the full strain tensor from elastic wavefields of arbitrary complexity.  相似文献   

6.
A modular borehole monitoring concept has been implemented to provide a suite of well‐based monitoring tools that can be deployed cost effectively in a flexible and robust package. The initial modular borehole monitoring system was deployed as part of a CO2 injection test operated by the Southeast Regional Carbon Sequestration Partnership near Citronelle, Alabama. The Citronelle modular monitoring system transmits electrical power and signals, fibre‐optic light pulses, and fluids between the surface and a reservoir. Additionally, a separate multi‐conductor tubing‐encapsulated line was used for borehole geophones, including a specialized clamp for casing clamping with tubing deployment. The deployment of geophones and fibre‐optic cables allowed comparison testing of distributed acoustic sensing. We designed a large source effort (>64 sweeps per source point) to test fibre‐optic vertical seismic profile and acquired data in 2013. The native measurement in the specific distributed acoustic sensing unit used (an iDAS from Silixa Ltd) is described as a localized strain rate. Following a processing flow of adaptive noise reduction and rebalancing the signal to dimensionless strain, improvement from repeated stacking of the source was observed. Conversion of the rebalanced strain signal to equivalent velocity units, via a scaling by local apparent velocity, allows quantitative comparison of distributed acoustic sensing and geophone data in units of velocity. We see a very good match of uncorrelated time series in both amplitude and phase, demonstrating that velocity‐converted distributed acoustic sensing data can be analyzed equivalent to vertical geophones. We show that distributed acoustic sensing data, when averaged over an interval comparable to typical geophone spacing, can obtain signal‐to‐noise ratios of 18 dB to 24 dB below clamped geophones, a result that is variable with noise spectral amplitude because the noise characteristics are not identical. With vertical seismic profile processing, we demonstrate the effectiveness of downgoing deconvolution from the large spatial sampling of distributed acoustic sensing data, along with improved upgoing reflection quality. We conclude that the extra source effort currently needed for tubing‐deployed distributed acoustic sensing vertical seismic profile, as part of a modular monitoring system, is well compensated by the extra spatial sampling and lower deployment cost as compared with conventional borehole geophones.  相似文献   

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.
Like most other industrial activities that affect the subsurface, hydraulic fracturing carries the risk of reactivating pre‐existing faults and thereby causing induced seismicity. In some regions, regulators have responded to this risk by imposing traffic light scheme‐type regulations, where fracture stimulation programs must be amended or shut down if events larger than a given magnitude are induced. Some sites may be monitored with downhole arrays and/or dense near‐surface arrays, capable of detecting very small microseismic events. However, such monitoring arrangements will not be logistically or economically feasible at all sites. Instead, operators are using small, sparse arrays of surface seismometers to meet their monitoring obligations. The challenge we address in this paper is to maximise the detection thresholds of such small, sparse, surface arrays so that they are capable of robustly identifying small‐magnitude events whose signal‐to‐noise ratios may be close to 1. To do this, we develop a beamforming‐and‐stacking approach, computing running short‐term/long‐term average functions for each component of each recorded trace (P, SH, and SV), time‐shifting these functions by the expected travel times for a given location, and performing a stack. We assess the effectiveness of this approach with a case study using data from a small surface array that recorded a multi‐well, multi‐stage hydraulic fracture stimulation in Oklahoma over a period of 8 days. As a comparison, we initially used a conventional event‐detection algorithm to identify events, finding a total of 17 events. In contrast, the beamforming‐and‐stacking approach identified a total of 155 events during this period (including the 17 events detected by the conventional method). The events that were not detected by the conventional algorithm had low‐signal‐to‐noise ratios to the extent that, in some cases, they would be unlikely to be identified even by manual analysis of the seismograms. We conclude that this approach is capable of improving the detection thresholds of small, sparse arrays and thus can be used to maximise the information generated when deployed to monitor industrial sites.  相似文献   

9.
Testing the ability of surface arrays to monitor microseismic activity   总被引:2,自引:0,他引:2  
Recently there has been much interest in the use of data from surface arrays in conjunction with migration‐based processing methods for passive seismic monitoring. In this study we use an example of this kind of data recorded whilst 18 perforation shots, with a variety of positions and propellant amounts, were detonated in the subsurface. As the perforation shots provide signals with known source positions and origin times, the analysis of these data is an invaluable opportunity to test the accuracy and ability of surface arrays to detect and locate seismic sources in the subsurface. In all but one case the signals from the perforation shots are not visible in the raw or preprocessed data. However, clear source images are produced for 12 of the perforation shots showing that arrays of surface sensors are capable of imaging microseismic events, even when the signals are not visible in individual traces. We find that point source locations are within typically 45 m (laterally) of the true shot location, however the depths are less well constrained (~150 m). We test the sensitivity of our imaging method to the signal‐to‐noise ratio in the data using signals embedded in realistic noise. We find that the position of the imaged shot location is quite insensitive to the level of added noise, the primary effect of increased noise being to defocus the source image. Given the migration approach, the array geometry and the nature of coherent noise during the experiment, signals embedded in noise with ratios ≥0.1 can be used to successfully image events. Furthermore, comparison of results from data and synthetic signals embedded in noise shows that, in this case, prestack corrections of traveltimes to account for near‐surface structure will not enhance event detectability. Although, the perforation shots have a largely isotropic radiation pattern the results presented here show the potential for the use of surface sensors in microseismic monitoring as a viable alternative to classical downhole methods.  相似文献   

10.
Seismic anisotropy which is common in shale and fractured rocks will cause travel-time and amplitude discrepancy in different propagation directions. For microseismic monitoring which is often implemented in shale or fractured rocks, seismic anisotropy needs to be carefully accounted for in source location and mechanism determination. We have developed an efficient finite-difference full waveform modeling tool with an arbitrary moment tensor source. The modeling tool is suitable for simulating wave propagation in anisotropic media for microseismic monitoring. As both dislocation and non-double-couple source are often observed in microseismic monitoring, an arbitrary moment tensor source is implemented in our forward modeling tool. The increments of shear stress are equally distributed on the staggered grid to implement an accurate and symmetric moment tensor source. Our modeling tool provides an efficient way to obtain the Green’s function in anisotropic media, which is the key of anisotropic moment tensor inversion and source mechanism characterization in microseismic monitoring. In our research, wavefields in anisotropic media have been carefully simulated and analyzed in both surface array and downhole array. The variation characteristics of travel-time and amplitude of direct P- and S-wave in vertical transverse isotropic media and horizontal transverse isotropic media are distinct, thus providing a feasible way to distinguish and identify the anisotropic type of the subsurface. Analyzing the travel-times and amplitudes of the microseismic data is a feasible way to estimate the orientation and density of the induced cracks in hydraulic fracturing. Our anisotropic modeling tool can be used to generate and analyze microseismic full wavefield with full moment tensor source in anisotropic media, which can help promote the anisotropic interpretation and inversion of field data.  相似文献   

11.
The hydrocarbon industry is moving increasingly towards tight sandstone and shale gas resources – reservoirs that require fractures to be produced economically. Therefore, techniques that can identify sets of aligned fractures are becoming more important. Fracture identification is also important in the areas of coal bed methane production, carbon capture and storage (CCS), geothermal energy, nuclear waste storage and mining. In all these settings, stress and pore pressure changes induced by engineering activity can generate or reactivate faults and fractures. P‐ and S‐waves are emitted by such microseismic events, which can be recorded on downhole geophones. The presence of aligned fracture sets generates seismic anisotropy, which can be identified by measuring the splitting of the S‐waves emitted by microseismic events. The raypaths of the S‐waves will have an arbitrary orientation, controlled by the event and geophone locations, meaning that the anisotropy system may only be partly illuminated by the available arrivals. Therefore to reliably interpret such splitting measurements it is necessary to construct models that compare splitting observations with modelled values, allowing the best fitting rock physics parameters to be determined. Commonly, splitting measurements are inverted for one fracture set and rock fabrics with a vertical axis of symmetry. In this paper we address the challenge of identifying multiple aligned fracture sets using splitting measured on microseismic events. We analyse data from the Weyburn CCS‐EOR reservoir, which is known to have multiple fracture sets, and from a hydraulic fracture stimulation, where it is believed that only one set is present. We make splitting measurements on microseismic data recorded on downhole geophone arrays. Our inversion technique successfully discriminates between the single and multiple fracture cases and in all cases accurately identifies the strikes of fracture sets previously imaged using independent methods (borehole image logs, core samples, microseismic event locations). We also generate a synthetic example to highlight the pitfalls that can be encountered if it is assumed that only one fracture set is present when splitting data are interpreted, when in fact more than one fracture set is contributing to the anisotropy.  相似文献   

12.
李稳  刘伊克  刘保金 《地球物理学报》2016,59(10):3869-3882
井下微震监测获得的地震记录往往包含大量的噪声,记录信噪比很低.有效地震信号的识别与提取是进行后续地震定位等工作之前需要优先解决的问题.经过研究发现,井下水压裂微地震信号具有稀疏分布的特征,而井下环境噪声则具有更多的Gaussian分布特征.为此,本文提出将图像处理领域适宜于稀疏分布信号降噪处理的稀疏码收缩方法应用于井下微震监测数据处理.为解决需要利用与待处理数据中有效信号成分具有相似分布特征的无噪信号序列估算正交基以及计算效率等问题,将原方法与小波变换理论相结合.即通过优选小波基函数作为正交基进行小波变换将信号分解为不同级的小波系数,利用稀疏码收缩方法中对稀疏编码施加的非线性收缩方式作为阈值准则对小波系数进行改造.通过多方面的数值实验证明了该方法在处理地震子波及井下微地震信号方面准确可靠.含噪记录经过处理后有效地震信号的到时、波形、时频谱特征等均能得到良好的识别和恢复.并且该方法具有很强的抗噪能力,当信噪比低至-20~-30db时,仍然能够发挥作用.在处理大量实际井下微震监测数据的过程中,面对多种复杂情况,本方法展现出了计算效率高、计算结果可靠、应用简单等优势,证明了其本身具有实际应用价值,值得进一步的研究和推广.  相似文献   

13.
Microseismic monitoring is an approach for mapping hydraulic fracturing. Detecting the accurate locations of microseismic events relies on an accurate velocity model. The one‐dimensional layered velocity model is generally obtained by model calibration from inverting perforation data. However, perforation shots may only illuminate the layers between the perforation shots and the recording receivers with limited raypath coverage in a downhole monitoring problem. Some of the microseismic events may occur outside of the depth range of these layers. To derive an accurate velocity model covering all of the microseismic events and locating events at the same time, we apply the cross double‐difference method for the simultaneous inversion of a velocity model and event locations using both perforation shots and microseismic data. The cross double‐difference method could provide accurate locations in both the relative and absolute sense, utilizing cross traveltime differences between P and S phases over different events. At the downhole monitoring scale, the number of cross traveltime differences is sufficiently large to constrain events locations and velocity model as well. In this study, we assume that the layer thickness is known, and velocities of P‐ and S‐wave are inverted. Different simultaneous inversion methods based on the Geiger's, double‐difference, and cross double‐difference algorithms have been compared with the same input data. Synthetic and field data experiments suggest that combining both perforation shots and microseismic data for the simultaneous cross double‐difference inversion of the velocity model and event locations is available for overcoming the trade‐offs in solutions and producing reliable results.  相似文献   

14.
地震检测与震相自动拾取研究   总被引:3,自引:2,他引:1       下载免费PDF全文
针对微震事件易受噪声干扰等特点,本文将STA/LTA方法和基于方差的AIC方法(var-AIC)相结合,在震相到时初步拾取的基础上,使用台站的德洛内(Delaunay)三角剖分及台站间最大走时差约束来减少噪声干扰的影响. 利用到时进行地震定位之后,根据台站预测到时,在设定的时间窗内对地震震相进行更精细的分析. 特别是针对微震事件信噪比低的特点,设计了基于偏振分析的拾取函数,根据窗内STA/LTA方法和var-AIC方法的拾取结果自动选择合适的值作为震相到时. 最后,对西昌流动地震台阵2013年304个单事件波形数据的分析处理和检验结果表明,本文方法较传统方法具有更高的地震事件检测能力和更高的震相拾取精度.   相似文献   

15.
Microseismic monitoring in the oil and gas industry commonly uses migration‐based methods to locate very weak microseismic events. The objective of this study is to compare the most popular migration‐based methods on a synthetic dataset that simulates a strike‐slip source mechanism event with a low signal‐to‐noise ratio recorded by surface receivers (vertical components). The results show the significance of accounting for the known source mechanism in the event detection and location procedures. For detection and location without such a correction, the ability to detect weak events is reduced. We show both numerically and theoretically that neglecting the source mechanism by using only absolute values of the amplitudes reduces noise suppression during stacking and, consequently, limits the possibility to retrieve weak microseismic events. On the other hand, even a simple correction to the data polarization used with otherwise ineffective methods can significantly improve detections and locations. A simple stacking of the data with a polarization correction provided clear event detection and location, but even better results were obtained for those data combined with methods that are based on semblance and cross‐correlation.  相似文献   

16.
For years, severe rockburst problems at the Lucky Friday mine in northern Idaho have been a persistent safety hazard and an impediment to production. An MP250 based microseismic monitoring system, which uses simple voltage threshold picking of first arrivals, has been used in this mine since 1973 to provide source locations and energy estimates of seismic events. Recently, interest has been expressed in developing a whole waveform microseismic monitoring system for the mine to provide more accurate source locations and information about source characteristics. For this study, we have developed a prototype whole-waveform microseismic monitoring system based on a 80386 computer equipped with a 50 kHz analog-digital convertor board. The software developed includes a data collection program, a data analysis program, and an event detection program. Whole-waveform data collected and analyzed using this system during a three-day test have been employed to investigate sources of error in the hypocenter location process and to develop an automatic phase picker appropriate for microseismic events.Comparison of hypocenter estimates produced by the MP250 system to those produced by the whole-waveform system shows that significant timing errors are common in the MP250 system and that these errors caused a large part of the scatter evident in the daily activity plots produced at the mine. Simulations and analysis of blast data show that analytical control over the solutions is strongly influenced by the array geometry. Within the geophone array, large errors in the velocity model or moderate timing errors may result in small changes in the solution, but outside the array, the solution is very sensitive to small changes in the data.Our whole-waveform detection program picks event onset times and determines event durations by analysis of a segmented envelope function (SEF) derived from the microseismic signal. The detection program has been tested by comparing its arrival time picks to those generated by human analysis of the data set. The program picked 87% of the channels that were picked by hand with a standard error of 0.75 milliseconds. Source locations calculated using times provided by our entire waveform detection program were similar to those calculated using hand-picked arrival times. In particular, they show far less scatter than source locations calculated using arrival times based on simple voltage threshold picking of first arrivals.  相似文献   

17.
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.  相似文献   

18.
We present results of processed microseismic events induced by hydraulic fracturing and detected using dual downhole monitoring arrays. The results provide valuable insight into hydraulic fracturing. For our study, we detected and located microseismic events and determined their magnitudes, source mechanisms and inverted stress field orientation. Event locations formed a distinct linear trend above the stimulated intervals. Source mechanisms were only computed for high‐quality events detected on a sufficient number of receivers. All the detected source mechanisms were dip‐slip mechanisms with steep and nearly horizontal nodal planes. The source mechanisms represented shear events and the non‐double‐couple components were very small. Such small, non‐double‐couple components are consistent with a noise level in the data and velocity model uncertainties. Strikes of inverted mechanisms corresponding to the nearly vertical fault plane are (within the error of measurements) identical with the strike of the location trend. Ambient principal stress directions were inverted from the source mechanisms. The least principal stress, σ3, was determined perpendicular to the strike of the trend of the locations, indicating that the hydraulic fracture propagated in the direction of maximum horizontal stress. Our analysis indicated that the source mechanisms observed using downhole instruments are consistent with the source mechanisms observed in microseismic monitoring arrays in other locations. Furthermore, the orientation of the inverted principal components of the ambient stress field is in agreement with the orientation of the known regional stress, implying that microseismic events induced by hydraulic fracturing are controlled by the regional stress field.  相似文献   

19.
Seismic site response analysis is commonly used to predict ground response due to local soil effects. An increasing number of downhole arrays are deployed to measure motions at the ground surface and within the soil profile and to provide a check on the accuracy of site response analysis models. Site response analysis models, however, cannot be readily calibrated to match field measurements. A novel inverse analysis framework, self-learning simulations (SelfSim), to integrate site response analysis and field measurements is introduced. This framework uses downhole array measurements to extract the underlying soil behavior and develops a neural network-based constitutive model of the soil. The resulting soil model, used in a site response analysis, provides correct ground response. The extracted cyclic soil behavior can be further enhanced using multiple earthquake events. The performance of the algorithm is successfully demonstrated using synthetically generated downhole array recordings.  相似文献   

20.
微震监测是直观评价压裂过程和压裂效果的有效手段.微震事件识别是微震监测的首要步骤.然而对于低信噪比微震监测数据,常规识别方法很难取得满意效果.基于微震事件在时频域中的稀疏性,本文提出利用Renyi熵值表示微震监测数据的时频稀疏程度,并以时频距离为约束条件,建立以低熵值的道数为判别阈值的目标函数.本文方法能在识别出微震事件的同时,恢复出较为清晰的微震事件.通过数值计算和对实际监测数据的测试,表明该方法对低信噪比的微震监测数据有较好的处理效果.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号