Microseismic event group positioning method and device, equipment and storage medium

Abstract

The invention discloses a micro-seismic event group positioning method, device and equipment and a storage medium, and relates to the technical field of tunnel engineering, and the method comprises the steps: obtaining a micro-seismic event set, and carrying out the pairing of micro-seismic events in the micro-seismic event set, and obtaining a micro-seismic event pair; clustering the micro-seismic event pairs to obtain micro-seismic event classes; when the number of the micro-seismic event classes is not a preset value, calculating the source spacing of each micro-seismic event pair in the micro-seismic event classes; and according to the source spacing, micro-seismic event group positioning is carried out. According to the invention, the influence of sensor array arrangement can be eliminated, and the positioning precision of the micro-seismic event group is improved.

Description

Technical Field

[0001] The present application relates to the field of tunnel engineering technology, and in particular to a method, device, equipment and storage medium for locating a microseismic event group. Background Art

[0002] As the number of deep tunnel projects increases, early warning of rock bursts is becoming more and more important, because rock burst disasters pose a serious threat to the lives of construction workers. Microseismic monitoring technology is currently the most effective means of early warning of rock bursts, among which microseismic event location is an important analysis step in microseismic monitoring and is closely related to the accuracy of rock burst warning.

[0003] At present, the methods widely used for microseismic event location are all based on the time difference principle, and the location of a single microseismic event is obtained by solving the time difference equation through different optimizations, such as the particle swarm method and the Newton downhill method. However, the time difference method positioning will be affected by the arrangement of the sensor array. However, in linear projects such as tunnels, the arrangement of the sensor array is deformed, and the waveform reception in the Z-axis direction is seriously insufficient, which will inevitably reduce the positioning accuracy of the time difference method in tunnels.

[0004] Therefore, how to improve the accuracy of microseismic event positioning is a problem that still needs to be solved.

[0005] The above contents are only used to assist in understanding the technical solution of the present invention and do not constitute an admission that the above contents are prior art. Summary of the invention

[0006] The main purpose of the present application is to provide a method, device, equipment and storage medium for locating a microseismic event group, aiming to solve the technical problem of how to improve the accuracy of microseismic event locating.

[0007] To achieve the above objectives, the present application proposes a method for locating a microseismic event group, the method comprising:

[0008] Acquire a microseismic event set, and pair the microseismic events in the microseismic event set to obtain microseismic event pairs;

[0009] Clustering the microseismic event pairs to obtain microseismic event classes;

[0010] When the number of the microseismic event class is not a preset value, calculating the source spacing of each pair of microseismic events in the microseismic event class;

[0011] According to the source spacing, the microseismic event group is located.

[0012] In one embodiment, when the number of the microseismic event class is not a preset value, the step of calculating the source spacing of each pair of microseismic event pairs in the microseismic event class includes:

[0013] When the number of the microseismic event class is not a preset value, performing secondary pairing on the microseismic events in the microseismic event class to obtain a target microseismic event pair;

[0014] Calculating the amplitude mean and mean square angular frequency of the target microseismic signal of the target microseismic event pair;

[0015] Obtaining a target amplitude segment according to the amplitude average value;

[0016] calculating the maximum correlation of microseismic event pairs in the microseismic event class according to the target amplitude segment;

[0017] The source spacing of the microseismic event pair is calculated based on the mean square angular frequency and the maximum correlation.

[0018] In one embodiment, the step of locating a microseismic event group according to the source spacing includes:

[0019] According to the source spacing, a source spacing mean and a source spacing standard deviation are calculated;

[0020] Calculating the relative position of the microseismic event group according to the mean value of the source spacing and the standard deviation of the source spacing;

[0021] Obtaining the actual positions, translational positions, circumcenter coordinates, and rotation angles of a preset number of microseismic events;

[0022] The microseismic event group is located according to the relative position, the actual position, the translation position, the circumcenter coordinates and the rotation angle.

[0023] In one embodiment, the step of acquiring a microseismic event set and pairing the microseismic events in the microseismic event set to obtain microseismic event pairs includes:

[0024] Acquire a microseismic event set, and obtain microseismic events and target sensors that record all microseismic events according to the microseismic event set;

[0025] According to the target sensor, the initial wave arrival time of the original microseismic signal is obtained;

[0026] According to the initial wave arrival time, the original microseismic signals are aligned to obtain a target microseismic signal;

[0027] The microseismic events are paired according to the target microseismic signals to obtain microseismic event pairs.

[0028] In one embodiment, the step of clustering the microseismic event pairs to obtain microseismic event classes includes:

[0029] Calculating the correlation of the microseismic event pairs;

[0030] According to the magnitude of the correlation, the microseismic events are clustered to obtain microseismic event classes.

[0031] In one embodiment, the step of calculating the correlation of the microseismic event pairs includes:

[0032] Obtaining a target microseismic signal according to the microseismic event pair;

[0033] Dividing the target microseismic signal into a plurality of calculation windows according to a preset division rule;

[0034] The correlation between the microseismic event pairs is calculated through the multiple calculation windows.

[0035] In one embodiment, after the step of locating a microseismic event group according to the source spacing, the method further includes:

[0036] A preset number of microseismic events are selected from the microseismic event group and added to the microseismic event class, and the step of calculating the source spacing of each pair of microseismic events in the microseismic event class when the number of microseismic events in the microseismic event class is not a preset value is returned.

[0037] In addition, to achieve the above-mentioned purpose, the present application also proposes a microseismic event group positioning device, the microseismic event group positioning device comprising:

[0038] A pairing module, used to obtain a microseismic event set, and pair the microseismic events in the microseismic event set to obtain microseismic event pairs;

[0039] A clustering module clusters the microseismic event pairs to obtain microseismic event classes;

[0040] A calculation module, when the number of the microseismic event class is not a preset value, calculates the source spacing of each pair of microseismic event pairs in the microseismic event class;

[0041] The positioning module locates the microseismic event group according to the source spacing.

[0042] In addition, to achieve the above-mentioned purpose, the present application also proposes a microseismic event group locating device, which includes: a memory, a processor, and a computer program stored in the memory and executable on the processor, and the computer program is configured to implement the steps of the microseismic event group locating method as described above.

[0043] In addition, to achieve the above-mentioned purpose, the present application also proposes a storage medium, which is a computer-readable storage medium, and a computer program is stored on the storage medium. When the computer program is executed by a processor, the steps of the microseismic event group locating method as described above are implemented.

[0044] In addition, to achieve the above-mentioned purpose, the present application also provides a computer program product, which includes a computer program, and when the computer program is executed by a processor, the steps of the microseismic event group positioning method described above are implemented.

[0045] The present invention obtains a microseismic event set, pairs the microseismic events in the microseismic event set, obtains microseismic event pairs, and then clusters the microseismic event pairs to obtain microseismic event classes. When the number of microseismic event classes is not a preset value, the source spacing of each pair of microseismic event pairs in the microseismic event class is calculated, and the microseismic event group is located according to the source spacing. Compared with the existing time difference method positioning, the present application can be unaffected by the arrangement of the sensor array, and can perform microseismic event group positioning, thereby improving the positioning accuracy. BRIEF DESCRIPTION OF THE DRAWINGS

[0046] The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the present application and, together with the description, serve to explain the principles of the present application.

[0047] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, for ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative labor.

[0048] Figure 1 A schematic diagram of a flow chart provided for the first embodiment of the method for locating a microseismic event group of the present application;

[0049] Figure 2 A schematic diagram of a flow chart provided for the second embodiment of the method for locating a microseismic event group of the present application;

[0050] Figure 3 This is a schematic diagram of dividing windows for microseismic events in the second embodiment of the method for locating a microseismic event group of the present application;

[0051] Figure 4 A schematic diagram of a flow chart provided for Embodiment 3 of the method for locating a microseismic event group of the present application;

[0052] Figure 5 Schematic diagram of projections of the relative positions of a microseismic event group in three different directions in a three-dimensional coordinate system in the third embodiment of the microseismic event group positioning method of the present application;

[0053] Figure 6 This is a schematic diagram of positioning results after the overall positioning results of a microseismic event group are translated in the third embodiment of the microseismic event group positioning method of the present application;

[0054] Figure 7This is a schematic diagram of the absolute position results of a microseismic event group in the third embodiment of the microseismic event group positioning method of the present application;

[0055] Figure 8 A schematic diagram of a flow chart provided for Embodiment 4 of the method for locating a microseismic event group of the present application;

[0056] Fig. 9 This is a schematic diagram of the microseismic signal alignment process of two microseismic events in the fourth embodiment of the microseismic event group positioning method of the present application;

[0057] Fig.10 This is a schematic diagram of a microseismic event pairing mode in the fourth embodiment of the microseismic event group positioning method of the present application;

[0058] Fig.11 A flowchart diagram of the fifth embodiment of the method for locating a microseismic event group of the present application is provided;

[0059] Fig.12 This is a schematic diagram of the correlation results of each event pair calculated by different sensors in the fifth embodiment of the microseismic event group positioning method of the present application;

[0060] Fig.13 This is a brief flowchart of an embodiment of a method for locating a microseismic event group of the present application;

[0061] Fig.14 This is a schematic diagram of the module structure of a microseismic event group positioning device according to an embodiment of the present application;

[0062] Fig.15 Schematic diagram of the device structure of the hardware operating environment involved in the microseismic event group positioning method in the embodiment of the present application.

[0063] The purpose, features and advantages of this application will be further described in conjunction with the embodiments and with reference to the accompanying drawings. DETAILED DESCRIPTION

[0064] It should be understood that the specific embodiments described herein are only used to explain the technical solutions of the present application and are not used to limit the present application.

[0065] In order to better understand the technical solution of the present application, a detailed description will be given below in conjunction with the accompanying drawings and specific implementation methods.

[0066] The main solution of the embodiment of the present application is: to obtain a microseismic event set, and pair the microseismic events in the microseismic event set to obtain microseismic event pairs, and then cluster the microseismic event pairs to obtain microseismic event classes, and when the number of the microseismic event classes is not a preset value, calculate the source spacing of each pair of microseismic event pairs in the microseismic event class, and locate the microseismic event group based on the source spacing.

[0067] At present, the methods widely used for microseismic event location are all based on the time difference principle, and the location of a single microseismic event is obtained by solving the time difference equation through different optimizations, such as particle swarm method, Newton downhill method, etc. However, the time difference method positioning will be affected by the arrangement of the sensor array, especially in linear projects such as tunnels, where the arrangement of the sensor array is deformed and the waveform reception in the Z-axis direction is seriously insufficient, which will inevitably reduce the positioning accuracy of the time difference method in tunnels.

[0068] The present application provides a solution that can locate a group of microseismic events without being affected by the arrangement of the sensor array, thereby improving the accuracy of locating the group of microseismic events.

[0069] Based on this, the present application embodiment provides a method for locating a microseismic event group, referring to Figure 1 , Figure 1 This is a flow chart of the first embodiment of the microseismic event group positioning method of the present application.

[0070] In this embodiment, the microseismic event group location method includes steps S10 to S40:

[0071] Step S10: obtaining a microseismic event set, and pairing the microseismic events in the microseismic event set to obtain microseismic event pairs;

[0072] It should be noted that the execution subject of this embodiment may be a computing service device with data processing, network communication and program running functions, such as a tablet computer, a personal computer, a mobile phone, etc., or an electronic device capable of realizing the above functions, a microseismic event group positioning device, etc. The following takes the microseismic event group positioning device as an example to illustrate this embodiment and the following embodiments.

[0073] It should be noted that the microseismic event set refers to a collection of tiny vibration phenomena caused by rock fracture or fluid disturbance.

[0074] It is understandable that the main purpose of pairing microseismic events is to identify microseismic events with similar characteristics or originating from the same geological process. Pairing methods include waveform matching method based on the similarity of microseismic waveform characteristics and arrival time difference method based on the time difference of microseismic events arriving at different monitoring stations.

[0075] Step S20: clustering the microseismic event pairs to obtain microseismic event classes;

[0076] It should be noted that there are many ways to cluster microseismic events, such as hierarchical clustering methods, density-based clustering methods, and grid-based clustering methods. Clustering microseismic event pairs is an important step in locating microseismic event groups. A microseismic event class is a collection of multiple microseismic events divided according to a certain feature.

[0077] Step S30: when the number of the microseismic event class is not a preset value, calculating the source spacing of each pair of microseismic event pairs in the microseismic event class;

[0078] It should be noted that the preset value is a value set according to the actual situation, which can be 1, 2 or other values; the source spacing refers to the distance between the sources of two microseismic events, which is an important indicator for evaluating the spatial distribution and mutual proximity of microseismic events. There are many methods for calculating the source spacing, such as using waveform matching and similarity measurement methods to indirectly estimate the source spacing, and simulating seismic waves through numerical simulation methods to obtain the source spacing.

[0079] Step S40: locating the microseismic event group according to the source spacing.

[0080] It should be noted that the actual positions of some microseismic events in the acquired microseismic event set are known. When the actual positions of some microseismic events are known, the positions of all microseismic events can be calculated through the source spacing of microseismic event pairs, thereby completing the positioning of the microseismic event group.

[0081] This embodiment obtains a microseismic event set, pairs the microseismic events in the microseismic event set, obtains microseismic event pairs, and then clusters the microseismic event pairs to obtain microseismic event classes. When the number of microseismic event classes is not a preset value, the source spacing of each pair of microseismic event pairs in the microseismic event class is calculated, and the microseismic event group is located according to the source spacing. The arrangement of the sensor array is eliminated, and the accuracy of microseismic event group location is improved.

[0082] Based on the above first embodiment, a second embodiment of the present application is proposed. In this embodiment, the same or similar contents as those of the above first embodiment can be referred to the above introduction, and will not be described in detail later. Please refer to Figure 2 , the step S30 includes steps S301 to S305:

[0083] Step S301: When the number of the microseismic event class is not a preset value, performing secondary pairing on the microseismic events in the microseismic event class to obtain a target microseismic event pair;

[0084] In this embodiment, the preset value is set to 1. When the number of microseismic event classes is not 1, it is necessary to perform secondary pairing on the microseismic events in the microseismic event class. First, the microseismic event class containing three microseismic events with known actual positions is calculated, and all microseismic events in the microseismic event class are secondary paired. For other microseismic event classes, it is necessary to wait for the positions of all microseismic events in the first microseismic event class to be calculated, and select three microseismic events from them to join, and then start secondary pairing, and so on, to complete secondary pairing for all microseismic event classes.

[0085] It is understandable that when the number of microseismic event classes is 1, no secondary pairing is required and the subsequent steps can be performed directly.

[0086] Step S302: Calculate the amplitude mean and mean square angular frequency of the target microseismic signal of the target microseismic event pair;

[0087] It should be noted that microseismic signals refer to elastic waves or stress waves released by microseismic events.

[0088] The specific calculation process is as follows: first, divide the calculation window, which can be referred to Figure 3 , starting from 70 milliseconds after the P wave arrival, the microseismic signal is evenly divided into four small calculation windows of 50 milliseconds in length. The P wave arrival time refers to the arrival time of the first wave of the microseismic signal. According to the P wave arrival time, the original waveform of the target microseismic event 100 milliseconds before the P wave of the two signals arrives is intercepted, and the amplitude mean is calculated. The calculation result of this embodiment is 1.86×10 -7 m / s.

[0089] Calculate the mean square angular frequency of the waveform, specifically:

[0090]

[0091] in, is the time derivative of the waveform a(t), is the mean square angular frequency.

[0092] Step S303: obtaining a target amplitude segment according to the amplitude average value;

[0093] In this embodiment, the target amplitude segment refers to the amplitude segment with an amplitude value greater than the mean amplitude. According to the maximum value of the waveform envelope of the microseismic signal (the waveform envelope refers to the line between the highest point and the lowest point of the waveform, which is used to represent the fluctuation range of the waveform) and the calculation window, the waveform amplitude value within the calculation window is less than the mean amplitude of 1.86×10 -7 The amplitude segments of m / s are deleted, and the waveform amplitude values ​​in the calculation window are retained to be greater than or equal to the amplitude mean 1.86×10 -7 The amplitude segment of m / s is the target amplitude segment.

[0094] Step S304: calculating the maximum correlation of the microseismic event pairs in the microseismic event class according to the target amplitude segment;

[0095] According to the target amplitude segment, the maximum correlation within each small calculation window for each microseismic event pair is calculated, specifically:

[0096]

[0097] Among them, a time calculation window is set from time t1 to time t2, ∈ represents the time stretch factor, R max (t1, t2) represents the maximum correlation between the two signal waveforms from time t1 to time t2.

[0098] Step S305: Calculate the source spacing of the microseismic event pair according to the mean square angular frequency and the maximum correlation.

[0099] According to the mean square angular frequency and the maximum correlation, the source spacing is calculated as follows:

[0100]

[0101] Where r is the source spacing, is the mean square angular frequency, and v is the P-wave velocity.

[0102] According to the source spacing calculated in each small calculation window for each microseismic event pair, the calculated source spacing results of the four small windows are 3.4470m, 3.5461m, 3.9632m and 3.3571m respectively.

[0103] This embodiment performs secondary pairing of microseismic events in the microseismic event class to obtain a target microseismic event pair when the number of the microseismic event class is not a preset value; calculates the amplitude mean and mean square angular frequency of the target microseismic signal of the target microseismic event pair; obtains a target amplitude segment according to the amplitude mean; calculates the maximum correlation of the microseismic event pair in the microseismic event class according to the target amplitude segment; and calculates the source spacing of the microseismic event pair according to the mean square angular frequency and the maximum correlation. This improves the accuracy of positioning a microseismic event group.

[0104] Based on the above first embodiment, a third embodiment of the present application is proposed. In this embodiment, the same or similar contents as those of the above first embodiment can be referred to the above introduction, and will not be described in detail later. Please refer to Figure 4 , the step S40 includes steps S401 to S404:

[0105] Step S401: Calculate the source spacing mean and source spacing standard deviation according to the source spacing;

[0106] In this embodiment, if the calculated results of the source spacing are 3.4470 m, 3.5461 m, 3.9632 m, and 3.3571 m, respectively, the mean of the source spacing is 3.5784 m, and the standard deviation of the source spacing is 0.4641 m.

[0107] Step S402: Calculating the relative position of the microseismic event group according to the mean value of the source spacing and the standard deviation of the source spacing;

[0108] In this embodiment, according to the mean value and standard deviation of the source spacing of each pair of microseismic events, the relative position of the microseismic event group is obtained by solving the following equation using the conjugate gradient method:

[0109]

[0110] Among them, μ(δ t ) represents the mean of the source spacing; σ(δt) represents the standard deviation of the source spacing. cwi and δ t They represent the source spacing calculated by the method and the real source spacing, respectively, c is a constant; N e is the number of microseismic events; e i =(x i ,y i , z i ) represents the coordinate position of microseismic event i; using δ t，ij represents the source distance between event i and event j; parameters a1 = 0.4661, a2 = 48.9697, a3 = 2.4693, a4 = 4.2467, a5 = 1,1619, b1 = 0.1441, b2 = 101.0376, b3 = 120.3864, b4 = 2.8430, b5 = 6.0823; the relative position of the microseismic event group is finally obtained, refer to Figure 5 The projections of the relative positions of the microseismic event group in three different directions in the three-dimensional coordinate system are shown. The CWI positioning in the figure represents the relative position.

[0111] Step S403: obtaining the actual positions, translation positions, circumcenter coordinates and rotation angles of a preset number of microseismic events;

[0112] In this embodiment, the preset number is 3, that is, the actual positions, translation positions, circumcenter coordinates, and rotation angles of three microseismic events are obtained;

[0113] It should be noted that the microseismic event group includes three microseismic events with known actual positions. According to the actual position coordinates, the three actual position coordinates are connected to form an absolute triangle, and the circumcenter coordinates of the absolute triangle are calculated. Similarly, according to the relative position coordinates, the three relative position coordinates are connected to form a relative triangle, and the circumcenter coordinates of the relative triangle are calculated. In this embodiment, the circumcenter of the absolute triangle is (10336.69, 0.98, -27.54), and the circumcenter of the relative triangle is (10341.89, 4.89, -31.85). Therefore, the overall positioning result of the microseismic event group needs to be translated by 5.20m, 3.91m and -4.31m along the x, y and z axes respectively; the positioning result after translation is referenced Figure 6 .

[0114] Step S404: positioning the microseismic event group according to the relative position, the actual position, the translation position, the circumcenter coordinates and the rotation angle.

[0115] According to the position of the microseismic event group after the translation, the microseismic event group is spatially rotated as a whole, and the rotation angle is solved by using the Lagrange multiplier method to solve the minimum point of the following equation:

[0116]

[0117] Where (x i ,y i , z i ) are the actual positions of the three known microseismic events, (x p i ,y p i , z p i ) is the relative position of the three known actual positions of the microseismic events after translation, (x0, y0, z0) is the circumcenter coordinate of the absolute triangle, θ x ,θ y ,θ z are the angles of rotation around the x, y and z axes respectively. In this embodiment, the microseismic event group as a whole needs to be rotated by -112°, 77° and -149° along the x, y and z axes respectively with the absolute circumcenter as the rotation center, and finally the absolute position result of the microseismic event group is obtained. Figure 7 ;

[0118] This embodiment calculates the mean value and the standard deviation of the source spacing according to the source spacing; calculates the relative position of the microseismic event group according to the mean value and the standard deviation of the source spacing; obtains the actual position, translation position, circumcenter coordinates and rotation angle of a preset number of microseismic events; and locates the microseismic event group according to the relative position, the actual position, the translation position, the circumcenter coordinates and the rotation angle. The process of locating the microseismic event group in this embodiment is not affected by the sensor array, thereby improving the accuracy of locating the microseismic event group.

[0119] Based on the above first embodiment, a fourth embodiment of the present application is proposed. In this embodiment, the same or similar contents as those of the above first embodiment can be referred to the above introduction, and will not be described in detail later. Please refer to Figure 8 , the step S10 includes steps S101 to S104:

[0120] Step S101: Acquire a microseismic event set, and obtain microseismic events and target sensors that record all microseismic events according to the microseismic event set;

[0121] It should be noted that in actual operation, there may be a situation where no target sensor has recorded all microseismic events. In this case, it is necessary to reselect the event set until at least one target sensor that has recorded all microseismic events is found.

[0122] Step S102: obtaining the initial wave arrival time of the original microseismic signal according to the target sensor;

[0123] It should be noted that the initial wave arrival time, that is, the time when the P wave (i.e. the first arrival wave in the seismic wave, also called compression wave or longitudinal wave) arrives at a certain observation point or sensor in a microseismic event.

[0124] The target sensor records the P-wave arrival time of all original microseismic signals, from which the first wave arrival time of the original microseismic signals can be obtained.

[0125] Step S103: aligning the original microseismic signals according to the initial wave arrival time to obtain a target microseismic signal;

[0126] It should be noted that signal alignment means moving all signals to the same time when their P waves arrive. Fig. 9 , showing the signal alignment process of two signals.

[0127] Step S104: Pairing the microseismic events according to the target microseismic signal to obtain a microseismic event pair.

[0128] In this embodiment, the microseismic events are paired in pairs, referring to Fig.10 , the total number of microseismic events input in this embodiment is 13, then event 1 is respectively e Pairing is performed to form a total of 12 event pairs, event 2 is paired with event 3, event 4, ..., event N respectively. e Pairing is performed to form a total of 11 event pairs, and so on, the last event pair consists of event N e -1 and event N e A total of 78 microseismic event pairs were formed.

[0129] This embodiment obtains a microseismic event set, obtains microseismic events and a target sensor that records all microseismic events according to the microseismic event set; obtains the initial wave arrival time of the original microseismic signal according to the target sensor; aligns the original microseismic signals according to the initial wave arrival time to obtain a target microseismic signal; and pairs the microseismic events according to the target microseismic signal to obtain a microseismic event pair. This embodiment obtains a target microseismic signal and a microseismic event pair through microseismic events, and the process is not affected by the sensor array, which facilitates the subsequent microseismic event group positioning and improves the accuracy of microseismic event group positioning.

[0130] Based on the above first embodiment, a fifth embodiment of the present application is proposed. In this embodiment, the same or similar contents as those of the above first embodiment can be referred to the above introduction, and will not be described in detail later. Please refer to Fig.11 , the step S20 includes steps S201 to S202:

[0131] Step S201: Calculate and obtain the correlation of the microseismic event pairs;

[0132] It should be noted that the correlation represents the degree of correlation between pairs of microseismic events. The higher the correlation, the closer the relationship between the pairs of microseismic events.

[0133] In a feasible manner, the step of calculating the correlation of the microseismic event pairs includes:

[0134] Obtaining a target microseismic signal according to the microseismic event pair;

[0135] Dividing the target microseismic signal into a plurality of calculation windows according to a preset division rule;

[0136] The correlation between the microseismic event pairs is calculated through the multiple calculation windows.

[0137] Specifically, a correlation calculation window is divided according to the microseismic event pair and the target microseismic signal of the microseismic event pair. The correlation calculation window is a window starting from the arrival time of the aligned signal P wave and lasting for 1000 milliseconds.

[0138] Furthermore, according to the calculation window, the correlation of each pair of microseismic events recorded by each sensor is calculated, specifically:

[0139]

[0140] Among them, R is the correlation of the microseismic event pair, and s1 and s2 are two microseismic events respectively.

[0141] For each event pair, take the average value of the calculation results of multiple sensors as the correlation, which can be referred to Fig.12 , Fig.12 This is the result of using different sensors to calculate the correlation of each event pair in this embodiment, as well as the average correlation result.

[0142] Step S202: clustering the microseismic events according to the magnitude of the correlation to obtain microseismic event classes.

[0143] First, the microseismic event pairs are sorted according to the size of the correlation, and the microseismic event pairs with a correlation less than 0.5 are deleted, and the microseismic event pairs with an average correlation greater than or equal to 0.5 are retained. The retained microseismic event pairs are clustered into microseismic events. The clustering principle is: if event a is related to event b, and event b is also related to event c, then events a, b, and c are related, and events a, b, and c are classified as the same type of events. Through the secondary clustering principle, all microseismic event pairs are clustered to obtain microseismic event classes.

[0144] This embodiment screens the microseismic event pairs by calculating the correlation of the microseismic event pairs, removes some microseismic event pairs with little correlation, and clusters the remaining microseismic event pairs to obtain microseismic event classes. This method can effectively improve the accuracy of subsequent microseismic event group positioning.

[0145] For example, to help understand the implementation process of the microseismic event group location method obtained by combining this embodiment with the above-mentioned embodiment 1, please refer to Fig.13 , Fig.13 A brief flow chart of a method for locating a microseismic event group is provided, specifically: a microseismic event set is input, and according to the microseismic event set, the sensors and the number of sensors that record all microseismic events are counted; when the number of sensors is less than 1, the microseismic event set is reselected; when the number of sensors is greater than or equal to 1, according to the sensors, the P-wave arrival times of all original microseismic signals recorded by the sensors are obtained; according to the P-wave arrival times, all microseismic signals are aligned; according to the aligned signals, the microseismic events are paired and the correlation is calculated; according to the correlation, the microseismic events are clustered; according to the microseismic event classes, the number of microseismic event classes is counted; when the number of microseismic event classes is equal to 1, the P-wave arrival times of each pair are calculated. Source spacing of microseismic event pairs; calculating the relative position of a microseismic event group based on the source spacing of the microseismic event pairs; calculating the absolute position of the microseismic event group based on the relative position of the microseismic event group; when the number of microseismic event classes is not equal to 1, selecting a microseismic event class, performing secondary pairing on the microseismic events therein, and then performing source spacing of microseismic event pairs, calculating the relative position of the microseismic event group based on the source spacing of the microseismic event pairs, calculating the absolute position of the microseismic event group based on the relative position of the microseismic event group, and then selecting three microseismic events from them to add to the next microseismic event class for secondary pairing and subsequent steps, until the absolute positions of the microseismic events in all microseismic event classes are calculated.

[0146] This can effectively improve the positioning accuracy of microseismic event groups in the case of abnormal sensor array layout.

[0147] It should be noted that the above examples are only used to understand the present application and do not constitute a limitation on the microseismic event group locating method of the present application. More simple transformations based on this technical concept are all within the protection scope of the present application.

[0148] This application also provides a microseismic event group positioning device, please refer to Fig.14 , the microseismic event group positioning device comprises:

[0149] A pairing module 10 is used to obtain a microseismic event set and pair the microseismic events in the microseismic event set to obtain microseismic event pairs;

[0150] A clustering module 20, configured to cluster the microseismic event pairs to obtain microseismic event classes;

[0151] A calculation module 30, configured to calculate the source spacing of each pair of microseismic events in the microseismic event class when the number of the microseismic event class is not a preset value;

[0152] The positioning module 40 is used to locate the microseismic event group according to the source spacing.

[0153] In one embodiment, the calculation module 30 is further used to perform secondary pairing of the microseismic events in the microseismic event class to obtain a target microseismic event pair when the number of the microseismic event class is not a preset value; calculate the amplitude mean and the mean square angular frequency of the target microseismic signal of the target microseismic event pair; obtain a target amplitude segment based on the amplitude mean; calculate the maximum correlation of the microseismic event pair in the microseismic event class based on the target amplitude segment; and calculate the source spacing of the microseismic event pair based on the mean square angular frequency and the maximum correlation.

[0154] In one embodiment, the positioning module 40 is also used to calculate the source spacing mean and the source spacing standard deviation based on the source spacing; calculate the relative position of the microseismic event group based on the source spacing mean and the source spacing standard deviation; obtain the actual position, translation position, circumcenter coordinates and rotation angle of a preset number of microseismic events; and locate the microseismic event group based on the relative position, the actual position, the translation position, the circumcenter coordinates and the rotation angle.

[0155] In one embodiment, the pairing module 10 is also used to obtain a microseismic event set, obtain microseismic events and a target sensor that records all microseismic events based on the microseismic event set; obtain the initial wave arrival time of the original microseismic signal based on the target sensor; align the original microseismic signals based on the initial wave arrival time to obtain a target microseismic signal; and pair the microseismic events based on the target microseismic signal to obtain a microseismic event pair.

[0156] In one embodiment, the clustering module 20 is further configured to calculate the correlation between the microseismic event pairs; and cluster the microseismic events according to the magnitude of the correlation to obtain microseismic event classes.

[0157] In one embodiment, the clustering module 20 is further used to obtain a target microseismic signal according to the microseismic event pair; divide the target microseismic signal into multiple calculation windows according to a preset division rule; and calculate the correlation of the microseismic event pair through the multiple calculation windows.

[0158] In one embodiment, the microseismic event group locating device also includes an adding module for selecting a preset number of microseismic events from the microseismic event group and adding them to the microseismic event class, and returning to the step of calculating the source spacing of each pair of microseismic events in the microseismic event class when the number of microseismic events in the microseismic event class is not a preset value.

[0159] This embodiment obtains a microseismic event set, pairs the microseismic events in the microseismic event set, obtains microseismic event pairs, and then clusters the microseismic event pairs to obtain microseismic event classes. When the number of microseismic event classes is not a preset value, the source spacing of each pair of microseismic event pairs in the microseismic event class is calculated, and the microseismic event group is located according to the source spacing. The arrangement of the sensor array is eliminated, and the accuracy of microseismic event group location is improved.

[0160] The microseismic event group positioning device provided by the present application adopts the microseismic event group positioning method in the above-mentioned embodiment, which can solve the technical problem of how to improve the accuracy of microseismic event positioning. Compared with the prior art, the beneficial effects of the microseismic event group positioning device provided by the present application are the same as the beneficial effects of the microseismic event group positioning method provided by the above-mentioned embodiment, and other technical features in the microseismic event group positioning device are the same as the features disclosed in the above-mentioned embodiment method, which will not be repeated here.

[0161] The present application provides a microseismic event group locating device, which includes: at least one processor; and a memory communicatively connected to the at least one processor; wherein the memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor so that the at least one processor can execute the microseismic event group locating method in the above-mentioned embodiment one.

[0162] Reference below Fig.15, which shows a schematic diagram of the structure of a microseismic event group positioning device suitable for implementing the embodiment of the present application. The microseismic event group positioning device in the embodiment of the present application may include but is not limited to mobile terminals such as mobile phones, laptop computers, digital broadcast receivers, PDAs (Personal Digital Assistants), PADs (Portable Application Descriptions), PMPs (Portable Media Players), vehicle-mounted terminals (such as vehicle-mounted navigation terminals), etc., and fixed terminals such as digital TVs, desktop computers, etc. Fig.15 The microseismic event group positioning device shown is only an example and should not bring any limitation to the functions and scope of use of the embodiments of the present application.

[0163] like Fig.15 As shown, the microseismic event group location device may include a processing device 1001 (such as a central processing unit, a graphics processor, etc.), which can perform various appropriate actions and processes according to a program stored in a read-only memory (ROM) 1002 or a program loaded from a storage device 1003 to a random access memory (RAM) 1004. Various programs and data required for the operation of the microseismic event group location device are also stored in the RAM 1004. The processing device 1001, the ROM 1002, and the RAM 1004 are connected to each other via a bus 1005. An input / output (I / O) interface 1006 is also connected to the bus. Typically, the following systems can be connected to the I / O interface 1006: input devices 1007 including, for example, a touch screen, a touchpad, a keyboard, a mouse, an image sensor, a microphone, an accelerometer, a gyroscope, etc.; output devices 1008 including, for example, a liquid crystal display (LCD), a speaker, a vibrator, etc.; storage devices 1003 including, for example, a magnetic tape, a hard disk, etc.; and communication devices 1009. The communication device 1009 can allow the microseismic event group location device to communicate with other devices wirelessly or by wire to exchange data. Although the microseismic event group location device with various systems is shown in the figure, it should be understood that it is not required to implement or have all the systems shown. More or fewer systems may be implemented or have instead.

[0164] In particular, according to the embodiments disclosed in the present application, the process described above with reference to the flowchart can be implemented as a computer software program. For example, the embodiments disclosed in the present application include a computer program product, which includes a computer program carried on a computer-readable medium, and the computer program includes a program code for executing the method shown in the flowchart. In such an embodiment, the computer program can be downloaded and installed from a network through a communication device, or installed from a storage device 1003, or installed from a ROM 1002. When the computer program is executed by the processing device 1001, the above-mentioned functions defined in the method of the embodiment disclosed in the present application are executed.

[0165] The microseismic event group positioning device provided by the present application adopts the microseismic event group positioning method in the above embodiment, which can solve the technical problem of how to improve the accuracy of microseismic event positioning. Compared with the prior art, the beneficial effects of the microseismic event group positioning device provided by the present application are the same as the beneficial effects of the microseismic event group positioning method provided by the above embodiment, and other technical features in the microseismic event group positioning device are the same as the features disclosed in the method of the previous embodiment, which will not be repeated here.

[0166] It should be understood that the various parts disclosed in this application can be implemented by hardware, software, firmware or a combination thereof. In the description of the above embodiments, specific features, structures, materials or characteristics can be combined in any one or more embodiments or examples in a suitable manner.

[0167] The above is only a specific implementation of the present application, but the protection scope of the present application is not limited thereto. Any person skilled in the art who is familiar with the present technical field can easily think of changes or substitutions within the technical scope disclosed in the present application, which should be included in the protection scope of the present application. Therefore, the protection scope of the present application should be based on the protection scope of the claims.

[0168] The present application provides a computer-readable storage medium having computer-readable program instructions (ie, computer programs) stored thereon, and the computer-readable program instructions are used to execute the microseismic event group location method in the above-mentioned embodiment.

[0169] The computer-readable storage medium provided in the present application may be, for example, a USB flash drive, but is not limited to electrical, magnetic, optical, electromagnetic, infrared, or semiconductor systems, systems or devices, or any combination of the above. More specific examples of computer-readable storage media may include, but are not limited to: an electrical connection with one or more wires, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the above. In this embodiment, the computer-readable storage medium may be any tangible medium containing or storing a program that can be used by or in combination with an instruction execution system, system or device. The program code contained on the computer-readable storage medium may be transmitted using any appropriate medium, including but not limited to: wires, optical cables, RF (Radio Frequency), etc., or any suitable combination of the above.

[0170] The computer-readable storage medium may be included in the microseismic event cluster locating device; or may exist independently without being assembled into the microseismic event cluster locating device.

[0171] The computer-readable storage medium carries one or more programs. When the one or more programs are executed by the microseismic event group locating device, the microseismic event group locating device: obtains a microseismic event set, and pairs the microseismic events in the microseismic event set to obtain microseismic event pairs;

[0172] Clustering the microseismic event pairs to obtain microseismic event classes;

[0173] When the number of the microseismic event class is not a preset value, calculating the source spacing of each pair of microseismic events in the microseismic event class;

[0174] According to the source spacing, the microseismic event group is located.

[0175] Computer program code for performing the operations of the present application may be written in one or more programming languages ​​or a combination thereof, including object-oriented programming languages ​​such as Java, Smalltalk, C++, and conventional procedural programming languages ​​such as "C" or similar programming languages. The program code may be executed entirely on the user's computer, partially on the user's computer, as a separate software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or server. In the case of a remote computer, the remote computer may be connected to the user's computer through any type of network, including a local area network (LAN) or a wide area network (WAN), or may be connected to an external computer (e.g., via the Internet using an Internet service provider).

[0176] The flow chart and block diagram in the accompanying drawings illustrate the possible architecture, function and operation of the system, method and computer program product according to various embodiments of the present application. In this regard, each square box in the flow chart or block diagram can represent a module, a program segment or a part of a code, and the module, the program segment or a part of the code contains one or more executable instructions for realizing the specified logical function. It should also be noted that in some alternative implementations, the functions marked in the square box can also occur in a sequence different from that marked in the accompanying drawings. For example, two square boxes represented in succession can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each square box in the block diagram and / or flow chart, and the combination of the square boxes in the block diagram and / or flow chart can be implemented with a dedicated hardware-based system that performs a specified function or operation, or can be implemented with a combination of dedicated hardware and computer instructions.

[0177] The modules involved in the embodiments described in this application may be implemented by software or hardware, wherein the name of the module does not constitute a limitation on the unit itself in some cases.

[0178] The readable storage medium provided in the present application is a computer-readable storage medium, which stores computer-readable program instructions (i.e., computer programs) for executing the above-mentioned microseismic event group locating method, and can solve the technical problem of how to improve the accuracy of microseismic event group locating. Compared with the prior art, the beneficial effects of the computer-readable storage medium provided in the present application are the same as the beneficial effects of the microseismic event group locating method provided in the above-mentioned embodiment, and will not be repeated here.

[0179] The above descriptions are only some embodiments of the present application, and are not intended to limit the patent scope of the present application. All equivalent structural changes made using the contents of the present application specification and drawings under the technical concept of the present application, or direct / indirect applications in other related technical fields are included in the patent protection scope of the present application.

Claims

1. A method for locating a microseismic event group, characterized in that: The microseismic event group positioning method comprises: Acquire a microseismic event set, and pair the microseismic events in the microseismic event set to obtain microseismic event pairs; Clustering the microseismic event pairs to obtain microseismic event classes; When the number of the microseismic event class is not a preset value, calculating the source spacing of each pair of microseismic events in the microseismic event class; According to the source spacing, the microseismic event group is located.

2. The method according to claim 1, characterized in that When the number of the microseismic event class is not a preset value, the step of calculating the source spacing of each pair of microseismic event pairs in the microseismic event class comprises: When the number of the microseismic event class is not a preset value, performing secondary pairing on the microseismic events in the microseismic event class to obtain a target microseismic event pair; Calculating the amplitude mean and mean square angular frequency of the target microseismic signal of the target microseismic event pair; Obtaining a target amplitude segment according to the amplitude average value; calculating the maximum correlation of microseismic event pairs in the microseismic event class according to the target amplitude segment; The source spacing of the microseismic event pair is calculated based on the mean square angular frequency and the maximum correlation.

3. The method according to claim 1, characterized in that The step of locating a microseismic event group according to the source spacing includes: According to the source spacing, a source spacing mean and a source spacing standard deviation are calculated; Calculating the relative position of the microseismic event group according to the mean value of the source spacing and the standard deviation of the source spacing; Obtaining the actual positions, translational positions, circumcenter coordinates, and rotation angles of a preset number of microseismic events; The microseismic event group is located according to the relative position, the actual position, the translation position, the circumcenter coordinates and the rotation angle.

4. The method according to claim 1, characterized in that The step of obtaining a microseismic event set and pairing the microseismic events in the microseismic event set to obtain microseismic event pairs comprises: Acquire a microseismic event set, and obtain microseismic events and target sensors that record all microseismic events according to the microseismic event set; According to the target sensor, the initial wave arrival time of the original microseismic signal is obtained; According to the initial wave arrival time, the original microseismic signals are aligned to obtain a target microseismic signal; The microseismic events are paired according to the target microseismic signals to obtain microseismic event pairs.

5. The method according to claim 1, characterized in that The step of clustering the microseismic event pairs to obtain microseismic event classes comprises: Calculating the correlation of the microseismic event pairs; According to the magnitude of the correlation, the microseismic events are clustered to obtain microseismic event classes.

6. The method according to claim 5, characterized in that The step of calculating the correlation of the microseismic event pairs comprises: Obtaining a target microseismic signal according to the microseismic event pair; Dividing the target microseismic signal into a plurality of calculation windows according to a preset division rule; The correlation between the microseismic event pairs is calculated through the multiple calculation windows.

7. The method according to claim 1, characterized in that After the step of locating a microseismic event group according to the source spacing, the method further includes: A preset number of microseismic events are selected from the microseismic event group and added to the microseismic event class, and the step of calculating the source spacing of each pair of microseismic events in the microseismic event class when the number of microseismic events in the microseismic event class is not a preset value is returned.

8. A microseismic event group positioning device, characterized in that: The device comprises: A pairing module, used to obtain a microseismic event set, and pair the microseismic events in the microseismic event set to obtain microseismic event pairs; A clustering module, used for clustering the microseismic event pairs to obtain microseismic event classes; A calculation module, configured to calculate the source spacing of each pair of microseismic event pairs in the microseismic event class when the number of the microseismic event class is not a preset value; A positioning module is used to locate the microseismic event group according to the source spacing.

9. A microseismic event group positioning device, characterized in that: The device comprises: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the computer program is configured to implement the steps of the microseismic event group location method according to any one of claims 1 to 7.

10. A storage medium, characterized in that: The storage medium is a computer-readable storage medium, and a computer program is stored on the storage medium. When the computer program is executed by a processor, the steps of the microseismic event group locating method according to any one of claims 1 to 7 are implemented.

Images