BSN矿山微震监测技术及其应用

深井高应力、高岩压诱发的冒顶、片帮、塌方、岩爆等地压灾害严重影响了矿山安全生产。随着高精度、高采样率传感器的研发和计算机技术的飞速发展,微震高精度实时监测技术不断完善,微震监测已成为矿山地压灾害监测预警的重要手段。通过设计、研发完全自主知识产权的BSN(BGRIMM Seismic Network)微震监测系统,实现了微震信号数据采集、多通道时钟同步、噪声压制、到时自动拾取、震源定位等功能,并结合统计地震学的相关方法,初步实现了基于微震活动性分析的岩体危险性评价。本文主要从矿山微震监测数据采集、预处理、定位、活动性分析等模块介绍BSN微震监测系统全流程,并对监测系统在广西某矿山的实际应用案例进行分析。     

三维地震与地面微地震联合校正方法

由于地面微地震监测台站布设在地表,会受到地表起伏、低降速带厚度和速度变化的影响,降低了微地震事件的识别准确度和定位精度,限制了地面微地震监测技术在复杂地表地区的应用.因此,将三维地震勘探技术的思路引入到地面微地震监测中,提出了三维地震与地面微地震联合校正方法,将油气勘探和开发技术更加紧密地结合在一起.根据三维地震数据和低降速带测量数据,通过约束层析反演方法建立精确的近地表速度模型,将地面微地震台站从起伏地表校正到高速层中的平滑基准面上,有效消除复杂近地表的影响.其次,根据射孔数据和声波测井速度信息,通过非线性反演方法建立最优速度模型,由于已经消除复杂近地表的影响,在进行速度模型优化时不需要考虑近地表的影响,因而建立的速度模型更加准确.最后,在精确速度模型的基础上,通过互相关方法求取剩余静校正量,进一步消除了复杂近地表和速度模型近似误差的影响.三维地震与地面微地震联合校正方法采用逐步校正的思路,能够有效消除复杂近地表的影响,提高微地震数据的品质和速度模型的精确度,保证了微地震事件的定位精度,具有良好的应用前景.     

Passive Seismic Monitoring of Mine-scale Geothermal Activity: A Trial at Lihir Open Pit Mine

The Lihir open pit mine in Papua New Guinea is located inside an old volcano where geothermal activity is strongly present. Outbursts of hot water and steam into the mining areas were a major safety concern. Passive seismic monitoring was carried out at the mine to investigate whether the geothermal activities could be detected and located using microseismic techniques in a mining environment. In this trial, sixteen triaxial geophones which can withstand temperature up to 200°C were used and installed in four deep boreholes inside the pit. The microseismic events were discriminated using the STA/LTA triggering criterion. During 6 weeks of monitoring, more than 17,000 events were recorded. Approximately 12% of the events showed harmonic vibration characteristics similar to those observed in other geothermal and volcanic areas, suggesting that the geothermal activity inside the pit was captured by the microseismic monitoring system. More than 75% of the events present both P and S waves and they were interpreted to be associated with rock fracturing due to stress release near the bottom of the pit. Many geothermal-type events were located in areas where shear events occurred, implying that the detected geothermal events were not far from the mining area below the pit and they may also be associated with mining. The borehole installation of the geophones significantly reduced the interference of mining noise and achieved good observation of the seismic events. However, equipment installation requires great attention as the geophones may be destroyed due to unexpected rising temperature within the boreholes.     

Application of Kalman Filtering Techniques for Microseismic Event Detection

—?Microseismic monitoring systems are generally installed in areas of induced seismicity caused by human activity. Induced seismicity results from changes in the state of stress which may occur as a result of excavation within the rock mass in mining (i.e., rockbursts), and changes in hydrostatic pressures and rock temperatures (e.g., during fluid injection or extraction) in oil exploitation, dam construction or fluid disposal. Microseismic monitoring systems determine event locations and important source parameters such as attenuation, seismic moment, source radius, static stress drop, peak particle velocity and seismic energy. An essential part of the operation of a microseismic monitoring system is the reliable detection of microseismic events. In the absence of reliable, automated picking techniques, operators rely upon manual picking. This is time-consuming, costly and, in the presence of background noise, very prone to error. The techniques described in this paper not only permit the reliable identification of events in cluttered signal environments they have also enabled the authors to develop reliable automated event picking procedures. This opens the way to use microseismic monitoring as a cost-effective production/operations procedure. It has been the experience of the authors that in certain noisy environments, the seismic monitoring system may trigger on and subsequently acquire substantial quantities of erroneous data, due to the high energy content of the ambient noise. Digital filtering techniques need to be applied on the microseismic data so that the ambient noise is removed and event detection simplified. The monitoring of seismic acoustic emissions is a continuous, real-time process and it is desirable to implement digital filters which can also be designed in the time domain and in real-time such as the Kalman Filter. This paper presents a real-time Kalman Filter which removes the statistically describable background noise from the recorded seismic traces.     

Management of Post-mining Large-scale Ground Failures: Blast Swarms Field Experiment for Calibration of Permanent Microseismic Early-warning Systems

In France, decades of coal and iron-ore mining have left extensive underground cavities beneath or in the vicinity of urban areas. This poses an environmental challenge for society. To ensure post-mining risk management and public safety, wherever remediation is not possible, numerous real-time microseismic monitoring systems are being installed. The objective is to detect remote rock mass fracturing processes, precursory events and acceleration phases for appropriate and timely action. Although no consistent collapse has occurred in any of the monitored areas yet, single 3-D probes record many microseismic events of very low amplitude which create difficulties in the quantitative data analysis. The development of specific quantitative processing has therefore become a major issue in our research work. For that purpose, a field experiment was carried out on six of the instrumented sites. It consisted of sequences of small blasts in mine pillars which were accurately controlled in terms of the location, orientation and energy of the explosive source. The data analysis was used to calibrate parameters (velocity model, 3-D sensor orientation, etc.) for reliable 3-D localization and to develop an empirical law to estimate the source energy from the sensor energy. This work now enables us to analyze real microseismic events with a considerably better level of accuracy and to obtain enough information and confidence to discuss these data in terms of site stability.     

Microseismic Full Waveform Modeling in Anisotropic Media with Moment Tensor Implementation

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.     

Tomographic inversion for microseismic source parameters in mining

We propose a new method for inverting source function of microseismic event induced in mining. The observed data from microseismic monitoring during mining are represented by a wave equation in a spherical coordinate system and then the data are transformed from the time-space domain to the time-slowness domain based on tomographic principle, from whichwe can obtain the signals related to the source in the time-slowness domain. Through analyzing the relationship between the signal located at the maximum energy and the source function, we derive the tomographic equations to compute the source function from the signals and to calculate the effective radiated energy based on the source function. Moreover, we fit the real amplitude spectrum of the source function computed from the observed data into the ?? ^?2 model based on the least squares principle and determine the zero-frequency level spectrum and the corner frequency, finally, the source rupture radius of the event is calculated and The synthetic and field examples demonstrate that the proposed tomographic inversion methods are reliable and efficient     

Results from microseismic monitoring,conventional instrumentation,and tomography surveys in the creation and thinning of a burst-prone sill pillar

Located in northern Québec, the Lac Shortt Mine was a small gold mine consisting of a thin subvertical orebody which was mined in three main phases. High stress and rockbursting conditions were experienced when ore was extracted in the upper zone between the surface and a depth of 500 metres during the first two phases of mining. Severe rockbursts were experienced in late 1989 near the shaft and in the footwall development following a deepening of the mine shaft to a depth of 830 m and partial development of footwall drift access for the third phase of mining (the mining of the lower zone starting at a depth of 830 m moving upward toward a depth of 500 m). A 16-channel Electrolab MP250 microseismic system, with a Queen's University Full-Waveform piggy-back system, was installed underground at the site due to these problems.It was expected that the thinning sill would be subjected to an ever-increasing load as the thickness of the 500 m sill pillar decreased in the face of the mining excavation from below. A monitoring program consisting of the microseismic monitoring system, a range of conventional geomechanics monitoring tools as well as the undertaking of periodic seismic tomography surveys to assess the ongoing state of stress and rock mass condition within the sill was therefore warranted.The anomalously high-magnitude stress field and the brittle rockmass created a situation in which rockmass failure was common and violent. In the creation and thinning of the sill pillar, the location of banded microseismic activity was crucial in tracing rockmass failure and the associated ground control problems. Reliable source-location determination enabled the identification of areas of stress increase. The movement of the rockmass failure front could be followed, and was responsible for stope dilution, footwall and orebody development deterioration, and caving.Source-mechanism analyses gave accurate double-couple solutions for approximately forty percent of these events having at least ten recognizable polarities. Results suggested movement along vertical north-south striking or vertical east-west striking features. Underground observation of damaged access points showed that vertical north-south striking joints were experiencing failure.The microseismic activity, which was consistently concentrated close to the southwest and northeast corners of current production stopes, could be explained by a stress field oriented obliquely to the strike of the orebody, as measured prior to shrinkage of the sill pillar byin situ stress measurements and observed borehole overbreaks. The orientations of theP andT axes for the microseismic activity further confirmed that the stress field oriented obliquely to strike.While an increase in compressional-wave velocity of 2.3 percent, corresponding to a measured stress increase of approximately 10 MPa could be measured by repeated tomographic surveys, it was relatively small and only a factor of two or so above the velocity measured uncertainty. The relative insensitivity of thein situ rock mass modulus to the applied stress is believed to be largely due to the rockmass discontinuities being relatively closed prior to stress increase, as substantiated by the small deformations seen by the extensometer and borehole camera. This situation existed because of the very high pre-mining stress level.The experimental demonstration that the rock could not absorb substantially increased load through the mechanism of discontinuity closure or tightening (which would be reflected in the modulus) may be evidence in itself of potentially burst-prone ground, such as encountered at Lac Shortt.     

基于向量扫描法的微地震震级确定方法

对地震记录中显著可视事件震级的确定,已有如里氏震级、矩震级、不同定义的震级之间的关联、震级与能量之间的关系等成熟方法。然而,由于微地震事件记录常淹没在背景噪声中,无法实施包括确定震级在内的传统定位。为比较不同类型不同地域的微震大小、研究微震机制及诱发原因、使用其安全预警等,我们外延里氏震级到微震范畴,定义了等效里氏微震震级:根据微破裂向量扫描输出的无量纲能量或最小信噪比,估计一定时空中单位时间单位台站所接收的在时空目标一点上的等效能量和震级。文中列出了目前我们在微震监测中发现的几类微震的大小规模。      

A strategy for automated analysis of passive microseismic data to image seismic anisotropy and fracture characteristics

Monitoring of induced seismicity is gaining importance in a broad range of industrial operations from hydrocarbon reservoirs to mining to geothermal fields. Such passive seismic monitoring mainly aims at identifying fractures, which is of special interest for safety and productivity reasons. By analysing shear‐wave splitting it is possible to determine the anisotropy of the rock, which may be caused by sedimentary layering and/or aligned fractures, which in turn offers insight into the state of stress in the reservoir. We present a workflow strategy for automatic and effective processing of passive microseismic data sets, which are ever increasing in size. The automation provides an objective quality control of the shear‐wave splitting measurements and is based on characteristic differences between the two independent eigenvalue and cross‐correlation splitting techniques. These differences are summarized in a quality index for each measurement, allowing identification of an appropriate quality threshold. Measurements above this threshold are considered to be of good quality and are used in further interpretation. We suggest an automated inversion scheme using rock physics theory to test for best correlation of the data with various combinations of fracture density, its strike and the background anisotropy. This fully automatic workflow is then tested on a synthetic and a real microseismic data set.     

A review of the microseismic focal mechanism research

Reservoir reconstructions implemented in unconventional oil and gas exploration usually adopt hydraulic fracturing techniques to inject high-pressure fluid into the reservoir and change its pore-fracture connection structure to enhance production. Hydraulic fracturing changes the reservoir stress and causes the rocks to crack, thus generating microseismic events.One important component of microseismic research is the source mechanism inversion. Through the research on the microseismic focal mechanism, information on the source mechanisms and in-situ stress status variations can be quantitatively revealed to effectively optimize the reservoir reconstruction design for increasing production. This paper reviews the recent progress in hydraulic fracturing induced microseismic focal mechanism research. We summarize their main principles and provide a detailed introduction of the research advances in source modeling, microseismic data synthesis, and focal mechanism inversion. We also discuss the challenges and limitations in the current microseismic focal mechanism research and propose prospects for future research ideas and directions.     

基于小波变换的大尺度岩体结构微震监测信号去噪方法研究

为将小波去噪方法应用于大尺度岩体结构微震监测信号的去噪研究,首先在MATLAB环境下进行仿真,验证了使用Symlet6小波进行小波去噪的可行性;利用4种自适应阈值规则对含噪信号进行去噪对比,结果表明4种阈值去噪后的信号在均方差较小的情况下都极大地提高了信号的信噪比,有效地去除了噪声,对不同的含噪信号,无偏似然原则阈值去...     

“交通灯”系统及其在诱发地震监测中的应用

随着页岩气开采、地热能源开采、CO2封存和水库蓄水等工业活动的开展,由此产生的诱发地震危害问题日益严重。为了对诱发地震进行有效监测和管控,地震学家设计了"交通灯"系统,当地震活动达到一定阈值时,相关人员可根据系统警报及时采取应对措施,减轻地震灾害。本文调研了近年来多个国家的"交通灯"系统研究和应用进展,介绍了其发展历史、设置原理及应用案例。"交通灯"系统的发展和完善应综合考虑震级、震动强度、b值和断层分布,以及建筑结构、城镇距离、人口密度、公众反应等社会因素。同时,在诱发地震频发区域,应建设密集的专用监测台网,提高微震检测和定位能力,并引入模板匹配和人工智能等新的自动处理方法,及时产出高完备性、高精度的微震目录,对微震的时空演化进行有效监测,形成时效性较强的"交通灯"系统。     

Cross double‐difference inversion for simultaneous velocity model update and microseismic event location

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.     

The Canada-Ontario-industry rockburst project

This project was started in 1985, after a particularly extensive period of rockburst activity in several Ontario mines. Three parties participated in the project: the Government of Canada, through CANMET, the Ontario Ministries of Labour and Northern Development and Mines, and the Ontario Mining Association.The first two years of the project have been mainly devoted to the design and installation of new seismic monitoring systems. It is intended to install three different types of monitoring systems at the four mining camps experiencing rockbursts (Red Lake, Elliot Lake, Sudbury, and Kirkland Lake). Seismograph units will be installed at each mining camp to obtain permanent records of the larger seismic events and their magnitude. Macroseismic systems are being installed around five mines (Campbell, Quirke, Strathcona, Creighton and Macassa). These systems consist of triaxial, strong-motion geophones with processing units for event detection and data digitization. Complete waveforms are captured to study first motion, peak particle velocity, seismic energy and spectral frequency. At present 13 mines in Ontario operate their own microseismic systems which are used exclusively for real time source location of seismic events.The instrumentation already installed and the present research activities at the mines are described in the paper.Presented at the Fred Leighton Memorial Workshop on Mining Induced Seismicity, Montreal 1987.     

基于稀疏分布特征的井下微地震信号识别与提取方法

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

Geotomographic imaging in the study of mining induced seismicity

Geotomographic imaging is a technique which allows seismic waves to be used to gain information about the internal structure of rock masses, in a way conceptually similar to medical CAT scanning. Traditional approaches to the study of mining-induced seismic phenomena have concentrated on using passive monitoring methods. This paper gives an overview of the developments in the acquisition, processing and interpretation of geotomographic data and outlines how images can be used in conjunction with passive techniques to study mining-induced seismicity.Presented at the Fred Leighton Memorial Workshop on Mining Induced Seismicity, Montreal 1987.     

Comparison of surface and borehole locations of induced seismicity

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.     

基于时频稀疏性分析法的低信噪比微震事件识别与恢复

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

Feasibility of borehole ambient noise interferometry for permanent reservoir monitoring

The analysis of seismic ambient noise acquired during temporary or permanent microseismic monitoring campaigns (e.g., improved/enhanced oil recovery monitoring, surveillance of induced seismicity) is potentially well suited for time‐lapse studies based on seismic interferometry. No additional data acquisition required, ambient noise processing can be automatized to a high degree, and seismic interferometry is very sensitive to small medium changes. Thus there is an opportunity for detection and monitoring of velocity variations in a reservoir at negligible additional cost and effort. Data and results are presented from an ambient noise interferometry study applied to two wells in a producing oil field in Romania. Borehole microseismic monitoring on three component geophones was performed for four weeks, concurrent with a water‐flooding phase for improved oil recovery from a reservoir in ca. 1 km depth. Both low‐frequency (2 Hz–50 Hz) P‐ and S‐waves propagating through the vertical borehole arrays were reconstructed from ambient noise by the virtual source method. The obtained interferograms clearly indicate an origin of the ambient seismic energy from above the arrays, thus suggesting surface activities as sources. It is shown that ambient noise from time periods as short as 30 seconds is sufficient to obtain robust interferograms. Sonic log data confirm that the vertical and horizontal components comprise first arrivals of P‐wave and S‐waves, respectively. The consistency and high quality of the interferograms throughout the entire observation period further indicate that the high‐frequency part (up to 100 Hz) represents the scattered wave field. The temporal variation of apparent velocities based on first‐arrival times partly correlates with the water injection rate and occurrence of microseismic events. It is concluded that borehole ambient noise interferometry in production settings is a potentially useful method for permanent reservoir monitoring due to its high sensitivity and robustness.