A natural gas hydrate top and bottom boundary AVO forward simulation method, device, equipment and medium

Abstract

The present invention discloses a method, apparatus, device, and medium for AVO forward modeling of natural gas hydrate top and bottom boundaries. The method comprises: obtaining set parameters in response to a first operation on a target object; obtaining drawing parameters input by the target object in response to a second operation on the target object, wherein the drawing parameters include seismic wavelet information; constructing AVO curves for the hydrate top and bottom boundaries based on the set parameters and the drawing parameters in response to a third operation on the target object; and obtaining AVO gathers based on the AVO curves for the hydrate top and bottom boundaries and the drawing parameters in response to a fourth operation on the target object. By invoking and implementing corresponding processing logic in response to relevant operations on the target object, the present invention can rapidly implement AVO feature analysis. Furthermore, unlike conventional AVO analysis that only analyzes the hydrate bottom boundary, embodiments of the present invention analyze both the top and bottom boundaries, resulting in more reliable hydrate evaluation results.

Description

Technical Field

[0001] The present invention relates to the field of data processing technology, and in particular to an AVO forward simulation method, device, equipment and medium for natural gas hydrate top and bottom boundaries. Background Art

[0002] Current AVO forward modeling approaches for hydrate reservoir bottom boundaries primarily focus on studying the hydrate BSR. Traditional methods typically require the use of large-scale commercial software, which is complex and expensive, limiting their usability and popularity. Specifically, existing technologies only perform AVO forward modeling on the BSR of the hydrate reservoir bottom boundary, resulting in low interpretation accuracy. Furthermore, traditional methods typically require the use of large-scale commercial software, which is complex and expensive, limiting their usability and popularity. Summary of the Invention

[0003] The present invention aims to at least partially address the limitations of related technologies. To this end, the present invention provides a method, apparatus, device, and medium for accurately performing AVO forward modeling of the top and bottom boundaries of natural gas hydrates.

[0004] In one aspect, an embodiment of the present invention provides an AVO forward modeling method for natural gas hydrate top and bottom boundaries, comprising the following steps:

[0005] In response to the first operation of the target object, setting parameters are obtained; the setting parameters include all types of petrophysical parameters and thickness data of the hydrate reservoir to be simulated and the upper and lower strata;

[0006] In response to the second operation of the target object, obtaining the drawing parameters input by the target object; the drawing parameters include seismic wavelet information;

[0007] In response to the third operation of the target object, constructing and obtaining the AVO curves of the top and bottom boundaries of the hydrate based on the set parameters and the drawing parameters;

[0008] In response to the fourth operation of the target object, an AVO gather is obtained based on the AVO curves of the top and bottom boundaries of the hydrate and the drawing parameters.

[0009] Optionally, the first operation includes a first input operation and a calculation confirmation operation; and obtaining the setting parameters in response to the first operation of the target object includes the following steps:

[0010] In response to a first input operation of the target object in the first interface region, acquiring thickness data of the target object and a preset number of rock physical parameters;

[0011] In response to a calculation confirmation operation of the target object in the second interface region, based on a preset number of rock physical parameters, a preset derivation relationship between elastic parameters and between elastic parameters and seismic wave parameters is called to complete all rock physical parameters to obtain all categories of rock physical parameters;

[0012] Among them, the full range of rock physical parameters includes compressional wave velocity, shear wave velocity, density, Poisson's ratio, bulk modulus, shear modulus, Young's modulus and Lame constant.

[0013] Optionally, the first interface area is provided with an input box for each parameter item in the set parameters; in response to a first input operation of the target object in the first interface area, obtaining thickness data and a preset number of rock physical parameters input by the target object includes the following steps:

[0014] In response to the target object inputting a first parameter in the thickness data input box, extracting information from the thickness data input box to obtain the thickness data input by the target object;

[0015] In response to the target object's second parameter input operation in the input boxes of multiple rock physical parameters, information is extracted from the input box where the second parameter input operation is performed to obtain a preset number of rock physical parameters input by the target object.

[0016] Optionally, in response to the second operation on the target object, obtaining the drawing parameters input by the target object includes the following steps:

[0017] In response to a second input operation of the target object in the third interface area, obtaining a drawing parameter input by the target object;

[0018] Among them, the drawing parameters include seismic wavelet information, sampling interval information and gain information; the seismic wavelet information includes the seismic wavelet main frequency and the seismic wavelet length; the third interface area is provided with an input box for each information item in the drawing parameters, and the second input operation includes the action of filling in information in the input box of each information item.

[0019] Optionally, in response to the third operation of the target object, constructing an AVO curve of the hydrate top and bottom boundaries based on the set parameters and the drawing parameters includes the following steps:

[0020] In response to the button click operation of the target object in the fourth interface area, the preset operation logic of the Zoplitz equation is called, and then the AVO curves of the top and bottom boundaries of the hydrate are constructed according to the set parameters and drawing parameters.

[0021] Optionally, the Zoplitz equation includes multiple simplified forms, including Aki-Richards, Shuey, and Hilterman, and the fourth interface area is provided with a calculation button for each simplified form; in response to the target object clicking a button in the fourth interface area, calling the preset calculation logic of the Zoplitz equation, and then constructing the AVO curve of the top and bottom boundaries of the hydrate according to the set parameters and the drawing parameters, including the following steps:

[0022] In response to the target object clicking the operation button of the target simplified formula in the fourth interface area, the operation logic of the target simplified formula is called to obtain the curve parameters according to the set parameters and the drawing parameters;

[0023] Among them, the curve parameters include the gradient, intercept and reflection coefficient of the top and bottom boundaries of the hydrate;

[0024] Data visualization is performed based on the numerical variation relationship between curve parameters and drawing parameters, and the AVO curves of the top and bottom boundaries of hydrates are drawn.

[0025] Optionally, in response to the fourth operation of the target object, obtaining an AVO gather based on the AVO curves of the top and bottom boundaries of the hydrate and the drawing parameters includes the following steps:

[0026] In response to a synthetic click operation of the target object in the fifth interface area, a convolution operation is performed on the AVO curves of the top and bottom boundaries of the hydrate and the seismic wavelet information;

[0027] Data visualization is performed based on the results of the convolution operation to draw the forward simulation image of the AVO gather.

[0028] On the other hand, an embodiment of the present invention provides an AVO forward simulation device for natural gas hydrate top and bottom boundaries, comprising:

[0029] The first module is configured to obtain setting parameters in response to a first operation of the target object; the setting parameters include all types of petrophysical parameters and thickness data of the hydrate reservoir to be simulated and its upper and lower strata;

[0030] The second module is configured to obtain drawing parameters input by the target object in response to a second operation of the target object; the drawing parameters include seismic wavelet information;

[0031] A third module is configured to construct an AVO curve of the top and bottom boundaries of the hydrate based on the set parameters and the drawing parameters in response to a third operation of the target object;

[0032] The fourth module is used to obtain an AVO gather based on the AVO curves of the top and bottom boundaries of the hydrate and the drawing parameters in response to the fourth operation of the target object.

[0033] On the other hand, an embodiment of the present invention provides an electronic device, including: a processor and a memory; the memory is used to store programs; the processor executes the program to implement the above-mentioned AVO forward simulation method for the top and bottom boundaries of natural gas hydrates.

[0034] On the other hand, an embodiment of the present invention provides a computer storage medium storing a program executable by a processor. When the program is executed by the processor, it is used to implement the above-mentioned AVO forward simulation method for the top and bottom boundaries of natural gas hydrates.

[0035] In an embodiment of the present invention, in response to a first operation on a target object, set parameters are obtained; the set parameters include all types of rock physical parameters and thickness data for the hydrate reservoir to be simulated and its upper and lower strata. In response to a second operation on the target object, drawing parameters input by the target object are obtained; the drawing parameters include seismic wavelet information. In response to a third operation on the target object, AVO curves for the top and bottom hydrate boundaries are constructed based on the set parameters and the drawing parameters. In response to a fourth operation on the target object, AVO gathers are drawn based on the AVO curves for the top and bottom hydrate boundaries and the drawing parameters. By invoking and implementing corresponding processing logic in response to relevant operations on the target object, the present invention can rapidly implement AVO feature analysis. Furthermore, unlike conventional AVO analysis that only analyzes the bottom hydrate boundary, the embodiment of the present invention analyzes both the top and bottom boundaries, resulting in more reliable hydrate evaluation results. BRIEF DESCRIPTION OF THE DRAWINGS

[0036] The accompanying drawings are used to provide a further understanding of the technical solution of the present invention and constitute a part of the specification. Together with the embodiments of the present invention, they are used to explain the technical solution of the present invention and do not constitute a limitation to the technical solution of the present invention.

[0037] Figure 1 This is a schematic diagram of an implementation environment for the AVO forward modeling method of natural gas hydrate top and bottom boundaries provided by an embodiment of the present invention;

[0038] Figure 2 1 is a flow chart of an AVO forward modeling method for the top and bottom boundaries of natural gas hydrates provided by an embodiment of the present invention;

[0039] Figure 3 A schematic diagram of the expanded flow of step S100 provided in an embodiment of the present invention;

[0040] Figure 4 A schematic diagram of the expanded flow of step S110 provided in an embodiment of the present invention;

[0041] Figure 5 A schematic diagram of a working interface for applying the AVO forward modeling method for the top and bottom boundaries of natural gas hydrates provided in an embodiment of the present invention;

[0042] Figure 6 A schematic diagram of the expansion process for constructing the AVO curves of the top and bottom boundaries of hydrates provided in an embodiment of the present invention;

[0043] Figure 7 A schematic diagram of the expanded flow of step S400 provided in an embodiment of the present invention;

[0044] Figure 8 A schematic diagram of an example of the AVO characteristics of the top and bottom boundaries of a hydrate provided in an embodiment of the present invention;

[0045] Figure 9 A schematic diagram of the structure of an AVO forward simulation device for the top and bottom boundaries of natural gas hydrates provided in an embodiment of the present invention;

[0046] Figure 10 A schematic structural diagram of an electronic device provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0047] In order to make the purpose, technical solutions and advantages of the present invention more clearly understood, the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention.

[0048] It should be noted that although the system diagrams illustrate functional module divisions and the flowcharts illustrate a logical sequence, in certain circumstances, the steps shown or described may be performed in a sequence that differs from the module divisions in the system or the sequence in the flowcharts. The terms "first / S100," "second / S200," and the like in the specification, claims, and drawings are used to distinguish similar objects and are not necessarily intended to describe a specific sequence or precedence.

[0049] References to "embodiments" in this disclosure mean that a particular feature, structure, or characteristic described in connection with the embodiment may be included in at least one embodiment of the disclosure. The appearance of this phrase in various places in the specification does not necessarily refer to the same embodiment, nor does it constitute an independent or alternative embodiment that is mutually exclusive of other embodiments. It is understood, both explicitly and implicitly, by those skilled in the art that the embodiments described herein may be combined with other embodiments.

[0050] In order to facilitate the understanding of the technical solution of the present invention, the proprietary technical terms involved in the technical features of the present invention are first explained:

[0051] BSR (Bottom Simulating Reflector) appears as a distinct reflector on a seismic profile, typically located at the bottom of a hydrate layer, hence the name "bottom-simulating reflector." BSR reflections are caused by the difference in acoustic velocity between the natural gas below the hydrate layer and the hydrate layer itself, resulting in seismic wave reflection at the interface between the two, forming a distinct reflector. BSR is often considered an indicator of hydrate presence in seismic exploration because the presence of a hydrate layer can significantly alter the subsurface formation velocity, creating the BSR reflection signature.

[0052] AVO (Amplitude Variation with Offset) technology is used to study how seismic reflection amplitude varies with the distance between the shot point and the receiver, i.e., the shot offset (or angle of incidence). This allows researchers to explore how the reflection coefficient response varies with offset (or angle of incidence), thereby determining the lithologic characteristics and physical properties of the media overlying and underlying the reflecting interface. AVO analysis enables geophysicists to better assess reservoir rock properties, including porosity, density, lithology, and fluid content. The theoretical basis of AVO is Zoeppritz's equation. Prestack seismic inversion, developed based on AVO theory, can predict formation elastic parameters from observed seismic data, making it an important quantitative seismic data interpretation technique.

[0053] It is understood that the AVO forward modeling method for the top and bottom boundaries of natural gas hydrates provided in the embodiments of the present invention can be applied to any computer device with data processing and computing capabilities, and this computer device can be various terminals or servers. When the computer device in the embodiment is a server, the server is an independent physical server, or a server cluster or distributed system composed of multiple physical servers, or a cloud server that provides basic cloud computing services such as cloud services, cloud databases, cloud computing, cloud functions, cloud storage, network services, cloud communications, middleware services, domain name services, security services, CDN (Content Delivery Network), and big data and artificial intelligence platforms. Optionally, the terminal is a smart phone, tablet computer, laptop computer, desktop computer, etc., but is not limited to this.

[0054] To facilitate understanding of the technical solutions of the present invention, the following are first explained regarding the technical features that may appear in the embodiments of the present invention:

[0055] like Figure 1 FIG. 1 is a schematic diagram of an implementation environment provided by an embodiment of the present invention. Figure 1The implementation environment includes at least one terminal 102 and a server 101. The terminal 102 and the server 101 can be connected to the network in a wireless or wired manner to complete data transmission and exchange.

[0056] Server 101 can be an independent physical server, or a server cluster or distributed system composed of multiple physical servers. It can also be a cloud server that provides basic cloud computing services such as cloud services, cloud databases, cloud computing, cloud functions, cloud storage, network services, cloud communications, middleware services, domain name services, security services, CDN (Content Delivery Network), as well as big data and artificial intelligence platforms.

[0057] In addition, server 101 can also be a node server in a blockchain network. Blockchain is a new application model of computer technologies such as distributed data storage, peer-to-peer transmission, consensus mechanism, and encryption algorithm.

[0058] The terminal 102 may be a smart phone, tablet computer, laptop computer, desktop computer, smart speaker, smart watch, etc., but is not limited thereto. The terminal 102 and the server 101 may be connected directly or indirectly via wired or wireless communication, which is not limited in this embodiment of the present invention.

[0059] Based on the example Figure 1 In the implementation environment shown, an embodiment of the present invention provides a natural gas hydrate top and bottom boundary AVO forward simulation method. The following uses the natural gas hydrate top and bottom boundary AVO forward simulation method applied to the server 101 as an example for description. It is understandable that the natural gas hydrate top and bottom boundary AVO forward simulation method can also be applied to the terminal 102. Specifically, in some specific application scenarios, the embodiment of the present invention can be integrated into office excel software and run on Figure 1 The implementation environment shown is used to realize the AVO forward simulation method of the top and bottom boundaries of natural gas hydrates.

[0060] Reference Figure 2 , Figure 2 The flowchart of the AVO forward simulation method for the top and bottom boundaries of natural gas hydrates provided by the embodiment of the present invention is applied to the server. The execution subject of the AVO forward simulation method for the top and bottom boundaries of natural gas hydrates can be any of the aforementioned computer devices (including servers or terminals). Figure 2 , the method comprises the following steps:

[0061] S100, obtaining setting parameters in response to a first operation of a target object;

[0062] The set parameters include all types of rock physical parameters and thickness data of the hydrate reservoir to be simulated and its upper and lower strata;

[0063] It should be noted that the first operation includes a first input operation and a calculation confirmation operation; in some embodiments, such as Figure 3 As shown, step S100 may include the following steps: S110, in response to a first input operation of the target object in the first interface area, obtaining the thickness data and a preset number of rock physical parameters input by the target object; S120, in response to a calculation confirmation operation of the target object in the second interface area, based on the preset number of rock physical parameters, calling the derivation relationship between the preset elastic parameters and between the elastic parameters and the seismic wave parameters to complete all rock physical parameters and obtain all categories of rock physical parameters; wherein, all categories of rock physical parameters include longitudinal wave velocity, shear wave velocity, density, Poisson's ratio, bulk modulus, shear modulus, Young's modulus and Lame constant.

[0064] In some embodiments, the first interface area is provided with an input box for each parameter item in the setting parameters; Figure 4 As shown, step S110 may include the following steps: S111, in response to the target object's first parameter input operation in the thickness data input box, extracting information from the thickness data input box to obtain the thickness data input by the target object; S112, in response to the target object's second parameter input operation in the input box of multiple rock physical parameters, extracting information from the input box performing the second parameter input operation to obtain a preset number of rock physical parameters input by the target object.

[0065] For example, in some specific implementations, when you open the office excel software, the following will appear: Figure 5 The working interface shown in the figure. Fill in the rock physical parameters of the hydrate reservoir and its upper and lower strata in the area (i.e. the first interface area). You only need to fill in any three (i.e. the preset number, you can choose at least three), such as P-wave velocity, S-wave velocity, and density. Click The "calculated elastic parameters" shown in the second interface area can be used to complete all other rock physical parameters through the derivation relationship between each elastic parameter and between the elastic parameters and seismic wave parameters. At the same time, fill in the thickness of the hydrate reservoir and its upper and lower strata. Among them, Figure 5 middle, The rock physical parameters involved in the region include vp (compressional wave velocity), vs (shear wave velocity), den (ρ, density), poisson_ratio (ν, Poisson's ratio, also expressed as σ), K (bulk modulus), mu (μ, shear modulus, also marked as G), E (Young's modulus) and lamda (λ, Lame constant); in addition, The setting parameters of the region also include thickness (i.e., thickness data of the hydrate reservoir and its upper and lower formations).

[0066] Specifically, the derived relationships between the rock physical parameters are as follows:

[0067] 1. Shear modulus (G, μ) and Young's modulus (E) and Poisson's ratio (ν): G = μ = E / [2(1 + ν)]; The shear modulus directly reflects the material's ability to resist shear deformation.

[0068] 2. Bulk modulus (K) and Young's modulus (E) and Poisson's ratio (ν): K = E / [3(1 - 2ν)]; Bulk modulus reflects the material's ability to resist volume compression.

[0069] 3. Lame constant (λ) and Young's modulus (E) and Poisson's ratio (ν): λ = Eν / [(1 + ν)(1 - 2ν)]; The physical meaning of the Lame constant λ itself is not as intuitive as K and G, but it is very useful in mathematical expressions.

[0070] 4. Poisson's ratio (ν) and bulk modulus (K) and shear modulus (G): ν = (3K - 2G) / [2(3K + G)]; Poisson's ratio is the ratio of the transverse contraction to the longitudinal extension of a material under uniaxial tension.

[0071] 5. Young's modulus (E) and bulk modulus (K) and shear modulus (G): E = 9KG / (3K + G); Young's modulus reflects the material's ability to resist uniaxial tension / compression deformation.

[0072] 5. Young's modulus (E) and Lame constant (λ, μ): E = μ(3λ + 2μ) / (λ + μ);

[0073] 6. Bulk modulus (K) and Lamé constants (λ, μ): K = λ + (2μ) / 3; this is the key formula relating K to λ and μ.

[0074] 7. Poisson's ratio (ν) and Lame constants (λ, μ): ν = λ / [2(λ + μ)]; relationship between wave velocity, elastic constants, and density.

[0075] 8. Longitudinal wave velocity (vp) and bulk modulus (K), shear modulus (G), and density (ρ): vp = √[(K + (4G) / 3) / ρ]. Longitudinal waves (P waves) are compressional waves that cause changes in the volume and shape of the medium as they propagate. Therefore, their velocity depends on K, which resists volume change, and G, which resists shape change. √ indicates the square root of the following content.

[0076] 9. P-wave velocity (vp) vs. Lame constants (λ, μ) and density (ρ): vp = √[(λ + 2μ) / ρ]; This is the most commonly used formula for expressing vp in terms of the Lame constant. From the relationship K = λ + (2μ) / 3, it can be deduced that this is equivalent to Equation 9: λ + 2μ = K + (4μ) / 3.

[0077] 10. Shear wave velocity (vs) and shear modulus (G, μ) and density (ρ): vs = √[G / ρ] = √[μ / ρ]. Shear waves (S waves) are shear waves. Their propagation only causes a change in the shape (shear deformation) of the medium, without causing a change in volume. Therefore, their velocity depends solely on the shear modulus G / μ, which resists shear deformation.

[0078] 11. Shear wave velocity (vs) and Young's modulus (E), Poisson's ratio (ν), and density (ρ): vs = √[E / (2ρ(1 +ν))] = √[G / ρ]; directly substitute into Equation 1 to obtain this.

[0079] Furthermore, there is a corollary and ratio: vp is always greater than vs: because K + 4G / 3 > G (both K and G are positive), so vp > vs. For common Earth rocks, vp / vs is typically between 1.6 and 2.0.

[0080] Relationship between Poisson's ratio ν and vp / vs: Combining equations 9 and 11, and using equation 4, we can obtain: (vp / vs)^2 = [K + (4G) / 3] / G = (K / G) + 4 / 3;

[0081] Substituting or transforming the relationship 4: ν = (3K - 2G) / [2(3K + G)] yields: (vp / vs)^2 = 2(1 - ν) / (1 - 2ν). This formula is crucial in seismology and rock physics, often used to calculate Poisson's ratio ν from observed vp and vs, or conversely, to estimate the wave velocity ratio. When ν approaches 0.5 (incompressible fluid), vp / vs approaches infinity (vs is 0).

[0082] The role of density ρ: The wave speed depends not only on the material's "stiffness" (elastic modulus) but also on its "inertia" (density). The greater the stiffness, the faster the wave speed; the greater the density, the slower the wave speed.

[0083] By configuring the derivation of the above-mentioned various relationship equations and the logic of the mutual relationship between the elastic parameters in the isotropic linear elastic medium in the office Excel software of the embodiment of the present invention, it is possible to complete all other rock physical parameters by using any three rock physical parameters.

[0084] S200, in response to a second operation on the target object, obtaining a drawing parameter input by the target object;

[0085] Among them, the drawing parameters include seismic wavelet information;

[0086] It should be noted that, in some embodiments, step S200 may include the following steps: in response to a second input operation of the target object in the third interface area, obtaining the drawing parameters input by the target object; wherein the drawing parameters include seismic wavelet information, sampling interval information and gain information; the seismic wavelet information includes the seismic wavelet main frequency and the seismic wavelet length; the third interface area is provided with an input box for each information item in the drawing parameters, and the second input operation includes an action of filling in information in the input box of each information item.

[0087] For example, in some specific embodiments, Figure 5 As shown, in The drawing parameters in the third interface area include the main frequency, length, sampling interval, and gain of the seismic wavelet, which are used for subsequent drawing. The drawing parameters also include parameters that help determine the drawing area and format (such as those on the left and below the frame).

[0088] S300, in response to a third operation on the target object, constructing an AVO curve of the top and bottom boundaries of the hydrate based on the set parameters and the drawing parameters;

[0089] It should be noted that, in some embodiments, step S300 may include the following steps: in response to a button click operation of the target object in the fourth interface area, calling the preset operation logic of the Zoplitz equation, and then constructing the AVO curve of the top and bottom boundaries of the hydrate according to the set parameters and drawing parameters.

[0090] In some embodiments, the Zoplitz equation includes multiple simplified forms, including Aki-Richards, Shuey, and Hilterman, and the fourth interface area is provided with operation buttons for each simplified form; Figure 6 As shown, in response to the target object clicking a button in the fourth interface area, the preset Zoplitz equation operation logic is called, and then the hydrate top and bottom boundary AVO curves are constructed according to the set parameters and drawing parameters. The following steps may be included: S311, in response to the target object clicking a calculation button for the target simplified formula in the fourth interface area, the calculation logic of the target simplified formula is called to calculate and obtain curve parameters according to the set parameters and drawing parameters; wherein the curve parameters include the gradient, intercept and reflection coefficient of the hydrate top and bottom boundaries; S312, data visualization processing is performed based on the numerical change relationship between the curve parameters and the drawing parameters, and the hydrate top and bottom boundary AVO curves are drawn.

[0091] The data visualization processing can be achieved by integrating a chart generation script into the office excel software of the embodiment of the present invention.

[0092] For example, in some specific embodiments, The three buttons "Aki-Richards1980", "Shuey1985" and "Hilterman1989" in the area (i.e. the fourth interface area) represent simplified forms of the Zoeppritz equation. Based on these simplified forms, the AVO curves of the top and bottom interfaces of the hydrate can be calculated for the set rock physical parameters.

[0093] In some specific application scenarios, the AVO curves of the hydrate top and bottom boundaries can be constructed based on the set parameters and drawing parameters as follows:

[0094] 1. Interface architecture:

[0095] The fourth interface area: integrates three calculation buttons, corresponding to three simplified forms of the Zoplitz equation:

[0096] Aki-Richards button, Shuey button, Hilterman button;

[0097] Underlying logic: Preset mathematical operation modules for various simplified formulas, dynamically generate AVO curves based on the user-entered parameters (elastic parameters of the hydrate layer / surrounding rock) and drawing parameters (incident angle range, step size).

[0098] 2. Supplementary explanation of the core algorithm principle:

[0099] (1) Aki-Richards simplified formula (1980) principle:

[0100] Based on the linear approximation of the Zoeppritz equation, the reflection coefficient Expressed as:

[0101]

[0102] in, 、 、 They are longitudinal wave velocity, shear wave velocity and density, 、 、 are the relative rates of change of longitudinal wave velocity, shear wave velocity and density respectively; coefficients a, b, c are the incidence angle Function:

[0103]

[0104] (2) Shuey Simplification (1985) Principle:

[0105] The reflection coefficient Reconstructed as intercept (normal incidence reflection coefficient) and gradient :

[0106]

[0107] in:

[0108]

[0109] Physical meaning: Reflects the difference in wave impedance, Reflects the change of Poisson's ratio. The time accuracy is higher.

[0110] (3) Hilterman simplified formula (1989) principle:

[0111] For gas sandstone optimization, an explicit expression of Poisson's ratio σ is introduced:

[0112]

[0113] in, is the value of the reflection coefficient at normal incidence, is the Poisson's ratio difference; characteristics:

[0114] It has outstanding fluid sensitivity (especially for natural gas) and has significant advantages in hydrocarbon detection in the medium angle range (10° to 40°).

[0115] 3. Technical solution implementation process:

[0116] Step S311: curve parameter generation: after the user clicks the target button (such as the Shuey button), the system calls the corresponding simplified calculation logic;

[0117] Step S312: Taking the reflection coefficient as an example, the AVO curve is visualized:

[0118] Input: curve parameters (such as reflection coefficient array R) + drawing parameters (incident angle range 0° to 40°, step size 1°)

[0119] Output: Dynamic drawing Follow The changing curve.

[0120] In actual application scenarios, such as Figure 5 For example, click the "Shuey1985" button to get - The various curves shown. Among them, represents the gradient and intercept, the red diamond represents the gradient and intercept corresponding to the top boundary of the hydrate reservoir, and the black square represents the gradient and intercept corresponding to the bottom boundary of the hydrate reservoir; represents the AVO curve, the red curve represents the AVO curve of the top boundary of the hydrate, and the black curve represents the AVO curve of the bottom boundary of the hydrate reservoir; represents the seismic wavelet, and the main frequency of the wavelet in the figure is 40 Hz; Represents the amplitude spectrum corresponding to the seismic wavelet, the horizontal axis is the frequency, and the vertical axis is the amplitude value corresponding to each frequency.

[0121] S400, in response to a fourth operation on the target object, obtaining an AVO gather based on the AVO curves of the top and bottom boundaries of the hydrate and the drawing parameters;

[0122] It should be noted that, in some embodiments, Figure 7 As shown, step S400 may include the following steps: S410, in response to a synthetic click operation of the target object in the fifth interface area, performing a convolution operation on the AVO curves of the top and bottom boundaries of the hydrate and the seismic wavelet information; S420, performing data visualization processing based on the result of the convolution operation, and drawing a forward simulation image of the AVO gather.

[0123] For example, in some embodiments, after the AVO curve is calculated, click Click the "Synthetic Record" button in the area (i.e. the fifth interface area). The internal program performs convolution operation based on the AVO curve and seismic wavelet, and finally obtains the corresponding AVO forward gather image ( and ),in is the pre-stack forward modeling gather that changes with the incident angle, the horizontal axis is the incident angle, is the gather corresponding to the incident angle of 0.

[0124] Specifically, the convolution operation formula for the AVO curve and seismic wavelet is: x(t) = w(t) * r(t), where x(t) represents the synthesized seismic record (i.e., the AVO curve), w(t) represents the seismic wavelet, and r(t) represents the reflection coefficient calculated using the above formula. The "*" indicates the convolution operation.

[0125] In order to explain the principle of the technical solution of the present invention in detail, the overall process of the present invention is described below in combination with some specific embodiments. It is easy to understand that the following is an explanation of the technical principle of the present invention and cannot be regarded as a limitation of the present invention.

[0126] First, it's important to note that traditional seismic reservoir analysis methods rely on gather stacking, a method with significant limitations. Compared to simply stacking seismic data, AVO technology can provide richer geological information, such as the variation of reflection coefficient with angle. This is crucial for identifying potential oil and gas reservoirs, assessing reserves, and selecting exploration targets. Therefore, AVO technology plays a crucial role in identifying the top and bottom interfaces of hydrates. AVO forward modeling can better understand the AVO characteristics of the top and bottom interfaces of hydrate layers, facilitating the use of prestack gather data to assess hydrate reservoir development.

[0127] In view of this, it is necessary to develop a fast AVO forward modeling device to clarify the AVO characteristics of the top and bottom interfaces of the hydrate reservoir. Using various characteristics as identification charts or standards, the top and bottom interfaces of the hydrate can be quickly and effectively characterized using seismic pre-stack data sets. At the same time, the hydrate saturation value can be quantitatively obtained, providing a data basis for subsequent hydrate resource calculations. To this end, an embodiment of the present invention provides an AVO forward modeling method for the top and bottom boundaries of natural gas hydrates that can be applied to office excel software. Specifically, in some practical application scenarios, the embodiment of the present invention can be integrated into office excel software to implement the following: Figure 5 The operation interface shown in the figure can be implemented as follows:

[0128] 1. Hardware and software preparation: a computer capable of running Office Excel software, and Office Excel software.

[0129] 2. Rock physical parameter setting:

[0130] Open the office excel software and you will see the following Figure 5 The working interface shown in the figure. Fill in the rock physical parameters of the hydrate reservoir and its upper and lower strata in the area. You only need to fill in any three, such as P-wave velocity, S-wave velocity, and density. Click Click the "Calculate Elastic Parameters" button to complete all other petrophysical parameters. At the same time, fill in the thickness of the hydrate reservoir and the strata above and below it.

[0131] 3. Drawing parameter settings:

[0132] exist Fill in the drawing parameters in the area shown, including the main frequency of the seismic wavelet, the length of the seismic wavelet, the sampling interval, the gain and other information, which will be used for subsequent drawing.

[0133] 4. Calculate the AVO curves of the top and bottom interfaces of the hydrate:

[0134] The three buttons in the area, "Aki-Richards1980", "Shuey1985" and "Hilterman1989", represent simplified forms of the Zoeppritz equation. Based on these simplified forms, the AVO curves of the top and bottom interfaces of the hydrate can be calculated for the set rock physical parameters.

[0135] For example, click the "Shuey1985" button to get - The various curves shown. Among them, represents the gradient and intercept, the red diamond represents the gradient and intercept corresponding to the top boundary of the hydrate reservoir, and the black square represents the gradient and intercept corresponding to the bottom boundary of the hydrate reservoir; represents the AVO curve, the red curve represents the AVO curve of the top boundary of the hydrate, and the black curve represents the AVO curve of the bottom boundary of the hydrate reservoir; represents the seismic wavelet, and the main frequency of the wavelet in the figure is 40 Hz; Represents the amplitude spectrum corresponding to the seismic wavelet, the horizontal axis is the frequency, and the vertical axis is the amplitude value corresponding to each frequency.

[0136] 5. AVO gather drawing:

[0137] After calculating the AVO curve, click Click the "Synthetic Record" button shown in the figure. The internal program performs convolution operation based on the AVO curve and seismic wavelet, and finally obtains the corresponding AVO forward gather image ( and ),in is the pre-stack forward modeling gather that changes with the incident angle, the horizontal axis is the incident angle, is the gather corresponding to the incident angle of 0.

[0138] The convolution formula for the AVO curve and seismic wavelet is: x(t) = w(t) * r(t), where x(t) represents the synthesized seismic record, w(t) represents the seismic wavelet, and r(t) represents the reflection coefficient calculated using the above formula. The "*" indicates the convolution operation.

[0139] Using the office excel software of the embodiment of the present invention, the information such as hydrate saturation, hydrate reservoir thickness, wavelet main frequency, free gas concentration, etc. is changed to perform AVO simulation of the top and bottom interfaces of the hydrate reservoir, as shown in FIG. Figure 8 As shown:

[0140] Figure 8In the boxes labeled A1, A2, and A3, the hydrate saturations are 20%, 40%, and 60%, respectively. The hydrate layer thickness is 50 m, the dominant wavelet frequency is 50 Hz, and the free gas concentration is 0. It can be seen that as hydrate saturation increases, the amplitude energy of the waveform at the upper and lower interfaces gradually increases. The AVO characteristic of the upper and lower interfaces is that the amplitude energy increases with increasing incident angle.

[0141] Figure 8 In the boxes labeled B1, B2, and B3, the hydrate layer thicknesses are 5 m, 10 m, and 20 m, respectively. The hydrate saturation is 40%, the wavelet dominant frequency is 50 Hz, and the free gas concentration is 0. It can be seen that as the hydrate layer thickness increases, the top boundary peak and bottom boundary trough gradually separate, and wavelet sidelobes appear between the peaks and troughs. The waveform energy at both the top and bottom boundaries increases with increasing incident angle. Furthermore, due to the tuning effect, the overall amplitude energy of the waveform at the upper and lower interfaces gradually decreases.

[0142] Figure 8 In the boxes labeled C1, C2, and C3, the wavelet dominant frequencies are 30 Hz, 40 Hz, and 50 Hz, respectively. The hydrate layer thickness is 10 m, the hydrate saturation is 40%, and the free gas concentration is 0. It can be seen that as the dominant wavelet frequency increases, the peak and trough events narrow, and the resolution increases.

[0143] Figure 8 In the boxes labeled D1, D2, and D3, the free gas concentrations in the formation below the hydrate are 0, 2%, and 4%, respectively; the hydrate saturation is 40%, the dominant wavelet frequency is 50 Hz, and the hydrate reservoir thickness is 50 m. It can be seen that as the free gas concentration increases, the energy of both the top and bottom wave peaks increases. Furthermore, the AVO characteristic of the upper and lower interfaces is that the amplitude energy increases with increasing incident angle.

[0144] In addition to the aforementioned cases, there are also cases where waveform amplitude decreases with increasing incident angle, which will not be discussed here. The above research clearly demonstrates the impact of changes in information such as hydrate saturation, hydrate reservoir thickness, wavelet dominant frequency, and free gas concentration on hydrate forward modeling gathers. Using these various features as identification charts or standards, seismic prestack gathers can be used to quickly and effectively characterize the top and bottom interfaces of hydrates. Furthermore, quantitative hydrate saturation values ​​can be obtained, providing a data foundation for subsequent hydrate resource calculations.

[0145] In summary, the embodiments of the present invention can be implemented using an AVO forward modeling program written in Microsoft Office Excel. Unlike traditional AVO analysis that only targets the bottom boundary of the hydrate, the present invention analyzes both the top and bottom boundaries. Furthermore, the embodiments of the present invention facilitate AVO simulation by varying information such as hydrate saturation, hydrate reservoir thickness, wavelet dominant frequency, and free gas concentration. Compared to existing technologies, the present invention offers at least the following advantages: Microsoft Office Excel is easy to use and, in terms of operation, allows for rapid AVO feature analysis. Furthermore, unlike traditional AVO analysis that only targets the bottom boundary of the hydrate, the present invention analyzes both the top and bottom boundaries, resulting in more reliable hydrate evaluation results.

[0146] On the other hand, Figure 9 As shown, an embodiment of the present invention provides an AVO forward simulation device 900 for the top and bottom boundaries of natural gas hydrates, which may include:

[0147] The first module 901 is configured to obtain setting parameters in response to a first operation on a target object; the setting parameters include all types of petrophysical parameters and thickness data of the hydrate reservoir to be simulated and its upper and lower strata;

[0148] The second module 902 is configured to obtain drawing parameters input by the target object in response to the second operation of the target object; the drawing parameters include seismic wavelet information;

[0149] The third module 903 is configured to construct an AVO curve of the top and bottom boundaries of the hydrate based on the set parameters and the drawing parameters in response to the third operation of the target object;

[0150] The fourth module 904 is configured to generate an AVO gather based on the AVO curves of the top and bottom boundaries of the hydrate and the drawing parameters in response to the fourth operation of the target object.

[0151] The contents of the method embodiments of the present invention are all applicable to the device embodiments. The functions specifically implemented by the device embodiments are the same as those of the above method embodiments, and the beneficial effects achieved are also the same as those achieved by the above method embodiments.

[0152] In another aspect, an embodiment of the present invention further provides an electronic device comprising a memory and a processor. The memory stores a computer program, and the processor executes the computer program to implement the aforementioned AVO forward modeling method for determining the top and bottom boundaries of natural gas hydrates. The electronic device can be any smart terminal, including a tablet computer and an in-vehicle computer.

[0153] It can be understood that the contents of the above method embodiments are applicable to the present device embodiments, the functions specifically implemented by the present device embodiments are the same as those of the above method embodiments, and the beneficial effects achieved are also the same as those achieved by the above method embodiments.

[0154] like Figure 10 As shown, Figure 10 The hardware structure of an electronic device 1000 according to another embodiment is shown. The electronic device 1000 includes:

[0155] The processor 1001 may be implemented as a general-purpose CPU (Central Processing Unit), a microprocessor, an application-specific integrated circuit (ASIC), or one or more integrated circuits, and is configured to execute relevant programs to implement the technical solutions provided by the embodiments of the present invention.

[0156] The memory 1002 can be implemented in the form of a read-only memory (ROM), a static storage device, a dynamic storage device, or a random access memory (RAM). The memory 1002 can store an operating system and other application programs. When the technical solutions provided in the embodiments of this specification are implemented through software or firmware, the relevant program code is stored in the memory 1002 and is called by the processor 1001 to execute the methods of the embodiments of the present invention.

[0157] Input / output interface 1003, used to implement information input and output;

[0158] Communication interface 1004, used to implement communication interaction between this device and other devices, which can be achieved through wired means (such as USB, network cable, etc.) or wireless means (such as mobile network, WiFi, Bluetooth, etc.);

[0159] Bus 1005 , which transmits information between various components of the device (e.g., processor 1001 , memory 1002 , input / output interface 1003 , and communication interface 1004 );

[0160] The processor 1001 , the memory 1002 , the input / output interface 1003 and the communication interface 1004 are connected to each other in communication within the device via a bus 1005 .

[0161] The electronic device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate, i.e., they may be located in one location or distributed across multiple network units. Some or all of these modules may be selected based on actual needs to achieve the objectives of this embodiment.

[0162] The contents of the method embodiments of the present invention are all applicable to the electronic device embodiments. The functions specifically implemented by the electronic device embodiments are the same as those of the above method embodiments, and the beneficial effects achieved are also the same as those achieved by the above method embodiments.

[0163] Another aspect of an embodiment of the present invention further provides a computer-readable storage medium, wherein the storage medium stores a program, and the program is executed by a processor to implement the above method.

[0164] It should be noted that the computer-readable medium described in the embodiments of the present invention may be a computer-readable signal medium or a computer-readable storage medium, or any combination thereof. A computer-readable storage medium may be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination thereof. More specific examples of computer-readable storage media may include, but are not limited to, an electrical connection having one or more conductors, a portable computer disk, a hard disk, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM), flash memory, optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination thereof. In the present invention, a computer-readable storage medium may be any tangible medium containing or storing a program that can be used by or in conjunction with an instruction execution system, apparatus, or device. In the present invention, a computer-readable signal medium may include a data signal propagated in baseband or as part of a carrier wave, which carries computer-readable program code. Such propagated data signals may take various forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination thereof. A computer-readable signal medium may also be any computer-readable medium other than a computer-readable storage medium that can transmit, propagate, or transport a program for use by or in connection with an instruction execution system, apparatus, or device. Program code embodied on a computer-readable medium may be transmitted using any suitable medium, including but not limited to wireless, wired, or any suitable combination thereof.

[0165] The contents of the method embodiments of the present invention are all applicable to the computer-readable storage medium embodiments. The functions specifically implemented by the computer-readable storage medium embodiments are the same as those of the above method embodiments, and the beneficial effects achieved are also the same as those achieved by the above method embodiments.

[0166] The present invention also discloses a computer program product or computer program, which includes computer instructions stored in a computer-readable storage medium. A processor of a computer device can read the computer instructions from the computer-readable storage medium and execute the computer instructions, causing the computer device to perform the above method.

[0167] The flowcharts and block diagrams in the accompanying drawings illustrate the possible implementation architecture, functions and operations of the systems, methods and computer program products according to various embodiments of the present invention. In this regard, each box in the flowchart or block diagram can represent a module, program segment, or a part of code, and the above-mentioned module, program segment, or a part of code contains one or more executable instructions for implementing the specified logical function. It should also be noted that in some alternative implementations, the functions marked in the box can also occur in an order different from that marked in the accompanying drawings. For example, two boxes represented in succession can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each box in the block diagram or flowchart, and the combination of boxes in the block diagram or flowchart, can be implemented with a dedicated hardware-based system that performs the specified function or operation, or can be implemented with a combination of dedicated hardware and computer instructions.

[0168] It should be noted that although several modules of the device for action execution are mentioned in the above detailed description, this division is not mandatory. In fact, according to an embodiment of the present invention, the features and functions of two or more modules or units described above can be concretized in one module or unit. Conversely, the features and functions of one module or unit described above can be further divided into multiple modules or units to be concretized.

[0169] Through the above description of the embodiments, it will be readily understood by those skilled in the art that the example embodiments described herein can be implemented via software or via a combination of software and necessary hardware. Therefore, the technical solutions according to the embodiments of the present invention can be embodied in the form of a software product, which can be stored on a non-volatile storage medium (such as a CD-ROM, USB flash drive, or mobile hard drive) or on a network and includes instructions for causing a computing device (such as a personal computer, server, touch terminal, or network device) to execute the methods according to the embodiments of the present invention.

[0170] In some optional embodiments, the function / operation mentioned in the block diagram may not occur in the order mentioned in the operation diagram. For example, depending on the function / operation involved, the two boxes shown in succession can actually be executed substantially simultaneously or the boxes can sometimes be executed in reverse order. In addition, the embodiment presented and described in the flow chart of the present invention is provided in an exemplary manner for the purpose of providing a more comprehensive understanding of the technology. The disclosed method is not limited to the operation and logic flow presented by the present invention. Optional embodiments are contemplated in which the order of the various operations is changed and the sub-operations described as a part of a larger operation are performed independently.

[0171] Furthermore, while the present invention has been described in the context of functional modules, it should be understood that, unless otherwise indicated, one or more of the functions and / or features may be integrated into a single physical device and / or software module, or one or more functions and / or features may be implemented in separate physical devices or software modules. It will also be understood that a detailed discussion of the actual implementation of each module is not necessary for understanding the present invention. More specifically, given the properties, functions, and internal relationships of the various functional modules in the devices disclosed in the present invention, the actual implementation of the module will be understood within the ordinary skill of an engineer. Therefore, a person skilled in the art will be able to implement the present invention as set forth in the claims using ordinary skill without undue experimentation. It will also be understood that the specific concepts disclosed are merely illustrative and are not intended to limit the scope of the present invention, which is determined by the full scope of the appended claims and their equivalents.

[0172] If a function is implemented as a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, or the portion that contributes to the prior art, or a portion of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions for causing a computer device (which can be a personal computer, server, or network device, etc.) to execute all or part of the steps of the various embodiments of the method of the present invention. The aforementioned storage medium includes various media that can store program code, such as USB flash drives, mobile hard drives, read-only memories (ROMs), random access memories (RAMs), magnetic disks, or optical disks.

[0173] The logic and / or steps represented in the flowcharts or otherwise described herein, for example, can be considered as a sequenced list of executable instructions for implementing the logical functions, and can be embodied in any computer-readable medium for use by, or in conjunction with, an instruction execution apparatus, device, or apparatus (e.g., a computer-based apparatus, a processor-included apparatus, or other apparatus that can fetch and execute instructions from, an instruction execution apparatus, device, or apparatus). For purposes of this specification, a "computer-readable medium" can be any device that can contain, store, communicate, propagate, or transport a program for use by, or in conjunction with, an instruction execution apparatus, device, or apparatus.

[0174] More specific examples (a non-exhaustive list) of computer-readable media include the following: an electrical connection with one or more wires (electronic devices), a portable computer disk cartridge (magnetic devices), random access memory (RAM), read-only memory (ROM), erasable and programmable read-only memory (EPROM or flash memory), fiber optic devices, and portable compact disk read-only memory (CDROM). In addition, the computer-readable medium may even be paper or other suitable medium on which the program is printed, since the program may be obtained electronically, for example, by optically scanning the paper or other medium and then editing, interpreting, or processing it in another suitable manner as necessary, and then storing it in a computer memory.

[0175] It should be understood that various components of the present invention may be implemented using hardware, software, firmware, or a combination thereof. In the aforementioned embodiments, multiple steps or methods may be implemented using software or firmware stored in a memory and executed by a suitable instruction execution device. For example, if implemented using hardware, as in another embodiment, any one of the following technologies known in the art or a combination thereof may be used: a discrete logic circuit having logic gate circuits for implementing logic functions on data signals, an application-specific integrated circuit having suitable combinational logic gate circuits, a programmable gate array (PGA), a field-programmable gate array (FPGA), etc.

[0176] Throughout this specification, reference to terms such as "one embodiment," "some embodiments," "examples," "specific examples," or "some examples" means that a specific feature, structure, material, or characteristic described in connection with that embodiment or example is included in at least one embodiment or example of the present invention. In this specification, schematic representations of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.

[0177] While embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that various changes, modifications, substitutions, and variations may be made to the embodiments without departing from the principles and spirit of the invention, and that the scope of the invention is defined by the claims and their equivalents.

[0178] The above is a specific description of the preferred implementation of the present invention, but the present invention is not limited to the embodiments. Those skilled in the art can make various equivalent modifications or substitutions without violating the spirit of the present invention. These equivalent modifications or substitutions are all included in the scope defined by the claims of the present invention.

Claims

1. A natural gas hydrate top and bottom boundary AVO forward simulation method, characterized by: The following steps are involved: In response to the first operation of the target object, setting parameters are obtained; the setting parameters include full-category rock physical parameters and thickness data of the hydrate reservoir to be simulated and the strata above and below it; wherein the full-category rock physical parameters include compressional wave velocity, shear wave velocity, density, Poisson's ratio, bulk modulus, shear modulus, Young's modulus, and Lame constant; The first operation includes a first input operation and a calculation confirmation operation; and obtaining the setting parameters in response to the first operation of the target object includes the following steps: In response to the first input operation of the target object in the first interface region, acquiring the thickness data and a preset number of rock physical parameters input by the target object; In response to the calculation confirmation operation of the target object in the second interface region, based on the preset number of rock physical parameters, calling the preset derivation relationships between the elastic parameters and between the elastic parameters and the seismic wave parameters to complete all the rock physical parameters, thereby obtaining the full range of rock physical parameters; In response to the second operation of the target object, obtaining a drawing parameter input by the target object; the drawing parameter includes seismic wavelet information; In response to a third operation on the target object, constructing an AVO curve of the top and bottom boundaries of the hydrate based on the set parameters and the drawing parameters; The step of constructing the AVO curves of the top and bottom boundaries of the hydrate based on the setting parameters and the drawing parameters in response to the third operation of the target object comprises the following steps: In response to a button click operation of the target object in the fourth interface area, calling the preset operation logic of the Zoplitz equation, and then constructing the hydrate top and bottom boundary AVO curve according to the set parameters and the drawing parameters; In response to the fourth operation of the target object, an AVO gather is obtained based on the hydrate top and bottom boundary AVO curves and the drawing parameters.

2. The AVO forward modeling method for natural gas hydrate top and bottom boundaries according to claim 1, characterized in that: The first interface area is provided with an input box for each parameter item in the setting parameters; The step of obtaining the thickness data and a preset number of rock physical parameters input by the target object in response to the first input operation of the target object in the first interface region comprises the following steps: In response to the target object inputting a first parameter into the thickness data input box, extracting information from the thickness data input box to obtain the thickness data input by the target object; In response to the target object's second parameter input operation in the input boxes of the plurality of rock physical parameters, information is extracted from the input box where the second parameter input operation is performed to obtain the preset number of rock physical parameters input by the target object.

3. The AVO forward modeling method for natural gas hydrate top and bottom boundaries according to claim 1, characterized in that: The step of obtaining the drawing parameters input by the target object in response to the second operation of the target object comprises the following steps: In response to a second input operation of the target object in the third interface area, acquiring the drawing parameter input by the target object; Among them, the drawing parameters include the seismic wavelet information, sampling interval information and gain information; the seismic wavelet information includes the seismic wavelet main frequency and the seismic wavelet length; the third interface area is provided with an input box for each information item in the drawing parameters, and the second input operation includes the action of filling in information in the input box of each information item.

4. The AVO forward modeling method for natural gas hydrate top and bottom boundaries according to claim 1, characterized in that: The Zoplitz equation includes multiple simplified formulas, including Aki-Richards, Shuey, and Hilterman. The fourth interface area is provided with a calculation button for each simplified formula. In response to the target object clicking a button in the fourth interface area, the preset Zoplitz equation calculation logic is called, and then the hydrate top and bottom boundary AVO curves are constructed according to the set parameters and the drawing parameters, including the following steps: In response to the target object clicking the operation button of the target simplified formula in the fourth interface area, calling the operation logic of the target simplified formula to obtain curve parameters according to the set parameters and the drawing parameters; The curve parameters include the gradient, intercept and reflection coefficient of the top and bottom boundaries of the hydrate; Data visualization processing is performed based on the numerical variation relationship between the curve parameters and the drawing parameters to draw the AVO curves of the top and bottom boundaries of the hydrate.

5. The AVO forward modeling method for natural gas hydrate top and bottom boundaries according to claim 1, characterized in that: The step of obtaining an AVO gather based on the AVO curves of the top and bottom boundaries of the hydrate and the drawing parameters in response to the fourth operation of the target object includes the following steps: In response to a synthetic click operation of the target object in the fifth interface region, performing a convolution operation on the hydrate top and bottom boundary AVO curves and the seismic wavelet information; Data visualization processing is performed based on the result of the convolution operation to draw a forward simulation image of the AVO gather.

6. A natural gas hydrate top and bottom boundary AVO forward simulation device, characterized by: The AVO forward modeling method for the top and bottom boundaries of natural gas hydrates according to any one of claims 1 to 5, wherein the device comprises: A first module is configured to obtain setting parameters in response to a first operation on a target object; the setting parameters include full-category rock physical parameters and thickness data of a hydrate reservoir to be simulated and its upper and lower strata; wherein the full-category rock physical parameters include compressional wave velocity, shear wave velocity, density, Poisson's ratio, bulk modulus, shear modulus, Young's modulus, and Lame constant; A second module is configured to obtain drawing parameters input by the target object in response to a second operation of the target object; the drawing parameters include seismic wavelet information; A third module is configured to construct an AVO curve of the top and bottom boundaries of the hydrate based on the set parameters and the drawing parameters in response to a third operation of the target object; The fourth module is used to obtain an AVO gather based on the hydrate top and bottom boundary AVO curves and the drawing parameters in response to a fourth operation of the target object.

7. An electronic device, characterized in that: including a processor and a memory; The memory is used to store programs; The processor executes the program to implement the method according to any one of claims 1 to 5.

8. A computer storage medium storing a program executable by a processor, characterized in that: The program executable by the processor is used to implement the method according to any one of claims 1 to 5 when executed by the processor.

Images