Blasting seismic wave simulation method, system, device and storage medium

Abstract

The present invention provides a method, system, device and storage medium for simulating blasting seismic waves, wherein the method comprises: obtaining the measured waveform of the blasting seismic wave, filtering the noise in the measured waveform, obtaining the zero point and extreme point of the waveform after filtering the noise, eliminating the redundant zero point, and using the waveform after eliminating the redundant as a sample waveform; fitting the half-wave frequency corresponding to each extreme point according to the zero point of the sample waveform, constructing the time history scatter points of the extreme point-half-wave frequency, bringing the time history scatter points into the frequency error function, solving the frequency attenuation parameter that makes the frequency error function take a minimum value, and fitting the frequency value of the wave based on the frequency attenuation parameter; bringing the extreme point into the intensity error function, solving the gamma intensity parameter that makes the intensity error function take a minimum value, and fitting the intensity value of the wave based on the gamma intensity parameter; bringing the frequency value and the intensity value into the simulation model to obtain the simulated blasting seismic wave. The present invention overcomes the inaccuracy of the traditional empirical formula by the method of theoretical derivation.

Description

Technical Field

[0001] The present invention relates to the field of waveform simulation technology, and in particular to a method, system, device and storage medium for simulating blasting seismic waves. Background Art

[0002] Seismic wave simulation is used to simulate the temporal or spatial variation of seismic waves, to correct collected waveform data, and to assist in assessing the impact of earthquakes on structures such as buildings, bridges, and roads.

[0003] In the related technologies, earthquake motion waves are often simulated, but there are few methods or programs for simulating blasting seismic waves. In the existing blasting seismic simulation technologies, most schemes directly rely on wave formulas and require a large number of on-site formulas to generate models. The above technical schemes have the following defects: lack of theoretical derivation process, resulting in a large difference between the simulated waves and the on-site waveforms; excessive calculation leads to low simulation efficiency.

[0004] Based on the above analysis of the development status of this technical field, the existing technology lacks a solution for extracting frequency-related parameters and intensity-related parameters by using theoretical derivation to solve extreme values ​​and zero points to simulate blasting seismic waves. Summary of the invention

[0005] The object of the present invention is to provide a method, system, device and storage medium for simulating blasting seismic waves, aiming to solve the above-mentioned problems in the prior art.

[0006] According to a first aspect of an embodiment of the present invention, there is provided a method for simulating blasting seismic waves, comprising:

[0007] Obtain the measured waveform of the blasting seismic wave, filter the noise in the measured waveform, obtain the zero point and extreme point of the waveform after the noise is filtered, eliminate the redundant zero points, and use the waveform after eliminating the redundant as the sample waveform;

[0008] According to the zero point of the sample waveform, the half-wave frequency corresponding to each extreme point is fitted, and the time-history scatter points of the extreme point-half-wave frequency are constructed. The time-history scatter points are brought into the frequency error function to solve the frequency attenuation parameter that makes the frequency error function take the minimum value, and the frequency value of the fitting wave is based on the frequency attenuation parameter;

[0009] Substitute the extreme value point into the intensity error function, solve the gamma intensity parameter that makes the intensity error function take the minimum value, and fit the intensity value of the wave based on the gamma intensity parameter;

[0010] The frequency and intensity values ​​are brought into the simulation model to obtain simulated blasting seismic waves.

[0011] According to a second aspect of an embodiment of the present invention, there is provided a blasting seismic wave simulation system, comprising:

[0012] The preprocessing module is used to obtain the measured waveform of the blasting seismic wave, filter the noise in the measured waveform, obtain the zero point and extreme point of the waveform after the noise is filtered, eliminate the redundant zero points, and use the waveform after eliminating the redundant as the sample waveform;

[0013] The frequency fitting module is used to fit the half-wave frequency corresponding to each extreme point according to the zero point of the sample waveform, construct the time-history scatter points of the extreme point-half-wave frequency, bring the time-history scatter points into the frequency error function, solve the frequency attenuation parameter that makes the frequency error function take the minimum value, and fit the frequency value of the wave based on the frequency attenuation parameter;

[0014] An intensity fitting module is used to bring the extreme value point into the intensity error function, solve the gamma intensity parameter that makes the intensity error function take the minimum value, and fit the intensity value of the wave based on the gamma intensity parameter;

[0015] The simulation module is used to bring the frequency value and intensity value into the simulation model to obtain the simulated blasting seismic wave.

[0016] According to a third aspect of an embodiment of the present invention, an electronic device is provided, comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the computer program is executed by the processor, the steps of the blasting seismic wave simulation method provided in the first aspect of the present disclosure are implemented.

[0017] According to a fourth aspect of an embodiment of the present invention, a computer-readable storage medium is provided, on which a program for implementing information transmission is stored. When the program is executed by a processor, the steps of the blasting seismic wave simulation method provided in the first aspect of the present disclosure are implemented.

[0018] The technical solution provided by the embodiment of the present invention includes the following beneficial effects: using an error function to represent the error between simulation and actual measurement, utilizing a method for solving the minimum value of the error function to obtain accurate coefficients of simulation frequency and intensity related formulas, and overcoming the inaccuracy of traditional empirical formulas by theoretical derivation methods; in the process of solving the minimum value, the number of zero points and extreme points is limited, which reduces the amount of calculation in the simulation process.

[0019] It is to be understood that the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the present disclosure. BRIEF DESCRIPTION OF THE DRAWINGS

[0020] In order to more clearly illustrate one or more embodiments of this specification or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments recorded in this specification. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative labor.

[0021] Figure 1 is a flow chart of a method for simulating blasting seismic waves according to an embodiment of the present invention;

[0022] Figure 2 is a schematic diagram of measured waveform loading according to an embodiment of the present invention;

[0023] Figure 3 is a schematic diagram of solving the frequency attenuation parameter of an embodiment of the present invention;

[0024] Figure 4 is a schematic diagram of frequency value fitting of an embodiment of the present invention;

[0025] Figure 5 is a schematic diagram of intensity value fitting of an embodiment of the present invention;

[0026] Figure 6 is a schematic diagram of blasting seismic wave simulation according to an embodiment of the present invention;

[0027] Figure 7 is a schematic diagram of a blasting seismic wave simulation system according to an embodiment of the present invention;

[0028] Figure 8 is a schematic diagram of an electronic device according to an embodiment of the present invention. DETAILED DESCRIPTION

[0029] In order to enable those skilled in the art to better understand the technical solutions in one or more embodiments of this specification, the following will be combined with the drawings in one or more embodiments of this specification to clearly and completely describe the technical solutions in one or more embodiments of this specification. Obviously, the described embodiments are only part of the embodiments of this specification, not all of the embodiments. Based on one or more embodiments of this specification, all other embodiments obtained by ordinary technicians in this field without creative work should fall within the scope of protection of this document.

[0030] Method Embodiment

[0031] According to an embodiment of the present invention, a method for simulating blasting seismic waves is provided. Figure 1 FIG. 1 is a flow chart of a method for simulating blasting seismic waves according to an embodiment of the present invention. Figure 1 As shown, the blasting seismic wave simulation method according to an embodiment of the present invention specifically includes:

[0032] In step S110, the measured waveform of the blasting seismic wave is obtained, the noise in the measured waveform is filtered, the zero points and extreme points of the waveform after the noise is filtered are obtained, the redundant zero points are eliminated, and the waveform after eliminating the redundant is used as a sample waveform, which specifically includes:

[0033] Export the measured wave data of the blasting seismic wave as a txt document, store the txt document in the default folder of matlab, compile the matlab program to load the waveform, and obtain the measured waveform. Figure 2 Schematic diagram of the measured waveform loading of an embodiment of the present invention. Figure 2 As shown in FIG. 1 , the process of fitting the waveform measured on site is demonstrated, and the moving average function is used to filter the noise in the measured waveform. The moving average function adopts exponentially weighted moving average.

[0034] The horizontal axis interval where the product of the waveform values ​​at adjacent moments is negative is taken as the zero value interval, that is, roughly find the possible position of the zero point, and the points where the product of the waveform values ​​at adjacent moments is negative are taken as a pair of candidate points. From each pair of candidate points, the point whose waveform value is closer to the zero value is selected as the zero point, that is, the point closer to the x-axis is selected as the zero point. After obtaining the coordinates of all zero points, it is easy to have zero points with very small intervals in this case, which may affect the parameter calculation; use the maximum value max function to obtain the extreme value point.

[0035] In order to remove the zero points with small time intervals, the attenuation characteristics of the seismic wave frequency are used to obtain the time interval of the first oscillation cycle in the waveform after filtering the noise, which is used as the constraint threshold. The oscillation cycle refers to the process in which the wave velocity decreases from zero value to zero value after reaching the extreme value over time. The time interval between each zero point and its adjacent zero point is calculated, and the zero points whose bilateral time intervals are both less than the constraint threshold are regarded as redundant zero points. The redundant zero points are removed from the waveform after filtering the noise, that is, the zero points that meet the conditions are accurately found; no other processing is done on the extreme value points to obtain the sample waveform;

[0036] In step S120, the half-wave frequencies corresponding to the extreme points are fitted according to the zero point of the sample waveform, and the time-history scatter points of the extreme point-half-wave frequency are constructed. The time-history scatter points are brought into the frequency error function, and the frequency attenuation parameter that makes the frequency error function take the minimum value is solved. The frequency value of the wave is fitted based on the frequency attenuation parameter, which specifically includes:

[0037] For the sine wave sinx, the period can be determined by determining the zero point. The time of one period is T, so the frequency of the sine wave is 2π / T, and therefore the half-wave frequency is 2π / T;

[0038] The specific calculation of the half-wave frequency is as follows: calculate the difference between the horizontal coordinates of the right adjacent zero point and the zero point as the calculated difference, calculate the division result of π and the calculated difference, that is, π / calculated difference, and use the division result as the half-wave frequency corresponding to the extreme point in the interval between the zero point and the right adjacent zero point;

[0039] Get the time value corresponding to the horizontal coordinate of the extreme point, and use the time value and the corresponding half-wave frequency value as the coordinate node pair of the time history scatter point, where the time value is the horizontal coordinate ti of the coordinate node pair, and the half-wave frequency value is the vertical coordinate f of the coordinate node pairi , the time-course scatter points are in the form of (t1, f1), (t2, f2)…, (tn, fn n ), the time scatter points can be reflected in the coordinate diagram to form a time course scatter diagram.

[0040] Substitute the horizontal and vertical coordinate values ​​of the time history scatter points into the frequency error function Q of formula 1:

[0041] Q=(f1-ηe -γt1 ) 2 +(f2-ηe -γt2 ) 2 +(f n -ηe -γtn ) 2 Formula 1;

[0042] Among them, f i represents the half-wave frequency, ti represents the time value corresponding to the extreme point, i represents the subscript of the time series, n represents the total number of time series of scattered points, ηe -γti represents the frequency attenuation formula, η and γ represent unknown frequency attenuation parameters, where η represents the initial angular frequency, γ represents the decline exponent, and e represents the base of the natural logarithm;

[0043] From the form of formula 1, we can understand that each item of Q is the error between the measured value at the moment and the standard theoretical value. When Q takes the minimum value, the overall error value is the smallest, and the time-history frequency to be simulated is closest to the measured curve. There is a minimum value on the right side of formula 1. When the coordinate value is substituted into formula 1, only the values ​​of η and γ on the right side of the equation are unknown. Adjusting the values ​​of η and γ can change the size of the Q value.

[0044] Using the method of solving binary equations, let the value of the first-order partial derivative of formula 1 with respect to η be 0, let the value of the first-order partial derivative of formula 1 with respect to γ ​​be 0, and solve the two partial derivative equations to obtain η and γ that minimize the frequency error function Q, and then convert ηe -γti The value of is taken as the frequency value of the seismic wave at time i. The specific process of solving is as follows. The first-order partial derivative of formula 1 with respect to η is shown in formula 3, and the first-order partial derivative of formula 1 with respect to γ ​​is shown in formula 4:

[0045]

[0046]

[0047] Considering the computing power of Matlab, the logarithm of each item in the brackets is taken when writing the program, which improves the computing efficiency without affecting the calculation results. Figure 3 Schematic diagram of solving the frequency attenuation parameter of an embodiment of the present invention. Figure 3As shown, the matlab code for solving the frequency attenuation parameters is shown;

[0048] After obtaining the frequency values ​​at all times, the time-frequency curve can be fitted. Figure 4 is a schematic diagram of frequency value fitting of an embodiment of the present invention, such as Figure 4 As shown, the comparison between the frequency value fitting and the measured waveform frequency is shown, and the fitting effect is good.

[0049] In step S130, the extreme value point is introduced into the intensity error function to solve the gamma intensity parameter that makes the intensity error function take the minimum value, and the intensity value of the wave is fitted based on the gamma intensity parameter, which specifically includes:

[0050] Substitute the horizontal and vertical coordinate values ​​of the extreme points in the sample waveform into the intensity error function Q2 in formula 2:

[0051]

[0052] Among them, Y i represents the value of the extreme point, ti represents the time value corresponding to the extreme point, i represents the subscript of the time series, N represents the total number of sample waveform time series, I0ti α-1 e -βti represents the gamma function, that is, the intensity envelope function in the field, α, β and I0 represent unknown gamma intensity parameters, a represents the gamma peak position parameter, β represents the gamma function shape position parameter, I0 represents the actual value of the gamma intensity, and e represents the base of the natural logarithm;

[0053] Let the first-order partial derivative of formula 2 with respect to a be 0, let the first-order partial derivative of formula 2 with respect to β be 0, let the first-order partial derivative of formula 2 with respect to I0 be 0, and solve the three partial derivative equations to obtain a, β, and I0 that minimize the intensity error function Q2. α-1 e -βti The value of is taken as the intensity value of the seismic wave at time i.

[0054] After obtaining the intensity values ​​at all times, the time-history intensity curve can be fitted. Figure 5 is a schematic diagram of the intensity value fitting of an embodiment of the present invention, such as Figure 5 As shown, the comparison between the intensity value fitting and the measured waveform intensity is shown, and the fitting effect is good.

[0055] In step S140, the frequency value and the intensity value are brought into the simulation model to obtain a simulated blasting seismic wave, which specifically includes:

[0056] Substituting the frequency value and the intensity value into the simulation model of Formula 5, the parameter expressions in Formula 5 are shown in Formula 6 and Formula 7:

[0057]

[0058]

[0059]

[0060] Among them, X(t) represents the amplitude of the simulated blasting seismic wave at time t, A(t) represents the intensity function of the model, which is used to control the intensity of the blasting seismic wave curve, x(t) represents the frequency function of the model, which is used to control the frequency of the blasting seismic wave curve, f(t) represents the impulse response curve of white noise, σ f represents the normalization curve, which is used to normalize the unnormalized f(t) to x(t), w(τ) represents Gaussian white noise, S represents the power spectrum density of white noise, and is mainly used to construct the normalization curve σ f , τ represents the relative time variable in the integration process, ξ represents the damping ratio, also known as wave impedance, ω(τ) represents the actual angular frequency, h represents the unit impact velocity response, and the time-varying angular frequency in h is given by ω(τ)=ηe -γτ , and convolve white noise, the unit pulse expression in the convolution operation is h[t-τ, λ(τ)], where λ(τ) represents the angular frequency expression relative to the time variable;

[0061] In the above formula, A(t) = I0t α-1 e -βt , or ω(τ)=ηe -γτ .

[0062] Export the final simulated blasting seismic wave, Figure 6 Schematic diagram of blasting seismic wave simulation according to an embodiment of the present invention. Figure 6 As shown, the effect of the final simulated seismic wave is shown. The embodiment of the present invention performs two simulations. The frequency value and intensity value solved by the present application are unique, but the white noise in the simulation model is random, so different results appear. Figure 6 Both sets of waveforms in may be waveforms under this geological condition. Preferably, the addresses corresponding to the seismic waves can be derived while deriving the blasting seismic waves.

[0063] In summary, in response to the existing problems, the blasting seismic wave simulation method invented this time cleverly uses the negative product of both sides of the zero point to roughly find the zero point, and uses the frequency attenuation characteristics of seismic waves to remove redundant zero points to reduce possible errors in parameter estimation; the error function is used to represent the error between simulation and measurement, and the method of solving the minimum value of the error function is used to obtain accurate coefficients of the simulation frequency and intensity related formulas, and the inaccuracy of traditional empirical formulas is overcome by theoretical derivation. The frequency attenuation formula and the intensity envelope function have high theoretical depth, which ensures the effectiveness of blasting seismic wave simulation; in the process of solving the minimum value, the number of zero points and extreme points is limited, which makes the amount of calculation in the simulation process small.

[0064] System Example

[0065] According to an embodiment of the present invention, a blasting seismic wave simulation system is provided. Figure 7 Schematic diagram of a blasting seismic wave simulation system according to an embodiment of the present invention. Figure 7 As shown, the blasting seismic wave simulation system according to an embodiment of the present invention specifically includes:

[0066] The preprocessing module 70 is used to obtain the measured waveform of the blasting seismic wave, filter the noise in the measured waveform, obtain the zero point and extreme point of the waveform after the noise is filtered, eliminate the redundant zero points, and use the waveform after eliminating the redundant as the sample waveform, which is specifically used for:

[0067] Use the moving average function to filter out noise in the measured waveform.

[0068] The points where the product of the waveform values ​​at adjacent moments is negative are taken as a pair of candidate points, and the point whose waveform value is closer to zero is selected from each pair of candidate points as the zero point;

[0069] Use the maximum function to find the extreme points.

[0070] The time interval of the first oscillation cycle in the waveform after noise filtering is obtained as the constraint threshold, and the time interval between each zero point and its adjacent zero point is calculated. The zero points whose bilateral time intervals are both less than the constraint threshold are regarded as redundant zero points, and the redundant zero points are removed from the waveform after noise filtering.

[0071] The frequency fitting module 72 is used to fit the half-wave frequency corresponding to each extreme point according to the zero point of the sample waveform, construct the time history scatter points of the extreme point-half-wave frequency, bring the time history scatter points into the frequency error function, solve the frequency attenuation parameter that makes the frequency error function take the minimum value, and fit the frequency value of the wave based on the frequency attenuation parameter, which is specifically used for:

[0072] Calculate the difference between the horizontal coordinates of the right adjacent zero point and the zero point as the calculated difference, calculate the division result of π and the calculated difference, and use the division result as the half-wave frequency corresponding to the extreme point in the interval between the zero point and the right adjacent zero point;

[0073] Get the time value corresponding to the horizontal coordinate of the extreme point, and use the time value and the corresponding half-wave frequency value as the coordinate node pair of the time history scatter point, where the time value is the horizontal coordinate ti of the coordinate node pair, and the half-wave frequency value is the vertical coordinate f of the coordinate node pair i .

[0074] Substitute the horizontal and vertical coordinate values ​​of the time history scatter points into the frequency error function Q of formula 1:

[0075] Q=(f1-ηe -γt1 ) 2 +(f2-ηe -γt2 ) 2 +(f n -ηe -γtn ) 2 Formula 1;

[0076] Among them, f i represents the half-wave frequency, ti represents the time value corresponding to the extreme point, i represents the subscript of the time series, n represents the total number of time series of scattered points, ηe -γti represents the frequency attenuation formula, η and γ represent unknown frequency attenuation parameters, η represents the initial angular frequency, and γ represents the decline exponent;

[0077] Let the value of the first-order partial derivative of formula 1 with respect to η be 0, let the value of the first-order partial derivative of formula 1 with respect to γ ​​be 0, combine the two partial derivative equations, and solve to obtain η and γ that minimize the frequency error function Q, and then ηe -γti The value of is taken as the frequency value of the seismic wave at moment i.

[0078] The intensity fitting module 74 is used to bring the extreme value point into the intensity error function, solve the gamma intensity parameter that makes the intensity error function take the minimum value, and fit the intensity value of the wave based on the gamma intensity parameter, specifically for:

[0079] Substitute the horizontal and vertical coordinate values ​​of the extreme points in the sample waveform into the intensity error function Q2 in formula 2:

[0080]

[0081] Among them, Y i represents the value of the extreme point, ti represents the time value corresponding to the extreme point, i represents the subscript of the time series, N represents the total number of sample waveform time series, I0ti α-1 e -βti represents the gamma function, α, β and I0 represent unknown gamma intensity parameters, α represents the gamma peak position parameter, β represents the gamma function shape position parameter, and I0 represents the actual value of the gamma intensity;

[0082] Let the first-order partial derivative of formula 2 with respect to α be 0, let the first-order partial derivative of formula 2 with respect to β be 0, let the first-order partial derivative of formula 2 with respect to I0 be 0, combine the three partial derivative equations, and solve to obtain α, β and I0 that make the intensity error function Q2 take the minimum value, and then I0ti a-1 e -βti The value of is taken as the intensity value of the seismic wave at time i.

[0083] The simulation module 76 is used to bring the frequency value and the intensity value into the simulation model to obtain the simulated blasting seismic wave.

[0084] In summary, in response to the existing problems, this invention has invented a blasting seismic wave simulation system, which cleverly uses the negative product of both sides of the zero point to roughly find the zero point, and uses the frequency attenuation characteristics of seismic waves to remove redundant zero points to reduce possible errors in parameter estimation; uses an error function to represent the error between simulation and actual measurement, and uses the method of solving the minimum value of the error function to obtain accurate coefficients of the simulation frequency and intensity related formulas, and overcomes the inaccuracy of traditional empirical formulas by theoretical deduction. The frequency attenuation formula and the intensity envelope function have a high theoretical depth, which ensures the effectiveness of blasting seismic wave simulation; in the process of solving the minimum value, the number of zero points and extreme points is limited, which makes the amount of calculation in the simulation process small.

[0085] Electronic device embodiment

[0086] Figure 8 8 is a schematic diagram of an electronic device according to an embodiment of the present invention. The electronic device 800 may include at least one processor 810 and a memory 820. The processor 810 may execute instructions stored in the memory 820. The processor 810 is connected to the memory 820 through a data bus. In addition to the memory 820, the processor 810 may also be connected to an input device 830, an output device 840, and a communication device 850 through a data bus.

[0087] The processor 810 may be any conventional processor, such as a commercially available CPU. The processor may also include a graphics processor (Graphic Process Unit, GPU), a field programmable gate array (Field Programmable Gate Array, FPGA), a system on chip (System on Chip, SOC), an application specific integrated circuit (Application Specific Integrated Circuit, ASIC) or a combination thereof.

[0088] The memory 820 may be implemented by any type of volatile or nonvolatile storage device or a combination thereof, such as static random access memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic memory, flash memory, magnetic disk or optical disk.

[0089] In the embodiment of the present disclosure, executable instructions are stored in the memory 820, and the processor 810 can read the executable instructions from the memory 820 and execute the instructions to implement all or part of the steps of any blasting seismic wave simulation method in the above exemplary embodiments.

[0090] Computer Readable Storage Medium Embodiments

[0091] In addition to the above-mentioned methods and devices, an exemplary embodiment of the present disclosure may also be a computer program product or a computer-readable storage medium storing the computer program product. The computer product includes computer program instructions, and the computer program instructions can be executed by a processor to implement all or part of the steps described in any of the blasting seismic wave simulation methods in the above-mentioned exemplary embodiments.

[0092] The computer program product may be written in any combination of one or more programming languages ​​to write program codes for performing the operations of the embodiments of the present application, including object-oriented programming languages ​​such as Java, C++, etc., and also conventional procedural programming languages ​​such as "C" language or similar programming languages ​​and scripting languages ​​(e.g., Python). The program code may be executed entirely on the user computing device, partially on the user computing device, as an independent software package, partially on the user computing device and partially on a remote computing device, or entirely on a remote computing device or server.

[0093] Computer readable storage media can adopt any combination of one or more readable media. The readable medium can be a readable signal medium or a readable storage medium. The readable storage medium can include, for example, but is not limited to, a system, device or device of electricity, magnetism, light, electromagnetic, infrared, or semiconductor, or any combination of the above. More specific examples of readable storage media include: a static random access memory (SRAM) with one or more wires electrically connected, an electrically erasable programmable read-only memory (EEPROM), an erasable programmable read-only memory (EPROM), a programmable read-only memory (PROM), a read-only memory (ROM), a magnetic memory, a flash memory, a magnetic disk or an optical disk, or any suitable combination of the above.

[0094] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or replace some or all of the technical features therein with equivalents. However, these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for simulating blasting seismic waves, characterized in that: include: Obtaining a measured waveform of a blasting seismic wave, filtering noise in the measured waveform, obtaining zero points and extreme points of the waveform after the noise is filtered, eliminating redundant zero points, and using the waveform after eliminating the redundant points as a sample waveform; According to the zero point of the sample waveform, the half-wave frequency corresponding to each extreme point is fitted, and the time-history scatter points of the extreme point-half-wave frequency are constructed. The time-history scatter points are brought into the frequency error function to solve the frequency attenuation parameter that makes the frequency error function take the minimum value, and the frequency value of the fitting wave is based on the frequency attenuation parameter; Substituting the extreme value point into the intensity error function, solving the gamma intensity parameter that makes the intensity error function take a minimum value, and fitting the intensity value of the wave based on the gamma intensity parameter; The frequency value and the intensity value are brought into a simulation model to obtain a simulated blasting seismic wave.

2. The method according to claim 1, characterized in that The filtering of the noise in the measured waveform specifically includes: using a moving average function to filter the noise in the measured waveform.

3. The method according to claim 1, characterized in that The step of obtaining the zero point and the extreme point of the waveform after filtering the noise specifically includes: The points where the product of the waveform values ​​at adjacent moments is negative are taken as a pair of candidate points, and the point whose waveform value is closer to zero is selected from each pair of candidate points as the zero point; The extreme point is obtained using a maximum value function.

4. The method according to claim 1, characterized in that The redundant zero points removed specifically include: The time interval of the first oscillation cycle in the waveform after noise filtering is obtained as the constraint threshold, the time interval between each zero point and its adjacent zero point is calculated, and the zero points whose bilateral time intervals are both less than the constraint threshold are taken as redundant zero points, and the redundant zero points are removed from the waveform after noise filtering.

5. The method according to claim 1, characterized in that The step of fitting the half-wave frequencies corresponding to the extreme points according to the zero point of the sample waveform to construct the time-history scatter points of the extreme point-half-wave frequency specifically includes: Calculate the difference between the horizontal coordinates of the right adjacent zero point and the zero point as a calculated difference, calculate the division result of π and the calculated difference, and use the division result as the half-wave frequency corresponding to the extreme point in the interval between the zero point and the right adjacent zero point; Get the time value corresponding to the horizontal coordinate of the extreme point, and use the time value and the corresponding half-wave frequency value as the coordinate node pair of the time history scatter point, wherein the time value is the horizontal coordinate ti of the coordinate node pair, and the half-wave frequency value is the vertical coordinate fi of the coordinate node pair i .

6. The method according to claim 5, characterized in that The step of bringing the time history scatter points into the frequency error function, solving the frequency attenuation parameter that makes the frequency error function take a minimum value, and fitting the frequency value of the wave based on the frequency attenuation parameter specifically includes: Substitute the horizontal and vertical coordinate values ​​of the time history scatter points into the frequency error function Q of formula 1: Q = (f1 - ηe -γt1 ) 2 + (f2 - ηe -γt2 ) 2 + (f n - ηe -γtn ) 2 Formula 1; Among them, f i represents the half-wave frequency, ti represents the time value corresponding to the extreme point, i represents the subscript of the time series, n represents the total number of time series of scattered points, ηe -γti represents the frequency attenuation formula, η and γ represent unknown frequency attenuation parameters, η represents the initial angular frequency, and γ represents the decline exponent; Let the value of the first-order partial derivative of formula 1 with respect to η be 0, let the value of the first-order partial derivative of formula 1 with respect to γ ​​be 0, combine the two partial derivative equations, and solve to obtain η and γ that minimize the frequency error function Q, and then ηe -γti The value of is taken as the frequency value of the seismic wave at moment i.

7. The method according to claim 1, characterized in that The step of bringing the extreme value point into the intensity error function, solving the gamma intensity parameter that makes the intensity error function take a minimum value, and fitting the intensity value of the wave based on the gamma intensity parameter specifically includes: Substitute the horizontal and vertical coordinate values ​​of the extreme points in the sample waveform into the intensity error function Q2 in formula 2: Among them, Y i represents the value of the extreme point, ti represents the time value corresponding to the extreme point, i represents the subscript of the time series, N represents the total number of sample waveform time series, I0ti α-1 e -βti represents the gamma function, α, β and I0 represent unknown gamma intensity parameters, α represents the gamma peak position parameter, β represents the gamma function shape position parameter, and I0 represents the actual value of the gamma intensity; Let the first-order partial derivative of formula 2 with respect to α be 0, let the first-order partial derivative of formula 2 with respect to β be 0, let the first-order partial derivative of formula 2 with respect to I0 be 0, combine the three partial derivative equations, and solve to obtain α, β and I0 that make the intensity error function Q2 take the minimum value, and then I0ti α-1 e -βti The value of is taken as the intensity value of the seismic wave at time i.

8. A blasting seismic wave simulation device, characterized in that: include: A preprocessing module is used to obtain the measured waveform of the blasting seismic wave, filter the noise in the measured waveform, obtain the zero point and extreme point of the waveform after the noise is filtered, eliminate the redundant zero points, and use the waveform after eliminating the redundant as a sample waveform; A frequency fitting module is used to fit the half-wave frequency corresponding to each extreme point according to the zero point of the sample waveform, construct a time-history scatter point of the extreme point-half-wave frequency, bring the time-history scatter point into the frequency error function, solve the frequency attenuation parameter that makes the frequency error function take a minimum value, and fit the frequency value of the wave based on the frequency attenuation parameter; An intensity fitting module, used for bringing the extreme value point into the intensity error function, solving the gamma intensity parameter that makes the intensity error function take a minimum value, and fitting the intensity value of the wave based on the gamma intensity parameter; The simulation module is used to bring the frequency value and the intensity value into the simulation model to obtain simulated blasting seismic waves.

9. An electronic device, characterized in that: include: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the computer program, when executed by the processor, implements the steps of the blasting seismic wave simulation method according to any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores an implementation program for information transmission, and when the program is executed by a processor, the steps of the blasting seismic wave simulation method according to any one of claims 1 to 7 are implemented.

Images