Seismic wavelet analysis-based seismic excitation parameter optimization method and related equipment

Abstract

The invention discloses a seismic excitation parameter optimization method based on seismic wavelet analysis and related equipment, and the method mainly employs a nonlinear finite element explicit dynamics numerical simulation technology, constructs a numerical model according to rock physical parameters of an area adjacent to an explosion source, characteristics of a near-surface low-speed-reduction layer structure and a seismic excitation charging structure, and carries out the optimization of seismic excitation parameters. Energy and pressure propagation characteristics after explosion excitation are simulated, seismic wavelets excited by different geological conditions and excitation parameters are obtained, field construction seismic excitation parameters are optimized on the basis, rock mechanics physical parameters, seismic geological conditions and site construction parameters are fully utilized, and the method is suitable for mass production. Seismic wavelets are obtained through nonlinear finite element explicit dynamics numerical simulation, conditions are created for optimizing seismic excitation parameters in complex work area construction, and the whole method is short in research period and high in efficiency.

Description

Technical Field

[0001] The invention belongs to the field of oil exploration seismic data acquisition, and specifically relates to a seismic excitation parameter optimization method based on seismic wavelet analysis and related equipment. Background Art

[0002] The use of explosives to stimulate is the main method for onshore oil and gas seismic exploration, which has the advantages of wide vibration frequency band and easy energy adjustment. Practice has proved that the lithology of the stimulus, the type of explosives, the amount of explosives and the coupling factors all have an important influence on the amplitude and frequency characteristics of the signal in the seismic record.

[0003] In order to improve the signal-to-noise ratio and resolution of seismic data, the current method mainly analyzes the characteristics of the collected seismic recording signals through field tests, so as to optimize the excitation parameters. This field test method obtains accurate data, but has the characteristics of long research cycle, high cost, low efficiency, and difficulty in achieving a single factor, which urgently needs to be improved.

[0004] (1) Concept of seismic wavelet

[0005] A seismic wavelet is a signal with a certain start time, limited energy and a certain duration. It is the basic unit in seismic records. It is generally believed that the seismic wave generated by the explosion in the rock is a sharp pulse with a very short duration. As the sharp pulse propagates in the viscoelastic medium, the high-frequency component of the sharp pulse decays quickly, and the waveform grows accordingly, forming a seismic wavelet. A seismic wavelet generally has a duration of 2 to 3 phases, and then propagates underground in the form of a seismic wavelet.

[0006] Seismic wavelet is a very critical concept in seismic exploration. In forward modeling problems, it is necessary to combine the wave equation or convolution model with the seismic wavelet to form forward simulation seismic data. In inversion and deconvolution problems, it is also necessary to extract a wavelet from the seismic trace. Different wavelets often have different effects on the inversion results.

[0007] Seismic exploration practice shows that the sharper the seismic wavelet, the wider its frequency band and the higher the resolution of seismic exploration; the characteristics of seismic wavelets excited by different construction parameters are also different. Near-surface geological conditions, charge structure, rock physical parameters and the type of coupling medium all have a certain influence on the characteristics of wavelets.

[0008] The acquisition of seismic wavelets through on-site field seismic excitation construction has high construction costs and long cycles, which is not conducive to conducting seismic exploration and research on complex structures. It is necessary to make improvements.

[0009] (2) Relationship between particle vibration and seismic wavelet

[0010] The practice of oil and gas exploration research has proved that seismic detectors respond to the vertical and horizontal displacements of ground particles, converting mechanical vibrations into electrical signals to obtain seismic waves ( Figure 1 ).

[0011] Research shows that the seismic waves collected by the detector are related to velocity-type wavelets.

[0012] The detector motion equation is as follows:

[0013]

[0014] The detector output voltage equation is as follows:

[0015]

[0016] Among them, X: displacement of the coil relative to the ground, Z: displacement of ground particle vibration, v: vibration velocity of ground particle, u: output voltage of the detector, the value of the output voltage represents the strength change of the seismic signal.

[0017] Conventional moving coil geophones measure the rate of change or velocity of displacement motion of ground particles, while geophones such as DSU-3 measure the secondary derivative of displacement motion, i.e. acceleration. After transforming the geophone output voltage equation into the frequency domain, the output voltage spectrum of the moving coil geophone is proportional to the vibration velocity spectrum of the ground particles. This relationship is the transfer function H(jW). Its equation is as follows:

[0018] U(jω)=H(jω)υ(jω)

[0019] In this equation, U(jw) is the frequency spectrum of the detector output voltage, H(jw) is the input function, and v(jw) is the frequency spectrum of the particle vibration velocity. This shows that the seismic waves collected by the detector are related to the velocity-type sub-wave.

[0020] (3) Feasibility of obtaining seismic wavelets through numerical simulation

[0021] Using nonlinear finite element numerical simulation software, a seismic material strength model and a geometric structure numerical model are established. Through numerical simulation of explosion impact, the vibration velocity information of underground particles near the explosion point is obtained, and the basic input signal (i.e., seismic wavelet) of the seismic exploration excitation link is obtained.

[0022] The amplitude and frequency characteristics of the seismic wavelet obtained through numerical simulation are related to the physical parameters (density, Poisson's ratio, elastic modulus, etc.) of the surrounding rock and seismic wave propagation medium excited by explosives, and are also related to the field construction method (such as the amount of explosives, drilling depth, and coupling of the medium in the well). Therefore, conducting relevant research on the numerical simulation wavelet of seismic excitation has certain guiding significance for quickly determining the excitation parameters of field construction and optimizing the construction method. Summary of the invention

[0023] The present invention provides a method for optimizing seismic excitation parameters based on seismic wavelet analysis and related equipment to solve the problems of long field test research cycle, high cost, low efficiency, and difficulty in achieving a single factor.

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

[0025] The method for optimizing seismic excitation parameters based on seismic wavelet analysis includes:

[0026] Obtain the basic physical parameters of rocks and establish the constitutive model of rocks according to the basic physical parameters of rocks;

[0027] Obtaining explosive parameters, and establishing a state equation according to the explosive parameters;

[0028] Establish numerical models based on equations of state and constitutive models;

[0029] Based on the numerical model, different parameters are simulated to obtain the vibration velocity information of the observed particle;

[0030] Perform amplitude and spectrum analysis on the observed point vibration velocity information to obtain the basic characteristics of the seismic wavelet;

[0031] The construction method is optimized based on the basic characteristics of seismic wavelets.

[0032] Preferably, the basic physical parameters of the rock are obtained by coring the rock in the analysis area and then performing mechanical deformation testing.

[0033] Preferably, the rock constitutive model adopts the Johnson-Cook constitutive relationship model.

[0034] Preferably, the constitutive model of the rock is specifically:

[0035]

[0036] Among them, σ is the dynamic yield stress of the material; A is the quasi-static yield stress of the material; B is the strain hardening modulus of the material; ε p is the effective plastic strain; ε p * is the normalized effective plastic strain rate, ε is the plastic strain rate, T H is the homologous temperature, C, n and m are material constants.

[0037] Preferably, the state equation is:

[0038]

[0039] Where p is the pressure; E is the internal energy per unit volume of the detonation products; V is the relative volume of the detonation products, that is, the ratio of the volume of the detonation products after the explosion to their initial volume; A, B, R1, R2 and ω are constants.

[0040] Preferably, the numerical model is divided into finite element meshes.

[0041] Preferably, based on the basic characteristics of seismic wavelets, the preferred construction method is specifically:

[0042] The energy growth trend of the seismic wavelet is obtained according to the basic characteristics of the seismic wavelet, the inflection point of the seismic excitation energy is obtained according to the energy growth trend of the seismic wavelet, and the construction method is optimized according to the inflection point of the seismic excitation energy.

[0043] The seismic excitation parameter optimization system based on seismic wavelet analysis includes:

[0044] Constitutive model acquisition module: obtains the basic physical parameters of rocks and establishes the constitutive model of rocks according to the basic physical parameters of rocks;

[0045] State equation establishment module: obtain explosive parameters and establish the state equation according to the explosive parameters;

[0046] Numerical model building module: build numerical models based on state equations and constitutive models;

[0047] Vibration velocity acquisition module: simulates different parameters based on the numerical model to obtain the vibration velocity information of the observed particle;

[0048] Vibration velocity analysis module: performs amplitude and spectrum analysis on the observed vibration velocity information to obtain the basic characteristics of the seismic wavelet;

[0049] Optimal construction module: Optimal construction method based on the basic characteristics of seismic wavelets.

[0050] A computer device comprises a memory, a processor and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the steps of the method for optimizing seismic excitation parameters based on seismic wavelet analysis are implemented.

[0051] A computer-readable storage medium stores a computer program, which, when executed by a processor, implements the steps of the method for optimizing seismic excitation parameters based on seismic wavelet analysis.

[0052] Compared with the prior art, the present invention has the following beneficial effects: the present invention provides a method for optimizing seismic excitation parameters based on seismic wavelet analysis, which mainly utilizes nonlinear finite element explicit dynamics numerical simulation technology, and constructs a numerical model according to the rock physical parameters of the area adjacent to the explosion source, the structural characteristics of the near-surface low-velocity reduction layer, and the seismic excitation charge structure, simulates the energy and pressure propagation characteristics after the explosion excitation, obtains seismic wavelets excited by different geological conditions and excitation parameters, and optimizes the field construction seismic excitation parameters on this basis, fully utilizes the rock mechanical and physical parameters, seismic geological conditions, and on-site construction parameters, obtains seismic wavelets through nonlinear finite element explicit dynamics numerical simulation, and creates conditions for optimizing seismic excitation parameters in complex construction areas, and the entire method has a short research cycle and high efficiency. BRIEF DESCRIPTION OF THE DRAWINGS

[0053] Figure 1 The internal structure of the existing detector;

[0054] Figure 2 The earthquake excitation model of the present invention;

[0055] Figure 3 It is the rock explosion pressure cloud map of the present invention;

[0056] Figure 4 This is a cross-sectional diagram of rock explosion pressure of the present invention;

[0057] Figure 5 This is a statistical diagram of particle vibration velocity (from left to right, the distance from the explosion point is from far to near);

[0058] Figure 6 It is the amplitude-frequency diagram of seismic wavelet;

[0059] Figure 7 This is a statistical diagram of the energy of seismic wavelets simulated numerically in the present invention;

[0060] Figure 8 It is a histogram of seismic wavelet energy numerically simulated by the present invention;

[0061] Fig. 9 This is a flow chart of the method for optimizing seismic excitation parameters based on seismic wavelet analysis of the present invention;

[0062] Fig.10 This is a block diagram of the seismic excitation parameter optimization system based on seismic wavelet analysis of the present invention. DETAILED DESCRIPTION

[0063] In order to make the purpose, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Generally, the components of the embodiments of the present invention described and shown in the drawings here can be arranged and designed in various different configurations.

[0064] Therefore, the following detailed description of the embodiments of the present invention provided in the accompanying drawings is not intended to limit the scope of the invention claimed for protection, but merely represents selected embodiments of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0065] It should be noted that similar reference numerals and letters denote similar items in the following drawings, and therefore, once an item is defined in one drawing, further definition and explanation thereof is not required in subsequent drawings.

[0066] In the description of the embodiments of the present invention, it should be noted that if the terms "upper", "lower", "horizontal", "inner", etc. indicate an orientation or positional relationship based on the orientation or positional relationship shown in the drawings, or the orientation or positional relationship in which the product of the invention is usually placed when in use, it is only for the convenience of describing the present invention and simplifying the description, and does not indicate or imply that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore cannot be understood as a limitation on the present invention. In addition, the terms "first", "second", etc. are only used to distinguish the description, and cannot be understood as indicating or implying relative importance.

[0067] In addition, if the term "horizontal" appears, it does not mean that the component must be absolutely horizontal, but can be slightly tilted. For example, "horizontal" only means that its direction is more horizontal than "vertical", which does not mean that the structure must be completely horizontal, but can be slightly tilted.

[0068] In the description of the embodiments of the present invention, it is also necessary to explain that, unless otherwise clearly specified and limited, the terms "set", "install", "connect", and "connect" should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be a direct connection, or it can be indirectly connected through an intermediate medium, or it can be the internal connection of two components. For ordinary technicians in this field, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

[0069] In order to enable those skilled in the art to better understand the technical solution of the present invention, the present invention will be further described in detail below with reference to the accompanying drawings.

[0070] like Fig. 9 As shown, the present invention provides a method for optimizing seismic excitation parameters based on seismic wavelet analysis, comprising:

[0071] S101 obtains basic physical parameters of rock and establishes a constitutive model of rock according to the basic physical parameters of rock;

[0072] S102 obtains explosive parameters and establishes a state equation according to the explosive parameters;

[0073] S103 establishes a numerical model based on the state equation and the constitutive model;

[0074] S104 simulates different parameters based on the numerical model to obtain vibration velocity information of the observed particle;

[0075] S105 performs amplitude and spectrum analysis on the observed point vibration velocity information to obtain the basic characteristics of the seismic wavelet;

[0076] S106 Optimize the construction method based on the basic characteristics of the seismic wavelet.

[0077] The basic physical parameters of the rock are obtained by coring the rock in the analysis area and then using mechanical deformation tests.

[0078] The rock constitutive model adopts the Johnson-Cook constitutive relationship model.

[0079] The specific constitutive model of rock is:

[0080]

[0081] Among them, σ is the dynamic yield stress of the material; A is the quasi-static yield stress of the material; B is the strain hardening modulus of the material; ε p is the effective plastic strain; ε p * is the normalized effective plastic strain rate, ε is the plastic strain rate, T H is the homologous temperature, C, n and m are material constants.

[0082] The state equation is:

[0083]

[0084] Where p is the pressure; E is the internal energy per unit volume of the detonation products; V is the relative volume of the detonation products, that is, the ratio of the volume of the detonation products after the explosion to their initial volume; A, B, R1, R2 and ω are constants.

[0085] The numerical model is meshed with finite elements.

[0086] Based on the basic characteristics of seismic wavelets, the preferred construction method is as follows:

[0087] The energy growth trend of the seismic wavelet is obtained according to the basic characteristics of the seismic wavelet, the inflection point of the seismic excitation energy is obtained according to the energy growth trend of the seismic wavelet, and the construction method is optimized according to the inflection point of the seismic excitation energy.

[0088] Example:

[0089] The present invention mainly realizes a method for extracting seismic wavelets based on finite element numerical simulation technology and optimizing field construction parameters through seismic wavelet attribute analysis. The method mainly utilizes nonlinear finite element explicit dynamic numerical simulation technology, constructs a numerical model based on the rock physical parameters of the area adjacent to the explosion source, the structural characteristics of the near-surface low-velocity reduction layer, and the seismic excitation charge structure, simulates the energy and pressure propagation characteristics after explosion excitation, obtains seismic wavelets excited by different geological conditions and excitation parameters, and optimizes field construction seismic excitation parameters on this basis.

[0090] The implementation method of the present invention is completed according to the following steps:

[0091] (1) Core the rock in the analysis area and obtain the basic physical parameters of the rock through mechanical deformation testing.

[0092] (2) Based on the shear resistance performance parameters of rock materials, a constitutive model of rock materials is established.

[0093] (3) Select an appropriate equation of state based on the parameters of the explosive.

[0094] (4) Establish a numerical model based on the near-surface geological conditions and charge structure, and divide the finite element grid.

[0095] (5) Simulate different parameters to perform excitation simulation and obtain the vibration velocity information of the observed particle.

[0096] (6) Perform amplitude and spectrum analysis on the vibration velocity information of the observed particles in the elastic deformation area to obtain the basic characteristics of the seismic wavelet.

[0097] (7) Analyze the seismic wavelets and characteristics of different excitation parameters and construction methods, and optimize the construction method.

[0098] Specific:

[0099] Step 1: In response to the blasting impact phenomenon in rocks, the Johnson-Cook constitutive relationship model is used to describe the stress-strain relationship of the surrounding rock. It is suitable for the strength properties of large strain and high strain rate materials and is suitable for the strong impact load problems caused by high-speed collisions and explosive detonations.

[0100]

[0101] Among them, σ is the dynamic yield stress of the material; A is the quasi-static yield stress of the material; B is the strain hardening modulus of the material; εp is the effective plastic strain; εp* is the normalized effective plastic strain rate, εp*=ε / ε0, ε is the plastic strain rate, ε0 is the critical strain rate, TH is the homologous temperature, and its expression is TH=(T-T0) / (T-Tm), T, T0 and Tm are the temperatures of the material during the deformation process, the reference temperature and the melting point temperature; C, n and m are material constants.

[0102] In order to represent the plastic yield of rock under huge shear stress, the volume failure model is applied. When the maximum principal stress tensile failure strain and the maximum shear strain exceed the limit (tensile strength), the material will handle the mesh failure in an erosion manner. (For example, Table 1)

[0103] Table 1 Physical and mechanical parameters of rock formations

[0104]

[0105] Step 2: For high-energy emulsion explosives, the JWL (Jones-Wilkins-Lee) state equation is used to describe the pressure, volume and energy characteristics of the gas products during the detonation process. Its expression is:

[0106]

[0107] Where p is the pressure; E is the internal energy per unit volume of the detonation product; V is the relative volume of the detonation product, that is, the ratio of the volume of the detonation product after the explosion to its initial volume; A, B, R1, R2 and ω are constants determined by experiments. (For example, Table 2)

[0108] Table 2 Explosive materials and their state equation parameters

[0109]

[0110] Step 3: Establish an excitation numerical model and start numerical simulation calculations.

[0111] According to the investigation results of the near-surface low-velocity layer, a near-surface digital model is established. For example, the P-wave velocity of the low-velocity layer is 700m / s, the P-wave velocity of the limestone is 3398m / s, and the thickness is 60m. Figure 2 The construction method adopts single well excitation, the depth is 16.6m, and the amount of explosive increases from 6kg to 20kg, as shown in Table 3.

[0112] Table 3 The detonation point is located 5m below the surface of the high-speed layer, and the numerical simulation of 6kg-20kg excited sub-wave

[0113] Serial number Excitation parameters 1 16.6m, 6kg 2 16.6m, 8kg 3 16.6m, 10kg 4 16.6m, 12kg 5 16.6m, 16kg 6 16.6m, 20kg

[0114] Step 4: Analyze the detonation wave pressure cloud diagram ( Figure 3-Figure 4 ), extract the vibration velocity information of the observed particles below the explosion point, and analyze the amplitude spectrum characteristics of the seismic wavelet.

[0115] From the amplitude-frequency spectrum of the seismic wavelet, it can be seen that as the observation point gradually moves away from the detonation point until it enters the elastic region, the seismic wavelet energy and frequency range gradually tend to be stable.

[0116] Step 5: Analyze the amplitude-frequency characteristics of the seismic wavelet under different excitation parameters and optimize the construction parameters. As shown in the statistical chart, when the amount of seismic excitation charge reaches 12kg, the energy growth trend of the seismic wavelet changes, and 12kg is the turning point of the seismic excitation energy.

[0117] Step 6: Modify the seismic excitation numerical model according to the field acquisition parameters (such as near-surface rock physical parameters, single well or combination well excitation, combination well spacing, explosive quantity, explosive type, filler below the explosive or well stuffing material above the explosive, wellhead aperture size, etc.), so as to quickly evaluate the amplitude-frequency characteristics of the seismic wavelet.

[0118] By using the above working steps, the amplitude-frequency characteristics of seismic wavelets with different excitation parameters can be quickly obtained, providing reference information for optimizing seismic acquisition excitation parameters and improving the quality of seismic data.

[0119] like Fig.10 As shown, the present invention also provides a seismic excitation parameter optimization system based on seismic wavelet analysis, comprising:

[0120] Constitutive model acquisition module: obtains the basic physical parameters of rocks and establishes the constitutive model of rocks according to the basic physical parameters of rocks;

[0121] State equation establishment module: obtain explosive parameters and establish the state equation according to the explosive parameters;

[0122] Numerical model building module: build numerical models based on state equations and constitutive models;

[0123] Vibration velocity acquisition module: simulates different parameters based on the numerical model to obtain the vibration velocity information of the observed particle;

[0124] Vibration velocity analysis module: performs amplitude and spectrum analysis on the observed vibration velocity information to obtain the basic characteristics of the seismic wavelet;

[0125] Optimal construction module: Optimal construction method based on the basic characteristics of seismic wavelets.

[0126] A terminal device is provided in one embodiment of the present invention. The terminal device of this embodiment includes: a processor, a memory, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the steps in the above-mentioned method embodiments are implemented. Alternatively, when the processor executes the computer program, the functions of the modules / units in the above-mentioned device embodiments are implemented.

[0127] The computer program may be divided into one or more modules / units, and the one or more modules / units are stored in the memory and executed by the processor to accomplish the present invention.

[0128] The terminal device may be a computing device such as a desktop computer, a notebook, a PDA, a cloud server, etc. The terminal device may include, but is not limited to, a processor and a memory.

[0129] The processor can be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc.

[0130] The memory may be used to store the computer program and / or module, and the processor implements various functions of the terminal device by running or executing the computer program and / or module stored in the memory and calling the data stored in the memory.

[0131] If the module / unit integrated in the terminal device is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the present invention implements all or part of the processes in the above-mentioned embodiment method, and can also be completed by instructing the relevant hardware through a computer program. The computer program can be stored in a computer-readable storage medium, and the computer program can implement the steps of the above-mentioned various method embodiments when executed by the processor. Among them, the computer program includes computer program code, and the computer program code can be in source code form, object code form, executable file or some intermediate form. The computer-readable medium may include: any entity or device capable of carrying the computer program code, recording medium, U disk, mobile hard disk, disk, optical disk, computer memory, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), electric carrier signal, telecommunication signal and software distribution medium. It should be noted that the content contained in the computer-readable medium can be appropriately increased or decreased according to the requirements of legislation and patent practice in the jurisdiction. For example, in some jurisdictions, according to legislation and patent practice, computer-readable media do not include electric carrier signals and telecommunication signals.

[0132] Although the embodiments of the present invention are described above in conjunction with the accompanying drawings, the present invention is not limited to the above-mentioned specific embodiments and application fields, and the above-mentioned specific embodiments are only illustrative and instructive, rather than restrictive. Under the guidance of the specification, a person skilled in the art can also make many forms without departing from the scope of protection of the claims of the present invention, all of which belong to the scope of protection of the present invention.

Claims

1. A method for optimizing seismic excitation parameters based on seismic wavelet analysis, characterized in that: include: Obtain the basic physical parameters of rocks and establish the constitutive model of rocks according to the basic physical parameters of rocks; Obtaining explosive parameters, and establishing a state equation according to the explosive parameters; Establish numerical models based on equations of state and constitutive models; Based on the numerical model, different parameters are simulated to obtain the vibration velocity information of the observed particle; Perform amplitude and spectrum analysis on the observed point vibration velocity information to obtain the basic characteristics of the seismic wavelet; The construction method is optimized based on the basic characteristics of seismic wavelets.

2. The method for optimizing seismic excitation parameters based on seismic wavelet analysis according to claim 1, characterized in that: The basic physical parameters of the rock are obtained by coring the rock in the analysis area and then using mechanical deformation testing.

3. The method for optimizing seismic excitation parameters based on seismic wavelet analysis according to claim 1, characterized in that: The rock constitutive model adopts the Johnson-Cook constitutive relationship model.

4. The method for optimizing seismic excitation parameters based on seismic wavelet analysis according to claim 3, characterized in that: The constitutive model of the rock is specifically: Among them, σ is the dynamic yield stress of the material; A is the quasi-static yield stress of the material; B is the strain hardening modulus of the material; ε p is the effective plastic strain; ε p * is the normalized effective plastic strain rate, ε is the plastic strain rate, T H is the homologous temperature, C, n and m are material constants.

5. The method for optimizing seismic excitation parameters based on seismic wavelet analysis according to claim 1, characterized in that: The state equation is: Where p is the pressure; E is the internal energy per unit volume of the detonation products; V is the relative volume of the detonation products, that is, the ratio of the volume of the detonation products after the explosion to their initial volume; A, B, R1, R2 and ω are constants.

6. The method for optimizing seismic excitation parameters based on seismic wavelet analysis according to claim 1, characterized in that: The numerical model is divided into finite element meshes.

7. The method for optimizing seismic excitation parameters based on seismic wavelet analysis according to claim 1, characterized in that: Based on the basic characteristics of seismic wavelets, the preferred construction method is as follows: The energy growth trend of the seismic wavelet is obtained according to the basic characteristics of the seismic wavelet, the inflection point of the seismic excitation energy is obtained according to the energy growth trend of the seismic wavelet, and the construction method is optimized according to the inflection point of the seismic excitation energy.

8. A seismic excitation parameter optimization system based on seismic wavelet analysis, characterized in that: include: Constitutive model acquisition module: obtains the basic physical parameters of rocks and establishes the constitutive model of rocks according to the basic physical parameters of rocks; State equation establishment module: obtain explosive parameters and establish the state equation according to the explosive parameters; Numerical model building module: build numerical models based on state equations and constitutive models; Vibration velocity acquisition module: simulates different parameters based on the numerical model to obtain the vibration velocity information of the observed particle; Vibration velocity analysis module: performs amplitude and spectrum analysis on the observed vibration velocity information to obtain the basic characteristics of the seismic wavelet; Optimal construction module: Optimal construction method based on the basic characteristics of seismic wavelets.

9. A computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the steps of the method for optimizing seismic excitation parameters based on seismic wavelet analysis as claimed in any one of claims 1 to 7 are implemented.

10. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the steps of the method for optimizing seismic excitation parameters based on seismic wavelet analysis as claimed in any one of claims 1 to 7 are implemented.

Images