CN114200540B - Method for dynamically predicting carbonate weathered shell karst reservoir development area - Google Patents
Method for dynamically predicting carbonate weathered shell karst reservoir development area Download PDFInfo
- Publication number
- CN114200540B CN114200540B CN202010985809.XA CN202010985809A CN114200540B CN 114200540 B CN114200540 B CN 114200540B CN 202010985809 A CN202010985809 A CN 202010985809A CN 114200540 B CN114200540 B CN 114200540B
- Authority
- CN
- China
- Prior art keywords
- target layer
- karst
- deposition
- discontinuities
- discontinuity
- 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
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000011161 development Methods 0.000 title claims abstract description 46
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 title claims abstract description 18
- 230000008021 deposition Effects 0.000 claims abstract description 64
- 238000012876 topography Methods 0.000 claims abstract description 43
- 238000005553 drilling Methods 0.000 claims abstract description 13
- 238000012360 testing method Methods 0.000 claims abstract description 11
- 230000015572 biosynthetic process Effects 0.000 claims description 26
- 230000015556 catabolic process Effects 0.000 claims description 11
- 238000006731 degradation reaction Methods 0.000 claims description 11
- 238000005516 engineering process Methods 0.000 claims description 7
- 230000002349 favourable effect Effects 0.000 claims description 4
- 238000011084 recovery Methods 0.000 claims description 4
- 238000002679 ablation Methods 0.000 claims description 3
- 230000009471 action Effects 0.000 claims description 3
- 230000035939 shock Effects 0.000 claims 1
- 239000011435 rock Substances 0.000 abstract description 11
- 230000008569 process Effects 0.000 abstract description 5
- 230000003872 anastomosis Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000009933 burial Methods 0.000 description 1
- 238000005137 deposition process Methods 0.000 description 1
- 235000021185 dessert Nutrition 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000002386 leaching Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a method for dynamically predicting a carbonate weathered shell karst reservoir development zone, which comprises the following steps: 1) Acquiring rock core, slice, logging, drilling, testing and three-dimensional seismic data; 2) Performing stratum division comparison, and determining a top boundary of a target layer, a bottom boundary of the target layer and a deposition discontinuity; 3) Explaining a top boundary of a target layer, a bottom boundary of the target layer and a deposition discontinuity; 4) Restoring karst paleo-topography (corresponding to multiple stages of the entire weathering stage, such as early, medium and end stages) between the bottom boundary of the target layer and adjacent depositional discontinuities, between adjacent depositional discontinuities, and between the top boundary of the target layer and adjacent depositional discontinuities; 5) And comprehensively predicting the development area of the weathered shell karst reservoir according to the recovered multiple karst palace appearances. According to the method, the evolution process of the karst palace and landform of the target layer in the whole weathering period is dynamically revealed, so that the prediction precision of the development area of the target layer weathering crust karst reservoir is improved.
Description
Technical Field
The invention relates to the technical field of karst reservoir distribution prediction in the petroleum industry, in particular to a method for dynamically predicting a carbonate weathered shell karst reservoir development area.
Background
Global oil and gas exploration practices have shown that karst is an important mechanism for the formation of high quality carbonate reservoirs. A plurality of large and medium-sized carbonate rock weathered karst reservoir oil and gas fields are found in a Tarim basin, an Erdos basin, a Sichuan basin and a plurality of foreign basins in China at present. Research shows that the distribution of karst reservoirs is closely related to the ancient karst landform, and the ancient karst slope and the ancient karst plateau edge are favorable development areas of the reservoirs, wherein karst hillocks are the areas where the reservoirs develop most. Therefore, accurately restoring karst palace is a key to carbonate rock weathered karst reservoir development area prediction.
In recent years, karst paleo-topography is mainly characterized by a layer sequence stratum method, a impression method, a residual thickness method and the like, so that a carbonate weathered shell karst reservoir development area is predicted. However, the stratum is degraded under the influence of the construction movement, and the karst paleo-landform restored by the residual thickness method has larger error with the actual paleo-landform. Meanwhile, the target layer may have a deposition break in the deposition process, and the stratum can also form a weathered crust karst reservoir after receiving the weathered leaching effect in the deposition break period. At present, the development area of the weathering crust karst reservoir predicts the karst paleo-topography of the whole target layer only to restore the karst paleo-topography of the whole target layer in a certain stage of the whole weathering period, is static, has relatively large error, and fails to dynamically reveal the karst paleo-topography evolution process of the target layer in a plurality of stages of the whole weathering period.
Therefore, in order to further improve the prediction precision, the invention provides a new thought of recovering the karst palace of a plurality of stages such as early stage (namely a first deposition intermittence period of a target layer), middle stage (namely a second deposition intermittence period of the target layer), end stage (namely a construction movement end stage) and the like of the weathering effect, and forms a new method for dynamically predicting the development area of the carbonate weathered shell karst reservoir.
Disclosure of Invention
Developing weathered karst reservoir predictive research is one of the carbonate reservoir research hotspots for a long time. The method for accurately predicting the development area of the weathered shell karst reservoir has very important significance for searching the dessert reservoir, improving the yield of a single well and saving the exploration investment. Therefore, in order to further improve the prediction precision of the distribution of the carbonate weathered crust karst reservoir, the invention provides a method for dynamically predicting the development area of the carbonate weathered crust karst reservoir.
In order to achieve the above purpose, the present invention adopts the following technical scheme:
a method of dynamically predicting carbonate weathered shell karst reservoir development zones, the method comprising the steps of:
1) Acquiring rock core, slice, logging, drilling, testing and three-dimensional seismic data;
2) Performing stratum division comparison, and determining a top boundary of a target layer, a bottom boundary of the target layer and a deposition discontinuity;
3) Explaining a top boundary of a target layer, a bottom boundary of the target layer and a deposition discontinuity;
4) Restoring karst paleo-topography (corresponding to multiple stages of the entire weathering stage, such as early, medium and end stages; when only one deposition discontinuity exists, the karst paleo-topography between adjacent deposition discontinuities does not exist, and the corresponding first and last stages of the entire weathering period are the same; of course, more than two deposition discontinuities may be present);
5) And comprehensively predicting the development area of the weathered shell karst reservoir according to the recovered multiple karst palace appearances.
According to the method, the evolution process of the karst palace and landform of the target layer in the whole weathering period is dynamically revealed, so that the prediction accuracy of the development area of the target layer weathering crust karst reservoir is improved.
Based on the method of the invention, preferably, the specific method for acquiring the rock core, the slice, the logging, the drilling, the testing and the three-dimensional seismic data in the step 1) is as follows: a) Visual rock sample characteristics such as lithology, pore size and the like are determined through photographing, visual observation and chemical analysis of the rock core; b) Obtaining data of mineral components, pore types and the like by observing the flakes under a microscope or an electron microscope; c) Obtaining information such as formation lithology composition, rock porosity, permeability, formation oil-gas property and the like by analyzing a logging curve; d) Analyzing formation lithology composition and hole seam development degree through rock scraps and drilling fluid usage degree obtained in the drilling process; e) Directly obtaining the characteristic of the stratum oil-gas by the oil test and gas test conditions; f) The plane spreading of the stratum is described through the explanation of the three-dimensional seismic data, and other information including the structural characteristics of the stratum, the burial depth of the stratum and the like is included;
based on the method of the invention, preferably, in the step 2), based on three-dimensional seismic data, the stratigraphic division comparison is carried out by adopting an interval stratigraphy method in combination with core, slice, logging, well drilling and test data.
Based on the method of the invention, preferably, the well earthquake combination in the step 2) is used for determining the top boundary of the target layer, the bottom boundary of the target layer and the deposition discontinuities by means of the analysis technology such as the same phase axis of the earthquake.
Preferably, the process of determining the top boundary of the destination layer, the bottom boundary of the destination layer and the deposition discontinuities by means of the seismic event isochronous analysis technology comprises the following steps of:
decomposing original seismic data into high-frequency data bodies and low-frequency data bodies, and calculating the inclination angle difference of the two data bodies to form an inclination angle difference data body; and comparing the inclination angle difference data body with the original seismic data, wherein the coincident seismic event is an equal time interface, namely a sedimentary discontinuity.
Based on the method of the invention, preferably, the seismic sequence interpretation technique is adopted in the step 3) to interpret the top boundary of the target layer, the bottom boundary of the target layer and the sedimentary discontinuities. The interpretation of the invention is needed to achieve fine interpretation, and the interpretation is performed on a line-by-line basis on three-dimensional seismic data.
Preferably, step 3) according to the method of the present invention comprises in particular: determining the seismic reflection characteristics of a top boundary of a target layer, a bottom boundary of the target layer and a deposition discontinuity according to the well earthquake fine calibration; and (3) comprehensively considering by adopting a seismic sequence interpretation technology according to a phase mark, an energy mark, a waveform mark and an interface reflection time difference stability comparison principle, and tracking and interpreting the top and bottom boundaries of a target layer and the layer positions of a deposition discontinuities.
Based on the method of the invention, preferably, in step 4), when recovering the karst paleo-topography between the top boundary of the target layer and the adjacent depositional discontinuities, the karst paleo-topography is characterized by the residual formation thickness for the regions of the formation that are not degraded;
for the region of the stratum suffering from the degradation, the karst paleo-topography is characterized according to the stratum deposit thickness of different stratum segments, which is calculated by the deposit inheritance in the same stratum.
The deposition discontinuities may include 1, 2, 3 or more based on the method of the present invention; preferably, the deposition discontinuities include a first deposition discontinuity and a second deposition discontinuity from bottom to top.
Preferably, the residual thickness between the top boundary of the target layer and the second deposition discontinuity is H, corresponding to the region of the formation not degraded, according to the method of the invention 3 The thickness of the stratum between the first deposition discontinuity and the second deposition discontinuity is H 2 ;
Corresponding to the region of the formation subject to degradation, the residual thickness between the top boundary of the target layer and the second deposition discontinuity is h 3 The formation degradation amount is Deltah 3 The thickness of the formation between the first deposition discontinuity and the second deposition discontinuity is h 2 ;
H is obtained by having deposit inheritance in the same layer according to the stratum deposit thickness of different layer segments 3 /H 2 =(h 3 +△h 3 )/h 2 Then (h) 3 +△h 3 )=H 3 ·h 2 /H 2 The method comprises the steps of carrying out a first treatment on the surface of the With (h) 3 +△h 3 ) To characterize the karst paleo-topography between the top boundary of the target layer and the second depositional discontinuity in the region of the formation subject to ablation.
Based on the method of the invention, preferably, the karst paleo-topography between the bottom boundary of the recovered target layer and the adjacent deposition discontinuities, between the adjacent deposition discontinuities and between the top boundary of the target layer and the adjacent deposition discontinuities correspond to the karst paleo-topography at the early stage, the middle stage and the end stage of the weathering respectively;
the step 5) specifically comprises the following steps: according to the principle that the paleo-topography high potential area is a favorable development area of the weathered crust karst reservoir, the weathered crust karst reservoir development areas of different stages are determined according to the karst paleo-topography recovery results of different stages in the whole weathered action period, on the basis, the weathered crust karst reservoir development areas of different stages are fused, and the distribution of the weathered crust karst reservoir development areas of the target layer is comprehensively predicted.
Compared with the static prediction method of the development area of the weathering crust karst reservoir in the prior art, the method fully considers the karst paleo-topography characteristics of multiple stages of the whole weathering period of the target layer, and simultaneously reduces the formation degradation amount (delta h in figure 2) 3 ) And the prediction accuracy of the development area of the weathered crust karst reservoir is greatly improved by recovering.
Drawings
FIG. 1 is a flow chart of a method of the present invention for dynamically predicting carbonate weathered shell karst reservoir development zones.
Fig. 2 is a schematic diagram of a karst paleo-geographic restoration technique at different stages of the weathering period.
FIG. 3 is a four-stage weathering early stage karst paleo-topography of a lamp assembly lamp recovered in an example.
Fig. 4 is a karst paleo-topography at the end of four-stage weathering of a lamp set lamp recovered in an example.
FIG. 5 is a plan view of a predicted four segment weathered karst reservoir development zone for a lamp array lamp according to an example.
Detailed Description
In order to more clearly illustrate the present invention, the present invention will be further described with reference to preferred embodiments. It is to be understood by persons skilled in the art that the following detailed description is illustrative and not restrictive, and that this invention is not limited to the details given herein.
As shown in fig. 1, the present invention provides a method for dynamically predicting carbonate weathered shell karst reservoir development zone, comprising the steps of:
1) Acquiring rock core, slice, logging, drilling, testing and three-dimensional seismic data;
2) Performing stratum division comparison, and determining a top boundary of a target layer, a bottom boundary of the target layer and a deposition discontinuity;
3) Explaining a top boundary of a target layer, a bottom boundary of the target layer and a deposition discontinuity;
4) Restoring karst paleo-topography (corresponding to multiple stages of the entire weathering stage, such as early, medium and end stages) between the bottom boundary of the target layer and adjacent depositional discontinuities, between adjacent depositional discontinuities, and between the top boundary of the target layer and adjacent depositional discontinuities;
5) And comprehensively predicting the development area of the weathered shell karst reservoir according to the recovered multiple karst palace appearances.
The present invention provides in part a preferred embodiment, the method comprising the steps of:
1) And acquiring core, slice, logging, drilling, testing and three-dimensional seismic data.
2) And (5) carrying out stratum division comparison, and determining a top boundary of the target layer, a bottom boundary of the target layer and a deposition discontinuity.
Based on three-dimensional seismic data, combining rock cores, slices, well logging, well drilling and test data, carrying out stratum division comparison by adopting a layer sequence stratigraphy method, combining well and earthquake, and determining a top boundary of a target layer, a bottom boundary of the target layer and a deposition discontinuity by means of an equal-time analysis technology (technical thought: decomposing original seismic data into high-frequency and low-frequency data bodies, calculating the inclination angle difference of the two data bodies to form an inclination angle difference data body, comparing the inclination angle difference data body with the original data body, and determining the coincident equal-time interfaces of the seismic phase bodies, namely the deposition discontinuities). As shown in FIG. 2, the deposition discontinuities include a first deposition discontinuity and a second deposition discontinuity from bottom to top.
3) The top boundary of the target layer, the bottom boundary of the target layer and the deposition discontinuities are explained.
Determining the seismic reflection characteristics of a top boundary of a target layer, a bottom boundary of the target layer and a deposition discontinuity according to the well earthquake fine calibration; by adopting a seismic sequence interpretation technology, the layers of the top, bottom and sedimentary discontinuities of the target layer are tracked and interpreted according to comparison principles such as phase marks (in-phase), energy marks (amplitude intensity and stability), waveform marks (waveform characteristic similarity), interface reflection time difference stability and the like.
4) Restoring karst paleo-topography (corresponding to multiple stages of the entire weathering stage, such as early, medium and end stages) between the bottom boundary of the target layer and adjacent depositional discontinuities, between adjacent depositional discontinuities, and between the top boundary of the target layer and adjacent depositional discontinuities;
as shown in fig. 2, an example includes two deposition discontinuities. In restoring the karst paleo-topography between the top boundary of the target layer and the adjacent sedimentary discontinuities (second sedimentary discontinuities), the residual layer thickness h of the target layer is used for the region of the stratum subject to degradation 3 Can not truly reflect the karst paleomorphology of the target layer at the end of weathering, can only use h 3 +△h 3 Can truly represent the karst paleo-topography. For the regions of the formation not degraded, the residual formation thickness H 3 Characterizing karst paleo-topography; for the region of the stratum suffering from the degradation, the karst paleo-topography is characterized according to the stratum deposit thickness of different stratum segments, which is calculated by the deposit inheritance in the same stratum.
Specifically, corresponding to the region of the formation not degraded, the residual thickness between the top boundary of the target layer and the second deposition discontinuity is H 3 (characterizing the end-stage paleomorphology of weathering), the thickness of the formation between the first and second depositional discontinuities is H 2 (characterizing mid-weathering paleo-topography), the formation thickness between the first sedimentary discontinuities and the bottom boundary of the target layer is H 1 (characterizing early paleomorphology of weathering). Corresponding to the region of the formation subject to degradation, a residual thickness between the top boundary of the target layer and the second deposition discontinuity of h 3 The formation degradation amount is Deltah 3 The thickness of the formation between the first deposition discontinuity and the second deposition discontinuity is h 2 (characterizing mid-weathering paleo-topography), a formation thickness between the first sedimentary discontinuities and the bottom boundary of the target layer of h 1 (characterizing early paleomorphology of weathering).
H is obtained by having deposit inheritance in the same layer according to the stratum deposit thickness of different layer segments 3 /H 2 =(h 3 +△h 3 )/h 2 Then (h) 3 +△h 3 )=H 3 ·h 2 /H 2 The method comprises the steps of carrying out a first treatment on the surface of the With (h) 3 +△h 3 ) To characterize the karst paleo-topography between the top boundary of the target layer and the second depositional discontinuity in the region of the formation subject to ablation.
5) And comprehensively predicting the development area of the weathered shell karst reservoir according to the recovered multiple karst palace appearances.
According to the principle that the paleo-topography high potential area is a favorable development area of the weathered crust karst reservoir, the weathered crust karst reservoir development areas of different stages are determined according to the karst paleo-topography recovery results of different stages in the whole weathered action period, on the basis, the weathered crust karst reservoir development areas of different stages are fused, and the distribution of the weathered crust karst reservoir development areas of the target layer is comprehensively predicted.
The embodiment of the invention adopts the existing static prediction method and the dynamic prediction method of the preferable scheme to comparatively predict the carbonate weathered shell karst reservoir development area aiming at the SC basin GS area. The weathering stage in this area includes a pre-stage and an end-stage.
Figure 3 shows a karst ancient plot of a certain stage (early stage) of recovery based on weathering. The high and low areas of the topography can be seen, and the areas with higher topography develop weathered crust karst reservoirs, so that the weathered crust karst reservoir distribution at a certain stage (early stage) based on the weathering can be determined. By comparing with the development degree of the reservoir in actual drilling in the graph, the anastomosis rate is lower: of the 14 wells drilled, 10 wells fit, with a fit rate of only 71%.
Figure 4 shows a karst ancient plot of a restoration at some stage (end stage) based on weathering. The high and low areas of topography can be seen, and the areas with higher topography develop weathered crust karst reservoirs, so that the distribution of the weathered crust karst reservoirs at a certain stage (end stage) of weathering can be determined. By comparing with the development degree of the reservoir in actual drilling in the graph, the anastomosis rate is lower: of the 14 wells drilled, 8 wells fit with only 57%.
FIG. 5 shows a weathered karst reservoir profile determined based on the dynamic restored karst paleo-topography of the present invention. In the figure, a dark gray region is an early karst reservoir development region, a light gray region is a terminal karst reservoir development region, and the dark gray region and the light gray region jointly form a target karst reservoir development region. By comparing with the development degree of the reservoir in actual drilling, the anastomosis rate is higher: of the 14 wells drilled, 13 wells were anastomosed, and the anastomosis rate reached 93%.
The exploration practice shows that the dynamic prediction method provided by the invention is feasible in the prediction of carbonate weathered shell karst reservoir development areas, and is worth being used as a reference in similar areas.
It should be understood that the foregoing examples of the present invention are provided merely for clearly illustrating the present invention and are not intended to limit the embodiments of the present invention, and that various other changes and modifications may be made therein by one skilled in the art without departing from the spirit and scope of the present invention as defined by the appended claims.
Claims (4)
1. A method of dynamically predicting carbonate weathered shell karst reservoir development zones, comprising the steps of:
1) Acquiring core, slice, logging, drilling, testing and three-dimensional seismic data,
2) Performing stratum division comparison, and determining a top boundary of a target layer, a bottom boundary of the target layer and a deposition discontinuity;
3) Explaining a top boundary of a target layer, a bottom boundary of the target layer and a deposition discontinuity;
4) Restoring karst paleo-topography between a bottom boundary of a target layer and an adjacent deposition discontinuity, between adjacent deposition discontinuities, and between a top boundary of the target layer and an adjacent deposition discontinuity by adopting a stratum thickness method;
5) Comprehensively predicting a development area of the weathered shell karst reservoir according to the recovered multiple karst palace appearances;
step 3) adopting a seismic sequence interpretation technology to interpret the top boundary of the target layer, the bottom boundary of the target layer and the deposition discontinuities;
the step 3) specifically comprises the following steps: determining the seismic reflection characteristics of a top boundary of a target layer, a bottom boundary of the target layer and a deposition discontinuity according to the well earthquake fine calibration; comprehensively considering according to phase marks, energy marks, waveform marks and interface reflection time difference stability comparison principles by adopting a seismic sequence interpretation technology, and tracking and interpreting the top and bottom boundaries of a target layer and the layer positions of a deposition discontinuities;
in the step 4), when the karst paleo-topography between the top boundary of the target layer and the adjacent deposition discontinuities is restored, the karst paleo-topography is characterized by the thickness of the residual stratum for the area of the stratum which is not degraded;
for the region of the stratum suffering from the degradation, representing the karst paleo-topography according to the stratum thickness calculated by the deposit inheritance of the stratum deposit thickness of different stratum sections in the same stratum;
the deposition discontinuities include a first deposition discontinuity and a second deposition discontinuity from bottom to top;
corresponding to the non-degraded area of the stratum, the residual thickness between the top boundary of the target layer and the second deposition gap is H 3 The thickness of the stratum between the first deposition discontinuity and the second deposition discontinuity is H 2 ;
Corresponding to the region of the formation subject to degradation, the residual thickness between the top boundary of the target layer and the second deposition discontinuity is h 3 The formation degradation amount is Deltah 3 The thickness of the formation between the first deposition discontinuity and the second deposition discontinuity is h 2 ;
H is obtained by having deposit inheritance in the same layer according to the stratum deposit thickness of different layer segments 3 /H 2 =(h 3 +△h 3 )/h 2 Then (h) 3 +△h 3 )=H 3 ·h 2 /H 2 The method comprises the steps of carrying out a first treatment on the surface of the With (h) 3 +△h 3 ) To characterize a karst paleo-topography between a top boundary of a target layer and a second depositional discontinuity of an area of the formation subject to ablation;
the karst paleo-topography between the restored bottom boundary of the target layer and the adjacent deposition discontinuities, between the adjacent deposition discontinuities and between the top boundary of the target layer and the adjacent deposition discontinuities correspond to the karst paleo-topography at the early stage, the middle stage and the end stage of the weathering respectively;
the step 5) specifically comprises the following steps: according to the principle that the paleo-topography high potential area is a favorable development area of the weathered crust karst reservoir, the weathered crust karst reservoir development areas of different stages are determined according to the karst paleo-topography recovery results of different stages in the whole weathered action period, on the basis, the weathered crust karst reservoir development areas of different stages are fused, and the distribution of the weathered crust karst reservoir development areas of the target layer is comprehensively predicted.
2. The method of claim 1, wherein the stratigraphic division comparison is performed in step 2) using a sequence stratigraphy method based on three-dimensional seismic data in combination with core, slice, log, well and test data.
3. The method of claim 2, wherein the top boundary of the destination layer, the bottom boundary of the destination layer and the depositional discontinuities are determined by means of the seismic event isochronous analysis technique in step 2) in conjunction with the well shocks.
4. The method of claim 3, wherein determining the top boundary of the destination layer, the bottom boundary of the destination layer, and the depositional discontinuities by means of seismic event isochronous analysis techniques comprises:
decomposing original seismic data into high-frequency data bodies and low-frequency data bodies, and calculating the inclination angle difference of the two data bodies to form an inclination angle difference data body; and comparing the inclination angle difference data body with the original seismic data, wherein the coincident seismic event is an equal time interface, namely a sedimentary discontinuity.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010985809.XA CN114200540B (en) | 2020-09-18 | 2020-09-18 | Method for dynamically predicting carbonate weathered shell karst reservoir development area |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010985809.XA CN114200540B (en) | 2020-09-18 | 2020-09-18 | Method for dynamically predicting carbonate weathered shell karst reservoir development area |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114200540A CN114200540A (en) | 2022-03-18 |
CN114200540B true CN114200540B (en) | 2024-01-30 |
Family
ID=80645310
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010985809.XA Active CN114200540B (en) | 2020-09-18 | 2020-09-18 | Method for dynamically predicting carbonate weathered shell karst reservoir development area |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114200540B (en) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109886457A (en) * | 2019-01-15 | 2019-06-14 | 中国石油大学(华东) | A method of judging the advantageous crust reservoirs of weathering development range of complex lithology buried hill |
CN110824557A (en) * | 2019-10-22 | 2020-02-21 | 中国石油天然气股份有限公司 | Quantitative karst ancient landform recovery method based on well-seismic combination |
CN111427087A (en) * | 2020-04-08 | 2020-07-17 | 中国石油天然气股份有限公司 | Method for determining carbonate reservoir position in ancient landform |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105626058B (en) * | 2015-12-30 | 2018-11-16 | 中国石油天然气股份有限公司 | A kind of method and device of determining reservoir karst development degree |
-
2020
- 2020-09-18 CN CN202010985809.XA patent/CN114200540B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109886457A (en) * | 2019-01-15 | 2019-06-14 | 中国石油大学(华东) | A method of judging the advantageous crust reservoirs of weathering development range of complex lithology buried hill |
CN110824557A (en) * | 2019-10-22 | 2020-02-21 | 中国石油天然气股份有限公司 | Quantitative karst ancient landform recovery method based on well-seismic combination |
CN111427087A (en) * | 2020-04-08 | 2020-07-17 | 中国石油天然气股份有限公司 | Method for determining carbonate reservoir position in ancient landform |
Non-Patent Citations (3)
Title |
---|
古地貌与不整合动态结合预测风化壳岩溶储集层分布-以四川盆地乐山—龙女寺古隆起灯影组为例;汤济广 等;《石油勘探与开发》;第40卷(第6期);第674-681页 * |
古地貌恢复在岩溶风化壳储层研究中的应用-以川中磨溪地区灯影组四段为例;朱茂 等;《海相油气地质》;第23卷(第4期);第87-95页 * |
地层剥蚀厚度恢复研究-以南黄海北部盆地东北凹陷为例;张海军 等;《中国石油勘探》(第5期);第21-25页 * |
Also Published As
Publication number | Publication date |
---|---|
CN114200540A (en) | 2022-03-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20150081265A1 (en) | Combining downhole fluid analysis and petroleum systems modeling | |
CN111046585B (en) | Shale gas dessert prediction method based on multiple linear regression analysis | |
CN108374657B (en) | Automatic well breakpoint identification method | |
CN109870721A (en) | A kind of method of sea area hydrate concentration prediction | |
CN113970796B (en) | Method for accurately recovering paleo-water depth of sedimentary basin | |
CN113608278A (en) | Sandstone-type uranium ore positioning method in red heterolayer of sedimentary basin | |
CN113093286B (en) | Inversion method for reservoir heterogeneity of cold spring development area | |
CN112505754B (en) | Method for collaborative partitioning sedimentary microfacies by well-seismic based on high-precision sequence grid model | |
CN114200540B (en) | Method for dynamically predicting carbonate weathered shell karst reservoir development area | |
CN111101935B (en) | Oil shale prediction method under few-well condition | |
Wang et al. | Seismic attributes for characterization and prediction of carbonate faulted karst reservoirs in the Tarim Basin, China | |
CN113311502B (en) | Method and device for identifying conventional oil layer and shale oil layer in shale layer system | |
CN110795513B (en) | Method for predicting distribution of river facies source storage ectopic type compact oil gas dessert area | |
US20240093597A1 (en) | Rapid Identification of Hydrodynamic Traps in Hydrocarbon Reservoirs | |
CN111708100B (en) | Deep thin layer oil and gas reservoir determination method and related device | |
Narimanov et al. | The Bahar oil and gas-condensate field in the South Caspian Basin | |
US11940398B2 (en) | Method to determine relative thermal maturity from porosities measured by quantitative imaging on uncleaned samples and GRI porosity measured on crushed cleaned samples | |
CN111983689B (en) | Near-source electromagnetic seismoelectric joint GR inversion method | |
CN112835098B (en) | Method and device for predicting energy storage coefficient of weathered-crust karst reservoir | |
US20230314646A1 (en) | Reservoir fluid mapping in mature fields | |
CN114200539A (en) | Karst ancient landform restoration method taking hydrothermal formation silicalite as isochronal datum plane | |
Sena et al. | Multiscale and multidisciplinary data-driven reservoir characterization of a fractured carbonate field in Kurdistan | |
Mendez et al. | Exploring La Luna Formation unconventional plays in the Maracaibo Lake Basin, through 3D modeling based on acoustic seismic inversion | |
Makra et al. | Classification of Miocene deposits using CPT data | |
Easow | Wellsite Real-Time Geochemical and Petrophysical Logging in North America Unconventionals |
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 |