CN107991714A - The quantitative approach recovered based on lake basin paleotopography - Google Patents
The quantitative approach recovered based on lake basin paleotopography Download PDFInfo
- Publication number
- CN107991714A CN107991714A CN201711213606.3A CN201711213606A CN107991714A CN 107991714 A CN107991714 A CN 107991714A CN 201711213606 A CN201711213606 A CN 201711213606A CN 107991714 A CN107991714 A CN 107991714A
- Authority
- CN
- China
- Prior art keywords
- thickness
- depth
- formation
- stratum
- lake basin
- 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.)
- Granted
Links
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
A kind of quantitative approach recovered based on lake basin paleotopography, using following steps:One:Establish depressed area high accuracy chronostratigraphic architecture;Two:Selection marker layer, determines two sections of depressed area residual stratum thickness;Three:It is actual formation thickness by the residual stratum thickness correction of target zone using mud stone sonic time difference;Four:Carry out stratigraphic dip correction;Five:Compaction treatment is carried out to the actual formation thickness of retrieving layer position, to determine embryo deposit period thickness;Six:By core analysis and thickness measure, recover paleao-water depth value at several well points, and according to paleao-water depth embryo deposit period thickness correspondence, formation thickness is converted into paleao-water depth distribution, i.e. quantitative reconstruction lake basin landform.The present invention not only solves the problem of stratigraphic dip, strata denudation and formation compaction effect are to thickness effect;Moreover, going back quantitatively characterizing lake basin landform, quantify paleao-water depth distribution, solved precision problem when recovering for lake basin landform, improve the recovery precision of lake basin landform.
Description
Technical field
The present invention relates to lake basin paleotopography, more particularly to a kind of quantitative approach recovered based on lake basin paleotopography.Belong to oil
Gas field exploration development technique field.
Background technology
The correlative study that lake basin landform is recovered originates from the fifties in last century, by the development in more than 60 years, the research
It is widely used in multiple subjects.In geological research, lake basin landform is recovered to play the effect being even more important, to not
Reservoir genesis, development models and the Distribution Pattern of synsedimentary type provide important guidance and reference function.It is in addition, raw in Gu
Very important effect is also played in thing, ancient geoaraply, paleotectonic research.However, at present, lake basin landform restoration methods exist
Problems, rest on qualitative stage mostly, therefore, cause quantitative study not carry out, this is also that lake basin landform recovers needs
The direction of development.
The seventies in last century, lake basin landform gradually play a significant role during oil-gas exploration and development.At present, often
Lake basin landform restoration methods include:Flattening method, sedimentary analysis method, impression method etc., difference of these methods in oil field
Development phase effectively directs oil-gas exploration and exploitation.Wherein, flattening method is mainly in combination with Based On Sequence Stratigraphic Theory and thing
Visit theoretical method, basic principle be using maximum flood surface or base level of deposition as object of reference, sequence calibration on the basis of, profit
Bottom surface subtracts top surface when being travelled with earthquake, and using the top surface after evening up as deposition lake level, corresponding bottom surface form, is the layer
Lake basin landform before sequence sedimentation;Sedimentary analysis method is that depositional environment and deposit are analyzed, according to deposition object space
Configuration relation and lithology combination feature etc. recover lake basin landform;Stamp rule be make use of participate in lake basin landform and superstratum it
Between " mirror ", pass through superstratum recover deposition period lake basin landform.
Three of the above method more or less be there are problems that in practical applications.Such as:Flattening method is relatively directly perceived true,
And it is easy to operate, but error is larger;Sedimentary analysis method needs to consider many factors, and workload is huge;Although impression method
Principle is simple but precision is not good enough.Moreover, the above method does not account for large dip angle strata recovers the important of generation to lake basin landform
Influence, in addition, lacking consideration to the influence factor aspect of paleao-water depth.
The content of the invention
It is a primary object of the present invention to overcome disadvantages mentioned above existing in the prior art, and provide a kind of based on the ancient ground of lake basin
The quantitative approach that shape is recovered, can not only be right that takes into account the influence of strata denudation, stratigraphic dip and differential compaction to thickness
Lake basin landform is recovered, and solves the problems, such as stratigraphic dip, strata denudation and formation compaction effect to thickness effect;Moreover,
Residual stratum thickness correction is converted to depth of water during primary deposit, i.e. quantitatively characterizing lake basin landform, has quantified paleao-water depth distribution,
Solve precision problem when recovering for lake basin landform, substantially increase the recovery precision of lake basin landform;Meanwhile realize
Transformation of the lake basin landform from qualitative to quantitative, makes restoration result more accurately and reliably.
The purpose of the present invention is what is realized by following technical scheme:
1st, a kind of quantitative approach recovered based on lake basin paleotopography, it is characterised in that:Using following steps:
The first step:Establish depressed area high accuracy chronostratigraphic architecture;
Second step:Selection marker layer, determines two sections of depressed area residual stratum thickness;
3rd step:It is actual formation thickness by the residual stratum thickness correction of target zone using mud stone sonic time difference;
4th step:Carry out stratigraphic dip correction;
5th step:Compaction treatment is carried out to the actual formation thickness of retrieving layer position, to determine embryo deposit period thickness;
6th step:By core analysis and thickness measure, recover paleao-water depth value at several well points, and according to paleao-water depth-just
Beginning deposition period thickness correspondence, paleao-water depth distribution, i.e. quantitative reconstruction lake basin landform are converted to by formation thickness.
The specific practice of the first step:It is a set of on identifying purpose stratum according to Theory of High-Resolution Sequence Stratigraphy
Stable mud stone section is maximum flood surface, determines sequence interface, and establishes the high-precision chronostratigraphic architecture of target zone, by maximum
Flood surface, determines sequence interface and is divided into two complete third-band sequence Q1, Q2.
The specific practice of the second step:Selection it is nearest, most stable of Q1 flood surfaces be reference lamina, and in seismic data
Layer position corresponding to the upper maximum flood surface Q1 in explanation two sections of depressed area residual stratum bottom surface and its inside, and it is true according to the following equation
Determine the formation thickness value between two sections of depressed area residual stratum bottom surface and reference lamina, i.e.,:Residual stratum thickness:
T3=T1-T2
Wherein, T3 represents residual stratum thickness, and T1 represents the depth of reference lamina on target zone, and T2 represents to recover destination
The depth of layer bottom interface.
The specific practice of 3rd step:The mud stone sound of the several wells in Yellow River mouth sag area is counted by well-log information
The ripple time difference and corresponding depth, fitting interval transit time and depth relationship Δ M=Δ Moe (- CH)
In formula:Δ Mo is the interval transit time value that earth's surface does not consolidate mud stone, and unit is:us/m;
C values are to represent normal compaction slope of a curve;
Δ M is:The mud stone interval transit time of any buried depth, unit are:us/m;H is:Mud stone buried depth, unit m;E is
Natural logrithm bottom;
(1) it will be extended down to outside the compaction curve of the following mud stone of plane of unconformity at Δ M=Δs Mo, be:Embryo deposit earth's surface, just
The deposition that begins the distance between earth's surface and plane of unconformity are erosion sediment thickness T4;
(2) it is target zone residual stratum thickness correction is as follows for the calculation formula of actual formation thickness:
T5=T3+T4
In formula;T5 is:Actual formation thickness;T3 represents residual stratum thickness;T4 is:Degrade formation thickness.
The specific practice of 4th step:
(1), according to selected depressed area top surface structural map, the distribution of selected depressed area plane stratigraphic dip is calculated;
(2), using geometric function method correction stratigraphic dip, dip correction is carried out to actual formation thickness;Assuming that T6 is stratum
Actual thickness, T5 is that drilled well bore bores the formation thickness met, α angles between T5 and T6, then T6=T5*cos α
In formula:α can be calculated on selected depressed area top surface structural map.
The specific practice of 5th step:
1. buried depth (H)~porosity (Φ) equation of depressed area depression sandstone-mud stone selected by fitting;From well-log information
The porosity value of the different depth of statistical number mouth well, sand, depth (H)~porosity of mud stone by exponential function fitting of distribution
(Φ) equation:
Ф=Ф0exp(-c*h) (1)
Wherein:H is the buried depth depth on stratum, m;The rock porosity that it is h in depth that Φ, which is,;Φ0It in depth is h=0 to be
Porosity;C is compacting factor, m-1;
The sand of fitting, the porosity of mud stone~depth relationship equation are as follows:
Mud stone:Φ=0.5947exp (- 0.000762*h) (2)
Sandstone:Φ=0.4906exp (- 0.000442*h) (3)
2. the sand shale thickness of statistical number mouth well respectively, and according to skeleton thickness iterative Integral Method mathematical model, pass through ground
Layer decompaction correction technology, is corrected to actual formation thickness (T6) formation thickness (T7) during embryo deposit;
The thickness of the sand shale of depressed area selected by statistics respectively, will be real using skeleton thickness iterative Integral Method mathematical model
Border formation thickness T6 is corrected to embryo deposit period formation thickness T7;The Computing Principle of skeleton thickness iterative Integral Method mathematical model
It is as follows, depth of stratum is calculated first as the weighting porosity at Z:
Ф (Z)=Ps Ф s (Z)+Pm Ф m (Z) (4)
In formula;Ps and Pm is the sandstone on stratum and the content of mud stone respectively, decimal;Ф s (Z) and Ф m (Z) are depth respectively
For the sandstone and the porosity of mud stone at Z;
Secondly, using skeleton thickness iterative Integral Method mathematical model, primary deposit period will be reverted to by formation thickness now
Formation thickness:
In formula:T7 is embryo deposit period formation thickness, and Ф (Z) is that depth is weighting porosity at Z, Z1 and Z2 difference
For the top circle on stratum and the depth on bottom circle;
Bring formula (2)-(4) into formula (5) and integrate, utilize above formula, you can obtain embryo deposit period formation thickness
T7;Then, further according to well log interpretation data, the percentage composition of depressed area sandstone selected by difference, then institute is recovered using the above method
The actual formation thickness of depressed area number mouthful well is selected, and counts formation compaction rate (compacting rate=embryo deposit formation thickness T7/ reality
Formation thickness corrects T6), the formation compaction rate of selected depressed area is generally in 1.8-2.0 scopes.
The skeleton thickness iterative Integral Method mathematical model has following premise when establishing:In whole deposition process, stratum
Skeleton volume do not change;And formation compaction process does not reverse.
The specific practice of 6th step:
1. by core analysis lithofacies characteristics, observation grain graininess, rock deposition construction and ground mark etc. come qualitative
Judge the scope of the depth of water;
2. being directed to lake depositional environment, " beach dam Thickness Method " Fine calibration depth of water is proposed;In the hydrodynamic force band of lake
The beach dam of a variety of origin causes of formation has been developed, has corresponded respectively to littoral dam, offshore dam and remote bank dam;First formed according to lithofacies characteristics and beach dam
Its origin cause of formation of position judgment, then, can be similar to the depth of water herein further according to the deposit thickness on beach dam of single phase time;
3. by the thickness after several Kou Jing beaches dam thickness correction, i.e.,:For paleao-water depth T8 herein;Several mouthfuls of well point depth of water T8 are counted
With the relation between initial formation thickness T7, and according to relevance formula, thickness conversion during by selected depressed area embryo deposit
For the quantitative distribution of the depth of water, it is:Quantitative reconstruction lake basin landform.
Beneficial effects of the present invention:Due to the adoption of the above technical scheme, that takes into account strata denudation, stratigraphic dip by the present invention
And influence of the differential compaction to thickness, lake basin landform can not only be recovered, solve stratigraphic dip, strata denudation and
The problem of formation compaction effect is to thickness effect;Moreover, depth of water when residual stratum thickness correction to be converted to primary deposit, i.e.,
Quantitatively characterizing lake basin landform, has quantified paleao-water depth distribution, solves precision problem when recovering for lake basin landform, significantly
Improve the recovery precision of lake basin landform;Meanwhile transformation of the lake basin landform from qualitative to quantitative is realized, make restoration result more
Accurately and reliably.
Brief description of the drawings
Fig. 1 is flow diagram of the present invention.
Fig. 2 is chronostratigraphic architecture schematic diagram of the present invention.
Fig. 3 corrects schematic diagram for stratigraphic dip of the present invention.
Fig. 4 is that sand shale porosity-depth relationship formula schematic diagram is fitted in decompaction correction of the present invention.
Fig. 5 recovers " sandstone thickness method " schematic diagram for paleao-water depth of the present invention.
Embodiment
As Figure 1-Figure 5, the present invention uses following steps:
The first step:Establish depressed area high accuracy chronostratigraphic architecture;
According to Theory of High-Resolution Sequence Stratigraphy, the mud stone section of a set of stabilization is general for most big lake on identifying purpose stratum
Face, determines sequence interface, and establishes the high-precision chronostratigraphic architecture of target zone, by maximum flood surface, determines sequence interface and draws
It is divided into two complete third-band sequence Q1, Q2;
Second step:Selection marker layer, determines two sections of depressed area residual stratum thickness;
Selection it is nearest, most stable of Q1 flood surfaces be reference lamina, and explanation two sections of depressed area remnants on seismic data
Stratum bottom surface and its layer position corresponding to internal maximum flood surface Q1, and two sections of depressed area residual stratum is determined according to the following equation
Formation thickness value between bottom surface and reference lamina, i.e.,:Residual stratum thickness:
T3=T1-T2
Wherein, T3 represents residual stratum thickness, and T1 represents the depth of reference lamina on target zone, and T2 represents to recover destination
The depth of layer bottom interface.
3rd step:It is actual formation thickness by the residual stratum thickness correction of target zone using mud stone sonic time difference;
By well-log information count the several wells in Yellow River mouth sag area (20 mouthfuls of wells of the present embodiment) mud stone interval transit time and
Corresponding depth, is fitted interval transit time and depth relationship
Δ M=Δ Moe (- CH)
In formula:Δ Mo is the interval transit time value that earth's surface does not consolidate mud stone, and unit is:us/m;
C values are to represent normal compaction slope of a curve;
Δ M is:The mud stone interval transit time of any buried depth, unit are:us/m;H is:Mud stone buried depth, unit m;E is
Natural logrithm bottom;
It will be extended down to outside the compaction curve of the following mud stone of plane of unconformity at Δ M=Δs Mo, be:Embryo deposit earth's surface, initially
It is erosion sediment thickness T4 to deposit the distance between earth's surface and plane of unconformity;
(2) it is target zone residual stratum thickness correction is as follows for the calculation formula of actual formation thickness:
T5=T3+T4
In formula;T5 is:Actual formation thickness;T3 represents residual stratum thickness;T4 is:Degrade formation thickness.
4th step:Carry out stratigraphic dip correction;
Due to the stratum observed and remained now, topography variation or tectonic activities when living through deposition
It is strong to influence, cause nowadays stratum form and well track to show as certain angle oblique.Therefore, well bore bores the formation thickness met
It is not the actual thickness on stratum, it is necessary to be corrected for the inclination effect on stratum.
1. as shown in figure 3, according to selected depressed area top surface structural map, the distribution of its plane stratigraphic dip is calculated, this
Embodiment is:Operated with professional treatment software.
2. using geometric function method correction stratigraphic dip, dip correction is carried out to actual formation thickness.Assuming that T6 is stratum
Actual thickness, T5 be drilled well bore bore meet formation thickness, α angles between T5 and T6, then
T6=T5*cos α
In formula:α can be calculated on selected depressed area top surface structural map.
5th step:Compaction treatment is carried out to the actual formation thickness of retrieving layer position, to determine embryo deposit period thickness;
Since the earth formation of different regions differs greatly, under the compacting of overlying formation pressure, different regions compacting feelings
Condition is different.Therefore, decompaction correction is one of committed step that lake basin landform is recovered.
(1) buried depth (H)~porosity (Φ) equation of depression sandstone-mud stone in depressed area selected by fitting;From well-log information
The porosity value of the different depth of statistical number mouth well (20 mouthfuls of wells of the present embodiment), sand, mud stone by exponential function fitting of distribution
Depth (H)~porosity (Φ) equation:
Ф=Ф0exp(-c*h) (1)
Wherein:H is the buried depth depth on stratum, m;The rock porosity that it is h in depth that Φ, which is,;Φ0It in depth is h=0 to be
Porosity;C is compacting factor, m-1。
As shown in figure 4, the Φ of the sedimentation type such as different lithology, different grain size0And c value differences are not larger, association index function,
The sand of fitting, the porosity of mud stone~depth relationship equation, it is as a result as follows:
Mud stone:Φ=0.5947exp (- 0.000762*h) (2)
Sandstone:Φ=0.4906exp (- 0.000442*h) (3)
(2) the sand shale thickness of 20 mouthfuls of wells is counted respectively, and according to skeleton thickness iterative Integral Method mathematical model, passes through ground
Layer decompaction correction technology, is corrected to actual formation thickness (T6) formation thickness (T7) during embryo deposit.
Skeleton thickness iterative Integral Method mathematical model has following premise when establishing:In whole deposition process, the bone on stratum
Frame body product does not change;And formation compaction process does not reverse.
The thickness of the sand shale of depressed area selected by statistics respectively, will be real using skeleton thickness iterative Integral Method mathematical model
Border formation thickness T6 is corrected to embryo deposit period formation thickness T7.The Computing Principle of skeleton thickness iterative Integral Method mathematical model
It is as follows, depth of stratum is calculated first as the weighting porosity at Z:
Ф (Z)=Ps Ф s (Z)+Pm Ф m (Z) (4)
In formula;Ps and Pm is the sandstone on stratum and the content of mud stone respectively, decimal;Ф s (Z) and Ф m (Z) are depth respectively
For the sandstone and the porosity of mud stone at Z.
Secondly, using skeleton thickness iterative Integral Method mathematical model, primary deposit period will be reverted to by formation thickness now
Formation thickness:
In formula:T7 is embryo deposit period formation thickness, and Ф (Z) is that depth is weighting porosity at Z, Z1 and Z2 difference
For the top circle on stratum and the depth on bottom circle.
Bring formula (2)-(4) into formula (5) and integrate, embryo deposit period formation thickness can be obtained using above formula
T7。
According to well log interpretation data, the percentage composition of depressed area sandstone selected by difference, then Huang has been recovered using the above method
The actual formation thickness of river mouth about 20 mouthfuls of wells of depression, and count formation compaction rate (compacting rate=embryo deposit formation thickness T7/ realities
Border formation thickness correction T6), the formation compaction rate of selected depressed area is generally in 1.8-2.0 or so.
6th step:By core analysis and thickness measure, recover paleao-water depth value at some well points, and according to paleao-water depth-just
Formation thickness is converted to paleao-water depth distribution, i.e. quantitative reconstruction lake basin landform by beginning deposition period thickness correspondence;
Formation thickness is only the qualitative description to lake basin landform, and the landform after paleao-water depth corrects is only most true, accurate
Lake basin landform.But quantitative reconstruction paleao-water depth is relatively difficult, for lake beach dam phase sedimentary characteristic, it is proposed that paleao-water depth recovers
New approaches:Paleao-water depth scope is determined using traditional lithofacies analysis method, it is quantitative further according to " the beach dam Thickness Method " for lake
Demarcate paleao-water depth, finally, according to paleao-water depth~embryo deposit when thickness (T7) correspondence, formation thickness is converted into Gu Shui
Deep quantitative distribution.
1. by core analysis lithofacies characteristics, observation grain graininess, rock deposition construction and ground mark etc. come qualitative
Judge the scope of the depth of water;Since from bank to the middle of a lake, the glutenite content of deposition is reduced, the clay content of deposition becomes more, so,
According to the scope of the particle of rock depth of water when that can reflect deposition.Such as:The depth of water is 1- when conglomerate, medium coarse sand corresponding deposition
10 meters, and the depth of water is 5-20m when packsand, siltstone corresponding deposition.Sedimentary structure of rock etc. is also influenced by the depth of water, such as
Parallel bedding, current bedding etc. are more developed in depth of water 5-20m environment.
2. as shown in figure 5, it is directed to lake depositional environment, it is proposed that " beach dam Thickness Method " Fine calibration depth of water;In lake
(surfing band, breaker zone, cleave band) has developed the beach dam of a variety of origin causes of formation in hydrodynamic force band, corresponds respectively to littoral dam, offshore
Dam and remote bank dam etc..First its origin cause of formation is judged according to lithofacies characteristics and beach dam forming position, then according to the deposition on beach dam of single phase time
Thickness can be similar to the depth of water herein (due to beach dam of single phase time from embryo deposit to deposition at the end of thickness be similar to beach dam and sink
Depth of water during product).
Exemplified by studying the 3D wells in area, pass through core observation:Beach dam deposition top bottom is directly contacted with exposed red
Mud stone (3180.5m), confirms that this section of rock core is in the surfing band of exposed environments once in a while;(channel-shaped is seen in lower part to sandstone based on cobble
Cross-bedding:3179.4m), sort it is poor, align, shale content it is higher, drop hydrochloric acid has no the reaction of obvious bubble, again
Confirm that this section is in the medium surfing band of hydrodynamic condition, according to the surfing of selected depressed area band depth of water empirical value, confirm the well
Depth of water scope at point is 1-6m.Due to a complete beach dam thickness in surfing band, it is similar to send out after formation compaction corrects
Depth of water when educating, so by the thickness (about 3m) after 3D Jing Tan dams thickness correction, herein as paleao-water depth T8.20 mouthfuls of wells are counted
Relation between point depth of water T8 and initial formation thickness T7, according to relevance formula, thickness during by selected depressed area embryo deposit
Degree is converted to the quantitative distribution of the depth of water, is:Quantitative reconstruction lake basin landform.
The above described is only a preferred embodiment of the present invention, not make limitation in any form to the present invention, it is all
It is any simple modification, equivalent change and modification made according to the technical spirit of the present invention to above example, still falls within
In the range of technical solution of the present invention.
Claims (8)
- A kind of 1. quantitative approach recovered based on lake basin paleotopography, it is characterised in that:Using following steps:The first step:Establish depressed area high accuracy chronostratigraphic architecture;Second step:Selection marker layer, determines two sections of depressed area residual stratum thickness;3rd step:It is actual formation thickness by the residual stratum thickness correction of target zone using mud stone sonic time difference;4th step:Carry out stratigraphic dip correction;5th step:Compaction treatment is carried out to the actual formation thickness of retrieving layer position, to determine embryo deposit period thickness;6th step:By core analysis and thickness measure, recover paleao-water depth value at several well points, and according to paleao-water depth-initial heavy Product period thickness correspondence, paleao-water depth distribution, i.e. quantitative reconstruction lake basin landform are converted to by formation thickness.
- 2. the quantitative approach according to claim 1 recovered based on lake basin paleotopography, it is characterised in that:The first step Specific practice:According to Theory of High-Resolution Sequence Stratigraphy, the mud stone section of a set of stabilization is most big lake on identifying purpose stratum General face, determines sequence interface, and establishes the high-precision chronostratigraphic architecture of target zone, by maximum flood surface, determines sequence interface It is divided into two complete third-band sequence Q1, Q2.
- 3. the quantitative approach according to claim 1 recovered based on lake basin paleotopography, it is characterised in that:The second step Specific practice:Selection it is nearest, most stable of Q1 flood surfaces be reference lamina, and explanation two sections of depressed area remnants on seismic data Stratum bottom surface and its layer position corresponding to internal maximum flood surface Q1, and two sections of depressed area residual stratum is determined according to the following equation Formation thickness value between bottom surface and reference lamina, i.e.,:Residual stratum thickness:T3=T1-T2Wherein, T3 represents residual stratum thickness, and T1 represents the depth of reference lamina on target zone, and T2 represents to recover purpose stratum bottom The depth at interface.
- 4. the quantitative approach according to claim 1 recovered based on lake basin paleotopography, it is characterised in that:3rd step Specific practice:The mud stone interval transit time of the several wells in Yellow River mouth sag area and corresponding depth are counted by well-log information, are fitted sound wave The time difference and depth relationship Δ M=Δ Moe (- CH)In formula:Δ Mo is the interval transit time value that earth's surface does not consolidate mud stone, and unit is:us/m;C values are to represent normal compaction slope of a curve;Δ M is:The mud stone interval transit time of any buried depth, unit are:us/m;H is:Mud stone buried depth, unit m;E is nature Logarithm bottom;(1) it will be extended down to outside the compaction curve of the following mud stone of plane of unconformity at Δ M=Δs Mo, be:Embryo deposit earth's surface, it is initial heavy Product the distance between earth's surface and plane of unconformity are erosion sediment thickness T4;(2) it is target zone residual stratum thickness correction is as follows for the calculation formula of actual formation thickness:T5=T3+T4In formula;T5 is:Actual formation thickness;T3 represents residual stratum thickness;T4 is:Degrade formation thickness.
- 5. the quantitative approach according to claim 1 recovered based on lake basin paleotopography, it is characterised in that:4th step Specific practice:(1), according to selected depressed area top surface structural map, the distribution of selected depressed area plane stratigraphic dip is calculated;(2), using geometric function method correction stratigraphic dip, dip correction is carried out to actual formation thickness;Assuming that T6 is the true of stratum Real thickness, T5 are that drilled well bore bores the formation thickness met, α angles between T5 and T6, then T6=T5*cos αIn formula:α can be calculated on selected depressed area top surface structural map.
- 6. the quantitative approach according to claim 1 recovered based on lake basin paleotopography, it is characterised in that:5th step Specific practice:1. buried depth (H)~porosity (Φ) equation of depressed area depression sandstone-mud stone selected by fitting;Counted from well-log information The porosity value of the different depth of several mouthfuls of wells, depth (H)~porosity (Φ) side of sand, mud stone by exponential function fitting of distribution Journey:Ф=Ф0exp(-c*h) (1)Wherein:H is the buried depth depth on stratum, m;The rock porosity that it is h in depth that Φ, which is,;Φ0It is in the hole that depth is h=0 Degree;C is compacting factor, m-1;The sand of fitting, the porosity of mud stone~depth relationship equation are as follows:Mud stone:Φ=0.5947exp (- 0.000762*h) (2)Sandstone:Φ=0.4906exp (- 0.000442*h) (3)2. the sand shale thickness of statistical number mouth well respectively, and according to skeleton thickness iterative Integral Method mathematical model, gone by stratum Compaction correction technology, is corrected to actual formation thickness (T6) formation thickness (T7) during embryo deposit;The thickness of the sand shale of depressed area selected by statistics respectively, will practically using skeleton thickness iterative Integral Method mathematical model Layer thickness T6 is corrected to embryo deposit period formation thickness T7;The Computing Principle of skeleton thickness iterative Integral Method mathematical model is such as Under, depth of stratum is calculated first as the weighting porosity at Z:Ф (Z)=PsФs(Z)+PmФm(Z) (4)In formula;Ps and Pm is the sandstone on stratum and the content of mud stone respectively, decimal;Ф s (Z) and Ф m (Z) are that depth is Z respectively The sandstone at place and the porosity of mud stone;Secondly, using skeleton thickness iterative Integral Method mathematical model, primary deposit period stratum will be reverted to by formation thickness now Thickness:<mrow> <mi>H</mi> <mn>7</mn> <mo>=</mo> <msubsup> <mo>&Integral;</mo> <mrow> <mi>Z</mi> <mn>1</mn> </mrow> <mrow> <mi>Z</mi> <mn>2</mn> </mrow> </msubsup> <mo>&lsqb;</mo> <mn>1</mn> <mo>-</mo> <mi>&Phi;</mi> <mrow> <mo>(</mo> <mi>Z</mi> <mo>)</mo> </mrow> <mo>&rsqb;</mo> <mi>d</mi> <mi>Z</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>In formula:T7 is embryo deposit period formation thickness, and Ф (Z) is that depth is weighting porosity at Z, and Z1 and Z2 are respectively The top circle of layer and the depth on bottom circle;Bring formula (2)-(4) into formula (5) and integrate, utilize above formula, you can obtain embryo deposit period formation thickness T7;So Afterwards, further according to well log interpretation data, the percentage composition of depressed area sandstone selected by difference, then selected depression is recovered using the above method The actual formation thickness of area's number mouthful well, and count formation compaction rate (compacting rate=embryo deposit formation thickness T7/ practically thickness Degree correction T6), the formation compaction rate of selected depressed area is generally in 1.8-2.0 scopes.
- 7. the quantitative approach according to claim 6 recovered based on lake basin paleotopography, it is characterised in that:The skeleton thickness Iterative Integral Method mathematical model has following premise when establishing:In whole deposition process, the skeleton volume on stratum does not change; And formation compaction process does not reverse.
- 8. the quantitative approach according to claim 6 recovered based on lake basin paleotopography, it is characterised in that:6th step Specific practice:1. by core analysis lithofacies characteristics, grain graininess, rock deposition construction and ground mark etc. are observed to qualitatively judge The scope of the depth of water;2. being directed to lake depositional environment, " beach dam Thickness Method " Fine calibration depth of water is proposed;Developed in the hydrodynamic force band of lake The beach dam of a variety of origin causes of formation, corresponds respectively to littoral dam, offshore dam and remote bank dam;First according to lithofacies characteristics and beach dam forming position Judge its origin cause of formation, then, the depth of water herein can be similar to further according to the deposit thickness on beach dam of single phase time;3. by the thickness after several Kou Jing beaches dam thickness correction, i.e.,:For paleao-water depth T8 herein;Counted several mouthfuls of well point depth of water T8 with just Relation between beginning formation thickness T7, and according to relevance formula, thickness during by selected depressed area embryo deposit are converted to water Deep quantitative distribution, is:Quantitative reconstruction lake basin landform.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711213606.3A CN107991714B (en) | 2017-11-28 | 2017-11-28 | Quantification method based on lake basin ancient terrain recovery |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711213606.3A CN107991714B (en) | 2017-11-28 | 2017-11-28 | Quantification method based on lake basin ancient terrain recovery |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107991714A true CN107991714A (en) | 2018-05-04 |
CN107991714B CN107991714B (en) | 2020-02-28 |
Family
ID=62032423
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711213606.3A Active CN107991714B (en) | 2017-11-28 | 2017-11-28 | Quantification method based on lake basin ancient terrain recovery |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107991714B (en) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109709302A (en) * | 2018-11-30 | 2019-05-03 | 中国海洋石油集团有限公司 | The method of clastic rock material resource system is differentiated based on multi-parameter comprehensive |
CN110187391A (en) * | 2019-06-20 | 2019-08-30 | 克拉玛依亚新石油科技有限公司 | A method of stratigraphic sequence interface is determined using interval transit time |
CN111045109A (en) * | 2018-10-11 | 2020-04-21 | 中国石油化工股份有限公司 | Method for recovering stratum denudation amount in basin |
CN111474601A (en) * | 2019-01-24 | 2020-07-31 | 中国石油天然气股份有限公司 | Method for predicting volcanic rock high-quality reservoir stage by stage |
CN111766630A (en) * | 2020-07-02 | 2020-10-13 | 中国地质大学(北京) | Restoration method for ancient landform of basin |
CN111983671A (en) * | 2019-05-23 | 2020-11-24 | 中国石油天然气股份有限公司 | Shallow water lake basin reservoir prediction method and device based on micro-ancient landform restoration |
CN113008921A (en) * | 2021-02-04 | 2021-06-22 | 新疆大学 | Method for determining ancient wind direction of basin through precipitation differentiation phenomena on two sides of lake center island |
CN113589365A (en) * | 2020-04-30 | 2021-11-02 | 中国石油化工股份有限公司 | Reservoir pinch-out line description method based on time-frequency domain information |
CN113625342A (en) * | 2021-08-04 | 2021-11-09 | 北京师范大学 | Method for determining denudation thickness of old stratum of clastic rock |
CN114526066A (en) * | 2022-01-24 | 2022-05-24 | 西南石油大学 | Method for compacting and recovering continental facies shale oil reservoir |
CN114624775A (en) * | 2020-12-11 | 2022-06-14 | 中国石油化工股份有限公司 | Comprehensive quantitative restoration method for ancient water depth of sedimentary lake basin |
CN117851748A (en) * | 2024-03-07 | 2024-04-09 | 中国石油大学(华东) | Method and system for calculating paleo-water depth of low exploration degree region |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104932031A (en) * | 2014-03-19 | 2015-09-23 | 中国石油化工股份有限公司 | Paleo-water-depth quantitative calculation method aiming at lake facies deposition |
CN105204069A (en) * | 2014-06-27 | 2015-12-30 | 中国石油化工股份有限公司 | Strata denudation amount recovery method |
CN105425316A (en) * | 2015-10-28 | 2016-03-23 | 中国石油化工股份有限公司 | Ancient landform restoration method based on impression method |
CN106094030A (en) * | 2016-08-24 | 2016-11-09 | 青岛海洋地质研究所 | A kind of method by seismic profile quantitative reconstruction lake basin maximum paleao-water depth |
-
2017
- 2017-11-28 CN CN201711213606.3A patent/CN107991714B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104932031A (en) * | 2014-03-19 | 2015-09-23 | 中国石油化工股份有限公司 | Paleo-water-depth quantitative calculation method aiming at lake facies deposition |
CN105204069A (en) * | 2014-06-27 | 2015-12-30 | 中国石油化工股份有限公司 | Strata denudation amount recovery method |
CN105425316A (en) * | 2015-10-28 | 2016-03-23 | 中国石油化工股份有限公司 | Ancient landform restoration method based on impression method |
CN106094030A (en) * | 2016-08-24 | 2016-11-09 | 青岛海洋地质研究所 | A kind of method by seismic profile quantitative reconstruction lake basin maximum paleao-water depth |
Non-Patent Citations (3)
Title |
---|
王彬 等: ""潜江凹陷潜江组地层剥蚀厚度恢复研究"", 《重庆石油高等专科学校学院》 * |
董刚 等: ""根据地层厚度恢复古水深的研究"", 《长江大学学报(自然科学版)》 * |
袁帅: ""博兴洼陷重点区块沙四上亚段滩坝薄互层砂体形成与分布预测"", 《中国优秀硕士学位论文数据库-基础科学辑》 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111045109A (en) * | 2018-10-11 | 2020-04-21 | 中国石油化工股份有限公司 | Method for recovering stratum denudation amount in basin |
CN111045109B (en) * | 2018-10-11 | 2023-06-02 | 中国石油化工股份有限公司 | Method for recovering stratum erosion amount in basin |
CN109709302A (en) * | 2018-11-30 | 2019-05-03 | 中国海洋石油集团有限公司 | The method of clastic rock material resource system is differentiated based on multi-parameter comprehensive |
CN111474601A (en) * | 2019-01-24 | 2020-07-31 | 中国石油天然气股份有限公司 | Method for predicting volcanic rock high-quality reservoir stage by stage |
CN111983671A (en) * | 2019-05-23 | 2020-11-24 | 中国石油天然气股份有限公司 | Shallow water lake basin reservoir prediction method and device based on micro-ancient landform restoration |
CN111983671B (en) * | 2019-05-23 | 2023-04-07 | 中国石油天然气股份有限公司 | Shallow water lake basin reservoir prediction method and device based on micro-ancient landform restoration |
CN110187391A (en) * | 2019-06-20 | 2019-08-30 | 克拉玛依亚新石油科技有限公司 | A method of stratigraphic sequence interface is determined using interval transit time |
CN110187391B (en) * | 2019-06-20 | 2022-02-01 | 克拉玛依亚新石油科技有限公司 | Method for determining stratum sequence interface by using acoustic wave time difference |
CN113589365A (en) * | 2020-04-30 | 2021-11-02 | 中国石油化工股份有限公司 | Reservoir pinch-out line description method based on time-frequency domain information |
CN111766630A (en) * | 2020-07-02 | 2020-10-13 | 中国地质大学(北京) | Restoration method for ancient landform of basin |
CN111766630B (en) * | 2020-07-02 | 2022-04-22 | 中国地质大学(北京) | Restoration method for ancient landform of basin |
CN114624775A (en) * | 2020-12-11 | 2022-06-14 | 中国石油化工股份有限公司 | Comprehensive quantitative restoration method for ancient water depth of sedimentary lake basin |
CN114624775B (en) * | 2020-12-11 | 2023-04-07 | 中国石油化工股份有限公司 | Comprehensive quantitative restoration method for ancient water depth of sedimentary lake basin |
CN113008921A (en) * | 2021-02-04 | 2021-06-22 | 新疆大学 | Method for determining ancient wind direction of basin through precipitation differentiation phenomena on two sides of lake center island |
CN113625342A (en) * | 2021-08-04 | 2021-11-09 | 北京师范大学 | Method for determining denudation thickness of old stratum of clastic rock |
CN114526066A (en) * | 2022-01-24 | 2022-05-24 | 西南石油大学 | Method for compacting and recovering continental facies shale oil reservoir |
CN114526066B (en) * | 2022-01-24 | 2024-06-11 | 西南石油大学 | Method for recovering compaction of land shale oil reservoir |
CN117851748A (en) * | 2024-03-07 | 2024-04-09 | 中国石油大学(华东) | Method and system for calculating paleo-water depth of low exploration degree region |
Also Published As
Publication number | Publication date |
---|---|
CN107991714B (en) | 2020-02-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107991714A (en) | The quantitative approach recovered based on lake basin paleotopography | |
CN106875471B (en) | Three-dimensional visual modeling method for coal-series water-containing or water-resisting layer | |
CN105373648B (en) | A kind of modeling method of sand-mud interbed type lithosomic body reservoir architecture | |
CN104632200B (en) | Fluvial facies multilayer and thick-layer oil reservoir fine modeling method | |
CN105425315B (en) | A kind of small scale depositional phase ancient landform inversion method in weak structural deformation area | |
CN106154351A (en) | A kind of evaluation method of low porosity permeability reservoir permeability | |
CN105822302A (en) | Oil-water distribution recognition method based on well-ground potential method | |
CN105137482B (en) | A kind of computational methods of the lithosomic body Gu gradient | |
CN105370269B (en) | One kind tracking control methods based on compaction correction chronostratigraphic interface | |
CN110056346B (en) | Oil reservoir three-dimensional original water saturation simulation method based on trend change function | |
CN102243678A (en) | Method for analyzing sand bodies in reservoirs based on inversion technique of sedimentary dynamics | |
CN106094030A (en) | A kind of method by seismic profile quantitative reconstruction lake basin maximum paleao-water depth | |
CN107066749B (en) | A kind of method of quantitative assessment Seam Roof And Floor capping performance | |
CN106570262A (en) | Reservoir configuration structure description method | |
CN104834003B (en) | Phase-controlled compression coefficient earthquake prediction method of unconventional tight gas reservoir | |
CN105868482A (en) | Method and device for calculating depositional stage micro-amplitude palaeohigh | |
CN110424956A (en) | Evaluation unit saves coefficient weights quantization assignment method in shale oil Resources calculation | |
CN115220101A (en) | Modeling method for deep clastic rock sand body sediment structure | |
Niu et al. | Hierarchical modeling method based on multilevel architecture surface restriction and its application in point-bar internal architecture of a complex meandering river | |
CN108267797B (en) | Method for restoring ancient landform of syngeneic normal fault development area | |
Xia et al. | The controlling factors of modern facies distributions in a half‐graben lacustrine rift basin: A case study from Hulun Lake, Northeastern China | |
Liu et al. | The Control Theory and Application for Well Pattern Optimization of Heterogeneous Sandstone Reservoirs | |
CN110428497A (en) | Braided stream training image generation method | |
Guderian et al. | Draugen field–successful reservoir management using 4d seismic | |
Siddiqui et al. | Facies characteristics and static reservoir connectivity of some siliciclastic tertiary outcrop successions in Bintulu and Miri, Sarawak, East Malaysia |
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 |