CN114609675A - Quantitative recovery method for carbonate rock stratum sedimentary micro-landform based on high-frequency cycle - Google Patents
Quantitative recovery method for carbonate rock stratum sedimentary micro-landform based on high-frequency cycle Download PDFInfo
- Publication number
- CN114609675A CN114609675A CN202210246166.6A CN202210246166A CN114609675A CN 114609675 A CN114609675 A CN 114609675A CN 202210246166 A CN202210246166 A CN 202210246166A CN 114609675 A CN114609675 A CN 114609675A
- Authority
- CN
- China
- Prior art keywords
- sequence
- sedimentary
- frequency
- logging
- microfacies
- 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.)
- Pending
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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a quantitative restoration method for carbonate stratum sedimentary micro-landform based on high-frequency convolution, and relates to the technical field of ancient landform restoration. The method is based on sequence stratigraphy as a theoretical basis, fully applies basic data, adopts a spectrum analysis technology and a multi-scale matching method, establishes a seismic-logging sequence stratigraphic framework, analyzes the type and the development characteristics of a sequence stratigraphic interface, determines an isochronous geologic body, and obtains a high-frequency cycle combination type and characteristics, thereby realizing the quantitative recovery of carbonate rock stratum sedimentary micro-topography, not only improving the scale and the precision of the sedimentary topography, but also providing powerful support for deepening geological understanding and the next exploration direction, and further reducing the exploration risk and the cost investment.
Description
Technical Field
The invention relates to the technical field of petroleum exploration and development, in particular to a quantitative recovery method for deep carbonate rock stratum sedimentary micro-geomorphology based on high-frequency gyrus.
Background
With the continuous maturity of oil and gas exploration theory, more and more examples prove that the sedimentary environment is a prerequisite for the development of a deep carbonate reservoir, and a sedimentary facies provides a material basis for the formation of the reservoir, and the sedimentary facies and the reservoir together control the development and the spatial distribution of the reservoir. Based on the research of predecessors, the sedimentary environment and sedimentary facies are closely related to sedimentary micro-landform, so that the restoration of the sedimentary micro-landform is the key point for understanding the geological significance of the lithofacies ancient geography and oil gas thereof in the region. Nowadays, for the restoration of deposited geomorphic units with large geomorphic altitude difference, the traditional restoration method is difficult to finely carve such micro geomorphic units, such as the Hui-shui method (Daijingyou et al 2005), the deposition analysis method (Jun Xing et al 2001; Wu Li Yan et al 2005), the sequence stratigraphy method (fan Tai Liang et al 1999; Wang Jiahao et al 2003), and the like. In addition, seismic data or well logging data are often used for landform restoration, and the obtained result is difficult to convince due to lack of constraint and verification. Therefore, the recovery of the sedimentary micro-topography of the deep carbonate rock stratum is to be solved urgently, and the accuracy of the recovery is directly related to the next exploration deployment and decision.
Disclosure of Invention
In view of the above technical problems, the present invention aims to overcome the defects of the prior art and provide a method for quantitatively recovering sedimentary micro-topography of a deep carbonate rock stratum based on high-frequency cycle.
The invention adopts the following technical scheme that:
s1, establishing a stratum lattice of a stratum-logging curve based on the logging data;
s2, establishing a sequence stratigraphic framework by utilizing sequence stratigraphy based on the stratigraphic framework and combining the characteristics of the ancient creatures, the rock types, the lithofacies and the well logging curve;
s3, establishing a seismic-sequence stratigraphic framework based on the sequence stratigraphic framework and in combination with seismic data by considering the relationship between the earthquake and the sequence stratigraphic framework;
s4, establishing a high-frequency convolution combination based on core and slice data, establishing a high-frequency sequence grid by combining an earthquake-sequence stratum grid, and determining the development characteristics of the isochronous geologic body;
s5, based on the development characteristics of the isochronous geologic body, obtaining the deposition thickness of different microfacies according to the relationship between the deposited microfeatures and the microfacies differentiation;
and S6, based on the deposition thicknesses of different microfacies, combining geological background and deposition microfacies development, drawing contour maps of different deposition microfacies with the thickness of an isochronous geologic body as a base map, restoring the deposition microfeature distribution map of the isochronous geologic body by using a sequence stratigraphic method, and finely depicting the multilevel topographic units.
The invention has the beneficial effects that: the invention takes sequence stratigraphy as a theoretical basis, fully applies basic data, adopts a spectrum analysis technology and a multi-scale matching method, establishes a seismic-logging sequence stratigraphic framework, analyzes the type and the development characteristics of a sequence stratigraphic interface, determines an isochronous geologic body, and obtains the type and the characteristics of a high-frequency cycle combination, thereby realizing the quantitative recovery of carbonate rock stratum sedimentary micro-landform, not only improving the scale and the precision of the sedimentary landform, but also providing powerful support for deepening geological understanding and the next exploration direction, and further reducing exploration risk and cost input.
Detailed Description
The technical solutions of the present invention will be described in detail below in order to clearly understand the technical features, objects, and advantages of the present invention, but the present invention is not limited to the practical scope of the present invention.
A quantitative recovery method for deep carbonate stratum sedimentary micro-topography based on high-frequency gyrus comprises the following steps:
s1, establishing a stratum lattice of a stratum-logging curve based on the logging data;
specifically, the matching relation of lithology, lithofacies and logging curves is determined by combining logging data and depth correction and core homing, and a lithology-electrical recognition model and a stratum framework of the stratum-logging curves are established.
S2, establishing a sequence stratigraphic framework by utilizing sequence stratigraphy based on the stratigraphic framework and combining the characteristics of the ancient creatures, the rock types, the lithofacies and the well logging curve;
specifically, based on a stratum framework, analyzing the characteristics of the ancient creatures, rock types, lithofacies and logging curves of the stratum, wherein the logging curves comprise logging curves obtained by logging by utilizing gamma rays, neutrons and sound waves, summarizing the matching relation of the ancient creatures, rocks, lithofacies and logging curves of the stratum corresponding to each other, forming and establishing identification marks with unified four attributes of the ancient creatures, rocks, lithofacies and logging curves of the stratum, and providing multi-type data support for accurate division and checking of the stratum;
processing the logging curve by adopting a spectrum analysis technology to enable the formation convolution characteristic to be more obvious on the spectrum analysis curve, wherein in the prior art, various spectrum analysis methods exist, such as a traditional FFT algorithm, a ZoomFFT method, YIp-ZOOM conversion and the like, and can be applied to the invention;
and on the basis of the processed logging curve, dividing the stratum sequence by adopting a layer sequence stratigraphy dividing rule and method, determining the structure of the layer sequence, and establishing a layer sequence stratum framework of the stratum-logging curve.
S3, establishing a seismic-sequence stratigraphic framework based on the sequence stratigraphic framework and in combination with seismic data by considering the relationship between the earthquake and the sequence stratigraphic framework;
specifically, based on a stratigraphic sequence stratigraphic framework, combined with the research results of predecessors, according to the stratigraphic, the logging curve and the presented identification mark of seismic attributes, the identification mark or the corresponding mark of the stratigraphic and the logging curve on a seismic section is selected by utilizing seismic data, the logging curve is superposed on the seismic section, and the precise coincidence of the stratigraphic on the logging and the seismic two identification layers is realized by adopting the technical means of software automatic matching and manual fine adjustment through the comparison of two types of multiple point positions (two or more), so that the fine calibration of the well seismic data is realized;
secondly, analyzing and obtaining the development characteristics of the sequence on the seismic section according to the sequence structure of the logging curve;
thirdly, according to the development characteristics of the sequence on the seismic section, combining a sequence stratum framework to establish a seismic-sequence stratum framework;
and finally, identifying and determining the spreading characteristics of the sequence interfaces in the transverse direction and the longitudinal direction according to the earthquake-sequence stratigraphic framework by combining the sequence stratigraphic principle and the interface identification marks thereof, judging the types of the interfaces of different sequence, further analyzing and obtaining the stacking sequence and the contact relation of the interfaces of different sequence in the vertical direction and the transverse direction, and determining the deposition process and the deposition environment.
S4, establishing a high-frequency convolution combination based on core and slice data, establishing a high-frequency sequence grid by combining an earthquake-sequence stratum grid, and determining the development characteristics of the isochronous geologic body;
specifically, firstly, based on core and slice data, according to a carbonate rock classification scheme (the classification scheme can refer to the following documents of limestone classification, Folk, 1962, carbonate deposit structure classification, R.J.Dunham, 1962, carbonate classification, Von Zeng Sho, 1982), the texture and development characteristics of carbonate rock are observed and identified mainly from three aspects of color, mineral components, structure and content thereof, the development types and development characteristics of rock are summarized, the development characteristics, the deposition processes and environments of the sequence in the vertical direction and the transverse direction are combined, the development characteristics of different rock types in the vertical direction are analyzed and obtained, and the high-frequency gyrus type and the stacking combination relationship are determined;
secondly, based on the earthquake-sequence stratigraphic framework and high-frequency cycle combination characteristics, selecting a well bore capable of representing the rock, sedimentary facies and stratigraphic characteristics as a representative well according to the geological characteristics of the region totally reflected by the stratigraphic, rock, logging curve and earthquake data and well bit plane distribution, performing spectral analysis and wavelet transformation on the logging curve representing the well, performing time-frequency spectral analysis on seismic data of a well side channel, combining the data characteristics of the stratigraphic, rock, logging curve and earthquake under four different scales, and adopting a multi-scale matching method of well-earthquake combination to establish the high-frequency cycle framework and determine the development characteristics and the transverse and longitudinal spread characteristics of a high-frequency cycle interface;
and finally, checking and correcting the high-frequency convolution interface based on the high-frequency convolution grid by combining the rock core data, the sedimentary facies and the identification marks of the previous convolution interface, determining the geologic body between the two high-frequency convolution interfaces as an isochronous geologic body, analyzing and obtaining the distribution range and the spreading characteristics in the transverse direction of the isochronous geologic body, and drawing an isochronous geologic body thickness plane distribution map.
S5, based on the development characteristics of the isochronous geologic body, obtaining the deposition thickness of different microfacies according to the relationship between the deposited microfeatures and the microfacies differentiation;
specifically, firstly, selecting an isochronous geologic body based on a high-frequency sequence grid, identifying the rock type of the isochronous geologic body by using a lithology-electrical property identification model, determining the sedimentary microfacies type of the isochronous geologic body according to sedimentary facies classification standards, constructing a sedimentary evolution mode by combining sedimentary geological backgrounds, analyzing and obtaining sedimentary microfacies types and characteristics mainly from four aspects of rock types, development positions, distribution ranges and thicknesses, and determining identification marks of different sedimentary microfacies;
secondly, establishing a logging facies interpretation model of different sedimentary microfacies based on the lithology-electrical property identification model and the sedimentary microfacies types and characteristics;
finally, the deposition thicknesses of different microphases are respectively counted according to the coring data and the logging interpretation model, wherein the coring section can be directly counted, and the coring single well can obtain the deposition thicknesses by using the logging interpretation model so as to draw the planar distribution map of the thickness of the different microphases.
S6, based on the deposition thicknesses of different microfacies, in combination with geological background and deposition microfacies development, drawing contour maps of different deposition microfacies with the thickness of an isochronous geologic body as a base map, recovering the deposition microfeatus distribution map of the isochronous geologic body by using a sequence stratigraphic method, and finely depicting multi-stage topographical units;
specifically, firstly, analyzing and obtaining microphase plane distribution characteristics from five aspects of microphase types, rock types, thicknesses, positions and ranges based on different sedimentary microphase thickness plane distribution maps and combining sedimentary microphase longitudinal development characteristics and tectonic development characteristics, and determining the overall sedimentary geological background and sedimentary environment development characteristics;
secondly, restoring the sedimentary micro-geomorphic distribution map of the selected isochronous geologic body by combining the overall sedimentary geologic background and the sedimentary environment development characteristics and utilizing a sedimentary micro-geomorphic restoration principle and adopting a sequence stratigraphy to obtain the overall sedimentary micro-geomorphic distribution map;
and finally, superposing the equal-time geologic body thickness plane distribution map and the deposition micro-landform distribution map, finely depicting the multi-stage landform units by combining the research result of the rock facies distribution rule of the predecessor and the landform unit division standard and basis, and summarizing the characteristics of the rock types, the deposition micro-facies and the distribution ranges of different deposition micro-landform units.
Through the steps, quantitative recovery of the sedimentary micro-landform of the carbonate rock stratum can be achieved, powerful support is provided for deepening geological understanding and the next exploration direction, and exploration risks and cost investment can be reduced.
The above description is only a preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the embodiments of the present invention are included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (7)
1. A quantitative recovery method for deep carbonate stratum sedimentary micro-topography based on high-frequency gyrus is characterized by comprising the following steps:
s1, establishing a stratum lattice of a stratum-logging curve based on the logging data;
s2, establishing a sequence stratigraphic framework by utilizing sequence stratigraphy based on the stratigraphic framework and combining the characteristics of the ancient creatures, the rock types, the lithofacies and the well logging curve;
s3, establishing a seismic-sequence stratigraphic framework based on the sequence stratigraphic framework and in combination with seismic data by considering the relationship between the earthquake and the sequence stratigraphic framework;
s4, establishing a high-frequency convolution combination based on core and slice data, establishing a high-frequency sequence grid by combining an earthquake-sequence stratum grid, and determining the development characteristics of the isochronous geologic body;
s5, based on the development characteristics of the isochronous geologic body, obtaining the deposition thickness of different microfacies according to the relationship between the deposited microfeatures and the microfacies differentiation;
and S6, based on the deposition thicknesses of different microfacies, combining geological background and deposition microfacies development, drawing contour maps of different deposition microfacies with the thickness of an isochronous geologic body as a base map, restoring the deposition microfeature distribution map of the isochronous geologic body by using a sequence stratigraphic method, and finely depicting the multilevel topographic units.
2. The method as claimed in claim 1, wherein the specific operation of S1 is to establish a lithology-electrical recognition model and a stratigraphic framework of stratigraphic-logging curves by determining matching relationships of lithology, lithofacies and logging curves through depth correction and core homing in combination with logging data.
3. The method of claim 1, wherein the specific operation of S2 is to analyze the ancient creatures, rock types, lithofacies and well log characteristics of the stratum based on the stratum lattice to clarify the coupling relationship; and establishing a sequence stratigraphic framework by adopting a classification rule and a method of sequence stratigraphy and combining a frequency spectrum analysis method.
4. The method of claim 1, wherein the specific operation of S3 is to perform fine well seismic calibration based on sequence stratigraphic framework in combination with seismic data, and analyze the relationship between the seismic and sequence stratigraphic framework; and establishing a seismic-sequence stratigraphic framework according to the relation between the earthquake and the sequence stratigraphic layer, and analyzing the stacking relation and the development characteristics of the sequence interface.
5. The method as claimed in claim 1, wherein the specific operation of S4 is that the high frequency convolution is determined by observing and identifying the organization and development characteristics of the core and slice data; the specific operation of establishing the high-frequency sequence grid is that based on the earthquake-sequence stratigraphic grid and the high-frequency cycle combination characteristics, a representative well is selected, the logging curve of the representative well is subjected to spectrum analysis and wavelet transformation, time-frequency spectrum analysis is carried out on well side channel earthquake data, a well-earthquake combined multi-scale matching method is adopted to establish the high-frequency sequence grid, and after the high-frequency sequence grid is established, the high-frequency sequence interface is corrected by combining rock core data, sedimentary facies and development cycle characteristics thereof, so that the development characteristics of the isochronous geologic body are determined.
6. The method as claimed in claim 1, wherein the specific operation of S5 is to select an isochronous geologic body based on a high frequency sequence grid, construct a sedimentary evolution model in combination with a sedimentary geological background, analyze the relationship between sedimentary microfeatures and sedimentary microfacies differentiation, clarify the developmental features and signatures of different sedimentary microfacies, and establish a conventional well log facies and an imaging well log facies of different sedimentary microfacies; secondly, respectively counting the deposition thickness of different microphases according to the coring data, the conventional logging phase and the imaging logging phase.
7. The method according to claim 1, wherein the specific operations of S6 are: the specific operation of drawing the sedimentary micro-geomorphic profile is to analyze the spreading characteristics of different microphase planes by combining the sedimentary microfacies longitudinal development characteristics and the tectonic development characteristics based on different sedimentary microphase thicknesses, and recover the sedimentary micro-geomorphic profile of the selected isochronous geological body by utilizing a sedimentary microfeature recovery principle and adopting a sequence stratigraphy; and on the basis of depositing the micro-landform distribution diagram, combining with planar rock spread distribution to finely depict the multi-stage landform units.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210246166.6A CN114609675A (en) | 2022-03-14 | 2022-03-14 | Quantitative recovery method for carbonate rock stratum sedimentary micro-landform based on high-frequency cycle |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210246166.6A CN114609675A (en) | 2022-03-14 | 2022-03-14 | Quantitative recovery method for carbonate rock stratum sedimentary micro-landform based on high-frequency cycle |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114609675A true CN114609675A (en) | 2022-06-10 |
Family
ID=81863890
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210246166.6A Pending CN114609675A (en) | 2022-03-14 | 2022-03-14 | Quantitative recovery method for carbonate rock stratum sedimentary micro-landform based on high-frequency cycle |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114609675A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115450611A (en) * | 2022-09-16 | 2022-12-09 | 中国地质大学（北京） | Deep carbonate rock sedimentary microphase analysis method based on random forest |
| CN115903047A (en) * | 2022-12-27 | 2023-04-04 | 中国地质调查局油气资源调查中心 | Method and device for identifying marine carbonate beach facies reservoir |
| WO2024087800A1 (en) * | 2022-10-27 | 2024-05-02 | 中国石油天然气股份有限公司 | Conventional well logging data-based method for identifying high-frequency cyclic carbonate rock |
-
2022
- 2022-03-14 CN CN202210246166.6A patent/CN114609675A/en active Pending
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115450611A (en) * | 2022-09-16 | 2022-12-09 | 中国地质大学（北京） | Deep carbonate rock sedimentary microphase analysis method based on random forest |
| WO2024087800A1 (en) * | 2022-10-27 | 2024-05-02 | 中国石油天然气股份有限公司 | Conventional well logging data-based method for identifying high-frequency cyclic carbonate rock |
| CN115903047A (en) * | 2022-12-27 | 2023-04-04 | 中国地质调查局油气资源调查中心 | Method and device for identifying marine carbonate beach facies reservoir |
| CN115903047B (en) * | 2022-12-27 | 2023-07-04 | 中国地质调查局油气资源调查中心 | Method and device for identifying sea carbonate beach reservoir |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109061765B (en) | Trap evaluation method for heterogeneous thin sandstone interbed reservoir | |
| US5838634A (en) | Method of generating 3-D geologic models incorporating geologic and geophysical constraints | |
| US7113869B2 (en) | Method for prospect identification in asset evaluation | |
| US7769546B2 (en) | Method for indexing a subsurface volume for the purpose of inferring geologic information | |
| CN114609675A (en) | Quantitative recovery method for carbonate rock stratum sedimentary micro-landform based on high-frequency cycle | |
| CN110824557B (en) | Quantitative karst ancient landform recovery method based on well-seismic combination | |
| US10732310B2 (en) | Seismic attributes derived from the relative geological age property of a volume-based model | |
| CN103993871A (en) | Method and device for processing well logging information of thin interbed stratums in standardization mode | |
| CN104678434A (en) | Method for predicting storage layer crack development parameters | |
| CN103775075A (en) | Method for identifying lithology of whole well section | |
| CN113031068B (en) | Reflection coefficient accurate base tracking prestack seismic inversion method | |
| CN105842733A (en) | Shale reservoir earthquake identification method | |
| CN112505754B (en) | Method for collaborative partitioning sedimentary microfacies by well-seismic based on high-precision sequence grid model | |
| CN112946782B (en) | Earthquake fine depicting method for dense oil-gas storage seepage body | |
| GB2583910A (en) | Method of analysing seismic data | |
| CN109283577B (en) | Seismic horizon calibration method | |
| CN107132574B (en) | Prediction method for lithologic combination of marine stratum clastic rocks | |
| CN114861515A (en) | Method, device, equipment and medium for calculating layer speed data volume | |
| Wethington et al. | Stratigraphic architecture of the Mississippian limestone through integrated electrofacies classification, Hardtner field area, Kansas and Oklahoma | |
| CN112946751A (en) | Method and system for predicting plane distribution of different reservoirs | |
| CN113484907B (en) | Method for predicting distribution on different types of reservoir planes | |
| CN115857047A (en) | Comprehensive prediction method for seismic reservoir | |
| CN116027454A (en) | Method for imaging contour map of sandstone thickness and intelligent drawing device | |
| CN110389381B (en) | Sand reservoir prediction method and device based on seismic attributes | |
| CN113514884A (en) | Compact sandstone reservoir prediction method |
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 |