Method and device for updating seismic wave velocity model, electronic equipment and storage medium

Abstract

The application provides a seismic wave velocity model updating method and device, electronic equipment and storage medium. The method comprises the following steps: obtaining microseismic events collected by a geophone module, wherein the microseismic events comprise occurrence coordinates of the microseismic events and cutting times of the microseismic events; obtaining P-wave arrival times and S-wave arrival times of the microseismic events according to the microseismic events; updating a pre-constructed seismic wave velocity model according to the occurrence coordinates, the cutting times, the P-wave arrival times and the S-wave arrival times; and predicting travel times of seismic waves through the updated seismic wave velocity model. The application can correct the pre-constructed seismic wave velocity model through inherent microseismic events generated in mining work, and is a safe, environmentally-friendly and simple velocity model correction method which can replace explosive sources.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of microseismic, in particular to a method and device for updating a seismic wave velocity model, an electronic device and a storage medium. BACKGROUND

[0002] At present, the downhole directional hole fracturing of the mine stratum is mostly carried out in the open hole by the back-off fracturing, and no perforation is made, so the perforation signal cannot be used to correct the velocity model; although the velocity model correction can be carried out based on the calibration gun by using explosives. However, with the more strict control of explosives and other fire attack products in the coal mine, the multiple shooting in the fracturing area not only pollutes the downhole environment, has high cost, and also has safety hazards, and does not meet the requirements of the green mine construction.

[0003] Therefore, the problems in the prior art need to be solved. SUMMARY

[0004] The present application provides a method and device for updating a seismic wave velocity model, an electronic device and a storage medium, to solve the defects of the velocity model correction based on the calibration gun in the prior art, and to realize simple, green and low-cost velocity model correction.

[0005] The present application provides a method for updating a seismic wave velocity model, comprising:

[0006] Acquiring a microseismic event collected by a geophone module, the microseismic event comprising an occurrence coordinate of the microseismic event and a cutting time of the microseismic event;

[0007] According to the microseismic event, acquiring a P-wave arrival time and an S-wave arrival time of the microseismic event;

[0008] According to the occurrence coordinate, the cutting time, the P-wave arrival time and the S-wave arrival time, updating a pre-constructed seismic wave velocity model;

[0009] Predicting a travel time of a seismic wave through the updated seismic wave velocity model.

[0010] According to the method for updating a seismic wave velocity model provided by the present application, according to the occurrence coordinate and the cutting time, the P-wave arrival time and the S-wave arrival time of the microseismic event are acquired, which specifically comprises:

[0011] Data preprocessing is performed on the microseismic event, and the preprocessing comprises denoising processing, filtering processing and gain correction processing;

[0012] First arrival picking is performed on the preprocessed microseismic event to acquire the P-wave arrival time and the S-wave arrival time of the microseismic event.

[0013] The application provides a method for updating a seismic wave velocity model, which comprises the following steps of: inputting the occurrence coordinates into the velocity model to obtain a predicted travel time of a seismic wave; obtaining a true travel time of the seismic wave according to the occurrence coordinates, the cutting time, the P-wave arrival time and the S-wave arrival time; and updating the velocity model according to the predicted travel time and the true travel time.

[0014] The application provides a method for updating a seismic wave velocity model, which comprises the following steps of: inputting the occurrence coordinates into the velocity model to obtain a predicted travel time of a seismic wave; obtaining a true travel time of the seismic wave according to the occurrence coordinates, the cutting time, the P-wave arrival time and the S-wave arrival time; and updating the velocity model according to the predicted travel time and the true travel time.

[0015] The application provides a method for updating a seismic wave velocity model, which comprises the following steps of: inputting the occurrence coordinates into the velocity model to obtain a predicted travel time of a seismic wave; obtaining a true travel time of the seismic wave according to the occurrence coordinates, the cutting time, the P-wave arrival time and the S-wave arrival time; and updating the velocity model according to the predicted travel time and the true travel time.

[0016] The application provides a method for updating a seismic wave velocity model, which comprises the following steps of: inputting the occurrence coordinates into the velocity model to obtain a predicted travel time of a seismic wave; obtaining a true travel time of the seismic wave according to the occurrence coordinates, the cutting time, the P-wave arrival time and the S-wave arrival time; and updating the velocity model according to the predicted travel time and the true travel time.

[0017] The application provides a method for updating a seismic wave velocity model, which comprises the following steps of: inputting the occurrence coordinates into the velocity model to obtain a predicted travel time of a seismic wave; obtaining a true travel time of the seismic wave according to the occurrence coordinates, the cutting time, the P-wave arrival time and the S-wave arrival time; and updating the velocity model according to the predicted travel time and the true travel time.

[0018] The application provides a method for updating a seismic wave velocity model, which comprises the following steps of: inputting the occurrence coordinates into the velocity model to obtain a predicted travel time of a seismic wave; obtaining a true travel time of the seismic wave according to the occurrence coordinates, the cutting time, the P-wave arrival time and the S-wave arrival time; and updating the velocity model according to the predicted travel time and the true travel time.

[0019] The application provides a method for updating a seismic wave velocity model, which comprises the following steps of: inputting the occurrence coordinates into the velocity model to obtain a predicted travel time of a seismic wave; obtaining a true travel time of the seismic wave according to the occurrence coordinates, the cutting time, the P-wave arrival time and the S-wave arrival time; and updating the velocity model according to the predicted travel time and the true travel time.

[0020] The application provides a method for updating a seismic wave velocity model, which comprises the following steps of: inputting the occurrence coordinates into the velocity model to obtain a predicted travel time of a seismic wave; obtaining a true travel time of the seismic wave according to the occurrence coordinates, the cutting time, the P-wave arrival time and the S-wave arrival time; and updating the velocity model according to the predicted travel time and the true travel time.

[0021] The application provides a method for updating a seismic wave velocity model, which comprises the following steps of: inputting the occurrence coordinates into the velocity model to obtain a predicted travel time of a seismic wave; obtaining a true travel time of the seismic wave according to the occurrence coordinates, the cutting time, the P-wave arrival time and the S-wave arrival time; and updating the velocity model according to the predicted travel time and the true travel time.

[0022] The application provides a method for updating a seismic wave velocity model, which comprises the following steps of: inputting the occurrence coordinates into the velocity model to obtain a predicted travel time of a seismic wave; obtaining a true travel time of the seismic wave according to the occurrence coordinates, the cutting time, the P-wave arrival time and the S-wave arrival time; and updating the velocity model according to the predicted travel time and the true travel time.

[0023] The application provides a method for updating a seismic wave velocity model, which comprises the following steps of: inputting the occurrence coordinates into the velocity model to obtain a predicted travel time of a seismic wave; obtaining a true travel time of the seismic wave according to the occurrence coordinates, the cutting time, the P-wave arrival time and the S-wave arrival time; and updating the velocity model according to the predicted travel time and the true travel time.

[0024] The application provides a method for updating a seismic wave velocity model, which comprises the following steps of: inputting the occurrence coordinates into the velocity model to obtain a predicted travel time of a seismic wave; obtaining a true travel time of the seismic wave according to the occurrence coordinates, the cutting time, the P-wave arrival time and the S-wave arrival time; and updating the velocity model according to the predicted travel time and the true travel time.

[0025] The application provides a method for updating a seismic wave velocity model, which comprises the following steps of: inputting the occurrence coordinates into the velocity model to obtain a predicted travel time of a seismic wave; obtaining a true travel time of the seismic wave according to the occurrence coordinates, the cutting time, the P-wave arrival time and the S-wave arrival time; and updating the velocity model according to the predicted travel time and the true travel time.

[0026] The application provides a method for updating a seismic wave velocity model, which comprises the following steps of: inputting the occurrence coordinates into the velocity model to obtain a predicted travel time of a seismic wave; obtaining a true travel time of the seismic wave according to the occurrence coordinates, the cutting time, the P-wave arrival time and the S-wave arrival time; and updating the velocity model according to the predicted travel time and the true travel time.

[0027] A prediction unit is configured to predict the travel time of the seismic wave by using the updated seismic wave velocity model.

[0028] The application further provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor implements the method for updating the seismic wave velocity model according to any one of the above when executing the program.

[0029] The application further provides a non-transitory computer readable storage medium, which stores a computer program, and the computer program is executable on a processor to implement the method for updating the seismic wave velocity model according to any one of the above.

[0030] The application further provides a computer program product, which includes a computer program, and the computer program is executable on a processor to implement the method for updating the seismic wave velocity model according to any one of the above.

[0031] The application provides a method and device for updating a seismic wave velocity model, an electronic device, and a storage medium. The method includes obtaining microseismic events collected by a geophone module, the microseismic events including occurrence coordinates of the microseismic events and cutting times of the microseismic events; obtaining P-wave arrival times and S-wave arrival times of the microseismic events according to the microseismic events; updating a pre-constructed seismic wave velocity model according to the occurrence coordinates, the cutting times, the P-wave arrival times, and the S-wave arrival times; and predicting travel times of seismic waves by using the updated seismic wave velocity model. The application can correct the pre-constructed seismic wave velocity model by using inherent microseismic events generated in mining work, instead of using explosives to generate seismic waves to correct the velocity model, and is a safe, environmentally friendly, and simple velocity model correction method. BRIEF DESCRIPTION OF DRAWINGS

[0032] In order to more clearly illustrate the technical solutions of the present application or the prior art, the following will briefly introduce the drawings needed in the embodiments or prior art description. Obviously, the drawings in the following description are some embodiments of the present application, and those skilled in the art can obtain other drawings according to these drawings without creative labor.

[0033] Figure 1 is a flowchart of the method for updating the seismic wave velocity model provided by the application;

[0034] Figure 2 is a schematic diagram of the geophone module provided by the application;

[0035] Figure 3 is a structural schematic diagram of the device for updating the seismic wave velocity model provided by the application;

[0036] Figure 4 is a structural schematic diagram of an electronic device provided by the present application. DETAILED DESCRIPTION

[0037] To make the objectives, technical solutions and advantages of the present application clearer, the technical solutions in the present application will be described clearly and completely below with reference to the drawings in the present application. Obviously, the described embodiments are part of, rather than all of, the embodiments of the present application. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without any creative work fall within the protection scope of the present application.

[0038] At present, the downhole directional hole fracturing of the mine stratum is mostly carried out in the open hole in a retreating fracturing mode, and there is no perforation, so the perforation signal cannot be used to correct the velocity model; although the velocity model correction can be carried out based on the calibration gun through explosive excitation. However, with the more stringent control of explosives and other fire attack products in coal mines, multiple shooting in the fracturing area not only pollutes the downhole environment, has high cost, and also has safety hazards, and does not meet the requirements of green mine construction.

[0039] In order to solve the above problems, the present application provides a method for updating a seismic wave velocity model, which includes but is not limited to steps 110-140:

[0040] Step 110, acquiring the microseismic event collected by the geophone module, the microseismic event including the occurrence coordinates of the microseismic event and the cutting time of the microseismic event.

[0041] In step 110, the microseismic event of the to-be-detected area needs to be collected by the geophone module to obtain the occurrence coordinates of the microseismic event and the cutting time of the microseismic event.

[0042] A geophone is a device or apparatus used for detecting and measuring signals. In seismic monitoring, a geophone is also called a seismograph and is used to measure the vibration of seismic waves. The geophone module is a group of geophones arranged in the to-be-detected area, as shown in FIG. Figure 2 By arranging the geophone module, the microseismic event of the to-be-detected area can be collected.

[0043] Microseismic events refer to very small magnitude seismic activity on Earth. They are usually caused by the breaking or movement of rocks in the Earth's crust, but the energy release is very small in scale. Microseismic events are often recorded in seismic monitoring networks. These events are of great significance in seismic research and can help understand the movement of the Earth's crust and the patterns of seismic activity. In the present invention, microseismic events refer to the microseismic events generated by the cutting of a coal mining machine. The microseismic events generated by the cutting of a coal mining machine are a kind of interference signal for hydraulic fracturing microseismic monitoring. However, the present invention can utilize the inherent microseismic events generated by the cutting of a coal mining machine, without the need for using explosives underground to generate seismic waves to correct the velocity model.

[0044] Step 120, according to the microseismic event, obtaining the P-wave arrival time and S-wave arrival time of the microseismic event.

[0045] In step 120, the P-wave arrival time and S-wave arrival time need to be obtained according to the obtained microseismic event. The P-wave arrival time and S-wave arrival time of the microseismic event are important parameters for measuring the propagation time of seismic waves. P-wave (P-wave) is a compressive wave that propagates along the medium, while S-wave (S-wave) is a shear wave that propagates along the medium. P-wave arrival time refers to the time from the occurrence of a seismic event to the first arrival of P-wave at a seismic detection point. It is the initial stage of the propagation of seismic waves in the Earth's crust and is usually the first wave to be recorded. The arrival time information of P-wave can be used to determine the occurrence time of the earthquake and the approximate location of the earthquake source. S-wave arrival time refers to the time from the occurrence of a seismic event to the first arrival of S-wave at a seismic detection point. S-wave propagates slower than P-wave, so it arrives later in the sequence of seismic wave arrival. The arrival time information of S-wave can help determine the magnitude of the earthquake and the more accurate location of the earthquake source. By measuring the P-wave and S-wave arrival times at multiple seismic detection points, triangulation or other seismic positioning algorithms can be used to determine the location and depth of the earthquake source. This is very important for seismic research and seismic monitoring, and can be used for the operation of the earthquake early warning system.

[0046] Step 130, updating the pre-constructed seismic wave velocity model according to the occurrence coordinates, the cutting time, the P-wave arrival time and the S-wave arrival time.

[0047] In step 130, the true travel time measured by steps 110 and 120 needs to be compared with the predicted travel time predicted by the pre-constructed seismic wave velocity model, and the Levenberg-Marquard (LM) algorithm is used to continuously update the velocity model until the difference between the observed travel time and the forward travel time is minimized.

[0048] Specifically, the ray tracing forward method of the present application obtains the travel time of the seismic source to each geophone, i.e., the travel time of the inherent microseismic event generated by the coal mining machine cutting to each geophone. The ray tracing forward method is a commonly used method for simulating seismic wave propagation, which can be used to calculate the travel time of seismic waves from the seismic source to each geophone. The specific steps include:

[0049] Building a velocity model: First, a seismic wave velocity model needs to be established to describe the velocity distribution inside the earth. This model is usually obtained based on seismic observation data or geophysical exploration data.

[0050] Launching a ray: Starting from the seismic source, a beam of rays (i.e., seismic wave paths) is launched in all directions. The path of the ray will bend, refract or reflect according to the speed changes of different media in the seismic wave velocity model.

[0051] Calculating travel time: According to the ray path and the velocity model, the travel time of each ray is calculated. The travel time depends on the length of the ray and the speed of the different media it passes through.

[0052] Summarizing travel time: For each geophone, the travel times of all rays arriving at that geophone are summarized to obtain the total travel time of the seismic source to each geophone.

[0053] The travel time information obtained by the ray tracing forward method can be used for seismic positioning, seismic waveform simulation, seismic imaging and seismic velocity model verification, etc. This method has wide application value in seismology research and earthquake monitoring.

[0054] Step 140, predicting the travel time of the seismic wave through the updated seismic wave velocity model.

[0055] In step 140, the travel time of the seismic wave can be predicted through the updated seismic wave velocity model. It can be understood that the updated seismic wave velocity model is more accurate than the un-updated seismic wave velocity model.

[0056] The application provides a seismic wave velocity model updating method and device, electronic equipment and storage medium. The method comprises the following steps: obtaining a microseismic event collected by a geophone module, wherein the microseismic event comprises occurrence coordinates of the microseismic event and a cutting time of the microseismic event; obtaining P-wave arrival time and S-wave arrival time of the microseismic event according to the microseismic event; updating a pre-constructed seismic wave velocity model according to the occurrence coordinates, the cutting time, the P-wave arrival time and the S-wave arrival time; and predicting travel time of a seismic wave through the updated seismic wave velocity model. The application can correct the pre-constructed seismic wave velocity model through inherent microseismic events generated in mining work without using explosives to generate seismic waves for correcting the velocity model, and is a safe, environmentally friendly and simple velocity model correction method that can replace explosive sources.

[0057] According to the application, the P-wave arrival time and the S-wave arrival time of the microseismic event are obtained according to the occurrence coordinates and the cutting time, and specifically comprise:

[0058] The microseismic event is preprocessed, and the preprocessing comprises denoising processing, filtering processing and gain correction processing.

[0059] The preprocessed microseismic event is first break picked to obtain the P-wave arrival time and the S-wave arrival time of the microseismic event.

[0060] In this embodiment, in order to obtain more accurate P-wave arrival time and S-wave arrival time of the microseismic event, the microseismic event data needs to be preprocessed first to reduce the influence of noise, and then first break picking is performed to determine the arrival time of the P-wave and the S-wave of the seismic wave.

[0061] The data preprocessing of the microseismic event generally comprises denoising processing, filtering processing and gain correction processing. These steps help to extract the microseismic signal and reduce the influence of noise.

[0062] Denoising processing: the target of denoising processing is to reduce the interference of background noise on the microseismic signal. Common denoising methods include median filtering, wavelet denoising and frequency domain filtering, etc. These methods can select the most suitable method to reduce the noise level according to the specific situation.

[0063] Filtering processing: filtering processing is used to selectively enhance the frequency components of the microseismic signal and reduce the influence of other frequency components. A band-pass filter is generally used to filter out low-frequency and high-frequency noise to highlight the frequency characteristics of the microseismic event. Common filtering methods include low-pass filtering, high-pass filtering and band-pass filtering.

[0064] Gain correction processing: Gain correction processing aims to adjust the amplitude of microseismic signals for better visualization and more accurate first arrival picking. Common gain correction methods include linear gain correction and logarithmic gain correction. These methods can be adjusted according to the dynamic range of the signal to highlight the microseismic signal and maintain appropriate contrast.

[0065] First arrival picking of pre-processed microseismic events can be done manually or automatically.

[0066] Manual first arrival picking: In manual first arrival picking, seismologists carefully observe the pre-processed microseismic signal and mark the arrival times of P and S waves. This usually involves visually inspecting the arrival time of the signal to determine the arrival times of P and S waves of seismic waves.

[0067] Automatic first arrival picking: Automatic first arrival picking uses computer algorithms to automatically identify and extract the arrival times of microseismic events. This can be achieved by applying signal processing and machine learning techniques. Automatic first arrival picking methods can improve efficiency and quickly extract first arrivals in large amounts of data.

[0068] The goal of first arrival picking is to determine the P and S wave arrival times of microseismic events, thereby providing information about the time and location of the earthquake. This is of great significance for earthquake monitoring and seismological research.

[0069] According to the method for updating the seismic wave velocity model provided by the application, the pre-constructed seismic wave velocity model is updated according to the occurrence coordinates, the cutting time, the P wave arrival time and the S wave arrival time, and specifically comprises:

[0070] Input the occurrence coordinates into the velocity model to obtain the predicted travel time of the seismic wave;

[0071] According to the occurrence coordinates, the cutting time, the P wave arrival time and the S wave arrival time, the true travel time of the seismic wave is obtained;

[0072] According to the predicted travel time and the true travel time, the velocity model is updated.

[0073] In this embodiment, the method for updating the seismic wave velocity model can be carried out according to the following steps:

[0074] Input the occurrence coordinates into the velocity model: According to the occurrence coordinates, input them as input parameters into the pre-constructed seismic wave velocity model. This can be used to calculate the predicted travel time of the seismic wave at this location.

[0075] Obtaining true travel times from hypocenter coordinates, cutting times, P-wave arrival times, and S-wave arrival times: Using the data of hypocenter coordinates, cutting times, P-wave arrival times, and S-wave arrival times, the true travel times of seismic waves can be calculated. These data can be obtained based on seismic observation data, seismic network or other related measurement methods.

[0076] Updating the velocity model: By comparing the predicted travel times and the true travel times, the accuracy and errors of the velocity model can be evaluated. According to the comparison results, the velocity model can be updated. A common method is to use an inversion algorithm to adjust the parameters in the velocity model by minimizing the difference between the predicted travel times and the true travel times.

[0077] This updating method allows the seismic wave velocity model to be improved according to actual observation data, thus more accurately describing the velocity structure of the Earth's interior. This is of great significance for seismological research, earthquake monitoring and earthquake prediction, and helps to improve the accuracy of earthquake location and seismic waveform simulation.

[0078] According to the updating method of the seismic wave velocity model provided by the present application, the velocity model is updated according to the predicted travel times and the true travel times, specifically including:

[0079] According to the predicted travel times and the true travel times, the velocity model is updated by Levenberg-Marquardt algorithm.

[0080] In this embodiment, the updating method of the seismic wave velocity model can use Levenberg-Marquardt algorithm to update the velocity model. Levenberg-Marquardt algorithm is a nonlinear least squares optimization algorithm commonly used for parameter estimation and model fitting. The specific updating steps are as follows:

[0081] First, according to the difference between the predicted travel times and the true travel times, a target function is constructed, which is used to measure the fitting degree of the velocity model. Then, Levenberg-Marquardt algorithm is used to minimize the target function. This algorithm will search for the optimal solution in the parameter space of the velocity model to make the target function reach the minimum value.

[0082] In the iterative process of Levenberg-Marquardt algorithm, the parameters of the velocity model will be updated according to the difference between the predicted travel times and the true travel times. Through iterative optimization, the velocity model can be gradually adjusted to better fit the actual observation data.

[0083] The Levenberg-Marquardt algorithm has the advantages of being able to handle nonlinear problems and having fast convergence characteristics. It is widely used in the updating and inversion of seismic wave velocity models to improve the description of the velocity structure of the Earth's interior. Using the Levenberg-Marquardt algorithm to update the velocity model can improve the accuracy of seismic wave propagation simulation and help improve the results of applications such as earthquake location, seismic waveform simulation, and earthquake prediction.

[0084] According to the updating method of the seismic wave velocity model provided by the application, the microseismic events collected by the geophone module are obtained, specifically including:

[0085] According to the preset time interval, the microseismic events collected by the geophone module are obtained.

[0086] In this embodiment, the microseismic events collected by the geophone module are obtained according to the preset time interval.

[0087] Specifically, the time interval is determined: according to the needs and actual conditions, a suitable time interval is set to represent the frequency at which microseismic events are desired to be obtained. The time interval can be determined according to the research needs, monitoring purposes, or requirements of the data acquisition system.

[0088] Setting up the data acquisition system: configuring the geophone module to collect data at the preset time interval. This can be done through corresponding software or hardware settings. Ensure that the geophone module collects data within each time interval.

[0089] Obtaining microseismic event data: after the end of the preset time interval, the collected microseismic event data from the geophone module is obtained. These data are usually stored in digital form and can include the amplitude, timestamp, and other related parameters of the microseismic signal.

[0090] According to the updating method of the seismic wave velocity model provided by the application, after predicting the travel time of the source to the geophone module through the updated seismic wave velocity model, the method further includes:

[0091] According to the occurrence coordinates, the cutting time, the P-wave arrival time, and the S-wave arrival time, pre-stack migration imaging is performed to obtain the migration imaging result map of the microseismic event.

[0092] In this embodiment, after predicting the travel time of the source to the geophone module, pre-stack migration imaging can be performed to obtain the migration imaging result map of the microseismic event. The specific steps are as follows:

[0093] Preparing pre-stack migration data: using the predicted travel time data of the source to the geophone module as input, and combining the occurrence coordinates, cutting time, P-wave arrival time, and S-wave arrival time, etc. information, preparing the data for pre-stack migration.

[0094] Performing pre-stack migration imaging: according to the pre-stack migration algorithm, the prepared data is processed to generate the migration imaging result map of the microseismic event. Pre-stack migration is a commonly used seismic imaging technique, which restores the underground seismic reflection interface and body wave scattering by propagating and stacking the data in the seismic wave velocity model.

[0095] Analyzing migration imaging results: by analyzing the migration imaging result map, information about the underground structure can be obtained, such as the position, shape and reflection intensity of the underground reflection interface. This helps geological interpretation, underground structure analysis and exploration target identification, etc.

[0096] By performing pre-stack migration imaging, the data of microseismic events can be converted into images of underground structures, providing more intuitive and detailed geological information. This can help seismologists and exploration geophysicists better understand underground structures, further improve seismic wave velocity models, and promote the development of seismic monitoring and underground resource exploration.

[0097] Reference Figure 3 The following describes the seismic wave velocity model updating device provided by the present application. The seismic wave velocity model updating device described below can be referred to in conjunction with the seismic wave velocity model updating method described above.

[0098] The present application also provides a seismic wave velocity model updating device, comprising:

[0099] The acquisition unit 310 is configured to acquire microseismic events collected by the geophone module, wherein the microseismic events include occurrence coordinates of the microseismic events and cutting times of the microseismic events.

[0100] The calculation unit 320 is configured to acquire P-wave arrival times and S-wave arrival times of the microseismic events according to the microseismic events.

[0101] The updating unit 330 is configured to update a pre-constructed seismic wave velocity model according to the occurrence coordinates, the cutting times, the P-wave arrival times and the S-wave arrival times.

[0102] The prediction unit 340 is configured to predict the travel time of the seismic wave by using the updated seismic wave velocity model.

[0103] Figure 4 An example of an electronic device is shown in the schematic diagram of the physical structure of the electronic device, as shown in Figure 4As shown, the electronic device can include a processor 410, a communications interface 420, a memory 430, and a communications bus 440, wherein the processor 410, the communications interface 420, and the memory 430 complete mutual communication through the communications bus 440. The processor 410 can invoke a logical instruction in the memory 430 to execute an updating method of a seismic wave velocity model, the method comprising:

[0104] Obtaining a microseismic event collected by a geophone module, the microseismic event comprising occurrence coordinates of the microseismic event and a cutting time of the microseismic event;

[0105] According to the microseismic event, obtaining P-wave arrival time and S-wave arrival time of the microseismic event;

[0106] According to the occurrence coordinates, the cutting time, the P-wave arrival time, and the S-wave arrival time, updating a pre-constructed seismic wave velocity model;

[0107] Predicting a travel time of a seismic wave through the updated seismic wave velocity model.

[0108] In addition, the logical instruction in the memory 430 described above can be implemented in the form of a software function unit and sold or used as an independent product, and can be stored in a computer-readable storage medium. Based on such understanding, the technical solutions of the present application essentially or the part that contributes to the prior art or part of the technical solutions can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a plurality of instructions to make a computer device (which can be a personal computer, a server, or a network device, etc.) execute all or part of the steps of the method described in various embodiments of the present application. The foregoing storage medium includes: a U disk, a mobile hard disk, a read-only memory (ROM, Read-Only Memory), a random access memory (RAM, Random Access Memory), a magnetic disk or an optical disk, and various media that can store program codes.

[0109] On the other hand, the present application also provides a computer program product, the computer program product comprising a computer program, the computer program being storable on a non-transitory computer-readable storage medium, and the computer program being executable by a processor, so that a computer can execute the updating method of the seismic wave velocity model provided by the above-mentioned methods, the method comprising:

[0110] Obtaining a microseismic event collected by a geophone module, the microseismic event comprising occurrence coordinates of the microseismic event and a cutting time of the microseismic event;

[0111] According to the microseismic event, P-wave arrival time and S-wave arrival time of the microseismic event are acquired;

[0112] According to the occurrence coordinates, the cutting time, the P-wave arrival time and the S-wave arrival time, a pre-constructed seismic wave velocity model is updated;

[0113] Travel time of a seismic wave is predicted through the updated seismic wave velocity model.

[0114] In another aspect, the application further provides a non-transitory computer readable storage medium, which stores a computer program, and the computer program is executed by a processor to implement the method for updating the seismic wave velocity model provided by the above method, and the method comprises:

[0115] Microseismic events collected by a receiver module are acquired, and the microseismic events comprise occurrence coordinates of the microseismic events and cutting time of the microseismic events;

[0116] According to the microseismic event, P-wave arrival time and S-wave arrival time of the microseismic event are acquired;

[0117] According to the occurrence coordinates, the cutting time, the P-wave arrival time and the S-wave arrival time, a pre-constructed seismic wave velocity model is updated;

[0118] Travel time of a seismic wave is predicted through the updated seismic wave velocity model.

[0119] The device embodiments described above are only schematic, wherein the units described as separate components may or may not be physically separate, and the components displayed as units may or may not be physical units, i.e., may be located in one place, or may be distributed on multiple network units. Part or all of the modules can be selected according to actual needs to achieve the purpose of the embodiment. Those skilled in the art can understand and implement without creative labor.

[0120] From the above description of the embodiments, those skilled in the art can clearly understand that the embodiments can be realized by means of software plus necessary universal hardware platforms, and of course can also be realized by hardware. Based on such understanding, the above technical solutions, essentially or in other words, the part that contributes to the prior art can be embodied in the form of a software product, which can be stored in a computer readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes a plurality of instructions to make a computer device (which can be a personal computer, a server, or a network device, etc.) execute the methods described in each embodiment or some parts of the embodiments.

[0121] It should be pointed out finally that the above embodiments are only used to illustrate the technical solutions of the present application, but not to limit the same; and although the present application has been described in detail with reference to the foregoing embodiments, it should be appreciated by those skilled in the art that the technical solutions recorded in the foregoing embodiments can be modified, or some technical features thereof can be replaced equivalently; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present application.

Claims

1. A method of updating a seismic wave velocity model, characterized by, The method comprises the following steps: acquiring microseismic events collected by a detector module, the microseismic events comprising occurrence coordinates of the microseismic events and cutting times of the microseismic events, the microseismic events being microseismic events generated by cutting of a coal mining machine; acquiring P-wave arrival times and S-wave arrival times of the microseismic events according to the microseismic events; updating a pre-constructed seismic wave velocity model according to the occurrence coordinates, the cutting times, the P-wave arrival times and the S-wave arrival times; predicting travel times of seismic waves through the updated seismic wave velocity model; the step of updating the pre-constructed seismic wave velocity model according to the occurrence coordinates, the cutting times, the P-wave arrival times and the S-wave arrival times specifically comprises the following steps: inputting the occurrence coordinates into the velocity model to obtain predicted travel times of seismic waves; obtaining true travel times of seismic waves according to the occurrence coordinates, the cutting times, the P-wave arrival times and the S-wave arrival times; updating the velocity model according to the predicted travel times and the true travel times.

2. The method of updating a seismic wave velocity model according to claim 1, characterized in that, the step of acquiring the P-wave arrival times and the S-wave arrival times of the microseismic events according to the occurrence coordinates and the cutting times specifically comprises the following steps: performing data preprocessing on the microseismic events, the preprocessing comprising denoising processing, filtering processing and gain correction processing; performing first arrival picking on the preprocessed microseismic events to obtain the P-wave arrival times and the S-wave arrival times of the microseismic events.

3. The method of updating a seismic wave velocity model according to claim 1, characterized by, the step of updating the velocity model according to the predicted travel times and the true travel times specifically comprises the following step: updating the velocity model through a Levenberg-Marquardt algorithm according to the predicted travel times and the true travel times.

4. The method of updating a seismic wave velocity model according to claim 1, characterized by, the step of acquiring the microseismic events collected by the detector module specifically comprises the following step: acquiring the microseismic events collected by the detector module at a preset time interval.

5. The method of updating a seismic wave velocity model according to claim 1, characterized by, after predicting travel times of seismic sources to the detector module through the updated seismic wave velocity model, the method further comprises the following step: performing pre-stack migration imaging according to the occurrence coordinates, the cutting times, the P-wave arrival times and the S-wave arrival times to obtain a migration imaging result map of the microseismic events.

6. An apparatus for updating a seismic wave velocity model, characterized by, The method comprises the following steps: an acquisition unit is configured to acquire microseismic events collected by a detector module, the microseismic events comprising occurrence coordinates of the microseismic events and cutting times of the microseismic events, the microseismic events being microseismic events generated by cutting of a coal mining machine; a calculation unit is configured to acquire P-wave arrival times and S-wave arrival times of the microseismic events according to the microseismic events; an updating unit is configured to update a pre-constructed seismic wave velocity model according to the occurrence coordinates, the cutting times, the P-wave arrival times and the S-wave arrival times; a prediction unit is configured to predict travel times of seismic waves through the updated seismic wave velocity model; the step of updating the pre-constructed seismic wave velocity model according to the occurrence coordinates, the cutting times, the P-wave arrival times and the S-wave arrival times specifically comprises the following steps: inputting the occurrence coordinates into the velocity model to obtain predicted travel times of seismic waves; obtaining true travel times of seismic waves according to the occurrence coordinates, the cutting times, the P-wave arrival times and the S-wave arrival times; The velocity model is updated based on the predicted travel times and the real travel times.

7. An electronic device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, The processor implements the method for updating the seismic wave velocity model as claimed in any one of claims 1 to 5 when executing the program.

8. A non-transitory computer-readable storage medium having stored thereon a computer program, characterized in that, The computer program implements the method for updating the seismic wave velocity model as claimed in any one of claims 1 to 5 when executed by the processor.

9. A computer program product comprising a computer program, characterized in that, The computer program implements the method for updating the seismic wave velocity model as claimed in any one of claims 1 to 5 when executed by the processor. The computer program implements the method for updating the seismic wave velocity model as claimed in any one of claims 1 to 5 when executed by the processor.

Images