CN107870374B - A kind of strata denudation evolution recovery method - Google Patents
A kind of strata denudation evolution recovery method Download PDFInfo
- Publication number
- CN107870374B CN107870374B CN201610851492.4A CN201610851492A CN107870374B CN 107870374 B CN107870374 B CN 107870374B CN 201610851492 A CN201610851492 A CN 201610851492A CN 107870374 B CN107870374 B CN 107870374B
- Authority
- CN
- China
- Prior art keywords
- denudation
- strata
- stratum
- plane
- recovery method
- 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 42
- 238000011084 recovery Methods 0.000 title claims abstract description 24
- 230000003628 erosive effect Effects 0.000 claims abstract description 12
- 230000000593 degrading effect Effects 0.000 claims description 10
- 230000004992 fission Effects 0.000 claims description 3
- 229910052586 apatite Inorganic materials 0.000 claims description 2
- 238000005259 measurement Methods 0.000 claims description 2
- VSIIXMUUUJUKCM-UHFFFAOYSA-D pentacalcium;fluoride;triphosphate Chemical compound [F-].[Ca+2].[Ca+2].[Ca+2].[Ca+2].[Ca+2].[O-]P([O-])([O-])=O.[O-]P([O-])([O-])=O.[O-]P([O-])([O-])=O VSIIXMUUUJUKCM-UHFFFAOYSA-D 0.000 claims description 2
- 230000001186 cumulative effect Effects 0.000 claims 1
- 230000008569 process Effects 0.000 abstract description 4
- 238000011160 research Methods 0.000 abstract description 4
- 238000004088 simulation Methods 0.000 abstract description 4
- 238000012545 processing Methods 0.000 abstract description 2
- 239000003208 petroleum Substances 0.000 abstract 1
- 239000010410 layer Substances 0.000 description 16
- 230000015572 biosynthetic process Effects 0.000 description 8
- 238000010586 diagram Methods 0.000 description 3
- 230000010429 evolutionary process Effects 0.000 description 3
- 235000019738 Limestone Nutrition 0.000 description 1
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 239000006028 limestone Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052698 phosphorus Inorganic materials 0.000 description 1
- 239000011574 phosphorus Substances 0.000 description 1
- 238000004540 process dynamic Methods 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V9/00—Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/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 kind of strata denudation evolution recovery methods.This method is to restore a kind of method of processing for stratum uplift and erosion process in petroleum geology industry, time started, end time, each Eroded Thickness data and a small amount of erosional rate model are degraded including be easier to obtain in one phase denudation of application, realize the prediction to the accumulative denudation of any time in the period is degraded, and each stratum paleostructure that underlies at any crucial moment in the phase is degraded according to accumulative erosion restoration, so that the temporal-spatial evolution analysis for stratum during geological work provides relatively reliable foundation.It is more nearly geology reality by the research purpose layer paleostructure form that this method obtains, while this method can also be provided rapidly for individual well simulation and degrade allocation plan.
Description
Technical field
The present invention relates to technical field of geological exploration, degrade more particularly, to stratum uplift during geological research dynamic
State evolution recovery method.
Background technique
Sedimentary formation uplift and erosion is a kind of most common geological evolution, carries out sedimentary formation erosion restoration and stripping
Process dynamics evolution research is lost, quickly, accurately restores in denudation, target zone paleostructure form is covered under any crucial moment,
It is of great significance to the analysis of sedimentary basin oil gas Dynamic accumulation.Although forefathers have done a large amount of work for the recovery of denudation
Make, also achieves achievement abundant, phenomenon is degraded to stratum uplift, preliminary discussion has also been made, but to denudation
Dynamic Evolution restores research and not yet carries out.Developing at present for strata denudation, there are mainly two types of methods for processing: a kind of method
Assuming that strata denudation rate is identical in denudation, i.e., the stratum of same thickness is peelled off in the identical time；Another method is false
If arbitrarily being degraded stratum, degrade that beginning, end time are all the same, i.e., in multilayer formation denudation, superstratum is degraded
Start lower layer's strata denudation after complete.In geological evolution, inside large range of sedimentary basin, it is related to multilayer formation and degrades
When, each region erosional rate is influenced by uplift size and process, in the plane and degrade shown as in evolutionary process it is non-thread
Property feature.Therefore, existing method is mainly insufficient are as follows: excessively idealizes, has in the recovery of multilayer formation denudation dynamic evolution
Body shows as 1) assuming in denudation, strata denudation rate is identical, 2) assume it is any degraded stratum, degrade beginning, terminate
Time is all the same.They with geology real process exist is not inconsistent, can not more it is careful it is accurate announcement degrade evolutionary process.
Summary of the invention
For the deficiency of the above-mentioned restoration methods of denudation in the prior art, the invention proposes a kind of new strata denudations
Evolution recovery method, comprising the following steps:
S10 obtains each layer Eroded Thickness involved in this denudation, adds up to each layer Eroded Thickness,
Obtain total denudation of this denudation；
S20, at the beginning of determining this denudation and the end time；
S30 obtains stratum any time denudation model according to stratum erosional rate model；
S40 carries out interpolation by the data to denudation models multiple in plane to obtain accumulative denudation body；
S50, at the beginning of determining each layer strata denudation according to accumulative denudation body and the end time；
S60 obtains the erosion stratum at the moment according to accumulative denudation body for any time in this denudation
Overall thickness is subtracted the stratum of the erosion overall thickness at the moment using each layer stratigraphic structure buried depth before degrading, obtains underliing for the moment
Each stratum paleostructure.
According to an embodiment of the invention, each layer stratum can be shelled using inverse distance weight in above-mentioned steps S10
Erosion amount is converted to unified grid data, then adds up, and obtains total denudation of this denudation.
According to an embodiment of the invention, the inverse distance weight assumes that n scatterplot in plane in above-mentioned steps S10
Data, each point attribute value are vi, carry out interpolation to these scatterplot data by following formula, obtain certain point attribute V in plane, wherein ri
It is the point at a distance from i-th of scatterplot,
According to an embodiment of the invention, inverse distance weight can be used this denudation in above-mentioned steps S20
Starting and end time is converted to unified grid data.
According to an embodiment of the invention, the strata denudation Rate Models can be by measuring phosphorus in above-mentioned steps S30
Certain point strata denudation rate curve in the plane that lime stone fission track data determine.
According to an embodiment of the invention, the strata denudation Rate Models can be by uplift rate in above-mentioned steps S30
Certain point strata denudation rate curve in the plane that Curve transform obtains.
According to an embodiment of the invention, in above-mentioned steps S40, it can be by using inverse distance weight to multiple in plane
The data of denudation model carry out interpolation to obtain accumulative denudation body.
According to an embodiment of the invention, the inverse distance weight is by following formula to n in plane in above-mentioned steps S40
A denudation model e1(t) ..., en(t) data carry out interpolation, obtain certain point in plane accumulative denudation body E (x, y,
T), wherein ri is the point at a distance from i-th of scatterplot,
According to an embodiment of the invention, above-mentioned accumulative denudation body E (x, y, t) is the monotonically increasing function about the time,
At the beginning of determining each layer strata denudation by linear search method and the end time.
According to an embodiment of the invention, above-mentioned strata denudation evolution recovery method, can with the following steps are included:
S70 exports each stratum paleotectonic map of any time in this denudation.
Compared with prior art, one or more embodiments of the invention can have following advantage:
Strata denudation evolution recovery method provided by the invention is that one kind can embody in strata denudation recovery evolutionary process,
The uplift for being more nearly geology real process degrades method.Using this method, can accurately and efficiently restore in the phase of degrading
Each stratum paleostructure that underlies at any crucial moment is during geological work, and the temporal-spatial evolution analysis on purpose stratum provides more
Add reliable foundation.Specifically,
(1) present invention being capable of degrading beginning, the end time, respectively degraded strata denudation by one phase denudation of input
Amount and part erosional rate model realize a phase multilayer differential denudation amount recovery process by numerical simulation；
(2) practical mistake of the phase multilayer differential denudation amount recovery process realized through the invention closer to Geological Evolution
Journey solves the problems, such as that simulation is more reasonable；
(3) the accumulative body-measure data that degrades that the present invention obtains can provide multilayer denudation for any single-point burial history reconstruction
Allocation plan can greatly improve individual well simulation precision.
(4) present invention may be equally applicable to the erosion restoration work that single layer degrades event.
Other features and advantages of the present invention will be illustrated in the following description, also, partly becomes from specification
It obtains it is clear that understand through the implementation of the invention.The objectives and other advantages of the invention can be by specification, right
Specifically noted structure is achieved and obtained in claim and attached drawing.
Detailed description of the invention
Attached drawing is used to provide further understanding of the present invention, and constitutes part of specification, with reality of the invention
It applies example and is used together to explain the present invention, be not construed as limiting the invention.In the accompanying drawings:
The flow chart of strata denudation evolution recovery method in Fig. 1 embodiment of the present invention；
Fig. 2 is the schematic diagram of strata denudation Rate Models in the embodiment of the present invention；
Fig. 3 is the stratum any time denudation obtained in the embodiment of the present invention according to the strata denudation Rate Models of Fig. 3
The schematic diagram of model；
Fig. 4 is to determine each layer strata denudation in multilayer differential denudation model according to accumulative denudation body in the embodiment of the present invention
At the beginning of and the end time schematic diagram.
Specific embodiment
Spirit of the invention be degraded using a phase be easier in event obtain always degrade the time started, at the end of
Between, each Eroded Thickness data and a small amount of erosional rate model, degrade dynamic evolution restoration methods using stratum uplift,
It realizes the prediction for adding up denudation to any time in the phase of degrading, and denudation is added up according to this, restore to throw oneself on the ground under any time
Layer paleostructure.
To make the object, technical solutions and advantages of the present invention clearer, below in conjunction with drawings and examples to the present invention
It is described in further detail.
As shown in Figure 1, strata denudation evolution recovery method mainly comprises the steps that in embodiment
S10 obtains each layer Eroded Thickness involved in this denudation, adds up to each layer Eroded Thickness,
Obtain total denudation of this denudation.
Wherein, input in a denudation involved each layer Eroded Thickness data, and preferably using instead away from
Each layer is degraded into data conversion into uniform grid data from weighting method, then adds up to each layer denudation, is degraded with obtaining this
Total denudation of event.
Specifically, inverse distance weight assumes that n scatterplot data in plane, and each point attribute value is vi, passes through following formula pair
These scatterplot data carry out interpolation, obtain certain point attribute V in plane, and wherein ri is the point at a distance from i-th of scatterplot,
S20, at the beginning of determining this denudation and the end time.
Similarly, it preferably will be converted into uniting with the end time at the beginning of this denudation using inverse distance weight
One grid data.
S30 obtains stratum any time denudation model according to stratum erosional rate model.
In the present embodiment, the erosional rate model is that certain point is degrading stage any time strata denudation speed in plane
Rate curve (as shown in Figure 2).The related data of rate curve can be obtained by measurement Apatite fission track data, can also be with
It is obtained by application uplift rate Curve transform.As shown in figure 3, can be by being integrated to stratum erosional rate curve, it will
It is converted to any time Eroded Thickness model.Such as following formula:
S40 carries out interpolation by the data to denudation models multiple in plane to obtain accumulative denudation body.
Accumulative denudations multiple in plane are obtained into accumulative denudation body by inverse distance weight interpolation.
If n denudation model is e in plane1(t) ..., en(t), it is carried out by the data of these denudation models slotting
It is accumulative denudation body Model E (x, y, t) to obtain certain point in plane for value.Likewise it is possible to be added using the anti-distance being shown below
Power method carries out interpolation to the data of denudation models multiple in plane, obtains the accumulative denudation body Model of certain point in plane.
In formula, ri is that certain is put at a distance from i-th of scatterplot in plane.
S50, at the beginning of determining each layer strata denudation according to accumulative denudation body and the end time.
Fig. 4 is a multilayer differential denudation model in the embodiment of the present invention.Wherein L1, L2, L3, L4, L5 indicate stratal surface；
H1, h2, h3, h4 respectively indicate L1-L2, L2-L3, L3-L4, L4-L5 formation thickness；E1, e2, e3, e4 indicate stratum primary
It degrades in event 4 and degrades time point (or moment).For the principle that the present invention is further explained, can select without loss of generality
Select any (x0, y0) point, it is assumed that once degrading event start time is e4, end time e1, the e3 moment degrade just to
L3 level (stratum L3-L4, L4-L5 of face covering has been degraded, and the stratum L2-L3 not yet starts to degrade), it is lucky at the e2 moment
It degrades to L2 level.In actual operation the known primary event of degrading start, the end time, i.e., known e4, e1, it is known that degraded
Each formation thickness that event influences, i.e. L1-L2, L2-L3, L3-L4, L4-L5 formation thickness is respectively h1, h2, h3, h4, then needs
It seeks being degraded degrading beginning and end time for each stratum of event influence, that such as seeks the stratum L2-L3 degrades beginning, knot
Beam time e3 and e2.And the two material times seek passing through respectively equation E (x0, y0, t)=h4+h3 and E (x0, y0,
T)=h4+h3+h2 is solved.In the present embodiment, from the figure 3, it may be seen that E (x, y, t) is the monotonic increasing function about the time, therefore can
To degrade starting and end time (as shown in Figure 4) using linear search method acquisition L2-L3 level.
S60 obtains the erosion stratum at the moment according to accumulative denudation body for any time in this denudation
Overall thickness is subtracted the stratum of the erosion overall thickness at the moment using each layer stratigraphic structure buried depth before degrading, obtains underliing for the moment
Each stratum paleostructure.
S70 exports any time underlying strata paleotectonic map in this denudation.
The above, specific implementation case only of the invention, scope of protection of the present invention is not limited thereto, any ripe
Those skilled in the art are known in technical specification of the present invention, modifications of the present invention or replacement all should be in the present invention
Protection scope within.For example, inverse distance-weighting described in the present invention is a kind of method of mathematical interpolation, not uniquely,
Plane and body interpolation algorithm can also be using interpolation algorithms such as Ke Lijin, sequence Gauss.
Claims (10)
1. a kind of strata denudation evolution recovery method, comprising the following steps:
S10 obtains each layer Eroded Thickness involved in this denudation, adds up to each layer Eroded Thickness, obtains
Total denudation of this denudation；
S20, at the beginning of determining this denudation and the end time；
S30 obtains stratum any time denudation model according to stratum erosional rate model；
S40 carries out interpolation by the data to denudation models multiple in plane to obtain accumulative denudation body；
S50, at the beginning of determining each layer strata denudation according to accumulative denudation body and the end time；
S60 obtains the stratum of the erosion total thickness at the moment according to accumulative denudation body for any time in this denudation
Degree is subtracted the stratum of the erosion overall thickness at the moment using each layer stratigraphic structure buried depth before degrading, obtains the various regions of underliing at the moment
Layer paleostructure.
2. strata denudation evolution recovery method according to claim 1, it is characterised in that:
In step S10, each layer Eroded Thickness is converted to by unified grid data using inverse distance weight, then into
Row is cumulative, obtains total denudation of this denudation.
3. strata denudation evolution recovery method according to claim 2, it is characterised in that:
In step S10, the inverse distance weight assumes that n scatterplot data in plane, and each point attribute value is vi, passes through following formula
Interpolation is carried out to these scatterplot data, obtains certain point attribute V in plane, wherein ri is the point at a distance from i-th of scatterplot,
4. strata denudation evolution recovery method according to claim 1, it is characterised in that:
In step S20, unification will be converted to the end time using inverse distance weight at the beginning of this denudation
Grid data.
5. strata denudation evolution recovery method according to claim 1, it is characterised in that:
In step S30, the strata denudation Rate Models are certain in the plane determined by measurement Apatite fission track data
Point strata denudation rate curve.
6. strata denudation evolution recovery method according to claim 1, it is characterised in that:
In step S30, the strata denudation Rate Models are certain point stratum strippings in the plane obtained by uplift rate Curve transform
Lose rate curve.
7. strata denudation evolution recovery method according to claim 1, it is characterised in that:
In step S40, interpolation is carried out by using data of the inverse distance weight to denudation models multiple in plane to obtain
Accumulative denudation body.
8. strata denudation evolution recovery method according to claim 7, it is characterised in that:
In step S40, the inverse distance weight is by following formula to n denudation model e in plane1(t) ..., en(t)
Data carry out interpolation, obtain the accumulative denudation body E (x, y, t) of certain point in plane, wherein ri be the point and i-th scatterplot away from
From,
Wherein eiIt (t) is i-th of denudation model in 1 to n denudation model in plane.
9. strata denudation evolution recovery method according to claim 8, it is characterised in that:
The accumulative denudation body E (x, y, t) is the monotonically increasing function about the time, with determining each layer by linear search method
Layer degrade at the beginning of and the end time.
10. strata denudation evolution recovery method according to claim 8, which is characterized in that further comprising the steps of:
S70 exports each stratum paleotectonic map that underlies of any time in this denudation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610851492.4A CN107870374B (en) | 2016-09-26 | 2016-09-26 | A kind of strata denudation evolution recovery method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610851492.4A CN107870374B (en) | 2016-09-26 | 2016-09-26 | A kind of strata denudation evolution recovery method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107870374A CN107870374A (en) | 2018-04-03 |
| CN107870374B true CN107870374B (en) | 2019-06-18 |
Family
ID=61751881
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610851492.4A Active CN107870374B (en) | 2016-09-26 | 2016-09-26 | A kind of strata denudation evolution recovery method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107870374B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108828685B (en) * | 2018-04-24 | 2019-05-10 | 北京师范大学 | A kind of determination method of formation erosion thickness |
| CN111624674A (en) * | 2020-06-09 | 2020-09-04 | 中国地质大学（北京） | Simulation system for evolution of deposition sequence formation |
| CN119045051B (en) * | 2023-05-29 | 2025-09-12 | 中国石油天然气股份有限公司 | A method for correcting stratum erosion using two-dimensional simulation |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104122602A (en) * | 2014-07-31 | 2014-10-29 | 中国石油大学(北京) | Method and device for acquiring erosion thickness of sedimentary basin |
| CN105044797A (en) * | 2015-08-06 | 2015-11-11 | 中国石油天然气股份有限公司 | Quantitative recovery method for carbonate rock stratum denudation amount |
| CN105204069A (en) * | 2014-06-27 | 2015-12-30 | 中国石油化工股份有限公司 | Strata denudation amount recovery method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2384304B (en) * | 2002-01-04 | 2003-12-03 | Nigel Allister Anstey | Method of distinguishing types of geologic sedimentation |
-
2016
- 2016-09-26 CN CN201610851492.4A patent/CN107870374B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105204069A (en) * | 2014-06-27 | 2015-12-30 | 中国石油化工股份有限公司 | Strata denudation amount recovery method |
| CN104122602A (en) * | 2014-07-31 | 2014-10-29 | 中国石油大学(北京) | Method and device for acquiring erosion thickness of sedimentary basin |
| CN105044797A (en) * | 2015-08-06 | 2015-11-11 | 中国石油天然气股份有限公司 | Quantitative recovery method for carbonate rock stratum denudation amount |
Non-Patent Citations (2)
| Title |
|---|
| Luca M. Abbühl et al..El Niño forcing on 10Be-based surface denudation rates in the northwestern Peruvian Andes?.《Geomorphology》.2010,第123卷第257-268页. |
| 江陵凹陷荆沙组地层剥蚀厚度恢复探讨;陈晓辉等;《长江大学学报（自然科学版）》;20101231;第7卷(第4期);第48-52页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107870374A (en) | 2018-04-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104765065B (en) | A method and device for determining paleogeomorphology of sedimentary strata | |
| CA2912626C (en) | 3d trap evaluation method of searching for oil-gas reservoir | |
| CN102748016A (en) | Geologic history period sandstone reservoir porosity evolution recovery method | |
| CN103077558A (en) | Method for establishing large-sized karst cave reservoir body distribution model of fractured-vuggy carbonate rock reservoir | |
| CN110320574B (en) | Method for describing thin-layer sand body based on gentle slope delta | |
| CN105425315B (en) | A kind of small scale depositional phase ancient landform inversion method in weak structural deformation area | |
| CN103293554B (en) | Statics correction method based on seismic data | |
| CN107870374B (en) | A kind of strata denudation evolution recovery method | |
| CN111308549B (en) | Variable-speed mapping method based on model inversion | |
| CN105607119B (en) | Near-surface model construction method and static correction value acquiring method | |
| Mickus et al. | Gravity analysis of the main Ethiopian rift | |
| CN102879820A (en) | Three-dimensional surface layer model construction method based on triangular grid | |
| CN115373024B (en) | Method and device for inverting passive land edge crustal structure based on stratum recording settlement | |
| CN104375172A (en) | Volcanic underlayer structural configuration correct imaging method | |
| CN111624651B (en) | Reservoir prediction method and device based on paleo-topography constraint | |
| Zhao et al. | Velocity structure heterogeneity and tectonic motion in and around the Tan–Lu fault of China | |
| CN106772584B (en) | A kind of fine volume elements attribute assignment method for taking geological tectonic conditions constraint into account | |
| CN103628866A (en) | Method for obtaining stratum pressure coefficients of stratum | |
| CN103837140B (en) | A Surveying and Mapping Method for Initial Back-arc Basin Continental Shelf Delineation | |
| CN104898181B (en) | The theoretical maximum wave height and cumulative frequency wave height of ancient times phase determines method | |
| CN109188525B (en) | A kind of acquisition methods and system of marine facies rammell buried depth data | |
| CN106908838A (en) | The method for building target area stratum inelastic attenuation quality factor three-dimensional model | |
| CN110488385A (en) | Method for making three-dimensional fractured-karst model | |
| CN108254779A (en) | Three-dimensional construction modeling method | |
| CN105205863B (en) | A kind of 3-dimensional digital basin modeling method for grading |
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 |