Advanced geological data calibration method and system based on blast wave reflection method

Abstract

The present invention provides an advanced geological data calibration method and system based on a blast wave reflection method, and relates to the technical field of data calibration. The present invention obtains raw data and geological data to be calibrated based on the blast wave reflection method, constructs a blast wave velocity analysis model according to the geological data, generates a formation velocity correction factor through the seismic wave velocity analysis model and the raw data to be calibrated, obtains an environmental interference signal, converts the environmental interference signal from the time domain to the frequency domain through Fourier transform, extracts the frequency domain characteristics of the environmental interference signal, analyzes the frequency domain characteristics of the environmental interference signal to obtain an environmental signal correction coefficient, generates a geological data correction factor according to the velocity correction factor and the environmental signal correction coefficient, and calibrates the geological data through the geological data correction factor. The present invention takes into account complex geological factors and surrounding environmental interference to improve the accuracy of calibration.

Description

Technical Field

[0001] The invention relates to the technical field of data calibration, and in particular to an advanced geological data calibration method and system based on a blasting wave reflection method. Background Art

[0002] At present, the blast wave reflection method is an important means of advanced geological prediction in underground projects such as tunnels and mines. By analyzing the propagation and reflection information of the blast wave in the rock layer, the geological conditions ahead can be predicted. However, during the propagation process, the blast wave is affected by factors such as the heterogeneity, porosity, and water content of the geological medium, and its propagation speed and energy will change, resulting in a decrease in the accuracy of the seismic data. In addition, the reflection signal of the blast wave is affected by many factors, such as the first arrival time of the blast wave, the propagation speed, and the Q value factor. The accurate acquisition and calibration of these parameters are crucial to improving the accuracy of earthquake prediction. At present, the processing of seismic signals mostly relies on traditional velocity models and empirical formulas, and it is impossible to effectively combine multiple geological parameters for comprehensive calibration, resulting in insufficient accuracy and reliability of seismic data. Especially in complex geological environments, traditional methods cannot fully reflect the propagation characteristics of blast waves under different geological conditions, resulting in large errors in the prediction results.

[0003] In the prior art, the publication number CN117805896A discloses a method and system for processing data for advanced geological prediction using the blast wave reflection method, which obtains the original data of the blast wave reflection method, performs processing preparation and data setting on the original data; and uses an intelligent analysis model to process the original data after processing preparation and data setting to obtain processed data. The prior art uses an intelligent analysis model to replace manual steps such as bandpass filtering, first arrival picking, and Q analysis in the processing of advanced geological prediction data using the blast wave reflection method to improve the efficiency and accuracy of data processing. However, the prior art still has defects. The prior art completely relies on the original data, lacks investigation of specific geological conditions and calibration of the original data, and the accuracy and reliability of seismic data are limited. Especially in complex geological environments, the processing method that relies solely on the original data cannot accurately reflect the propagation characteristics of the blast wave under different geological conditions, and is prone to produce large prediction errors, affecting the accuracy of geological prediction and the scientific nature of decision-making.

[0004] The above information disclosed in this Background section is only for enhancement of understanding of the background of the present disclosure and therefore it may contain information that does not constitute the prior art that is already known to one of ordinary skill in the art. Summary of the invention

[0005] The object of the present invention is to provide a method and system for advanced geological data calibration based on blasting wave reflection method to solve the problems raised in the above background technology.

[0006] To achieve the above object, the present invention provides the following technical solutions:

[0007] A method for calibrating advanced geological data based on blast wave reflection method, the specific steps include:

[0008] Step 1: Based on the blast wave reflection method, the original data and geological data to be calibrated are obtained; the original data include the formation direct wave velocity and the blast wave first arrival time in the direct wave adjustment stage; the geological data include the formation geological thickness, elastic modulus, bulk modulus, shear modulus and density;

[0009] Step 2: Construct a blast wave velocity analysis model based on geological data to obtain the propagation velocity of the blast wave in the formation, and obtain the propagation time of the direct wave based on the propagation velocity of the blast wave in the formation;

[0010] Step 3: The difference between the first arrival time of the blast wave and the propagation time of the direct wave is used as the calculation time error, and the calculation time error, the propagation velocity of the blast wave in the formation and the direct wave velocity of the formation are analyzed to obtain the formation velocity adjustment amount, and the formation velocity correction factor is obtained according to the formation velocity adjustment amount;

[0011] Step 4: Obtain environmental interference signals, convert the environmental interference signals from the time domain to the frequency domain through Fourier transform, extract the frequency domain characteristics of the environmental interference signals, and analyze the frequency domain characteristics of the environmental interference signals to obtain environmental signal correction coefficients, wherein the environmental interference signals include vibration signals generated by drilling and noise signals generated by blasting;

[0012] Step 5: Generate a geological data correction factor according to the velocity correction factor and the environmental signal correction coefficient, and calibrate the geological data using the geological data correction factor.

[0013] Furthermore, the specific logic for obtaining the propagation velocity of the blast wave in the formation is: obtaining geological data of the area to be calibrated, constructing a blast wave velocity analysis model based on the geological data, and obtaining the propagation velocity of the blast wave in the formation; the specific formula for obtaining the propagation velocity of the blast wave in the formation is:

[0014]

[0015] Among them, v i is the propagation speed of the explosion wave in the i-th layer, E i is the elastic modulus of the i-th layer, G i is the shear modulus of the i-th layer, K i is the bulk modulus of the i-th layer, ρ i is the density of the i-th layer, and i is the index of the stratum in the area to be analyzed.

[0016] Furthermore, the specific logic for obtaining the propagation time of the direct wave is: the propagation time of the direct wave in the blast wave velocity analysis model is calculated using the ray tracing method, and the specific formula is:

[0017]

[0018] Among them, T S is the propagation time of the direct wave, v i is the propagation speed of the explosion wave in the i-th layer, d i is the geological thickness of the i-th layer, N is the number of stratigraphic layers, and i is the stratigraphic layer index.

[0019] Furthermore, the specific logic for obtaining the velocity adjustment amount is: subtracting the first arrival time of the blast wave from the propagation time of the direct wave to obtain the calculated time error, and analyzing the calculated time error, the propagation velocity of the blast wave in the formation and the direct wave velocity of the formation to obtain the velocity adjustment amount; the specific formula for obtaining the calculated time error is:

[0020]

[0021] in, To calculate the time error, is the first arrival time of the blast wave, T S is the propagation time of the direct wave;

[0022] The specific formula for obtaining the speed adjustment is:

[0023]

[0024] in, is the speed adjustment of the i-th layer, v i is the propagation speed of the explosion wave in the i-th layer, d i is the geological thickness of the i-th layer, V i is the direct wave velocity of the i-th layer, N is the number of layers, and i is the index of the layer number.

[0025] Furthermore, the specific logic for obtaining the formation velocity correction factor is: obtaining the formation velocity correction factor according to the formation velocity adjustment amount, and the specific formula for obtaining the formation velocity correction factor is:

[0026]

[0027] Where Sp is the formation velocity correction factor.

[0028] Furthermore, the vibration signal generated by drilling and the noise signal generated by blasting are converted from the time domain to the frequency domain through Fourier transform, the frequency domain characteristics of the vibration signal generated by drilling and the noise signal generated by blasting are extracted, and the frequency domain characteristics of the vibration signal generated by drilling and the noise signal generated by blasting are analyzed to obtain the environmental signal correction coefficient, wherein the frequency domain characteristics include the center frequency and half-height width of the spectrum of the vibration signal generated by drilling, and the power spectrum density, spectrum entropy and main frequency of the noise signal generated by blasting; the specific logic for generating the environmental signal correction coefficient is:

[0029]

[0030] Where Ae is the environmental signal correction coefficient, Ae is the environmental signal correction coefficient, f 0 The center frequency of the vibration signal generated by drilling, P is the half-height width of the spectrum of the vibration signal generated by drilling; 0 is the power spectral density of the noise signal generated by the explosion, f 1 is the main frequency of the noise signal generated by the explosion, SA is the spectrum entropy of the noise signal generated by the explosion;

[0031] Furthermore, the specific logic for calibrating the geological data is as follows: a geological data correction factor is generated according to the velocity correction factor and the environmental signal correction coefficient, and the geological data is calibrated by the geological data correction factor; the specific formula for obtaining the geological data correction factor is as follows:

[0032] Dx=Sp*Ae+Ae+tanAe

[0033] Where Dx is the geological data correction factor, Sp is the formation velocity correction factor, and Ae is the environmental signal correction coefficient;

[0034] The specific logic for calibrating geological data through geological data correction factors is as follows:

[0035] x′=x(1-Dx%)

[0036] Among them, x′ is the corrected geological data, and x is the original geological data.

[0037] The present invention further provides an advanced geological data calibration system based on the blast wave reflection method, the system is used to implement the advanced geological data calibration method based on the blast wave reflection method, specifically comprising:

[0038] A data acquisition module, used to acquire raw data and geological data to be calibrated based on the blast wave reflection method; the raw data includes the formation direct wave velocity and the blast wave first arrival time in the direct wave adjustment stage; the geological data includes the formation geological thickness, elastic modulus, bulk modulus, shear modulus and density;

[0039] The speed-time analysis module is used to construct a blast wave velocity analysis model based on geological data, obtain the propagation velocity of the blast wave in the formation, and obtain the propagation time of the direct wave based on the propagation velocity of the blast wave in the formation;

[0040] The velocity calibration module is used to calculate the time error by taking the difference between the first arrival time of the blast wave and the propagation time of the direct wave, analyze the calculated time error, the propagation velocity of the blast wave in the formation and the direct wave velocity of the formation to obtain the formation velocity adjustment amount, and obtain the formation velocity correction factor according to the formation velocity adjustment amount;

[0041] A signal analysis module is used to obtain environmental interference signals, convert the environmental interference signals from the time domain to the frequency domain through Fourier transform, extract the frequency domain characteristics of the environmental interference signals, and analyze the frequency domain characteristics of the environmental interference signals to obtain environmental signal correction coefficients, wherein the environmental interference signals include vibration signals generated by drilling and noise signals generated by blasting;

[0042] The data calibration module is used to generate a geological data correction factor according to a velocity correction factor and an environmental signal correction coefficient, and calibrate the geological data through the geological data correction factor.

[0043] Compared with the prior art, the present invention has the following beneficial effects:

[0044] The present invention establishes a blast wave velocity analysis model, accurately calculates the propagation velocity of the blast wave in the stratum, generates a stratum velocity correction factor through the velocity model, corrects the geological data, and reduces the error caused by geological heterogeneity.

[0045] The present invention also takes into account the interference of complex environmental factors on the blasting wave, and quantifies the interference through a formula to generate an environmental signal correction coefficient. The geological data is calibrated by the environmental signal correction coefficient, which can make the calibration result more accurate. BRIEF DESCRIPTION OF THE DRAWINGS

[0046] Figure 1 It is a schematic diagram of the overall method flow of the present invention.

[0047] Figure 2 It is a schematic diagram of the overall system structure of the present invention. DETAILED DESCRIPTION

[0048] In order to make the objectives, technical solutions and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with specific embodiments.

[0049] It should be noted that, unless otherwise defined, the technical terms or scientific terms used in the present invention should be understood by people with ordinary skills in the field to which the present invention belongs. The words "first", "second" and similar words used in the present invention do not indicate any order, quantity or importance, but are only used to distinguish different components. "Include" or "comprise" and similar words mean that the elements or objects appearing before the word include the elements or objects listed after the word and their equivalents, without excluding other elements or objects. "Connect" or "connected" and similar words are not limited to physical or mechanical connections, but may include electrical connections, whether direct or indirect. "Up", "down", "left", "right" and the like are only used to indicate relative positional relationships. When the absolute position of the described object changes, the relative positional relationship may also change accordingly.

[0050] Example:

[0051] See also Figure 1 , the present invention provides a technical solution:

[0052] A method for calibrating advanced geological data based on blast wave reflection method, the specific steps include:

[0053] Step 1: Based on the blast wave reflection method, the original data and geological data to be calibrated are obtained; the original data include the formation direct wave velocity and the blast wave first arrival time in the direct wave adjustment stage; the geological data include the formation geological thickness, elastic modulus, bulk modulus, shear modulus and density;

[0054] Step 2: Construct a blast wave velocity analysis model based on geological data to obtain the propagation velocity of the blast wave in the formation, and obtain the propagation time of the direct wave based on the propagation velocity of the blast wave in the formation;

[0055] The specific logic for obtaining the propagation velocity of the blast wave in the formation is: obtain the geological data of the area to be analyzed, build a blast wave velocity analysis model based on the geological data, and obtain the propagation velocity of the blast wave in the formation; the specific formula for obtaining the propagation velocity of the blast wave in the formation is:

[0056]

[0057] Among them, v i is the propagation speed of the explosion wave in the i-th layer, E i is the elastic modulus of the i-th layer, G i is the shear modulus of the i-th layer, K i is the bulk modulus of the i-th layer, ρ i is the density of the i-th layer. The propagation speed v of the blast wave in the i-th layer i The direct wave can be quantified from a theoretical level, providing an important basis for the adjustment of the direct wave.

[0058] The specific logic for obtaining the propagation time of the direct wave is: the propagation time of the direct wave in the blast wave velocity analysis model is calculated using the ray tracing method, and the specific formula is:

[0059]

[0060] Among them, T S is the propagation time of the direct wave, v i is the propagation speed of the blast wave at layer i, d i is the geological thickness of the i-th layer, N is the number of layers, and i is the index of the layer number. The propagation time of the direct wave is T S It is an important parameter to determine whether the direct wave needs to be adjusted. If the direct wave propagation time error is within the preset range, there is no need to adjust the direct wave velocity; if the direct wave propagation time error exceeds the preset range, the direct wave velocity needs to be adjusted until the direct wave propagation time error is within the preset range.

[0061] Step 3: The difference between the first arrival time of the blast wave and the propagation time of the direct wave is used as the calculation time error, and the calculation time error, the propagation velocity of the blast wave in the formation and the direct wave velocity of the formation are analyzed to obtain the formation velocity adjustment amount, and the formation velocity correction factor is obtained according to the formation velocity adjustment amount;

[0062] The specific logic for obtaining the velocity adjustment is: subtract the first arrival time of the blast wave from the propagation time of the direct wave to obtain the calculated time error, and analyze the calculated time error, the propagation velocity of the blast wave in the formation and the direct wave velocity in the formation to obtain the velocity adjustment; the specific formula for obtaining the calculated time error is:

[0063]

[0064] in, To calculate the time error, T M is the first arrival time of the blast wave, T 1 is the propagation time of the direct wave; calculate the time error It is used to measure the error between the constructed velocity model and the original data. The larger the value, the greater the degree of calibration required.

[0065] The specific formula for obtaining the speed adjustment is:

[0066]

[0067] in, is the speed adjustment of the i-th layer, v i is the propagation speed of the blast wave at layer i, d iis the geological thickness of the i-th layer, V i is the direct wave velocity of the i-th layer, N is the number of layers, and i is the index of the layer number.

[0068] The specific logic for obtaining the formation velocity correction factor is: obtaining the formation velocity correction factor according to the formation velocity adjustment amount, and the specific formula for obtaining the formation velocity correction factor is:

[0069]

[0070] Among them, Sp is the formation velocity correction factor. The formation velocity correction factor Sp reflects the influence of each layer velocity on the geological data. The larger its value, the greater the influence of each layer velocity on the geological data, and the greater the calibration required. The generation of this factor can provide an important basis for the calibration of geological data.

[0071] Step 4: Obtain environmental interference signals, convert the environmental interference signals from the time domain to the frequency domain through Fourier transform, extract the frequency domain characteristics of the environmental interference signals, and analyze the frequency domain characteristics of the environmental interference signals to obtain environmental signal correction coefficients, wherein the environmental interference signals include vibration signals generated by drilling and noise signals generated by blasting;

[0072] The vibration signal generated by drilling and the noise signal generated by blasting are converted from the time domain to the frequency domain by Fourier transform, and the frequency domain characteristics of the vibration signal generated by drilling and the noise signal generated by blasting are extracted. The frequency domain characteristics of the vibration signal generated by drilling and the noise signal generated by blasting are analyzed to obtain the environmental signal correction coefficient. The frequency domain characteristics include the center frequency and half-height width of the spectrum of the vibration signal generated by drilling, and the power spectrum density, spectrum entropy and main frequency of the noise signal generated by blasting. The specific logic for generating the environmental signal correction coefficient is:

[0073]

[0074] Where Ae is the environmental signal correction coefficient, Ae is the environmental signal correction coefficient, f 0 The center frequency of the vibration signal generated by drilling, P is the half-height width of the spectrum of the vibration signal generated by drilling; 0 is the power spectral density of the noise signal generated by the explosion, f 1 is the main frequency of the noise signal generated by blasting, SA is the spectrum entropy of the noise signal generated by blasting; the environmental signal correction coefficient Ae reflects the degree of interference of vibration and noise interference in the environment on the accuracy of geological data. The larger its value, the greater the interference, the lower the accuracy of geological data, and the greater the calibration required. The generation of this coefficient can provide an important basis for the calibration of geological data; It reflects the comprehensive interference degree of the vibration signal generated by drilling on the accuracy of geological data. The larger the value, the greater the interference degree of the vibration signal generated by drilling on the geological data. 0 *f 1 *SA reflects the comprehensive interference degree of the noise signal generated by drilling and blasting on the accuracy of geological data. The larger the value, the greater the interference degree of the vibration signal generated by drilling and the noise signal generated by blasting on the geological data;

[0075] Step 5: Generate a geological data correction factor according to the velocity correction factor and the environmental signal correction coefficient, and calibrate the geological data using the geological data correction factor.

[0076] The specific logic for calibrating geological data is as follows: the geological data correction factor is generated according to the velocity correction factor and the environmental signal correction coefficient, and the geological data is calibrated by the geological data correction factor; the specific formula for obtaining the geological data correction factor is as follows:

[0077] Dx=Sp*Ae+Ae+tanAe

[0078] Among them, Dx is the geological data correction factor, Sp is the formation velocity correction factor, and Ae is the environmental signal correction coefficient; the geological data correction factor Dx reflects the degree to which the geological data needs to be corrected. The larger its value, the greater the amplitude of the geological data calibration.

[0079] The specific logic for calibrating geological data through geological data correction factors is as follows:

[0080] x′=x(1-Dx%)

[0081] Among them, x′ is the corrected geological data, and x is the original geological data.

[0082] See also Figure 2 The present invention further provides an advanced geological data calibration system based on the blast wave reflection method, the system is used to implement the advanced geological data calibration method based on the blast wave reflection method, specifically comprising:

[0083] A data acquisition module, used to acquire raw data and geological data to be calibrated based on the blast wave reflection method; the raw data includes the formation direct wave velocity and the blast wave first arrival time in the direct wave adjustment stage; the geological data includes the formation geological thickness, elastic modulus, bulk modulus, shear modulus and density;

[0084] The speed-time analysis module is used to construct a blast wave velocity analysis model based on geological data, obtain the propagation velocity of the blast wave in the formation, and obtain the propagation time of the direct wave based on the propagation velocity of the blast wave in the formation;

[0085] The velocity calibration module is used to calculate the time error by taking the difference between the first arrival time of the blast wave and the propagation time of the direct wave, analyze the calculated time error, the propagation velocity of the blast wave in the formation and the direct wave velocity of the formation to obtain the formation velocity adjustment amount, and obtain the formation velocity correction factor according to the formation velocity adjustment amount;

[0086] A signal analysis module is used to obtain environmental interference signals, convert the environmental interference signals from the time domain to the frequency domain through Fourier transform, extract the frequency domain characteristics of the environmental interference signals, and analyze the frequency domain characteristics of the environmental interference signals to obtain environmental signal correction coefficients, wherein the environmental interference signals include vibration signals generated by drilling and noise signals generated by blasting;

[0087] The data calibration module is used to generate a geological data correction factor according to a velocity correction factor and an environmental signal correction coefficient, and calibrate the geological data through the geological data correction factor.

[0088] The above formulas are all dimensionless and numerical calculations. The formula is a formula for the most recent real situation obtained by collecting a large amount of data and performing software simulation. The preset parameters in the formula are set by technicians in this field according to actual conditions.

[0089] The above embodiments may be implemented in whole or in part by software, hardware, firmware or any other combination thereof. When implemented using software, the above embodiments may be implemented in whole or in part in the form of a computer program product. Those skilled in the art may appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein may be implemented by electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are performed by hardware or software methods depends on the specific application and design constraints of the technical solution.

[0090] The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, and may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0091] The above description is only a specific implementation manner of the present application, but the protection scope of the present application is not limited thereto. Any technician familiar with the technical field can easily think of changes or substitutions within the technical scope disclosed in the present application, which should be included in the protection scope of the present application.

Claims

1. A method for advanced geological data calibration based on blast wave reflection method, characterized in that: The specific steps include: Step 1: Based on the blast wave reflection method, the original data and geological data to be calibrated are obtained; the original data include the formation direct wave velocity and the blast wave first arrival time in the direct wave adjustment stage; the geological data include the formation geological thickness, elastic modulus, bulk modulus, shear modulus and density; Step 2: Construct a blast wave velocity analysis model based on geological data to obtain the propagation velocity of the blast wave in the formation, and obtain the propagation time of the direct wave based on the propagation velocity of the blast wave in the formation; Step 3: The difference between the first arrival time of the blast wave and the propagation time of the direct wave is used as the calculation time error, and the calculation time error, the propagation velocity of the blast wave in the formation and the direct wave velocity of the formation are analyzed to obtain the formation velocity adjustment amount, and the formation velocity correction factor is obtained according to the formation velocity adjustment amount; Step 4: Obtain environmental interference signals, convert the environmental interference signals from the time domain to the frequency domain through Fourier transform, extract the frequency domain characteristics of the environmental interference signals, and analyze the frequency domain characteristics of the environmental interference signals to obtain environmental signal correction coefficients, wherein the environmental interference signals include vibration signals generated by drilling and noise signals generated by blasting; Step 5: Generate a geological data correction factor according to the velocity correction factor and the environmental signal correction coefficient, and calibrate the geological data using the geological data correction factor.

2. The method for advanced geological data calibration based on blast wave reflection method according to claim 1, characterized in that: The specific logic for obtaining the propagation velocity of the blast wave in the formation is: obtain the geological data of the area to be calibrated, build a blast wave velocity analysis model based on the geological data, and obtain the propagation velocity of the blast wave in the formation; the specific formula for obtaining the propagation velocity of the blast wave in the formation is: Among them, v i is the propagation speed of the explosion wave in the i-th layer, E i is the elastic modulus of the i-th layer, G i is the shear modulus of the i-th layer, K i is the bulk modulus of the i-th layer, ρ i is the density of the i-th layer, and i is the index of the stratum in the area to be analyzed.

3. The method for advanced geological data calibration based on blast wave reflection method according to claim 1, characterized in that: The specific logic for obtaining the propagation time of the direct wave is: the propagation time of the direct wave in the blast wave velocity analysis model is calculated using the ray tracing method, and the specific formula is: Among them, T S is the propagation time of the direct wave, v i is the propagation speed of the explosion wave in the i-th layer, d i is the geological thickness of the i-th layer, N is the number of stratigraphic layers, and i is the stratigraphic layer index.

4. The method for advanced geological data calibration based on blast wave reflection method according to claim 1, characterized in that: The specific logic for obtaining the velocity adjustment is: subtract the first arrival time of the blast wave from the propagation time of the direct wave to obtain the calculated time error, and analyze the calculated time error, the propagation velocity of the blast wave in the formation and the direct wave velocity in the formation to obtain the velocity adjustment; the specific formula for obtaining the calculated time error is: in, To calculate the time error, T M is the first arrival time of the blast wave, T S is the propagation time of the direct wave; The specific formula for obtaining the speed adjustment is: in, is the speed adjustment of the i-th layer, v i is the propagation speed of the explosion wave in the i-th layer, d i is the geological thickness of the i-th layer, V i is the direct wave velocity of the i-th layer, N is the number of layers, and i is the index of the layer number.

5. The method for advanced geological data calibration based on blast wave reflection method according to claim 4 is characterized in that: The specific logic for obtaining the formation velocity correction factor is: obtaining the formation velocity correction factor according to the formation velocity adjustment amount, and the specific formula for obtaining the formation velocity correction factor is: Where Sp is the formation velocity correction factor.

6. The method for advanced geological data calibration based on blast wave reflection method according to claim 1, characterized in that: The specific logic for obtaining the environmental signal correction coefficient is as follows: the vibration signal generated by drilling and the noise signal generated by blasting are converted from the time domain to the frequency domain through Fourier transform, the frequency domain characteristics of the vibration signal generated by drilling and the noise signal generated by blasting are extracted, and the frequency domain characteristics of the vibration signal generated by drilling and the noise signal generated by blasting are analyzed to obtain the environmental signal correction coefficient, wherein the frequency domain characteristics include the center frequency and half-height width of the spectrum of the vibration signal generated by drilling, and the power spectrum density, spectrum entropy and main frequency of the noise signal generated by blasting; the specific logic for generating the environmental signal correction coefficient is as follows: Among them, Ae is the environmental signal correction coefficient, f0 is the center frequency of the vibration signal generated by drilling, is the half-height width of the spectrum of the vibration signal generated by drilling; P0 is the power spectral density of the noise signal generated by blasting, f1 is the main frequency of the noise signal generated by blasting, and SA is the spectral entropy of the noise signal generated by blasting.

7. The method for advanced geological data calibration based on blast wave reflection method according to claim 1, characterized in that: The specific logic for calibrating geological data is as follows: the geological data correction factor is generated according to the velocity correction factor and the environmental signal correction coefficient, and the geological data is calibrated by the geological data correction factor; the specific formula for obtaining the geological data correction factor is as follows: Dx=Sp*Ae+Ae+tanAe Where Dx is the geological data correction factor, Sp is the formation velocity correction factor, and Ae is the environmental signal correction coefficient; The specific logic for calibrating geological data through geological data correction factors is as follows: x′=x(1-Dx%) Among them, x′ is the corrected geological data, and x is the original geological data.

8. An advanced geological data calibration system based on seismic wave reflection method, characterized in that: The system is used to implement the advanced geological data calibration method based on seismic wave reflection method as described in any one of claims 1 to 7, specifically comprising: A data acquisition module, used to acquire raw data and geological data to be calibrated based on the blast wave reflection method; the raw data includes the formation direct wave velocity and the blast wave first arrival time in the direct wave adjustment stage; the geological data includes the formation geological thickness, elastic modulus, bulk modulus, shear modulus and density; The speed-time analysis module is used to construct a blast wave velocity analysis model based on geological data, obtain the propagation velocity of the blast wave in the formation, and obtain the propagation time of the direct wave based on the propagation velocity of the blast wave in the formation; The velocity calibration module is used to calculate the time error by taking the difference between the first arrival time of the blast wave and the propagation time of the direct wave, analyze the calculated time error, the propagation velocity of the blast wave in the formation and the direct wave velocity of the formation to obtain the formation velocity adjustment amount, and obtain the formation velocity correction factor according to the formation velocity adjustment amount; A signal analysis module is used to obtain environmental interference signals, convert the environmental interference signals from the time domain to the frequency domain through Fourier transform, extract the frequency domain characteristics of the environmental interference signals, and analyze the frequency domain characteristics of the environmental interference signals to obtain environmental signal correction coefficients, wherein the environmental interference signals include vibration signals generated by drilling and noise signals generated by blasting; The data calibration module is used to generate a geological data correction factor according to a velocity correction factor and an environmental signal correction coefficient, and calibrate the geological data through the geological data correction factor.

Images