CN114859411A - Method and system for determining microseism magnitude - Google Patents
Method and system for determining microseism magnitude Download PDFInfo
- Publication number
- CN114859411A CN114859411A CN202210459414.5A CN202210459414A CN114859411A CN 114859411 A CN114859411 A CN 114859411A CN 202210459414 A CN202210459414 A CN 202210459414A CN 114859411 A CN114859411 A CN 114859411A
- Authority
- CN
- China
- Prior art keywords
- target monitoring
- monitoring point
- seismic
- micro
- determining
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000012544 monitoring process Methods 0.000 claims abstract description 172
- 230000005855 radiation Effects 0.000 claims abstract description 43
- 238000004364 calculation method Methods 0.000 claims description 12
- 238000004590 computer program Methods 0.000 claims description 4
- 238000005336 cracking Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 3
- 238000003860 storage Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 239000003245 coal Substances 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 239000003129 oil well Substances 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 238000006467 substitution reaction 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. analysis, for interpretation, for correction
- G01V1/288—Event detection in seismic signals, e.g. microseismics
-
- 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. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/12—Signal generation
- G01V2210/123—Passive source, e.g. microseismics
- G01V2210/1232—Earthquakes
-
- 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/63—Seismic attributes, e.g. amplitude, polarity, instant phase
Abstract
The application provides a method and a system for determining a microseismic magnitude, wherein the method comprises the following steps: carrying out vector scanning on the micro-seismic target monitoring points to obtain the fracture radiation energy corresponding to the target monitoring points; calculating the lower limit of the signal-to-noise ratio corresponding to the target monitoring point based on the rupture radiation energy corresponding to the target monitoring point; determining the equivalent micro-seismic energy corresponding to the target monitoring point according to the lower limit of the signal-to-noise ratio corresponding to the target monitoring point; and determining the micro-seismic magnitude corresponding to the target monitoring point according to the micro-seismic equivalent energy corresponding to the target monitoring point. According to the technical scheme, the micro-earthquake magnitude can be effectively and accurately determined through the micro-earthquake equivalent energy corresponding to the determined target monitoring point.
Description
Technical Field
The application relates to the technical field of micro-seismic monitoring, in particular to a method and a system for determining a micro-seismic magnitude.
Background
The determination of magnitude of a seismic event is typically considered significant in seismic records, and well-established methods exist, such as the richter magnitude, the moment magnitude, and the association between differently defined magnitudes, the relationship between magnitude and energy, and so forth. The determination of these magnitude levels is based on the energy released by the earthquake, or the energy is proportional to the square of the maximum amplitude of the recorded seismic event, or the fault movement moment; this requires specific measurements of the amplitude, period, diastole, etc. of the recorded seismic events and the calculation of the magnitude from the characteristics of the recording instrument, where the richter magnitude has been applied in the seismic world for nearly a hundred years.
However, because the microseisms are small, the salient amplitude of the visible seismic event can only be recorded close to the target of monitoring, and such proximity monitoring is often difficult; in order to expand the monitoring range of the microseismic monitoring station network, a large number of microseismic monitoring stations are more than several hundred meters away from a seismic source and reach two to three kilometers, so that the microseismic event records are completely submerged in background noise or even if barely visible, the microseismic magnitude cannot be effectively and accurately determined due to low signal-to-noise ratio. Therefore, it is highly desirable to find a method for determining microseismic magnitude that can communicate with existing systems for estimating the magnitude of the richter.
Disclosure of Invention
The application provides a method and a system for determining a microseismic magnitude, which are used for at least solving the technical problem that the microseismic magnitude cannot be effectively and accurately determined in the related technology.
The embodiment of the first aspect of the application provides a method for determining a microseismic magnitude, which comprises the following steps:
carrying out vector scanning on the micro-seismic target monitoring points to obtain the fracture radiation energy corresponding to the target monitoring points;
calculating the lower limit of the signal-to-noise ratio corresponding to the target monitoring point based on the rupture radiation energy corresponding to the target monitoring point;
determining the equivalent micro-seismic energy corresponding to the target monitoring point according to the lower limit of the signal-to-noise ratio corresponding to the target monitoring point;
and determining the micro-seismic magnitude corresponding to the target monitoring point according to the micro-seismic equivalent energy corresponding to the target monitoring point.
In a second aspect, an embodiment of the present application provides a system for determining a microseismic magnitude, the system including:
the scanning module is used for carrying out vector scanning on the micro-seismic target monitoring points to obtain the fracture radiation energy corresponding to the target monitoring points;
the calculation module is used for calculating the lower limit of the signal-to-noise ratio corresponding to the target monitoring point based on the rupture radiation energy corresponding to the target monitoring point;
the first determining module is used for determining the micro-seismic equivalent energy corresponding to the target monitoring point according to the lower limit of the signal-to-noise ratio corresponding to the target monitoring point;
and the second determination module is used for determining the micro-seismic magnitude corresponding to the target monitoring point according to the micro-seismic equivalent energy corresponding to the target monitoring point.
A computer-readable storage medium is provided in an embodiment of the third aspect of the present application, on which a computer program is stored, which when executed by a processor, implements the method as described in the first aspect above.
The technical scheme provided by the embodiment of the application at least has the following beneficial effects:
the application provides a method and a system for determining a microseismic magnitude, wherein the method comprises the following steps: carrying out vector scanning on the micro-seismic target monitoring points to obtain the fracture radiation energy corresponding to the target monitoring points; calculating the lower limit of the signal-to-noise ratio corresponding to the target monitoring point based on the rupture radiation energy corresponding to the target monitoring point; determining the equivalent micro-seismic energy corresponding to the target monitoring point according to the lower limit of the signal-to-noise ratio corresponding to the target monitoring point; and determining the micro-seismic magnitude corresponding to the target monitoring point according to the micro-seismic equivalent energy corresponding to the target monitoring point. According to the technical scheme, the micro-earthquake magnitude can be effectively and accurately determined through the micro-earthquake equivalent energy corresponding to the determined target monitoring point.
Additional aspects and advantages of the present application will be set forth in part in the description which follows and, in part, will be obvious from the description, or may be learned by practice of the present application.
Drawings
The foregoing and/or additional aspects and advantages of the present application will become apparent and readily appreciated from the following description of the embodiments, taken in conjunction with the accompanying drawings of which:
FIG. 1 is a flow chart of a method for determining a microseismic magnitude provided in accordance with one embodiment of the present application;
FIG. 2 is a block diagram of a system for determining microseismic magnitudes provided in accordance with one embodiment of the present application.
FIG. 3 is a block diagram of a calculation module in a microseismic magnitude determination system provided in accordance with one embodiment of the present application;
FIG. 4 is a block diagram of a first determination module in a system for determining a microseismic magnitude provided in accordance with one embodiment of the present application.
Detailed Description
Reference will now be made in detail to embodiments of the present application, examples of which are illustrated in the accompanying drawings, wherein like or similar reference numerals refer to the same or similar elements or elements having the same or similar function throughout. The embodiments described below with reference to the drawings are exemplary and intended to be used for explaining the present application and should not be construed as limiting the present application.
The application provides a method and a system for determining a microseismic magnitude, wherein the method comprises the following steps: carrying out vector scanning on the micro-seismic target monitoring points to obtain the fracture radiation energy corresponding to the target monitoring points; calculating the lower limit of the signal-to-noise ratio corresponding to the target monitoring point based on the rupture radiation energy corresponding to the target monitoring point; determining the equivalent micro-seismic energy corresponding to the target monitoring point according to the lower limit of the signal-to-noise ratio corresponding to the target monitoring point; and determining the micro-seismic magnitude corresponding to the target monitoring point according to the micro-seismic equivalent energy corresponding to the target monitoring point. According to the technical scheme, the micro-earthquake magnitude can be effectively and accurately determined through the micro-earthquake equivalent energy corresponding to the determined target monitoring point.
A method and system for determining a microseismic magnitude according to an embodiment of the present application will be described with reference to the accompanying drawings.
Example one
Fig. 1 is a flowchart of a method for determining a microseismic magnitude according to an embodiment of the present application, as shown in fig. 1, the method includes:
step 1: and carrying out vector scanning on the micro-seismic target monitoring points to obtain the fracture radiation energy corresponding to the target monitoring points.
In an embodiment of the present disclosure, the calculation formula of the corresponding rupture radiation energy of the target monitoring point is as follows:
in the formula, C k The cracking radiation energy corresponding to a target monitoring point k is obtained, M is the number of stations in a space area corresponding to the target monitoring point k, L is the number of recording sampling points of each time period window of the stations, f ij A recorded signal pointing to a target monitor point k for the ith station jth recorded sample,
P k the total energy of the recorded signal for the target monitoring point k.
Step 2: and calculating the lower limit of the signal-to-noise ratio corresponding to the target monitoring point based on the rupture radiation energy corresponding to the target monitoring point.
In an embodiment of the present disclosure, the step 2 specifically includes:
step 2-1: acquiring rupture radiation energy corresponding to each point in a space region corresponding to the target monitoring point, and determining the value of the minimum rupture radiation energy in the rupture radiation energy corresponding to each point in the space region;
it should be noted that the coordinates corresponding to each point in the spatial region corresponding to the target monitoring point are(x, y, z), wherein x 1 ≤x≤x 2 ，y 1 ≤y≤y 2 ，z 1 ≤z≤z 2 Determining the cracking radiation energy corresponding to each point in the range, and finding out the minimum value C min 。
Step 2-2: determining the difference value between the rupture radiation energy corresponding to the target monitoring point and the minimum rupture radiation energy value;
step 2-3: and respectively taking the difference values as the lower limits of the signal-to-noise ratios corresponding to the target monitoring points.
It should be noted that, in order to ensure the accuracy of the lower limit of the signal-to-noise ratio corresponding to the target monitoring point, the following conditions need to be satisfied:
(1) each seismic station should be at a quantitatively determined quiet location;
(2) using a special instrument suitable for monitoring microseisms, in particular whose detector should have a low natural frequency (e.g.. ltoreq.7 Hz);
(3) scanning superposition must take into account the shear fracture characteristics, generally abandoning the superposition of longitudinal waves of much smaller amplitude, while using transverse waves that carry much greater energy when arriving at a distance;
(4) using a minimum number of scanning stations M that is greater than or equal to a statistical significance;
(5) efficient denoising, including removing interference from surface underground machines, remote seismic, non-monitoring microseismic targets.
And step 3: and determining the equivalent micro-seismic energy corresponding to the target monitoring point according to the lower limit of the signal-to-noise ratio corresponding to the target monitoring point.
In an embodiment of the present disclosure, the step 3 specifically includes:
step 3-1: acquiring a background noise energy value corresponding to the target monitoring point, and determining a unit energy value of the background noise;
wherein the calculation formula of the unit energy value of the background noise is as follows:
P J,k ＝P k /G 2 /V vms 2 /M/T sec
in the formula, P J,k The unit energy value of the background noise corresponding to the target monitoring point kG is the magnification of the seismic instrument, V vms The voltage of the speed unit output by the seismic instrument, M is the number of stations in a space area corresponding to the target monitoring point, T sec Is the scan duration.
Step 3-2: determining the signal-to-noise ratio of unit time length corresponding to the target monitoring point according to the signal-to-noise ratio lower limit corresponding to the target monitoring point;
the calculation formula of the signal-to-noise ratio of the unit time corresponding to the target monitoring point is as follows:
in the formula, R J,k Is the signal-to-noise ratio, R, of unit time length corresponding to the target monitoring point k k To the lower limit of the signal-to-noise ratio, T sec Is the scan duration, and T is the microseismic period.
Step 3-3: and determining the equivalent micro-seismic energy corresponding to the target monitoring point according to the unit energy value of the background noise and the signal-to-noise ratio of unit time length.
Wherein, the calculation formula of the micro-seismic equivalent energy corresponding to the target monitoring point is as follows:
E J,k ＝P J,k ·R J,k
in the formula, E J,k Is the micro-seismic equivalent energy, P, corresponding to the target monitoring point k J,k Is the unit energy value, R, of the background noise corresponding to the target monitoring point k J,k The signal-to-noise ratio of the unit time length corresponding to the target monitoring point k.
And 4, step 4: and determining the micro-seismic magnitude corresponding to the target monitoring point according to the micro-seismic equivalent energy corresponding to the target monitoring point.
In an embodiment of the present disclosure, the step 4 specifically includes:
the determining the micro-seismic magnitude corresponding to the target monitoring point according to the micro-seismic equivalent energy corresponding to the target monitoring point comprises the following steps:
substituting the micro-seismic equivalent energy corresponding to the target monitoring point into a relational expression of an E-Richter magnitude M, and solving the relational expression to obtain the micro-seismic magnitude corresponding to the target monitoring point;
wherein the relation of the E-Richter magnitude M is log 10 E J,k 11.8+1.5A, wherein A is the microseismic magnitude corresponding to the target monitoring point, and the energy unit is erg-10 -7 J。
The method for determining the magnitude of a microseism provided by the embodiment of the present disclosure is specifically described with reference to the above method, and includes:
the method provided by the invention is used for respectively determining the unit energy value and the seismic level of background noise corresponding to the collapse of the coal mine goaf, the pretreatment of oil well hot water injection (the top plate is fractured in dozens of meters of rock mass in advance), the fracture of the coal mine top plate and the fracture of the oil well, and the method is shown in the table 1:
it can be seen that the determined values in the table correspond substantially to the magnitude range of the near monitor contrast in the well, it being noted that the equivalent energy levels here are to equate microseismic effects in the scan domain and time period to a point in space.
The method provided by the implementation can also be applied to geological disaster monitoring of a mining area, and when an event with a small earthquake (the earthquake magnitude A is more than or equal to 0) is monitored, the disaster often occurs; therefore, such monitoring should be developed for early warning of slight shock. The microseismic which most possibly causes underground mines and ground geological disasters at present is a strong microseismic with A being more than or equal to-1.5. The inducement and the occurrence condition of the event A is more than or equal to-1.5 are clarified, and the time-space correlation of strong microseisms and small microseisms is the key point of the microseismic monitoring of the mining area. Meanwhile, the method provided by the embodiment can also be applied to oil and gas production, and an oil and gas reservoir, a production water injection air cavity, a reservoir oil gas front edge, an oil reservoir geological state and a fracturing process state are closely related to microseismic activities including the size of microseismic (group), and the microseismic size needs to be considered in the microseismic monitoring.
In summary, the method for determining the magnitude of a microseism provided by the embodiment includes: carrying out vector scanning on the micro-seismic target monitoring points to obtain the fracture radiation energy corresponding to the target monitoring points; calculating the lower limit of the signal-to-noise ratio corresponding to the target monitoring point based on the rupture radiation energy corresponding to the target monitoring point; determining the equivalent micro-seismic energy corresponding to the target monitoring point according to the lower limit of the signal-to-noise ratio corresponding to the target monitoring point; and determining the micro-seismic magnitude corresponding to the target monitoring point according to the micro-seismic equivalent energy corresponding to the target monitoring point. According to the technical scheme, the micro-earthquake magnitude can be effectively and accurately determined through the micro-earthquake equivalent energy corresponding to the determined target monitoring point.
Example two
Fig. 2 is a block diagram of a system for determining a microseismic magnitude according to an embodiment of the present application, and as shown in fig. 4, the system may include:
the scanning module 100 is used for performing vector scanning on the micro-seismic target monitoring points to obtain fracture radiation energy corresponding to the target monitoring points;
a calculating module 200, configured to calculate, based on the burst radiation energy corresponding to the target monitoring point, a signal-to-noise ratio lower limit corresponding to the target monitoring point;
the first determining module 300 is configured to determine microseism equivalent energy corresponding to the target monitoring point according to a signal-to-noise ratio lower limit corresponding to the target monitoring point;
a second determining module 400, configured to determine, according to the micro-seismic equivalent energy corresponding to the target monitoring point, a micro-seismic magnitude corresponding to the target monitoring point.
Specifically, the calculation formula of the fracture radiation energy corresponding to the target monitoring point is as follows:
in the formula, C k The cracking radiation energy corresponding to a target monitoring point k is obtained, M is the number of stations in a space area corresponding to the target monitoring point k, L is the number of recording sampling points of each time period window of the stations, f ij A recorded signal pointing to a target monitor point k for the ith station jth recorded sample,
P k the total energy of the recorded signal for the target monitoring point k.
In an embodiment of the present disclosure, as shown in fig. 3, the calculating module 200 includes:
a first determining unit 201, configured to obtain the corresponding rupture radiation energy in each point in the spatial region corresponding to the target monitoring point, and determine a value of the minimum rupture radiation energy in the corresponding rupture radiation energy in the spatial region;
the second determining unit 202 is configured to determine a difference between the fracture radiation energy corresponding to the target monitoring point and the minimum fracture radiation energy value;
and the signal-to-noise ratio unit 203 is configured to take the difference values as signal-to-noise ratio lower limits corresponding to the target monitoring points, respectively.
In the disclosed embodiment, as shown in fig. 4, the first determining module 300 includes:
a third determining unit 301, configured to obtain an energy value of background noise corresponding to the target monitoring point, and determine a unit energy value of the background noise;
a fourth determining unit 302, configured to determine a signal-to-noise ratio of a unit duration corresponding to the target monitoring point according to a signal-to-noise ratio lower limit corresponding to the target monitoring point;
a fifth determining unit 303, configured to determine microseismic equivalent energy corresponding to the target monitoring point according to the unit energy value of the background noise and the signal-to-noise ratio of the unit time length.
Further, the calculation formula of the unit energy value of the background noise is as follows:
P J,k ＝P k /G 2 /V vms 2 /M/T sec
in the formula, P J,k Is the unit energy value of the background noise corresponding to the target monitoring point k, G is the magnification factor of the seismic instrument, V vms In units of velocity output by seismic instrumentsVoltage, M is the number of stations in the space region corresponding to the target monitoring point, T sec Is the scan duration.
Further, the calculation formula of the signal-to-noise ratio of the target monitoring point in unit time duration is as follows:
in the formula, R J,k Is the signal-to-noise ratio, R, of unit time length corresponding to the target monitoring point k k To the lower limit of the signal-to-noise ratio, T sec Is the scan duration, and T is the microseismic period.
Further, the calculation formula of the micro-seismic equivalent energy corresponding to the target monitoring point is as follows:
E J,k ＝P J,k ·R J,k
in the formula, E J,k Micro-seismic equivalent energy, P, corresponding to target monitoring point k J,k Is the unit energy value, R, of the background noise corresponding to the target monitoring point k J,k The signal-to-noise ratio of the unit time length corresponding to the target monitoring point k.
In an embodiment of the present disclosure, the second determining module 400 is specifically configured to:
substituting the micro-seismic equivalent energy corresponding to the target monitoring point into a relational expression of an E-Richter magnitude M, and solving the relational expression to obtain the micro-seismic magnitude corresponding to the target monitoring point;
wherein the relation of the E-Richter magnitude M is log 10 E J,k 11.8+1.5A, wherein A is the micro-seismic magnitude corresponding to the target monitoring point.
In summary, the present application provides a system for determining a microseismic magnitude, the system comprising: the scanning module, the calculating module, the first determining module and the second determining module can effectively and accurately determine the earthquake magnitude of the micro earthquake through the determined micro earthquake equivalent energy corresponding to the target monitoring point.
Example 3
In order to implement the above embodiments, the present disclosure also proposes a computer-readable storage medium.
The present embodiment provides a computer device having a computer program stored thereon, wherein the computer program, when executed by a processor, implements the method of embodiment 1.
In the description herein, reference to the description of the term "one embodiment," "some embodiments," "an example," "a specific example," or "some examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the application. In this specification, the schematic representations of the terms used above are not necessarily intended to refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Furthermore, various embodiments or examples and features of different embodiments or examples described in this specification can be combined and combined by one skilled in the art without contradiction.
Any process or method descriptions in flow charts or otherwise described herein may be understood as representing modules, segments, or portions of code which include one or more executable instructions for implementing steps of a custom logic function or process, and alternate implementations are included within the scope of the preferred embodiment of the present application in which functions may be executed out of order from that shown or discussed, including substantially concurrently or in reverse order, depending on the functionality involved, as would be understood by those reasonably skilled in the art of the present application.
Although embodiments of the present application have been shown and described above, it is understood that the above embodiments are exemplary and should not be construed as limiting the present application, and that variations, modifications, substitutions and alterations may be made to the above embodiments by those of ordinary skill in the art within the scope of the present application.
Claims (10)
1. A method of determining a microseismic magnitude, the method comprising:
carrying out vector scanning on the micro-seismic target monitoring points to obtain the fracture radiation energy corresponding to the target monitoring points;
calculating the lower limit of the signal-to-noise ratio corresponding to the target monitoring point based on the rupture radiation energy corresponding to the target monitoring point;
determining the equivalent micro-seismic energy corresponding to the target monitoring point according to the lower limit of the signal-to-noise ratio corresponding to the target monitoring point;
and determining the micro-seismic magnitude corresponding to the target monitoring point according to the micro-seismic equivalent energy corresponding to the target monitoring point.
2. The method of claim 1, wherein the burst radiation energy corresponding to the target monitoring point is calculated as follows:
in the formula, C k The cracking radiation energy corresponding to a target monitoring point k is obtained, M is the number of stations in a space area corresponding to the target monitoring point k, L is the number of recording sampling points of each time period window of the stations, f ij A recorded signal pointing to a target monitor point k for the ith station jth recorded sample,
P k the total energy of the recorded signal for the target monitoring point k.
3. The method of claim 1, wherein said calculating a lower signal-to-noise ratio limit for the target monitoring point based on the corresponding target monitoring point's fracture radiation energy comprises:
acquiring rupture radiation energy corresponding to each point in a space region corresponding to the target monitoring point, and determining the value of the minimum rupture radiation energy in the rupture radiation energy corresponding to each point in the space region;
determining the difference value between the rupture radiation energy corresponding to the target monitoring point and the minimum rupture radiation energy value;
and respectively taking the difference values as the lower limits of the signal-to-noise ratios corresponding to the target monitoring points.
4. The method of claim 1, wherein determining the corresponding microseismic equivalent energy of the target monitoring point according to the corresponding lower signal-to-noise ratio limit of the target monitoring point comprises:
acquiring a background noise energy value corresponding to the target monitoring point, and determining a unit energy value of the background noise;
determining the signal-to-noise ratio of unit time length corresponding to the target monitoring point according to the signal-to-noise ratio lower limit corresponding to the target monitoring point;
and determining the equivalent micro-seismic energy corresponding to the target monitoring point according to the unit energy value of the background noise and the signal-to-noise ratio of unit time length.
5. The determination method of claim 4, wherein the energy per unit value of the background noise is calculated as follows:
P J,k ＝P k /G 2 /V vms 2 /M/T sec
in the formula, P J,k Is the unit energy value of the background noise corresponding to the target monitoring point k, G is the magnification factor of the seismic instrument, V vms The voltage of the speed unit output by the seismic instrument, M is the number of stations in a space area corresponding to the target monitoring point, T sec Is the scan duration.
6. The method of claim 4, wherein the snr per unit time corresponding to the target monitoring point is calculated as follows:
in the formula, R J,k Signal noise of unit time length corresponding to target monitoring point kRatio, R k To the lower limit of the signal-to-noise ratio, T sec Is the scan duration, and T is the microseismic period.
7. The method of claim 4, wherein the microseismic equivalent energy for the target survey point is calculated as follows:
E J,k ＝P J,k ·R J,k
in the formula, E J,k Is the micro-seismic equivalent energy, P, corresponding to the target monitoring point k J,k Is the unit energy value, R, of the background noise corresponding to the target monitoring point k J,k The signal-to-noise ratio of the unit time length corresponding to the target monitoring point k.
8. The method for determining as claimed in claim 7, wherein said determining the micro-seismic magnitude corresponding to the target monitoring point according to the micro-seismic equivalent energy corresponding to the target monitoring point comprises:
substituting the micro-seismic equivalent energy corresponding to the target monitoring point into a relational expression of an E-Richter magnitude M, and solving the relational expression to obtain the micro-seismic magnitude corresponding to the target monitoring point;
wherein the relation of the E-Richter magnitude M is log 10 E J,k 11.8+1.5A, wherein A is the micro-seismic magnitude corresponding to the target monitoring point.
9. A microseismic magnitude determination system, comprising:
the scanning module is used for carrying out vector scanning on the micro-seismic target monitoring points to obtain the fracture radiation energy corresponding to the target monitoring points;
the calculation module is used for calculating the lower limit of the signal-to-noise ratio corresponding to the target monitoring point based on the rupture radiation energy corresponding to the target monitoring point;
the first determining module is used for determining the micro-seismic equivalent energy corresponding to the target monitoring point according to the lower limit of the signal-to-noise ratio corresponding to the target monitoring point;
and the second determination module is used for determining the micro-seismic magnitude corresponding to the target monitoring point according to the micro-seismic equivalent energy corresponding to the target monitoring point.
10. A computer-readable storage medium, on which a computer program is stored, which program, when being executed by a processor, carries out the determination method according to any one of claims 1 to 8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210459414.5A CN114859411A (en) | 2022-04-27 | 2022-04-27 | Method and system for determining microseism magnitude |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210459414.5A CN114859411A (en) | 2022-04-27 | 2022-04-27 | Method and system for determining microseism magnitude |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114859411A true CN114859411A (en) | 2022-08-05 |
Family
ID=82633369
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210459414.5A Pending CN114859411A (en) | 2022-04-27 | 2022-04-27 | Method and system for determining microseism magnitude |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114859411A (en) |
-
2022
- 2022-04-27 CN CN202210459414.5A patent/CN114859411A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Parker et al. | Active‐source seismic tomography at the Brady geothermal field, Nevada, with dense nodal and fiber‐optic seismic arrays | |
| Murphy et al. | Advanced seismic analyses of the source characteristics of the 2006 and 2009 North Korean nuclear tests | |
| Mendecki et al. | A guide to routine seismic monitoring in mines | |
| US20090299637A1 (en) | Continuous Reservoir Monitoring for Fluid Pathways Using Microseismic Data | |
| McNutt et al. | Seismic and infrasonic monitoring | |
| Zeng et al. | Surface microseismic monitoring of hydraulic fracturing of a shale‐gas reservoir using short‐period and broadband seismic sensors | |
| Husen et al. | Induced seismicity during the construction of the Gotthard Base Tunnel, Switzerland: hypocenter locations and source dimensions | |
| Gok et al. | Effect of soil conditions on predicted ground motion: Case study from Western Anatolia, Turkey | |
| NO341717B1 (en) | Stacking of seismic noise data to analyze microseismic events | |
| Murphy et al. | Exploitation of the IMS and other data for a comprehensive, advanced analysis of the North Korean nuclear tests | |
| CN113153430B (en) | Roadway surrounding rock damage acoustic emission positioning and wave velocity imaging monitoring and catastrophe early warning method | |
| Zollo et al. | Space and time behavior of seismic activity at Mt. Vesuvius volcano, southern Italy | |
| Butcher et al. | Evaluating rock mass disturbance within open-pit excavations using seismic methods: A case study from the Hinkley Point C nuclear power station | |
| US20080112263A1 (en) | System and method for determining seismic event location | |
| Dai et al. | Focal mechanism determination for microseismic events and its application to the left bank slope of the Baihetan hydropower station in China | |
| D’Alessandro et al. | Montana Seismic Network Performance: an evaluation through the SNES method | |
| JP2007040949A (en) | System for predicting seismic motion utilizing real-time seismic information | |
| Shashidhar et al. | Relative locations of an earthquake sequence recorded during June 2017 on the Koyna–Warna borehole seismic network of western India | |
| CN114859411A (en) | Method and system for determining microseism magnitude | |
| Lim et al. | Horizontal seismic wave at ground surface from transfer function based on ambient noise | |
| Kong et al. | Evaluation of earthquake signal characteristics for early warning | |
| Macpherson et al. | Using local infrasound to estimate seismic velocity and earthquake magnitudes | |
| Zollo et al. | Source characterization for earthquake early warning | |
| Wang et al. | Magnitude-scaling relationships based on initial P-wave information in the Xinjiang region, China | |
| Peng et al. | A τ c magnitude estimation of the 20 April 2013 Lushan earthquake, Sichuan, China |
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 |