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

Abstract

The application provides a seismic wave field analytical simulation method, comprising the following steps: constructing a seabed interface analytical equation according to a displacement potential function of a reflection / transmission wave and a boundary condition; obtaining a first function relationship; obtaining a second function relationship; transforming the seabed interface analytical equation by using the first function relationship and the second function relationship to obtain a function relationship formula of a reflection / transmission coefficient and an influence factor, and determining the function relationship formula as a third function relationship; obtaining a plurality of influence factors, and obtaining the reflection / transmission coefficients corresponding to the influence factors based on the third function relationship; and obtaining seabed seismic wave field data through forward modeling based on the one-to-one corresponding influence factors and reflection / transmission coefficients. The seabed seismic wave field obtained by the above seismic wave field analytical simulation method is highly similar to the real data, and can provide a theoretical model for the prediction and attenuation of seabed multiple waves and ghost waves.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of seismic data processing, and particularly relates to a seismic wave field analytical simulation method and device, computer equipment and a storage medium. BACKGROUND

[0002] In marine seismic exploration, the reflection and transmission of seismic waves at the seabed interface are more complex than those at the solid-solid interface, because the two sides of the seabed interface are liquid medium and solid medium respectively. On the one hand, the seabed interface is a velocity discontinuity surface, and the wave impedance difference between the two sides of the interface is large. Seismic waves behave as acoustic wave propagation in the upper water layer and as solid layer propagation in the lower layer, resulting in complex calculation of reflection and transmission coefficients at the seabed interface. On the other hand, the reflection / transmission coefficients of seabed interfaces with different medium properties vary with the incident angle in different ways. Not only do the reflection / transmission coefficients of hard and soft seabed interfaces have different variation characteristics, but also the size of the longitudinal and transverse wave velocities of the solid layer affects the reflection and transmission coefficient characteristics of the seabed interface.

[0003] Studying the reflection / transmission coefficient characteristics of the seabed can help us better understand the medium properties of the seabed, and is of great significance to the processing and interpretation of seabed multi-component seismic data. Since the shear modulus is zero in the water layer and there is shear wave in the seabed solid layer, the seabed interface can be regarded as a secondary transverse wave source. Using the seabed reflection / transmission coefficient to obtain the transverse wave velocity information of the seabed solid layer is conducive to the study of seabed conversion wave problems. In addition, compared with land, the multiple waves and ghost waves in the water layer are very complex in marine seismic exploration. Studying the reflection / transmission coefficient of the seabed interface can help us 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 coefficients, they cannot analytically discuss the influence of factors on reflection / transmission coefficients, and lack of understanding of the law of the change of reflection / transmission coefficients with the change of influencing factors. SUMMARY

[0005] Therefore, it is necessary to provide a seismic wave field analytical simulation method, device, computer equipment and storage medium to solve the above technical problems.

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

[0007] constructing a seabed interface analytical equation according to the displacement potential function of the reflection / transmission wave and the boundary condition;

[0008] 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;

[0009] obtaining a function relationship between the boundary condition and the influence factor, and determining a second function relationship;

[0010] transforming 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 a third function relationship;

[0011] obtaining a plurality of influence factors, and obtaining the reflection / transmission coefficient corresponding to each of the influence factors based on the third function relationship;

[0012] performing forward modeling to obtain seabed seismic wave field data based on the one-to-one correspondence between the influence factor and the reflection / transmission coefficient.

[0013] In an embodiment, the influence factor includes an incident angle and a seabed medium parameter.

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

[0015] obtaining the seabed medium parameter, and obtaining a function relationship between the reflection / transmission coefficient and the incident angle based on the third function relationship, and determining a fourth function relationship;

[0016] obtaining an incident angle, and obtaining the reflection / transmission coefficient corresponding to each of the incident angles based on the fourth function relationship.

[0017] In an embodiment, the seabed medium parameter includes a wave impedance data set and a solid interface transmission angle data set.

[0018] In an 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 set includes solid interface P-wave transmission angle and solid interface S-wave transmission angle.

[0020] In an embodiment, the transmission coefficient includes P-wave transmission coefficient and converted wave transmission coefficient.

[0021] In an embodiment, the third function relationship is:

[0022]

[0023] In the formula, R PP is the reflection coefficient, T PP is the P-wave transmission coefficient, and T PSFor 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 incidence angle of the downgoing P wave at the seabed liquid interface, i2 is the transmission angle of the transmitted P wave at the seabed solid interface, and j2 is the transmission angle of the transmitted S wave at the seabed solid interface; G is the seabed analytic variable, and G is related to I1, I2, I3, i1, i2 and j2 respectively.

[0024] In one embodiment, the function relationship between the displacement potential function of the reflected / transmitted wave and the displacement amplitude of the reflected / transmitted wave is obtained, and the step of determining the first function relationship comprises:

[0025] The function relationship between the displacement potential function of the reflected / transmitted wave and the displacement amplitude of the reflected / transmitted wave is obtained, and the fifth function relationship is determined;

[0026] The displacement potential function of the reflected / transmitted wave in the seabed interface analytic equation is converted into the displacement amplitude of the reflected / transmitted wave by using the fifth function relationship;

[0027] The function relationship between the displacement amplitude of the reflected / transmitted wave and the reflected / transmitted coefficient is obtained, and the sixth function relationship is determined;

[0028] The displacement amplitude of the reflected / transmitted wave in the seabed interface analytic equation is converted into the reflected / transmitted coefficient by using the sixth function relationship.

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

[0030] The model building module builds the seabed interface analytic equation according to the displacement potential function of the reflected / transmitted wave and the boundary condition;

[0031] The first function relationship acquisition module obtains the function relationship between the displacement potential function of the reflected / transmitted wave and the reflected / transmitted coefficient, and determines the first function relationship;

[0032] The second function relationship acquisition module obtains the function relationship between the boundary condition and the influence factor, and determines the second function relationship;

[0033] The data conversion module converts the seabed interface analytic equation by using the first function relationship and the second function relationship, obtains the function relationship between the reflected / transmitted coefficient and the influence factor, and determines the third function relationship;

[0034] The data acquisition module acquires a plurality of influence factors, and acquires the reflected / transmitted coefficient corresponding to each influence factor based on the third function relationship;

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

[0036] A computer device comprises a memory and a processor, the memory stores a computer program, and the processor implements the following steps when executing the computer program:

[0037] An analytical equation of a seabed interface is constructed according to a displacement potential function of a reflection / transmission wave and a boundary condition;

[0038] A function relationship between the displacement potential function of the reflection / transmission wave and a reflection / transmission coefficient is obtained, and is determined as a first function relationship;

[0039] A function relationship between the boundary condition and an influence factor is obtained, and is determined as a second function relationship;

[0040] The first function relationship and the second function relationship are used to transform the analytical equation of the seabed interface, a function relationship between the reflection / transmission coefficient and the influence factor is obtained, and is determined as a third function relationship;

[0041] A plurality of influence factors are obtained, and a reflection / transmission coefficient corresponding to each of the influence factors is obtained based on the third function relationship;

[0042] Seabed seismic wave field data are obtained through forward simulation based on the one-to-one correspondence between the influence factors and the reflection / transmission coefficients.

[0043] A computer readable storage medium stores a computer program, and the computer program is executed by a processor to implement the following steps:

[0044] An analytical equation of a seabed interface is constructed according to a displacement potential function of a reflection / transmission wave and a boundary condition;

[0045] A function relationship between the displacement potential function of the reflection / transmission wave and a reflection / transmission coefficient is obtained, and is determined as a first function relationship;

[0046] A function relationship between the boundary condition and an influence factor is obtained, and is determined as a second function relationship;

[0047] The first function relationship and the second function relationship are used to transform the analytical equation of the seabed interface, a function relationship between the reflection / transmission coefficient and the influence factor is obtained, and is determined as a third function relationship;

[0048] A plurality of influence factors are obtained, and a reflection / transmission coefficient corresponding to each of the influence factors is obtained based on the third function relationship;

[0049] Based on the one-to-one correspondence between the influence factor and the reflection / transmission coefficient, the seabed seismic wave field data is obtained through forward simulation.

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

[0051] The reflection / transmission wave includes reflection wave and transmission wave.

[0052] The reflection / transmission wave of the present application includes but is not limited to the reflection / transmission wave obtained by the downward P wave incident of the seabed liquid-solid interface.

[0053] The boundary condition can be the boundary condition of the downward P wave incident of the seabed liquid-solid interface.

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

[0055] 1) According to the displacement potential function of the reflection / transmission wave and the boundary condition, the seabed interface analytical equation is constructed, and the seabed interface analytical equation is transformed to construct the function relationship between the reflection / transmission coefficient and the influence factor, so that the change rule of the reflection / transmission coefficient with the change of the influence factor can be found, and the influence of the influence factor on the reflection / transmission coefficient can also be quantitatively discussed.

[0056] 2) Based on the one-to-one correspondence between the influence factor and the reflection / transmission coefficient, the seabed seismic wave field data is obtained through forward simulation. By comparing the seabed seismic wave field data obtained through forward simulation with the actual seabed seismic wave field data, the accuracy of the third function relationship can be verified, and 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 DRAWINGS

[0057] Figure 1 It is a flowchart of the seismic wave field analytical simulation method in an embodiment;

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

[0059] Figure 3 It is a schematic diagram of the boundary condition in an embodiment;

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

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

[0062] Figure 4c is a function image of the incident angle of the hard seafloor liquid-solid interface and the amplitude value of the reflection coefficient simulated based on the third function relationship in an embodiment;

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

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

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

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

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

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

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

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

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

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

[0073] Figure 7 is an internal structure diagram of a computer device in an embodiment. DETAILED DESCRIPTION

[0074] In order to make the purposes, technical solutions and advantages of the present application clearer, the present application will be further described in detail below with reference to the drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and not used to limit the present application.

[0075] Embodiment one

[0076] In this embodiment, as shown in the following formula, a seismic wave field analytical simulation method is provided, which comprises: Figure 1

[0077] In step 110, according to the displacement potential function of the reflected / transmitted wave and the boundary condition, a seabed interface analytical equation is constructed.

[0078] In this embodiment, the boundary condition is a constraint condition of the displacement potential function of the reflected / transmitted wave. In various embodiments, the reflected / transmitted wave comprises a reflected wave and a transmitted wave, and accordingly, the displacement potential function of the reflected / transmitted wave comprises a displacement potential function of the reflected wave and a displacement potential function of the transmitted wave.

[0079] In this embodiment, as shown in the following formula, the horizontal line between fluid and solid represents a seabed liquid-solid interface; fluid represents a seabed liquid medium, i.e., seawater, ρ1 is the density of the seabed liquid medium, and α1 is the P-wave velocity of seawater; solid represents a seabed solid medium, including a hard seabed medium and a soft seabed medium, ρ2 is the density of the seabed solid medium, and α2 and β2 are the P-wave and S-wave velocities of the seabed solid medium, respectively. Figure 3

[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] represents the displacement amplitude of the transmitted S-wave in the solid medium. In this embodiment, the boundary condition is expressed by the following calculation formula:

[0086]

[0087]

[0088] ​​In the formula, φ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, and are all material-related quantities that appear in the strain-stress relationship.

[0089] In the embodiment, the analytic equation about the seabed interface of the reflected wave is constructed according to the displacement potential function of the reflected wave and the boundary condition.

[0090] In the embodiment, the analytic equation about the seabed interface of the transmitted wave is constructed according to the displacement potential function of the transmitted wave and the boundary condition.

[0091] In step 120, the function relationship between the displacement potential function of the reflected / transmitted wave and the reflection / transmission coefficient is obtained, and is determined as the first function relationship.

[0092] In the embodiment, the analytic equation about the seabed interface is constructed based on the displacement potential function of the reflected / transmitted wave. To obtain the function relationship between the reflection / transmission coefficient and the influence factor, the displacement potential function of the reflected / transmitted wave needs to be converted into the reflection / transmission coefficient. Therefore, the function relationship between the displacement potential function of the reflected / transmitted wave and the reflection / transmission coefficient needs to be obtained, and is determined as the first function relationship.

[0093] In the embodiment, the function relationship between the displacement potential function of the reflected wave and the reflection coefficient is obtained, and is determined as the first function relationship about the reflected wave.

[0094] In the embodiment, the function relationship between the displacement potential function of the transmitted wave and the transmission coefficient is obtained, and is determined as the first function relationship about the transmitted wave.

[0095] In step 130, the function relationship between the boundary condition and the influence factor is obtained, and is determined as the second function relationship.

[0096] In the embodiment, the analytic equation about the seabed interface is constructed based on the boundary condition. To obtain the function relationship between the reflection / transmission coefficient and the influence factor, the boundary condition needs to be converted into the influence factor. Therefore, the function relationship between the boundary condition and the influence factor needs to be obtained, and is determined as the second function relationship.

[0097] In step 140, the first function relationship and the second function relationship are used to convert the analytic equation about the seabed interface, to obtain the function relationship between the reflection / transmission coefficient and the influence factor, and is determined as the third function relationship.

[0098] In the embodiment, the first function relationship about the reflected wave and the second function relationship are used to obtain the third function relationship about the reflected wave.

[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: Acquire multiple influencing factors, and based on the third functional relationship, acquire 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 subsequently used as data for 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 factors on the reflection / transmission coefficients 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 wavefield 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 includes:

[0108] Step 151: 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;

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

[0110] In the above embodiments, the seabed medium parameters are set as constants, the function relationship between the reflection / transmission coefficient and the incidence angle is obtained, and then the law that the reflection / transmission coefficient changes with the incidence angle is obtained.

[0111] The seabed medium parameters are different in different sea areas. For example, the seabed medium parameters can be divided into hard seabed medium parameters and soft seabed medium parameters, and the parameters are different between different hard seabed media and different soft seabed media. 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 an embodiment, that is, the seabed medium parameters of the corresponding sea area are obtained in each embodiment, which are substituted into the third function relationship to find the law that the reflection / transmission coefficient changes with the incidence angle in the sea area. The influence of the seabed medium parameters on the reflection / transmission coefficient can also be obtained by comparing at least two embodiments.

[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 function relationship between the reflection / transmission coefficient and the influence factors, the more types of influence factors, the more the actual reflection / transmission coefficient and the influence factors can be simulated, 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 impedance, solid P-wave impedance, and solid S-wave impedance.

[0115] In this embodiment, in the seismic wave field based on the seabed liquid-solid interface, the liquid P-wave impedance, the solid P-wave impedance, and the solid S-wave impedance all affect the reflection / transmission coefficient. Therefore, the function relationship between the reflection / transmission coefficient and the influence factors is constructed based on the liquid P-wave impedance, the solid P-wave impedance, and the solid S-wave impedance, so that the simulated seismic wave field is highly similar to the actual seismic wave field. Theoretical models are provided for the prediction and attenuation of seabed multiple waves and ghost waves, which can be used for the simulation and prediction of seabed interface seismic wave fields, and the accuracy of seismic data denoising is significantly improved.

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

[0117] In the embodiment, in the seismic wave field based on the liquid-solid interface of seabed, the solid interface P wave transmission angle and the solid interface S wave transmission angle both have influence on the reflection / transmission coefficient, therefore, the function relationship between the reflection / transmission coefficient and the influence 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, a theoretical model is provided for the prediction and attenuation of seabed multiple waves and ghost waves, and the model can be used for the simulation and prediction of the seabed interface seismic wave field, and the precision of seismic data denoising is significantly improved.

[0118] The P wave transmission coefficient is an index for evaluating the ability of seismic wave to penetrate strata in seismic exploration, and represents the ratio of the amplitude of the seismic P wave reaching a receiving point after penetrating a certain interface to the amplitude of the incident wave to the interface. The S wave transmission coefficient is the SV wave transmission coefficient, and the SV wave transmission coefficient can be used to evaluate the propagation and conversion of SV wave in seismic exploration. The P wave transmission coefficient and the S wave transmission coefficient have different meanings for seismic wave. In order to improve the fidelity of the simulated seismic wave field, in one embodiment, the transmission coefficient includes the P wave transmission coefficient and the S wave transmission coefficient.

[0119] In the embodiment, the P wave transmission coefficient and the S wave transmission coefficient can make the related data of the seismic wave field more abundant, and effectively improve the accuracy of the seismic wave field obtained through subsequent forward simulation.

[0120] The reflection / transmission coefficient is influenced by many influence factors, and the wide range of influence factors involved in the research of the relationship between the reflection / transmission coefficient and the influence factors 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 function relationship is the following calculation formula:

[0121]

[0122] In the formula, R PP is the reflection coefficient, T PP is the P wave transmission coefficient, T PS is the S 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 downgoing P wave at the seabed liquid interface, i2 is the transmission angle of the transmitted P wave at the seabed solid interface, and j2 is the transmission angle of the transmitted S wave at the seabed solid interface.

[0123] In the formula, G is the seabed analytic variable, and the specific form is

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

[0125] In the embodiment, I1 = p1a1, I2 = p2a2, I3 = p2b2, wherein p1 is the density of the liquid medium in the sea, a1 is the P-wave velocity of the sea water; p2 is the density of the solid medium in the sea, a2 is the P-wave velocity of the solid medium in the sea, and b2 is the S-wave velocity of the solid medium in the sea.

[0126] In one embodiment, the function relationship between the displacement potential function of the reflected / transmitted wave and the reflection / transmission coefficient is obtained, and the step of determining the first function relationship comprises:

[0127] Step 121: obtaining a function relationship between the displacement potential function of the reflected / transmitted wave and the displacement amplitude of the reflected / transmitted wave, and determining the fifth function relationship;

[0128] Step 122: converting the displacement potential function of the reflected / transmitted wave in the sea bottom interface analytical equation into the displacement amplitude of the reflected / transmitted wave by using the fifth function relationship;

[0129] Step 123: obtaining a function relationship between the displacement amplitude of the reflected / transmitted wave and the reflection / transmission coefficient, and determining the sixth function relationship;

[0130] Step 124: converting the displacement amplitude of the reflected / transmitted wave in the sea bottom interface analytical equation into the reflection / transmission coefficient by using the sixth function relationship.

[0131] In the embodiment, the displacement potential function of the reflected / transmitted wave in the sea bottom analytical equation is converted into the reflection / transmission coefficient through the displacement amplitude of the reflected / transmitted wave, and then the function relationship between the reflection / transmission coefficient and the influence factor is obtained to obtain the rule that the reflection / transmission coefficient changes with the influence factor.

[0132] It should be understood that, although Figure 1 the steps in the flowchart of FIG. 1 are shown in sequence according to the direction of the arrows, these steps are not necessarily executed in sequence according to the direction of the arrows. Unless otherwise specified herein, the execution of these steps is not strictly limited in sequence, and these steps can be executed in other sequences. Moreover, Figure 1 at least part of the steps in FIG. 1 can include multiple sub-steps or multiple stages, which are not necessarily executed at the same time, but can be executed at different times, and the execution sequence of these sub-steps or stages is not necessarily sequential, but can be executed in rotation or alternation with other steps or sub-steps or stages of other steps.

[0133] Embodiment Two

[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 functions 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, and λ1, μ1, λ2, and μ2 are all Lame constants, which are material-related quantities that appear in the strain-stress relationship.

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

[0139]

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

[0141] Based on the obtained reflection / transmission coefficient of the liquid-solid interface of the hard seabed, the seabed seismic wavefield 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 solid medium of the hard seabed, the P wave velocity is α2 = 4000 m / s, the S wave velocity is β2 = 1920 m / s, and the density of the solid layer is ρ2 = 2460 kg / m 3 , using the following calculation formula:

[0143]

[0144] Obtaining reflection / transmission coefficients corresponding to the respective influencing factors;

[0145] Based on the one-to-one correspondence between the influencing factors and the reflection / transmission coefficients, seabed 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 seawater, the sound wave velocity in seawater has little interference with the P-wave velocity and 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 exhibit a variation characteristic with two critical transmission angles as the dividing points, namely the critical angle for P-wave transmission through the hard seabed interface and the critical angle for 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, such as 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, so I will not go into details later. Figure 5a As shown in the figure, Phase angle is the phase angle. The same is true in other figures, which will not be described in detail later.

[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 function relationship, Figure 4a It reflects the change of the 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 function 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 function 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 at 1 at incident angles greater than 50°.

[0153] like Figure 4b As shown in FIG, 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 in the figure, the amplitude value of the converted wave transmission coefficient reaches a minimum value when the incidence angle is 22°, and the amplitude value of the converted wave transmission coefficient is 0 when the incidence angle is greater than or equal to 50°.

[0155] As shown in the figure, the amplitude value of the converted wave transmission coefficient reaches a minimum value when the incidence angle is 22°, and the amplitude value of the converted wave transmission coefficient is 0 when the incidence angle is greater than or equal to 50°. Figure 4d As shown in the figure, the phase angle of the reflection coefficient is 0 when the incidence angle is less than 22°, the phase angle of the reflection coefficient generally shows a downward trend when the incidence angle is 22° to 50°, and the downward amplitude is small, the phase angle of the reflection coefficient generally shows a trend of first decreasing and then increasing when the incidence angle is 50° to 90°, and the change trend is large. That is, the phase angle of the reflection coefficient produces obvious turning points when the incidence angle is 22° and 50°.

[0156] As shown in the figure, the amplitude value of the converted wave transmission coefficient reaches a minimum value when the incidence angle is 22°, and the amplitude value of the converted wave transmission coefficient is 0 when the incidence angle is greater than or equal to 50°. Figure 4e As shown in the figure, the amplitude value of the converted wave transmission coefficient reaches a minimum value when the incidence angle is 22°, and the amplitude value of the converted wave transmission coefficient is 0 when the incidence angle is greater than or equal to 50°.

[0157] Figure 4f As shown in the figure, the phase angle of the converted wave transmission coefficient is constant when the incidence angle is 0° to 22°, the amplitude value of the converted wave transmission coefficient first decreases and then increases until 0 when the incidence angle is 22° to 50°, and the amplitude value of the converted wave transmission coefficient is 0 when the incidence angle is greater than or equal to 50°. That is, the phase angle of the converted wave transmission coefficient produces obvious turning points when the incidence angle is 22° and 50°.

[0158] The longitudinal wave transmission coefficient T PP is 0 at the P-wave critical angle, and then the transmitted P-wave becomes a non-uniform surface wave propagating along the seabed interface and the amplitude exponentially decays with depth as the incidence angle continues to increase; as the incidence angle further increases, T PS is 0 at the S-wave critical angle, and then the S-wave becomes a non-uniform surface wave propagating along the seabed interface. The reflection coefficient Rpp and the transmission coefficients T PP and T PS all show a change characteristic with two transmission critical angles as dividing points.

[0159] Figure 6 FIG. 1 is a comparison diagram of the seabed seismic wave field data obtained by forward modeling and the actual seabed seismic wave field data in an embodiment, wherein the left side figure Synthetics is the seabed seismic wave field data obtained by forward modeling, and the right side 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 a horizontal velocity model is adopted, as Figure 6 ​As shown, the trend and amplitude variation characteristics of the synthetic seabed seismic wave field obtained by forward modeling of the plurality of 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 trend and amplitude of different order multiples are similar to the actual data, and the multiples with an incidence angle close to the critical angle of the longitudinal wave or greater than the critical angle of the transverse wave are seismic waves with relatively strong amplitudes. On the other hand, the synthetic seabed seismic wave field record simulates the non-uniform wave of the actual seabed data, which is similar to the refracted wave but has a higher speed.

[0161] The high degree of closeness between the seabed seismic wave field data obtained by forward modeling and the actual seabed seismic wave field data proves that the plurality of data sets obtained by the seismic wave field analytical simulation method can effectively simulate the seabed liquid-solid interface seismic wave field. That is, the denoising accuracy of the seabed seismic wave field data obtained by the above method is good.

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

[0163] Embodiment Three

[0164] In this embodiment, a seismic wave field analytical simulation method based on a 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 seabed interface analytical equation is constructed according to the boundary conditions of the downgoing P-wave incident on the seabed liquid-solid interface and the displacement potential functions on both sides of the seabed liquid-solid interface; wherein the boundary conditions are expressed by the following calculation formula:

[0166]

[0167] In the formula, φ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 appear in the strain-stress relationship.

[0168] The seabed interface analytical equation is solved to obtain an analytical expression of the seabed liquid-solid interface reflection / transmission coefficient; wherein the analytical expression of the seabed liquid-solid interface reflection / transmission coefficient is as follows:

[0169]

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

[0171] The seabed seismic wave field data is obtained by forward modeling according to the obtained soft seabed liquid-solid interface reflection / transmission coefficient.

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

[0173]

[0174] The reflection / transmission coefficients corresponding to each of the influence factors are obtained.

[0175] Based on the one-to-one correspondence between the influence factors and the reflection / transmission coefficients, seabed seismic wave field data are obtained through forward modeling.

[0176] The density of soft seabed solid medium is relatively small, the S-wave velocity propagating in soft seabed solid medium is relatively low, the S-wave velocity at soft seabed liquid-solid interface is generally lower than the sound wave velocity in seawater, the sound wave velocity in seawater has a relatively large interference on the S-wave velocity at soft seabed liquid-solid interface, and the S-wave critical angle does not exist at soft seabed liquid-solid interface; and the P-wave velocity at soft seabed liquid-solid interface is generally higher than the sound wave velocity in seawater, and the sound wave velocity in seawater has a relatively small interference on the P-wave velocity at soft seabed liquid-solid interface, so the variation characteristics of the reflection coefficient and the transmission coefficient with the incident angle at soft seabed liquid-solid interface are mainly controlled by the transmission P-wave critical angle. In the embodiment, the transmission P-wave critical angle at soft seabed interface is different from the transmission P-wave critical angle at hard seabed interface. In other embodiments, the transmission P-wave critical angle at soft seabed interface is the same as the transmission P-wave critical angle at hard seabed interface.

[0177] According to

[0178]

[0179] At the P-wave critical angle, the transmission angle i2 of the transmission P-wave is 90°, and the transmission coefficient is 0. As shown in Figure 5b and Figure 5c , the incident angle at which the transmission coefficient is 0 is the P-wave critical angle. When the incident angle reaches the P-wave critical angle, the transmission P-wave becomes a non-uniform surface wave.

[0180] Specifically, Figure 5d is a function image of the incident angle and the reflection coefficient phase angle of the soft seabed liquid-solid interface simulated based on the third function relationship, Figure 5d indicating the variation of the reflection coefficient phase angle of the soft seabed liquid-solid interface with the incident angle; Figure 5e is a function image of the incident angle and the P-wave transmission coefficient phase angle of the soft seabed liquid-solid interface simulated based on the third function relationship; Figure 5fis a function image of the incident angle of the soft seabed liquid-solid interface based on the third function relationship simulation and the phase angle of the conversion wave transmission coefficient; Figure 5e and Figure 5f both represent the change of the phase angle of the transmission coefficient with the incident angle.

[0181] As shown in Figure 5a to Figure 5c , after the incident angle reaches the P-wave critical angle, the amplitude values of the reflection coefficient, the longitudinal wave transmission coefficient and the conversion wave transmission coefficient all change obviously, that is, the amplitude values of the reflection / transmission coefficients all change obviously at the P-wave critical angle.

[0182] As shown in Figure 5c , when the incident angle is greater than the P-wave critical angle, the phase angle of the conversion wave transmission coefficient presents the characteristics of first increasing and then decreasing.

[0183] As shown in Figure 5d and Figure 5f , after the incident angle reaches the P-wave critical angle, the phase angles of the reflection coefficient and the conversion wave transmission coefficient both change obviously, that is, the phase angles of the reflection / transmission coefficients both change obviously at the P-wave critical angle.

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

[0185] As shown in Figure 4c and Figure 5c , compared with the amplitude value change range of the T PS in the hard seabed medium, the amplitude value change range of the conversion wave transmission coefficient T PS in the soft seabed medium is smaller.

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

[0187] Embodiment Four

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

[0189] The model construction module 210 constructs a seabed interface analytical equation according to the displacement potential function of the reflection / transmission wave and the boundary condition;

[0190] The first function relationship acquisition module 220 acquires a function relationship between a displacement potential function of the reflected / transmitted wave and a reflection / transmission coefficient, and determines the function relationship as a first function relationship.

[0191] The second function relationship acquisition module 230 acquires a function relationship between the boundary condition and the influence factor, and determines the function relationship as a second function relationship.

[0192] The data conversion module 240 converts the sea bottom interface analytical equation by using the first function relationship and the second function relationship, to obtain a function relationship formula between the reflection / transmission coefficient and the influence factor, and determines the function relationship formula as a third function relationship.

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

[0194] The forward simulation module 260 obtains the sea bottom seismic wave field data through forward simulation based on the one-to-one corresponding influence factors and reflection / transmission coefficients.

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

[0196] The medium parameter acquisition unit is configured to acquire the sea bottom medium parameter, obtain a function relationship between the reflection / transmission coefficient and the incidence angle based on the third function relationship, and determine the function relationship as a fourth function relationship.

[0197] The incidence angle acquisition unit is configured to acquire the incidence angle, and acquire a reflection / transmission coefficient corresponding to each of the incidence angles based on the fourth function relationship.

[0198] In the embodiment, the influence factor includes the incidence angle and the sea bottom medium parameter.

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

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

[0201] In an 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 an embodiment, the transmission coefficient includes a P-wave transmission coefficient and a converted wave transmission coefficient.

[0203] In an embodiment, the third function relationship is:

[0204] In an embodiment, the third function relationship is:

[0205] wherein R PP is the reflection coefficient, T PP is the P-wave transmission coefficient, T PS is the converted 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 incidence angle of the down-going P-wave at the seabed liquid interface, i2 is the transmission angle of the transmitted P-wave at the seabed solid interface, and j2 is the transmission angle of the transmitted S-wave at the seabed solid interface; and G is the seabed analytic variable, which is related to I1, I2, I3, i1, i2, and j2, respectively.

[0206] In one embodiment, the first function relationship obtaining module comprises:

[0207] a first function relationship obtaining unit configured to obtain a function relationship between the displacement potential function of the reflected / transmitted wave and the displacement amplitude of the reflected / transmitted wave, and determine the function relationship as a fifth function relationship.

[0208] a first function relationship converting unit configured to convert, by using the fifth function relationship, the displacement potential function of the reflected / transmitted wave in the seabed interface analytic equation into the displacement amplitude of the reflected / transmitted wave.

[0209] a second function relationship obtaining unit configured to obtain a function relationship between the displacement amplitude of the reflected / transmitted wave and the reflection / transmission coefficient, and determine the function relationship as a sixth function relationship.

[0210] a second function relationship converting unit configured to convert, by using the sixth function relationship, the displacement amplitude of the reflected / transmitted wave in the seabed interface analytic equation into the reflection / transmission coefficient.

[0211] The specific limitations of the seismic wave field analytic simulation device can refer to the limitations of the seismic wave field analytic simulation method in the foregoing, and will not be described herein. Each unit in the above seismic wave field analytic simulation device can be realized by software, hardware, and combinations thereof, in whole or in part. 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 in the computer device in the form of software, so as to be called and executed by the processor to perform the operations corresponding to each unit.

[0212] Embodiment Five

[0213] In this embodiment, a computer device is provided. Its internal structure diagram can be as shown in Figure 7As shown in the figure. 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, 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 for storing influence 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 deployed with application software. The computer program is executed by the processor to implement a seismic wave field analytical simulation method. 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 overlaid on the display screen, or a key, trackball or touchpad arranged on the shell of the computer device, or an external keyboard, touchpad or mouse, etc.

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

[0215] In one embodiment, a computer device is provided, comprising a memory and a processor, the memory stores a computer program, and the processor executes the computer program to implement the following steps:

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

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

[0218] Step 130, obtaining a functional relationship between the boundary condition and the influence factor, and determining it as a second functional relationship.

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

[0220] Step 150, obtaining a plurality of influence factors, and based on the third functional relationship, obtaining a reflection / transmission coefficient corresponding to each of the influence factors.

[0221] Step 160, based on the one-to-one correspondence between the influence factor and the reflection / transmission coefficient, the seabed seismic wave field data is obtained by forward simulation.

[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 solid interface P-wave transmission angle and 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 converted 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, I3 is the wave impedance of S wave in the seabed solid medium, i1 is the incidence angle of downgoing P wave at the seabed liquid interface, i2 is the transmission angle of transmitted P wave at the seabed solid interface, and j2 is the transmission angle of transmitted S wave at the seabed solid interface.

[0229] Wherein, G is the seabed analytic variable, and 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, wherein ρ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 following steps are also implemented:

[0233] Step 121, the function relationship between the displacement potential function of the reflected / transmitted wave and the displacement amplitude of the reflected / transmitted wave is obtained, and the fifth function relationship is determined.

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

[0235] Step 123, obtaining the function relationship between the displacement amplitude of the reflected / transmitted wave and the reflection / transmission coefficient, and determining the sixth function relationship;

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

[0237] Embodiment six

[0238] A computer readable storage medium is provided, and a computer program is stored on the computer readable storage medium. The computer program is executed by a processor to implement the following steps:

[0239] Step 110, constructing a seabed interface analytical equation according to the displacement potential function of the reflected / transmitted wave and the boundary condition.

[0240] Step 120, obtaining the function relationship between the displacement potential function of the reflected / transmitted wave and the reflection / transmission coefficient, and determining the first function relationship.

[0241] Step 130, obtaining the function relationship between the boundary condition and the influence factor, and determining the second function relationship.

[0242] Step 140, converting the seabed interface analytical equation by using the first function relationship and the second function relationship, obtaining the function relationship between the reflection / transmission coefficient and the influence factor, and determining the third function relationship.

[0243] Step 150, obtaining a plurality of influence factors, and obtaining the reflection / transmission coefficient corresponding to each influence factor based on the third function relationship.

[0244] Step 160, obtaining seabed seismic wave field data by forward modeling based on the one-to-one corresponding influence factor and reflection / transmission coefficient.

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

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

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

[0248] In one embodiment, the transmission coefficients include 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] wherein R PP is a reflection coefficient, T PP is a longitudinal wave transmission coefficient, T PS is a converted wave transmission coefficient, I1 is the wave impedance of a P wave in a seabed liquid medium, I2 is the wave impedance of a P wave in a seabed solid medium, and I3 is the wave impedance of an S wave in the seabed solid medium; i1 is the incidence angle of a downward P wave at a seabed liquid interface, i2 is the transmission angle of a transmitted P wave at a seabed solid interface, and j2 is the transmission angle of a transmitted S wave at the seabed solid interface.

[0252] wherein G is a seabed analytic variable, and the specific form is

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

[0254] In the present embodiment, I1 = p1a1, I2 = p2a2, and I3 = p2b2, wherein p1 is the density of the seabed liquid medium, a1 is the P wave velocity of seawater; p2 is the density of the seabed solid medium, a2 is the P wave velocity of the seabed solid medium, and b2 is the S wave velocity of the seabed solid medium.

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

[0256] Step 121: obtaining a functional relationship between the displacement potential function of the reflected / transmitted wave and the displacement amplitude of the reflected / transmitted wave, and determining the functional relationship as a fifth functional relationship;

[0257] Step 122: using the fifth functional relationship to convert the displacement potential function of the reflected / transmitted wave in the seabed interface analytic equation into the displacement amplitude of the reflected / transmitted wave;

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

[0259] Step 124: using the sixth functional relationship to convert the displacement amplitude of the reflected / transmitted wave in the seabed interface analytic 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. The computer program can be stored in a non-volatile computer readable storage medium, and when the computer program is executed, the processes of the above-mentioned embodiments of the methods can be included. Any reference to memory, storage, databases, or other media in the embodiments provided by the present 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 but not limitation, RAM is available in many forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDR SDRAM), 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), etc.

[0261] The technical features of the above embodiments can be combined in any way. In order to make the description simple, not all possible combinations of the technical features in the above embodiments are described, but as long as the combinations of the technical features do not exist, they should be considered as the scope of the present application.

[0262] The above embodiments only express several implementation manners of the present application, and the description is more specific and detailed, but it should not be understood as a limitation on the scope of the patent. It should be pointed out that for ordinary skilled in the art, without departing from the concept of the present application, a number of modifications and improvements can be made, which are all within the scope of the present application. Therefore, the scope of the patent protection of the present application should be subject to the appended 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 a 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, seabed 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 includes: Acquiring the seabed medium parameters, and obtaining a functional relationship between the reflection / transmission coefficient and the incident angle based on the third functional relationship, and determining the 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, wherein 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: Where 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 waves in the seabed liquid medium, I2 is the wave impedance of P waves in the seabed solid medium, and I3 is the wave impedance of S waves 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, which 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 a 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 includes: 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 relationship 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; Obtaining a functional relationship between the displacement amplitude of the reflection / transmission wave and the reflection / transmission coefficient, and determining the relationship 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, which acquires a functional relationship between the displacement potential function of the reflection / transmission wave and the reflection / transmission coefficient, and determines the functional relationship as a first functional relationship; A second functional relationship acquisition module is configured to acquire a functional relationship between the boundary condition and the influencing factor and determine the functional relationship as a second functional relationship; A data conversion module is configured to transform the seabed interface analytical equation 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 determine the functional relationship as a third functional relationship; A data acquisition module is configured to acquire a plurality of influencing factors and, based on the third functional relationship, acquire a reflection / transmission coefficient corresponding to each of the influencing factors; The forward simulation module obtains seabed 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