Method, device, equipment and storage medium for picking up arrival time of S-wave of microseismic signal

Abstract

The present invention relates to the technical field of tunnel engineering, and discloses a method, device, equipment and storage medium for picking up the arrival time of the S wave of microseismic signals. The method includes: acquiring microseismic signals under a single microseismic event, and obtaining an S-wave arrival time solution signal and an auxiliary microseismic signal according to the microseismic signals; determining a target phase difference according to the S-wave arrival time solution signal and the auxiliary microseismic signal; determining an instantaneous phase difference according to the reconstructed signals of the S-wave arrival time solution signal and the auxiliary microseismic signal; picking up the arrival time of the S wave of the S-wave arrival time solution signal according to the target phase difference and the instantaneous phase difference; By the above method, the S-wave arrival time solution signal and the auxiliary microseismic signal are selected from the microseismic signals, then the target phase difference and the instantaneous phase difference are determined respectively, and then the arrival time of the S wave of the S-wave arrival time solution signal is determined according to the target phase difference and the instantaneous phase difference, so as to effectively improve the accuracy and applicability of picking up the S-wave arrival time.

Description

Technical Field

[0001] The present invention relates to the technical field of tunnel engineering, and particularly to a method, device, equipment and storage medium for picking up the arrival time of S-wave of microseismic signals. Background Art

[0002] Microseismic monitoring is a necessary safety monitoring means for high-stress tunnel engineering. In tunnel engineering, microseismic monitoring is usually carried out in units of the heading face. The monitoring array generally consists of 6-8 sensors and moves closely behind the heading face in an arrangement of 2-3 rows. Engineering practice shows that the spatial form of the monitoring array in tunnel engineering is approximately planar and the sensors in the same row are relatively close to each other. Only using the arrival time of P-wave of microseismic signals to solve the location of microseismic events of rock fracture, the error of the solved result is large and unstable. In addition, the microseismic monitoring in tunnel engineering is at the scale of hundreds of meters, and the distance between the sensor and the seismic source generally does not exceed 150m. When the S-wave arrives, the P-wave has not ended yet, and the initial arrival of the S-wave is severely and complexly intertwined with the tail of the P-wave. Eventually, the above method for picking up the arrival time of the S-wave has poor applicability in tunnel engineering.

[0003] The above content is only used to assist in understanding the technical solution of the present invention, and does not represent an admission that the above content is prior art. Summary of the Invention

[0004] The main object of the present invention is to provide a method, device, equipment and storage medium for picking up the arrival time of S-wave of microseismic signals, aiming to solve the technical problem of low accuracy and applicability of picking up the arrival time of S-wave in the prior art.

[0005] To achieve the above object, the present invention provides a method for picking up the arrival time of S-wave of microseismic signals, and the method for picking up the arrival time of S-wave of microseismic signals includes the following steps:

[0006] Obtain the microseismic signal under a single microseismic event, and obtain the S-wave arrival time solving signal and the auxiliary microseismic signal according to the microseismic signal;

[0007] Determine the target phase difference according to the S-wave arrival time solving signal and the auxiliary microseismic signal;

[0008] Determine the instantaneous phase difference according to the reconstructed signals of the S-wave arrival time solving signal and the auxiliary microseismic signal;

[0009] Pick up the arrival time of the S-wave of the S-wave arrival time solving signal according to the target phase difference and the instantaneous phase difference.

[0010] Optionally, the obtaining the microseismic signal under a single microseismic event, and obtaining the S-wave arrival time solving signal and the auxiliary microseismic signal according to the microseismic signal includes:

[0011] Obtain the microseismic signal under a single microseismic event, and determine the excitation time of the microseismic signal;

[0012] Sort the microseismic signals according to the time according to a preset sorting rule;

[0013] Name the sorted microseismic signals in sequence according to a preset naming rule;

[0014] Select the S-wave arrival time solution signal and the auxiliary microseismic signal from the named microseismic signals.

[0015] Optionally, the determining the target phase difference according to the S-wave arrival time solution signal and the auxiliary microseismic signal includes:

[0016] Perform a Fourier transform on the S-wave arrival time solution signal, and pick up the main frequency of the signal according to the transformation result;

[0017] Calculate the wavelength of the S-wave according to the main frequency of the signal;

[0018] Obtain the propagation speed of the P-wave in the rock medium;

[0019] Obtain the P-wave arrival time of the S-wave arrival time solution signal and the P-wave arrival time of the auxiliary microseismic signal;

[0020] Calculate the signal propagation distance difference according to the propagation speed of the P-wave in the rock medium, the P-wave arrival time of the S-wave arrival time solution signal, and the P-wave arrival time of the auxiliary microseismic signal;

[0021] Calculate the target phase difference according to the wavelength of the S-wave and the signal propagation distance difference.

[0022] Optionally, the determining the instantaneous phase difference according to the S-wave arrival time solution signal and the reconstructed signal of the auxiliary microseismic signal includes:

[0023] Perform a global empirical mode decomposition on the S-wave arrival time solution signal and the auxiliary microseismic signal respectively to obtain a certain number of empirical modes;

[0024] Reconstruct the signals of the certain number of empirical modes by setting the frequency range respectively to obtain the reconstructed signals of the S-wave arrival time solution signal and the auxiliary microseismic signal;

[0025] Perform a Hilbert transform on the reconstructed signals of the S-wave arrival time solution signal and the auxiliary microseismic signal respectively to obtain the instantaneous phase of the reconstructed signal of the S-wave arrival time solution signal and the instantaneous phase of the reconstructed signal of the auxiliary microseismic signal;

[0026] Perform a difference calculation on the instantaneous phase of the reconstructed signal of the S-wave arrival time solution signal and the instantaneous phase of the reconstructed signal of the auxiliary microseismic signal to obtain the instantaneous phase difference.

[0027] Optionally, picking up the arrival time of the S wave of the signal for solving the S wave arrival time based on the target phase difference and the instantaneous phase difference includes:

[0028] Calculating a target parameter based on the target phase difference and the instantaneous phase difference;

[0029] Calculating the period of the S wave based on the main frequency of the signal;

[0030] Calculating the intensity of the instantaneous phase difference based on the target parameter and the period of the S wave;

[0031] Picking up the arrival time of the S wave of the signal for solving the S wave arrival time according to the intensity of the instantaneous phase difference through a target decision rule.

[0032] Optionally, picking up the arrival time of the S wave of the signal for solving the S wave arrival time according to the intensity of the instantaneous phase difference through a target decision rule includes:

[0033] Obtaining the end time of the microseismic waveform;

[0034] Generating a target time range based on the end time of the microseismic waveform and the arrival time of the P wave of the signal for solving the S wave arrival time;

[0035] Picking up the first significant peak that appears within the target time range according to the intensity of the instantaneous phase difference through a target decision rule;

[0036] Finding the corresponding time according to the first significant peak, and taking the time as the arrival time of the S wave of the signal for solving the S wave arrival time.

[0037] Optionally, after picking up the arrival time of the S wave of the signal for solving the S wave arrival time based on the target phase difference and the instantaneous phase difference, it further includes:

[0038] If there are other signals that need to pick up the arrival time of the S wave, then encapsulating each step of picking up the arrival time of the S wave of the signal for solving the S wave arrival time above to obtain a signal S wave arrival time picking strategy;

[0039] Picking up the arrival time of the S wave of the signal according to the signal S wave arrival time picking strategy.

[0040] In addition, to achieve the above object, the present invention also proposes a microseismic signal S wave arrival time picking device, and the microseismic signal S wave arrival time picking device includes:

[0041] An acquisition module, configured to acquire a microseismic signal under a single microseismic event, and obtain a signal for solving the S wave arrival time and an auxiliary microseismic signal according to the microseismic signal;

[0042] A determination module, configured to determine a target phase difference according to the S-wave arrival time solution signal and the auxiliary microseismic signal;

[0043] The determination module is further configured to determine an instantaneous phase difference according to the reconstructed signals of the S-wave arrival time solution signal and the auxiliary microseismic signal;

[0044] A picking module, configured to pick up the arrival time of the S-wave of the S-wave arrival time solution signal according to the target phase difference and the instantaneous phase difference.

[0045] In addition, to achieve the above object, the present invention further provides a microseismic signal S-wave arrival time picking device, where the microseismic signal S-wave arrival time picking device includes: a memory, a processor, and a microseismic signal S-wave arrival time picking program stored on the memory and executable on the processor, and the microseismic signal S-wave arrival time picking program is configured to implement the microseismic signal S-wave arrival time picking method as described above.

[0046] In addition, to achieve the above object, the present invention further provides a storage medium, on which a microseismic signal S-wave arrival time picking program is stored, and when the microseismic signal S-wave arrival time picking program is executed by a processor, it implements the microseismic signal S-wave arrival time picking method as described above.

[0047] The microseismic signal S-wave arrival time picking method provided by the present invention obtains a microseismic signal under a single microseismic event, and obtains an S-wave arrival time solution signal and an auxiliary microseismic signal according to the microseismic signal; determines a target phase difference according to the S-wave arrival time solution signal and the auxiliary microseismic signal; determines an instantaneous phase difference according to the reconstructed signals of the S-wave arrival time solution signal and the auxiliary microseismic signal; picks up the arrival time of the S-wave of the S-wave arrival time solution signal according to the target phase difference and the instantaneous phase difference; through the above method, the S-wave arrival time solution signal and the auxiliary microseismic signal are selected from the microseismic signal, then the target phase difference and the instantaneous phase difference are respectively determined, and then the arrival time of the S-wave of the S-wave arrival time solution signal is determined according to the target phase difference and the instantaneous phase difference, so as to effectively improve the accuracy and applicability of picking up the S-wave arrival time. Description of the Drawings

[0048] Figure 1 It is a schematic structural diagram of a microseismic signal S-wave arrival time picking device in the hardware operating environment related to the embodiment solution of the present invention;

[0049] Figure 2 It is a schematic flowchart of the first embodiment of the microseismic signal S-wave arrival time picking method of the present invention;

[0050] Figure 3 It is a wave train diagram of microseismic signals received by each sensor in an embodiment of the microseismic signal S-wave arrival time picking method of the present invention;

[0051] Figure 4 It is the spectrogram after Fourier transform of an embodiment of the method for picking up the arrival time of S-wave of microseismic signals in the present invention;

[0052] Figure 5 It is the schematic diagram of empirical mode of an embodiment of the method for picking up the arrival time of S-wave of microseismic signals in the present invention;

[0053] Figure 6 It is the schematic diagram of instantaneous phase of an embodiment of the method for picking up the arrival time of S-wave of microseismic signals in the present invention;

[0054] Figure 7 It is the schematic diagram of instantaneous phase difference of an embodiment of the method for picking up the arrival time of S-wave of microseismic signals in the present invention;

[0055] Figure 8 It is the schematic flow chart of the second embodiment of the method for picking up the arrival time of S-wave of microseismic signals in the present invention;

[0056] Figure 9 It is the curve graph of target parameters of an embodiment of the method for picking up the arrival time of S-wave of microseismic signals in the present invention;

[0057] Figure 10 It is the curve graph of the intensity of instantaneous phase difference of an embodiment of the method for picking up the arrival time of S-wave of microseismic signals in the present invention;

[0058] Figure 11 It is the schematic overall flow chart of an embodiment of the method for picking up the arrival time of S-wave of microseismic signals in the present invention;

[0059] Figure 12 It is the schematic diagram of functional modules of the first embodiment of the device for picking up the arrival time of S-wave of microseismic signals in the present invention.

[0060] The realization, functional features and advantages of the object of the present invention will be further described with reference to the embodiments and the accompanying drawings. Specific Embodiments

[0061] It should be understood that the specific embodiments described herein are only used to explain the present invention and are not used to limit the present invention.

[0062] Refer to Figure 1 , Figure 1 It is the schematic diagram of the structure of the device for picking up the arrival time of S-wave of microseismic signals in the hardware operating environment related to the solution of the embodiment of the present invention.

[0063] Such as Figure 1As shown in the figure, the microseismic signal S-wave arrival time picking device may include: a processor 1001, such as a Central Processing Unit (CPU), a communication bus 1002, a user interface 1003, a network interface 1004, and a memory 1005. Among them, the communication bus 1002 is used to realize the connection and communication between these components. The user interface 1003 may include a display screen and an input unit such as a keyboard. Optionally, the user interface 1003 may further include a standard wired interface and a wireless interface. The network interface 1004 may optionally include a standard wired interface and a wireless interface (such as a Wireless-Fidelity (Wi-Fi) interface). The memory 1005 may be a high-speed Random Access Memory (RAM) or a stable Non-Volatile Memory (NVM), such as a disk memory. Optionally, the memory 1005 may also be a storage device independent of the aforementioned processor 1001.

[0064] Those skilled in the art can understand that Figure 1 the structure shown in does not constitute a limitation on the microseismic signal S-wave arrival time picking device, and it may include more or fewer components than shown, or combine some components, or have different component arrangements.

[0065] As Figure 1 shown, the memory 1005, as a storage medium, may include an operating system, a network communication module, a user interface module, and a microseismic signal S-wave arrival time picking program.

[0066] In Figure 1 the microseismic signal S-wave arrival time picking device shown, the network interface 1004 is mainly used for data communication with a network integrated platform workstation; the user interface 1003 is mainly used for data interaction with users; the processor 1001 and the memory 1005 in the microseismic signal S-wave arrival time picking device of the present invention may be provided in the microseismic signal S-wave arrival time picking device. The microseismic signal S-wave arrival time picking device calls the microseismic signal S-wave arrival time picking program stored in the memory 1005 through the processor 1001 and executes the microseismic signal S-wave arrival time picking method provided by the embodiments of the present invention.

[0067] Based on the above hardware structure, an embodiment of the microseismic signal S-wave arrival time picking method of the present invention is proposed.

[0068] Referring to Figure 2 , Figure 2 is a schematic flowchart of the first embodiment of the microseismic signal S-wave arrival time picking method of the present invention.

[0069] In the first embodiment, the method for picking up the arrival time of the S-wave of the microseismic signal includes the following steps:

[0070] Step S10: Obtain the microseismic signal under a single microseismic event, and obtain the S-wave arrival time solving signal and the auxiliary microseismic signal according to the microseismic signal.

[0071] It should be noted that the execution subject of this embodiment is a device for picking up the arrival time of the S-wave of the microseismic signal, and it can also be other devices that can achieve the same or similar functions, such as an S-wave arrival time picking controller, etc. This embodiment does not limit this. In this embodiment, the S-wave arrival time picking controller is taken as an example for illustration.

[0072] It should be understood that the microseismic signal refers to multiple microseismic signals under a single microseismic event, that is, the microseismic signal can be a set of microseismic event signals. The S-wave arrival time solving signal refers to the signal for which the arrival time of the S-wave needs to be solved, and the auxiliary microseismic signal refers to the microseismic signal used to assist in calculating the arrival time of the S-wave.

[0073] Further, step S10 includes: obtaining the microseismic signal under a single microseismic event, and determining the time when the microseismic signal is excited; sorting the microseismic signals according to the time according to a preset sorting rule; naming the sorted microseismic signals in sequence according to a preset naming rule; and selecting the S-wave arrival time solving signal and the auxiliary microseismic signal from the named microseismic signals.

[0074] It can be understood that when the microseismic signal under a single microseismic event is obtained, the time when each channel of the microseismic signal is excited is recorded at the same time. The preset sorting rule refers to the rule for sorting multiple microseismic signals, and this preset sorting rule can be a time sequence sorting rule. The preset naming rule refers to the rule for naming the sorted microseismic signals, that is, the named microseismic signals can be W i (i = 1, 2,..., n). Specifically, the first excited microseismic signal is named W 1 ; the second excited microseismic signal is named W 2 ; and so on. The last excited signal is named W n . Then, the S-wave arrival time solving signal and the auxiliary microseismic signal are selected from the named microseismic signals. For example, the S-wave arrival time solving signal is W i , and the auxiliary microseismic signal is W i-1 Reference Figure 3 , Figure 3 is the wave train diagram of the microseismic signal received by each sensor, specifically including the microseismic signals W 1 to W 6 . There are a total of 6 microseismic signals, and the excitation times in sequence are the 54th, 56th, 59th, 71st, 73rd, and 76th sampling points. The above microseismic signals are in the order of being excited from first to last, atFigure 3 They are arranged in sequence from top to bottom.

[0075] Step S20: Determine the target phase difference according to the S-wave arrival time solving signal and the auxiliary microseismic signal.

[0076] Further, step S20 includes: performing Fourier transform on the S-wave arrival time solving signal, and picking up the main frequency of the signal according to the transform result; calculating the wavelength of the S-wave according to the main frequency of the signal; obtaining the propagation speed of the P-wave in the rock medium; obtaining the P-wave arrival time of the S-wave arrival time solving signal and the P-wave arrival time of the auxiliary microseismic signal; calculating the signal propagation distance difference according to the propagation speed of the P-wave in the rock medium, the P-wave arrival time of the S-wave arrival time solving signal, and the P-wave arrival time of the auxiliary microseismic signal; calculating the target phase difference according to the wavelength of the S-wave and the signal propagation distance difference.

[0077] It should be understood that the main frequency of the signal refers to the main frequency of the S-wave arrival time solving signal, and the signal propagation distance difference refers to the distance difference between the propagation of the S-wave arrival time solving signal and the auxiliary microseismic signal. After obtaining the propagation speed of the P-wave in the rock medium, the P-wave arrival time of the S-wave arrival time solving signal, and the P-wave arrival time of the auxiliary microseismic signal, calculate the signal propagation distance difference, specifically:

[0078]

[0079] Among them, PDD represents the signal propagation distance difference. represents the P-wave arrival time of the S-wave arrival time solving signal. represents the P-wave arrival time of the auxiliary microseismic signal, V p represents the propagation speed of the P-wave in the rock medium.

[0080] It can be understood that after calculating the signal propagation distance difference, it is judged whether the signal propagation distance difference is less than the wavelength of the S-wave. If so, the target phase difference is calculated in combination with the wavelength of the S-wave, specifically:

[0081] α e = PDD / l i · 2π;

[0082] Among them, α e represents the target phase difference, PDD represents the signal propagation distance difference, l i represents the wavelength of the S-wave.

[0083] It should be understood that referring to Figure 4 , Figure 4 is the spectrogram after Fourier transform. Specifically, the abscissa corresponding to the maximum amplitude point is 302. At this time, the abscissa 302 is the S-wave arrival time solving signal W 3The main frequency is 302 Hz, and then combined with the wave velocity of the S wave, which is 3000 m / s, the arrival time of the S wave is calculated to solve for signal W. 3 The wavelength l of the S wave 3 is 9.93 m, and the S wave period T 3 is 3.31 ms. The auxiliary microseismic signal W is selected 2 , and based on the arrival time of the P wave and the arrival time of the S wave of the auxiliary microseismic signal W 2 to solve for signal W. 3 The arrival time of the P wave of signal W is calculated, and the signal propagation distance difference PDD between the two signals is 2.5 m. Then, when the signal propagation distance difference PDD is less than the wavelength l of the S wave 3 , combined with the wavelength l of the S wave 3 the target phase difference α is calculated e to be π / 2.

[0084] Step S30: Determine the instantaneous phase difference according to the reconstruction signal of the S wave arrival time solution signal and the auxiliary microseismic signal.

[0085] It should be understood that the instantaneous phase difference refers to the difference between the instantaneous phase of the reconstruction signal of the S wave arrival time solution signal and the instantaneous phase of the reconstruction signal of the auxiliary microseismic signal.

[0086] Further, step S30 includes: respectively performing global empirical mode decomposition on the S wave arrival time solution signal and the auxiliary microseismic signal to obtain a certain number of empirical modes; respectively performing signal reconstruction on the certain number of empirical modes by setting frequency ranges to obtain the reconstruction signals of the S wave arrival time solution signal and the auxiliary microseismic signal; respectively performing Hilbert transform on the reconstruction signals of the S wave arrival time solution signal and the auxiliary microseismic signal to obtain the instantaneous phase of the reconstruction signal of the S wave arrival time solution signal and the instantaneous phase of the reconstruction signal of the auxiliary microseismic signal; calculating the difference between the instantaneous phase of the reconstruction signal of the S wave arrival time solution signal and the instantaneous phase of the reconstruction signal of the auxiliary microseismic signal to obtain the instantaneous phase difference.

[0087] It can be understood that after obtaining the S wave arrival time solution signal and the auxiliary microseismic signal, global empirical mode decomposition is respectively performed on the S wave arrival time solution signal and the auxiliary microseismic signal to obtain a certain number of empirical modes, then noise filtering is performed on the certain number of empirical modes by setting frequency ranges, and then the certain number of empirical modes within the set frequency range are added together for signal reconstruction. Then, Hilbert transform is respectively performed on the S wave arrival time solution signal and the auxiliary microseismic signal according to the reconstruction signal, and then the difference is calculated between the instantaneous phase of the reconstruction signal of the S wave arrival time solution signal and the instantaneous phase of the reconstruction signal of the auxiliary microseismic signal to obtain the instantaneous phase difference. Refer to Figure 5 , Figure 5It is a schematic diagram of empirical mode. Specifically: select the empirical modes with the main frequency between 100 Hz and 350 Hz. The empirical modes within the frequency range are IMF3 to IMF6. Add IMF3 to IMF6 to obtain the signal W for solving the arrival time of the S wave. 3 The reconstructed signal W' 3 According to the process of obtaining the reconstructed signal W' 3 Perform the same decomposition and reconstruction on the auxiliary microseismic signal W 2 Finally, obtain the reconstructed signal W' 2 of the auxiliary microseismic signal W 2 Then, according to the reconstructed signal W' 2 and the reconstructed signal W' 3 Perform Hilbert transform on the signal for solving the arrival time of the S wave and the auxiliary microseismic signal respectively to obtain the instantaneous phase of the reconstructed signal. Refer to Figure 6 , Figure 6 It is a schematic diagram of the instantaneous phase. Specifically: the instantaneous phase continuously jumps between -π and π with time. Then, perform a difference between the instantaneous phase of the reconstructed signal W' 2 of the signal for solving the arrival time of the S wave and the instantaneous phase of the reconstructed signal W' 3 of the auxiliary microseismic signal to obtain the instantaneous phase difference Δα(t). Refer to Figure 7 , Figure 7 It is a schematic diagram of the instantaneous phase difference. Specifically: the instantaneous phase difference continuously jumps between -2π and 2π with time.

[0088] Step S40: Pick up the arrival time of the S wave of the signal for solving the arrival time of the S wave according to the target phase difference and the instantaneous phase difference.

[0089] It can be understood that after obtaining the target phase difference and the instantaneous phase difference, pick up the arrival time of the S wave of the signal for solving the arrival time of the S wave according to the target determination rule based on the target phase difference and the instantaneous phase difference.

[0090] Further, after step S40, it further includes: if there are other signals for which the arrival time of the S wave needs to be picked up, encapsulate the above steps of picking up the arrival time of the S wave of the signal for solving the arrival time of the S wave to obtain the signal S wave arrival time picking strategy; pick up the arrival time of the S wave of the signal according to the signal S wave arrival time picking strategy.

[0091] It should be understood that the S-wave arrival time picking strategy of the signal refers to the strategy of picking the S-wave arrival time of the microseismic signal in the tunnel project based on the instantaneous phase difference. After picking the S-wave arrival time and solving for the arrival time of the S-wave of the signal, it is continued to determine whether there are other signals for which the S-wave arrival time needs to be picked. If so, the various steps of picking the S-wave arrival time and solving for the arrival time of the S-wave of the signal are encapsulated as the S-wave arrival time picking strategy of the signal, and then the S-wave arrival time picking strategy of the signal is used to pick the S-wave arrival time of the signal until the S-wave arrival times of all signals are picked.

[0092] In this embodiment, by obtaining the microseismic signal under a single microseismic event, the S-wave arrival time solving signal and the auxiliary microseismic signal are obtained according to the microseismic signal; the target phase difference is determined according to the S-wave arrival time solving signal and the auxiliary microseismic signal; the instantaneous phase difference is determined according to the reconstructed signals of the S-wave arrival time solving signal and the auxiliary microseismic signal; the S-wave arrival time of the S-wave arrival time solving signal is picked according to the target phase difference and the instantaneous phase difference; in the above manner, the S-wave arrival time solving signal and the auxiliary microseismic signal are selected from the microseismic signal, then the target phase difference and the instantaneous phase difference are respectively determined, and then the S-wave arrival time of the S-wave arrival time solving signal is determined according to the target phase difference and the instantaneous phase difference, so as to effectively improve the accuracy and applicability of picking the S-wave arrival time.

[0093] In one embodiment, as Figure 8 described, based on the first embodiment, the second embodiment of the microseismic signal S-wave arrival time picking method of the present invention is proposed, and the step S40 includes:

[0094] Step S401, calculate the target parameter according to the target phase difference and the instantaneous phase difference.

[0095] It should be understood that after obtaining the target phase difference and the instantaneous phase difference, the target parameter is calculated using the target phase difference and the instantaneous phase difference, specifically:

[0096] i(t) = (|Δα(t) - α e |) 2 + (2π - |α e - Δα(t)|) 2 ;

[0097] wherein, i(t) represents the target parameter, α e represents the target phase difference, and Δα(t) represents the instantaneous phase difference.

[0098] It can be understood that referring to Figure 9 Figure 9 is the curve graph of the target parameter, specifically: as time goes by, the value of the target parameter is constantly changing, that is, between time 50 - 250, the value of the target parameter fluctuates between 15 - 55.​

[0099] Step S402: Calculate the period of the S-wave based on the main frequency of the signal.

[0100] Step S403: Calculate the instantaneous phase difference intensity according to the target parameter and the period of the S-wave.

[0101] It should be understood that after obtaining the target parameter and the period of the S-wave, the instantaneous phase difference intensity is calculated using the target parameter and the period of the S-wave. Specifically:

[0102]

[0103] where IPDI represents the instantaneous phase difference intensity, T i represents the period of the S-wave, and i(t) represents the target parameter.

[0104] Step S404: Pick up the arrival time of the S-wave of the signal by solving the arrival time of the S-wave according to the instantaneous phase difference intensity through the target determination rule.

[0105] It can be understood that the target determination rule refers to the rule for determining the first significant peak that appears within the target time range, that is, using the target determination rule to pick up the arrival time of the S-wave according to the instantaneous phase difference intensity to solve the arrival time of the S-wave of the signal.

[0106] Further, step S404 includes: obtaining the end time of the microseismic waveform; generating a target time range by solving the arrival time of the P-wave of the signal based on the end time of the microseismic waveform and the arrival time of the S-wave; picking up the first significant peak that appears within the target time range according to the instantaneous phase difference intensity through the target determination rule; finding the corresponding time according to the first significant peak, and taking the time as the arrival time of the S-wave of the signal for solving the arrival time of the S-wave.

[0107] It should be understood that after obtaining the end time of the microseismic waveform, a target time range is generated in combination with the arrival time of the P-wave of the signal for solving the arrival time of the S-wave. For example, if the arrival time of the P-wave of the signal for solving the arrival time of the S-wave is the 59th sampling point and the end time of the microseismic waveform is the 200th sampling point, then the 59th sampling point to the 200th sampling point is defined as the target time range, and then the first significant peak that appears within the target time range is picked up according to the instantaneous phase difference intensity through the target determination rule, and then the time corresponding to the first significant peak is used as the arrival time of the S-wave of the signal for solving the arrival time of the S-wave. Refer to Figure 10 , Figure 10It is an instantaneous phase difference intensity curve graph. Specifically: after generating the target time range U, pick up the first significant peak that appears within the target time range according to the instantaneous phase difference intensity, and then determine that the time corresponding to this first significant peak is 98, that is, time 98 is the arrival time of the S wave of the S wave arrival time solving signal, so that the arrival time of the S wave of the S wave arrival time solving signal can be picked up more conveniently, and when the tail of the P wave intersects with the head of the S wave and the arrival time of the S wave takeoff point is picked up, it provides more accurate data information for the calculation of microseismic positioning.

[0108] It can be understood that referring to Figure 11 , Figure 11 is a schematic diagram of the overall process. Specifically: name the microseismic signals received by each sensor in the order of excitation time as W i (i = 1, 2,..., n). When i = 2, perform FFT transformation on the S wave arrival time solving signal, then calculate the signal main frequency of the S wave arrival time solving signal, and then calculate the period and wavelength of the S wave according to the signal main frequency. Then calculate the signal propagation distance difference PDD between the S wave arrival time solving signal and the auxiliary microseismic signal. When it is determined that the signal propagation distance difference PDD is less than the wavelength of the S wave, calculate the target phase difference in combination with the wavelength of the S wave, and then decompose and reconstruct the S wave arrival time solving signal and the auxiliary microseismic signal based on the ensemble empirical mode decomposition (EEMD) to obtain the reconstructed signal W' 3 of the S wave arrival time solving signal and the reconstructed signal W' 2 of the auxiliary microseismic signal. Then perform Hilbert transform on the reconstructed signal W' 3 of the S wave arrival time solving signal and the reconstructed signal W 2 ' of the auxiliary microseismic signal. Then subtract the instantaneous phase of the reconstructed signal W' 2 of the S wave arrival time solving signal from the instantaneous phase of the reconstructed signal W' 3 of the auxiliary microseismic signal to obtain the instantaneous phase difference Δα(t), and then calculate the instantaneous phase difference intensity IPID. Then generate the target time range according to the microseismic waveform end time and the P wave arrival time of the S wave arrival time solving signal, and determine the arrival time of the S wave of the S wave arrival time solving signal. Then continue to judge whether there are other signals that need to pick up the S wave arrival time. If so, i = i + 1, and use the S wave arrival time picking strategy of the signal to pick up the arrival time of the S wave of the signal until the arrival times of the S waves of all signals are picked up.

[0109] This embodiment calculates a target parameter based on the target phase difference and the instantaneous phase difference; calculates the period of the S-wave according to the main frequency of the signal; calculates the instantaneous phase difference intensity according to the target parameter and the period of the S-wave; and picks up the arrival time of the S-wave of the signal by using a target determination rule based on the instantaneous phase difference intensity. In this way, the target parameter is calculated by using the target phase difference and the instantaneous phase difference, then the instantaneous phase difference intensity is calculated in combination with the period of the S-wave, and then the arrival time of the S-wave of the signal is picked up by using the target determination rule based on the instantaneous phase difference intensity, so as to effectively improve the accuracy of picking up the arrival time of the S-wave, and further effectively improve the accuracy of microseismic positioning.

[0110] In addition, an embodiment of the present invention further provides a storage medium, on which a program for picking up the arrival time of the S-wave of a microseismic signal is stored. When the program for picking up the arrival time of the S-wave of the microseismic signal is executed by a processor, the steps of the method for picking up the arrival time of the S-wave of the microseismic signal as described above are implemented.

[0111] Since this storage medium adopts all the technical solutions of the above-mentioned all embodiments, it has at least all the beneficial effects brought by the technical solutions of the above-mentioned embodiments, which will not be elaborated herein one by one.

[0112] In addition, referring to Figure 12 , an embodiment of the present invention further provides a device for picking up the arrival time of the S-wave of a microseismic signal. The device for picking up the arrival time of the S-wave of the microseismic signal includes:

[0113] An acquisition module 10, configured to acquire a microseismic signal under a single microseismic event, and obtain an S-wave arrival time solution signal and an auxiliary microseismic signal according to the microseismic signal.

[0114] A determination module 20, configured to determine a target phase difference according to the S-wave arrival time solution signal and the auxiliary microseismic signal.

[0115] The determination module 20 is further configured to determine an instantaneous phase difference according to the reconstructed signals of the S-wave arrival time solution signal and the auxiliary microseismic signal.

[0116] A picking-up module 30, configured to pick up the arrival time of the S-wave of the S-wave arrival time solution signal according to the target phase difference and the instantaneous phase difference.

[0117] In this embodiment, by acquiring the microseismic signal under a single microseismic event, the S-wave arrival time solving signal and the auxiliary microseismic signal are obtained based on the microseismic signal; the target phase difference is determined according to the S-wave arrival time solving signal and the auxiliary microseismic signal; the instantaneous phase difference is determined according to the reconstructed signals of the S-wave arrival time solving signal and the auxiliary microseismic signal; the arrival time of the S-wave of the S-wave arrival time solving signal is picked up according to the target phase difference and the instantaneous phase difference; in the above manner, the S-wave arrival time solving signal and the auxiliary microseismic signal are selected from the microseismic signal, then the target phase difference and the instantaneous phase difference are respectively determined, and then the arrival time of the S-wave of the S-wave arrival time solving signal is determined according to the target phase difference and the instantaneous phase difference, so as to effectively improve the accuracy and applicability of picking up the S-wave arrival time.

[0118] It should be noted that the above-described workflow is only illustrative and does not limit the protection scope of the present invention. In actual applications, those skilled in the art can select some or all of them according to actual needs to achieve the purpose of the solution of this embodiment, and no limitation is made here.

[0119] In addition, for the technical details not described in detail in this embodiment, reference can be made to the method for picking up the S-wave arrival time of the microseismic signal provided in any embodiment of the present invention, which will not be elaborated here.

[0120] In one embodiment, the acquisition module 10 is further configured to acquire the microseismic signal under a single microseismic event and determine the time when the microseismic signal is excited; sort the microseismic signals according to the time according to a preset sorting rule; name the sorted microseismic signals in sequence according to a preset naming rule; select the S-wave arrival time solving signal and the auxiliary microseismic signal from the named microseismic signals.

[0121] In one embodiment, the determination module 20 is further configured to perform Fourier transform on the S-wave arrival time solving signal, pick up the signal main frequency according to the transformation result; calculate the wavelength of the S-wave according to the signal main frequency; acquire the propagation speed of the P-wave in the rock medium; acquire the P-wave arrival time of the S-wave arrival time solving signal and the P-wave arrival time of the auxiliary microseismic signal; calculate the signal propagation distance difference according to the propagation speed of the P-wave in the rock medium, the P-wave arrival time of the S-wave arrival time solving signal, and the P-wave arrival time of the auxiliary microseismic signal; calculate the target phase difference according to the wavelength of the S-wave and the signal propagation distance difference.

[0122] In one embodiment, the determining module 20 is further configured to perform empirical mode decomposition on the S-wave arrival time solving signal and the auxiliary microseismic signal respectively to obtain a number of empirical modes; perform signal reconstruction on the number of empirical modes by setting frequency ranges respectively to obtain the reconstructed signals of the S-wave arrival time solving signal and the auxiliary microseismic signal; perform Hilbert transform on the reconstructed signals of the S-wave arrival time solving signal and the auxiliary microseismic signal respectively to obtain the instantaneous phase of the reconstructed signal of the S-wave arrival time solving signal and the instantaneous phase of the reconstructed signal of the auxiliary microseismic signal; perform a difference calculation on the instantaneous phase of the reconstructed signal of the S-wave arrival time solving signal and the instantaneous phase of the reconstructed signal of the auxiliary microseismic signal to obtain an instantaneous phase difference.

[0123] In one embodiment, the picking module 30 is further configured to calculate a target parameter according to the target phase difference and the instantaneous phase difference; calculate the period of the S-wave according to the signal main frequency; calculate the instantaneous phase difference intensity according to the target parameter and the period of the S-wave; pick the arrival time of the S-wave of the S-wave arrival time solving signal according to the target determination rule based on the instantaneous phase difference intensity.

[0124] In one embodiment, the picking module 30 is further configured to obtain the end time of the microseismic waveform; generate a target time range according to the end time of the microseismic waveform and the arrival time of the P-wave of the S-wave arrival time solving signal; pick the first significant peak that appears within the target time range according to the target determination rule based on the instantaneous phase difference intensity; find the corresponding time according to the first significant peak and use the time as the arrival time of the S-wave of the S-wave arrival time solving signal.

[0125] In one embodiment, the picking module 30 is further configured to, if there are other signals for which the arrival time of the S-wave needs to be picked, package the above steps for picking the arrival time of the S-wave of the S-wave arrival time solving signal to obtain a signal S-wave arrival time picking strategy; pick the arrival time of the S-wave of the signal according to the signal S-wave arrival time picking strategy.

[0126] Other embodiments or implementation methods of the microseismic signal S-wave arrival time picking device of the present invention can refer to the above method embodiments and will not be repeated here.

[0127] In addition, it should be noted that in this text, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, such that a process, method, article or system comprising a series of elements not only includes those elements but also includes other elements not expressly listed, or further includes elements inherent to such process, method, article or system. Without further limitation, an element defined by the statement "comprising one..." does not exclude the existence of additional identical elements in the process, method, article or system comprising such element.

[0128] The serial numbers of the above embodiments of the present invention are only for description and do not represent the superiority or inferiority of the embodiments.

[0129] Through the description of the above embodiments, those skilled in the art can clearly understand that the above embodiment methods can be implemented by means of software plus a necessary general hardware platform. Of course, they can also be implemented by hardware, but in many cases the former is a better implementation. Based on such an understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product is stored in a storage medium (such as Read Only Memory (ROM) / RAM, magnetic disk, optical disc), and includes several instructions to enable a terminal device (which can be a mobile phone, computer, integrated platform workstation, or network device, etc.) to execute the methods described in the various embodiments of the present invention.

[0130] The above are only the preferred embodiments of the present invention and do not limit the patent scope of the present invention accordingly. Any equivalent structural or equivalent process transformation made by using the description of the present invention and the drawings, or directly or indirectly applied in other related technical fields, shall be equally included in the patent protection scope of the present invention.

Claims

1. A method for picking up the arrival time of S - wave of micro - seismic signals, characterized in that, the method for picking up the arrival time of S - wave of micro - seismic signals includes the following steps: Obtain micro - seismic signals under a single micro - seismic event, and obtain S - wave arrival - time solving signals and auxiliary micro - seismic signals according to the micro - seismic signals; Determine the target phase difference according to the S - wave arrival - time solving signals and the auxiliary micro - seismic signals; Determine the instantaneous phase difference according to the reconstruction signals of the S - wave arrival - time solving signals and the auxiliary micro - seismic signals; Pick up the arrival time of the S - wave of the S - wave arrival - time solving signals according to the target phase difference and the instantaneous phase difference.

2. The method for picking up the arrival time of S - wave of micro - seismic signals according to claim 1, characterized in that, the step of obtaining micro - seismic signals under a single micro - seismic event and obtaining S - wave arrival - time solving signals and auxiliary micro - seismic signals according to the micro - seismic signals includes: Obtain micro - seismic signals under a single micro - seismic event, and determine the excitation time of the micro - seismic signals; Sort the micro - seismic signals according to the time according to a preset sorting rule; Name the sorted micro - seismic signals in sequence according to a preset naming rule; Select S - wave arrival - time solving signals and auxiliary micro - seismic signals from the named micro - seismic signals.

3. The method for picking up the arrival time of S - wave of micro - seismic signals according to claim 1, characterized in that, the step of determining the target phase difference according to the S - wave arrival - time solving signals and the auxiliary micro - seismic signals includes: Perform Fourier transform on the S - wave arrival - time solving signals, and pick up the signal main frequency according to the transform result; Calculate the wavelength of the S - wave according to the signal main frequency; Obtain the propagation speed of P - wave in rock media; Obtain the P - wave arrival time of the S - wave arrival - time solving signals and the P - wave arrival time of the auxiliary micro - seismic signals; Calculate the signal propagation distance difference according to the propagation speed of P - wave in rock media, the P - wave arrival time of the S - wave arrival - time solving signals and the P - wave arrival time of the auxiliary micro - seismic signals; Calculate the target phase difference according to the wavelength of the S - wave and the signal propagation distance difference.

4. The method for picking up the arrival time of S - wave of micro - seismic signals according to claim 1, characterized in that, the step of determining the instantaneous phase difference according to the reconstruction signals of the S - wave arrival - time solving signals and the auxiliary micro - seismic signals includes: Perform global empirical mode decomposition on the S - wave arrival - time solving signals and the auxiliary micro - seismic signals respectively to obtain a number of empirical modes; Reconstruct signals for the number of empirical modes respectively by setting frequency ranges to obtain the reconstruction signals of the S - wave arrival - time solving signals and the auxiliary micro - seismic signals; Perform Hilbert transform on the reconstruction signals of the S - wave arrival - time solving signals and the auxiliary micro - seismic signals respectively to obtain the instantaneous phase of the reconstruction signal of the S - wave arrival - time solving signals and the instantaneous phase of the reconstruction signal of the auxiliary micro - seismic signals; Perform difference calculation on the instantaneous phase of the reconstruction signal of the S - wave arrival - time solving signals and the instantaneous phase of the reconstruction signal of the auxiliary micro - seismic signals to obtain the instantaneous phase difference.

5. The method for picking up the arrival time of S - wave of micro - seismic signals according to claim 1, characterized in that, the step of picking up the arrival time of the S - wave of the S - wave arrival - time solving signals according to the target phase difference and the instantaneous phase difference includes: Calculate the target parameter based on the target phase difference and the instantaneous phase difference; Calculate the period of the S-wave according to the main frequency of the signal; Calculate the instantaneous phase difference intensity based on the target parameter and the period of the S-wave; Pick up the arrival time of the S-wave of the signal by using the target decision rule according to the instantaneous phase difference intensity to solve the arrival time of the S-wave of the signal.

6. The method for picking up the arrival time of the S-wave of the microseismic signal according to claim 5, characterized in that, The step of picking up the arrival time of the S-wave of the signal by using the target decision rule according to the instantaneous phase difference intensity to solve the arrival time of the S-wave of the signal includes: Obtain the end time of the microseismic waveform; Generate a target time range according to the end time of the microseismic waveform and the arrival time of the P-wave of the signal to solve the arrival time of the S-wave; Pick up the first significant peak that appears within the target time range by using the target decision rule according to the instantaneous phase difference intensity; Find the corresponding time according to the first significant peak, and use the time as the arrival time of the S-wave of the signal to solve the arrival time of the S-wave.

7. The method for picking up the arrival time of the S-wave of the microseismic signal according to any one of claims 1 to 6, characterized in that, After picking up the arrival time of the S-wave of the signal by using the target phase difference and the instantaneous phase difference to solve the arrival time of the S-wave of the signal, it further includes: If there are other signals that need to pick up the arrival time of the S-wave, then encapsulate the steps of picking up the arrival time of the S-wave of the signal to solve the arrival time of the S-wave of the signal, and obtain the signal S-wave arrival time picking strategy; Pick up the arrival time of the S-wave of the signal according to the signal S-wave arrival time picking strategy.

8. A device for picking up the arrival time of the S-wave of the microseismic signal, characterized in that, The device for picking up the arrival time of the S-wave of the microseismic signal includes: An acquisition module, configured to acquire a microseismic signal under a single microseismic event, and obtain a signal for solving the arrival time of the S-wave and an auxiliary microseismic signal according to the microseismic signal; A determination module, configured to determine a target phase difference according to the signal for solving the arrival time of the S-wave and the auxiliary microseismic signal; The determination module is further configured to determine an instantaneous phase difference according to the reconstructed signal of the signal for solving the arrival time of the S-wave and the auxiliary microseismic signal; A picking module, configured to pick up the arrival time of the S-wave of the signal for solving the arrival time of the S-wave according to the target phase difference and the instantaneous phase difference.

9. A device for picking up the arrival time of the S-wave of the microseismic signal, characterized in that, The device for picking up the arrival time of the S-wave of the microseismic signal includes: a memory, a processor, and a microseismic signal S-wave arrival time picking program stored on the memory and executable on the processor, and the microseismic signal S-wave arrival time picking program is configured to implement the microseismic signal S-wave arrival time picking method according to any one of claims 1 to 7.

10. A storage medium, characterized in that, A microseismic signal S-wave arrival time picking program is stored on the storage medium, and when the microseismic signal S-wave arrival time picking program is executed by a processor, it implements the microseismic signal S-wave arrival time picking method according to any one of claims 1 to 7.

Images