CN111025392B - Coal rock body fracturing crack real-time rapid monitoring and evaluation method utilizing microseismic signals - Google Patents
Coal rock body fracturing crack real-time rapid monitoring and evaluation method utilizing microseismic signals Download PDFInfo
- Publication number
- CN111025392B CN111025392B CN201911378413.2A CN201911378413A CN111025392B CN 111025392 B CN111025392 B CN 111025392B CN 201911378413 A CN201911378413 A CN 201911378413A CN 111025392 B CN111025392 B CN 111025392B
- Authority
- CN
- China
- Prior art keywords
- microseismic
- curve
- waveform
- time window
- fracturing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000011435 rock Substances 0.000 title claims abstract description 98
- 239000003245 coal Substances 0.000 title claims abstract description 93
- 238000012544 monitoring process Methods 0.000 title claims abstract description 43
- 238000011156 evaluation Methods 0.000 title claims description 15
- 238000000034 method Methods 0.000 claims abstract description 39
- 230000000977 initiatory effect Effects 0.000 claims abstract description 31
- 230000008569 process Effects 0.000 claims description 21
- 238000005553 drilling Methods 0.000 claims description 20
- 230000035515 penetration Effects 0.000 claims description 6
- 238000001514 detection method Methods 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 3
- 239000002243 precursor Substances 0.000 claims description 3
- 238000012935 Averaging Methods 0.000 claims description 2
- 230000008859 change Effects 0.000 abstract description 9
- 230000007246 mechanism Effects 0.000 abstract description 4
- 238000004364 calculation method Methods 0.000 abstract description 3
- 239000012530 fluid Substances 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 6
- 238000013461 design Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 108091034117 Oligonucleotide Proteins 0.000 description 1
- 241000695274 Processa Species 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000009412 basement excavation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000005507 spraying Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/34—Displaying seismic recordings or visualisation of seismic data or attributes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/646—Fractures
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
A method for rapidly monitoring and evaluating a coal rock body fracturing crack in real time by using a microseismic signal comprises the steps of firstly, calculating microseismic waveform characteristic parameters such as channel average waveform number, time window maximum amplitude sum, channel average pulse number, time window waveform average duration, time window waveform average dominant frequency and the like in any time window; on the basis of the above, respectively obtainingA curve,A curve,A curve,Curves andand according to the size of the longitudinal coordinate value of the curve and the curve change trend, monitoring and evaluating the initiation, expansion and communication of the main fracture of the coal rock body fracturing and the initiation of the branch fracture without carrying out complicated microseismic positioning and seismic source mechanism inversion calculation.
Description
Technical Field
The invention belongs to the technical field of fracture monitoring, relates to monitoring of coal rock body fractures, and particularly relates to a real-time rapid monitoring and evaluation method of coal rock body fractures by using microseismic signals.
Background
In recent years, the fracturing technology is widely applied to the aspects of oil gas and hot dry rock development, mine dynamic disaster prevention and control and the like. The fracturing technology is characterized in that high-pressure fluid is continuously injected into the coal rock body through a fracturing well or a drilling hole, so that the coal rock is fractured and a fracturing crack is generated; and high-pressure fluid is continuously injected, and the cracks can be expanded towards the inside of the coal rock body to form a crack network system mainly comprising main cracks and wing-shaped branch cracks. The rapid real-time monitoring and evaluation of fracture initiation, expansion and expansion of the fracturing fracture have important significance for improving the fracturing effect, optimizing the fracturing scheme design, preventing and treating fracturing-induced dynamic disasters and the like.
The most common coal rock fracturing monitoring method on site is to monitor and evaluate fracturing fractures by observing a pressure curve, however, the coal rock stratum on site is complex and changeable, the pressure curve can only indicate the pressure in a fracturing pipeline, and the pressure in the fractures and at the tips of the fractures in the fracturing process cannot be truly reflected; and the pressure difference with the fracturing crack and the pressure at the tip of the crack is larger, when the fracturing area has special structures such as a latent fault, a hole, a goaf and the like, a pressure curve can be sharply reduced, but the fracturing crack does not crack or expand at the moment, and misjudgment is easily caused in actual operation. Before and after the on-site coal rock body fracturing, checking holes are usually arranged on the periphery of a fracturing drilling well or a drilling hole, and the fracturing process is monitored by comparing the drilling cutting quantity of the checking holes before and after fracturing, the water content of the coal rock body, the stress and the like. Firstly, the method can only concentrate on a certain point in the potential fracturing influence range, cannot continuously monitor all point positions of the whole fracturing area, and cannot continuously monitor in real time; secondly, the method mainly carries out point evaluation on the fracturing affected zone through stress and water content change, and cannot carry out real-time continuous monitoring on fracture initiation and expansion of the fracturing fracture. The fracturing effect is evaluated by monitoring the change of the oil and gas yield after fracturing, and because oil and gas extraction can be carried out after drainage and exhaust for a certain time after fracturing, the method cannot carry out real-time continuous monitoring on the fracturing process; in addition, the method can only indirectly evaluate the fracturing effect through the oil gas yield change after fracturing, and cannot carry out real-time continuous monitoring evaluation on the fracture initiation and expansion process of the fracturing fracture.
In recent years, microseismic technology is adopted to monitor coal and rock fracturing at home and abroad. At present, the microseismic monitoring technology is mainly used for positioning and monitoring the coal and rock body fracturing process by a seismic source positioning method. However, due to the large computation amount of microseismic positioning, the theoretical and technical requirements on a monitoring and computing system and practitioners are high, the monitoring efficiency is low, the monitoring cost is high, and the real-time monitoring is difficult to realize; more importantly, the positioning result of the micro-seismic source fractured by the coal rock mass only can reflect the form of the fractured crack, and the fracture initiation and expansion processes of the main crack and the branch crack of the fractured crack are difficult to be rapidly monitored in real time.
In summary, a method for rapidly monitoring the fracture initiation and propagation of a coal rock body fracture is urgently needed so as to provide a basis for real-time adjustment of fracture parameters in the fracturing implementation process.
Disclosure of Invention
The invention aims to provide a method for rapidly monitoring and evaluating a coal rock body fracturing crack in real time by using a microseismic signal, which has simple steps and is rapid in real time.
In order to achieve the purpose, the technical scheme adopted by the invention is as follows: a real-time rapid monitoring and evaluation method for coal rock body fracturing cracks by using microseismic signals comprises the following specific steps:
(1) microseismic signal background noise of coal rock mass fracturing environment and induced original microseismic signal acquisition
N microseismic sensors are arranged at the periphery of a coal rock fracturing well or a fracturing drill hole and are sequentially marked as No. 1,2,3, …, i, … and N microseismic sensors (channels). Before the fracturing operation begins, the background noise of microseismic signals of the fracturing environment in a period of time is collected, and the average value of the absolute value of the amplitude of the background noise is calculatedAnd frequency distribution rules, providing input data for fracturing to induce effective microseismic signal detection; and starting coal rock mass fracturing, and continuously acquiring original micro-seismic signals induced in the whole coal rock mass fracturing process in real time.
(2) Effective microseismic waveform detection induced by coal rock body fracturing and characteristic parameter extraction thereof
Selecting a reasonable time window (time window for short) according to the actual situation of coal rock mass fracturing, recording the time length (duration for short) of the time window as L, and sequentially dividing original microseismic signals monitored in the whole coal rock mass fracturing process into M (j is 1,2,3, …, M) time windows with the time length of L according to the time sequence;
scanning and detecting the original microseismic signals monitored by each sensor in all time windows by adopting a long-short time window method based on instantaneous frequency and energy, and automatically identifying effective microseismic waveforms (effective microseismic waveform signals) induced in the whole process of coal rock fracturing; on the basis, the arrival time, the ending time, the waveform duration, the maximum amplitude and the main frequency of each effective microseismic waveform are respectively calculated; defining the number of effective microseismic waveforms monitored by each microseismic sensor in the jth (j is 1,2,3, …, M) time window as the waveform number, using letter W to represent the waveform number, and recording the waveform number of the ith microseismic sensor in the jth time window asOn the basis, the number of effective microseismic waveforms monitored by each microseismic sensor in any time window is recorded as the average waveform number of channels in the time window, and the average waveform number is usedIndicating that the average waveform number of the channel in the jth time window is usedIt is shown that,calculated by the following formula:
calculating to obtain each microseismic sensorMaximum amplitude A of the effective microseismic waveform monitored by the devicemax(ii) a Calculating to obtain the duration T of the effective microseismic waveform monitored by each microseismic sensord(ii) a Calculating to obtain the main frequency F of the effective microseismic waveform monitored by each microseismic sensorp(ii) a Averaging absolute values of background noise amplitudes of microseismic signalsAs a threshold value, respectively calculating the times of wave crests exceeding the threshold value in the effective microseismic waveform monitored by each microseismic sensor, and defining the times as the pulse number H;
recording the maximum amplitude sum of effective microseismic waveforms monitored by all microseismic sensors in any time window as the maximum amplitude sum of the time window, and usingIndicating that the maximum amplitude sum of the time windows in the jth time window isIt is shown that,calculated by the following formula:
in the above formula (2), i is 1,2,3, …, N,representing the maximum amplitude of the kth effective microseismic waveform of the ith microseismic sensor in the jth time window;
recording the average pulse number monitored by each microseismic sensor in any time window as the average pulse number of the channel in the time window, and using the average pulse numberIndicates that the channel in the jth time window is averagedFor the number of pulsesIt is shown that,calculated by the following formula:
in the above formula (3), i is 1,2,3, …, N,representing the number of pulses in the kth effective microseismic waveform of the ith microseismic sensor in the jth time window;
recording the average duration of all effective microseismic waveforms monitored by all microseismic sensors in any time window as the average duration of the time window waveform, and usingIndicating that the average duration of the waveform in the jth time window isIt is shown that,calculated by the following formula:
in the above formula (4), i is 1,2,3, …, N, representing the k-th effective microseismic waveform of the i-th microseismic sensor in the j-th time windowA duration of time;
recording the average main frequency of all effective microseismic waveforms monitored by all sensors in any time window as the average main frequency of the time window waveform, and using the average main frequencyIndicating that the average dominant frequency of the waveform in the jth time windowIt is shown that,calculated by the following formula:
in the above formula (5), i is 1,2,3, …, N, representing the dominant frequency of the kth effective microseismic waveform of the ith microseismic sensor in the jth time window.
(3) Real-time rapid monitoring and evaluation of coal rock body fracturing crack based on micro-seismic waveform characteristic parameters
Respectively drawing the average waveform number of the channelsMaximum amplitude sum of time windowNumber of channel average pulsesAverage duration of time window waveformTime window waveform average dominant frequencyGraphs of time window number j (which is essentially time) and are separately labeledA curve,A curve,A curve,Curves anda curve;
according toA curve,A curve,A curve,Curves andand monitoring and evaluating the cracking, expansion, penetration and branch crack initiation of the main crack of the coal rock body fracturing by the size of the longitudinal coordinate value of the curve and the change trend of each curve.
When in useCurves andthe slope of the curve simultaneously presents a turning point which is obviously increased, anCurves andthe curve shows local sudden drop, andwhen the curve is locally suddenly increased, the fracture initiation of the coal rock fracturing main fracture is shown;
thereafter, whenCurve or line ofThe slope of any one of the curves increases; orThe slope of the curve increases and at the same timeWhen the slope of the curve is reduced, the coal rock mass fracturing is shown to have branch crack initiation;
when the main fracture of the coal rock body fracturing is initiated,a curve,Curves andthe curve continues to increase; and the average duration of the time window waveformThe duration is greater than the average duration of the time window waveform of the fracture initiation stage of the main fractureNamely, it isTime window waveform average dominant frequencyTime window waveform average main frequency continuously smaller than fracture initiation stage of main fractureNamely, it isIndicating that the coal rock fracturing main crack enters a stable expansion stage;
when in useA curve,A curve,The curves are simultaneously subjected to global sudden drop and are dropped to the global minimum value and are continuously positioned at the minimum value; and isThe curve is abruptly increased to a global maximum,the curve suddenly drops to the global minimum value, which shows that the coal rock fracturing main fracture penetrates through the existing fault, large fracture or roadway section and the like, and the continuous fracturing is meaningless, so that the curve can be used as a microseismic signal precursor characteristic for finishing the fracturing operation.
The coal rock body fracturing fracture is generated under the action of high-pressure fluid. Compared with the microseismic waveform induced by the coal rock body fracture and crack caused by stress without participation of fluid, the microseismic waveform induced by the fracture and crack has unique characteristics in amplitude frequency, time frequency, duration, pulse number and other parameters due to the participation of high-pressure fluid. In addition, a main fracture and a branch fracture can be fractured in the coal rock fracturing process, and different microseismic waveform signals can also be generated when the main fracture and the branch fracture are initiated; the induced microseismic waveform characteristic parameters of the cracks are different at different stages of initiation, expansion, penetration and the like due to stress, fluid pressure, fluid flow state and the like. The invention systematically monitors the initiation, expansion, penetration and the like of the coal-rock body fracturing fracture in real time and rapidly according to the difference characteristics of microseismic waveform parameters induced by the coal-rock body fracturing for the first time, and mainly has the following beneficial effects:
(1) the invention utilizes the microseism waveform characteristic parameters induced by coal rock body fracturing and the change characteristics of the microseism waveform characteristic parameters along with time to realize real-time and rapid monitoring and evaluation of different stages of initiation, expansion, penetration and the like of coal rock body fracturing cracks, and is suitable for monitoring and evaluation technologies of various fracturing cracks of coal rock body hydraulic fracturing, carbon dioxide gas fracturing and the like in the fields of oil gas, mines, dry and hot rocks and the like; in addition, the method can also be applied to monitoring and evaluating coal rock stratum fractures caused by excavation of underground engineering such as mines, tunnels and the like.
(2) The traditional coal rock body fracturing micro-seismic monitoring method needs complex and large-computation-amount computing work such as micro-seismic source positioning, micro-seismic wave velocity imaging, a seismic source mechanism and the like, and the traditional micro-seismic monitoring time cost and the economic cost are very high; in addition, the analysis and interpretation of monitoring results such as microseismic positioning, wave velocity imaging, seismic source mechanisms and the like are obscure and difficult to understand, and the traditional microseismic monitoring has high requirements on professional theories and technologies of related practitioners. The method for rapidly monitoring the coal rock body fracturing fracture in real time by utilizing the micro-seismic waveform signal characteristics has the advantages of small calculated amount, high calculating speed and simple and understandable result analysis and explanation; the method not only greatly reduces the calculation workload, the economic cost and the time cost of microseismic monitoring, but also greatly reduces the difficulty of analyzing and explaining the monitoring result, thereby reducing the professional technical requirements on related operators.
(3) The coal rock body fracturing induced microseismic waveform characteristic parameters and the change characteristics of the characteristic parameters along with time, which are obtained by the invention, can provide reference and comparison for positioning of a fracturing fracture microseismic seismic source, inversion of a seismic source mechanism and the like, and have important significance for real-time adjustment of fracturing parameters in the design and implementation processes of coal rock body fracturing schemes.
Drawings
FIG. 1 is a flow chart of the method of the present invention;
FIG. 2 is a schematic diagram of characteristic parameters of microseismic waveforms induced by coal-rock body fracturing: (a) inducing effective microseismic wave forms by coal rock body fracturing; (b) a spectrogram corresponding to the effective microseismic waveform;
fig. 3 is an effective microseismic waveform induced by coal rock body fracturing monitored by the microseismic sensor No. i (i ═ 8) of the jth (j ═ 29) time window actually monitored in the process of coal mine underground coal rock body drilling fracturing: (a) 9 effective microseismic waveforms induced by coal rock mass drilling fracturing; (b)1# effective microseismic waveform enlarged view; (c) a time-frequency curve corresponding to the 1# effective microseismic waveform;
Detailed Description
The invention is described in further detail below with reference to the figures and specific examples.
As shown in fig. 1A flow chart of a real-time rapid monitoring and evaluation method for coal rock body fracturing cracks by using microseismic signals. FIG. 2 is a schematic diagram of characteristic parameters of microseismic waveform induced by coal-rock body fracturing, and it can be seen from the diagram that the average value of the absolute values of the amplitudes of the background noise isThe arrival time of effective microseismic wave shape induced by fracturing is taEnd time teThe waveform has a duration of Td=te-ta(ii) a The pulse number H of the waveform is 12; the maximum amplitude of the waveform is Amax(ii) a The dominant frequency of the waveform is Fp。
Fig. 3 shows effective microseismic waveforms induced by coal rock body fracturing monitored by a microseismic sensor No. i (8) in a jth (j (29)) time window of coal rock body drilling fracturing under a coal mine. As can be seen from fig. 3(a), the time window is 120 seconds long, and the number 8 microseismic sensors in the time window monitor 9 effective microseismic waveforms induced by fracturing. FIG. 3(b) is an enlarged view of the 1# effective microseismic waveform in the time window, FIG. 3(c) is a time-frequency curve corresponding to the 1# effective microseismic waveform, and it can be seen from FIGS. 3(b) and (c) that the average value of the absolute values of the background noise amplitudes is the average valueThe background noise frequency is mainly distributed around 15Hz, and the main frequency of the 1# effective microseismic waveform is 109 Hz.
The coal mine underground coal rock mass drilling fracturing micro-seismic monitoring is carried out by arranging 12 micro-seismic sensors in total, namely N is 12; the microseismic sensors are numbered as sensors No. 1-12 in sequence.
According to the coal rock body fracturing crack micro-seismic monitoring process, before the coal rock body drilling and fracturing of the coal mine, micro-seismic signal background noise test is carried out, and the average value of the absolute value of the amplitude of the background noise is calculatedThe frequency of the background noise is mainly concentrated on about 15Hz, and the amplitude and the frequency information of the background noise provide input data for detecting effective microseismic signals induced by the coal-rock body fracturingAccordingly.
According to the actual situation of the coal mine underground coal rock body drilling fracturing microseismic monitoring, a reasonable time window is determined, the time duration of the time window L is 2 minutes and 120 seconds, and original microseismic signals monitored by coal rock body fracturing are sequentially divided into M (j is 1,2,3, …, M) time windows with the time duration of L according to the time sequence.
And scanning and detecting microseismic signals monitored by each sensor in all time windows by adopting a long-time window method based on instantaneous frequency and energy, automatically identifying effective microseismic waveforms induced in the whole process of coal rock mass drilling fracturing, and respectively calculating the arrival time, the ending time, the waveform duration, the maximum amplitude and the main frequency of each effective microseismic waveform on the basis. Calculating to obtain the number of effective microseismic waveforms monitored by each microseismic sensor in each time window, namely calculating to obtain the number of effective microseismic waveforms monitored by the ith microseismic sensor in the jth time window, and recording the number asAnd setting the average value of the absolute values of the amplitudes of the background noise as a threshold value, and calculating to obtain the pulse number in each effective microseismic waveform duration.
Calculating the average waveform number of the channels in the jth time window according to the following formula
Calculating the sum of the maximum amplitudes of the effective microseismic waveforms monitored by all the sensors in the jth time window according to the following formula, namely the sum of the maximum amplitudes of the instant windows
The average pulse number monitored by each sensor in the jth time window is calculated according to the following formula, namely the average pulse number of the channels
Calculating the average duration of all effective microseismic waveforms monitored by all sensors in the jth time window according to the following formula, namely the average duration of the waveform of the time window
Calculating the average main frequency of all effective microseismic waveforms monitored by all sensors in the jth time window according to the following formula, and calculating the average main frequency of the waveform of the real-time window
Respectively drawing the channel average waveform number according to the calculation resultsMaximum amplitude sum of time windowNumber of channel average pulsesAverage duration of time window waveformTime window waveform average dominant frequencyThe time window number j is plotted and recorded asA curve,A curve,A curve,Curves andthe curves, the respective graphs are shown in detail in fig. 4.
According to the detailed characteristics of microseismic waveform signals and the change along with time given in figure 1, the fracture initiation, expansion, penetration and branch fracture initiation of the main fracture of the coal rock body fracture are monitored and evaluated, and the specific process and the monitoring and evaluation results are as follows:
(1) in the context of figure 4(a),the slope of curve OA is gw,1The slope of the curve in the AB section is gw,2The slope of the curve at the BC section is gw,3The slope of the curve in the CD segment is gw,4Slope of curve in DE section is gw,5The slope of the EF-segment curve is gw,6And is andthe slope of each section of the curve satisfies the following condition: gw,5<gw,6<gw,1<gw,2<gw,3<gw,4(ii) a In FIG. 4(b)The slope of curve OA is gH,1The slope of the curve in the AB section is gH,2The slope of the BD segment curve is gH,3Slope of curve in DE section is gH,4The slope of the EF-segment curve is gH,5And is andthe slope of each section of the curve satisfies the following condition: gH,4<gH,5<gH,1<gH,3<gH,2。
(2) When in useCurves andthe slope of the curve simultaneously presents a turning point which is obviously increased, anCurves andthe curve shows local sudden drop, andwhen the curve is locally suddenly increased, the fracture initiation of the coal rock fracturing main fracture is shown. As can be seen from fig. 4, in the a-th time window, i.e. at point a on each curve,curves andthe slope of the curve increases simultaneously, and this timeCurves andthe curve is locally suddenly dropped, and the local sudden drop occurs,the curve is locally suddenly increased, so that the main fracture of the coal rock drilling fracture is initiated.
(3) After the main crack of the coal rock mass drilling and fracturing is initiated, the main crack is treatedCurve or line ofThe slope of any one curve in the curves is increased; orThe slope of the curve increases and at the same timeWhen the slope of the curve is reduced, the coal rock mass fracturing is shown to have branch crack initiation. As can be seen from FIG. 4, at the b-th time window, i.e., atCurves andwhen the point B is on the curve, the point B,the slope of the curve increases and at the same timeThe slope of the curve is reduced, which indicates that the coal rock drilling fracture has branch fracture initiation at the moment. In addition, at the time of the c-th time window,the slope of the curve increases, and the branch crack initiation occurs again in the coal rock drilling and fracturing process.
(4) When coal is usedAfter the main fracture of the rock mass fracturing is initiated,a curve,Curves andthe curve continues to increase; and the average duration of the time window waveformThe duration is greater than the average duration of the time window waveform of the fracture initiation stage of the main fractureNamely, it isTime window waveform average dominant frequencyTime window waveform average main frequency continuously smaller than fracture initiation stage of main fractureNamely, it is And indicating that the coal rock fracturing main crack enters a stable expansion stage. As can be seen from FIG. 4, the main fracture of the coal rock drilling fracture starts at the a-th time window until the d-th time window, and in the processA curve,Curves andthe curve continues to increase; and in the process, the average duration of the time window waveformIs always longer than the average duration time of the time window waveform of the fracture initiation stage of the main fractureNamely, it isTime window waveform average dominant frequencyThe average main frequency of the waveform of the time window always smaller than the fracture initiation stage of the main fractureNamely, it isThis indicates that the coal-rock body borehole fracture main fracture is in a stable propagation stage from the a-th time window to the d-th time window.
(5) When in useA curve,A curve,The curves simultaneously have global sudden drop and drop to a global minimum value, and then continuously fluctuate near the minimum value; and isThe curve abruptly increases to a global maximum, anThe curve suddenly drops to the global minimum value, which shows that the main fracture of the coal rock mass fracturing is through and the fracturing is continued to lose significance, and the curve can be used as the precursor characteristic of the microseismic signal for finishing the fracturing operation. As can be seen from fig. 4, after the d-th time window,a curve,A curve,The curve also shows a sudden sharp decrease and drops to a global minimum in a very short time, i.e. at the e-th time window, and thereafter continues at the minimum until the fracturing is finished; and isThe curve rapidly increases after the d-th time window to the global maximum at the e-th time window,and the curve rapidly drops to the global minimum value of the time window e after the time window d, which shows that the main fracture of the coal-rock body drilling fracture starts to accelerate and expand after the time window d until the main fracture of the coal-rock body drilling fracture is communicated with the time window e. The coal mine underground field observation also shows that the water spraying phenomenon occurs on the roadway roof around the fracturing drill hole in the e-th time window, which shows that the main fracturing crack of the coal rock body drill hole penetrates through the broken zone of the roadway roof, and simultaneously shows that the monitoring and evaluation method provided by the invention is feasible and has important significance for the real-time adjustment of the fracturing parameters in the design and implementation processes of the coal rock body fracturing scheme.
The above description is only of the preferred embodiments of the present invention, and it should be noted that: it will be apparent to those skilled in the art that various modifications and adaptations can be made without departing from the principles of the invention and these are intended to be within the scope of the invention.
Claims (1)
1. A real-time rapid monitoring and evaluation method for coal rock body fracturing cracks by using microseismic signals is characterized by comprising the following steps: the method comprises the following specific steps:
(1) microseismic signal background noise of coal rock mass fracturing environment and induced original microseismic signal acquisition
Arranging N microseismic sensors around a coal rock fracturing well or a fracturing drilling hole, sequentially recording as No. 1,2,3, L, i, L and N microseismic sensors, collecting background noise of microseismic signals of a fracturing environment for a period of time before the fracturing operation starts, and calculating to obtain the average value of the absolute value of the amplitude of the background noiseAnd the frequency distribution rule thereof, and provides input data for the detection of effective microseismic signals induced by fracturing; starting coal rock mass fracturing, and continuously acquiring original micro-seismic signals induced in the whole coal rock mass fracturing process in real time;
(2) effective microseismic waveform detection induced by coal rock body fracturing and characteristic parameter extraction thereof
Selecting a reasonable time window according to the actual situation of coal rock body fracturing, recording the time length of the time window as L, and sequentially dividing the original microseismic signals monitored in the whole coal rock body fracturing process into M (j is 1,2,3, L, M) time windows with the time length of L according to the time sequence;
scanning and detecting original microseismic signals monitored by each sensor in all time windows by adopting a long-time window method based on instantaneous frequency and energy to automatically identify effective microseismic waveforms induced in the whole process of coal rock fracturing, respectively calculating the arrival time, the ending time, the waveform duration, the maximum amplitude and the dominant frequency of each effective microseismic waveform on the basis, defining the number of the effective microseismic waveforms monitored by each microseismic sensor in the jth (j is 1,2,3, L, M) time window as a waveform number, representing the waveform number by using a letter W, and recording the waveform number of the ith microseismic sensor in the jth time window as the waveform numberOn the basis, the number of effective microseismic waveforms monitored by each microseismic sensor in the jth time window is recorded as the average waveform number of channels in the jth time window for useIs shown to beCalculated by the following formula:
calculating to obtain the maximum amplitude A of the effective microseismic waveform monitored by each microseismic sensormax(ii) a Calculating to obtain the duration T of the effective microseismic waveform monitored by each microseismic sensord(ii) a Calculating to obtain the main frequency F of the effective microseismic waveform monitored by each microseismic sensorp(ii) a Averaging absolute values of background noise amplitudes of microseismic signalsAs a threshold value, respectively calculating the times of wave crests exceeding the threshold value in the effective microseismic waveform monitored by each microseismic sensor, and defining the times as the pulse number H;
recording the maximum amplitude sum of effective microseismic waveforms monitored by all microseismic sensors in the jth time window as the maximum amplitude sum of the time window, and using the sumIs shown to beCalculated by the following formula:
in the above formula (2), i is 1,2,3, L, N, k is 1,2,3, L, representing the maximum amplitude of the kth effective microseismic waveform of the ith microseismic sensor in the jth time window;
recording the average pulse number monitored by each microseismic sensor in the jth time window as the average pulse number of the channel, and using the average pulse numberIs shown to beCalculated by the following formula:
in the above formula (3), i is 1,2,3, L, N, k is 1,2,3, L, representing the number of pulses in the kth effective microseismic waveform of the ith microseismic sensor in the jth time window;
recording the average duration of all effective microseismic waveforms monitored by all microseismic sensors in the jth time window as the average duration of the time window waveform, and using the average durationIs shown to beCalculated by the following formula:
in the above formula (4), i is 1,2,3, L, N, k is 1,2,3, L, representing the duration of the kth effective microseismic waveform of the ith microseismic sensor in the jth time window;
recording the average main frequency of all effective microseismic waveforms monitored by all microseismic sensors in the jth time window as the average main frequency of the time window waveform, and using the average main frequencyIs shown to beCalculated by the following formula:
in the above formula (5), i is 1,2,3, L, N, k is 1,2,3, L, representing the main frequency of the kth effective microseismic waveform of the ith microseismic sensor in the jth time window;
(3) real-time rapid monitoring and evaluation of coal rock body fracturing crack based on micro-seismic waveform characteristic parameters
Respectively drawing the average waveform number of the channelsMaximum amplitude sum of time windowNumber of channel average pulsesAverage duration of time window waveformTime window waveform average dominant frequencyThe time window number j is plotted and recorded asA-j curve,A-j curve,A-j curve,The j curve and-j-curve;
according toA-j curve,A-j curve,A-j curve,The j curve andand (3) monitoring and evaluating the fracture initiation, expansion, penetration and branch fracture initiation of the coal rock body fracturing main fracture by the longitudinal coordinate value of the j curve and the variation trend of each curve:
when in useThe j curve andthe slope of the j-curve simultaneously presents a turning point which is significantly increased, andthe j curve andthe curve-j shows a local dip, andwhen the curve-j is locally suddenly increased, the fracture initiation of the coal rock fracturing main fracture is indicated;
thereafter, whenThe j curve or-an increase in the slope of any one of the j curves; orThe slope of the-j curve increases and at the same timeWhen the slope of the-j curve is reduced, indicating that the coal rock mass fracturing has branch crack initiation;
when the main fracture of the coal rock body fracturing is initiated,a-j curve,The j curve andthe j curve increases continuously; and the average duration of the time window waveformThe duration is greater than the average duration of the time window waveform of the fracture initiation stage of the main fractureNamely, it isTime window waveform average dominant frequencyTime window waveform average main frequency continuously smaller than fracture initiation stage of main fractureNamely, it isIndicating that the coal rock fracturing main crack enters a stable expansion stage;
when in useA-j curve,A-j curve,The j curve is simultaneously subjected to global sudden drop and drops to the global minimum value and is continuously positioned at the minimum value; and isThe curve j is abruptly increased to a global maximum,and (4) suddenly reducing the curve-j to a global minimum value, indicating that the coal rock mass fracturing main fracture penetrates through the existing fault, large fracture or roadway section, continuing fracturing and losing significance, and taking the fracture as a microseismic signal precursor characteristic for finishing fracturing operation.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911378413.2A CN111025392B (en) | 2019-12-27 | 2019-12-27 | Coal rock body fracturing crack real-time rapid monitoring and evaluation method utilizing microseismic signals |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911378413.2A CN111025392B (en) | 2019-12-27 | 2019-12-27 | Coal rock body fracturing crack real-time rapid monitoring and evaluation method utilizing microseismic signals |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111025392A CN111025392A (en) | 2020-04-17 |
CN111025392B true CN111025392B (en) | 2020-07-24 |
Family
ID=70196380
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911378413.2A Active CN111025392B (en) | 2019-12-27 | 2019-12-27 | Coal rock body fracturing crack real-time rapid monitoring and evaluation method utilizing microseismic signals |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111025392B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111691871A (en) * | 2020-06-12 | 2020-09-22 | 鞍钢集团矿业有限公司 | Crack propagation monitoring method for hydrofracturing ore rock pretreatment test |
CN111636859B (en) * | 2020-07-09 | 2022-08-16 | 中煤科工集团重庆研究院有限公司 | Coal rock while-drilling self-identification method based on micro-fracture wave detection |
CN113050159B (en) * | 2021-03-23 | 2021-11-16 | 中国矿业大学 | Coal rock hydraulic fracturing crack micro-seismic positioning and propagation mechanism monitoring method |
CN113703038B (en) * | 2021-08-31 | 2024-05-07 | 中煤科工集团重庆研究院有限公司 | Automatic microseismic signal acquisition and identification method and system |
CN117289344B (en) * | 2023-11-24 | 2024-01-30 | 北京科技大学 | Quick coal rock destabilization damage judgment method based on seismic source spatial distribution |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016134210A1 (en) * | 2015-02-20 | 2016-08-25 | Schlumberger Technology Corporation | Microseismic sensitivity analysis and scenario modelling |
CN106918841A (en) * | 2017-05-05 | 2017-07-04 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | A kind of microseism fracturing fracture analysis method based on multi-factor comprehensive |
CN107728200A (en) * | 2017-09-29 | 2018-02-23 | 中国石油化工股份有限公司 | Ground micro-seismic fracturing fracture dynamic spread method of real-time |
CN108873081A (en) * | 2018-05-31 | 2018-11-23 | 湖南继善高科技有限公司 | A kind of main seam net of oil gas pressure break stitches four-dimensional electromagnetism method of real-time and system |
CN109595036A (en) * | 2018-12-14 | 2019-04-09 | 北京矿冶科技集团有限公司 | A kind of method for early warning of Mine Ground Pressure Disaster |
-
2019
- 2019-12-27 CN CN201911378413.2A patent/CN111025392B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016134210A1 (en) * | 2015-02-20 | 2016-08-25 | Schlumberger Technology Corporation | Microseismic sensitivity analysis and scenario modelling |
CN106918841A (en) * | 2017-05-05 | 2017-07-04 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | A kind of microseism fracturing fracture analysis method based on multi-factor comprehensive |
CN107728200A (en) * | 2017-09-29 | 2018-02-23 | 中国石油化工股份有限公司 | Ground micro-seismic fracturing fracture dynamic spread method of real-time |
CN108873081A (en) * | 2018-05-31 | 2018-11-23 | 湖南继善高科技有限公司 | A kind of main seam net of oil gas pressure break stitches four-dimensional electromagnetism method of real-time and system |
CN109595036A (en) * | 2018-12-14 | 2019-04-09 | 北京矿冶科技集团有限公司 | A kind of method for early warning of Mine Ground Pressure Disaster |
Non-Patent Citations (1)
Title |
---|
微震监测技术及其在煤矿的应用现状与展望;李楠;《煤炭学报》;20170731;第42卷;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111025392A (en) | 2020-04-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111025392B (en) | Coal rock body fracturing crack real-time rapid monitoring and evaluation method utilizing microseismic signals | |
CN113050159B (en) | Coal rock hydraulic fracturing crack micro-seismic positioning and propagation mechanism monitoring method | |
Zhu et al. | Interpretation of the extent of hydraulic fracturing for rockburst prevention using microseismic monitoring data | |
Dai et al. | Microseismic early warning of surrounding rock mass deformation in the underground powerhouse of the Houziyan hydropower station, China | |
CN104390537B (en) | A kind of side slope pre split Blasting Excavation damage control method based on blasting vibration test | |
CN204462405U (en) | A kind of rock burst omen early warning system based on acoustic emission | |
WO2019161593A1 (en) | Monitoring and early warning method for electromagnetic radiation and underground sound of coal and rock dynamic disaster hazard | |
CN101581234B (en) | Comprehensive underground test method for deformation and damage of terranes of mining top plate and mining bottom plate of coal bed | |
Li et al. | Characteristics of microseismic waveforms induced by hydraulic fracturing in coal seam for coal rock dynamic disasters prevention | |
NO334626B1 (en) | Procedure for monitoring fracture | |
CN106353792A (en) | Method suitable for positioning hydraulic fracturing micro-seismic source | |
Zhang et al. | A new monitoring-while-drilling method of large diameter drilling in underground coal mine and their application | |
CN105334548A (en) | Geological forecasting method for tunnel construction in karst area | |
RU2467171C1 (en) | Method of diagnosing dangerous situations in deep mining and forecasting parameters of fissuring zones formed by fracturing | |
CN107479098B (en) | Same-well micro-seismic monitoring process in hydraulic fracturing process | |
CN112731525B (en) | Intelligent prediction method for stability of surrounding rock of roadway based on synchronous monitoring of microseismic and electromagnetic radiation | |
Liu et al. | Application of surface–downhole combined microseismic monitoring technology in the Fuling shale gas field and its enlightenment | |
CN110259442B (en) | Coal measure stratum hydraulic fracturing fracture horizon identification method | |
Zhang et al. | Vibration events in underground heading face and useful index for rock burst monitoring | |
Ren et al. | Microseismic signals in heading face of tengdong coal mine and their application for rock burst monitoring | |
CN106032750B (en) | Geological logging instrument based on drilling energy spectrum | |
Zhang et al. | A multi-channel verification index to improve distinguish accuracy of target signals in rock burst monitoring of heading face | |
Stoeckhert et al. | Hydraulic fractures in discontinuous, anisotropic and heterogeneous rock-a lab study | |
CN106869904B (en) | A method of Rock Damage state is determined in real time using drilling machine operating parameter is in situ | |
Cheng et al. | Vibration behavior during underground drilling based on an innovative measurement method and the application |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |