Seismic wave field analysis simulation method and device, computer equipment and storage medium

Abstract

The invention provides a seismic wave field analytical simulation method, which comprises the following steps of: constructing a seabed interface analytical equation according to a displacement potential function and a boundary condition of reflected / transmitted waves; obtaining a first function relationship; obtaining a second function relationship; converting the seabed interface analytical equation by using the first function relationship and the second function relationship to obtain a function relationship between the reflection / transmission coefficient and the influence factor, and determining the function relationship as a third function relationship; a plurality of influence factors are obtained, and reflection / transmission coefficients corresponding to the influence factors are obtained based on the third function relation; and obtaining seabed seismic wave field data through forward modeling based on the influence factors and the reflection / transmission coefficients which are in one-to-one correspondence. According to the seismic wave field analysis simulation method, the obtained seabed seismic wave field is highly similar to real data, and a theoretical model can be provided for prediction and attenuation of seabed multiples and ghost waves.

Description

Technical Field

[0001] The present invention relates to the technical field of seismic data processing, and in particular to a seismic wave field analysis simulation method, device, computer equipment and storage medium. Background Art

[0002] In marine seismic exploration, since the two sides of the seabed interface are liquid and solid media respectively, the reflection and transmission of seismic waves at this interface are more complicated than those at the solid-solid interface. On the one hand, the seabed interface is a velocity discontinuity surface for seismic waves, and the wave impedance on both sides of the interface is very different. Seismic waves show the characteristics of acoustic wave propagation in the seawater layer of the upper interface and the characteristics of solid layer propagation at the lower interface, which makes the calculation of reflection and transmission coefficients at the seabed interface very complicated. On the other hand, the reflection / transmission coefficients of seabed interfaces with different medium properties vary with the incident angle. Not only do the reflection / transmission coefficients of hard and soft seabed interfaces have different characteristics, but the longitudinal and transverse wave velocities of the solid layer will also affect the reflection and transmission coefficient characteristics of the seabed interface.

[0003] Studying the characteristics of the seafloor reflection / transmission coefficient can help us better understand the properties of the seafloor medium, and is of great significance for the processing and interpretation of seafloor multi-component seismic data. Since the shear modulus in the seawater layer is zero, and there are shear waves in the seafloor solid layer, the seafloor interface can be regarded as a secondary shear wave source. Using the seafloor reflection / transmission coefficient to obtain the shear wave velocity information of the seafloor solid layer is conducive to the study of seafloor conversion wave problems. In addition, compared with land, marine seismic exploration, the multiple waves and ghost waves in the seawater layer are very complex. The study of the reflection / transmission coefficient of the seafloor interface can better understand the characteristics of multiple waves and ghost waves in theory, which is conducive to the prediction and attenuation of multiple waves and accompanying ghost waves in marine seismic data processing.

[0004] Although numerical algorithms can calculate complex seabed models in solving the nonlinear problem of seabed reflection / transmission coefficient, they cannot analytically discuss the impact of influencing factors on the reflection / transmission coefficient, and lack an understanding of the law of how the reflection / transmission coefficient changes with changes in influencing factors. Summary of the invention

[0005] Based on this, it is necessary to provide a seismic wave field analysis simulation method, device, computer equipment and storage medium to address the above technical problems.

[0006] A seismic wave field analysis simulation method, comprising:

[0007] According to the displacement potential function and boundary conditions of the reflected / transmitted waves, the analytical equation of the seabed interface is constructed;

[0008] Obtaining a functional relationship between a displacement potential function of the reflection / transmission wave and a reflection / transmission coefficient, and determining the relationship as a first functional relationship;

[0009] Obtaining a functional relationship between the boundary condition and the influencing factor, and determining it as a second functional relationship;

[0010] The first functional relationship and the second functional relationship are used to transform the seabed interface analytical equation to obtain a functional relationship between the reflection / transmission coefficient and the influencing factor, which is determined as a third functional relationship;

[0011] Acquire multiple influencing factors, and acquire the reflection / transmission coefficient corresponding to each of the influencing factors based on the third functional relationship;

[0012] Based on the one-to-one correspondence between the influencing factors and the reflection / transmission coefficients, forward simulation is performed to obtain seafloor seismic wave field data.

[0013] In one embodiment, the influencing factors include incident angle and seabed medium parameters;

[0014] The step of obtaining a plurality of influencing factors and obtaining a reflection / transmission coefficient corresponding to each of the influencing factors based on the third functional relationship comprises:

[0015] Acquire the seabed medium parameters, and obtain the functional relationship between the reflection / transmission coefficient and the incident angle based on the third functional relationship, and determine it as a fourth functional relationship;

[0016] The incident angle is obtained, and based on the fourth functional relationship, the reflection / transmission coefficient corresponding to each incident angle is obtained.

[0017] In one embodiment, the seabed medium parameters include a wave impedance data set and a solid interface transmission angle data set.

[0018] In one embodiment, the wave impedance data set includes liquid P-wave wave impedance, solid P-wave wave impedance and solid S-wave wave impedance;

[0019] The solid interface transmission angle data group includes a solid interface P-wave transmission angle and a solid interface S-wave transmission angle.

[0020] In one embodiment, the transmission coefficient includes a longitudinal wave transmission coefficient and a converted wave transmission coefficient.

[0021] In one embodiment, the third functional relationship is:

[0022]

[0023] In the formula, R PP is the reflection coefficient, T PP is the longitudinal wave transmission coefficient, T PSis the conversion wave transmission coefficient, I1 is the wave impedance of P wave in the seabed liquid medium, I2 is the wave impedance of P wave in the seabed solid medium, and I3 is the wave impedance of S wave in the seabed solid medium. i1 is the incident angle of the downward P wave at the seabed liquid interface, i2 is the transmission angle of the P wave transmitted by the seabed solid interface, and j2 is the transmission angle of the S wave transmitted by the seabed solid interface; G is the seabed analytical variable, and G is related to I1, I2, I3, i1, i2 and j2 respectively.

[0024] In one embodiment, the step of obtaining the functional relationship between the displacement potential function of the reflection / transmission wave and the reflection / transmission coefficient and determining the first functional relationship includes:

[0025] Acquire a functional relationship between a displacement potential function of the reflection / transmission wave and a displacement amplitude of the reflection / transmission wave, and determine it as a fifth functional relationship;

[0026] Using the fifth functional relationship, the displacement potential function of the reflection / transmission wave in the seabed interface analytical equation is converted into the displacement amplitude of the reflection / transmission wave;

[0027] Acquire a functional relationship between the displacement amplitude of the reflection / transmission wave and the reflection / transmission coefficient, and determine it as a sixth functional relationship;

[0028] The sixth functional relationship is used to convert the displacement amplitude of the reflection / transmission wave in the seabed interface analytical equation into the reflection / transmission coefficient.

[0029] A seismic wave field analysis simulation device, comprising:

[0030] Construct a model module and construct the analytical equation of the seabed interface based on the displacement potential function and boundary conditions of the reflection / transmission waves;

[0031] A first functional relationship acquisition module is used to acquire the functional relationship between the displacement potential function of the reflection / transmission wave and the reflection / transmission coefficient, and determine it as a first functional relationship;

[0032] A second functional relationship acquisition module is used to acquire the functional relationship between the boundary condition and the influencing factor and determine it as a second functional relationship;

[0033] A data conversion module, using the first functional relationship and the second functional relationship, transforms the seabed interface analytical equation to obtain a functional relationship between the reflection / transmission coefficient and the influencing factor, and determines it as a third functional relationship;

[0034] A data acquisition module, which acquires a plurality of influencing factors, and acquires a reflection / transmission coefficient corresponding to each of the influencing factors based on the third functional relationship;

[0035] The forward simulation module obtains the seafloor seismic wave field data through forward simulation based on the one-to-one correspondence between the influencing factors and the reflection / transmission coefficients.

[0036] A computer device comprises a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the following steps are implemented:

[0037] According to the displacement potential function and boundary conditions of the reflected / transmitted waves, the analytical equation of the seabed interface is constructed;

[0038] Obtaining a functional relationship between a displacement potential function of the reflection / transmission wave and a reflection / transmission coefficient, and determining the relationship as a first functional relationship;

[0039] Obtaining a functional relationship between the boundary condition and the influencing factor, and determining it as a second functional relationship;

[0040] The first functional relationship and the second functional relationship are used to transform the seabed interface analytical equation to obtain a functional relationship between the reflection / transmission coefficient and the influencing factor, which is determined as a third functional relationship;

[0041] Acquire multiple influencing factors, and acquire the reflection / transmission coefficient corresponding to each of the influencing factors based on the third functional relationship;

[0042] Based on the one-to-one correspondence between the influencing factors and the reflection / transmission coefficients, seafloor seismic wave field data is obtained through forward simulation.

[0043] A computer-readable storage medium stores a computer program, which, when executed by a processor, implements the following steps:

[0044] According to the displacement potential function and boundary conditions of the reflected / transmitted waves, the analytical equation of the seabed interface is constructed;

[0045] Obtaining a functional relationship between a displacement potential function of the reflection / transmission wave and a reflection / transmission coefficient, and determining the relationship as a first functional relationship;

[0046] Obtaining a functional relationship between the boundary condition and the influencing factor, and determining it as a second functional relationship;

[0047] The first functional relationship and the second functional relationship are used to transform the seabed interface analytical equation to obtain a functional relationship between the reflection / transmission coefficient and the influencing factor, which is determined as a third functional relationship;

[0048] Acquire multiple influencing factors, and acquire the reflection / transmission coefficient corresponding to each of the influencing factors based on the third functional relationship;

[0049] Based on the one-to-one correspondence between the influencing factors and the reflection / transmission coefficients, seafloor seismic wave field data is obtained through forward simulation.

[0050] The seabed seismic wave field refers to the seismic wave field based on the liquid-solid interface of the seabed.

[0051] The reflected / transmitted waves include the reflected waves and the transmitted waves.

[0052] The reflected / transmitted waves of the present application include but are not limited to the reflected / transmitted waves obtained by the incidence of downward P waves at the seabed liquid-solid interface.

[0053] The boundary condition may be a boundary condition for the incidence of a downward P wave at the seafloor liquid-solid interface.

[0054] Compared with the prior art, the seismic wave field analysis simulation method provided by the present invention has the following advantages or beneficial effects:

[0055] 1) According to the displacement potential function and boundary conditions of the reflection / transmission wave, the analytical equation of the seabed interface is constructed, and the analytical equation of the seabed interface is transformed, and the functional relationship between the reflection / transmission coefficient and the influencing factor is constructed, so as to find the law of the change of the reflection / transmission coefficient with the change of the influencing factor, and the influence of the influencing factor on the reflection / transmission coefficient can be quantitatively discussed.

[0056] 2) Based on the one-to-one correspondence between the influencing factors and the reflection / transmission coefficients, seabed seismic wave field data are obtained through forward simulation. The seabed seismic wave field data obtained through forward simulation are compared with the actual seabed seismic wave field data to verify the accuracy of the third functional relationship. The seabed seismic wave field data obtained through forward simulation can provide a theoretical model for the prediction and attenuation of seabed multiple waves and ghost waves. BRIEF DESCRIPTION OF THE DRAWINGS

[0057] Figure 1 is a flow chart of a seismic wave field analytical simulation method in an embodiment;

[0058] Figure 2 is a structural block diagram of a seismic wave field analysis simulation device in an embodiment;

[0059] Figure 3 is a schematic diagram of boundary conditions in an embodiment;

[0060] Figure 4a is a functional graph of the incident angle and the reflection coefficient amplitude value of the liquid-solid interface of the hard seabed simulated based on the third functional relationship in an embodiment;

[0061] Figure 4b is a functional graph of the incident angle and the amplitude of the longitudinal wave transmission coefficient of the hard seabed liquid-solid interface simulated based on the third functional relationship in an embodiment;

[0062] Figure 4c is a functional image of the incident angle of the liquid-solid interface of the hard seabed and the amplitude value of the converted wave transmission coefficient simulated based on the third functional relationship in an embodiment;

[0063] Figure 4d is a functional image of the incident angle and the reflection coefficient phase angle of the liquid-solid interface of the hard seabed simulated based on the third functional relationship in an embodiment;

[0064] Figure 4e is a functional image of the incident angle and the phase angle of the longitudinal wave transmission coefficient of the hard seabed liquid-solid interface simulated based on the third functional relationship in an embodiment;

[0065] Figure 4f is a functional image of the incident angle and the phase angle of the converted wave transmission coefficient of the hard seabed liquid-solid interface simulated based on the third functional relationship in an embodiment;

[0066] Figure 5a is a functional graph of the incident angle and the reflection coefficient amplitude value of the soft seabed liquid-solid interface simulated based on the third functional relationship in an embodiment;

[0067] Figure 5b is a functional image of the incident angle and the amplitude of the longitudinal wave transmission coefficient of the soft seabed liquid-solid interface simulated based on the third functional relationship in an embodiment;

[0068] Figure 5c is a functional image of the incident angle of the soft seabed liquid-solid interface and the amplitude value of the converted wave transmission coefficient simulated based on the third functional relationship in an embodiment;

[0069] Figure 5d is a functional image of the incident angle and the reflection coefficient phase angle of the soft seabed liquid-solid interface simulated based on the third functional relationship in an embodiment;

[0070] Figure 5e is a functional image of the incident angle and the phase angle of the longitudinal wave transmission coefficient of the soft seabed liquid-solid interface simulated based on the third functional relationship in an embodiment;

[0071] Figure 5f is a functional image of the incident angle and the phase angle of the converted wave transmission coefficient of the soft seabed liquid-solid interface simulated based on the third functional relationship in an embodiment;

[0072] Figure 6 is a comparison diagram of seabed seismic wave field data obtained by forward modeling in an embodiment and actual seabed seismic wave field data;

[0073] Figure 7 FIG. 4 is a diagram showing the internal structure of a computer device in one embodiment. DETAILED DESCRIPTION

[0074] In order to make the purpose, technical solution and advantages of the present application more clearly understood, the present application is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and are not used to limit the present application.

[0075] Embodiment 1

[0076] In this embodiment, Figure 1 As shown, a seismic wave field analytical simulation method is provided, which includes:

[0077] Step 110, constructing the seabed interface analytical equation according to the displacement potential function of the reflection / transmission wave and the boundary conditions.

[0078] In this embodiment, the boundary condition is a constraint condition of the displacement potential function of the reflection / transmission wave. In each embodiment, the reflection / transmission wave includes the reflection wave and the transmission wave, and accordingly, the displacement potential function of the reflection / transmission wave includes the displacement potential function of the reflection wave and the displacement potential function of the transmission wave.

[0079] In this embodiment, if Figure 3 As shown in the figure, the horizontal line between fluid and solid represents the liquid-solid interface of the seabed; fluid represents the seabed liquid medium, that is, seawater, ρ1 is the density of the seabed liquid medium, and α1 is the P-wave velocity of seawater; solid represents the seabed solid medium, including hard seabed medium and soft seabed medium, ρ2 is the density of the seabed solid medium, α2 is the P-wave velocity of the seabed solid medium, and β2 is the S-wave velocity of the seabed solid medium.

[0080] represents the displacement amplitude of the incident wave in the liquid medium,

[0081] represents the displacement amplitude of the incident P wave in the solid medium,

[0082] represents the displacement amplitude of the incident S-wave in the solid medium,

[0083] represents the displacement amplitude of the reflected wave in the liquid medium,

[0084] represents the displacement amplitude of the transmitted P wave in the solid medium,

[0085] It represents the displacement amplitude of the transmitted S wave in the solid medium.

[0086] In this embodiment, the boundary condition is expressed by the following calculation formula:

[0087]

[0088] Where φ1 is the displacement potential function of the reflected / transmitted wave, φ2 is the displacement potential function of the incident wave, z is the depth, λ1, μ1, λ2 and μ2 are all Lame constants, which are material-related quantities that will appear in the strain-stress relationship.

[0089] In this embodiment, an analytical equation of the seabed interface for the reflected wave is constructed based on the displacement potential function and boundary conditions of the reflected wave.

[0090] In this embodiment, an analytical equation of the seabed interface for the transmission wave is constructed based on the displacement potential function and boundary conditions of the transmission wave.

[0091] Step 120: Obtain a functional relationship between the displacement potential function of the reflection / transmission wave and the reflection / transmission coefficient, and determine it as a first functional relationship.

[0092] In this embodiment, the analytical equation of the seabed interface is constructed based on the displacement potential function of the reflection / transmission wave. To obtain the functional relationship between the reflection / transmission coefficient and the influencing factor, it is necessary to convert the displacement potential function of the reflection / transmission wave into the reflection / transmission coefficient. Therefore, it is necessary to obtain the functional relationship between the displacement potential function of the reflection / transmission wave and the reflection / transmission coefficient and determine it as the first functional relationship.

[0093] In this embodiment, a functional relationship between a displacement potential function of the reflected wave and a reflection coefficient is obtained and determined as a first functional relationship with respect to the reflected wave.

[0094] In this embodiment, a functional relationship between the displacement potential function and the transmission coefficient of the transmitted wave is obtained and determined as a first functional relationship with respect to the transmitted wave.

[0095] Step 130: Obtain a functional relationship between the boundary condition and the influencing factor, and determine it as a second functional relationship.

[0096] In this embodiment, the analytical equation of the seabed interface is constructed based on the boundary conditions. To obtain the functional relationship between the reflection / transmission coefficient and the influencing factor, the boundary condition needs to be converted into the influencing factor. Therefore, it is necessary to obtain the functional relationship between the boundary condition and the influencing factor and determine it as the second functional relationship.

[0097] Step 140, using the first functional relationship and the second functional relationship, transforming the seabed interface analytical equation to obtain a functional relationship between the reflection / transmission coefficient and the influencing factor, and determining it as a third functional relationship.

[0098] In this embodiment, the third functional relationship regarding the reflected wave is obtained by using the first functional relationship regarding the reflected wave and the second functional relationship.

[0099] In this embodiment, the first functional relationship regarding the transmitted wave and the second functional relationship are used to obtain the third functional relationship regarding the transmitted wave.

[0100] Step 150: Obtain multiple influencing factors, and based on the third functional relationship, obtain the reflection / transmission coefficient corresponding to each of the influencing factors.

[0101] In this embodiment, the reflection / transmission coefficient corresponding to each of the influencing factors is obtained to obtain multiple data groups, which are used as data for subsequent forward simulation to obtain seabed seismic wave field data; and by observing each influencing factor and its corresponding reflection / transmission coefficient, the law of change of the reflection / transmission coefficient with the change of the influencing factor can be obtained.

[0102] In this embodiment, the influence of the influencing factor on the reflection / transmission coefficient can also be quantitatively discussed.

[0103] Step 160: Based on the one-to-one correspondence between the influencing factors and the reflection / transmission coefficients, obtain seafloor seismic wave field data through forward modeling.

[0104] In this embodiment, based on the one-to-one correspondence between the influencing factors and the reflection / transmission coefficients, seabed seismic wavefield data is obtained through forward simulation. The seabed seismic wavefield data obtained through forward simulation is compared with the actual seabed seismic wavefield data to verify the accuracy of the third functional relationship. The seabed seismic wavefield data obtained through forward simulation can provide a theoretical model for the prediction and attenuation of seabed multiple waves and ghost waves.

[0105] There are many influencing factors that affect the reflection / transmission coefficient. When observing the relationship between the reflection / transmission coefficient and the influencing factors, some of the influencing factors can be set as constants in order to observe the influence of other influencing factors on the reflection / transmission coefficient, so as to find the relationship between the reflection / transmission coefficient and the influencing factors.

[0106] In one embodiment, the influencing factors include incident angle and seabed medium parameters;

[0107] The step of obtaining a plurality of influencing factors and obtaining a reflection / transmission coefficient corresponding to each of the influencing factors based on the third functional relationship comprises:

[0108] Step 151, acquiring the seabed medium parameters, and obtaining the functional relationship between the reflection / transmission coefficient and the incident angle based on the third functional relationship, and determining it as a fourth functional relationship;

[0109] Step 152, obtaining the incident angle, and based on the fourth functional relationship, obtaining the reflection / transmission coefficient corresponding to each incident angle.

[0110] In the above embodiment, after setting the seabed medium parameters as constants, the functional relationship between the reflection / transmission coefficient and the incident angle is obtained, and then the law of the reflection / transmission coefficient changing with the incident angle is obtained.

[0111] The seabed medium parameters in different sea areas are different, for example, they can be divided into hard seabed medium and soft seabed medium, and the parameters between different hard seabed media are different, and the parameters between different soft seabed media are different. At least two embodiments can be set, that is, at least one hard seabed medium and at least one soft seabed medium can be selected as embodiments respectively, that is, each embodiment obtains the seabed medium parameters of the corresponding sea area, and substitutes them into the third functional relationship to find out the law of the reflection / transmission coefficient of the sea area with the change of the incident angle. The comparison between at least two embodiments can also derive the influence of the seabed medium parameters on the reflection / transmission coefficient.

[0112] In one embodiment, the seabed medium parameters include a wave impedance data set and a solid interface transmission angle data set.

[0113] In this embodiment, in the seismic wave field based on the seabed liquid-solid interface, the reflection / transmission coefficient is affected by many factors. Therefore, in the functional relationship between the reflection / transmission coefficient and the influencing factors, the more types of influencing factors there are, the more likely it is that the simulation can be close to the actual law of the reflection / transmission coefficient and the influencing factors, that is, the simulated seismic wave field is highly similar to the actual seismic wave field.

[0114] In one embodiment, the wave impedance data set includes liquid P-wave wave impedance, solid P-wave wave impedance and solid S-wave wave impedance.

[0115] In this embodiment, in the seismic wave field based on the liquid-solid interface of the seabed, the liquid P-wave impedance, solid P-wave impedance and solid S-wave impedance all affect the reflection / transmission coefficient. Therefore, a functional relationship between the reflection / transmission coefficient and the influencing factor is constructed based on the liquid P-wave impedance, solid P-wave impedance and solid S-wave impedance, so that the simulated seismic wave field is highly similar to the actual seismic wave field, providing a theoretical model for the prediction and attenuation of seabed multiple waves and ghost waves, which can be used for the simulation and prediction of the seabed interface seismic wave field, and significantly improving the accuracy of seismic data denoising.

[0116] In one embodiment, the solid interface transmission angle data set includes a solid interface P-wave transmission angle and a solid interface S-wave transmission angle.

[0117] In this embodiment, in the seismic wave field based on the seabed liquid-solid interface, the solid interface P-wave transmission angle and the solid interface S-wave transmission angle both affect the reflection / transmission coefficient. Therefore, a functional relationship between the reflection / transmission coefficient and the influencing factor is constructed based on the solid interface P-wave transmission angle and the solid interface S-wave transmission angle, so that the simulated seismic wave field is highly similar to the actual seismic wave field, providing a theoretical model for the prediction and attenuation of seabed multiple waves and ghost waves, which can be used for the simulation and prediction of the seabed interface seismic wave field, and significantly improving the accuracy of seismic data denoising.

[0118] The longitudinal wave transmission coefficient is an indicator for evaluating the ability of seismic waves to penetrate the strata in seismic exploration. It indicates the ratio of the amplitude of a seismic longitudinal wave that reaches a receiving point after passing through a certain interface when propagating in an underground medium to the amplitude of the wave incident on the interface. The converted wave transmission coefficient refers to the SV wave transmission coefficient. The SV wave transmission coefficient refers to the SV wave transmission coefficient after the seismic SV wave propagates through the underground medium to the receiving point. The SV wave transmission coefficient can be used to evaluate the propagation and conversion of SV waves in seismic exploration. The longitudinal wave transmission coefficient and the converted wave transmission coefficient have different meanings for seismic waves. In order to improve the fidelity of the simulated seismic wave field, in one embodiment, the transmission coefficient includes the longitudinal wave transmission coefficient and the converted wave transmission coefficient.

[0119] In this embodiment, the combination of the longitudinal wave transmission coefficient and the converted wave transmission coefficient can enrich the relevant data of the seismic wave field, and effectively improve the accuracy of the seismic wave field obtained through subsequent forward simulation.

[0120] The reflection / transmission coefficient is affected by many influencing factors. When studying the law between the reflection / transmission coefficient and the influencing factors, the extensiveness of the influencing factors involved is positively correlated with the accuracy of the reflection / transmission coefficient. In order to make the reflection / transmission coefficient obtained by the above method highly close to the actual data, in one embodiment, the third functional relationship is the following calculation formula:

[0121]

[0122] In the formula, R PP is the reflection coefficient, T PP is the longitudinal wave transmission coefficient, T PS is the conversion wave transmission coefficient, I1 is the wave impedance of P wave in the seabed liquid medium, I2 is the wave impedance of P wave in the seabed solid medium, and I3 is the wave impedance of S wave in the seabed solid medium. i1 is the incident angle of the downward P wave at the seabed liquid interface, i2 is the transmission angle of the P wave transmitted by the seabed solid interface, and j2 is the transmission angle of the S wave transmitted by the seabed solid interface;

[0123] Among them, G is the seafloor analytical variable, the specific form is

[0124] G=2sin j2sin2j2cos i1(I2cos j2-I3cos i2)-I1cos i2-I2cos i1.

[0125] In this embodiment, I1=ρ1α1, I2=ρ2α2, I3=ρ2β2, where ρ1 is the density of the seabed liquid medium, α1 is the P-wave velocity of seawater; ρ2 is the density of the seabed solid medium, α2 is the P-wave velocity of the seabed solid medium, and β2 is the S-wave velocity of the seabed solid medium.

[0126] In one embodiment, the step of obtaining the functional relationship between the displacement potential function of the reflection / transmission wave and the reflection / transmission coefficient and determining the first functional relationship includes:

[0127] Step 121, obtaining a functional relationship between a displacement potential function of the reflection / transmission wave and a displacement amplitude of the reflection / transmission wave, and determining the fifth functional relationship;

[0128] Step 122, using the fifth functional relationship, converting the displacement potential function of the reflection / transmission wave in the seabed interface analytical equation into the displacement amplitude of the reflection / transmission wave;

[0129] Step 123, obtaining a functional relationship between the displacement amplitude of the reflection / transmission wave and the reflection / transmission coefficient, and determining it as a sixth functional relationship;

[0130] Step 124, using the sixth functional relationship, converting the displacement amplitude of the reflection / transmission wave in the seabed interface analytical equation into the reflection / transmission coefficient.

[0131] In this embodiment, the displacement potential function of the reflection / transmission wave in the seabed analytical equation is converted into the reflection / transmission coefficient through the displacement amplitude of the reflection / transmission wave, and then the functional relationship between the reflection / transmission coefficient and the influencing factor is obtained, so as to obtain the law of change of the reflection / transmission coefficient with the change of the influencing factor.

[0132] It should be understood that although Figure 1 The steps in the flowchart are shown in sequence as indicated by the arrows, but these steps are not necessarily executed in the order indicated by the arrows. Unless otherwise specified in this document, there is no strict order restriction for the execution of these steps, and these steps can be executed in other orders. Moreover, Figure 1 At least part of the steps may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily executed at the same time, but can be executed at different times. The execution order of these sub-steps or stages is not necessarily sequential, but can be executed in turn or alternately with other steps or at least part of the sub-steps or stages of other steps.

[0133] Embodiment 2

[0134] In this embodiment, a method for analyzing and simulating seismic wave fields based on the liquid-solid interface of the hard seabed is provided. The liquid-solid interface of the hard seabed is located between the liquid medium and the hard seabed medium. The specific implementation is as follows:

[0135] The analytical equation of the seabed interface is constructed based on the boundary conditions of the downward P wave incident on the seabed liquid-solid interface and the displacement potential function on both sides of the seabed liquid-solid interface; the boundary conditions are expressed by the following calculation formula:

[0136]

[0137] Where φ1 is the displacement potential function of the reflected / transmitted wave, φ2 is the displacement potential function of the incident wave, z is the depth, λ1, μ1, λ2 and μ2 are all Lame constants, which are material-related quantities that will appear in the strain-stress relationship.

[0138] Solve the seabed interface analytical equation to obtain the analytical formula of the seabed liquid-solid interface reflection / transmission coefficient; among them, the analytical formula of the seabed liquid-solid interface reflection / transmission coefficient is as follows:

[0139]

[0140] The reflection / transmission coefficient of the liquid-solid interface at the hard seabed is obtained according to the analytical formula of the reflection / transmission coefficient of the liquid-solid interface at the seabed;

[0141] According to the obtained reflection / transmission coefficient of the liquid-solid interface of the hard seabed, the seabed seismic wave field data are obtained through forward simulation.

[0142] In this embodiment, the seabed medium parameters used are hard seabed medium parameters, where the P wave velocity in the seawater layer is α1=1490m / s, and the seawater density ρ1 is 1025kg / m 3 In the hard seabed solid medium, the P wave velocity is α2 = 4000 m / s, the S wave velocity is β2 = 1920 m / s, and the solid layer density is ρ2 = 2460 kg / m 3 , using the following calculation formula:

[0143]

[0144] Obtaining the reflection / transmission coefficient corresponding to each of the influencing factors;

[0145] Based on the one-to-one correspondence between the influencing factors and the reflection / transmission coefficients, seafloor seismic wave field data is obtained through forward simulation.

[0146] At the hard seabed interface, since the P-wave velocity and S-wave velocity in the solid medium of the hard seabed are both greater than the sound wave velocity in the seawater, the sound wave velocity in the seawater has little interference with the P-wave velocity and the S-wave velocity in the solid medium of the hard seabed, so that the reflected P-wave, transmitted P-wave and transmitted S-wave all follow the law obtained by the above method, that is, the reflection coefficient R PP And the transmission coefficient T PP and T PS Both of them show the characteristics of variation with two critical transmission angles as the dividing points. The two critical transmission angles include the critical angle of P wave transmission through the hard seabed interface and the critical angle of S wave transmission through the hard seabed interface.

[0147] According to the following calculation formula:

[0148]

[0149] The critical angle for P wave transmission through the hard seabed interface is 22°, and the critical angle for S wave transmission through the hard seabed interface is 50°.

[0150] Specifically, Figure 4a As shown in the figure, Incidence angle is the incident angle, degrees is the degree, and Amplitude is the amplitude value. The same is true for other figures, which will not be repeated in the following text. Figure 5a As shown in the figure, Phase angle is the phase angle. The same is true for other figures, which will not be elaborated in the following text.

[0151] Figure 4a The function graph of the incident angle and the reflection coefficient amplitude value of the liquid-solid interface of the hard seabed simulated based on the third functional relationship, Figure 4a It reflects the change of reflection coefficient amplitude with the incident angle. Figure 4b It is a function image of the incident angle and the amplitude of the longitudinal wave transmission coefficient of the liquid-solid interface of the hard seabed simulated based on the third functional relationship. Figure 4c It is a function image of the incident angle and the converted wave transmission coefficient of the liquid-solid interface of the hard seabed simulated based on the third functional relationship. Figure 4b and Figure 4c It reflects the change of the amplitude of the transmission coefficient with the incident angle.

[0152] from Figure 4a It can be seen that the reflection coefficient R PP The amplitude reaches 1 at incident angles of 22° and 50°, and remains 1 at incident angles greater than 50°.

[0153] like Figure 4b As shown, when the incident angle is greater than or equal to 22°, the amplitude value of the longitudinal wave transmission coefficient is 0.

[0154] like Figure 4cAs shown, when the incident angle is 22°, the amplitude value of the converted wave transmission coefficient reaches a minimum value, and when the incident angle is greater than or equal to 50°, the amplitude value of the converted wave transmission coefficient is 0.

[0155] like Figure 4d As shown in the figure, when the incident angle is less than 22°, the reflection coefficient phase angle is 0. When the incident angle is between 22° and 50°, the reflection coefficient phase angle generally shows a downward trend, and the decline is small. When the incident angle is between 50° and 90°, the reflection coefficient phase angle generally shows a trend of first decreasing and then increasing, and the change trend is large. That is, when the incident angle is 22° and 50°, the reflection coefficient phase angle has a significant turning point.

[0156] like Figure 4e As shown, the longitudinal wave transmission coefficient is always 0, that is, the longitudinal wave transmission coefficient does not change with the change of the incident angle.

[0157] like Figure 4f As shown in the figure, when the incident angle is 0° to 22°, the phase angle of the converted wave transmission coefficient remains unchanged. When the incident angle is 22° to 50°, the amplitude of the converted wave transmission coefficient first decreases and then increases until it reaches 0. When the incident angle is greater than or equal to 50°, the amplitude of the converted wave transmission coefficient is 0. That is, when the incident angle is 22° and 50°, the phase angle of the converted wave transmission coefficient has an obvious turning point.

[0158] The longitudinal wave transmission coefficient T at the P-wave critical angle PP 0, and then as the incident angle continues to increase, the transmitted P wave becomes an inhomogeneous surface wave, which propagates along the seabed interface and its amplitude decays exponentially with depth; as the incident angle increases further, at the critical angle of the S wave T PS is 0, and then the S wave becomes an inhomogeneous surface wave propagating along the seabed interface. PP and T PS They all show the characteristics of change with two transmission critical angles as the dividing points.

[0159] Figure 6 It is a comparison diagram of the seabed seismic wave field data obtained by forward simulation in an embodiment and the actual seabed seismic wave field data, wherein the left figure Synthetics is the seabed seismic wave field data obtained by forward simulation, and the right figure Realdata is the actual seabed seismic wave field data. In this embodiment, the actual seabed seismic wave field data is the real data of a certain marine work area.

[0160] In this embodiment, the seabed liquid-solid interface is a hard seabed interface, and the P wave velocity of the solid layer is greater than the S wave velocity of the seawater layer. Although the horizontal velocity model is used, Figure 6As shown in the figure, the trends and amplitude variation characteristics of the seafloor seismic wave fields synthesized by forward modeling of multiple data sets obtained by the above method are highly similar to the real data of a certain marine work area. On the one hand, the trends and amplitudes of multiple waves of different orders are similar to the actual data, and the multiple waves with incident angles close to the critical angle of longitudinal waves or greater than the critical angle of transverse waves appear as seismic waves with stronger amplitudes. On the other hand, the records of the synthetic seafloor seismic wave fields simulate the medium non-uniform waves of the actual seafloor data, which are similar to refracted waves but with higher speeds.

[0161] The seafloor seismic wave field data obtained by forward modeling are very close to the actual seafloor seismic wave field data, proving that the multiple data sets obtained by the seismic wave field analytical simulation method can effectively simulate the seismic wave field of the seafloor liquid-solid interface. In other words, the seafloor seismic wave field data obtained by the above method has a good denoising accuracy.

[0162] Furthermore, the seismic wave field of the seafloor liquid-solid interface simulated by the above method can provide a theoretical model for the prediction and attenuation of seafloor multiple waves and ghost waves, significantly improving the accuracy of seismic data denoising.

[0163] Embodiment 3

[0164] In this embodiment, a seismic wave field analysis simulation method based on the soft seabed liquid-solid interface is provided. The soft seabed liquid-solid interface is located between the liquid medium and the soft seabed medium. The specific implementation is as follows:

[0165] The analytical equation of the seabed interface is constructed based on the boundary conditions of the downward P wave incident on the seabed liquid-solid interface and the displacement potential function on both sides of the seabed liquid-solid interface; the boundary conditions are expressed by the following calculation formula:

[0166]

[0167] Where φ1 is the displacement potential function of the reflected / transmitted wave, φ2 is the displacement potential function of the incident wave, z is the depth, λ1, μ1, λ2 and μ2 are all Lame constants, which are material-related quantities that will appear in the strain-stress relationship.

[0168] Solve the seabed interface analytical equation to obtain the analytical formula of the seabed liquid-solid interface reflection / transmission coefficient; among them, the analytical formula of the seabed liquid-solid interface reflection / transmission coefficient is as follows:

[0169]

[0170] The reflection / transmission coefficient of the liquid-solid interface of the soft seabed is obtained according to the analytical formula of the reflection / transmission coefficient of the liquid-solid interface of the seabed;

[0171] According to the reflection / transmission coefficient of the liquid-solid interface of the soft seabed, the seabed seismic wave field data are obtained through forward simulation.

[0172] In this embodiment, the seabed medium parameters used are soft seabed medium parameters, where the P wave velocity in the seawater layer is α1=1490m / s, and the seawater density ρ1 is 1025kg / m 3 , the P wave velocity in the soft seabed solid medium is α2 = 1800 m / s, the S wave velocity in the soft seabed solid medium is β2 = 522 m / s, and the density of the soft seabed solid medium is ρ2 = 1600 kg / m 3 , using the following calculation formula:

[0173]

[0174] Obtaining the reflection / transmission coefficient corresponding to each of the influencing factors;

[0175] Based on the one-to-one correspondence between the influencing factors and the reflection / transmission coefficients, seafloor seismic wave field data is obtained through forward simulation.

[0176] The density of the soft seabed solid medium is relatively small, and the S wave velocity propagating in the soft seabed solid medium is relatively low. Generally, the S wave velocity of the soft seabed liquid-solid interface is lower than the sound wave velocity in seawater. The sound wave velocity in seawater has a greater interference with the S wave velocity of the soft seabed liquid-solid interface. At the soft seabed liquid-solid interface, the S wave critical angle does not exist; and the P wave velocity of the soft seabed liquid-solid interface is usually higher than the sound wave velocity in seawater. The sound wave velocity in seawater has a smaller interference with the P wave velocity of the soft seabed liquid-solid interface. Therefore, at the soft seabed liquid-solid interface, the change characteristics of the reflection coefficient and the transmission coefficient with the incident angle are mainly controlled by the critical angle of the transmitted P wave. In this embodiment, the critical angle of the soft seabed interface for transmitting P waves is different from the critical angle of the hard seabed interface for transmitting P waves. In some other embodiments, the critical angle of the soft seabed interface for transmitting P waves is the same as the critical angle of the hard seabed interface for transmitting P waves.

[0177] according to

[0178]

[0179] At the critical angle of the P wave, the transmission angle i2 of the transmitted P wave is 90°, and the transmission coefficient is 0. Figure 5b and Figure 5c As shown in the figure, the incident angle when the transmission coefficient is 0 is the P-wave critical angle. When the incident angle reaches the P-wave critical angle, the transmitted P-wave becomes an inhomogeneous surface wave.

[0180] Specifically, Figure 5d It is a function image of the incident angle and the phase angle of the reflection coefficient of the soft seabed liquid-solid interface simulated based on the third functional relationship. Figure 5d It represents the change of phase angle of reflection coefficient of liquid-solid interface of soft seabed with incident angle; Figure 5e It is a functional image of the incident angle and the phase angle of the longitudinal wave transmission coefficient of the soft seabed liquid-solid interface simulated based on the third functional relationship; Figure 5fIt is a functional image of the incident angle and the phase angle of the converted wave transmission coefficient of the soft seabed liquid-solid interface simulated based on the third functional relationship; Figure 5e and Figure 5f Both represent the change of the transmission coefficient phase angle with the incident angle.

[0181] like Figure 5a to Figure 5c As shown, after the incident angle reaches the P-wave critical angle, the reflection coefficient amplitude, longitudinal wave transmission coefficient amplitude and converted wave transmission coefficient amplitude all undergo obvious changes, that is, the reflection / transmission coefficient amplitude all undergo an obvious turning point at the P-wave critical angle.

[0182] like Figure 5c As shown in the figure, when the incident angle is greater than the P-wave critical angle, the phase angle of the converted wave transmission coefficient increases first and then decreases.

[0183] like Figure 5d and Figure 5f As shown, after the incident angle reaches the P-wave critical angle, the reflection coefficient phase angle and the converted wave transmission coefficient phase angle both undergo obvious changes, that is, the reflection / transmission coefficient phase angles both undergo an obvious turning point at the P-wave critical angle.

[0184] like Figure 4e and Figure 5e As shown, the phase angle of the longitudinal wave transmission coefficient of the liquid-solid interface of the hard seabed and the phase angle of the longitudinal wave transmission coefficient of the liquid-solid interface of the soft seabed are always 0 and do not change with the change of the incident angle.

[0185] like Figure 4c and Figure 5c As shown, compared with T in hard seabed media PS The amplitude of the change in the value of the conversion wave transmission coefficient T in the soft seabed medium PS The amplitude value changes slightly.

[0186] like Figure 4c and Figure 5c As shown in the figure, when the incident angle is the critical angle of P wave, the conversion wave transmission coefficient T in the hard seabed medium is PS The amplitude value and the conversion wave transmission coefficient T in the soft seabed medium PS The amplitude values ​​are all the minimum value 0.

[0187] Embodiment 4

[0188] In this embodiment, Figure 2 As shown, a seismic wave field analysis simulation device is provided, comprising:

[0189] A model building module 210 is used to build an analytical equation of the seafloor interface according to the displacement potential function and boundary conditions of the reflection / transmission wave;

[0190] A first functional relationship acquisition module 220 is used to acquire a functional relationship between a displacement potential function of the reflection / transmission wave and a reflection / transmission coefficient, and determine the functional relationship as a first functional relationship;

[0191] A second functional relationship acquisition module 230 acquires the functional relationship between the boundary condition and the influencing factor and determines it as a second functional relationship;

[0192] The data conversion module 240 converts the seabed interface analytical equation by using the first functional relationship and the second functional relationship to obtain a functional relationship between the reflection / transmission coefficient and the influencing factor, and determines the third functional relationship;

[0193] The data acquisition module 250 acquires a plurality of influencing factors, and acquires a reflection / transmission coefficient corresponding to each of the influencing factors based on the third functional relationship;

[0194] The forward modeling module 260 obtains the seafloor seismic wave field data through forward modeling based on the one-to-one correspondence between the influencing factors and the reflection / transmission coefficients.

[0195] In one embodiment, the data acquisition module 250 includes:

[0196] A medium parameter acquisition unit, used to acquire the seabed medium parameter, obtain the functional relationship between the reflection / transmission coefficient and the incident angle based on the third functional relationship, and determine it as a fourth functional relationship;

[0197] The incident angle acquisition unit is used to acquire the incident angle, and based on the fourth functional relationship, acquire the reflection / transmission coefficient corresponding to each incident angle.

[0198] In this embodiment, the influencing factors include the incident angle and seabed medium parameters.

[0199] In one embodiment, the medium parameter acquisition unit includes a wave impedance acquisition subunit and a solid interface transmission angle acquisition subunit.

[0200] In one embodiment, the wave impedance acquisition subunit includes a liquid P-wave wave impedance acquisition subunit, a solid P-wave wave impedance acquisition subunit and a solid S-wave wave impedance acquisition subunit.

[0201] In one embodiment, the solid interface transmission angle acquisition subunit includes a solid interface P-wave transmission angle acquisition subunit and a solid interface S-wave transmission angle acquisition subunit.

[0202] In one embodiment, the transmission coefficient includes a longitudinal wave transmission coefficient and a converted wave transmission coefficient.

[0203] In one embodiment, the third functional relationship is:

[0204]

[0205] In the formula, R PP is the reflection coefficient, T PP is the longitudinal wave transmission coefficient, T PS is the conversion wave transmission coefficient, I1 is the wave impedance of P wave in the seabed liquid medium, I2 is the wave impedance of P wave in the seabed solid medium, and I3 is the wave impedance of S wave in the seabed solid medium. i1 is the incident angle of the downward P wave at the seabed liquid interface, i2 is the transmission angle of the P wave transmitted by the seabed solid interface, and j2 is the transmission angle of the S wave transmitted by the seabed solid interface; G is the seabed analytical variable, and G is related to I1, I2, I3, i1, i2 and j2 respectively.

[0206] In one embodiment, the first functional relationship acquisition module includes:

[0207] A first functional relationship acquisition unit is used to acquire a functional relationship between a displacement potential function of the reflection / transmission wave and a displacement amplitude of the reflection / transmission wave, and determine the relationship as a fifth functional relationship;

[0208] A first functional relationship conversion unit converts the displacement potential function of the reflection / transmission wave in the seabed interface analytical equation into the displacement amplitude of the reflection / transmission wave by using the fifth functional relationship;

[0209] A second functional relationship acquisition unit is configured to acquire a functional relationship between the displacement amplitude of the reflection / transmission wave and the reflection / transmission coefficient, and determine the functional relationship as a sixth functional relationship;

[0210] The second functional relationship conversion unit converts the displacement amplitude of the reflection / transmission wave in the seabed interface analytical equation into the reflection / transmission coefficient using the sixth functional relationship.

[0211] The specific definition of the seismic wave field analysis simulation device can be found in the definition of the seismic wave field analysis simulation method above, which will not be repeated here. Each unit in the above seismic wave field analysis simulation device can be implemented in whole or in part by software, hardware and a combination thereof. The above units can be embedded in or independent of the processor in the computer device in the form of hardware, or can be stored in the memory of the computer device in the form of software, so that the processor can call and execute the operations corresponding to the above units.

[0212] Embodiment 5

[0213] In this embodiment, a computer device is provided. Its internal structure diagram can be shown as follows: Figure 7As shown. The computer device includes a processor, a memory, a network interface, a display screen and an input device connected through a system bus. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program, and the non-volatile storage medium is deployed with a database, which is used to store influencing factors. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The network interface of the computer device is used to communicate with other computer devices that deploy application software. When the computer program is executed by the processor, a method for analyzing and simulating seismic wave fields is implemented. The display screen of the computer device can be a liquid crystal display screen or an electronic ink display screen, and the input device of the computer device can be a touch layer covered on the display screen, or a key, trackball or touchpad set on the computer device housing, or an external keyboard, touchpad or mouse, etc.

[0214] Those skilled in the art will understand that Figure 7 The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than those shown in the figure, or combine certain components, or have a different arrangement of components.

[0215] In one embodiment, a computer device is provided, including a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the following steps are implemented:

[0216] Step 110, constructing the seabed interface analytical equation according to the displacement potential function of the reflection / transmission wave and the boundary conditions.

[0217] Step 120: Obtain a functional relationship between the displacement potential function of the reflection / transmission wave and the reflection / transmission coefficient, and determine it as a first functional relationship.

[0218] Step 130: Obtain a functional relationship between the boundary condition and the influencing factor, and determine it as a second functional relationship.

[0219] Step 140, using the first functional relationship and the second functional relationship, transforming the seabed interface analytical equation to obtain a functional relationship between the reflection / transmission coefficient and the influencing factor, and determining it as a third functional relationship.

[0220] Step 150: Obtain multiple influencing factors, and based on the third functional relationship, obtain the reflection / transmission coefficient corresponding to each of the influencing factors.

[0221] Step 160: Based on the one-to-one correspondence between the influencing factors and the reflection / transmission coefficients, obtain seafloor seismic wave field data through forward modeling.

[0222] In one embodiment, the seabed medium parameters include a wave impedance data set and a solid interface transmission angle data set.

[0223] In one embodiment, the wave impedance data set includes liquid P-wave wave impedance, solid P-wave wave impedance and solid S-wave wave impedance.

[0224] In one embodiment, the solid interface transmission angle data set includes a solid interface P-wave transmission angle and a solid interface S-wave transmission angle.

[0225] In one embodiment, the transmission coefficient includes a longitudinal wave transmission coefficient and a converted wave transmission coefficient.

[0226] In one embodiment, the third functional relationship is the following calculation formula:

[0227]

[0228] In the formula, R PP is the reflection coefficient, T PP is the longitudinal wave transmission coefficient, T PS is the conversion wave transmission coefficient, I1 is the wave impedance of P wave in the seabed liquid medium, I2 is the wave impedance of P wave in the seabed solid medium, and I3 is the wave impedance of S wave in the seabed solid medium. i1 is the incident angle of the downward P wave at the seabed liquid interface, i2 is the transmission angle of the P wave transmitted by the seabed solid interface, and j2 is the transmission angle of the S wave transmitted by the seabed solid interface;

[0229] Among them, G is the seafloor analytical variable, the specific form is

[0230] G=2sin j2sin2j2cos i1(I2cos j2-I3cos i2)-I1cos i2-I2cos i1.

[0231] In this embodiment, I1=ρ1α1, I2=ρ2α2, I3=ρ2β2, where ρ1 is the density of the seabed liquid medium, α1 is the P-wave velocity of seawater; ρ2 is the density of the seabed solid medium, α2 is the P-wave velocity of the seabed solid medium, and β2 is the S-wave velocity of the seabed solid medium.

[0232] In one embodiment, when the processor executes the computer program, the processor further implements the following steps:

[0233] Step 121, obtaining a functional relationship between a displacement potential function of the reflection / transmission wave and a displacement amplitude of the reflection / transmission wave, and determining the fifth functional relationship;

[0234] Step 122, using the fifth functional relationship, converting the displacement potential function of the reflection / transmission wave in the seabed interface analytical equation into the displacement amplitude of the reflection / transmission wave;

[0235] Step 123, obtaining a functional relationship between the displacement amplitude of the reflection / transmission wave and the reflection / transmission coefficient, and determining it as a sixth functional relationship;

[0236] Step 124, using the sixth functional relationship, converting the displacement amplitude of the reflection / transmission wave in the seabed interface analytical equation into the reflection / transmission coefficient.

[0237] Embodiment 6

[0238] A computer-readable storage medium is provided, on which a computer program is stored, and when the computer program is executed by a processor, the following steps are implemented:

[0239] Step 110, constructing the seabed interface analytical equation according to the displacement potential function of the reflection / transmission wave and the boundary conditions.

[0240] Step 120: Obtain a functional relationship between the displacement potential function of the reflection / transmission wave and the reflection / transmission coefficient, and determine it as a first functional relationship.

[0241] Step 130: Obtain a functional relationship between the boundary condition and the influencing factor, and determine it as a second functional relationship.

[0242] Step 140, using the first functional relationship and the second functional relationship, transforming the seabed interface analytical equation to obtain a functional relationship between the reflection / transmission coefficient and the influencing factor, and determining it as a third functional relationship.

[0243] Step 150: Obtain multiple influencing factors, and based on the third functional relationship, obtain the reflection / transmission coefficient corresponding to each of the influencing factors.

[0244] Step 160: Based on the one-to-one correspondence between the influencing factors and the reflection / transmission coefficients, obtain seafloor seismic wave field data through forward modeling.

[0245] In one embodiment, the seabed medium parameters include a wave impedance data set and a solid interface transmission angle data set.

[0246] In one embodiment, the wave impedance data set includes liquid P-wave wave impedance, solid P-wave wave impedance and solid S-wave wave impedance.

[0247] In one embodiment, the solid interface transmission angle data set includes a solid interface P-wave transmission angle and a solid interface S-wave transmission angle.

[0248] In one embodiment, the transmission coefficient includes a longitudinal wave transmission coefficient and a converted wave transmission coefficient.

[0249] In one embodiment, the third functional relationship is the following calculation formula:

[0250]

[0251] In the formula, R PP is the reflection coefficient, T PP is the longitudinal wave transmission coefficient, T PS is the conversion wave transmission coefficient, I1 is the wave impedance of P wave in the seabed liquid medium, I2 is the wave impedance of P wave in the seabed solid medium, and I3 is the wave impedance of S wave in the seabed solid medium. i1 is the incident angle of the downward P wave at the seabed liquid interface, i2 is the transmission angle of the P wave transmitted by the seabed solid interface, and j2 is the transmission angle of the S wave transmitted by the seabed solid interface;

[0252] Among them, G is the seafloor analytical variable, the specific form is

[0253] G=2sin j2sin2j2cos i1(I2cosj2-I3cos i2)-I1cos i2-I2cos i1.

[0254] In this embodiment, I1=ρ1α1, I2=ρ2α2, I3=ρ2β2, where ρ1 is the density of the seabed liquid medium, α1 is the P-wave velocity of seawater; ρ2 is the density of the seabed solid medium, α2 is the P-wave velocity of the seabed solid medium, and β2 is the S-wave velocity of the seabed solid medium.

[0255] In one embodiment, when the computer program is executed by a processor, the following steps are implemented:

[0256] Step 121, obtaining a functional relationship between a displacement potential function of the reflection / transmission wave and a displacement amplitude of the reflection / transmission wave, and determining the fifth functional relationship;

[0257] Step 122, using the fifth functional relationship, converting the displacement potential function of the reflection / transmission wave in the seabed interface analytical equation into the displacement amplitude of the reflection / transmission wave;

[0258] Step 123, obtaining a functional relationship between the displacement amplitude of the reflection / transmission wave and the reflection / transmission coefficient, and determining it as a sixth functional relationship;

[0259] Step 124, using the sixth functional relationship, converting the displacement amplitude of the reflection / transmission wave in the seabed interface analytical equation into the reflection / transmission coefficient.

[0260] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, storage, database or other media used in the embodiments provided in this application can include non-volatile and / or volatile memory. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM) or flash memory. Volatile memory can include random access memory (RAM) or external cache memory. As an illustration and not limitation, RAM is available in many forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).

[0261] The technical features of the above embodiments may be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0262] The above-mentioned embodiments only express several implementation methods of the present application, and the descriptions thereof are relatively specific and detailed, but they cannot be understood as limiting the scope of the invention patent. It should be pointed out that, for a person of ordinary skill in the art, several variations and improvements can be made without departing from the concept of the present application, and these all belong to the protection scope of the present application. Therefore, the protection scope of the patent of the present application shall be subject to the attached claims.

Claims

1. A seismic wave field analytical simulation method, characterized in that: include: According to the displacement potential function and boundary conditions of the reflected / transmitted waves, the analytical equation of the seabed interface is constructed; Obtaining a functional relationship between a displacement potential function of the reflection / transmission wave and a reflection / transmission coefficient, and determining the relationship as a first functional relationship; Obtaining a functional relationship between the boundary condition and the influencing factor, and determining it as a second functional relationship; The first functional relationship and the second functional relationship are used to transform the seabed interface analytical equation to obtain a functional relationship between the reflection / transmission coefficient and the influencing factor, which is determined as a third functional relationship; Acquire multiple influencing factors, and acquire the reflection / transmission coefficient corresponding to each of the influencing factors based on the third functional relationship; Based on the one-to-one correspondence between the influencing factors and the reflection / transmission coefficients, seafloor seismic wave field data is obtained through forward simulation.

2. The method according to claim 1, characterized in that The influencing factors include incident angle and seabed medium parameters; The step of obtaining a plurality of influencing factors and obtaining a reflection / transmission coefficient corresponding to each of the influencing factors based on the third functional relationship comprises: Acquire the seabed medium parameters, and obtain the functional relationship between the reflection / transmission coefficient and the incident angle based on the third functional relationship, and determine it as a fourth functional relationship; The incident angle is obtained, and based on the fourth functional relationship, the reflection / transmission coefficient corresponding to each incident angle is obtained.

3. The method according to claim 2, characterized in that The seabed medium parameters include a wave impedance data set and a solid interface transmission angle data set.

4. The method according to claim 3, characterized in that The wave impedance data set includes liquid P-wave wave impedance, solid P-wave wave impedance and solid S-wave wave impedance; The solid interface transmission angle data group includes a solid interface P-wave transmission angle and a solid interface S-wave transmission angle.

5. The method according to claim 1, characterized in that The transmission coefficient includes a longitudinal wave transmission coefficient and a converted wave transmission coefficient.

6. The method according to claim 5, characterized in that The third functional relationship is: In the formula, R PP is the reflection coefficient, T PP is the longitudinal wave transmission coefficient, T PS is the conversion wave transmission coefficient, I1 is the wave impedance of P wave in the seabed liquid medium, I2 is the wave impedance of P wave in the seabed solid medium, and I3 is the wave impedance of S wave in the seabed solid medium. i1 is the incident angle of the downward P wave at the seabed liquid interface, i2 is the transmission angle of the P wave transmitted by the seabed solid interface, and j2 is the transmission angle of the S wave transmitted by the seabed solid interface; G is the seabed analytical variable, and G is related to I1, I2, I3, i1, i2 and j2 respectively.

7. The method according to claim 1, characterized in that The step of obtaining the functional relationship between the displacement potential function of the reflection / transmission wave and the reflection / transmission coefficient and determining the functional relationship as a first functional relationship comprises: Acquire a functional relationship between a displacement potential function of the reflection / transmission wave and a displacement amplitude of the reflection / transmission wave, and determine it as a fifth functional relationship; Using the fifth functional relationship, the displacement potential function of the reflection / transmission wave in the seabed interface analytical equation is converted into the displacement amplitude of the reflection / transmission wave; Acquire a functional relationship between the displacement amplitude of the reflection / transmission wave and the reflection / transmission coefficient, and determine it as a sixth functional relationship; The sixth functional relationship is used to convert the displacement amplitude of the reflection / transmission wave in the seabed interface analytical equation into the reflection / transmission coefficient.

8. A seismic wave field analysis simulation device, characterized in that: include: Construct a model module and construct the analytical equation of the seabed interface based on the displacement potential function and boundary conditions of the reflection / transmission waves; A first functional relationship acquisition module is used to acquire the functional relationship between the displacement potential function of the reflection / transmission wave and the reflection / transmission coefficient, and determine it as a first functional relationship; A second functional relationship acquisition module is used to acquire the functional relationship between the boundary condition and the influencing factor and determine it as a second functional relationship; A data conversion module, using the first functional relationship and the second functional relationship, transforms the seabed interface analytical equation to obtain a functional relationship between the reflection / transmission coefficient and the influencing factor, and determines it as a third functional relationship; A data acquisition module, which acquires a plurality of influencing factors, and acquires a reflection / transmission coefficient corresponding to each of the influencing factors based on the third functional relationship; The forward simulation module obtains the seafloor seismic wave field data through forward simulation based on the one-to-one correspondence between the influencing factors and the reflection / transmission coefficients.

9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, wherein: When the processor executes the computer program, the steps of the method according to any one of claims 1 to 7 are implemented.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 7 are implemented.

Images