Microseismic signal P wave arrival time picking method, device, equipment and storage medium

Abstract

The present invention relates to the field of tunnel engineering technology, and discloses a method, device, equipment and storage medium for picking up the P-wave arrival time of a microseismic signal. The method comprises: acquiring a microseismic signal under a microseismic event, and obtaining a P-wave arrival time solution signal and an auxiliary microseismic signal according to the microseismic signal; determining an ideal phase difference according to the P-wave arrival time solution signal and the auxiliary microseismic signal; determining an actual phase difference according to the P-wave arrival time solution signal and the auxiliary microseismic signal; picking up the arrival time of the P wave of the P-wave arrival time solution signal according to the ideal phase difference and the actual phase difference; in the above manner, selecting the P-wave arrival time solution signal and the auxiliary microseismic signal from the microseismic signal under the microseismic event, and then determining the ideal phase difference and the actual phase difference respectively, and then picking up the arrival time of the P wave of the P-wave arrival time solution signal according to the ideal phase difference and the actual phase difference, thereby effectively improving the accuracy of picking up the P-wave arrival time.

Description

Technical Field

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

[0002] Typical deep geological hazards such as rockbursts and stress-induced landslides are characterized by high intensity, high frequency, and significant damage. Microseismic monitoring, currently one of the most direct and effective technologies for monitoring and warning of stress-induced disasters, has become an essential monitoring tool for the safety of deep tunnel projects. Accurately locating the microseismic source of rock fractures is a prerequisite for early warning of stress-induced disasters, and the accuracy of picking up the P-wave of the microseismic signal is the key factor affecting the error in microseismic source location.

[0003] Currently, the commonly used methods for picking up the P-wave arrival time of microseismic signals are the long-short window ratio method (STA / LTA) and the Akaike Information Criterion method (AIC). However, these methods are only applicable to scenarios with a high signal-to-noise ratio (SNR) of microseismic signals. However, in tunnel engineering, microseismic monitoring equipment is arranged near the tunnel face. The rock fracture microseismic signals are severely interfered with by various types of construction noise, and the SNR of microseismic signals is often low. If these methods are continued to be used to pick up the P-wave arrival time, the accuracy of the picked P-wave arrival time will inevitably be low.

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

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

[0006] To achieve the above object, the present invention provides a method for picking up the arrival time of a microseismic signal P wave, the method comprising the following steps:

[0007] Acquire a microseismic signal under a microseismic event, and obtain a P-wave arrival time solution signal and an auxiliary microseismic signal according to the microseismic signal;

[0008] Determining an ideal phase difference based on the P-wave arrival time solution signal and the auxiliary microseismic signal;

[0009] determining an actual phase difference based on the P-wave arrival time solution signal and the auxiliary microseismic signal;

[0010] The arrival time of the P wave of the P wave arrival time solution signal is picked up according to the ideal phase difference and the actual phase difference.

[0011] Optionally, acquiring a microseismic signal under a microseismic event and obtaining a P-wave arrival time solution signal and an auxiliary microseismic signal according to the microseismic signal includes:

[0012] Acquire a microseismic signal under a microseismic event, and determine a time when a sensor corresponding to the microseismic signal is excited;

[0013] sorting the microseismic signals according to the time in a preset time sequence;

[0014] Name the sorted microseismic signals in sequence according to the preset naming rules;

[0015] The P-wave arrival time solution signal and auxiliary microseismic signal are determined based on the named microseismic signal.

[0016] Optionally, determining the P-wave arrival time solution signal and the auxiliary microseismic signal according to the named microseismic signal includes:

[0017] Select the P-wave arrival time solution signal from the named microseismic signals;

[0018] Performing Fourier transform on the P-wave arrival time solution signal, and picking up the signal main frequency according to the transformation result;

[0019] Calculating the wavelength of the P wave based on the main frequency of the signal;

[0020] An auxiliary microseismic signal is selected from the named microseismic signals according to the wavelength of the P wave.

[0021] Optionally, the microseismic signal under the microseismic event includes a current microseismic signal and multiple other microseismic signals;

[0022] The selecting of the auxiliary microseismic signal from the named microseismic signal according to the wavelength of the P wave includes:

[0023] respectively calculating straight-line distances between the sensor corresponding to the current microseismic signal and the multiple sensors corresponding to the other multiple microseismic signals;

[0024] determining a minimum straight-line distance among a plurality of said straight-line distances;

[0025] When the minimum straight-line distance is smaller than the wavelength of the P wave, the microseismic signal corresponding to the minimum straight-line distance is used as an auxiliary microseismic signal.

[0026] Optionally, determining the ideal phase difference according to the P-wave arrival time solution signal and the auxiliary microseismic signal includes:

[0027] Calculating the straight-line distance between the sensor corresponding to the P-wave arrival time solution signal and the sensor corresponding to the auxiliary microseismic signal;

[0028] Obtain the propagation velocity of P waves in geological media;

[0029] The ideal phase difference is calculated based on the straight-line distance, propagation velocity, and signal main frequency between the sensor corresponding to the P-wave arrival time solution signal and the sensor corresponding to the auxiliary microseismic signal.

[0030] Optionally, determining the actual phase difference according to the P-wave arrival time solution signal and the auxiliary microseismic signal includes:

[0031] Performing Hilbert transform on the P-wave arrival time solution signal and the auxiliary microseismic signal respectively to obtain the instantaneous phase of the P-wave arrival time solution signal and the instantaneous phase of the auxiliary microseismic signal;

[0032] The instantaneous phase of the P-wave arrival signal and the instantaneous phase of the auxiliary microseismic signal are subtracted to obtain an actual phase difference.

[0033] Optionally, picking up the arrival time of the P wave of the P wave arrival time solution signal according to the ideal phase difference and the actual phase difference includes:

[0034] Calculate the period of the P wave based on the dominant frequency of the signal;

[0035] Calculating instantaneous phase difference intensity according to the period of the P wave, the ideal phase difference, and the actual phase difference;

[0036] The maximum instantaneous phase difference intensity is selected from the multiple instantaneous phase difference intensities, and the time corresponding to the maximum instantaneous phase difference intensity is used as the arrival time of the P wave of the P wave arrival time solution signal.

[0037] In addition, to achieve the above-mentioned purpose, the present invention further proposes a device for picking up the arrival time of a microseismic signal P wave, the device comprising:

[0038] An acquisition module is used to acquire a microseismic signal under a microseismic event, and obtain a P-wave arrival time solution signal and an auxiliary microseismic signal according to the microseismic signal;

[0039] a determination module, configured to determine an ideal phase difference based on the P-wave arrival time solution signal and the auxiliary microseismic signal;

[0040] The determination module is further configured to determine an actual phase difference based on the P-wave arrival time solution signal and the auxiliary microseismic signal;

[0041] A picking module is used to pick up the arrival time of the P wave of the P wave arrival solution signal according to the ideal phase difference and the actual phase difference.

[0042] In addition, to achieve the above-mentioned purpose, the present invention also proposes a microseismic signal P-wave arrival time picking device, which includes: a memory, a processor, and a microseismic signal P-wave arrival time picking program stored on the memory and runnable on the processor, and the microseismic signal P-wave arrival time picking program is configured to implement the microseismic signal P-wave arrival time picking method described above.

[0043] In addition, to achieve the above-mentioned purpose, the present invention also proposes a storage medium, on which a microseismic signal P wave arrival time picking program is stored. When the microseismic signal P wave arrival time picking program is executed by a processor, the microseismic signal P wave arrival time picking method described above is implemented.

[0044] The present invention proposes a method for picking up the P-wave arrival time of microseismic signals, which obtains the microseismic signals under the microseismic events, obtains the P-wave arrival time solution signal and the auxiliary microseismic signals according to the microseismic signals; determines the ideal phase difference according to the P-wave arrival time solution signal and the auxiliary microseismic signals; determines the actual phase difference according to the P-wave arrival time solution signal and the auxiliary microseismic signals; picks up the arrival time of the P-wave of the P-wave arrival time solution signal according to the ideal phase difference and the actual phase difference; through the above-mentioned method, the P-wave arrival time solution signal and the auxiliary microseismic signals are selected from the microseismic signals under the microseismic events, and then the ideal phase difference and the actual phase difference are determined respectively, and then the arrival time of the P-wave of the P-wave arrival time solution signal is picked up according to the ideal phase difference and the actual phase difference, thereby effectively improving the accuracy of picking up the P-wave arrival time. BRIEF DESCRIPTION OF THE DRAWINGS

[0045] Figure 1 1 is a schematic diagram of the structure of a microseismic signal P-wave arrival time picking device in the hardware operating environment involved in an embodiment of the present invention;

[0046] Figure 2 1. It is a flow chart of the first embodiment of the method for picking up the arrival time of microseismic signals P waves according to the present invention;

[0047] Figure 3 A schematic diagram of sensor distribution in accordance with an embodiment of a method for picking up the arrival time of microseismic signals P waves according to the present invention;

[0048] Figure 4 Schematic diagram of instantaneous phase of an embodiment of a method for picking up the arrival time of microseismic signal P waves according to the present invention;

[0049] Figure 5 Schematic diagram of actual phase difference of an embodiment of a method for picking up the arrival time of microseismic signals P waves according to the present invention;

[0050] Figure 6 A schematic diagram of the instantaneous phase difference intensity of an embodiment of a method for picking up the arrival time of microseismic signals P waves according to the present invention;

[0051] Figure 7 This is a flow chart of a second embodiment of a method for picking up the arrival time of microseismic signals P waves according to the present invention;

[0052] Figure 8 A microseismic signal wave train diagram according to an embodiment of a method for picking up the arrival time of a microseismic signal P wave according to the present invention;

[0053] Figure 9 This is an amplitude-frequency diagram of an embodiment of a method for picking up the arrival time of microseismic signals P waves according to the present invention;

[0054] Figure 10 This is a schematic diagram of the overall flow of an embodiment of a method for picking up the arrival time of microseismic signals P waves according to the present invention;

[0055] Figure 11 Schematic diagram of the functional modules of the first embodiment of the microseismic signal P-wave arrival time pickup device of the present invention.

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

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

[0058] Reference Figure 1 , Figure 1 This is a schematic diagram of the structure of a microseismic signal P-wave arrival time picking device in the hardware operating environment involved in an embodiment of the present invention.

[0059] like Figure 1 As shown, the microseismic signal P-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. The communication bus 1002 is used to enable communication between these components. The user interface 1003 may include a display and an input unit, such as a keyboard. Optionally, the user interface 1003 may also include a standard wired interface or a wireless interface. The network interface 1004 may optionally include a standard wired interface or 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 drive. Optionally, the memory 1005 may be a storage device independent of the processor 1001.

[0060] Those skilled in the art will understand that Figure 1 The structure shown in the figure does not constitute a limitation on the microseismic signal P-wave arrival time picking device, and may include more or fewer components than shown in the figure, or combine certain components, or arrange the components differently.

[0061] like Figure 1 As 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 P-wave arrival time picking program.

[0062] exist Figure 1 In the microseismic signal P-wave arrival time picking device shown, the network interface 1004 is mainly used for data communication with the network integration platform workstation; the user interface 1003 is mainly used for data interaction with the user; the processor 1001 and the memory 1005 in the microseismic signal P-wave arrival time picking device of the present invention can be set in the microseismic signal P-wave arrival time picking device, and the microseismic signal P-wave arrival time picking device calls the microseismic signal P-wave arrival time picking program stored in the memory 1005 through the processor 1001 and executes the microseismic signal P-wave arrival time picking method provided by the embodiment of the present invention.

[0063] Based on the above hardware structure, an embodiment of a method for picking the arrival time of a microseismic signal P wave of the present invention is proposed.

[0064] Reference Figure 2 , Figure 2 FIG. 1 is a flow chart of a first embodiment of a method for picking up the arrival time of microseismic signals P waves according to the present invention.

[0065] In a first embodiment, the method for picking the arrival time of a microseismic signal P wave comprises the following steps:

[0066] Step S10: Acquire a microseismic signal under a microseismic event, and obtain a P-wave arrival time solution signal and an auxiliary microseismic signal according to the microseismic signal.

[0067] It should be noted that the executor of this embodiment is a microseismic signal P-wave arrival time picking device, and may also be other devices that can achieve the same or similar functions, such as a P-wave arrival time picking controller, etc. This embodiment does not limit this. In this embodiment, the P-wave arrival time picking controller is used as an example for explanation.

[0068] 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 P-wave arrival time solution signal refers to the signal that needs to solve the P-wave arrival time, and the auxiliary microseismic signal refers to the microseismic signal used to assist in calculating the P-wave arrival time.

[0069] It should be emphasized that the reference Figure 3 , Figure 3This is a schematic diagram of sensor distribution, including section 1, section 2, excavated section, tunnel face and unexcavated section. The sensors are named , and the sensor Set around section 1, sensor Set around section 2.

[0070] Step S20: determining an ideal phase difference based on the P-wave arrival time solution signal and the auxiliary microseismic signal.

[0071] It can be understood that after the P-wave arrival time solution signal and the auxiliary microseismic signal are selected, the ideal phase difference is determined using the P-wave arrival time solution signal and the auxiliary microseismic signal.

[0072] Furthermore, step S20 includes: calculating the straight-line distance between the sensor corresponding to the P-wave arrival time solution signal and the sensor corresponding to the auxiliary microseismic signal; obtaining the propagation velocity of the P-wave in the geological medium; and calculating the ideal phase difference based on the straight-line distance, propagation velocity, and signal main frequency between the sensor corresponding to the P-wave arrival time solution signal and the sensor corresponding to the auxiliary microseismic signal.

[0073] It should be understood that the straight-line distance refers to the distance between the sensor corresponding to the P-wave arrival signal and the sensor corresponding to the auxiliary microseismic signal, and the signal dominant frequency refers to the dominant frequency of the P-wave arrival signal. After obtaining the straight-line distance, the ideal phase difference is calculated by combining the propagation speed of the P-wave in the geological medium and the signal dominant frequency. Specifically, it is:

[0074] ;

[0075] in, represents the ideal phase difference, D represents the straight-line distance between the sensor corresponding to the P-wave arrival signal and the sensor corresponding to the auxiliary microseismic signal, Indicates the signal main frequency, Indicates the propagation speed of P waves in geological media. For example, D is 31, is 128, is 6500, then is 220°.

[0076] Step S30: determining the actual phase difference according to the P-wave arrival time solution signal and the auxiliary microseismic signal.

[0077] It should be understood that the actual phase difference refers to the actual difference between the instantaneous phase of the P-wave arrival time solution signal and the instantaneous phase of the auxiliary microseismic signal.

[0078] Furthermore, step S30 includes: performing Hilbert transform on the P-wave arrival time solution signal and the auxiliary microseismic signal respectively to obtain the instantaneous phase of the P-wave arrival time solution signal and the instantaneous phase of the auxiliary microseismic signal; performing difference calculation on the instantaneous phase of the P-wave arrival time solution signal and the instantaneous phase of the auxiliary microseismic signal to obtain an actual phase difference.

[0079] It can be understood that after obtaining the P-wave arrival time solution signal and the auxiliary microseismic signal, the P-wave arrival time solution signal and the auxiliary microseismic signal are subjected to Hilbert transform respectively to obtain the instantaneous phase of the P-wave arrival time solution signal and the instantaneous phase of the auxiliary microseismic signal. Figure 4 , Figure 4 This is a schematic diagram of the instantaneous phase. Specifically, as time goes by, the instantaneous phase of the P-wave arrival signal and the instantaneous phase of the auxiliary microseismic signal vary between -π and π. Then, the actual phase difference is calculated based on the instantaneous phase of the P-wave arrival signal and the instantaneous phase of the auxiliary microseismic signal, which is specifically:

[0080] ;

[0081] in, Indicates the actual phase difference, Indicates the instantaneous phase of the signal when the P wave arrives, Represents the instantaneous phase of the auxiliary microseismic signal, refer to Figure 5 , Figure 5 Schematic diagram of the actual phase difference. Specifically, as time goes by, the actual phase difference changes between -2π and 2π.

[0082] Step S40: Pick up the P-wave arrival time of the P-wave solution signal according to the ideal phase difference and the actual phase difference.

[0083] It can be understood that after obtaining the ideal phase difference and the actual phase difference, the arrival time of the P wave of the signal is solved by picking up the P wave arrival time according to the ideal phase difference and the actual phase difference.

[0084] Furthermore, step S40 includes: calculating the period of the P wave based on the main frequency of the signal; calculating the instantaneous phase difference intensity based on the period of the P wave, the ideal phase difference and the actual phase difference; selecting the maximum instantaneous phase difference intensity from multiple instantaneous phase difference intensities, and using the time corresponding to the maximum instantaneous phase difference intensity as the arrival time of the P wave of the P wave arrival solution signal.

[0085] It should be understood that after obtaining the signal main frequency, the period of the P wave is calculated based on the signal main frequency, specifically:

[0086] ;

[0087] in, represents the period of the P wave, Indicates the main frequency of the signal. For example, if the main frequency of the signal is 128 Hz, the period of the P wave is 7.813 ms.

[0088] It can be understood that after obtaining the period of the P wave, the instantaneous phase difference intensity is calculated by combining the ideal phase difference and the actual phase difference, specifically:

[0089] ;

[0090] in, represents the instantaneous phase difference intensity, Indicates the actual phase difference, Indicates the ideal phase difference, reference Figure 6 , Figure 6 Schematic diagram of instantaneous phase difference intensity. Specifically, as time goes by, the instantaneous phase difference intensity is constantly changing, and between time 500-800, the instantaneous phase difference intensity changes between 800-1600.

[0091] It should be understood that after obtaining multiple instantaneous phase difference intensities, the maximum instantaneous phase difference intensity is selected from the multiple instantaneous phase difference intensities, and then the time corresponding to the maximum instantaneous phase difference intensity is used as the arrival time of the P wave of the P-wave arrival solution signal. For example, the time corresponding to the maximum instantaneous phase difference intensity is t3, and at this time t3 is used as the arrival time of the P wave of the P-wave arrival solution signal.

[0092] This embodiment obtains the microseismic signal under the microseismic event, obtains the P-wave arrival time solution signal and the auxiliary microseismic signal according to the microseismic signal; determines the ideal phase difference according to the P-wave arrival time solution signal and the auxiliary microseismic signal; determines the actual phase difference according to the P-wave arrival time solution signal and the auxiliary microseismic signal; picks up the arrival time of the P wave of the P-wave arrival time solution signal according to the ideal phase difference and the actual phase difference; through the above method, the P-wave arrival time solution signal and the auxiliary microseismic signal are selected from the microseismic signal under the microseismic event, and then the ideal phase difference and the actual phase difference are determined respectively, and then the arrival time of the P wave of the P-wave arrival time solution signal is picked up according to the ideal phase difference and the actual phase difference, thereby effectively improving the accuracy of picking up the P-wave arrival time.

[0093] In one embodiment, if Figure 7 Based on the first embodiment, a second embodiment of the method for picking the arrival time of microseismic signals P waves of the present invention is proposed. Step S10 includes:

[0094] Step S101: Acquire a microseismic signal under a microseismic event, and determine the time when a sensor corresponding to the microseismic signal is excited.

[0095] It should be understood that, when a microseismic signal is obtained under a microseismic event, the time when the sensor corresponding to the microseismic signal is excited is determined, referring to Figure 8 , Figure 8 It is a microseismic signal wave train diagram, which specifically includes four sensors, and the microseismic signals under microseismic events are obtained through the four sensors.

[0096] Step S102: sorting the microseismic signals according to the time in a preset time sequence.

[0097] It can be understood that the preset time sequence refers to the order in which the microseismic signals are sorted according to the time relationship. The preset time sequence can be a chronological order, that is, the microseismic signals are arranged from top to bottom according to the preset time sequence.

[0098] Step S103: naming the sorted microseismic signals in sequence according to a preset naming rule.

[0099] It should be understood that the preset naming rule refers to the rule for naming the sorted microseismic signals. For example, the four sensors are named 、 、 as well as , and the sensor 、 、 as well as The microseismic signal is Figure 8 shown.

[0100] Step S104: determining a P-wave arrival time solution signal and an auxiliary microseismic signal based on the named microseismic signal.

[0101] It is understandable that after obtaining the named microseismic signal, the P-wave arrival time solution signal and the auxiliary microseismic signal are determined according to the named microseismic signal.

[0102] Furthermore, step S104 includes: selecting a P-wave arrival time solution signal from the named microseismic signals; performing Fourier transform on the P-wave arrival time solution signal, and picking up the signal main frequency according to the transformation result; calculating the wavelength of the P wave according to the signal main frequency; and selecting an auxiliary microseismic signal from the named microseismic signals according to the wavelength of the P wave.

[0103] It should be understood that after obtaining the named microseismic signal, a P-wave arrival time solution signal is selected from the named microseismic signal. For example, the P-wave arrival time solution signal is Then, the P-wave arrival time solution signal is Fourier transformed to obtain the transformation result, and the amplitude-frequency diagram is generated according to the change result. Figure 9 , Figure 9The amplitude-frequency diagram is as follows: the horizontal axis represents the frequency, the vertical axis represents the amplitude, and the amplitude ranges from 0 to 4.5×10 -6 .

[0104] It can be understood that after obtaining the amplitude-frequency diagram, the horizontal coordinate corresponding to the maximum amplitude point in the amplitude-frequency diagram is picked up as the main frequency of the signal. For example, the main frequency of the signal is 128HZ, and then the wavelength of the P wave is calculated based on the main frequency of the signal, specifically:

[0105] ;

[0106] in, represents the wavelength of the P wave, represents the propagation speed of P waves in geological media, Indicates the main frequency of the signal.

[0107] Furthermore, the microseismic signal under the microseismic event includes the current microseismic signal and multiple other microseismic signals; the selection of the auxiliary microseismic signal from the named microseismic signals according to the wavelength of the P wave includes: respectively calculating the straight-line distances between the sensor corresponding to the current microseismic signal and multiple sensors corresponding to the multiple other microseismic signals; determining the minimum straight-line distance among the multiple straight-line distances; when the minimum straight-line distance is less than the wavelength of the P wave, using the microseismic signal corresponding to the minimum straight-line distance as the auxiliary microseismic signal.

[0108] It is understandable that after obtaining the microseismic signal under the microseismic event, since the microseismic signal under the microseismic event includes the current microseismic signal and multiple other microseismic signals, it is necessary to calculate the straight-line distance between the sensor corresponding to the current microseismic signal and the multiple sensors corresponding to the multiple other microseismic signals, for example, to calculate the microseismic signal Corresponding sensors and microseismic signals 、 Corresponding linear distance of the sensor , that is, sensor Respectively with the sensor and sensors Straight-line distance , and then select the minimum straight-line distance from multiple straight-line distances, for example, =33m, =31m, the minimum straight-line distance selected at this time is 31m, and then determine whether the minimum straight-line distance is less than the wavelength of the P wave. If so, then select As an auxiliary microseismic signal, and renamed .

[0109] It should be understood that reference Figure 10 , Figure 10 The overall process diagram is as follows: the microseismic signals received by each sensor are named according to the excitation time. ,exist When the P wave arrives, the signal is Fourier transformed, the main frequency of the signal is picked up according to the transformation result, and the wavelength and period of the P wave are calculated, and then the microseismic signal is sequentially Corresponding sensors and microseismic signals … Corresponding linear distance of the sensor and , where superscripts 1 and 2 indicate that the two sensors are in different groups or the same group, and , then take the straight-line distance and Then determine whether the minimum straight-line distance is less than the wavelength of the P wave. If so, the microseismic signal corresponding to the lowest straight-line distance is used as the auxiliary microseismic signal. Then, the Hilbert transform is performed on the P-wave arrival time solution signal and the auxiliary microseismic signal to obtain the instantaneous phase of the P-wave arrival time solution signal and the instantaneous phase of the auxiliary microseismic signal. Then, the instantaneous phase difference between the instantaneous phase of the P-wave arrival time solution signal and the instantaneous phase of the auxiliary microseismic signal is calculated to obtain the instantaneous phase difference. Then, the maximum instantaneous phase difference intensity is selected from multiple instantaneous phase difference intensities, and the time corresponding to the maximum instantaneous phase difference intensity is used as the arrival time of the P wave to solve the signal P wave arrival time Then continue to judge whether there are other signals that need to be picked up when the P wave arrives. If so, then , and pick up the arrival time of the P wave of the signal according to the above steps until the arrival time of the P wave of all signals is picked up.

[0110] This embodiment obtains the microseismic signal under the microseismic event to determine the time when the sensor corresponding to the microseismic signal is excited; sorts the microseismic signal according to the time in a preset time sequence; names the sorted microseismic signals in sequence according to a preset naming rule; and determines the P-wave arrival time solution signal and the auxiliary microseismic signal based on the named microseismic signal; in the above manner, the time when the sensor corresponding to the microseismic signal under the microseismic event is excited is determined, and then the microseismic signals are sorted according to the preset time sequence, and the sorted microseismic signals are named, and then the named microseismic signals are used to determine the P-wave arrival time solution signal and the auxiliary microseismic signal, thereby effectively improving the accuracy of determining the P-wave arrival time solution signal and the auxiliary microseismic signal.

[0111] In addition, an embodiment of the present invention further proposes a storage medium, on which a microseismic signal P wave arrival time picking program is stored. When the microseismic signal P wave arrival time picking program is executed by a processor, the steps of the microseismic signal P wave arrival time picking method described above are implemented.

[0112] Since the storage medium adopts all the technical solutions of all the above embodiments, it has at least all the beneficial effects brought by the technical solutions of the above embodiments, which will not be described one by one here.

[0113] In addition, refer to Figure 11 The embodiment of the present invention further provides a device for picking up the arrival time of a microseismic signal P wave, the device comprising:

[0114] The acquisition module 10 is used to acquire microseismic signals under microseismic events, and obtain P-wave arrival time solution signals and auxiliary microseismic signals according to the microseismic signals.

[0115] The determination module 20 is configured to determine an ideal phase difference based on the P-wave arrival time solution signal and the auxiliary microseismic signal.

[0116] The determination module 20 is further configured to determine an actual phase difference based on the P-wave arrival time solution signal and the auxiliary microseismic signal.

[0117] The picking module 30 is used to pick up the arrival time of the P wave of the P wave arrival time solution signal according to the ideal phase difference and the actual phase difference.

[0118] This embodiment obtains the microseismic signal under the microseismic event, obtains the P-wave arrival time solution signal and the auxiliary microseismic signal according to the microseismic signal; determines the ideal phase difference according to the P-wave arrival time solution signal and the auxiliary microseismic signal; determines the actual phase difference according to the P-wave arrival time solution signal and the auxiliary microseismic signal; picks up the arrival time of the P wave of the P-wave arrival time solution signal according to the ideal phase difference and the actual phase difference; through the above method, the P-wave arrival time solution signal and the auxiliary microseismic signal are selected from the microseismic signal under the microseismic event, and then the ideal phase difference and the actual phase difference are determined respectively, and then the arrival time of the P wave of the P-wave arrival time solution signal is picked up according to the ideal phase difference and the actual phase difference, thereby effectively improving the accuracy of picking up the P-wave arrival time.

[0119] It should be noted that the workflow described above is merely illustrative and does not limit the scope of protection of the present invention. In practical applications, technicians in this field can select part or all of it according to actual needs to achieve the purpose of the embodiment scheme, and no limitation is made here.

[0120] In addition, for technical details not fully described in this embodiment, reference can be made to the method for picking the arrival time of the microseismic signal P wave provided in any embodiment of the present invention, and will not be repeated here.

[0121] In one embodiment, the acquisition module 10 is further used to acquire microseismic signals under microseismic events, determine the time when the sensor corresponding to the microseismic signal is excited; sort the microseismic signals according to the time in a preset time sequence; name the sorted microseismic signals in sequence according to a preset naming rule; and determine the P-wave arrival time solution signal and the auxiliary microseismic signal based on the named microseismic signals.

[0122] In one embodiment, the acquisition module 10 is further used to select a P-wave arrival time solution signal from the named microseismic signal; perform Fourier transform on the P-wave arrival time solution signal, and pick up the signal main frequency according to the transformation result; calculate the wavelength of the P wave according to the signal main frequency; and select an auxiliary microseismic signal from the named microseismic signal according to the wavelength of the P wave.

[0123] In one embodiment, the acquisition module 10 is further used to respectively calculate the straight-line distances between the sensor corresponding to the current microseismic signal and the multiple sensors corresponding to the other multiple microseismic signals; determine the minimum straight-line distance among the multiple straight-line distances; and when the minimum straight-line distance is less than the wavelength of the P wave, use the microseismic signal corresponding to the minimum straight-line distance as an auxiliary microseismic signal.

[0124] In one embodiment, the determination module 20 is further used to calculate the straight-line distance between the sensor corresponding to the P-wave arrival time solution signal and the sensor corresponding to the auxiliary microseismic signal; obtain the propagation velocity of the P wave in the geological medium; and calculate the ideal phase difference based on the straight-line distance, propagation velocity, and signal main frequency between the sensor corresponding to the P-wave arrival time solution signal and the sensor corresponding to the auxiliary microseismic signal.

[0125] In one embodiment, the determination module 20 is further configured to perform Hilbert transform on the P-wave arrival time solution signal and the auxiliary microseismic signal, respectively, to obtain the instantaneous phase of the P-wave arrival time solution signal and the instantaneous phase of the auxiliary microseismic signal; and perform difference calculation on the instantaneous phase of the P-wave arrival time solution signal and the instantaneous phase of the auxiliary microseismic signal to obtain an actual phase difference.

[0126] In one embodiment, the pickup module 30 is also used to calculate the period of the P wave based on the main frequency of the signal; calculate the instantaneous phase difference intensity based on the period of the P wave, the ideal phase difference and the actual phase difference; select the maximum instantaneous phase difference intensity from the multiple instantaneous phase difference intensities, and use the time corresponding to the maximum instantaneous phase difference intensity as the arrival time of the P wave of the P wave arrival solution signal.

[0127] Other embodiments or implementation methods of the microseismic signal P wave arrival time pickup device of the present invention can refer to the above-mentioned method embodiments, which are not repeated here.

[0128] In addition, it should be noted that, in this document, the terms "comprises," "includes," or any other variations thereof are intended to encompass non-exclusive inclusion, such that a process, method, article, or system comprising a series of elements includes not only those elements but also other elements not explicitly listed, or elements inherent to such process, method, article, or system. In the absence of further limitations, an element defined by the phrase "comprising a ..." does not exclude the presence of other identical elements in the process, method, article, or system comprising the element.

[0129] The serial numbers of the above embodiments of the present invention are for description only and do not represent the advantages or disadvantages of the embodiments.

[0130] Through the above description of the embodiments, those skilled in the art will clearly understand that the methods of the above embodiments can be implemented using software plus the necessary general-purpose hardware platform. Of course, hardware can also be used, but in many cases the former is a more preferred implementation method. Based on this understanding, the technical solution of the present invention, or the portion 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, a magnetic disk, or an optical disk) and includes a number of instructions for enabling a terminal device (which can be a mobile phone, a computer, an integrated platform workstation, or a network device, etc.) to execute the methods described in the various embodiments of the present invention.

[0131] The above are only preferred embodiments of the present invention and are not intended to limit the patent scope of the present invention. Any equivalent structure or equivalent process transformation made using the contents of the present invention description and drawings, or directly or indirectly applied in other related technical fields, are also included in the patent protection scope of the present invention.

Claims

1. A method for picking up the arrival time of microseismic signal P wave, characterized in that: The method for picking up the arrival time of microseismic signal P wave comprises the following steps: Acquire a microseismic signal under a microseismic event, and obtain a P-wave arrival time solution signal and an auxiliary microseismic signal according to the microseismic signal; Determining an ideal phase difference based on the P-wave arrival time solution signal and the auxiliary microseismic signal; determining an actual phase difference based on the P-wave arrival time solution signal and the auxiliary microseismic signal; The arrival time of the P wave of the P wave arrival time solution signal is picked up according to the ideal phase difference and the actual phase difference.

2. The method for picking up the arrival time of microseismic signal P wave according to claim 1, characterized in that: The step of acquiring a microseismic signal under a microseismic event and obtaining a P-wave arrival time solution signal and an auxiliary microseismic signal according to the microseismic signal includes: Acquire a microseismic signal under a microseismic event, and determine a time when a sensor corresponding to the microseismic signal is excited; sorting the microseismic signals according to the time in a preset time sequence; Name the sorted microseismic signals in sequence according to the preset naming rules; The P-wave arrival time solution signal and auxiliary microseismic signal are determined based on the named microseismic signal.

3. The method for picking up the arrival time of microseismic signal P wave according to claim 2, characterized in that: Determining the P-wave arrival time solution signal and the auxiliary microseismic signal according to the named microseismic signal includes: Select the P-wave arrival time solution signal from the named microseismic signals; Performing Fourier transform on the P-wave arrival time solution signal, and picking up the signal main frequency according to the transformation result; Calculating the wavelength of the P wave based on the main frequency of the signal; An auxiliary microseismic signal is selected from the named microseismic signals according to the wavelength of the P wave.

4. The method for picking up the arrival time of microseismic signal P wave according to claim 3, characterized in that: The microseismic signal under the microseismic event includes the current microseismic signal and multiple other microseismic signals; The selecting of the auxiliary microseismic signal from the named microseismic signal according to the wavelength of the P wave includes: respectively calculating straight-line distances between the sensor corresponding to the current microseismic signal and the multiple sensors corresponding to the other multiple microseismic signals; determining a minimum straight-line distance among a plurality of said straight-line distances; When the minimum straight-line distance is smaller than the wavelength of the P wave, the microseismic signal corresponding to the minimum straight-line distance is used as an auxiliary microseismic signal.

5. The method for picking up the arrival time of microseismic signal P wave according to claim 1, characterized in that: Determining the ideal phase difference according to the P-wave arrival time solution signal and the auxiliary microseismic signal includes: Calculating the straight-line distance between the sensor corresponding to the P-wave arrival time solution signal and the sensor corresponding to the auxiliary microseismic signal; Obtain the propagation velocity of P waves in geological media; The ideal phase difference is calculated based on the straight-line distance, propagation velocity, and signal main frequency between the sensor corresponding to the P-wave arrival time solution signal and the sensor corresponding to the auxiliary microseismic signal.

6. The method for picking up the arrival time of microseismic signal P wave according to claim 1, characterized in that: The determining of the actual phase difference according to the P-wave arrival time solution signal and the auxiliary microseismic signal includes: Performing Hilbert transform on the P-wave arrival time solution signal and the auxiliary microseismic signal respectively to obtain the instantaneous phase of the P-wave arrival time solution signal and the instantaneous phase of the auxiliary microseismic signal; The instantaneous phase of the P-wave arrival signal and the instantaneous phase of the auxiliary microseismic signal are subtracted to obtain an actual phase difference.

7. The method for picking up the arrival time of microseismic signals P waves according to any one of claims 1 to 6, characterized in that: The step of picking up the P-wave arrival time of the P-wave arrival time solution signal according to the ideal phase difference and the actual phase difference comprises: Calculate the period of the P wave based on the dominant frequency of the signal; Calculating instantaneous phase difference intensity according to the period of the P wave, the ideal phase difference, and the actual phase difference; The maximum instantaneous phase difference intensity is selected from the multiple instantaneous phase difference intensities, and the time corresponding to the maximum instantaneous phase difference intensity is used as the arrival time of the P wave of the P wave arrival time solution signal.

8. A device for picking up the arrival time of microseismic signal P wave, characterized in that: The microseismic signal P wave arrival time picking device comprises: An acquisition module is used to acquire a microseismic signal under a microseismic event, and obtain a P-wave arrival time solution signal and an auxiliary microseismic signal according to the microseismic signal; a determination module, configured to determine an ideal phase difference based on the P-wave arrival time solution signal and the auxiliary microseismic signal; The determination module is further configured to determine an actual phase difference based on the P-wave arrival time solution signal and the auxiliary microseismic signal; A picking module is used to pick up the arrival time of the P wave of the P wave arrival solution signal according to the ideal phase difference and the actual phase difference.

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

10. A storage medium, characterized in that: The storage medium stores a microseismic signal P wave arrival time picking program, which, when executed by the processor, implements the microseismic signal P wave arrival time picking method according to any one of claims 1 to 7.

Images