CN112965108B - Method and system for determining vertical closure of trap cap layer - Google Patents
Method and system for determining vertical closure of trap cap layer Download PDFInfo
- Publication number
- CN112965108B CN112965108B CN202110160884.7A CN202110160884A CN112965108B CN 112965108 B CN112965108 B CN 112965108B CN 202110160884 A CN202110160884 A CN 202110160884A CN 112965108 B CN112965108 B CN 112965108B
- Authority
- CN
- China
- Prior art keywords
- mudstone
- cover layer
- mudstone cover
- target interval
- post
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Abstract
The invention provides a method and a system for determining vertical closure of a trap cover layer. The method for determining the vertical closure of the trap cover layer comprises the following steps: acquiring post-stack seismic data volume, mudstone cover layer parameters and layers of a target interval; dividing the target interval into multiple types of mudstone cover layers according to the parameters of the mudstone cover layers; creating a low-frequency interpolation model according to reservoir inversion sensitive parameter curves corresponding to the horizons and various mudstone cover layers; generating mudstone cover layer thickness data of the target interval according to the low-frequency interpolation model, the post-stack seismic data body and cut-off values corresponding to various mudstone cover layers; and determining the vertical closure of the trap cover layer according to the mudstone cover layer thickness data. The method can accurately determine the vertical closure of the trap cap layer, so as to guide exploration and deployment and improve the success rate of the exploratory well.
Description
Technical Field
The invention relates to the technical field of reservoir geology, in particular to a method and a system for determining vertical closure of a trap cap layer.
Background
The cap layer is a protective layer which is positioned above the reservoir layer and can seal the reservoir layer to prevent oil gas in the reservoir layer from escaping upwards, and the vertical sealing property of the cap layer is used for determining the vertical sealing capacity of the trapped cap layer, so that whether the cap layer can seal the oil gas or not is directly related to preventing the oil gas from escaping upwards. Accurate prediction of trap cap vertical closure has an important impact on the reservoir analysis and overall evaluation of target traps. For the trap to be drilled, the vertical closure of the target trap cover layer is accurately predicted before drilling, so that the drilling success rate of the exploratory well is greatly improved.
The vertical closure of the trap cover layer is developed mainly aiming at the deposition environment, the thickness of a argillaceous single layer and the closure of faults, the main parameters affecting the vertical closure of the trap cover layer include porosity, permeability, argillaceous content, cover layer thickness, displacement pressure and the like, wherein the cover layer argillite cover layer thickness is one of the key parameters for evaluating the vertical closure of the trap cover layer, and exploration practices prove that the argillite cover layer is often not a pure argillite cover layer, the cover layer is often thick-layer argillite and thin-layer sandstone or sand-argillite interbed deposition, only lithology inversion is carried out on the trap cover combination, the lithology characteristics of the trap cover layer are finely carved, the argillite cover layer on the reservoir layer is comprehensively evaluated, and the influence of sandstone in the argillite cover layer is removed, so that the overall closure of the cover layer can be accurately evaluated. The conventional technical scheme at present cannot accurately determine the vertical closure of the trap cap layer, so that the success rate of drilling the structural trap during exploration and deployment is low.
Disclosure of Invention
The embodiment of the invention mainly aims to provide a method and a system for determining the vertical closure of a trap cover layer, so as to accurately determine the vertical closure of the trap cover layer, guide exploration deployment and improve the success rate of a exploratory well.
In order to achieve the above objective, an embodiment of the present invention provides a method for determining vertical closure of a trap cap layer, including:
acquiring post-stack seismic data volume, mudstone cover layer parameters and layers of a target interval;
dividing the target interval into multiple types of mudstone cover layers according to the parameters of the mudstone cover layers;
creating a low-frequency interpolation model according to reservoir inversion sensitive parameter curves corresponding to the horizons and various mudstone cover layers;
generating mudstone cover layer thickness data of the target interval according to the low-frequency interpolation model, the post-stack seismic data body and cut-off values corresponding to various mudstone cover layers;
and determining the vertical closure of the trap cover layer according to the mudstone cover layer thickness data.
The embodiment of the invention also provides a system for determining the vertical closure of the trap cover layer, which comprises the following steps:
the acquisition unit is used for acquiring post-stack seismic data volumes, mudstone cover layer parameters and horizons of the target interval;
the dividing unit is used for dividing the target interval into multiple types of mudstone cover layers according to the mudstone cover layer parameters;
the low-frequency interpolation model creation unit is used for creating a low-frequency interpolation model according to reservoir inversion sensitive parameter curves corresponding to the horizons and various mudstone cover layers;
the mudstone cover layer thickness unit is used for generating mudstone cover layer thickness data of the target layer section according to the low-frequency interpolation model, the post-stack seismic data body and cut-off values corresponding to various mudstone cover layers;
and the vertical sealing unit is used for determining the vertical sealing of the trap cover layer according to the mudstone cover layer thickness data.
The embodiment of the invention also provides computer equipment, which comprises a memory, a processor and a computer program stored on the memory and running on the processor, wherein the processor realizes the steps of the method for determining the vertical closure of the trap cover layer when executing the computer program.
The embodiment of the invention also provides a computer readable storage medium, wherein a computer program is stored on the computer readable storage medium, and the computer program realizes the steps of the method for determining the vertical closure of the trap cover layer when being executed by a processor.
According to the method and the system for determining the vertical closure of the trap cap layer, a target layer section is divided into a plurality of types, then a low-frequency interpolation model is created according to a horizon and a reservoir inversion sensitive parameter curve corresponding to each type, then mudstone cap layer thickness data of the target layer section are generated according to the low-frequency interpolation model, a post-stack seismic data body and cut-off values corresponding to each type, and finally the vertical closure of the trap cap layer is accurately determined according to the mudstone cap layer thickness data, so that exploration deployment is guided, and the success rate of exploratory well detection is improved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are needed in the description of the embodiments will be briefly described below, it being obvious that the drawings in the following description are only some embodiments of the present invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a method of determining vertical closure of a containment cap layer in accordance with an embodiment of the present invention;
FIG. 2 is a flow chart of a method of determining vertical closure of a containment cap layer in accordance with another embodiment of the present invention;
FIG. 3 is a flowchart of S102 in an embodiment of the invention;
FIG. 4 is a flowchart of S103 in an embodiment of the present invention;
FIG. 5 is a flowchart of S104 in an embodiment of the invention;
FIG. 6 is a preferred plot of reservoir inversion sensitivity parameters in an embodiment of the invention;
FIG. 7 is a schematic illustration of a mudstone cover standard plate in an embodiment of the invention;
FIG. 8 is a schematic diagram of a reservoir inversion sensitivity parameter for a class I mudstone cap in an embodiment of the invention;
FIG. 9 is a schematic diagram of a reservoir inversion sensitivity parameter for a class II mudstone cap layer in accordance with an embodiment of the present invention;
FIG. 10 is a schematic diagram of a reservoir inversion sensitivity parameter for a class III mudstone cap in an embodiment of the invention;
FIG. 11 is a schematic diagram of a low frequency interpolation model of a class I mudstone cover layer in accordance with an embodiment of the present invention;
FIG. 12 is a cross-sectional view of a post-stack inverted lithology data volume of a class I mudstone cap in an embodiment of the present invention;
FIG. 13 is a schematic view of a mudstone cap layer thickness profile of a class I mudstone cap layer in accordance with an embodiment of the present invention;
FIG. 14 is a schematic view of the vertical effective mudstone cap thickness of a target interval in an embodiment of the invention;
FIG. 15 is a block diagram of a vertical closure determination system for a containment cap in accordance with an embodiment of the present invention;
fig. 16 is a block diagram showing the structure of a computer device in the embodiment of the present invention.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Those skilled in the art will appreciate that embodiments of the invention may be implemented as a system, apparatus, device, method, or computer program product. Accordingly, the present disclosure may be embodied in the following forms, namely: complete hardware, complete software (including firmware, resident software, micro-code, etc.), or a combination of hardware and software.
In view of the fact that the vertical closure of the trap cap layer cannot be accurately determined by the conventional technical scheme at present, the success rate of drilling the structural trap during exploration and deployment is low, the embodiment of the invention provides a method for determining the vertical closure of the trap cap layer, so that the vertical closure of the trap cap layer is accurately determined, quantitative evaluation is carried out on the effectiveness of the trap cap, then drilling is carried out by optimizing the trap with good trap closure during exploration and deployment, and the success rate of a exploratory well is improved. The present invention will be described in detail with reference to the accompanying drawings.
FIG. 1 is a flow chart of a method for determining the vertical closure of a containment cap in an embodiment of the invention. FIG. 2 is a flow chart of a method for determining the vertical closure of a containment cap in accordance with another embodiment of the present invention. As shown in fig. 1-2, the method for determining the vertical closure of the trap cover layer comprises the following steps:
s101: and acquiring post-stack seismic data volume, mudstone cover layer parameters and layers of the target interval.
The well data in fig. 2 comprise layering data and logging data, and the layering data is calibrated and finely construed to obtain the horizons of the target intervals.
The logging data includes acoustic curves, gamma curves, and density curves. FIG. 6 is a preferred plot of reservoir inversion sensitivity parameters in an embodiment of the invention. As shown in fig. 6, the abscissa in fig. 6 is gamma in units of API; the ordinate is percent, in%. The well logging data of the target interval of the research area are subjected to histogram analysis and intersection chart analysis, and as can be seen from fig. 6, the target interval with the gamma value larger than 63 is mudstone and the target interval with the gamma value smaller than or equal to 63 is sandstone, so that the mudstone and the sandstone can be distinguished through a gamma curve, and the gamma curve is a reservoir inversion sensitive parameter curve.
S102: the target intervals are divided into multiple types of mudstone cover layers according to the mudstone cover layer parameters.
Fig. 3 is a flowchart of S102 in an embodiment of the present invention. As shown in fig. 3, S102 includes:
s201: and establishing a mudstone cover standard plate according to the parameters of the mudstone cover.
In one embodiment, a standard shale overburden pattern for each target interval may be created based on formation characteristics for different intervals.
Fig. 7 is a schematic illustration of a mudstone cover standard plate in an embodiment of the invention. As shown in FIG. 7, the abscissa of FIG. 7 is the specific surface area in m 2 ·g -1 The method comprises the steps of carrying out a first treatment on the surface of the The ordinate is the breakthrough pressure in Mpa. In specific implementation, parameters (including breakthrough pressure, specific surface area and permeability) of a mudstone cover layer of a target interval sampling point can be projected into a coordinate system, then the sampling points are divided into three types according to the permeability (i is the sampling point with the permeability of less than 300md, ii is the sampling point with the permeability of more than or equal to 300md and less than or equal to 600md, and iii is the sampling point with the permeability of more than 600 md), and finally a dividing curve for dividing the three types of sampling points into three areas is generated to establish a mudstone cover layer standard pattern. As shown in fig. 7, the porosities of the sampling points located on the two dividing curves are 10% and 20%, respectively.
S202: the target intervals are divided into multiple types of mudstone cover layers according to a mudstone cover layer standard plate.
In practice, mudstone cover layers can be classified into class I, class II and class III according to standard patterns of the mudstone cover layers. The lower limit of the thickness of the effective mudstone cover layer of each mudstone cover layer can be determined by counting the thicknesses of the mudstone cover layers in the various mudstone cover layers, and the lower limit of the thickness of the effective mudstone cover layer of each mudstone cover layer is used as a cut-off value corresponding to each mudstone cover layer. For example, the lower limit (cutoff value) of the effective mudstone cap thickness of the class i mudstone cap is 5m, the lower limit (cutoff value) of the effective mudstone cap thickness of the class ii mudstone cap is 10m, and the lower limit (cutoff value) of the effective mudstone cap thickness of the class iii mudstone cap is 20m.
S103: and creating a low-frequency interpolation model according to the reservoir inversion sensitive parameter curves corresponding to the horizons and various mudstone cover layers.
Fig. 4 is a flowchart of S103 in the embodiment of the present invention. As shown in fig. 4, S103 includes:
s301: and coarsening a corresponding reservoir inversion sensitive parameter curve according to the cutoff value.
FIG. 8 is a schematic representation of a reservoir inversion sensitivity parameter for a class I mudstone cap in an embodiment of the invention. FIG. 9 is a schematic diagram of a reservoir inversion sensitivity parameter for a class II mudstone cap in an embodiment of the invention. FIG. 10 is a schematic representation of a reservoir inversion sensitivity parameter for a class III mudstone cap in an embodiment of the invention. As shown in fig. 8-10, the target interval is divided into J-i, J-ii and J-iii according to the formation depth, the reservoir inversion sensitivity parameter curve is a gamma curve GR, the abscissa is Gamma (GR), the unit is API, the ordinate is depth, and the unit is m.
The higher the cut-off value corresponding to the mudstone cover layer type is, the lower the sampling rate of the coarsened gamma curve is. As shown in fig. 8-10, the class i mudstone cover layer has a cutoff value of 5m, and the coarsened gamma curve in fig. 8 has a sampling rate of 2m; the cutoff value of the II-type mudstone cover layer is 10m, and the sampling rate of the gamma curve after coarsening in FIG. 8 is 1m; the class iii mudstone cover has a cutoff of 20m and the gamma curve after coarsening in fig. 8 has a sampling rate of 0.5m.
S302: and creating a low-frequency interpolation model according to the coarsened reservoir inversion sensitive parameter curves corresponding to the horizons and various mudstone cover layers.
FIG. 11 is a schematic diagram of a low frequency interpolation model of a class I mudstone cover layer in an embodiment of the invention. As shown in fig. 11, the abscissa of fig. 11 is distance in km; the ordinate is depth in m. RL-1, RC-2, RC-1, RN-2, RN-3A and RN-1 are all well logs for drilling a class I mudstone cover layer.
In the specific implementation, taking the class I mudstone cover layer as an example, a horizon can be loaded into inversion software, a simple structure model of a target interval is created, the stratum contact relationship is equal division, and a low-frequency interpolation model is created by adopting a local weighted interpolation method under the constraint of a coarsened reservoir inversion sensitive parameter curve corresponding to the class I mudstone cover layer.
S104: and generating mudstone cover layer thickness data of the target interval according to the low-frequency interpolation model, the post-stack seismic data body and cut-off values corresponding to various mudstone cover layers.
Fig. 5 is a flowchart of S104 in the embodiment of the present invention. As shown in fig. 5, S104 includes:
s401: and generating a post-stack inversion lithology data body according to the low-frequency interpolation model and the post-stack seismic data body.
FIG. 12 is a cross-sectional view of a post-stack inverted lithology data volume of a class I mudstone cap in an embodiment of the invention. As shown in fig. 12, the abscissa of fig. 12 is distance in km; the ordinate is depth in m.
For example, based on analysis of the well data of the research area, the post-stack inversion lithology data body and the low-frequency interpolation model, the lithology inversion software system is used for setting the effective sample number, the high-frequency range and the target sampling rate, and performing argillaceous content inversion on the I-type mudstone cover layer, the II-type mudstone cover layer and the III-type mudstone cover layer to generate the post-stack inversion lithology data body.
S402: and generating mudstone cover thickness data of the target interval according to the cut-off values corresponding to various mudstone covers and the post-stack inversion lithology data body.
FIG. 13 is a schematic view of a mudstone cover layer thickness profile of a class I mudstone cover layer in accordance with an embodiment of the present invention. As shown in fig. 13, the abscissa of fig. 13 is distance in km; the ordinate is depth in m.
In specific implementation, cut-off values can be set on inversion software, and coarsened mudstone cover layer thickness data in the vertical direction of the post-stack inversion lithology data body of the I-type mudstone cover layer, the II-type mudstone cover layer and the III-type mudstone cover layer can be respectively obtained.
S105: and determining the vertical closure of the trap cover layer according to the mudstone cover layer thickness data.
FIG. 14 is a schematic view of the vertical effective mudstone cap thickness of a target interval in an embodiment of the invention. As shown in fig. 14, the abscissa of fig. 14 is distance in km; the ordinate is distance in km. In FIG. 14, there are multiple well logs (e.g., S-1, RNW-1, RNN-1, RL-1, RC-2, RC-1, RN-2, RN-3A, RNNE-1, RN-1, RNE-1, TE-1, R-1A, etc.) of various mudstone cover layers encountered by the drill, where the beginning sequence number of L in FIG. 14 represents the main line number and the beginning sequence number of T represents the track number. Mudstone cover thickness data (vertical effective mudstone cover thickness) for the target interval can be obtained from fig. 14.
In the specific implementation, thickness data of the class I mudstone cover layer, the class II mudstone cover layer and the class III mudstone cover layer in the target interval of the research area can be counted respectively to obtain a schematic diagram of the vertical effective mudstone cover layer thickness of the target interval of the research area. According to the schematic diagram of the vertical effective mudstone cover layer thickness in the target interval or the upper coating section of the target interval, the vertical closure of the mudstone cover layer to be drilled and closed can be determined, and under the condition that other conditions are the same, the trap with good vertical closure of the cover layer is preferably selected as a drilling deployment target, so that exploration deployment is guided, and the success rate of a exploratory well is improved.
The method for determining the vertical closure of the trap cover layer shown in fig. 1 may be performed by a computer. As can be seen from the flow shown in fig. 1, the method for determining the vertical closure of the trap cap layer in the embodiment of the invention divides the target layer segment into a plurality of types, creates a low-frequency interpolation model according to the reservoir inversion sensitive parameter curves corresponding to the layers and the types, generates mudstone cap layer thickness data of the target layer segment according to the low-frequency interpolation model, the post-stack seismic data body and the cut-off values corresponding to the types, and finally accurately determines the vertical closure of the trap cap layer according to the mudstone cap layer thickness data, thereby guiding exploration deployment and improving the success rate of exploratory wells.
The specific flow of the embodiment of the invention is as follows:
1. and acquiring post-stack seismic data volumes, mudstone cover layer parameters and layering data of the target interval.
2. And calibrating and finely constructing and explaining the layering data to obtain the horizon of the target interval.
3. And establishing a mudstone cover standard plate according to the parameters of the mudstone cover.
4. The target intervals are divided into multiple types of mudstone cover layers according to a mudstone cover layer standard plate.
5. And coarsening a corresponding reservoir inversion sensitive parameter curve according to the cutoff value.
6. And creating a low-frequency interpolation model according to the coarsened reservoir inversion sensitive parameter curves corresponding to the horizons and various mudstone cover layers.
7. And generating a post-stack inversion lithology data body according to the low-frequency interpolation model and the post-stack seismic data body.
8. And generating mudstone cover thickness data of the target interval according to the cut-off values corresponding to various mudstone covers and the post-stack inversion lithology data body.
9. And determining the vertical closure of the trap cover layer according to the mudstone cover layer thickness data.
In summary, the method for determining the vertical closure of the trap cap layer in the embodiment of the invention coarsens inversion sensitive parameters of a reservoir according to different scales by counting the obtained standard pattern plates of the target layer section mudstone cap layer, further obtains coarsened lithologic inversion data bodies, can accurately determine the vertical closure of the trap mudstone cap layer, and is favorable for drilling the trap cap layer under the condition of equal other conditions during exploration and deployment, thereby improving the drilling success rate.
Based on the same inventive concept, the embodiment of the invention also provides a system for determining the vertical closure of the trap cover layer, and because the principle of solving the problem of the system is similar to that of the method for determining the vertical closure of the trap cover layer, the implementation of the system can refer to the implementation of the method, and the repetition is omitted.
FIG. 15 is a block diagram of a vertical closure determination system for a containment cap in accordance with an embodiment of the present invention. As shown in fig. 15, the trap cap layer vertical closure determination system includes:
the acquisition unit is used for acquiring post-stack seismic data volumes, mudstone cover layer parameters and horizons of the target interval;
the dividing unit is used for dividing the target interval into multiple types of mudstone cover layers according to the mudstone cover layer parameters;
the low-frequency interpolation model creation unit is used for creating a low-frequency interpolation model according to reservoir inversion sensitive parameter curves corresponding to the horizons and various mudstone cover layers;
the mudstone cover layer thickness unit is used for generating mudstone cover layer thickness data of the target layer section according to the low-frequency interpolation model, the post-stack seismic data body and cut-off values corresponding to various mudstone cover layers;
and the vertical sealing unit is used for determining the vertical sealing of the trap cover layer according to the mudstone cover layer thickness data.
In one of these embodiments, the mudstone cover thickness unit is specifically for:
generating a post-stack inversion lithology data body according to the low-frequency interpolation model and the post-stack seismic data body;
and generating mudstone cover thickness data of the target interval according to the cut-off values corresponding to various mudstone covers and the post-stack inversion lithology data body.
In one embodiment, the dividing unit is specifically configured to:
establishing a mudstone cover standard chart according to the parameters of the mudstone cover;
the target intervals are divided into multiple types of mudstone cover layers according to a mudstone cover layer standard plate.
In one embodiment, the low frequency interpolation model creation unit is specifically configured to:
coarsening a corresponding reservoir inversion sensitive parameter curve according to the cutoff value;
and creating a low-frequency interpolation model according to the coarsened reservoir inversion sensitive parameter curves corresponding to the horizons and various mudstone cover layers.
In summary, the vertical closure determination system of the trap cap layer of the embodiment of the invention divides a target layer section into a plurality of types, then creates a low-frequency interpolation model according to the inversion sensitive parameter curves of the layers and the reservoirs corresponding to the types, then generates mudstone cap layer thickness data of the target layer section according to the low-frequency interpolation model, the post-stack seismic data body and the cut-off values corresponding to the types, and finally accurately determines the vertical closure of the trap cap layer according to the mudstone cap layer thickness data, thereby guiding exploration deployment and improving the success rate of exploratory well detection.
The embodiment of the invention also provides a concrete implementation mode of the computer equipment, which can realize all the steps in the method for determining the vertical closure of the trap cover layer in the embodiment. Fig. 16 is a block diagram of a computer device in an embodiment of the present invention, and referring to fig. 16, the computer device specifically includes:
a processor 1601 and a memory 1602.
The processor 1601 is configured to invoke a computer program in the memory 1602, where the processor executes the computer program to implement all the steps in the method for determining vertical closure of a trap cap in the foregoing embodiment, for example, the processor executes the computer program to implement the following steps:
acquiring post-stack seismic data volume, mudstone cover layer parameters and layers of a target interval;
dividing the target interval into multiple types of mudstone cover layers according to the parameters of the mudstone cover layers;
creating a low-frequency interpolation model according to reservoir inversion sensitive parameter curves corresponding to the horizons and various mudstone cover layers;
generating mudstone cover layer thickness data of the target interval according to the low-frequency interpolation model, the post-stack seismic data body and cut-off values corresponding to various mudstone cover layers;
and determining the vertical closure of the trap cover layer according to the mudstone cover layer thickness data.
In summary, the computer equipment of the embodiment of the invention divides the target interval into a plurality of types, creates a low-frequency interpolation model according to the horizon and the reservoir inversion sensitive parameter curves corresponding to the types, generates mudstone cap layer thickness data of the target interval according to the low-frequency interpolation model, the post-stack seismic data body and the cut-off values corresponding to the types, and finally accurately determines the vertical closure of the trap layer according to the mudstone cap layer thickness data, thereby guiding exploration deployment and improving the success rate of exploratory wells.
The embodiment of the present invention further provides a computer readable storage medium capable of implementing all the steps of the method for determining vertical closure of a trap cap in the above embodiment, where the computer readable storage medium stores a computer program, and the computer program when executed by a processor implements all the steps of the method for determining vertical closure of a trap cap in the above embodiment, for example, the processor implements the following steps when executing the computer program:
acquiring post-stack seismic data volume, mudstone cover layer parameters and layers of a target interval;
dividing the target interval into multiple types of mudstone cover layers according to the parameters of the mudstone cover layers;
creating a low-frequency interpolation model according to reservoir inversion sensitive parameter curves corresponding to the horizons and various mudstone cover layers;
generating mudstone cover layer thickness data of the target interval according to the low-frequency interpolation model, the post-stack seismic data body and cut-off values corresponding to various mudstone cover layers;
and determining the vertical closure of the trap cover layer according to the mudstone cover layer thickness data.
In summary, the computer readable storage medium of the embodiment of the invention divides the target interval into a plurality of types, creates a low-frequency interpolation model according to the horizon and the reservoir inversion sensitive parameter curves corresponding to the types, then generates mudstone cap layer thickness data of the target interval according to the low-frequency interpolation model, the post-stack seismic data body and the cut-off values corresponding to the types, and finally accurately determines the vertical closure of the trap layer according to the mudstone cap layer thickness data, thereby guiding exploration and deployment and improving the success rate of exploratory wells.
The foregoing description of the embodiments has been provided for the purpose of illustrating the general principles of the invention, and is not meant to limit the scope of the invention, but to limit the invention to the particular embodiments, and any modifications, equivalents, improvements, etc. that fall within the spirit and principles of the invention are intended to be included within the scope of the invention.
Those of skill in the art will further appreciate that the various illustrative logical blocks (illustrative logical block), units, and steps described in connection with the embodiments of the invention may be implemented by electronic hardware, computer software, or combinations of both. To clearly illustrate this interchangeability of hardware and software, various illustrative components (illustrative components), elements, and steps have been described above generally in terms of their functionality. Whether such functionality is implemented as hardware or software depends upon the particular application and design requirements of the overall system. Those skilled in the art may implement the described functionality in varying ways for each particular application, but such implementation is not to be understood as beyond the scope of the embodiments of the present invention.
The various illustrative logical blocks, or units, or devices described in the embodiments of the invention may be implemented or performed with a general purpose processor, a digital signal processor, an Application Specific Integrated Circuit (ASIC), a field programmable gate array or other programmable logic device, discrete gate or transistor logic, discrete hardware components, or any combination thereof designed to perform the functions described. A general purpose processor may be a microprocessor, but in the alternative, the general purpose processor may be any conventional processor, controller, microcontroller, or state machine. A processor may also be implemented as a combination of computing devices, e.g., a digital signal processor and a microprocessor, a plurality of microprocessors, one or more microprocessors in conjunction with a digital signal processor core, or any other similar configuration.
The steps of a method or algorithm described in connection with the embodiments disclosed herein may be embodied directly in hardware, in a software module executed by a processor, or in a combination of the two. A software module may be stored in RAM memory, flash memory, ROM memory, EPROM memory, EEPROM memory, registers, hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art. In an example, a storage medium may be coupled to the processor such that the processor can read information from, and write information to, the storage medium. In the alternative, the storage medium may be integral to the processor. The processor and the storage medium may reside in an ASIC, which may reside in a user terminal. In the alternative, the processor and the storage medium may reside as distinct components in a user terminal.
In one or more exemplary designs, the above-described functions of embodiments of the present invention may be implemented in hardware, software, firmware, or any combination of the three. If implemented in software, the functions may be stored on a computer-readable medium or transmitted as one or more instructions or code on the computer-readable medium. Computer readable media includes both computer storage media and communication media that facilitate transfer of computer programs from one place to another. A storage media may be any available media that can be accessed by a general purpose or special purpose computer. For example, such computer-readable media may include, but is not limited to, RAM, ROM, EEPROM, CD-ROM or other optical disk storage, magnetic disk storage or other magnetic storage devices, or any other medium that may be used to carry or store program code in the form of instructions or data structures and other data structures that may be read by a general or special purpose computer, or a general or special purpose processor. Further, any connection is properly termed a computer-readable medium, e.g., if the software is transmitted from a website, server, or other remote source via a coaxial cable, fiber optic cable, twisted pair, digital Subscriber Line (DSL), or wireless such as infrared, radio, and microwave, and is also included in the definition of computer-readable medium. The disks (disks) and disks (disks) include compact disks, laser disks, optical disks, DVDs, floppy disks, and blu-ray discs where disks usually reproduce data magnetically, while disks usually reproduce data optically with lasers. Combinations of the above may also be included within the computer-readable media.
Claims (4)
1. A method for determining the vertical closure of a trap cap layer, comprising:
acquiring post-stack seismic data volume, mudstone cover layer parameters and layers of a target interval;
dividing the target interval into multiple types of mudstone cover layers according to the mudstone cover layer parameters;
creating a low-frequency interpolation model according to the reservoir inversion sensitive parameter curves corresponding to the layers and various mudstone cover layers;
generating mudstone cover layer thickness data of a target layer section according to the low-frequency interpolation model, the post-stack seismic data volume and cutoff values corresponding to various mudstone cover layers;
determining the vertical closure of the trap cap layer according to the mudstone cap layer thickness data;
generating mudstone cap layer thickness data of the target interval according to the low-frequency interpolation model, the post-stack seismic data volume and cutoff values corresponding to various mudstone cap layers comprises:
generating a post-stack inversion lithology data body according to the low-frequency interpolation model and the post-stack seismic data body;
generating mudstone cover thickness data of a target interval according to cut-off values corresponding to various mudstone covers and the post-stack inversion lithology data body;
dividing the target interval into multiple types of mudstone cover layers according to the mudstone cover layer parameters comprises:
establishing a mudstone cover standard chart according to the mudstone cover parameters;
dividing the target interval into multiple types of mudstone cover layers according to the standard mudstone cover layer drawing;
the step of establishing the mudstone cover standard plate according to the mudstone cover parameters comprises the following steps:
projecting mudstone cover layer parameters of sampling points of a target interval into a coordinate system, and dividing the sampling points into three types according to permeability; wherein the mudstone cover parameters include breakthrough pressure, specific surface area and permeability;
the creating of the low-frequency interpolation model according to the reservoir inversion sensitive parameter curves corresponding to the horizons and various mudstone cover layers comprises the following steps:
coarsening a corresponding reservoir inversion sensitive parameter curve according to the cutoff value;
and creating a low-frequency interpolation model according to the coarsened reservoir inversion sensitive parameter curves corresponding to the horizons and various mudstone cover layers.
2. A trap cap layer vertical closure determination system, comprising:
the acquisition unit is used for acquiring post-stack seismic data volumes, mudstone cover layer parameters and horizons of the target interval;
the dividing unit is used for dividing the target interval into multiple types of mudstone cover layers according to the mudstone cover layer parameters;
the low-frequency interpolation model creation unit is used for creating a low-frequency interpolation model according to the reservoir inversion sensitive parameter curves corresponding to the layers and various mudstone cover layers;
the mudstone cover layer thickness unit is used for generating mudstone cover layer thickness data of a target layer section according to the low-frequency interpolation model, the post-stack seismic data body and cut-off values corresponding to various mudstone cover layers;
the vertical sealing unit is used for determining the vertical sealing of the trap cover layer according to the mudstone cover layer thickness data;
the mudstone cover layer thickness unit is specifically used for:
generating a post-stack inversion lithology data body according to the low-frequency interpolation model and the post-stack seismic data body;
generating mudstone cover thickness data of a target interval according to cut-off values corresponding to various mudstone covers and the post-stack inversion lithology data body;
the dividing unit is specifically configured to:
establishing a mudstone cover standard chart according to the mudstone cover parameters;
dividing the target interval into multiple types of mudstone cover layers according to the standard mudstone cover layer drawing;
projecting mudstone cover layer parameters of sampling points of a target interval into a coordinate system, and dividing the sampling points into three types according to permeability; wherein the mudstone cover parameters include breakthrough pressure, specific surface area and permeability;
the low-frequency interpolation model creation unit is specifically configured to:
coarsening a corresponding reservoir inversion sensitive parameter curve according to the cutoff value;
and creating a low-frequency interpolation model according to the coarsened reservoir inversion sensitive parameter curves corresponding to the horizons and various mudstone cover layers.
3. A computer device comprising a memory, a processor and a computer program stored on the memory and running on the processor, characterized in that the processor implements the steps of the trap vertical closure determination method of claim 1 when the computer program is executed.
4. A computer readable storage medium having stored thereon a computer program, characterized in that the computer program when executed by a processor realizes the steps of the method for determining the vertical closure of a trap cap as claimed in claim 1.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110160884.7A CN112965108B (en) | 2021-02-05 | 2021-02-05 | Method and system for determining vertical closure of trap cap layer |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110160884.7A CN112965108B (en) | 2021-02-05 | 2021-02-05 | Method and system for determining vertical closure of trap cap layer |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112965108A CN112965108A (en) | 2021-06-15 |
| CN112965108B true CN112965108B (en) | 2023-08-22 |
Family
ID=76274404
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110160884.7A Active CN112965108B (en) | 2021-02-05 | 2021-02-05 | Method and system for determining vertical closure of trap cap layer |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112965108B (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104063626A (en) * | 2014-07-11 | 2014-09-24 | 中国石油大学(华东) | Predication method of outer-source entrapment fullness degree |
| RU2646528C1 (en) * | 2016-12-07 | 2018-03-05 | Федеральное государственное бюджетное учреждение науки Институт физики Земли им. О.Ю. Шмидта Российской академии наук (ИФЗ РАН) | Method of searching for mineral resources on the shelf of seas coated by ice |
| CN108680956A (en) * | 2018-01-08 | 2018-10-19 | 中国石油大港油田勘探开发研究院 | A kind of oil rich subdepression mature exploration area entirety exploitation method |
| CN109187302A (en) * | 2018-05-29 | 2019-01-11 | 中国石油大学(华东) | A kind of effective seal rock THICKNESS CALCULATION method based on cap rock heterogeneity |
| CN109709607A (en) * | 2018-12-07 | 2019-05-03 | 中国石油天然气股份有限公司 | A kind of prediction thin sandstone reservoirs thickness approach and device |
| CN111077578A (en) * | 2018-10-22 | 2020-04-28 | 中国石油天然气股份有限公司 | Rock stratum distribution prediction method and device |
| CN111680386A (en) * | 2020-04-17 | 2020-09-18 | 中国石油天然气股份有限公司 | Method and device for analyzing sealing performance of brittle cover layer of broken block trap |
| WO2021008630A1 (en) * | 2019-07-18 | 2021-01-21 | 中国石油天然气股份有限公司 | Method, apparatus, and system for quantitative analysis of sealing strength of trap fault |
-
2021
- 2021-02-05 CN CN202110160884.7A patent/CN112965108B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104063626A (en) * | 2014-07-11 | 2014-09-24 | 中国石油大学(华东) | Predication method of outer-source entrapment fullness degree |
| RU2646528C1 (en) * | 2016-12-07 | 2018-03-05 | Федеральное государственное бюджетное учреждение науки Институт физики Земли им. О.Ю. Шмидта Российской академии наук (ИФЗ РАН) | Method of searching for mineral resources on the shelf of seas coated by ice |
| CN108680956A (en) * | 2018-01-08 | 2018-10-19 | 中国石油大港油田勘探开发研究院 | A kind of oil rich subdepression mature exploration area entirety exploitation method |
| CN109187302A (en) * | 2018-05-29 | 2019-01-11 | 中国石油大学(华东) | A kind of effective seal rock THICKNESS CALCULATION method based on cap rock heterogeneity |
| CN111077578A (en) * | 2018-10-22 | 2020-04-28 | 中国石油天然气股份有限公司 | Rock stratum distribution prediction method and device |
| CN109709607A (en) * | 2018-12-07 | 2019-05-03 | 中国石油天然气股份有限公司 | A kind of prediction thin sandstone reservoirs thickness approach and device |
| WO2021008630A1 (en) * | 2019-07-18 | 2021-01-21 | 中国石油天然气股份有限公司 | Method, apparatus, and system for quantitative analysis of sealing strength of trap fault |
| CN111680386A (en) * | 2020-04-17 | 2020-09-18 | 中国石油天然气股份有限公司 | Method and device for analyzing sealing performance of brittle cover layer of broken block trap |
Non-Patent Citations (1)
| Title |
|---|
| 多参数岩性地震反演在盐岩地下储气库泥岩夹层分布中的应用;董娉婷等;《物探与化探》;20130430;第37卷(第02期);318-322+327 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112965108A (en) | 2021-06-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103603659B (en) | Utilize the method for Using Conventional Logs identification Complicated Pore Structures fluid type of reservoir through | |
| CN111722277B (en) | Method and device for describing broken solution trap for broken solution oil and gas reservoir exploration | |
| CN110632651B (en) | Method and device for quantitatively analyzing fault closure | |
| WO2013151455A1 (en) | Method for determining the water encroachment factor and the influx composition of an oil well | |
| CN105986813A (en) | Quasi tight reservoir rapid evaluation method and quasi tight reservoir multi-index evaluation method | |
| CN107272065A (en) | The disconnected solution profile testing method of carbonate rock | |
| CN106979006A (en) | The determination method and apparatus of strata pressure | |
| CN107861162B (en) | Microelectrode logging data-based natural crack identification method and system | |
| Aljuboori et al. | The impact of diagenesis precipitation on fracture permeability in naturally fractured carbonate reservoirs | |
| CN112965108B (en) | Method and system for determining vertical closure of trap cap layer | |
| CN115877464B (en) | Lithology recognition method and device, computer equipment and storage medium | |
| CN110308497B (en) | Method for fusing different types of reservoir body models of fracture-cavity oil reservoir | |
| CN110910267A (en) | Method and device for determining marine carbonate rock interlayer | |
| CN115522912A (en) | Cable adsorption card risk quantitative evaluation method, system and medium based on big data | |
| CN113126155B (en) | Sandstone reservoir prediction method aiming at strong reflection influence distributed among coal rocks | |
| Aslanyan et al. | Assessing macroscopic dynamic permeability through pressure and noise analysis | |
| CN110208861B (en) | Prediction method and device for constructing soft coal development area | |
| CN113534248A (en) | Method, device and system for quantitatively analyzing closure of trap fault | |
| CN114428365A (en) | Gas-containing characteristic-based tight sandstone gas reservoir logging identification method and device | |
| CN113294146A (en) | Method and apparatus for calculating effective porosity of asphaltene-containing reservoir | |
| CN115903026B (en) | Method, equipment and medium for analyzing composite sand body configuration | |
| US20240044239A1 (en) | A method of and apparatus for optimizing a drilling strategy | |
| CN117008185A (en) | Stratigraphic unconformity trap identification method and device, storage medium and electronic equipment | |
| RU2225020C1 (en) | Method of geophysical prospecting to determine oil productivity of fractured argillacous collectors in space between wells | |
| CN113052356A (en) | Oil well single well productivity prediction method and device, electronic equipment and storage medium |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |