Concealed river channel identification method and device, computer device and storage medium

Abstract

The application provides a hidden river channel identification method and device, computer equipment and a storage medium, which comprises the following steps: obtaining well logging data; performing AVO forward simulation based on the well logging data to obtain a simulation result; analyzing the AVO response characteristics of a sandstone section based on the simulation result; selecting data within a first preset incidence angle range for stacking according to the AVO response characteristics of the sandstone section to obtain first seismic data; selecting seismic data within a second preset incidence angle range for stacking to obtain second seismic data; extracting the maximum wave trough attributes of the first seismic data and the second seismic data; comparing the maximum wave trough attributes of the first seismic data and the second seismic data to obtain a comparison result. Based on the AVO forward simulation of a real drilled well, the AVO response characteristics of the sandstone section of the real drilled well are analyzed, the advantage incidence angle information for better depicting the hidden river channel is selected, the pre-stack gather stacking within the advantage incidence angle range is performed, and the purpose of depicting the hidden river channel is achieved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of oil exploration, and in particular to a subtle channel identification method and device, computer equipment and a storage medium. BACKGROUND

[0002] At present, there are many channel delineation methods, but they are all based on the high wave impedance difference between the channel sandstone and the surrounding rock, so that a strong amplitude wave trough reflection is generated in the seismic profile. For the subtle channel, due to the special rock physical property, the wave trough amplitude energy is weak and comparable to the surrounding rock in the seismic profile, so the identification of the subtle channel is extremely difficult, and there is no targeted technology for identifying this type of subtle channel. SUMMARY

[0003] Therefore, it is necessary to provide a subtle channel identification method, device, computer equipment and storage medium aiming at the above technical problems.

[0004] A subtle channel identification method comprises the following steps:

[0005] Obtaining well logging data;

[0006] Based on the well logging data, AVO forward simulation is performed to obtain a simulation result;

[0007] Based on the simulation result, the AVO response characteristics of the sandstone section are analyzed and obtained;

[0008] According to the AVO response characteristics of the sandstone section, seismic data in a first preset incidence angle range is selected for stacking to obtain first seismic data;

[0009] Seismic data in a second preset incidence angle range is selected for stacking to obtain second seismic data, wherein the second preset incidence angle range is different from the first preset incidence angle range;

[0010] The maximum wave trough attribute of the first seismic data and the maximum wave trough attribute of the second seismic data are extracted, and the maximum wave trough attribute of the first seismic data and the maximum wave trough attribute of the second seismic data are compared to obtain a comparison result.

[0011] In one embodiment, the step of obtaining well logging data comprises:

[0012] Obtaining the well logging data and a well logging interpretation result;

[0013] The step of analyzing and obtaining the AVO response characteristics of the sandstone section based on the simulation result comprises:

[0014] Based on the comparison result of the simulation result and the well logging interpretation result, the AVO response characteristics of the sandstone section are analyzed and obtained.

[0015] In one of the embodiments, the logging data includes P-wave velocity, S-wave velocity, density, mineral content, porosity and gas saturation.

[0016] In one of the embodiments, the second preset incidence angle range is larger than the first preset incidence angle range.

[0017] In one of the embodiments, the first preset incidence angle range is 0-20°.

[0018] In one of the embodiments, the second preset incidence angle range is 0-45°.

[0019] A hidden channel identification device, comprising:

[0020] A logging data acquisition module is configured to acquire logging data.

[0021] A forward simulation module is configured to perform AVO forward simulation based on the logging data to obtain a simulation result.

[0022] A response feature acquisition module is configured to analyze and obtain sandstone section AVO response features based on the simulation result.

[0023] A first seismic data obtaining module is configured to select seismic data within a first preset incidence angle range based on the sandstone section AVO response features to perform stacking and obtain first seismic data.

[0024] A second seismic data obtaining module is configured to select seismic data within a second preset incidence angle range to perform stacking and obtain second seismic data, wherein the second preset incidence angle range is different from the first preset incidence angle range.

[0025] A seismic data comparison module is configured to extract maximum trough attributes of the first seismic data and maximum trough attributes of the second seismic data, compare the maximum trough attributes of the first seismic data and the maximum trough attributes of the second seismic data, and obtain a comparison result.

[0026] In one of the embodiments, the logging data acquisition module is further configured to acquire the logging data and logging interpretation results.

[0027] The response feature acquisition module is further configured to analyze and obtain sandstone section AVO response features based on a comparison result of the simulation result and the logging interpretation results.

[0028] A computer device, comprising a memory and a processor, wherein the memory stores a computer program, and the processor implements the following steps when executing the computer program:

[0029] Acquire logging data.

[0030] Based on the logging data, AVO forward simulation is performed to obtain a simulation result;

[0031] Based on the simulation result, a sandstone section AVO response characteristic is analyzed and obtained;

[0032] According to the sandstone section AVO response characteristic, seismic data in a first preset incidence angle range is selected for stacking to obtain first seismic data;

[0033] Seismic data in a second preset incidence angle range is selected for stacking to obtain second seismic data, wherein the second preset incidence angle range is different from the first preset incidence angle range;

[0034] Maximum trough attributes of the first seismic data and maximum trough attributes of the second seismic data are extracted, and the maximum trough attributes of the first seismic data and the maximum trough attributes of the second seismic data are compared to obtain a comparison result.

[0035] A computer readable storage medium has a computer program stored thereon, and the computer program is executed by a processor to implement the following steps:

[0036] Logging data is acquired;

[0037] Based on the logging data, AVO forward simulation is performed to obtain a simulation result;

[0038] Based on the simulation result, a sandstone section AVO response characteristic is analyzed and obtained;

[0039] According to the sandstone section AVO response characteristic, seismic data in a first preset incidence angle range is selected for stacking to obtain first seismic data;

[0040] Seismic data in a second preset incidence angle range is selected for stacking to obtain second seismic data, wherein the second preset incidence angle range is different from the first preset incidence angle range;

[0041] Maximum trough attributes of the first seismic data and maximum trough attributes of the second seismic data are extracted, and the maximum trough attributes of the first seismic data and the maximum trough attributes of the second seismic data are compared to obtain a comparison result.

[0042] The above-mentioned hidden channel identification method, device, computer equipment and storage medium are based on AVO forward simulation of a real drilled well, analyze the AVO response characteristic of the sandstone section of the real drilled well, select the advantage incidence angle information that can better depict the hidden channel, perform pre-stack gather stacking in the advantage incidence angle range, and then achieve the purpose of depicting the hidden channel. BRIEF DESCRIPTION OF DRAWINGS

[0043] Figure 1 a flowchart of a method for identifying a hidden river channel in an embodiment;

[0044] Figure 2 a structural block diagram of an apparatus for identifying a hidden river channel in an embodiment;

[0045] Figure 3 an internal structure diagram of a computer device in an embodiment;

[0046] Figure 4 a flowchart of a method for identifying a hidden river channel in another embodiment;

[0047] Figure 5 a schematic diagram of an AVO forward simulation result of a real drilled well in an embodiment;

[0048] Figure 6 a schematic diagram of a plane attribute based on 0-20° incident angle stacked seismic data in an embodiment;

[0049] Figure 7 a schematic diagram of a plane attribute based on 0-45° incident angle stacked seismic data in an embodiment. DETAILED DESCRIPTION

[0050] 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 do not limit the present application.

[0051] Embodiment One

[0052] In this embodiment, as shown in the accompanying drawings, a method for identifying a hidden river channel is provided, which comprises: Figure 1

[0053] Step 110, acquiring logging data.

[0054] In this step, the logging data of a real drilled well is collected to obtain typical well logging data.

[0055] In an embodiment, the logging data comprises P-wave velocity, S-wave velocity, density, mineral content, porosity and gas saturation.

[0056] Step 120, based on the logging data, performing AVO forward simulation to obtain a simulation result.

[0057] In this embodiment, AVO forward simulation is performed based on known logging data. In an embodiment, AVO forward simulation is performed based on known well P-wave velocity, S-wave velocity, density, porosity and gas saturation to obtain a forward simulation result, i.e. a simulation result. ​

[0058] As Figure 5 shown in the figure, the AVO forward simulation result of the real drilled well, Figure 5 The left side is the AVO forward trace set, and the right side is the AVO response characteristics of the target layer

[0059] Step 130, based on the simulation result, the AVO response characteristics of the sandstone section are analyzed.

[0060] In this embodiment, the simulation result is analyzed to obtain the AVO response characteristics of the sandstone section.

[0061] Specifically, the result of the forward simulation is analyzed, and the AVO response characteristics of the target layer section are typical AVO of the fourth type, that is, the target layer section is a low-impedance gas-bearing sandstone, and when the Poisson's ratio increases from top to bottom, the AVO absolute value shows a decreasing trend. The amplitude characteristics of this type of AVO are: strong trough reflection at zero incidence angle, and the trough energy gradually decreases with the increase of the incidence angle, and the polarity reverses when the incidence angle reaches a certain range. Because this type of AVO has polarity reversal at large incidence angle, when all incidence angle stacking is performed, the peak energy of large incidence angle will offset part of the trough energy, resulting in weakening of the trough energy, and finally the amplitude energy of this type of channel is weak, which is difficult to identify on conventional post-stack data.

[0062] Step 140, according to the AVO response characteristics of the sandstone section, seismic data in a first preset incidence angle range is selected for stacking to obtain first seismic data.

[0063] In this step, according to the law of the AVO response characteristics of the sandstone section, data in the first preset incidence angle range is selected for stacking to obtain partial incidence angle stacking seismic data A, which is the first seismic data, so as to maximize the amplitude energy of the channel sandstone. In one embodiment, the first preset incidence angle range is 0-20°.

[0064] Step 150, seismic data in a second preset incidence angle range is selected for stacking to obtain second seismic data, wherein the second preset incidence angle range is different from the first preset incidence angle range.

[0065] In this step, data in the second preset incidence angle range is stacked to obtain full-stacking seismic data B, which is the second seismic data. In this embodiment, the second preset incidence angle range is greater than the first preset incidence angle range. In one embodiment, the second preset incidence angle range is 0-45°.

[0066] Step 160, the maximum trough attribute of the first seismic data and the maximum trough attribute of the second seismic data are extracted, and the maximum trough attribute of the first seismic data and the maximum trough attribute of the second seismic data are compared to obtain a comparison result.

[0067] In this embodiment, the maximum wave trough attributes extracted from the first seismic data and the second seismic data are obtained.

[0068] As shown in Figure 5 , the maximum wave trough attributes of the first seismic data and the second seismic data along the target layer are extracted, Figure 6 the planar attributes of the first preset incident angle stacked seismic data are obtained, Figure 7 the planar attributes of the second preset incident angle stacked seismic data are obtained. From Figure 6 and Figure 7 , it can be seen that only a weak reflection feature of the channel can be seen in the planar attributes of the second incident angle stacked seismic data, but the development of the channel can be clearly seen in the planar attributes of the first incident angle stacked seismic data, and the channel has drilled sand bodies, verifying the advantages and effectiveness of the present application in depicting the channel.

[0069] In one embodiment, the step of obtaining the logging data comprises: obtaining the logging data and a logging interpretation result; and the step of analyzing the AVO response characteristics of the sandstone section based on the simulation result comprises: analyzing the AVO response characteristics of the sandstone section based on a comparison result of the simulation result and the logging interpretation result.

[0070] In this embodiment, the AVO response characteristics of the sandstone section are obtained by comparing the simulation result with the logging interpretation result of the known well and analyzing based on the comparison result.

[0071] Embodiment Two

[0072] The purpose of the present application is to analyze the AVO characteristics of the drilled well, select the advantageous incident angle that can highlight the energy of the hidden channel for stacking, and then realize the depiction of the hidden channel.

[0073] In this embodiment, a hidden channel identification method based on small incident angle stacking is proposed, which is based on the AVO forward simulation of the drilled well, analyzes the AVO characteristics, selects the advantageous incident angle information that can highlight the hidden channel, and stacks the advantageous incident angle to achieve the purpose of depicting the hidden channel.

[0074] Please refer to Figure 4 , the steps of the hidden channel identification method based on small incident angle stacking are:

[0075] Step 1: Collect the logging data of a typical well, including P-wave velocity, S-wave velocity, density, mineral content, porosity, gas saturation and logging interpretation result, and perform AVO forward simulation based on the P-wave velocity, S-wave velocity, density, porosity and gas saturation of the known well;

[0076] Step 2: Analyze the AVO response characteristics of the sandstone section by comparing the results of the forward simulation with the well logging interpretation results of the known well;

[0077] Step 3: Select data within a certain range of incident angles based on the AVO response characteristics of the sandstone section to obtain partial incident angle stacking seismic data A, and ensure that the amplitude energy of the channel sandstone is maximized;

[0078] Step 4: Stack all the data within the range of incident angles to obtain full stacking seismic data B;

[0079] Step 5: Extract the maximum wave trough attribute along the target layer for seismic data A and B, and compare and analyze the effect of channel depiction.

[0080] The basic principle of the hidden channel recognition method based on small incident angle stacking is as follows: based on the AVO forward simulation of the drilled well, by analyzing the AVO characteristics of the drilled well, the advantage incident angle information that can better depict the hidden channel is selected, and the pre-stack gather within the advantage incident angle range is stacked, so as to achieve the purpose of depicting the hidden channel. Based on the forward simulation of the drilled well, the selected advantage incident angle information is highly targeted, which can effectively depict the specific AVO characteristic channel in the actual work area, and the method is highly targeted.

[0081] In the hidden channel recognition method based on small incident angle stacking:

[0082] According to step 1, the well logging data of the typical well includes P-wave velocity, S-wave velocity, density, mineral content, porosity, gas saturation and well logging interpretation results, and the AVO forward simulation is carried out based on the P-wave velocity, S-wave velocity, density, porosity and gas saturation of the known well, Figure 5 The AVO forward simulation results of the drilled well are shown in the figure, the left side is the AVO forward gather, and the right side is the AVO response characteristics corresponding to the target layer (red solid line).

[0083] According to step 2, the results of the forward simulation are analyzed, and the AVO response characteristics of the target layer section are typical type IV AVO, that is, the target layer section is low impedance gas-bearing sandstone, and the Poisson's ratio increases from top to bottom, and the AVO absolute value shows a decreasing trend. The amplitude characteristics of this type of AVO are strong trough reflection at zero incident angle, and the trough energy gradually decreases with the increase of incident angle, and the polarity reverses when the incident angle reaches a certain range. Since this type of AVO has polarity reversal at large incident angle, when all incident angle stacking is performed, the peak energy of large incident angle will offset part of the trough energy, resulting in weakening of the trough energy, and finally the amplitude energy of this type of channel is weak, which is difficult to identify on the conventional post-stack data.

[0084] Based on step 3, according to the law of AVO response characteristics of sandstone segment, data in the range of 0-20° incident angle is selected to obtain partial incident angle stacking seismic data A to ensure the maximum amplitude energy of the channel sandstone;

[0085] Based on step 4, data in the range of 0-45° incident angle is stacked to obtain full stacking seismic data B;

[0086] Based on step 5, the maximum wave trough attribute of the seismic data A and B along the target layer is extracted, Figure 6 the planar attribute of the 0-20° incident angle stacking seismic data, Figure 7 the planar attribute of the 0-45° incident angle stacking seismic data. From Figure 6 and Figure 7 It can be seen that in the planar attribute of the 0-45° incident angle stacking seismic data, only a weak reflection feature of the channel can be seen, but in the planar attribute of the 0-20° incident angle stacking seismic data, the development of the channel can be clearly seen, and the channel has drilled sand body by the actual drilling, verifying the advantages and effectiveness of the channel depicted by the present application.

[0087] It should be understood that, although each step in the flowchart of Figure 1 is displayed in sequence according to the direction of the arrow, these steps are not necessarily executed in sequence according to the direction of the arrow. 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 may 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.

[0088] Embodiment three

[0089] In this embodiment, as shown in Figure 2 , a hidden channel identification device is provided, comprising:

[0090] The logging data acquisition module 210 is configured to acquire logging data;

[0091] The forward simulation module 220 is configured to perform AVO forward simulation based on the logging data to obtain a simulation result;

[0092] The response feature acquisition module 230 is configured to analyze and obtain AVO response characteristics of a sandstone segment based on the simulation result;

[0093] The first seismic data obtaining module 240 is configured to obtain first seismic data by stacking seismic data in a first preset incident angle range according to the AVO response feature of the sandstone section;

[0094] The second seismic data obtaining module 250 is configured to obtain second seismic data by stacking seismic data in a second preset incident angle range;

[0095] The seismic data comparison module 260 is configured to extract a maximum trough attribute of the first seismic data and a maximum trough attribute of the second seismic data, and compare the maximum trough attribute of the first seismic data with the maximum trough attribute of the second seismic data to obtain a comparison result.

[0096] In an embodiment, the logging data obtaining module is further configured to obtain the logging data and a logging interpretation result.

[0097] The response feature obtaining module is further configured to analyze the AVO response feature of the sandstone section based on the comparison result of the simulation result and the logging interpretation result.

[0098] In an embodiment, the logging data comprises a P-wave velocity, a S-wave velocity, a density, a mineral content, a porosity and a gas saturation.

[0099] In an embodiment, the second preset incident angle range is greater than the first preset incident angle range.

[0100] In an embodiment, the first preset incident angle range is 0-20°.

[0101] In an embodiment, the second preset incident angle range is 0-45°.

[0102] The specific limitations of the hidden channel identification device can refer to the limitations of the hidden channel identification method described above, and will not be repeated here. Each unit in the hidden channel identification device described above can be realized by software, hardware and combinations thereof in whole or in part. The above-mentioned units can be embedded in or independent of the processor in the computer device in hardware form, or can be stored in the memory in the computer device in software form, so as to call and execute the operations corresponding to each unit by the processor.

[0103] Embodiment Four

[0104] In this embodiment, a computer device is provided. Its internal structure diagram can be as follows Figure 4As 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 well logging data. 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 hidden channel identification 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.

[0105] Those skilled in the art can understand that, Figure 3 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.

[0106] In one embodiment, a computer device is provided, comprising a memory and a processor, the memory storing a computer program, and the processor executing the computer program to implement the following steps:

[0107] Obtaining well logging data;

[0108] Based on the well logging data, performing AVO forward simulation to obtain simulation results;

[0109] Based on the simulation results, the AVO response characteristics of the sandstone section are analyzed and obtained;

[0110] According to the AVO response characteristics of the sandstone section, seismic data in a first preset incidence angle range is selected for stacking to obtain first seismic data;

[0111] Seismic data in a second preset incidence angle range is selected for stacking to obtain second seismic data, wherein the second preset incidence angle range is different from the first preset incidence angle range;

[0112] The maximum trough attribute of the first seismic data and the maximum trough attribute of the second seismic data are extracted, and the maximum trough attribute of the first seismic data and the maximum trough attribute of the second seismic data are compared to obtain a comparison result.

[0113] In one embodiment, the processor, when executing the computer program, also implements the following steps:

[0114] obtaining the logging data and logging interpretation results;

[0115] based on the comparison results of the simulation results and the logging interpretation results, obtaining the AVO response characteristics of the sandstone section.

[0116] In one embodiment, the logging data includes P-wave velocity, S-wave velocity, density, mineral content, porosity and gas saturation.

[0117] In one embodiment, the second preset incidence angle range is greater than the first preset incidence angle range.

[0118] In one embodiment, the first preset incidence angle range is 0-20°.

[0119] In one embodiment, the second preset incidence angle range is 0-45°.

[0120] Embodiment five

[0121] In this embodiment, a computer readable storage medium is provided, which stores a computer program. When the computer program is executed by a processor, the following steps are implemented:

[0122] obtaining the logging data and logging interpretation results;

[0123] based on the logging data, performing AVO forward simulation to obtain simulation results;

[0124] based on the simulation results, obtaining the AVO response characteristics of the sandstone section;

[0125] selecting seismic data within a first preset incidence angle range according to the AVO response characteristics of the sandstone section to obtain first seismic data;

[0126] selecting seismic data within a second preset incidence angle range to obtain second seismic data, wherein the second preset incidence angle range is different from the first preset incidence angle range;

[0127] extracting the maximum trough attribute of the first seismic data and the maximum trough attribute of the second seismic data, and comparing the maximum trough attribute of the first seismic data with the maximum trough attribute of the second seismic data to obtain comparison results.

[0128] In one embodiment, the processor, when executing the computer program, also implements the following steps:

[0129] obtaining the logging data and logging interpretation results;

[0130] Based on a comparison result of the simulation result and the well logging interpretation result, a sandstone segment AVO response characteristic is analyzed.

[0131] In one embodiment, the well logging data includes P-wave velocity, S-wave velocity, density, mineral content, porosity, and gas saturation.

[0132] In one embodiment, the second preset incidence angle range is greater than the first preset incidence angle range.

[0133] In one embodiment, the first preset incidence angle range is 0-20°.

[0134] In one embodiment, the second preset incidence angle range is 0-45°.

[0135] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by a computer program instructing related hardware, and the computer program can be stored in a non-volatile computer readable storage medium. When the computer program is executed, it can include the processes of the above-mentioned embodiments. Any reference to memory, storage, database or other medium used 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 various 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).

[0136] The technical features of the above embodiments can be combined in any way. To make the description concise, 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 contradict, they should be considered within the scope of the present disclosure.

[0137] The above-described embodiments are merely illustrative of several embodiments of the present application, which are described in more detail and in a specific manner, but should not be construed as limiting the scope of the patent. It should be noted that for those skilled in the art, several modifications and improvements can be made without departing from the concept of the present application, and these all belong to the protection scope of the present application. Therefore, the protection scope of the patent of the present application should be subject to the appended claims.

Claims

1. A method for identifying a concealed river course, characterized by, The method comprises the following steps: obtaining well logging data; performing AVO forward simulation based on the well logging data to obtain simulation results; analyzing sandstone segment AVO response characteristics based on the simulation results; selecting seismic data within a first preset incidence angle range according to the sandstone segment AVO response characteristics to obtain first seismic data; selecting seismic data within a second preset incidence angle range to obtain second seismic data, wherein the second preset incidence angle range is different from the first preset incidence angle range; extracting the maximum trough attribute of the first seismic data and the maximum trough attribute of the second seismic data, and comparing the maximum trough attribute of the first seismic data with the maximum trough attribute of the second seismic data to obtain a comparison result.

2. The method of claim 1, wherein, The step of obtaining well logging data comprises: obtaining the well logging data and well logging interpretation results; The step of analyzing sandstone segment AVO response characteristics based on the simulation results comprises: analyzing sandstone segment AVO response characteristics based on the comparison result of the simulation results and the well logging interpretation results.

3. The method of claim 1, wherein, The well logging data comprises P-wave velocity, S-wave velocity, density, mineral content, porosity and gas saturation.

4. The method according to any one of claims 1 to 3, characterized in that, The second preset incidence angle range is greater than the first preset incidence angle range.

5. The method of claim 4, wherein, The first preset incidence angle range is 0-20°.

6. The method of claim 4, wherein, The second preset incidence angle range is 0-45°.

7. A concealed watercourse identification device, characterized by, The method comprises the following steps: a well logging data obtaining module for obtaining well logging data; a forward simulation module for performing AVO forward simulation based on the well logging data to obtain simulation results; a response characteristic obtaining module for analyzing sandstone segment AVO response characteristics based on the simulation results; a first seismic data obtaining module for selecting seismic data within a first preset incidence angle range according to the sandstone segment AVO response characteristics to obtain first seismic data; a second seismic data obtaining module for selecting seismic data within a second preset incidence angle range to obtain second seismic data, wherein the second preset incidence angle range is different from the first preset incidence angle range; a seismic data comparison module for extracting the maximum trough attribute of the first seismic data and the maximum trough attribute of the second seismic data, and comparing the maximum trough attribute of the first seismic data with the maximum trough attribute of the second seismic data to obtain a comparison result.

8. The apparatus of claim 7, wherein, The well logging data obtaining module is further used for obtaining the well logging data and well logging interpretation results; The response characteristic obtaining module is further used for analyzing sandstone segment AVO response characteristics based on the comparison result of the simulation results and the well logging interpretation results. 9.A computer device, comprising a memory and a processor, wherein the memory stores a computer program, and the computer device is configured to perform the method according to any one of claims 1-8 when the computer program is executed by the processor. The processor executes the computer program to realize the steps of the method in any one of claims 1 to 6.

10. A computer-readable storage medium having stored thereon a computer program, characterized in that, The computer program is executed by the processor to realize the steps of the method in any one of claims 1 to 6.

Images