Real-time monitoring method and device for micro-seismic of mining land

By synchronizing the downhole sensors with the surface sensors and processing the microseismic waveforms, the coordinates of the earthquake source were determined and iteratively located, which solved the problem of low accuracy in downhole microseismic monitoring and improved the accuracy of rockburst prediction and forecasting.

CN121784822APending Publication Date: 2026-04-03CCTEG COAL MINING RES INST +1
View PDF 0 Cites 0 Cited by

Patent Information

Authority / Receiving Office
CN · China
Patent Type
Applications(China)
Current Assignee / Owner
Filing Date
2025-11-18
Publication Date
2026-04-03

AI Technical Summary

Technical Problem

The existing technology for monitoring microseismic activity in mines is limited to a single mode, resulting in incomplete perception of the vibration field in and around the mine, which reduces the accuracy of microseismic location and affects the research and prevention of rockbursts.

Method used

By synchronizing downhole sensors with surface sensors to acquire sensor data, estimating the source coordinates based on microseismic waveforms, and performing microseismic positioning iterations, the positioning accuracy is improved.

Benefits of technology

It has effectively improved the accuracy of microseismic location and promoted the development of rockburst prediction, forecasting and prevention technologies.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN121784822A_ABST
    Figure CN121784822A_ABST
Patent Text Reader

Abstract

The invention discloses a mining ground micro-seismic real-time monitoring method and device. The mining ground micro-seismic real-time monitoring method comprises the steps that time service synchronization is conducted on an underground sensor and a ground sensor; sensing data collected by an underground sensor and a ground sensor are obtained; based on a microseismic waveform in the sensing data, determining an estimated seismic source coordinate; and performing micro-seismic positioning iteration based on the estimated seismic source coordinates to obtain a target seismic source position. According to the method, the estimated seismic source coordinates are determined based on the sensing data acquired by the underground sensor and the ground sensor, and the micro-seismic positioning iteration is performed based on the estimated seismic source coordinates to obtain the target seismic source position, so that the micro-seismic positioning precision is effectively improved, and the method plays a positive role in promoting the development of rock burst prediction and forecast and prevention and treatment technologies.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] This invention relates to the field of microseismic monitoring technology, and in particular to a method and device for real-time monitoring of microseismic activity in mining areas. Background Technology

[0002] During underground coal mining, the original stress balance is disrupted after the coal seam is extracted. Stress is redistributed and acts on the surrounding coal and rock mass and the roadway support structure, thus creating mine pressure. When the mine pressure borne by the surrounding rock exceeds the support resistance provided by the support structure, it can lead to roadway deformation, floor heave, spalling, or even complete collapse, easily triggering serious safety accidents and threatening miners' lives and normal mine production. Therefore, mine monitoring is necessary to determine the coordinates of seismic sources and prevent safety accidents.

[0003] In related technologies, monitoring underground microseismic activity using a single monitoring mode results in incomplete perception of the overall vibration field of the mine and its surrounding geological environment, thereby reducing the accuracy of microseismic location and affecting the research and prevention of mine rockbursts. Summary of the Invention

[0004] The present invention aims to at least partially solve one of the technical problems in the related art.

[0005] To address this, the present invention proposes a method and device for real-time monitoring of microseismic activity in mining areas. This method can determine the estimated source coordinates based on sensor data collected by downhole and surface sensors, and perform iterative microseismic positioning based on the estimated source coordinates to obtain the target source location. This effectively improves the accuracy of microseismic positioning and will play a positive role in promoting the development of rockburst prediction, forecasting and prevention technologies.

[0006] To achieve the above objectives, the present invention provides a method for real-time monitoring of microseismic activity in mining areas, comprising: Synchronize the downhole sensors with the surface sensors; Acquire the sensing data collected by the downhole sensor and the surface sensor; Based on the microseismic waveforms in the sensor data, the estimated source coordinates are determined; Based on the estimated source coordinates, microseismic location iterations are performed to obtain the target source location.

[0007] The real-time monitoring method for microseismic activity in mining areas according to embodiments of the present invention may also have the following additional technical features: In one embodiment of the present invention, the time synchronization between the downhole sensor and the surface sensor includes: The reference time is obtained through a navigation satellite receiver; Calculate the time reference deviation of each sensor based on the aforementioned reference time; Each of the sensors adjusts its own clock based on the time base deviation.

[0008] In one embodiment of the present invention, determining the estimated source coordinates based on the microseismic waveform in the sensing data includes: Acquire the micro-vibration waveform from the sensor data; The micro-vibration waveform is preprocessed to obtain the processed micro-vibration waveform; Based on the processed micro-vibration waveform, the waveform gradient corresponding to each sensor is determined. Based on the waveform gradient, the relative distance between the seismic source and the waveform source observed by each of the sensors is determined; The source depth corresponding to each sensor is calculated based on the relative distance; The estimated source coordinates are determined based on the relative distance and the source depth.

[0009] In one embodiment of the present invention, the step of performing microseismic location iteration based on the estimated source coordinates to obtain the target source location includes: Obtain the iteration parameters; The waveform gradient residual for estimating the source coordinates is calculated based on the wave velocity. Determine whether the waveform gradient residual meets the stopping iteration condition; If it is determined that the waveform gradient residual does not meet the stopping iteration condition, the estimated source coordinates and wave velocity are updated based on the iteration parameters to obtain the updated source coordinates and updated wave velocity. Based on the updated source coordinates and the updated wave velocity, repeat the above steps to determine the source coordinates that satisfy the stopping iteration condition or the source coordinates that satisfy the number of iterations as the target source location.

[0010] To achieve the above objectives, another aspect of the present invention provides a real-time monitoring device for microseismic activity in mining areas, the device comprising: The synchronization module is used to synchronize the downhole sensors with the surface sensors. The acquisition module is used to acquire the sensing data collected by the downhole sensor and the surface sensor; The first determining module is used to determine the estimated source coordinates based on the micro-seismic waveform in the sensing data; The second determining module is used to perform microseismic positioning iterations based on the estimated seismic source coordinates to obtain the target seismic source location.

[0011] In one embodiment of the present invention, the synchronization module is specifically used for: The reference time is obtained through a navigation satellite receiver; Calculate the time reference deviation of each sensor based on the aforementioned reference time; Each of the sensors adjusts its own clock based on the time base deviation.

[0012] In one embodiment of the present invention, the first determining module is specifically used for: Acquire the micro-vibration waveform from the sensor data; The micro-vibration waveform is preprocessed to obtain the processed micro-vibration waveform; Based on the processed micro-vibration waveform, the waveform gradient corresponding to each sensor is determined. Based on the waveform gradient, the relative distance between the seismic source and the waveform source observed by each of the sensors is determined; The source depth corresponding to each sensor is calculated based on the relative distance; The estimated source coordinates are determined based on the relative distance and the source depth.

[0013] In one embodiment of the present invention, the second determining module is specifically used for: Obtain the iteration parameters; The waveform gradient residual for estimating the source coordinates is calculated based on the wave velocity. Determine whether the waveform gradient residual meets the stopping iteration condition; If it is determined that the waveform gradient residual does not meet the stopping iteration condition, the estimated source coordinates and wave velocity are updated based on the iteration parameters to obtain the updated source coordinates and updated wave velocity. Based on the updated source coordinates and the updated wave velocity, repeat the above steps to determine the source coordinates that satisfy the stopping iteration condition or the source coordinates that satisfy the number of iterations as the target source location.

[0014] The real-time microseismic monitoring method and apparatus for mining areas according to embodiments of the present invention synchronizes downhole sensors and surface sensors; acquires sensor data collected by downhole and surface sensors; determines estimated source coordinates based on microseismic waveforms in the sensor data; and performs iterative microseismic positioning based on the estimated source coordinates to obtain the target source location. Therefore, the present invention can determine estimated source coordinates based on sensor data collected by downhole and surface sensors, and perform iterative microseismic positioning based on the estimated source coordinates to obtain the target source location, thereby effectively improving the accuracy of microseismic positioning and playing a positive role in promoting the development of rockburst prediction, forecasting, and prevention technologies.

[0015] Additional aspects and advantages of the invention will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. Attached Figure Description

[0016] The above and / or additional aspects and advantages of the present invention will become apparent and readily understood from the following description of the embodiments taken in conjunction with the accompanying drawings, wherein: Figure 1 This is a flowchart illustrating a method for real-time monitoring of microseismic activity in mining areas according to an embodiment of the present invention. Figure 2 This is a schematic diagram of a real-time monitoring device for microseismic activity in mining areas according to an embodiment of the present invention. Detailed Implementation

[0017] It should be noted that, unless otherwise specified, the embodiments and features described in the present invention can be combined with each other. The present invention will now be described in detail with reference to the accompanying drawings and embodiments.

[0018] To enable those skilled in the art to better understand the present invention, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort should fall within the scope of protection of the present invention.

[0019] The following description, with reference to the accompanying drawings, describes a method and apparatus for real-time monitoring of microseismic activity in mining areas according to embodiments of the present invention.

[0020] Figure 1 This is a flowchart illustrating the real-time monitoring method for microseismic activity in mining areas according to an embodiment of the present invention.

[0021] like Figure 1 As shown, the method may include the following steps: Step 101: Synchronize the downhole sensors with the surface sensors.

[0022] In one embodiment of the present invention, the method for synchronizing downhole sensors and surface sensors may include the following steps: Step 1011: Obtain the reference time through the navigation satellite receiver; Step 1012: Calculate the time reference deviation of each sensor based on the reference time; Step 1013: Each sensor adjusts its own clock based on the time base deviation.

[0023] In one embodiment of the present invention, the navigation satellite may include either GPS or BeiDou.

[0024] In one embodiment of the present invention, the downhole sensor and the surface sensor are respectively equipped with independent clocks, and each sensor can determine the time reference deviation of each sensor by the difference between the independent clock and the reference time.

[0025] In one embodiment of the present invention, after each sensor obtains the corresponding time reference deviation through the above steps, it can adjust its own clock based on the time reference deviation so as to unify the time of the sensor and the navigation satellite receiver.

[0026] Step 102: Acquire sensor data collected by downhole sensors and surface sensors.

[0027] In one embodiment of the present invention, the above-mentioned sensing data may include micro-vibration waveforms.

[0028] Step 103: Determine the estimated source coordinates based on the microseismic waveforms in the sensor data.

[0029] In one embodiment of the present invention, after obtaining the sensing data through the above steps, the estimated source coordinates can be determined based on the microseismic waveforms in the sensing data.

[0030] In one embodiment of the present invention, the method for determining the estimated source coordinates based on microseismic waveforms in sensor data may include the following steps: Step 1031: Obtain the micro-vibration waveform from the sensor data.

[0031] In one embodiment of the present invention, after obtaining the micro-vibration waveform in the sensing data through the above steps, the data corresponding to the obtained micro-vibration waveform can be numbered.

[0032] Step 1032: Preprocess the micro-vibration waveform to obtain the processed micro-vibration waveform.

[0033] In one embodiment of the present invention, the method for preprocessing the micro-vibration waveform to obtain the processed micro-vibration waveform may include: performing noise reduction processing and redundancy removal processing on the micro-vibration waveform to obtain the processed micro-vibration waveform.

[0034] Specifically, in one embodiment of the present invention, the above-mentioned denoising process can be wavelet threshold denoising (such as db4 wavelet, retaining the 10Hz-1kHz frequency band). Also, in one embodiment of the present invention, the above-mentioned redundancy removal process can include truncating the effective time period before and after the P-wave's first arrival point to avoid S-wave interference.

[0035] Step 1033: Based on the processed micro-vibration waveform, determine the waveform gradient corresponding to each sensor.

[0036] In one embodiment of the present invention, after obtaining the processed micro-vibration waveform through the above steps, the waveform gradient corresponding to each sensor can be determined based on the processed micro-vibration waveform.

[0037] In one embodiment of the present invention, the waveform gradient corresponding to each sensor can be determined based on a first formula, wherein the first formula is:

[0038] Wherein, the waveform x of the k-th sensor k (t) Calculate the time derivative, taking the gradient value G at the initial arrival time tp of wave P. k .

[0039] Step 1034: Based on the waveform gradient, determine the relative distance between the seismic source and the waveform sources observed by each sensor.

[0040] In one embodiment of the present invention, after determining the waveform gradient through the above steps, the relative distance between the seismic source and the waveform sources observed by each sensor can be determined based on the waveform gradient.

[0041] Specifically, in one embodiment of the present invention, the method for determining the relative distance between the seismic source and the waveform sources observed by each sensor based on the waveform gradient may include: determining the relative distance between the seismic source and the waveform sources observed by each sensor based on the waveform gradient using a second formula, wherein the second formula is: Δr k =α×v p ×Δt k / (v p -v s ) Where, Δr k The relative distance between the waveform sources observed by the k-th sensor and the source is α, where α is the rock stratum wave velocity correction coefficient, and v p For P-wave velocity, v s Let Δt be the S-wave velocity. k Let be the travel time of the P-wave of the k-th sensor relative to the reference reference.

[0042] Furthermore, in one embodiment of the present invention, the waveform gradient G corresponding to each sensor can be obtained through the above steps. k Then, the travel time difference Δt can be calculated. k , Δt k = G k -G p0 Among them, G p0 It is the earliest arriving P-wave gradient value.

[0043] Step 1035: Calculate the source depth corresponding to each sensor based on the relative distance.

[0044] In one embodiment of the present invention, after obtaining the relative distance between the waveform sources observed by each sensor through the above steps, the source depth corresponding to each sensor can be calculated based on the relative distance.

[0045] In one embodiment of the present invention, the method for calculating the source depth of each sensor based on relative distance may include: calculating the source depth of each sensor based on relative distance using a third formula, wherein the third formula is:

[0046] Among them, h k Let k be the source depth of the seismic source. For ground reference depth, Let be the horizontal projected distance between the k-th sensor and the seismic source.

[0047] Step 1036: Determine the estimated source coordinates based on the relative distance and source depth.

[0048] In one embodiment of the present invention, after determining the relative distance and the source depth through the above steps, the estimated source coordinates can be determined based on the relative distance and the source depth.

[0049] In one embodiment of the present invention, the method for determining and estimating the source coordinates based on relative distance and source depth may include: obtaining the horizontal coordinates of the source by intersecting the relative distances of the waveforms observed by each sensor using the least squares method. The source depth h of each sensor k The average value is determined as the vertical coordinate. The estimated source coordinates are obtained. ).

[0050] Step 104: Perform microseismic location iterations based on the estimated source coordinates to obtain the target source location.

[0051] In one embodiment of the present invention, after determining the estimated source coordinates through the above steps, the target source location can be obtained by performing microseismic positioning iterations based on the estimated source coordinates.

[0052] In one embodiment of the present invention, the method for obtaining the target seismic source location by performing microseismic location iteration based on estimated source coordinates may include the following steps: Step 1041: Obtain the iteration parameters.

[0053] In one embodiment of the present invention, the set iteration parameters can be obtained. In one embodiment of the present invention, the iteration parameters may include an iteration step size λ = 2m, an iteration count n = 100, and a speed step size v = 0.05m / s.

[0054] Step 1042: Calculate and estimate the waveform gradient residual of the source coordinates based on wave velocity.

[0055] In one embodiment of the present invention, the method for estimating the waveform gradient residual of the seismic source coordinates based on wave velocity may include: calculating the waveform gradient residual for estimating the seismic source coordinates based on the wave velocity using a fourth formula, wherein the fourth formula is:

[0056] In one embodiment of the present invention, en is the waveform gradient parameter of the nth iteration. When the P-wave arrives at the k-th sensor, To estimate the source coordinates (x) n , z n ) and wave speed The theory of calculation will be available by then. , where (x k , z k () represents the coordinates of the k-th sensor.

[0057] Step 1043: Determine whether the waveform gradient residual meets the stopping iteration condition.

[0058] In one embodiment of the present invention, after obtaining the waveform gradient residual through the above steps, it can be determined whether the waveform gradient residual meets the stopping iteration condition.

[0059] In one embodiment of the present invention, the waveform gradient residual is used to determine whether it meets the stopping iteration condition using a fifth formula, wherein the fifth formula is: ε≤| |≤δ Wherein, the gradient tolerance δ = 2m and the error tolerance ε = 0.1m, if the waveform gradient residual satisfies the fifth formula above, then the waveform gradient residual is determined to meet the stopping iteration condition; otherwise, the waveform gradient residual is determined to not meet the stopping iteration condition.

[0060] In one embodiment of the present invention, if it is determined that the waveform gradient residual does not meet the stopping iteration condition, the source coordinates that meet the stopping iteration condition are determined as the target source location.

[0061] Step 1044: If it is determined that the waveform gradient residual does not meet the stopping iteration condition, then the source coordinates and wave velocity are updated based on the iteration parameters to obtain the updated source coordinates and the updated wave velocity.

[0062] In one embodiment of the present invention, if it is determined that the waveform gradient residual does not meet the stopping iteration condition, the estimated source coordinates and wave velocity are updated based on the iteration parameters to obtain the updated source coordinates and the updated wave velocity.

[0063] Specifically, in one embodiment of the present invention, the method for updating the estimated source coordinates and wave velocity based on iterative parameters to obtain the updated source coordinates and wave velocity may include: updating the estimated source coordinates and wave velocity based on iterative parameters using a sixth formula to obtain the updated source coordinates and wave velocity, wherein the sixth formula is: ,

[0064]

[0065] Among them, the above (x) n+1 , z n+1 () represents the updated source coordinates. This is the updated wave speed.

[0066] Step 1045: Based on the updated source coordinates and the updated wave velocity, repeat the above steps to determine the source coordinates that meet the stopping iteration condition or the source coordinates that meet the number of iterations as the target source location.

[0067] In one embodiment of the present invention, after obtaining the updated source coordinates and the updated wave velocity through the above steps, steps 1042 to 1044 can be repeated, and the source coordinates that meet the stopping iteration condition or the source coordinates that meet the iteration number can be determined as the target source location.

[0068] The real-time microseismic monitoring method for mining areas according to embodiments of the present invention synchronizes downhole sensors and surface sensors in terms of time; acquires sensor data collected by downhole and surface sensors; determines estimated source coordinates based on microseismic waveforms in the sensor data; and performs iterative microseismic positioning based on the estimated source coordinates to obtain the target source location. Therefore, the present invention can determine estimated source coordinates based on sensor data collected by downhole and surface sensors, and perform iterative microseismic positioning based on the estimated source coordinates to obtain the target source location, thereby effectively improving the accuracy of microseismic positioning and playing a positive role in promoting the development of rockburst prediction, forecasting, and prevention technologies.

[0069] Figure 2 This is a schematic diagram of the structure of the real-time monitoring device for microseismic activity in mining areas according to an embodiment of the present invention.

[0070] like Figure 2 As shown, the device may include: Synchronization module 201 is used to synchronize downhole sensors with surface sensors. The acquisition module 202 is used to acquire sensor data collected by downhole sensors and surface sensors; The first determining module 203 is used to determine the estimated source coordinates based on the microseismic waveforms in the sensor data. The second determining module 204 is used to perform microseismic location iteration based on the estimated source coordinates to obtain the target source location.

[0071] In one embodiment of the present invention, the synchronization module 201 is specifically used for: The reference time is obtained through a navigation satellite receiver; Calculate the time reference deviation of each sensor based on the reference time; Each sensor adjusts its own clock based on the time base deviation.

[0072] In one embodiment of the present invention, the first determining module 203 is specifically used for: Acquire micro-vibration waveforms from sensor data; The micro-vibration waveform is preprocessed to obtain the processed micro-vibration waveform; Based on the processed micro-vibration waveform, the waveform gradient corresponding to each sensor is determined. Based on the waveform gradient, the relative distance between the seismic source and the waveform sources observed by each sensor is determined; The source depth corresponding to each sensor is calculated based on the relative distance; Based on the relative distance and focal depth, the estimated focal coordinates are determined.

[0073] In one embodiment of the present invention, the second determining module 204 is specifically used for: Obtain the iteration parameters; Waveform gradient residuals for estimating source coordinates based on wave velocity calculation; Determine whether the waveform gradient residual meets the stopping iteration condition; If it is determined that the waveform gradient residual does not meet the stopping iteration condition, the source coordinates and wave velocity are updated based on the iteration parameters to obtain the updated source coordinates and the updated wave velocity. Based on the updated source coordinates and the updated wave velocity, repeat the above steps to determine the source coordinates that meet the stopping iteration condition or the source coordinates that meet the number of iterations as the target source location.

[0074] The real-time microseismic monitoring device for mining areas according to this invention synchronizes the timing of downhole sensors and surface sensors; acquires sensor data collected by both downhole and surface sensors; determines estimated seismic source coordinates based on the microseismic waveforms in the sensor data; and performs iterative microseismic positioning based on the estimated seismic source coordinates to obtain the target seismic source location. Therefore, this invention can determine estimated seismic source coordinates based on sensor data collected by downhole and surface sensors, and perform iterative microseismic positioning based on these coordinates to obtain the target seismic source location, thereby effectively improving the accuracy of microseismic positioning and playing a positive role in promoting the development of rockburst prediction, forecasting, and prevention technologies.

[0075] In this specification, the use of terms such as "one embodiment," "some embodiments," "example," "specific example," or "some examples," etc., refers to a specific feature, structure, material, or characteristic described in connection with that embodiment or example, which is included in at least one embodiment or example of the invention. In this specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in one or more embodiments or examples. Moreover, without contradiction, those skilled in the art can combine and integrate the different embodiments or examples described in this specification, as well as the features of different embodiments or examples.

[0076] Furthermore, the terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include at least one of that feature. In the description of this invention, "a plurality of" means at least two, such as two, three, etc., unless otherwise explicitly specified.

Claims

1. A method for real-time monitoring of microseismic activity in mining areas, characterized in that, include: Synchronize the downhole sensors with the surface sensors; Acquire the sensing data collected by the downhole sensor and the surface sensor; Based on the microseismic waveforms in the sensor data, the estimated source coordinates are determined; Based on the estimated source coordinates, microseismic location iterations are performed to obtain the target source location.

2. The method as described in claim 1, characterized in that, The time synchronization between downhole sensors and surface sensors includes: The reference time is obtained through a navigation satellite receiver; Calculate the time reference deviation of each sensor based on the aforementioned reference time; Each of the sensors adjusts its own clock based on the time base deviation.

3. The method as described in claim 1, characterized in that, The step of determining the estimated source coordinates based on the microseismic waveforms in the sensing data includes: Acquire the micro-vibration waveform from the sensor data; The micro-vibration waveform is preprocessed to obtain the processed micro-vibration waveform; Based on the processed micro-vibration waveform, the waveform gradient corresponding to each sensor is determined. Based on the waveform gradient, the relative distance between the seismic source and the waveform source observed by each of the sensors is determined; The source depth corresponding to each sensor is calculated based on the relative distance; The estimated source coordinates are determined based on the relative distance and the source depth.

4. The method as described in claim 1, characterized in that, The step of performing microseismic location iterations based on the estimated source coordinates to obtain the target source location includes: Obtain the iteration parameters; The waveform gradient residual for estimating the source coordinates is calculated based on the wave velocity. Determine whether the waveform gradient residual meets the stopping iteration condition; If it is determined that the waveform gradient residual does not meet the stopping iteration condition, the estimated source coordinates and wave velocity are updated based on the iteration parameters to obtain the updated source coordinates and updated wave velocity. Based on the updated source coordinates and the updated wave velocity, repeat the above steps to determine the source coordinates that satisfy the stopping iteration condition or the source coordinates that satisfy the number of iterations as the target source location.

5. A real-time monitoring device for microseismic activity in mining areas, characterized in that, include: The synchronization module is used to synchronize the downhole sensors with the surface sensors. The acquisition module is used to acquire the sensing data collected by the downhole sensor and the surface sensor; The first determining module is used to determine the estimated source coordinates based on the micro-seismic waveform in the sensing data; The second determining module is used to perform microseismic positioning iterations based on the estimated source coordinates to obtain the target source location.

6. The apparatus as claimed in claim 5, characterized in that, The synchronization module is specifically used for: The reference time is obtained through a navigation satellite receiver; Calculate the time reference deviation of each sensor based on the aforementioned reference time; Each of the sensors adjusts its own clock based on the time base deviation.

7. The apparatus as claimed in claim 5, characterized in that, The first determining module is specifically used for: Acquire the micro-vibration waveform from the sensor data; The micro-vibration waveform is preprocessed to obtain the processed micro-vibration waveform; Based on the processed micro-vibration waveform, the waveform gradient corresponding to each sensor is determined. Based on the waveform gradient, the relative distance between the seismic source and the waveform source observed by each of the sensors is determined; The source depth corresponding to each sensor is calculated based on the relative distance; The estimated source coordinates are determined based on the relative distance and the source depth.

8. The apparatus as claimed in claim 5, characterized in that, The second determining module is specifically used for: Obtain the iteration parameters; The waveform gradient residual for estimating the source coordinates is calculated based on the wave velocity. Determine whether the waveform gradient residual meets the stopping iteration condition; If it is determined that the waveform gradient residual does not meet the stopping iteration condition, the estimated source coordinates and wave velocity are updated based on the iteration parameters to obtain the updated source coordinates and updated wave velocity. Based on the updated source coordinates and the updated wave velocity, repeat the above steps to determine the source coordinates that satisfy the stopping iteration condition or the source coordinates that satisfy the number of iterations as the target source location.

9. An electronic device, comprising: At least one processor; as well as A memory communicatively connected to the at least one processor; wherein, The memory stores instructions that can be executed by the at least one processor to enable the at least one processor to perform the method of any one of claims 1-4.

10. A computer storage medium, wherein, The computer storage medium stores computer-executable instructions; when executed by a processor, the computer-executable instructions can implement the method described in any one of claims 1-4.