A method, device, medium and equipment for rapid detection of reservoir gas cap risk

Abstract

The present invention relates to a method, device, medium, and equipment for rapid detection of reservoir gas cap risk. The method comprises the following steps: determining the seismic response characteristics of the gas layer through rock physics analysis; determining whether the amplitude attribute anomaly of the reservoir seismic data is consistent with the seismic response characteristics of the gas layer; extracting the superimposed amplitude values ​​of contour lines in different time domains; delineating the gas-bearing boundaries of the depth domain contour lines to determine the gas-bearing area, and obtaining the gas cap value by calculation. The present invention obtains the superimposed amplitude values ​​of the contour lines in different time domains through a three-step method, so as to facilitate conversion to the depth domain to delineate the gas-bearing boundaries and calculate the gas cap value, thereby enabling rapid prediction of the reservoir gas cap risk based on the size of the gas cap value.

Description

Technical Field

[0001] The present invention relates to the technical field of oil and gas field development, and in particular to a method, device, medium and equipment for rapid detection of reservoir gas cap risk. Background Art

[0002] Fluid testing is a further advancement of reservoir prediction during the oil and gas field development phase. The hydrocarbon content of the reservoir has a significant impact on reservoir development methods and well pattern design. This is especially true when a gas cap is present in the reservoir. The gas cap, which represents the geological reserves of the gas layer, and its size directly determine the oilfield development approach. Currently, fluid testing is primarily based on prestack seismic inversion, which requires a high level of seismic and well logging data. The prediction process is complex and time-consuming, and it is only applicable to specific types of oil and gas reservoirs. This restricts the widespread application of fluid prediction technology during the development phase, making it difficult for geophysical fluid testing to fully guide oilfield development methods and well pattern design. Summary of the Invention

[0003] The purpose of the present invention is to provide a method, device, medium and equipment for rapid detection of reservoir gas cap risk, so as to solve the problem that the current fluid detection method in the prior art is complex, time-consuming and has limited applicability in the process of predicting reservoir gas cap risk.

[0004] To achieve the above object, the present invention adopts the following technical solutions:

[0005] The present invention provides a method for rapid detection of reservoir gas cap risk, comprising:

[0006] Based on post-stack seismic data, the seismic response characteristics of the gas layer are clarified through rock physical analysis;

[0007] Interpret reservoir layers and extract the amplitude of reservoir seismic data to determine if the amplitude attribute anomalies of reservoir seismic data are consistent with the seismic response characteristics of the gas layer;

[0008] Seismic flat point amplitude enhancement is used to stack the amplitude of reservoir seismic data along the interpreted time domain contour lines of the reservoir plane, and the stacked amplitude values ​​of different time domain contour lines are extracted;

[0009] The extracted superimposed amplitude values ​​are converted into depth values ​​corresponding to the contour lines in different time domains. The gas-bearing area is determined by delineating the gas-bearing boundaries of the depth domain contour lines, and the gas cap value is obtained by calculation.

[0010] Furthermore, the method also includes a method for stacking the amplitude of reservoir seismic data along the interpreted reservoir layer time domain contour line using seismic flat point amplitude enhancement:

[0011] Based on post-stack pure wave seismic data, the top surface of the reservoir time domain is traced and interpreted, and the top surface time grid of the reservoir time domain is constructed;

[0012] The time domain contour lines of the reservoir top surface time grid are used as superposition units, and the amplitude values ​​of the pure wave seismic data samples corresponding to the same time domain contour lines are superimposed;

[0013] The transformation relationship between the superimposed amplitude value and the time domain contour line is constructed, and the superimposed amplitude values ​​corresponding to different time domain contour lines are counted.

[0014] Furthermore, it also includes a method for constructing a time grid on the top surface of the reservoir time domain: tracking seismic peaks, troughs or zero-value points during the interpretation of the top surface layer of the reservoir time domain, and generating a time grid on the top surface of the reservoir with preset grid parameters.

[0015] Furthermore, the method also includes a method for superimposing the amplitude values ​​of pure wave seismic data samples corresponding to the same time domain contour lines: superimposing the amplitude values ​​of pure wave seismic data samples corresponding to the same time domain contour lines using an amplitude superposition formula, the amplitude superposition formula being:

[0016] A(t0)=∑ x,y D(x,y,t=t0) / N x,y ;

[0017] Where A(t0) is the superimposed amplitude value of the time domain contour line value t0, D(x,y,t=t0) is the amplitude value of the three-dimensional pure wave seismic data at the sample point position, x is the main survey line direction, y is the connecting survey line direction, t is the time direction, N x,y is the total number of sample points involved in the accumulation.

[0018] Furthermore, the method also includes constructing a transformation relationship between the superimposed amplitude value and the time domain contour line, and counting the superimposed amplitude values ​​corresponding to different time domain contour lines:

[0019] Based on the obtained superimposed amplitude values, a structural map showing the transformation of the superimposed amplitude values ​​with the time domain contour lines at the reservoir level is constructed, with the time domain contour lines as the horizontal axis and the superimposed amplitude values ​​as the vertical axis. The superimposed amplitude values ​​corresponding to different time domain contour lines are statistically analyzed based on the order of the time domain contour lines from shallow to deep, and the time domain contour line corresponding to the amplitude mutation position is predicted to be the gas-bearing boundary.

[0020] Furthermore, it also includes a method for delineating the gas-bearing boundary in the reservoir top depth domain:

[0021] Based on the reservoir top surface structural map, the extracted stacking amplitude values ​​are assigned to all sample points corresponding to the corresponding depth domain contour lines;

[0022] Interpolate the stacked amplitude values ​​of the contour lines in each depth domain to obtain the stacked amplitude value plane attributes corresponding to the reservoir top surface structure map;

[0023] According to the boundary between gas layer and non-gas layer in the plane attributes of superimposed amplitude value, the gas-bearing boundary of the depth domain contour line is delineated.

[0024] Furthermore, the area above the boundary is the superimposed amplitude value characteristic of the reaction gas layer, and the area below the boundary is the superimposed amplitude value characteristic reflecting the non-gas layer.

[0025] Based on the above-mentioned method for rapid detection of reservoir gas cap risk, the present invention further provides an analysis device, comprising:

[0026] The first processing unit is used to determine the seismic response characteristics of the gas layer through rock physical analysis based on the post-stack seismic data;

[0027] The second processing unit is used to interpret the reservoir layer and extract the amplitude of the reservoir seismic data, and determine whether the amplitude attribute anomaly of the reservoir seismic data is consistent with the seismic response characteristics of the gas layer;

[0028] The third processing unit is used to superimpose the amplitude of the reservoir seismic data along the interpreted time domain contour line of the reservoir layer using seismic flat point amplitude enhancement, and extract the superimposed amplitude values ​​of different time domain contour lines;

[0029] The fourth processing unit is used to convert the extracted superimposed amplitude values ​​into depth domain contour values ​​corresponding to different time domain contour lines through time-depth conversion, determine the gas-bearing area by delineating the gas-bearing boundaries of the depth domain contour lines, and obtain the gas cap value by calculation.

[0030] Based on the above-mentioned method for rapid detection of reservoir gas cap risk, the present invention also provides a computer-readable storage medium having a computer program stored thereon. When the computer program is executed by a processor, the steps of the above-mentioned method for rapid detection of reservoir gas cap risk are implemented.

[0031] Based on the above-mentioned method for rapid detection of reservoir gas cap risks, the present invention also provides a computer device, including a memory, a processor, and a computer program stored in the memory and runnable on the processor. When the processor executes the computer program, the steps of the above-mentioned method for rapid detection of reservoir gas cap risks are implemented.

[0032] The present invention adopts the above technical solution, which has the following beneficial effects:

[0033] To achieve the above object, the present invention adopts the following technical solutions:

[0034] A three-step method is used to clarify the seismic response characteristics of the gas layer - determine whether the amplitude attribute anomaly is consistent with the seismic response characteristics of the gas layer - and extract the superimposed amplitude values ​​of the contour lines in different time domains. The extracted superimposed amplitude values ​​are converted into depth by time-to-depth conversion to obtain the depth domain contour line values ​​corresponding to the contour lines in different time domains. The gas-bearing boundaries of the depth domain contour lines are circled to determine the gas-bearing area, so that the gas cap value, that is, the geological reserves of the gas layer, can be obtained through corresponding calculations. Therefore, the gas cap risk of the reservoir can be quickly predicted based on the size of the gas cap value. The process is simple and can effectively play the guiding role of geophysical fluid detection in oilfield development methods and well pattern design. BRIEF DESCRIPTION OF THE DRAWINGS

[0035] Various other advantages and benefits will become apparent to those skilled in the art upon reading the detailed description of the preferred embodiment below. The accompanying drawings are for illustration purposes only and are not to be considered as limiting the present invention. Throughout the drawings, the same reference numerals are used to denote the same components. In the drawings:

[0036] Figure 1 This is a schematic flow chart of a method for rapid detection of reservoir gas cap risk provided by an embodiment of the present invention;

[0037] Figure 2 This is a schematic diagram of a petrophysical analysis histogram of a reservoir in the Bohai A oilfield according to a method for rapid detection of reservoir gas cap risk provided by an embodiment of the present invention;

[0038] Figure 3 Schematic diagram of amplitude attributes of seismic data of the Bohai A oilfield reservoir in a method for rapid detection of reservoir gas cap risk provided by an embodiment of the present invention;

[0039] Figure 4 This is a schematic diagram of the superposition of seismic data amplitude along the time domain contour lines in a method for rapid detection of reservoir gas cap risk provided by an embodiment of the present invention;

[0040] Figure 5 This is a schematic diagram of the change of superimposed amplitude values ​​versus time domain contour lines in the Bohai A oilfield according to a method for rapid detection of reservoir gas cap risk provided by an embodiment of the present invention;

[0041] Figure 6 This is a schematic diagram of the gas cap boundary predicted for the Bohai A oilfield reservoir according to a method for rapid detection of reservoir gas cap risk provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0042] Exemplary embodiments of the present invention will be described in more detail below with reference to the accompanying drawings. Although exemplary embodiments of the present invention are shown in the accompanying drawings, it should be understood that the present invention can be implemented in various forms and should not be limited by the embodiments described herein. Rather, these embodiments are provided to enable a more thorough understanding of the present invention and to fully convey the scope of the present invention to those skilled in the art.

[0043] Because the current fluid detection methods are complex, time-consuming, and have limited applicability in predicting reservoir gas cap risks. The present invention provides a method for rapid detection of reservoir gas cap risks, including clarifying the seismic response characteristics of gas layers through rock physical analysis; determining whether the amplitude attribute anomaly of reservoir seismic data is consistent with the seismic response characteristics of gas layers; extracting the superimposed amplitude values ​​of contour lines in different time domains; delineating the gas-bearing boundaries of contour lines in the depth domain to determine the gas-bearing area, and obtaining the gas cap value by calculation. The present invention obtains the superimposed amplitude values ​​of contour lines in different time domains through a three-step method, so as to facilitate conversion to the depth domain to delineate the gas-bearing boundaries and calculate the gas cap value, thereby realizing rapid prediction of reservoir gas cap risk according to the size of the gas cap value.

[0044] The scheme of the present invention is described in detail below through examples.

[0045] Example

[0046] like Figure 1 As shown, the present invention provides a method for rapid detection of reservoir gas cap risk, comprising:

[0047] Based on post-stack seismic data and rock physics analysis, the seismic response characteristics of the gas layer are clarified through rock physics analysis;

[0048] Based on the existence of amplitude response anomalies, the top and bottom surfaces of the reservoir are interpreted, and the amplitude of the reservoir seismic data is extracted to determine whether the amplitude attribute anomalies of the reservoir seismic data are consistent with the seismic response characteristics of the gas layer;

[0049] Seismic flat point amplitude enhancement is used to stack the amplitude of reservoir seismic data along the interpreted time domain contour lines of the reservoir layer, further amplifying the amplitude difference between gas layers and non-gas layers, and extracting the superimposed amplitude values ​​of different time domain contour lines;

[0050] The extracted superimposed amplitude values ​​are converted into depth values ​​corresponding to the contour lines in different time domains. The gas-bearing area is determined by delineating the gas-bearing boundaries of the depth domain contour lines, and the gas cap value is obtained by calculation.

[0051] Furthermore, the rapid detection method for reservoir gas cap risk also includes a method of using seismic flat point amplitude enhancement to superimpose the amplitude of reservoir seismic data along the interpreted time domain contour lines of the reservoir plane:

[0052] Based on post-stack pure wave seismic data, the top surface of the reservoir time domain is traced and interpreted, and the top surface time grid of the reservoir time domain is constructed;

[0053] Write a self-developed program code, input the post-stack pure wave seismic data and the constructed reservoir time domain top surface time grid, use the time domain contour line of the reservoir time domain top surface time grid as the superposition unit, and superimpose the amplitude values ​​of the pure wave seismic data sample points corresponding to the same time domain contour line;

[0054] The transformation relationship between the superimposed amplitude value and the time domain contour line is constructed, and the superimposed amplitude values ​​corresponding to different time domain contour lines are counted.

[0055] Furthermore, the rapid detection method for reservoir gas cap risk also includes a method for constructing a time grid of the reservoir top surface in the time domain: tracking seismic peaks, troughs or zero-value points in the process of interpreting the top surface layer of the reservoir time domain, and generating a reservoir top surface time grid with preset 10*10 grid parameters.

[0056] Furthermore, the rapid detection method for reservoir gas cap risk also includes a method for superimposing the amplitude values ​​of pure wave seismic data sample points corresponding to the same time domain contour line: the amplitude values ​​of pure wave seismic data sample points corresponding to the same time domain contour line are superimposed using the amplitude superposition formula, and the amplitude superposition formula is:

[0057] A(t0)=∑ x,y D(x,y,t=t0) / N x,y ;

[0058] Where A(t0) is the superimposed amplitude value of the time domain contour line value t0, D(x,y,t=t0) is the amplitude value of the three-dimensional pure wave seismic data at the sample point position, and t=t0, x is the main survey line direction, y is the connecting survey line direction, t is the time direction, and N x,y It plays a normalizing role for the total number of sample points involved in the accumulation.

[0059] Furthermore, the rapid detection method for reservoir gas cap risk also includes constructing a transformation relationship between the superimposed amplitude value and the time domain contour line, and a method for counting the superimposed amplitude values ​​corresponding to different time domain contour lines:

[0060] Based on the obtained superimposed amplitude values, a structural map showing the transformation of the superimposed amplitude values ​​with the time domain contour lines at the reservoir level is constructed, with the time domain contour lines as the horizontal axis and the superimposed amplitude values ​​as the vertical axis. The superimposed amplitude values ​​corresponding to different time domain contour lines are statistically analyzed based on the order of the time domain contour lines from shallow to deep, and the time domain contour line corresponding to the amplitude mutation position is predicted to be the gas-bearing boundary.

[0061] Among them, according to the time-depth conversion formula of the study area, the extracted superimposed amplitude values ​​are converted through time-depth conversion to obtain the depth domain contour line values ​​corresponding to different time domain contour lines, especially the depth domain values ​​corresponding to the predicted time domain gas-bearing boundary values.

[0062] Furthermore, the rapid detection method for reservoir gas cap risk also includes a method for delineating the gas-bearing boundary in the reservoir top depth domain:

[0063] Based on the reservoir top surface structural map, the extracted stacking amplitude values ​​are assigned to all sample points corresponding to the corresponding depth domain contour lines, so that the same depth domain contour lines correspond to the same stacking amplitude values, while the stacking amplitude values ​​corresponding to different depth domain contour lines are different.

[0064] Interpolate the stacked amplitude values ​​of the contour lines in each depth domain to obtain the stacked amplitude value plane attributes corresponding to the reservoir top surface structure map;

[0065] Based on the boundary between gas and non-gas layers in the plane properties of the superimposed amplitude value, the gas-bearing boundary of the depth domain contour is delineated. The area above the boundary reflects the superimposed amplitude value characteristics of the gas layer, while the area below the boundary reflects the superimposed amplitude value characteristics of the non-gas layer.

[0066] Specifically, the method of the present invention is used to implement the implementation method of reservoir gas cap prediction: taking the Shahejie Formation reservoir of Bohai A Oilfield as an example, the specific steps are:

[0067] S1, based on post-stack seismic data and rock physics analysis. Figure 2 As shown in the figure, rock physical analysis confirmed that the gas layer in oilfield A has low impedance characteristics, and the overlying surrounding rock also has low impedance characteristics. The reflection interface formed by the gas layer and the overlying surrounding rock corresponds to the seismic weak reflection response characteristics.

[0068] S2, based on the existence of amplitude response anomalies, interpret the reservoir top and bottom surfaces and extract the amplitude of reservoir seismic data. Figure 3 As shown, the potential gas layer corresponds to the weak seismic amplitude anomaly in the amplitude attribute, which is consistent with the weak seismic reflection response characteristics of the gas layer obtained from rock physics analysis. Step S3 can be carried out;

[0069] S3, using earthquake flat point amplitude enhancement as a means. Figure 4 As shown in the figure, the amplitude of the seismic data is superimposed along the time domain contour line of the interpreted top surface of the A oil field reservoir, and the amplitude superposition attribute of the reservoir seismic data along the time domain contour line is obtained, which further amplifies the amplitude difference between the gas layer and the non-gas layer. Further statistics are made on the amplitude superposition of the A oil field reservoir seismic data along the time domain contour line, and the superposition amplitude values ​​of different time domain contour lines are as follows: Figure 5As shown, the time domain contour line corresponding to the amplitude mutation position is 2025ms, which is predicted to be the gas top boundary. The area shallower than the time domain contour line is considered to be the gas layer.

[0070] S4, such as Figure 6 As shown in the figure, through time-depth conversion, the time domain contour line of 2025ms in oil field A corresponds to a depth of 2240m. Based on this, the gas-bearing boundary is delineated in the depth domain structural map of the reservoir top surface. The gas-bearing area is 0.6km2, and the calculated gas cap size is 258,000 cubic meters.

[0071] The rapid detection method for reservoir gas cap risk of the present invention adopts a three-step method of clarifying the seismic response characteristics of the gas layer - determining whether the amplitude attribute anomaly is consistent with the seismic response characteristics of the gas layer - extracting the superimposed amplitude values ​​of the contour lines in different time domains. The gas-bearing boundary of the depth domain contour line can be delineated to determine the gas-bearing area, so that the gas cap value can be obtained through corresponding calculations, thereby realizing a rapid prediction of the reservoir gas cap risk according to the size of the gas cap value. The process is simple and can effectively play the guiding role of geophysical fluid detection in oil field development methods and well network and well type design.

[0072] Based on the above-mentioned method for rapid detection of reservoir gas cap risk, the present invention further provides an analysis device, comprising:

[0073] The first processing unit is used to determine the seismic response characteristics of the gas layer through rock physical analysis based on the post-stack seismic data;

[0074] The second processing unit is used to interpret the reservoir layer and extract the amplitude of the reservoir seismic data, and determine whether the amplitude attribute anomaly of the reservoir seismic data is consistent with the seismic response characteristics of the gas layer;

[0075] The third processing unit is used to superimpose the amplitude of the reservoir seismic data along the interpreted time domain contour line of the reservoir layer using seismic flat point amplitude enhancement, and extract the superimposed amplitude values ​​of different time domain contour lines;

[0076] The fourth processing unit is used to convert the extracted superimposed amplitude values ​​into depth domain contour values ​​corresponding to different time domain contour lines through time-depth conversion, determine the gas-bearing area by delineating the gas-bearing boundaries of the depth domain contour lines, and obtain the gas cap value by calculation.

[0077] Based on the above-mentioned method for rapid detection of reservoir gas cap risk, the present invention also provides a computer-readable storage medium having a computer program stored thereon. When the computer program is executed by a processor, the steps of the above-mentioned method for rapid detection of reservoir gas cap risk are implemented.

[0078] Based on the above-mentioned method for rapid detection of reservoir gas cap risks, the present invention also provides a computer device, including a memory, a processor, and a computer program stored in the memory and runnable on the processor. When the processor executes the computer program, the steps of the above-mentioned method for rapid detection of reservoir gas cap risks are implemented.

[0079] The present invention is described in terms of flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to specific embodiments. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as a combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 a process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0080] These computer program instructions may also be stored in a computer readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 a process or multiple processes and / or boxes Figure 1 The function specified in one or more boxes.

[0081] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operational steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing the instructions executed on the computer or other programmable device for implementing the process. Figure 1 a process or multiple processes and / or boxes Figure 1 A step that specifies a function in one or more boxes.

[0082] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the various embodiments of the present invention.

Claims

1. A method for rapid detection of reservoir gas cap risk, characterized in that: The reservoir gas cap risk rapid detection method comprises: Based on post-stack seismic data, the seismic response characteristics of the gas layer are clarified through rock physical analysis; Interpret reservoir layers and extract the amplitude of reservoir seismic data to determine if the amplitude attribute anomalies of reservoir seismic data are consistent with the seismic response characteristics of the gas layer; Seismic flat point amplitude enhancement is used to stack the amplitude of reservoir seismic data along the interpreted time domain contour lines of the reservoir layer, and the stacked amplitude values ​​of different time domain contour lines are extracted; The extracted superimposed amplitude values ​​are converted into depth domain contour values ​​corresponding to different time domain contour lines through time-depth conversion. The gas-bearing area is determined by delineating the gas-bearing boundaries of the depth domain contour lines, and the gas cap value is obtained by calculation. Also included is a method for stacking the amplitude of reservoir seismic data along the interpreted time domain contour lines of the reservoir plane using seismic flat point amplitude enhancement: Based on post-stack pure wave seismic data, the top surface of the reservoir time domain is traced and interpreted, and the top surface time grid of the reservoir time domain is constructed; The time domain contour lines of the reservoir top surface time grid are used as superposition units, and the amplitude values ​​of the pure wave seismic data samples corresponding to the same time domain contour lines are superimposed; Construct the transformation relationship between the superimposed amplitude value and the time domain contour line, and count the superimposed amplitude values ​​corresponding to different time domain contour lines; It also includes methods for delineating gas-bearing boundaries: Based on the reservoir top surface structural map, the extracted stacking amplitude values ​​are assigned to all sample points corresponding to the corresponding depth domain contour lines; Interpolate the stacked amplitude values ​​of the contour lines in each depth domain to obtain the stacked amplitude value plane attributes corresponding to the reservoir top surface structure map; According to the depth domain contour values ​​corresponding to the time domain contour values ​​of the predicted gas-bearing boundary, the boundary between the gas layer and the non-gas layer in the plane attributes of the superimposed amplitude value is determined, and the gas-bearing boundary of the depth domain contour line is delineated.

2. A method for rapid detection of reservoir gas cap risk according to claim 1, characterized in that: It also includes a method for constructing a time grid on the top surface of a reservoir time domain: tracking seismic peaks, troughs or zero-value points during the interpretation of the top surface layer of the reservoir time domain, and generating a time grid on the top surface of the reservoir with preset grid parameters.

3. A rapid detection method for reservoir gas cap risk according to claim 1, characterized in that: The method also includes a method for superimposing the amplitude values ​​of pure wave seismic data samples corresponding to the same time domain contour line: the amplitude values ​​of pure wave seismic data samples corresponding to the same time domain contour line are superimposed using an amplitude superposition formula, and the amplitude superposition formula is: ； in, The time domain contour value is The superimposed amplitude value, is the amplitude value of the three-dimensional pure wave seismic data at the sample point position, The main survey line direction, For the direction of the communication line, is the time direction, is the total number of sample points involved in the accumulation.

4. A method for rapid detection of reservoir gas cap risk according to claim 1, characterized in that: It also includes the method of constructing the transformation relationship between the superimposed amplitude value and the time domain contour line, and counting the superimposed amplitude values ​​corresponding to different time domain contour lines: Based on the obtained superimposed amplitude values, a structural map showing the transformation of the superimposed amplitude values ​​with the time domain contour lines at the reservoir level is constructed, with the time domain contour lines as the horizontal axis and the superimposed amplitude values ​​as the vertical axis. The superimposed amplitude values ​​corresponding to different time domain contour lines are statistically analyzed based on the order of the time domain contour lines from shallow to deep, and the time domain contour line corresponding to the amplitude mutation position is predicted to be the gas-bearing boundary.

5. A rapid detection method for reservoir gas cap risk according to claim 4, characterized in that: The area above the gas-bearing boundary is the superimposed amplitude value characteristic of the reaction gas layer, and the area below the gas-bearing boundary is the superimposed amplitude value characteristic reflecting the non-gas layer.

6. An analysis device using the rapid detection method for reservoir gas cap risk according to any one of claims 1 to 5, characterized in that: The analysis device comprises: The first processing unit is used to determine the seismic response characteristics of the gas layer through rock physical analysis based on the post-stack seismic data; The second processing unit is used to interpret the reservoir layer and extract the amplitude of the reservoir seismic data, and determine whether the amplitude attribute anomaly of the reservoir seismic data is consistent with the seismic response characteristics of the gas layer; The third processing unit is used to superimpose the amplitude of the reservoir seismic data along the interpreted time domain contour line of the reservoir layer using seismic flat point amplitude enhancement, and extract the superimposed amplitude values ​​of different time domain contour lines; The fourth processing unit is used to convert the extracted superimposed amplitude values ​​into depth domain contour values ​​corresponding to different time domain contour lines through time-depth conversion, determine the gas-bearing area by delineating the gas-bearing boundaries of the depth domain contour lines, and obtain the gas cap value by calculation.

7. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method for rapid detection of reservoir gas cap risk according to any one of claims 1 to 5 are implemented.

8. A computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the computer program, the steps of the method for rapid detection of reservoir gas cap risk according to any one of claims 1 to 5 are implemented.

Images