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 PDF

Info

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
Application number
CN201911378413.2A
Other languages
Chinese (zh)
Other versions
CN111025392A (en
Inventor
李楠
孙炜辰
陈栋
张志博
房柳林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Mining and Technology CUMT
Original Assignee
China University of Mining and Technology CUMT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Mining and Technology CUMT filed Critical China University of Mining and Technology CUMT
Priority to CN201911378413.2A priority Critical patent/CN111025392B/en
Publication of CN111025392A publication Critical patent/CN111025392A/en
Application granted granted Critical
Publication of CN111025392B publication Critical patent/CN111025392B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/34Displaying seismic recordings or visualisation of seismic data or attributes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/64Geostructures, e.g. in 3D data cubes
    • G01V2210/646Fractures

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 obtaining
Figure DDA0002341630870000011
A curve,
Figure DDA0002341630870000012
A curve,
Figure DDA0002341630870000013
A curve,
Figure DDA0002341630870000014
Curves and
Figure DDA0002341630870000015
and 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

Coal rock body fracturing crack real-time rapid monitoring and evaluation method utilizing microseismic signals
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 calculated
Figure BDA0002341630850000021
And 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 as
Figure BDA0002341630850000031
On 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 used
Figure BDA0002341630850000032
Indicating that the average waveform number of the channel in the jth time window is used
Figure BDA0002341630850000033
It is shown that,
Figure BDA0002341630850000034
calculated by the following formula:
Figure BDA0002341630850000035
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 signals
Figure BDA0002341630850000036
As 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 using
Figure BDA0002341630850000037
Indicating that the maximum amplitude sum of the time windows in the jth time window is
Figure BDA0002341630850000038
It is shown that,
Figure BDA0002341630850000039
calculated by the following formula:
Figure BDA00023416308500000310
in the above formula (2), i is 1,2,3, …, N,
Figure BDA00023416308500000311
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 number
Figure BDA00023416308500000312
Indicates that the channel in the jth time window is averagedFor the number of pulses
Figure BDA00023416308500000313
It is shown that,
Figure BDA00023416308500000314
calculated by the following formula:
Figure BDA00023416308500000315
in the above formula (3), i is 1,2,3, …, N,
Figure BDA00023416308500000316
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 using
Figure BDA0002341630850000041
Indicating that the average duration of the waveform in the jth time window is
Figure BDA0002341630850000042
It is shown that,
Figure BDA0002341630850000043
calculated by the following formula:
Figure BDA0002341630850000044
in the above formula (4), i is 1,2,3, …, N,
Figure BDA0002341630850000045
Figure BDA0002341630850000046
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 frequency
Figure BDA0002341630850000047
Indicating that the average dominant frequency of the waveform in the jth time window
Figure BDA0002341630850000048
It is shown that,
Figure BDA0002341630850000049
calculated by the following formula:
Figure BDA00023416308500000410
in the above formula (5), i is 1,2,3, …, N,
Figure BDA00023416308500000411
Figure BDA00023416308500000412
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 channels
Figure BDA00023416308500000413
Maximum amplitude sum of time window
Figure BDA00023416308500000414
Number of channel average pulses
Figure BDA00023416308500000415
Average duration of time window waveform
Figure BDA00023416308500000416
Time window waveform average dominant frequency
Figure BDA00023416308500000417
Graphs of time window number j (which is essentially time) and are separately labeled
Figure BDA00023416308500000418
A curve,
Figure BDA00023416308500000419
A curve,
Figure BDA00023416308500000420
A curve,
Figure BDA00023416308500000421
Curves and
Figure BDA00023416308500000422
a curve;
according to
Figure BDA00023416308500000423
A curve,
Figure BDA00023416308500000424
A curve,
Figure BDA00023416308500000425
A curve,
Figure BDA00023416308500000426
Curves and
Figure BDA00023416308500000427
and 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 use
Figure BDA0002341630850000051
Curves and
Figure BDA0002341630850000052
the slope of the curve simultaneously presents a turning point which is obviously increased, an
Figure BDA0002341630850000053
Curves and
Figure BDA0002341630850000054
the curve shows local sudden drop, and
Figure BDA0002341630850000055
when the curve is locally suddenly increased, the fracture initiation of the coal rock fracturing main fracture is shown;
thereafter, when
Figure BDA0002341630850000056
Curve or line of
Figure BDA0002341630850000057
The slope of any one of the curves increases; or
Figure BDA0002341630850000058
The slope of the curve increases and at the same time
Figure BDA0002341630850000059
When 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,
Figure BDA00023416308500000510
a curve,
Figure BDA00023416308500000511
Curves and
Figure BDA00023416308500000512
the curve continues to increase; and the average duration of the time window waveform
Figure BDA00023416308500000513
The duration is greater than the average duration of the time window waveform of the fracture initiation stage of the main fracture
Figure BDA00023416308500000514
Namely, it is
Figure BDA00023416308500000515
Time window waveform average dominant frequency
Figure BDA00023416308500000516
Time window waveform average main frequency continuously smaller than fracture initiation stage of main fracture
Figure BDA00023416308500000517
Namely, it is
Figure BDA00023416308500000518
Indicating that the coal rock fracturing main crack enters a stable expansion stage;
when in use
Figure BDA00023416308500000519
A curve,
Figure BDA00023416308500000520
A curve,
Figure BDA00023416308500000521
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 is
Figure BDA00023416308500000522
The curve is abruptly increased to a global maximum,
Figure BDA00023416308500000523
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;
fig. 4 is a change curve of each microseismic waveform characteristic parameter induced by coal-rock mass fracturing drilling along with a time window sequence number: (a)
Figure BDA0002341630850000061
curve, (b)
Figure BDA0002341630850000062
Curve, (c)
Figure BDA0002341630850000063
Curve (d)
Figure BDA0002341630850000064
Curve, (e)
Figure BDA0002341630850000065
Curve line.
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 is
Figure BDA0002341630850000071
The 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 value
Figure BDA0002341630850000072
The 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 calculated
Figure BDA0002341630850000073
The 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 as
Figure BDA0002341630850000074
And 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
Figure BDA0002341630850000075
Figure BDA0002341630850000081
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
Figure BDA0002341630850000082
Figure BDA0002341630850000083
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
Figure BDA0002341630850000084
Figure BDA0002341630850000085
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
Figure BDA0002341630850000086
Figure BDA0002341630850000087
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
Figure BDA0002341630850000088
Figure BDA0002341630850000089
Respectively drawing the channel average waveform number according to the calculation results
Figure BDA00023416308500000810
Maximum amplitude sum of time window
Figure BDA00023416308500000811
Number of channel average pulses
Figure BDA00023416308500000812
Average duration of time window waveform
Figure BDA00023416308500000813
Time window waveform average dominant frequency
Figure BDA00023416308500000814
The time window number j is plotted and recorded as
Figure BDA00023416308500000815
A curve,
Figure BDA00023416308500000816
A curve,
Figure BDA00023416308500000817
A curve,
Figure BDA00023416308500000818
Curves and
Figure BDA0002341630850000091
the 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),
Figure BDA0002341630850000092
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 and
Figure BDA0002341630850000093
the 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)
Figure BDA0002341630850000094
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 and
Figure BDA0002341630850000095
the slope of each section of the curve satisfies the following condition: gH,4<gH,5<gH,1<gH,3<gH,2
(2) When in use
Figure BDA0002341630850000096
Curves and
Figure BDA0002341630850000097
the slope of the curve simultaneously presents a turning point which is obviously increased, an
Figure BDA0002341630850000098
Curves and
Figure BDA0002341630850000099
the curve shows local sudden drop, and
Figure BDA00023416308500000910
when 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,
Figure BDA00023416308500000911
curves and
Figure BDA00023416308500000912
the slope of the curve increases simultaneously, and this time
Figure BDA00023416308500000913
Curves and
Figure BDA00023416308500000914
the curve is locally suddenly dropped, and the local sudden drop occurs,
Figure BDA00023416308500000915
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 treated
Figure BDA00023416308500000916
Curve or line of
Figure BDA00023416308500000917
The slope of any one curve in the curves is increased; or
Figure BDA00023416308500000918
The slope of the curve increases and at the same time
Figure BDA00023416308500000919
When 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 and
Figure BDA00023416308500000921
when the point B is on the curve, the point B,
Figure BDA00023416308500000922
the slope of the curve increases and at the same time
Figure BDA00023416308500000923
The 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,
Figure BDA00023416308500000924
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,
Figure BDA00023416308500000925
a curve,
Figure BDA00023416308500000926
Curves and
Figure BDA00023416308500000927
the curve continues to increase; and the average duration of the time window waveform
Figure BDA0002341630850000101
The duration is greater than the average duration of the time window waveform of the fracture initiation stage of the main fracture
Figure BDA0002341630850000102
Namely, it is
Figure BDA0002341630850000103
Time window waveform average dominant frequency
Figure BDA0002341630850000104
Time window waveform average main frequency continuously smaller than fracture initiation stage of main fracture
Figure BDA0002341630850000105
Namely, it is
Figure BDA0002341630850000106
Figure BDA0002341630850000107
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 process
Figure BDA0002341630850000108
A curve,
Figure BDA0002341630850000109
Curves and
Figure BDA00023416308500001010
the curve continues to increase; and in the process, the average duration of the time window waveform
Figure BDA00023416308500001011
Is always longer than the average duration time of the time window waveform of the fracture initiation stage of the main fracture
Figure BDA00023416308500001012
Namely, it is
Figure BDA00023416308500001013
Time window waveform average dominant frequency
Figure BDA00023416308500001014
The average main frequency of the waveform of the time window always smaller than the fracture initiation stage of the main fracture
Figure BDA00023416308500001015
Namely, it is
Figure BDA00023416308500001016
This 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 use
Figure BDA00023416308500001017
A curve,
Figure BDA00023416308500001018
A curve,
Figure BDA00023416308500001019
The curves simultaneously have global sudden drop and drop to a global minimum value, and then continuously fluctuate near the minimum value; and is
Figure BDA00023416308500001020
The curve abruptly increases to a global maximum, an
Figure BDA00023416308500001021
The 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,
Figure BDA00023416308500001022
a curve,
Figure BDA00023416308500001023
A curve,
Figure BDA00023416308500001024
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 is
Figure BDA00023416308500001025
The curve rapidly increases after the d-th time window to the global maximum at the e-th time window,
Figure BDA00023416308500001026
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 noise
Figure FDA0002521928810000011
And 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 number
Figure FDA0002521928810000012
On 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 use
Figure FDA0002521928810000013
Is shown to be
Figure FDA0002521928810000014
Calculated by the following formula:
Figure FDA0002521928810000015
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 signals
Figure FDA0002521928810000021
As 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 sum
Figure FDA0002521928810000022
Is shown to be
Figure FDA0002521928810000023
Calculated by the following formula:
Figure FDA0002521928810000024
in the above formula (2), i is 1,2,3, L, N, k is 1,2,3, L,
Figure FDA0002521928810000025
Figure FDA0002521928810000026
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 number
Figure FDA0002521928810000027
Is shown to be
Figure FDA0002521928810000028
Calculated by the following formula:
Figure FDA0002521928810000029
in the above formula (3), i is 1,2,3, L, N, k is 1,2,3, L,
Figure FDA00025219288100000210
Figure FDA00025219288100000211
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 duration
Figure FDA00025219288100000212
Is shown to be
Figure FDA00025219288100000213
Calculated by the following formula:
Figure FDA00025219288100000214
in the above formula (4), i is 1,2,3, L, N, k is 1,2,3, L,
Figure FDA00025219288100000215
Figure FDA00025219288100000216
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 frequency
Figure FDA00025219288100000217
Is shown to be
Figure FDA00025219288100000218
Calculated by the following formula:
Figure FDA0002521928810000031
in the above formula (5), i is 1,2,3, L, N, k is 1,2,3, L,
Figure FDA0002521928810000032
Figure FDA0002521928810000033
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 channels
Figure FDA0002521928810000034
Maximum amplitude sum of time window
Figure FDA0002521928810000035
Number of channel average pulses
Figure FDA0002521928810000036
Average duration of time window waveform
Figure FDA0002521928810000037
Time window waveform average dominant frequency
Figure FDA0002521928810000038
The time window number j is plotted and recorded as
Figure FDA0002521928810000039
A-j curve,
Figure FDA00025219288100000310
A-j curve,
Figure FDA00025219288100000311
A-j curve,
Figure FDA00025219288100000312
The j curve and
Figure FDA00025219288100000313
-j-curve;
according to
Figure FDA00025219288100000314
A-j curve,
Figure FDA00025219288100000315
A-j curve,
Figure FDA00025219288100000316
A-j curve,
Figure FDA00025219288100000317
The j curve and
Figure FDA00025219288100000318
and (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 use
Figure FDA00025219288100000319
The j curve and
Figure FDA00025219288100000320
the slope of the j-curve simultaneously presents a turning point which is significantly increased, and
Figure FDA00025219288100000321
the j curve and
Figure FDA00025219288100000322
the curve-j shows a local dip, and
Figure FDA00025219288100000323
when the curve-j is locally suddenly increased, the fracture initiation of the coal rock fracturing main fracture is indicated;
thereafter, when
Figure FDA00025219288100000324
The j curve or
Figure FDA00025219288100000325
-an increase in the slope of any one of the j curves; or
Figure FDA00025219288100000326
The slope of the-j curve increases and at the same time
Figure FDA00025219288100000327
When 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,
Figure FDA00025219288100000328
a-j curve,
Figure FDA00025219288100000329
The j curve and
Figure FDA00025219288100000330
the j curve increases continuously; and the average duration of the time window waveform
Figure FDA00025219288100000331
The duration is greater than the average duration of the time window waveform of the fracture initiation stage of the main fracture
Figure FDA00025219288100000332
Namely, it is
Figure FDA00025219288100000333
Time window waveform average dominant frequency
Figure FDA00025219288100000334
Time window waveform average main frequency continuously smaller than fracture initiation stage of main fracture
Figure FDA00025219288100000335
Namely, it is
Figure FDA00025219288100000336
Indicating that the coal rock fracturing main crack enters a stable expansion stage;
when in use
Figure FDA00025219288100000337
A-j curve,
Figure FDA00025219288100000338
A-j curve,
Figure FDA00025219288100000339
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 is
Figure FDA0002521928810000041
The curve j is abruptly increased to a global maximum,
Figure FDA0002521928810000042
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.
CN201911378413.2A 2019-12-27 2019-12-27 Coal rock body fracturing crack real-time rapid monitoring and evaluation method utilizing microseismic signals Active CN111025392B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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