A method and system for determining microseismic magnitude

Abstract

The present application proposes a method and system for determining the magnitude of a microseism, the method comprising: performing vector scanning on a microseismic target monitoring point to obtain the rupture radiation energy corresponding to the target monitoring point; 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 microseismic 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 determining the microseismic magnitude corresponding to the target monitoring point according to the microseismic equivalent energy corresponding to the target monitoring point. The technical solution proposed in the present application can effectively and accurately determine the magnitude of a microseism by determining the microseismic equivalent energy corresponding to the target monitoring point.

Description

Technical Field

[0001] The present application relates to the field of microseismic monitoring technology, and in particular to a method and system for determining the magnitude of a microseismic event. Background Art

[0002] There are already mature methods for determining the magnitude of significant visible earthquake events in general earthquake records, such as the Richter magnitude, moment magnitude, and the correlation between magnitudes with different definitions, the relationship between magnitude and energy, etc. These magnitudes are determined based on the energy released by the earthquake, or the energy is proportional to the square of the maximum amplitude of the recorded earthquake event or the moment of fault dislocation; this requires specific measurements of the amplitude, period, dislocation, etc. of the recorded earthquake event, and the magnitude is calculated based on the characteristics of the recording instrument. The Richter magnitude has been used in the earthquake community for nearly a hundred years.

[0003] However, since microseismic events are so small, only when one is close to the monitoring target can the prominent amplitude of visible earthquake events be recorded, and this kind of proximity monitoring is often difficult. In order to expand the monitoring range of the microseismic monitoring network, a large number of microseismic monitoring stations are more than a few hundred meters away from the earthquake source, and it is common to reach two or three kilometers. In this case, the microseismic event records are either completely submerged in the background noise, or even if they are barely visible, the microseismic magnitude cannot be effectively and accurately determined due to the low signal-to-noise ratio. Therefore, we urgently need to find a method to determine the magnitude of microseismic events that can communicate with the existing Richter magnitude estimation system. Summary of the invention

[0004] The present application provides a method and system for determining the magnitude of a microseism, so as to at least solve the technical problem in the related art that the magnitude of a microseism cannot be determined effectively and accurately.

[0005] The first embodiment of the present application provides a method for determining the magnitude of a microseismic event, the method comprising:

[0006] Performing vector scanning on the microseismic target monitoring point to obtain the rupture radiation energy corresponding to the target monitoring point;

[0007] Calculating a lower limit of a signal-to-noise ratio corresponding to the target monitoring point based on the rupture radiation energy corresponding to the target monitoring point;

[0008] Determine the microseismic 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;

[0009] The microseismic magnitude corresponding to the target monitoring point is determined according to the microseismic equivalent energy corresponding to the target monitoring point.

[0010] The second embodiment of the present application provides a system for determining the magnitude of a microseismic event, the system comprising:

[0011] A scanning module is used to perform vector scanning on a microseismic target monitoring point to obtain the rupture radiation energy corresponding to the target monitoring point;

[0012] A calculation module, used for calculating a lower limit of a signal-to-noise ratio corresponding to the target monitoring point based on the rupture radiation energy corresponding to the target monitoring point;

[0013] A first determination module is used to determine the microseismic 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;

[0014] The second determination module is used to determine the microseismic magnitude corresponding to the target monitoring point according to the microseismic equivalent energy corresponding to the target monitoring point.

[0015] The third aspect of the present application proposes a computer-readable storage medium, on which a computer program is stored. When the program is executed by a processor, the method described in the first aspect is implemented.

[0016] The technical solution provided by the embodiments of the present application brings at least the following beneficial effects:

[0017] The present application proposes a method and system for determining the magnitude of a microseism, wherein the method includes: performing vector scanning on a microseismic target monitoring point to obtain the rupture radiation energy corresponding to the target monitoring point; 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 microseismic 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 determining the microseismic magnitude corresponding to the target monitoring point according to the microseismic equivalent energy corresponding to the target monitoring point. The technical solution proposed in the present application can effectively and accurately determine the magnitude of a microseism by determining the microseismic equivalent energy corresponding to the target monitoring point.

[0018] Additional aspects and advantages of the present application will be given in part in the description below, and in part will become apparent from the description below, or will be learned through the practice of the present application. BRIEF DESCRIPTION OF THE DRAWINGS

[0019] The above and / or additional aspects and advantages of the present application will become apparent and easily understood from the following description of the embodiments in conjunction with the accompanying drawings, in which:

[0020] Figure 1 A flowchart of a method for determining a microseismic magnitude according to an embodiment of the present application;

[0021] Figure 2 A structural diagram of a system for determining microseismic magnitude according to an embodiment of the present application.

[0022] Figure 3A structural diagram of a calculation module in a system for determining microseismic magnitude according to an embodiment of the present application;

[0023] Figure 4 This is a structural diagram of a first determination module in a system for determining microseismic magnitude according to an embodiment of the present application. DETAILED DESCRIPTION

[0024] Embodiments of the present application are described in detail below, and examples of the embodiments are shown in the accompanying drawings, wherein the same or similar reference numerals throughout represent the same or similar elements or elements having the same or similar functions. The embodiments described below with reference to the accompanying drawings are exemplary and are intended to be used to explain the present application, and should not be construed as limiting the present application.

[0025] The present application proposes a method and system for determining the magnitude of a microseism, wherein the method includes: performing vector scanning on a microseismic target monitoring point to obtain the rupture radiation energy corresponding to the target monitoring point; 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 microseismic 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 determining the microseismic magnitude corresponding to the target monitoring point according to the microseismic equivalent energy corresponding to the target monitoring point. The technical solution proposed in the present application can effectively and accurately determine the magnitude of a microseism by determining the microseismic equivalent energy corresponding to the target monitoring point.

[0026] A method and system for determining the magnitude of a microseism in an embodiment of the present application will be described below with reference to the accompanying drawings.

[0027] Embodiment 1

[0028] Figure 1 A flowchart of a method for determining a microseismic magnitude according to an embodiment of the present application is shown in FIG. Figure 1 As shown, the method includes:

[0029] Step 1: Perform vector scanning on the microseismic target monitoring point to obtain the rupture radiation energy corresponding to the target monitoring point.

[0030] In the embodiment of the present disclosure, the calculation formula of the rupture radiation energy corresponding to the target monitoring point is as follows:

[0031]

[0032] In the formula, C k is the rupture radiation energy corresponding to the target monitoring point k, M is the number of stations in the spatial area corresponding to the target monitoring point k, L is the number of recorded sample points in each time window of the station, and f ijis the recorded signal of the jth recorded sample point at the i-th station pointing to the target monitoring point k, P k is the total energy of the recorded signal of the target monitoring point k.

[0033] Step 2: Calculate 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.

[0034] In the embodiment of the present disclosure, step 2 specifically includes:

[0035] Step 2-1: obtaining the rupture radiation energy corresponding to each point in the spatial region corresponding to the target monitoring point, and determining the minimum rupture radiation energy value among the rupture radiation energies corresponding to each point in the spatial region;

[0036] It should be noted that the coordinates of each point in the spatial area corresponding to the target monitoring point are (x, y, z), where x1≤x≤x2, y1≤y≤y2, z1≤z≤z2. The rupture radiation energy corresponding to each point is determined within this range, and the minimum value C is found. min .

[0037] Step 2-2: determining the difference between the rupture radiation energy corresponding to the target monitoring point and the minimum rupture radiation energy value;

[0038] Step 2-3: The difference is used as the lower limit of the signal-to-noise ratio corresponding to the target monitoring point.

[0039] It should be noted that in order to ensure the accuracy of the lower limit of the signal-to-noise ratio corresponding to the above target monitoring points, the following conditions must be met:

[0040] (1) Each seismic station should be located in a quantitatively determined quiet location;

[0041] (2) Use specialized instruments suitable for monitoring microseisms, especially geophones with low natural frequencies (e.g. ≤7 Hz);

[0042] (3) Scanning stacking must take into account the shear fracture characteristics. Generally, the longitudinal waves with much smaller amplitudes are abandoned and the shear waves with much greater energy are used when reaching far distances.

[0043] (4) Use a minimum number of scanning stations M that is greater than or equal to a statistically significant number;

[0044] (5) Effective denoising, including removing interference from ground and underground machinery, teleseismic, and non-monitored microseismic targets.

[0045] Step 3: Determine the microseismic 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.

[0046] In the embodiment of the present disclosure, step 3 specifically includes:

[0047] Step 3-1: Obtain the background noise energy value corresponding to the target monitoring point, and determine the unit energy value of the background noise;

[0048] The calculation formula of the unit energy value of the background noise is as follows:

[0049] P J,k =P k / G 2 / V vms 2 / M / T sec

[0050] Where P J,k is the unit energy value of the background noise corresponding to the target monitoring point k, G is the magnification of the seismic instrument, V vms is the voltage in velocity units output by the seismic instrument, M is the number of stations in the spatial area corresponding to the target monitoring point, T sec The scan duration.

[0051] Step 3-2: determining the signal-to-noise ratio per unit time length corresponding to the target monitoring point according to the signal-to-noise ratio lower limit corresponding to the target monitoring point;

[0052] The calculation formula of the signal-to-noise ratio per unit time corresponding to the target monitoring point is as follows:

[0053]

[0054] In the formula, R J,k is the signal-to-noise ratio per unit time corresponding to the target monitoring point k, R k is the lower limit of the signal-to-noise ratio, T sec is the scanning time, and T is the microseismic period.

[0055] Step 3-3: Determine the 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 per unit time length.

[0056] The calculation formula of the microseismic equivalent energy corresponding to the target monitoring point is as follows:

[0057] E J,k =P J,k ·R J,k

[0058] In the formula, E J,k is the microseismic equivalent energy corresponding to the target monitoring point k, P J,k is the unit energy value of the background noise corresponding to the target monitoring point k, R J,k is the signal-to-noise ratio per unit time corresponding to the target monitoring point k.

[0059] Step 4: Determine the microseismic magnitude corresponding to the target monitoring point according to the microseismic equivalent energy corresponding to the target monitoring point.

[0060] In the embodiment of the present disclosure, step 4 specifically includes:

[0061] Determining the microseismic magnitude corresponding to the target monitoring point according to the microseismic equivalent energy corresponding to the target monitoring point includes:

[0062] Substituting the microseismic equivalent energy corresponding to the target monitoring point into the relationship between E and Richter magnitude M, and solving the relationship to obtain the microseismic magnitude corresponding to the target monitoring point;

[0063] The relationship between E and Richter magnitude M is log 10 E J,k =11.8+1.5A, where A is the microseismic magnitude corresponding to the target monitoring point, and the energy unit is erg=10 -7 J.

[0064] In combination with the above method, a method for determining the magnitude of a microseismic event provided by an embodiment of the present disclosure is specifically described, including:

[0065] The method provided by the present invention is used to determine the unit energy value and magnitude of the background noise corresponding to the collapse of the coal mine goaf, the collapse of the coal mine goaf, the pretreatment of hot water injection of the oil well (the rock mass of several tens of meters of the roof has been fractured in advance), the fracturing of the coal mine roof and the fracturing of the oil well, as shown in Table 1:

[0066]

[0067] It can be concluded that the determined values ​​in the table basically correspond to the magnitude range of the adjacent monitoring comparison in the well. It should be noted that the equivalent energy level here is to equate the microseismic effects within the scanning domain and time period to a point in space.

[0068] The method provided in this embodiment can also be used in the monitoring of geological disasters in mining areas. When an event above a small earthquake (magnitude A≥0) is monitored, it is often a disaster; therefore, this type of monitoring should be developed into microseismic early warning. At present, the microseisms that are most likely to cause underground mines and ground geological disasters are strong microseisms of A≥-1.5. It is the focus of microseismic monitoring in mining areas to clarify the causes and conditions for the occurrence of events A≥-1.5, and the temporal and spatial correlation between strong microseisms and smaller microseisms. At the same time, the method provided in this embodiment can also be applied to oil and gas production. Oil and gas reservoirs, production water injection cavities and oil and gas fronts, the geological state of oil reservoirs, and the state of the fracturing process are all closely related to microseismic activities, including the size of microseismic (groups), and the size of microseismic needs to be statistically considered in these microseismic monitoring.

[0069] In summary, the present embodiment provides a method for determining the magnitude of a microseism, the method comprising: performing vector scanning on a microseismic target monitoring point to obtain the rupture radiation energy corresponding to the target monitoring point; 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 microseismic 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 determining the microseismic magnitude corresponding to the target monitoring point according to the microseismic equivalent energy corresponding to the target monitoring point. The technical solution proposed in the present application can effectively and accurately determine the magnitude of a microseism by determining the microseismic equivalent energy corresponding to the target monitoring point.

[0070] Embodiment 2

[0071] Figure 2 FIG. 1 is a structural diagram of a system for determining microseismic magnitude according to an embodiment of the present application, such as Figure 4 As shown, the system may include:

[0072] The scanning module 100 is used to perform vector scanning on the microseismic target monitoring point to obtain the rupture radiation energy corresponding to the target monitoring point;

[0073] A calculation module 200, configured to calculate a lower limit of a signal-to-noise ratio corresponding to the target monitoring point based on the rupture radiation energy corresponding to the target monitoring point;

[0074] A first determination module 300 is used to determine the microseismic 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;

[0075] The second determination module 400 is used to determine the microseismic magnitude corresponding to the target monitoring point according to the microseismic equivalent energy corresponding to the target monitoring point.

[0076] Specifically, the calculation formula of the rupture radiation energy corresponding to the target monitoring point is as follows:

[0077]

[0078] In the formula, C k is the rupture radiation energy corresponding to the target monitoring point k, M is the number of stations in the spatial area corresponding to the target monitoring point k, L is the number of recorded sample points in each time window of the station, and f ij is the recorded signal of the jth recorded sample point at the i-th station pointing to the target monitoring point k, P k is the total energy of the recorded signal of the target monitoring point k.

[0079] In the embodiments of the present disclosure, Figure 3 As shown, the calculation module 200 includes:

[0080] The first determining unit 201 is used to obtain the fracture radiation energy corresponding to each point in the spatial region corresponding to the target monitoring point, and determine the minimum fracture radiation energy value among the fracture radiation energies corresponding to each point in the spatial region;

[0081] A second determining unit 202 is used to determine the difference between the rupture radiation energy corresponding to the target monitoring point and the minimum rupture radiation energy value;

[0082] The signal-to-noise ratio unit 203 is used to use the difference as the lower limit of the signal-to-noise ratio corresponding to the target monitoring point.

[0083] In the embodiments of the present disclosure, Figure 4 As shown, the first determination module 300 includes:

[0084] The third determining unit 301 is used to obtain the background noise energy value corresponding to the target monitoring point and determine the unit energy value of the background noise;

[0085] A fourth determining unit 302 is configured to determine the signal-to-noise ratio per unit time length corresponding to the target monitoring point according to the signal-to-noise ratio lower limit corresponding to the target monitoring point;

[0086] The fifth determining unit 303 is used to determine the 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 per unit time length.

[0087] Furthermore, the calculation formula of the unit energy value of the background noise is as follows:

[0088] P J,k =P k / G 2 / V vms 2 / M / T sec

[0089] Where P J,k is the unit energy value of the background noise corresponding to the target monitoring point k, G is the magnification of the seismic instrument, V vms is the voltage in velocity units output by the seismic instrument, M is the number of stations in the spatial area corresponding to the target monitoring point, T sec The scan duration.

[0090] Furthermore, the calculation formula of the signal-to-noise ratio per unit time corresponding to the target monitoring point is as follows:

[0091]

[0092] In the formula, R J,k is the signal-to-noise ratio per unit time corresponding to the target monitoring point k, Rk is the lower limit of the signal-to-noise ratio, T sec is the scanning time, and T is the microseismic period.

[0093] Furthermore, the calculation formula of the microseismic equivalent energy corresponding to the target monitoring point is as follows:

[0094] E J,k =P J,k ·R J,k

[0095] In the formula, E J,k is the microseismic equivalent energy corresponding to the target monitoring point k, P J,k is the unit energy value of the background noise corresponding to the target monitoring point k, R J,k is the signal-to-noise ratio per unit time corresponding to the target monitoring point k.

[0096] In the embodiment of the present disclosure, the second determining module 400 is specifically used to:

[0097] Substituting the microseismic equivalent energy corresponding to the target monitoring point into the relationship between E and Richter magnitude M, and solving the relationship to obtain the microseismic magnitude corresponding to the target monitoring point;

[0098] The relationship between E and Richter magnitude M is log 10 E J,k =11.8+1.5A, where A is the microseismic magnitude corresponding to the target monitoring point.

[0099] In summary, the present application proposes a system for determining the magnitude of a microseism, which includes: a scanning module, a calculation module, a first determination module and a second determination module. The microseismic magnitude can be effectively and accurately determined by determining the microseismic equivalent energy corresponding to the target monitoring point.

[0100] Example 3

[0101] In order to implement the above embodiments, the present disclosure also proposes a computer-readable storage medium.

[0102] The computer device provided in this embodiment stores a computer program, and when the computer program is executed by a processor, the method in Embodiment 1 is implemented.

[0103] In the description of this specification, the description with reference to the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" etc. means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present application. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or more embodiments or examples in a suitable manner. In addition, those skilled in the art may combine and combine the different embodiments or examples described in this specification and the features of the different embodiments or examples, without contradiction.

[0104] Any process or method description in a flowchart or otherwise described herein may be understood to represent a module, fragment or portion of code comprising one or more executable instructions for implementing the steps of a custom logical function or process, and the scope of the preferred embodiments of the present application includes alternative implementations in which functions may not be performed in the order shown or discussed, including performing functions in a substantially simultaneous manner or in the reverse order depending on the functions involved, which should be understood by technicians in the technical field to which the embodiments of the present application belong.

[0105] Although the embodiments of the present application have been shown and described above, it can be understood that the above embodiments are exemplary and cannot be understood as limitations on the present application. Ordinary technicians in this field can change, modify, replace and modify the above embodiments within the scope of the present application.

Claims

1. A method for determining the magnitude of a microseismic event, characterized in that: The method comprises: Performing vector scanning on the microseismic target monitoring point to obtain the rupture radiation energy corresponding to the target monitoring point; Calculating a lower limit of a signal-to-noise ratio corresponding to the target monitoring point based on the rupture radiation energy corresponding to the target monitoring point; Determine the microseismic 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; Determine the microseismic magnitude corresponding to the target monitoring point according to the microseismic equivalent energy corresponding to the target monitoring point; Wherein, determining the microseismic 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 includes: Obtaining a background noise energy value corresponding to the target monitoring point, and determining a unit energy value of the background noise; Determine the signal-to-noise ratio per unit time length corresponding to the target monitoring point according to the signal-to-noise ratio lower limit corresponding to the target monitoring point; Determine the microseismic equivalent energy corresponding to the target monitoring point according to the product of the unit energy value of the background noise and the signal-to-noise ratio per unit time length; The 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 includes: Acquire the rupture radiation energy corresponding to each point in the spatial region corresponding to the target monitoring point, and determine the minimum rupture radiation energy value among the rupture radiation energies corresponding to each point in the spatial region; Determine the difference between the rupture radiation energy corresponding to the target monitoring point and the minimum rupture radiation energy value; The differences are respectively used as the lower limits of the signal-to-noise ratio corresponding to the target monitoring points.

2. The determination method according to claim 1, characterized in that: The calculation formula of the rupture radiation energy corresponding to the target monitoring point is as follows: In the formula, C k is the rupture radiation energy corresponding to the target monitoring point k, M is the number of stations in the spatial area corresponding to the target monitoring point k, L is the number of recorded sample points in each time window of the station, and f ij is the recorded signal of the jth recorded sample point at the i-th station pointing to the target monitoring point k, P k is the total energy of the recorded signal of the target monitoring point k.

3. The determination method according to claim 1, characterized in that: 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 Where P J,k is the unit energy value of the background noise corresponding to the target monitoring point k, G is the magnification of the seismic instrument, V vms is the voltage in velocity units output by the seismic instrument, M is the number of stations in the spatial area corresponding to the target monitoring point, T sec is the scanning time, P k is the total energy of the recorded signal of the target monitoring point k.

4. The determination method according to claim 1, characterized in that: The calculation formula of the signal-to-noise ratio per unit time corresponding to the target monitoring point is as follows: In the formula, R J,k is the signal-to-noise ratio per unit time corresponding to the target monitoring point k, R k is the lower limit of the signal-to-noise ratio, T sec is the scanning time, and T is the microseismic period.

5. The determination method according to claim 1, characterized in that: The calculation formula of the microseismic 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 microseismic equivalent energy corresponding to the target monitoring point k, P J,k is the unit energy value of the background noise corresponding to the target monitoring point k, R J,k is the signal-to-noise ratio per unit time corresponding to the target monitoring point k.

6. The determination method according to claim 5, characterized in that: Determining the microseismic magnitude corresponding to the target monitoring point according to the microseismic equivalent energy corresponding to the target monitoring point includes: Substituting the microseismic equivalent energy corresponding to the target monitoring point into the relationship between E and Richter magnitude M, and solving the relationship to obtain the microseismic magnitude corresponding to the target monitoring point; The relationship between E and Richter magnitude M is log 10 E J,k =11.8+1.5A, where A is the microseismic magnitude corresponding to the target monitoring point.

7. A system for determining microseismic magnitude, characterized in that: The system comprises: A scanning module is used to perform vector scanning on a microseismic target monitoring point to obtain the rupture radiation energy corresponding to the target monitoring point; A calculation module, used for calculating a lower limit of a signal-to-noise ratio corresponding to the target monitoring point based on the rupture radiation energy corresponding to the target monitoring point; A first determination module is used to determine the microseismic 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; A second determination module is used to determine the microseismic magnitude corresponding to the target monitoring point according to the microseismic equivalent energy corresponding to the target monitoring point; Wherein, determining the microseismic 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 includes: Obtaining a background noise energy value corresponding to the target monitoring point, and determining a unit energy value of the background noise; Determine the signal-to-noise ratio per unit time length corresponding to the target monitoring point according to the signal-to-noise ratio lower limit corresponding to the target monitoring point; Determine the microseismic equivalent energy corresponding to the target monitoring point according to the product of the unit energy value of the background noise and the signal-to-noise ratio per unit time length; The 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 includes: Acquire the rupture radiation energy corresponding to each point in the spatial region corresponding to the target monitoring point, and determine the minimum rupture radiation energy value among the rupture radiation energies corresponding to each point in the spatial region; Determine the difference between the rupture radiation energy corresponding to the target monitoring point and the minimum rupture radiation energy value; The differences are respectively used as the lower limits of the signal-to-noise ratio corresponding to the target monitoring points.

8. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the program is executed by a processor, the determination method according to any one of claims 1 to 6 is implemented.

Images