CN106443772B - One kind going diapir original stratigraphic thickness restoration methods - Google Patents
One kind going diapir original stratigraphic thickness restoration methods Download PDFInfo
- Publication number
- CN106443772B CN106443772B CN201610784447.1A CN201610784447A CN106443772B CN 106443772 B CN106443772 B CN 106443772B CN 201610784447 A CN201610784447 A CN 201610784447A CN 106443772 B CN106443772 B CN 106443772B
- Authority
- CN
- China
- Prior art keywords
- stratum
- thickness
- salt
- rock
- diapir
- 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 28
- 235000002639 sodium chloride Nutrition 0.000 claims abstract description 60
- 150000003839 salts Chemical class 0.000 claims abstract description 58
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 19
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 claims abstract description 12
- 239000011780 sodium chloride Substances 0.000 claims abstract description 12
- 238000005553 drilling Methods 0.000 claims abstract description 4
- 239000011435 rock Substances 0.000 claims description 41
- 238000006243 chemical reaction Methods 0.000 claims description 11
- 230000004044 response Effects 0.000 claims description 3
- 238000005755 formation reaction Methods 0.000 abstract description 16
- 229910052602 gypsum Inorganic materials 0.000 description 7
- 239000010440 gypsum Substances 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 238000011084 recovery Methods 0.000 description 4
- 239000004576 sand Substances 0.000 description 4
- 230000007480 spreading Effects 0.000 description 4
- 230000008021 deposition Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 208000035126 Facies Diseases 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
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. for interpretation or for event detection
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses one kind to go diapir original stratigraphic thickness restoration methods, belongs to basin analysis and oil-gas geology field.This method mainly includes the following steps: that a. obtains formation thickness on the basis of high accuracy surface modelling data interpretation；B. rock salt plane distributing scope is determined using seismic data and drilling data；C. the adjacent earth formations volume that rock salt stratum is not plastically deformed with co-extensive is calculated, and finds out its ratio relation；D. using adjacent earth formations thickness spread as constraint condition, restore the original thickness on diapir stratum.The present invention predicts the stratum restoration methods relative to forefathers, and precision of prediction with higher, invention is applied in the multiple salt diapir development areas of In Jianghan Basin, practical.
Description
Technical Field
The invention relates to the field of basin analysis and oil-gas geology, in particular to a method for recovering the thickness of an original stratum of a diapir is removed.
Background
Salt structure is a structure type which is closely related to oil gas accumulation, and salt structure research is also one of hot spots and difficult problems of oil gas-containing basin structure research in the world at present. The development of salt layers in different degrees from the frigid system to the fourth system in the crust, the deformation of salt bodies and the influence thereof on the distribution of sedimentary facies zones and sand bodies can form different structural encirclements and formation-lithologic encirclement. Most of the over 150 salty basins worldwide contain abundant oil and gas resources, such as the famous gulf of mexico, the north sea, the bos gulf, the northern europe, and the continental africa. The research on salt structure in foreign countries is earlier developed, and great progress is made in the past decades, and the research results prove that the salt structure research has wide application prospect for oil and gas exploration and development and promote the exploration breakthrough of salt structure related oil and gas reservoirs. China also develops a multi-stage salt rock stratum and various salt structures widely in basins such as Tarim, Bohai Bay, Sichuan and Jianghan, but the research on the salt structures in China is relatively less and is still in a qualitative description stage. With the increasing demand for oil and gas, the research on salt formation is urgently needed.
Under the action of different geological conditions, rock salt generates plastic flow, structural deformation characteristics of surrounding stratums are influenced, salt structural deformation patterns with different forms can be formed finally, and the complex interaction process of the rock salt and an overburden layer in the salt structural growth process is reflected. Salt structures generally evolve from less extensive and mature integrated structures to more extensive and mature unconformity structures, gradually forming an invasive pioneering body. The original formation thickness recovery is typically accomplished by adding the residual thickness to the denuded thickness to approximately reflect the original spread of the formation. In a rigid stratum development area, the original spreading state of the stratum can be effectively recovered by using a traditional original stratum thickness recovery method, in a salt diapir development area, certain limiting conditions are useless for the salt diapir, the salt stratum is easy to generate plastic flow in the process of being stressed, and the salt structure generally forms three-dimensional flow of salt rocks, so that the inflow and outflow amount in the deformation process is difficult to know, and the thickness is difficult to control by using a conventional original thickness recovery method due to the influence of factors such as the deformation irregularity and the like, so that a large error exists, even an error exists.
Just because it is often difficult and uncertain to restore the original thickness of the salt diapir development area, in the process of the former restoration, the thickness abnormality caused by diapir effect is usually hidden deliberately, or the diapir structure is treated by adopting a thickness sharing mode, but the salt rock stratum is not deposited with equal thickness when being deposited, so the salt structure is restored on a small-scale basis, the restoration is continued, and the three-dimensional data of the area is analyzed fully.
Disclosure of Invention
The invention aims to solve the technical problem of providing a method for restoring the thickness of an original stratum by removing a diapir.
The technical scheme adopted by the invention for solving the technical problems is as follows: a method for recovering the thickness of an original stratum of a diapir is disclosed, which comprises the following steps:
1) obtaining stratum horizon interpretation stratum thickness according to high-precision three-dimensional seismic data interpretation;
2) determining the plane distribution range of the salt rock stratum by using the drilling data and the seismic data;
3) calculating the volume of the salt rock stratum and the adjacent stratum in the same range as the salt rock stratum according to the plane distribution range of the salt rock stratum; selecting a formation closest to the adjacent formation without plastic deformation;
4) restoring the original thickness of the diaplasis stratum by using the thickness distribution of the adjacent stratum without plastic deformation as a constraint condition;
let the original thickness of the rock salt formation be H1The thickness of the adjacent non-salt formation is H2Volume of rock salt formation is V1Volume of adjacent strata is V2，
Then H1＝H2*V1/V2。
The invention has the following beneficial effects: the invention provides a method for recovering the thickness of an original stratum of a diapir removing machine, namely a method for recovering the thickness of the original stratum of the diapir removing machine under geological constraint conditions. The invention is applied to stratum thickness recovery of salt diapir development areas of basin in Jianghuang, and has strong practicability.
Drawings
The invention will be further described with reference to the accompanying drawings and examples, in which:
FIG. 1 is a flow chart of a method of an embodiment of the present invention;
FIG. 2 is a schematic diagram of the restoration of original formation thickness by the method of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail with reference to the following embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
As shown in fig. 1 and fig. 2, the method comprises the following specific steps:
1) obtaining stratum horizon interpretation according to high-precision three-dimensional seismic data interpretation, wherein the density of seismic interpretation grids is generally more than 20m multiplied by 20m for ensuring the precision; in this example 50m × 100 m;
2) determining the distribution range of the rock salt plane by using logging data and seismic data;
the well logging response identification means that the well logging curves of the gypsum rock show great difference, and is a basic basis for identifying the gypsum rock. Particularly, the GR curve can directly identify the gypsum rock at a drilling position, and the gypsum rock and the thickness can be accurately determined; the seismic data are mainly used for identifying the gypsum rock in space, the density and seismic wave propagation speed of the gypsum rock are obviously different from those of sand and mudstone, and a stronger seismic reflection interface can be formed by generating a larger wave impedance difference with surrounding rock, so that the plane distribution range of the gypsum rock can be accurately identified and depicted by analyzing seismic response characteristics and utilizing seismic attributes; and the distribution of the gypsum-rock is predicted by integrating the well-seismic calibration and the well-seismic calibration information, and the area of the gypsum-rock is calculated with the accuracy rate of 99%.
3) And calculating the volume of the salt rock stratum and the adjacent stratum in the same range. The distribution range of the salt rock stratum is obtained in the step 2, the volume of the salt rock stratum in the diapir development area is calculated, and meanwhile, the volume of the adjacent stratum in the same range is calculated. The volume calculation was performed in the software buffer, in order: and (3) exporting horizon data in the geoframe, carrying out time-depth conversion after gridding, subtracting the grid data subjected to the time-depth conversion to obtain stratum thickness data, importing the area obtained in the step (2), and calculating the volume of the stratum.
Time depth conversion formula: 299.11t2+902.93t +63.137(t in seconds and y in meters)
Taking the four lower strata in the laten north area as an example, the time-depth conversion table of the T7 interface is:
TABLE 1 time-depth conversion table for T7 interface in the northwest of the basin in Jianghuang province
The time-depth transition table of the T6' interface is:
TABLE 2 time-depth conversion table for T6' interface in the northwest of the basin in Jianghuang province
The thickness data of the four potential formations is:
TABLE 2 thickness gauge of four submerged areas in the deep north of the basin in Jianghan
And importing the thickness data into a buffer, and adding boundary conditions to obtain the volume of the stratum.
4) And recovering the original thickness of the diaplasis stratum by using the thickness distribution of the adjacent stratum without plastic deformation as a constraint condition. If no large tectonic activity occurs during the deposition period of the salt rock stratum and the adjacent stratum, the spreading state of the adjacent stratum is approximately consistent with that of the salt rock stratum, so that the volume ratio of the salt rock stratum to the adjacent stratum is the thickness ratio of the salt rock stratum to the adjacent stratum under the condition of the same plane projection area. Although the rock salt flows through the shaping and the plane spreading state and the state of the rock salt during the sedimentation period are greatly changed due to the later diapir effect, the total volume of the rock salt is not changed all the time. In the process of forming the diapir structure, the plane spreading state of the non-salt stratum will not change greatly, so the original distribution state of the salt rock stratum can be reversely deduced by the distribution state of the non-salt stratum at present. Obtaining the volume V of the salt rock stratum through the step 31And volume V of adjacent formation2(the upper stratum or the lower stratum of the salt rock stratum is selected according to the construction background and the lithology of the stratum), and the ratio relation n is obtained as V1/V2That is, the volume of the salt rock stratum in the same range is n times of the volume of the adjacent non-salt system stratum, and under the condition that the original sedimentary grid and the plane projection area are the same, the original thickness H of the salt rock stratum1I.e. the thickness H of the non-salt stratum2N times of (i.e. H)1＝n×H2。
Example 1:
the invention is applied to the back inclined belt of the queen yard in the submarine north area. The method mainly comprises the following steps:
according to the three-dimensional seismic data of the potential north area, time depth data of T7, T6 'and T6 interfaces are obtained, then a time depth conversion formula of a standard well fitting basin is selected, the time depth is converted into the real depth, and the thicknesses of the strata of the potential four lower sections and the potential four upper sections are obtained by utilizing T7-T6' and T6-T6. Obtaining the distribution range of the salt rock through well data and seismic data, and obtaining the volume V of the salt rock stratum in the buffer1Comprises the following steps: 2.93829123203X 1011m3Volume V of the same range of upper four strata2Is 1.13133438322 multiplied by 1011m3Therefore, the ratio relationship is:
n＝V1/V2＝2.597190782505
the thickness H of the upper section of the submerged formation2Enlargement by a factor of n, i.e. H1＝H2’＝n×H2Re-imaging to obtain a volume of 2.91694838194 × 1011m3. And replacing the thickness data of the diapir development area with the new thickness data to re-map. After the original stratum thickness is recovered by the method, the salt rock stratum remains unchanged before and after, and the stratum distribution state is consistent with the upper section of the submarine four, thereby basically reflecting the original distribution state during stratum deposition.
Example 2:
the invention is applied to Jiangling sunken temple areas. The method mainly comprises the following steps:
according to three-dimensional seismic data of Jiangling areas, time depth data of T10, T9 and T8 interfaces are obtained, then a time depth conversion formula of a standard well fitting basin is selected, the time depth is converted into real depth, and stratum thicknesses of a sand city group and a new ditch group are obtained by utilizing T10-T9 and T9-T8. Obtaining the distribution range of the salt rock according to well data and seismic data, and obtaining the volume V of the stratums of the sand city in the buffer1Comprises the following steps: 6.27624469709X 1011m3Volume V of new trench group formation in the same range2Is 3.35257327128 multiplied by 1011m3Therefore, the ratio relationship is:
n＝V1/V2＝1.8720678682419
the thickness H of the new groove group2Enlargement by a factor of n, i.e. H1＝H2’＝n×H2Re-imaging to obtain a volume of 6.14363586425 × 1011m3. And replacing the thickness data of the diapir development area with the new thickness data to re-map. After the original stratum thickness is recovered by the method, the salt rock stratum remains unchanged before and after the salt rock stratum, and the original distribution state during stratum deposition is basically reflected.
It will be understood that modifications and variations can be made by persons skilled in the art in light of the above teachings and all such modifications and variations are intended to be included within the scope of the invention as defined in the appended claims.
Claims (1)
1. A method for recovering the thickness of an original stratum of a diapir is characterized by comprising the following steps:
1) obtaining a stratum horizon explanation according to the high-precision three-dimensional seismic data explanation;
2) determining the plane distribution range of the salt rock stratum by using the drilling data and the seismic data; the method comprises the following specific steps: forecasting the distribution of the salt rock by analyzing seismic response characteristics and calibrating well-to-well seismic, calculating the area of the distribution, and accurately identifying and depicting the plane distribution range of the salt rock stratum;
3) respectively calculating the volume of the salt rock stratum and the volume of the adjacent stratum in the same range with the salt rock stratum according to the plane distribution range of the salt rock stratum; selecting a formation closest to the adjacent formation without plastic deformation;
the method comprises the following specific steps: exporting horizon data in the geoframe, carrying out time-depth conversion after gridding, subtracting the grid data subjected to the time-depth conversion to obtain stratum thickness data, importing the area obtained in the step 2), and calculating the volume of the stratum; wherein,
time depth conversion formula: y =299.11t2+902.93t +63.137, t being time and y being depth;
4) restoring the original thickness of the diaplasis stratum by using the thickness distribution of the adjacent stratum without plastic deformation as a constraint condition;
let the original thickness of the rock salt formation be H1The thickness of the adjacent non-salt formation is H2Volume of rock salt formation is V1Volume of adjacent strata is V2，
Then H1= H2*V1/ V2。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610784447.1A CN106443772B (en) | 2016-08-31 | 2016-08-31 | One kind going diapir original stratigraphic thickness restoration methods |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610784447.1A CN106443772B (en) | 2016-08-31 | 2016-08-31 | One kind going diapir original stratigraphic thickness restoration methods |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106443772A CN106443772A (en) | 2017-02-22 |
| CN106443772B true CN106443772B (en) | 2018-12-04 |
Family
ID=58091252
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610784447.1A Active CN106443772B (en) | 2016-08-31 | 2016-08-31 | One kind going diapir original stratigraphic thickness restoration methods |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106443772B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111045109B (en) * | 2018-10-11 | 2023-06-02 | 中国石油化工股份有限公司 | Method for recovering stratum erosion amount in basin |
| CN112241578B (en) * | 2019-07-18 | 2022-10-04 | 中国石油天然气股份有限公司 | Salt structure style analysis method, device and system |
| CN113607920B (en) * | 2021-07-30 | 2023-10-24 | 中海石油(中国)有限公司深圳分公司 | Reconstruction analysis method, experimental device and medium for sedimentary basin by magma bottom wall |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| UA80043C2 (en) * | 2005-11-15 | 2007-08-10 | Ukrainian State Geol Explorati | Method for predicting hydrocarbon deposits |
| RU2484503C1 (en) * | 2011-12-12 | 2013-06-10 | Закрытое акционерное общество "Актуальная геология" | Method for helium survey on water bodies |
| CN105204069B (en) * | 2014-06-27 | 2018-08-17 | 中国石油化工股份有限公司 | A kind of Eroded Thickness restoration methods |
| CN105425316A (en) * | 2015-10-28 | 2016-03-23 | 中国石油化工股份有限公司 | Ancient landform restoration method based on impression method |
| CN105388526B (en) * | 2015-10-28 | 2020-10-23 | 中国石油化工股份有限公司 | Ancient geological structure recovery method for continental lake basin |
-
2016
- 2016-08-31 CN CN201610784447.1A patent/CN106443772B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN106443772A (en) | 2017-02-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7337069B2 (en) | Method of simulating the sedimentary deposition in a basin respecting the thicknesses of the sedimentary sequences | |
| EP3253946B1 (en) | Seismic attributes derived from the relative geological age property of a volume-based model | |
| Adeoti et al. | Static reservoir modeling using well log and 3-D seismic data in a KN field, offshore Niger Delta, Nigeria | |
| CN109143399B (en) | Method for identifying carbonate rock sequence interface | |
| Harris | The role of geology in reservoir simulation studies | |
| CN105386756A (en) | Method for applying dependent variables to calculate porosity of brittle formation | |
| CN106443772B (en) | One kind going diapir original stratigraphic thickness restoration methods | |
| Trippetta et al. | Carbonate-ramp reservoirs modelling best solutions: Insights from a dense shallow well database in Central Italy | |
| Dawuda et al. | Geologic modeling and ensemble-based history matching for evaluating co2 sequestration potential in point bar reservoirs | |
| Alhusseini et al. | Estimation of initial oil in place for Buzurgan Oil Field by using volumetric method and reservoir simulation method | |
| Whitaker et al. | Simulating flow through forward sediment model stratigraphies: insights into climatic control of reservoir quality in isolated carbonate platforms | |
| Yong et al. | Application of stratigraphic-sedimentological forward modeling of sedimentary processes to predict high-quality reservoirs within tight sandstone | |
| Reda et al. | Hydrocarbon reservoir characterization in the challenging structural setting of southern gulf of Suez: synergistic approach of well log analyses and 2D seismic data interpretation | |
| CUI et al. | A volumetric model for evaluating tight sandstone gas reserves in the Permian Sulige gas field, Ordos Basin, Central China | |
| Sunmonu et al. | Reservoir characterization and by-passed pay analysis of Philus Field in Niger Delta, Nigeria | |
| CN113885096B (en) | High-frequency layer sequence division and small layer comparison method and device, electronic equipment and medium | |
| CN114779333A (en) | Method and system for determining location of hydrocarbon reservoir | |
| Peijs-van Hilten et al. | Heterogeneity modeling and geopseudo upscaling applied to waterflood performance prediction of an incised valley reservoir: Countess YY Pool, southern Alberta, Canada | |
| CN111708100B (en) | Deep thin layer oil and gas reservoir determination method and related device | |
| Airen et al. | Reservoir Hydrocarbon Volumetric Analysis of Sapele Deep Field, Niger Delta, Southern Nigeria | |
| Adaeze et al. | Facies Modelling and Petrophysical Properties of X-Field, Onshore, Niger Delta, Nigeria | |
| CN117849874B (en) | Dolomite pore type thin reservoir prediction method based on reservoir cause geological model | |
| Shahbazi | Hybrid Strategy for Porosity Distribution Mapping on a Heterogeneous Reservoir Using Artificial Intelligence Method | |
| El Saadi | A geo-cellular model of a multiphase hydrocarbon field incorporating fault-seal analysis: Bredasdorp Basin, South Africa | |
| Franklin et al. | Development of an infill well program to maximise economic return from the Iagifu-Hedinia Field: part I. Integrated structural, stratigraphic, and reservoir attribute modelling as input to reservoir simulation and well targeting |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20191230 Address after: 100728 Beijing, Chaoyangmen, North Street, No. 22, No. Co-patentee after: Exploration and Development Research Institute of SINOPEC Jianghan Oilfield Branch Company Patentee after: China Petrochemical Co., Ltd. Address before: 430223, No. 18, Garden Road, East Lake hi tech Zone, Hubei, Wuhan Patentee before: Exploration and Development Research Institute of SINOPEC Jianghan Oilfield Branch Company |
|
| TR01 | Transfer of patent right |