CN106526679A - Method for recovering non-uniform erosion thickness through utilizing cycle analysis method - Google Patents
Method for recovering non-uniform erosion thickness through utilizing cycle analysis method Download PDFInfo
- Publication number
- CN106526679A CN106526679A CN201710015647.5A CN201710015647A CN106526679A CN 106526679 A CN106526679 A CN 106526679A CN 201710015647 A CN201710015647 A CN 201710015647A CN 106526679 A CN106526679 A CN 106526679A
- Authority
- CN
- China
- Prior art keywords
- thickness
- unconformity
- stratum
- erosion
- sedimentation rate
- 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
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/614—Synthetically generated data
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geology (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Medicinal Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- Food Science & Technology (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Acoustics & Sound (AREA)
- Geophysics (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
The invention provides a method for recovering non-uniform erosion thickness through utilizing a cycle analysis method and belongs to the technical field of geological exploration. The method provided by the invention comprises the following steps: calculating a compaction coefficient and an original thickness of a stratum before compaction through decompaction correction; calculating an average stratum deposition speed of each deposition period and drawing a deposition rate curve; determining a missed curve section according to the deposition rate curve; reckoning a period corresponding to the missed curve section and transverse coordinates corresponding to a wave peak and a wave trough of the missed curve section; carrying out frequency-divisional processing on the deposition rate curve to obtain a low-frequency wavelet and a high-frequency wavelet, reckoning longitudinal coordinates corresponding to the wave peak and the wave trough of the missed curve section; and recovering the complete deposition rate curve according to frequency spectrum analysis, and calculating the non-uniform erosion thickness through complete deposition rate curve integral. According to the method provided by the invention, a stratum deposition-erosion state of a whole geological historical period can be effectively recovered, the erosion thickness of different erosion stratum layers is recovered, and the method is high in accuracy and wide in applicable range.
Description
Technical field
The invention belongs to technical field of geological exploration, more particularly to a kind of recover unconformity using Cyclic Analysis method and degrade thickness
The method of degree.
Background technology
Recover the important content that unconformity erosion sediment thickness is basin-mountain frame work analysis and Hydrocarbon Formation Reservoirs research, its result
Accuracy is to the evaluation of Thermal Evolution of Source Rocks, reservoir development and transformation, Gas Accumulation and preservation condition with important impact.
The method for recovering unconformity erosion sediment thickness at present is more, including mud stone sonic time difference, vitrinite reflectance rate, bag
Wrap up in body homogenization temperature method, and based on the Analysis of Earthquake Trend method of seismic interpretation, adjacent bed thickness ratio method, refer to layer thickness variation
Rate method etc..However, practice finds that these methods still suffer from some shortcomings：(1) as above-mentioned existing method is by setting up one
Plant mathematical relationship to realize the quantitative reconstruction to erosion sediment thickness, and the foundation of this mathematical relationship often has nonuniqueness, from
And cause erosion sediment thickness value that there is multi-solution, and error is larger；(2) above-mentioned existing method has respective applicable elements and limitation
Property, for example：Vitrinite reflectance rate needs more test data, and is affected fairly obvious by incident heat；Based on seismic interpretation
Analysis of Earthquake Trend rule need to set up research area degrade scope it is less on the basis of, for degrade scope cover the whole district ground
Area does not apply to；(3) as unconformity is probably that a series of strata degrade product, it is also possible to which multiple series of strata degrade product
Thing, when denudation degree is larger, multiple series of strata degrade to form a unconformity, is only capable of recovery using above-mentioned existing method a set of not
The erosion sediment thickness of integration, i.e., many series of strata erosion sediment thickness sums, it is impossible to which the erosion sediment thickness of different series of strata is recovered.
Therefore, how to provide a kind of degree of accuracy it is high, applied widely and can recover difference degrade series of strata erosion sediment thickness
Method is the current technical barrier be badly in need of and solving.
The content of the invention
Method error for above-mentioned existing recovery unconformity erosion sediment thickness of the invention is big, the scope of application has limitation
Property and the technical problem that the erosion sediment thickness of different series of strata cannot be recovered, propose a kind of to recover not using Cyclic Analysis method
The method for integrating erosion sediment thickness, can effectively reduce sedimentation-degrade state of the research area in the whole geology period of history, can
Recover the erosion sediment thickness that difference degrades series of strata, and the method degree of accuracy is high, applied widely.
In order to achieve the above object, the technical solution used in the present invention is：
A kind of method that utilization Cyclic Analysis method recovers unconformity erosion sediment thickness, comprises the following steps：
Sedimentation rate Drawing of Curve：Research area is chosen, the geologic information in binding area is calculated by decompaction correction and pressed
Real coefficient and the original thickness on the front stratum of compacting, determine the initial time of each sedimentation with the end of according to stratigraphical time table
Between, the average stratum sedimentation rate of each deposition phase is calculated, sedimentation rate curve is drawn；
The spectrum analysis of sedimentation rate curve：Disappearance curved section is determined according to the cycle rule of sedimentation rate curve；According to
The periodic quantity of sedimentation rate curve before and after disappearance curved section, calculate disappearance curved section corresponding cycle, disappearance curved section crest and
The corresponding abscissa of trough；Scaling down processing is carried out to sedimentation rate curve and obtains correspondence lower frequency wavelet and high-frequency wavelet, according to scarce
Lose curved section before and after corresponding lower frequency wavelet and high-frequency wavelet crest and trough ordinate, calculate disappearance curved section crest and
The corresponding ordinate of trough；
Calculate unconformity erosion sediment thickness：Complete sedimentation rate curve is recovered according to result of spectrum analysis, unconformity is calculated
Erosion sediment thickness, the computing formula of the unconformity erosion sediment thickness are as follows：
Wherein, HeFor unconformity erosion sediment thickness, unit：m；F (t) is deposited intact rate curve function；HocIt is in deposition
Original thickness F (t) integrations obtained in time period；HrFor residual thickness now, unit：m；C is total compacting factor；t1
To degrade sedimentation time started, unit：Ma；t2To degrade sedimentation end time, unit：Ma；
Preferably, in sedimentation rate Drawing of Curve step, the computing formula of the original thickness for being compacted front stratum
For：
H0=Hs/Cs+Hm/Cm
Wherein, H0For the original thickness before formation compaction, unit：m；HsFor sandstone thickness, unit：m；CsFor sandstone compaction
Coefficient, is calculated using decompaction correction method and is obtained；HmFor mud stone thickness, unit：m；CmFor mudstone compacting coefficient, using going pressure
Real bearing calibration is calculated and is obtained；
When unconformity erosion sediment thickness is calculated, the computing formula of total compacting factor is：
Wherein, C is total compacting factor；HsFor sandstone thickness, unit：m；CsFor sandstone compaction coefficient, using decompaction correction
Method is calculated and is obtained；HmFor mud stone thickness, unit：m；CmFor mudstone compacting coefficient, calculated using decompaction correction method and obtained.
Preferably, the geologic information includes strata division data, stratigraphical time table, and log data or porosity
Data.
Preferably, the decompaction correction method is empirical formula method or porosity-lithology method；When the geology is provided
When material includes log data, the decompaction correction method adopts empirical formula method；When the geologic information includes that porosity is provided
During material, the decompaction correction method adopts porosity-lithology method.
Preferably, after unconformity erosion sediment thickness step is calculated, also including using strata division data or existing recovery
Obtained by the method validation of unconformity erosion sediment thickness the step of unconformity erosion sediment thickness accuracy.
Preferably, carrying out concretely comprising the following steps for Accuracy Verification using strata division data：It is complete heavy using what is recovered
Product rate curve, is integrated by the corresponding sedimentation rate curved section in stratum to degrading and is not degraded stratum
Original thickness, by compacting be calculated do not degrade after the compacting on stratum in formation thickness, with strata division data now
Thickness degree is contrasted, and verifies accuracy.
Preferably, the existing method for recovering unconformity erosion sediment thickness is anti-selected from mud stone sonic time difference, vitrinite
Penetrate rate method, Evidence from Apatite Fission Track Analysis, material balance method, inclusion enclave homogenization temperature method, sedimentary ratio analysis, sandstone porosity method,
Any one in wave theory approach, Analysis of Earthquake Trend method, adjacent bed thickness ratio method and reference layer thickness change method.
Compared with prior art, advantages of the present invention and good effect are：
1st, the present invention can effectively reduce sedimentation-degrade state of the research area in the whole geology period of history, recover
Sedimentation original thickness in any one time range, degrades the erosion sediment thickness of series of strata, Ke Yiwei so as to recover difference
Important foundation is established in the research such as migration direction of structural evolution, oil gas in the different periods of history, advantageously reduces oil-gas exploration wind
Danger；
It is 2nd, of the invention merely with basic datas such as strata division data, stratigraphical time table, log data and porosity data,
Requirement to data is simple and data volume that need is less, eliminates because experimental data quantity is few or machine error, sampling difference
Etc. the error brought, accuracy is improve；
3rd, the present invention for research area degrade scope or strata pressure etc. without it is specifically limited with requirement, the scope of application is more
For extensive；
4th, the present invention does not rely on experiment, it is not necessary to the support of special experimental facilities, simple to operate and cost is lower.
Description of the drawings
Fig. 1 is recovered the flow process of the method for unconformity erosion sediment thickness by the utilization Cyclic Analysis method that the embodiment of the present invention is provided
Figure；
Sedimentation rate broken line graph before the B202 wells recovery that Fig. 2 is provided by the embodiment of the present invention 1；
Sedimentation rate curve map before the B202 wells recovery that Fig. 3 is provided by the embodiment of the present invention 1；
Fig. 4 recovers the corresponding lower frequency wavelet curve of front sedimentation rate curve by the B202 wells that the embodiment of the present invention 1 is provided
Figure；
Fig. 5 recovers the corresponding high-frequency wavelet curve of front sedimentation rate curve by the B202 wells that the embodiment of the present invention 1 is provided
Figure；
Fig. 6 is by the complete sedimentation rate broken line graph of B202 wells that the embodiment of the present invention 1 is provided；
Fig. 7 is by the complete sedimentation rate curve map of B202 wells that the embodiment of the present invention 1 is provided；
The corresponding lower frequency wavelet curve map of B202 well deposited intact rate curves that Fig. 8 is provided by the embodiment of the present invention 1；
The corresponding high-frequency wavelet curve map of B202 well deposited intact rate curves that Fig. 9 is provided by the embodiment of the present invention 1.
Specific embodiment
Technical scheme in the embodiment of the present invention will be clearly and completely described below, it is clear that described enforcement
Example is only a part of embodiment of the invention, rather than the embodiment of whole.Based on the embodiment in the present invention, this area is common
The every other embodiment obtained under the premise of creative work is not made by technical staff, belongs to the model of present invention protection
Enclose.
Embodiments provide a kind of method that utilization Cyclic Analysis method recovers unconformity erosion sediment thickness, its flow chart
As shown in figure 1, comprising the following steps：
S1 sedimentation rate Drawing of Curve：Research area is chosen, the geologic information in binding area is calculated by decompaction correction
Compacting factor and the original thickness on the front stratum of compacting, determine initial time and the end of each sedimentation according to stratigraphical time table
Time, the average stratum sedimentation rate of each deposition phase is calculated, draw sedimentation rate curve.
In this step, it should be noted that due to there is compacting to make the stratum of first deposition on the stratum of rear deposition
With the original thickness for being compacted front stratum being obtained by decompaction correction, using the calculated average stratum of this original thickness
Sedimentation rate is closer to the deposition phase real sedimentation speed.The computing formula of the average stratum sedimentation rate is：
V=H0/(t4-t3) (1)
In formula (1), V is the average stratum sedimentation rate of each deposition phase, unit：m/Ma；H0It is original before for formation compaction
Thickness, unit：m；t3For each sedimentation time started deposition phase, unit：Ma；t4For each sedimentation end time deposition phase,
Unit：Ma.
In addition it is also necessary to explanation, the drafting of sedimentation rate curve can be carried out using Matlab softwares, with the time as horizontal stroke
Coordinate, average stratum sedimentation rate are that ordinate draws sedimentation rate broken line graph, are further obtained by cubic spline interpolation method
To sedimentation rate curve map.
The spectrum analysis of S2 sedimentation rate curves：Disappearance curved section is determined according to the cycle rule of sedimentation rate curve；Root
According to the periodic quantity of sedimentation rate curve before and after disappearance curved section, disappearance curved section corresponding cycle, disappearance curved section crest are calculated
With the corresponding abscissa of trough；Scaling down processing is carried out to sedimentation rate curve, the corresponding lower frequency wavelet of sedimentation rate curve is obtained
And high-frequency wavelet, disappearance curved section crest and the corresponding ordinate of trough are calculated according to lower frequency wavelet and high-frequency wavelet.
In this step, it should be noted that as the deposition on stratum is with certain cyclicity, concrete manifestation
For：The sedimentation rate on stratum fluctuates near sinusoidal, and its cycle, crest and trough numerical value are presented certain gradual change law.So
And, degrading for stratum can cause deposition rate data to lack, and then affect this cycle rule of sedimentation rate curve, therefore,
Can determine that curved section is as lacked with one section of inconsistent curved section of cycle rule in front and back in sedimentation rate curve.Additionally, knot
Close areal geology data and geophysical response (i.e. log mutation, identification of earthquake unconformity etc.), it may be verified that disappearance curve
The accuracy of section correspondence time interval.
Meanwhile, cycle of sedimentation rate curve, crest and trough numerical value are also presented certain linear gradient rule, that is, deposit
Cycle of rate curve, crest and trough numerical value taper off or progressive law.And, the disappearance curved section corresponding time must be situated between
Between anterioposterior curve, i.e. time-constrain principle.Therefore, using this linear gradient rule and time-constrain principle, according to scarce
The periodic quantity of sedimentation rate curve before and after mistake curved section, can extrapolate the periodicity of disappearance curved section disappearance, can further calculate
Go out to lack curved section corresponding cycle, disappearance curved section crest and the corresponding abscissa of trough.
In addition, although disappearance curved section crest and the corresponding ordinate of trough (i.e. sedimentation rate) numerical value also can be carried out in the same manner
Calculate, but scaling down processing is carried out to sedimentation rate curve using decomposed and reconstituted thought, the fluctuation change of the high-frequency wavelet for obtaining can
The cycle rule of reflection deposition, the fluctuation change of lower frequency wavelet can reflect the approximate trend of deposition, thus calculate the data for obtaining
More accuracy, is conducive to improving the accuracy of spectrum analysis.Specifically, scaling down processing is carried out to sedimentation rate curve and obtains low
Frequency wavelet and high-frequency wavelet, read the corresponding crest value of two wavelets respectively, and the two superposition is obtained total crest numerical value Yc,
Trough numerical value Yt can be obtained in the same manner.By the method calculate respectively disappearance curved section before and after the corresponding crest of sedimentation rate curve and
Trough numerical value, recycles linear gradient rule, you can extrapolate disappearance curved section crest and the corresponding ordinate of trough.
Also, it should be noted that scaling down processing can be completed in Matlab softwares, concretely comprise the following steps：By (i-N+1)
Secondary N points take the lower frequency wavelet that average method obtains sedimentation rate curve, sedimentation rate curve is deducted the lower frequency wavelet and is obtained
High-frequency wavelet.Wherein, N points take average method and are specially：For any one-dimensional data x, solid timing window is chosenFor
Each sampled point x in datai, progressively sliding window averaged, the process of asking for are represented by：
In formula (2), i is data point number；tiTime point corresponding to i-th data point；N is data sampling point
Number；J is cyclic variable；Δ t is time sampling interval.
S3 calculates unconformity erosion sediment thickness：Complete sedimentation rate curve is recovered according to result of spectrum analysis, is calculated not whole
Erosion sediment thickness is closed, the computing formula of the unconformity erosion sediment thickness is as follows：
In this step, disappearance curved section corresponding cycle, disappearance curved section crest and the ripple for being obtained according to spectrum analysis
The corresponding abscissa of paddy and ordinate, you can recover complete sedimentation rate curve, are sedimentation rate maximum at its medium wave peak
In period, be the period that sedimentation rate is minimum at trough, and the curve can truly reflect the cycle change of sedimentation rate and sedimentary evolution
Process.Further, by integrating to deposited intact rate curve, you can obtain the original thickness not by compaction, profit
Compacting calculating is carried out to the thickness with the total compacting factor for calculating acquisition before, you can obtain the original thickness on stratum, and degrade
Thickness is the original thickness on stratum and the difference of residual thickness now.Therefore, unconformity can be calculated using formula (3) and degrades thickness
Degree.It should be noted that parameter t in formula (3)1And t2Be for a certain deposition start and end time for specifically degrading stratum,
Thus the method that the utilization Cyclic Analysis method of the present embodiment offer recovers unconformity erosion sediment thickness can degrade stratum for difference
Individually calculated, and then the erosion sediment thickness that difference degrades series of strata can be recovered.
In a preferred embodiment, in sedimentation rate Drawing of Curve step, the original thickness for being compacted front stratum
Computing formula is：
H0=Hs/Cs+Hm/Cm (4)
When unconformity erosion sediment thickness is calculated, the computing formula of total compacting factor is：
In the present embodiment, the computing formula and total compacting factor of the original thickness for being compacted front stratum are further defined
Computing formula, it should be noted that as sandstone and/or mud stone are included in stratum, and sandstone is different from the compacting factor of mud stone,
Accordingly, it would be desirable to carry out decompaction correction to sandstone and mud stone respectively.
In a preferred embodiment, the geologic information includes strata division data, stratigraphical time table, and log data
Or porosity data.The present embodiment further defines the geologic information combination of optimum, it is to be understood that the present embodiment is not
It is confined to above-mentioned cited geologic information species, the conventional geology money of other that can also be well known to those skilled in the art
Material.
In a preferred embodiment, the decompaction correction method is empirical formula method or porosity-lithology method；Work as institute
When stating geologic information and including log data, the decompaction correction method adopts empirical formula method；When the geologic information includes
During porosity data, the decompaction correction method adopts porosity-lithology method.The present embodiment further defines compacting
Bearing calibration, it should be noted that those skilled in the art can be specifically chosen according to geology data type and stratum rock signature
Suitable decompaction correction method.And, the present embodiment is not limited to above-mentioned cited decompaction correction method, can be with
Other the conventional decompaction correction methods being well known to those skilled in the art.
In addition it is also necessary to explanation, using the formula of empirical formula method calculating compacting factor is：
C'=b-alnH (6)
In formula (6), when sandstone compaction coefficient is calculated, C ' is Cs；When mudstone compacting coefficient is calculated, C ' is Cm；a、
B is constant, for sandstone takes a=0.08, b=1.37, for mud stone takes a=0.1, b=1.46；H is buried depth, unit：
m。
Compacting factor is calculated using the porosity-lithology method, the i.e. relation using different lithology porosity with depth
Formula be：
In formula (7), when sandstone compaction coefficient is calculated, C " it is Cs；When mudstone compacting coefficient is calculated, C " it is Cm；H
For buried depth, unit：m；It is the porosity at H for depth；For earth's surface initial porosity.
In a preferred embodiment, after unconformity erosion sediment thickness step is calculated, also include using strata division data or
Obtained by the existing method validation for recovering unconformity erosion sediment thickness the step of unconformity erosion sediment thickness accuracy, by the verification step
Can ensure that the accuracy of gained unconformity erosion sediment thickness.
In a preferred embodiment, concretely comprising the following steps for Accuracy Verification is carried out using strata division data：Using recovery
Deposited intact rate curve, be integrated by the corresponding sedimentation rate curved section in stratum to degrading and do not shelled
The original thickness on erosion stratum, by being compacted be calculated after the compacting for not degrading stratum in formation thickness, with strata division data
Formation thickness now contrasted, verify accuracy.In the present embodiment, the formation thickness after the compacting for not degrading stratum
With the error size of formation thickness now, you can the accuracy of checking gained unconformity erosion sediment thickness, this verification method are not required to
Want the support of additive method and data, simple to operate, workload is little and low cost.
In a preferred embodiment, it is described it is existing recover unconformity erosion sediment thickness method selected from mud stone sonic time difference,
Vitrinite reflectance rate, Evidence from Apatite Fission Track Analysis, material balance method, inclusion enclave homogenization temperature method, sedimentary ratio analysis, sandstone hole
It is arbitrary in porosity method, wave theory approach, Analysis of Earthquake Trend method, adjacent bed thickness ratio method and reference layer thickness change method
Kind.In the present embodiment, it should be noted that those skilled in the art can be degraded according to geologic information and existing recovery unconformity
The applicable elements of thickness approach, specifically chosen suitable existing recovery unconformity erosion sediment thickness method are verified.And, this reality
The verification method for applying example employing is not limited to above-mentioned cited recovery unconformity erosion sediment thickness method, can also be this area
Other common methods known to technical staff.
The utilization Cyclic Analysis method provided to become apparent from introducing in detail the embodiment of the present invention is recovered unconformity and is degraded
The method of thickness, is described below in conjunction with specific embodiment.
Embodiment 1
In Song-liao basin Chang Ling rift dragon and phoenix Mountain area, a bite typical case well B202 is chosen, the well location is in dragon and phoenix Mountain area
Southeast gentle slope belt, the Cretaceous Strata that drilling well is got into are followed successively by open fire group, Sifangtai group, Nenjiang group, Yao Jia groups, green grass or young crops from top to bottom
Mountain pass group, spring head group, step on Lou Ku groups, battalion city group, Shahe subgroup (not drilling).
The unconformity erosion sediment thickness of B202 wells is recovered, step is as follows：
S1 sedimentation rate Drawing of Curve：
Strata division data, stratigraphical time table and log data with reference to B202 wells, carry out pressure by empirical formula method
It is real to correct, compacting factor is calculated using formula (6), calculate the original thickness on the front stratum of compacting using formula (4), according to stratigraphical time table
Determine initial time and the end time of each sedimentation, the flat of each stratum deposition period of the well is calculated using formula (1)
Equal sedimentation speed, data are as shown in table 1.
The average stratum sedimentation rate statistics of each stratum deposition period of 1 B202 wells of table
The average stratum sedimentation rate statistics of each stratum deposition period of 1 (Continued) B202 well of table
According to the result of calculation of table 1, using Matlab softwares with the time as abscissa, average stratum sedimentation rate be vertical seat
Mark, draws the sedimentation rate broken line graph of B202 wells, as shown in Figure 2.Sedimentation rate curve is obtained by cubic spline interpolation method
Figure, as shown in Figure 3.
The spectrum analysis of S2 sedimentation rate curves：
According to the cycle rule of sedimentation rate curve in Fig. 3, it is apparent that lacking between the 7th point and the 8th point
Data are lost, it is determined that the corresponding time interval of disappearance curved section is 109.835Ma-118Ma.Can according to the spectrum analysis of Fig. 2 and Fig. 3
Know, lack curved section before (i.e. in the 94.25Ma-109.835Ma times) curve be 2 cycles, its average period T1=
7.7925Ma；Disappearance curved section after (i.e. in the 118Ma-135.1429Ma times) curve be 2.75 cycles, its average week
Phase T2=6.234Ma.Thus can obtain, the cycle of sedimentation rate curve is as the age is from newly to the rule that tapers off always.In addition, according to
In Fig. 3, the cycle rule of sedimentation rate curve is visible, and disappearance curved section should be 1.25 cycles, that is, lack a crest and one
Trough, therefore, the average period for lacking curved section is 6.532Ma, according to linear gradient and time-constrain principle, lacks curved section
The corresponding abscissa of crest (i.e. time value) is 113.0605Ma, and the corresponding abscissa of trough is 116.286Ma.
Scaling down processing is carried out to the sedimentation rate curve in Fig. 3, lower frequency wavelet is extracted using formula (2), as shown in figure 4, will
Sedimentation rate curve deducts the lower frequency wavelet and obtains high-frequency wavelet, as shown in Figure 5.Read disappearance from Fig. 4 and Fig. 5 respectively bent
Lower frequency wavelet and the corresponding crest value of high-frequency wavelet before line segment, the two superposition obtain lacking the crest value Yc before curved section1
=81.79, obtain in the same manner lacking the valley value Yt before curved section1=43.87；In the same manner, obtain lack curved section after ripple
Peak value Yc2=148.68, valley value Yt2=39.01.There is linear gradient to advise using the crest and valley value of sedimentation rate curve
Rule, it is 115.235m/Ma that can extrapolate the disappearance corresponding ordinate (i.e. sedimentation rate) of curved section crest, the corresponding vertical seat of trough
It is designated as 41.44m/Ma.
S3 calculates unconformity erosion sediment thickness：
According to above result of spectrum analysis, will (113.0605,115.235) and (116.286,41.44) two data points
It is added in the data of Fig. 2, the deposited intact speed broken line graph of drafting is as shown in fig. 6, corresponding deposited intact rate profile
As shown in fig. 7, corresponding lower frequency wavelet and high-frequency wavelet curve map are respectively as shown in Figure 8, Figure 9.The ground that B202 wells degrade
To step on Lou Ku groups (112Ma-124Ma), it is 0.7207 to be calculated using formula (5) and step on the corresponding total compacting factor C of Lou Ku groups to layer,
By C=0.7207, t1=112Ma, t2=124Ma, Hr=318.8m is substituted in formula (3), is calculated the unconformity for stepping on Lou Ku groups
Erosion sediment thickness HeFor 379.69m.
S4 method validations：
Verified using strata division data, using the deposited intact rate curve for recovering, by degrading
The corresponding sedimentation rate curved section in stratum be integrated the original thickness for not degraded stratum, by compacting be calculated not
Formation thickness after the compacting on stratum is degraded, data are as shown in table 2.
Table 2 utilizes B202 well strata division data verification result statistical forms
After the compacting on stratum will not degraded, formation thickness is carried out with the formation thickness now obtained from strata division data
Contrast, as shown in table 2, in addition to Qingshankou group stratum, the error rate on remaining stratum is respectively less than 10%.Drilled according to Song-liao basin construction
The pertinent literature for changing history understands that Qingshankou group stratum occurred strata denudation, caused by the error exactly strata denudation on the stratum.
As can be seen here, it is reliable, effective that the utilization Cyclic Analysis method that the present embodiment is provided recovers the method for unconformity erosion sediment thickness.
Embodiment 2
With not being both for embodiment 1, the present embodiment is using mud stone sonic time difference checking gained unconformity erosion sediment thickness
Accuracy, in addition to method validation step, remaining step and its data obtained result are with embodiment 1.Due to mud stone sonic time difference
For the method for existing conventional recovery unconformity erosion sediment thickness, here is omitted.The present embodiment calculates the B202 wells of gained and steps on
The unconformity erosion sediment thickness of Lou Ku groups is 379.69m, adopts mud stone sonic time difference result of calculation for 408.19m, the two error
Rate is 6.98%, less than 10%, it can be seen that, the utilization Cyclic Analysis method that the present embodiment is provided recovers unconformity erosion sediment thickness
Method is reliable, effective.
Embodiment 3
The unconformity erosion sediment thickness of Song-liao basin Chang Ling rift dragon and phoenix Mountain area B204 wells is recovered, except method validation
Outside step, with embodiment 1, the unconformity erosion sediment thickness for being calculated B204 wells is 418.09m to remaining step, and the present embodiment is adopted
The accuracy of inclusion enclave homogenization temperature method checking gained unconformity erosion sediment thickness.As inclusion enclave homogenization temperature method is existing conventional
Recovery unconformity erosion sediment thickness method, here is omitted.The result calculated using first kind inclusion enclave homogenization temperature method
For 388.58m, error rate is 7.60%；The result that Equations of The Second Kind inclusion enclave homogenization temperature method is calculated is adopted for 446.15m, error rate
For 6.29%.As can be seen here, it is reliable that the utilization Cyclic Analysis method that the present embodiment is provided recovers the method for unconformity erosion sediment thickness
, effectively.
Claims (7)
1. a kind of method that utilization Cyclic Analysis method recovers unconformity erosion sediment thickness, it is characterised in that comprise the following steps：
Sedimentation rate Drawing of Curve：Research area is chosen, the geologic information in binding area calculates compacting system by decompaction correction
The original thickness on front stratum is counted and is compacted, initial time and the end time of each sedimentation are determined according to stratigraphical time table,
The average stratum sedimentation rate of each deposition phase is calculated, sedimentation rate curve is drawn；
The spectrum analysis of sedimentation rate curve：Disappearance curved section is determined according to the cycle rule of sedimentation rate curve；According to disappearance
The periodic quantity of sedimentation rate curve before and after curved section, calculates disappearance curved section corresponding cycle, disappearance curved section crest and trough
Corresponding abscissa；Scaling down processing is carried out to sedimentation rate curve and obtains correspondence lower frequency wavelet and high-frequency wavelet, it is bent according to disappearance
The ordinate of the crest and trough of corresponding lower frequency wavelet and high-frequency wavelet before and after line segment, calculates disappearance curved section crest and trough
Corresponding ordinate；
Calculate unconformity erosion sediment thickness：Complete sedimentation rate curve is recovered according to result of spectrum analysis, unconformity is calculated and is degraded
Thickness, the computing formula of the unconformity erosion sediment thickness are as follows：
Wherein, HeFor unconformity erosion sediment thickness, unit：m；F (t) is deposited intact rate curve function；HocIt is in sedimentation time
Original thickness F (t) integrations obtained in section；HrFor residual thickness now, unit：m；C is total compacting factor；t1For stripping
Erosion sedimentation time started, unit：Ma；t2To degrade sedimentation end time, unit：Ma.
2. the method that utilization Cyclic Analysis method according to claim 1 recovers unconformity erosion sediment thickness, it is characterised in that
In sedimentation rate Drawing of Curve step, the computing formula of the original thickness for being compacted front stratum is：
H0=Hs/Cs+Hm/Cm
Wherein, H0For the original thickness before formation compaction, unit：m；HsFor sandstone thickness, unit：m；CsFor sandstone compaction coefficient,
Calculated using decompaction correction method and obtained；HmFor mud stone thickness, unit：m；CmFor mudstone compacting coefficient, using decompaction correction
Method is calculated and is obtained；
When unconformity erosion sediment thickness is calculated, the computing formula of total compacting factor is：
Wherein, C is total compacting factor；HsFor sandstone thickness, unit：m；CsFor sandstone compaction coefficient, using decompaction correction method
Calculate and obtain；HmFor mud stone thickness, unit：m；CmFor mudstone compacting coefficient, calculated using decompaction correction method and obtained.
3. the method that utilization Cyclic Analysis method according to claim 1 recovers unconformity erosion sediment thickness, it is characterised in that：Institute
Stating geologic information includes strata division data, stratigraphical time table, and log data or porosity data.
4. the method that utilization Cyclic Analysis method according to claim 3 recovers unconformity erosion sediment thickness, it is characterised in that：Institute
Decompaction correction method is stated for empirical formula method or porosity-lithology method；When the geologic information includes log data, institute
State decompaction correction method and adopt empirical formula method；When the geologic information includes porosity data, the decompaction correction
Method adopts porosity-lithology method.
5. the method that utilization Cyclic Analysis method according to claim 1 recovers unconformity erosion sediment thickness, it is characterised in that：
After calculating unconformity erosion sediment thickness step, also include the method using strata division data or existing recovery unconformity erosion sediment thickness
The step of checking gained unconformity erosion sediment thickness accuracy.
6. the method that utilization Cyclic Analysis method according to claim 5 recovers unconformity erosion sediment thickness, it is characterised in that profit
Concretely comprising the following steps for Accuracy Verification is carried out with stratum individual-layer data：Using the deposited intact rate curve for recovering, by not
The corresponding sedimentation rate curved section in stratum for degrading is integrated the original thickness for not degraded stratum, by compacting meter
After the compacting on stratum is not degraded in calculation, formation thickness now in formation thickness, with strata division data is contrasted, and is verified
Accuracy.
7. the method that utilization Cyclic Analysis method according to claim 5 recovers unconformity erosion sediment thickness, it is characterised in that：Institute
State the existing method for recovering unconformity erosion sediment thickness and be selected from mud stone sonic time difference, vitrinite reflectance rate, apatite fission footpath
Mark method, material balance method, inclusion enclave homogenization temperature method, sedimentary ratio analysis, sandstone porosity method, wave theory approach, seismic trend
Any one in analytic approach, adjacent bed thickness ratio method and reference layer thickness change method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710015647.5A CN106526679B (en) | 2017-01-10 | 2017-01-10 | A kind of method for recovering unconformity erosion sediment thickness using Cyclic Analysis method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710015647.5A CN106526679B (en) | 2017-01-10 | 2017-01-10 | A kind of method for recovering unconformity erosion sediment thickness using Cyclic Analysis method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106526679A true CN106526679A (en) | 2017-03-22 |
| CN106526679B CN106526679B (en) | 2018-05-25 |
Family
ID=58335346
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710015647.5A Expired - Fee Related CN106526679B (en) | 2017-01-10 | 2017-01-10 | A kind of method for recovering unconformity erosion sediment thickness using Cyclic Analysis method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106526679B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108957555A (en) * | 2017-05-19 | 2018-12-07 | 中国石油化工股份有限公司 | Sedimentation process planar periodic analysis method |
| CN110263763A (en) * | 2019-07-09 | 2019-09-20 | 兰州大学 | Method based on lithology midpoint Sine Interpolation model analysis High-resolution Cyclostratigraphy |
| CN113777657A (en) * | 2021-09-07 | 2021-12-10 | 中国石油大学(华东) | Slope area well-connected profile seismic facies characterization industrialization process based on differential compaction recovery |
| CN113970796A (en) * | 2020-07-23 | 2022-01-25 | 中国石油化工股份有限公司 | Method for accurately recovering ancient water depth of sedimentary basin |
| CN113985476A (en) * | 2020-07-27 | 2022-01-28 | 中国石油天然气股份有限公司 | Stratum ancient structure recovery method, device and equipment and readable storage medium |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103424773A (en) * | 2012-05-25 | 2013-12-04 | 中国石油化工股份有限公司 | S Stratification position flattening method-based ancient landform restoration method |
| WO2014056541A1 (en) * | 2012-10-11 | 2014-04-17 | Omnisens S.A. | Method for monitoring soil erosion around buried devices, instrumented device and system implementing the 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 |
| CN105093311A (en) * | 2015-06-29 | 2015-11-25 | 成都理工大学 | Stratum denudation thickness measurement method for superposition basin multi-phase difference upheaval region |
| CN105204069A (en) * | 2014-06-27 | 2015-12-30 | 中国石油化工股份有限公司 | Strata denudation amount recovery method |
| CN105242328A (en) * | 2015-10-28 | 2016-01-13 | 中国石油大学(北京) | Paleo-thermal lithosphere thickness determination method and device |
| CN106291695A (en) * | 2015-05-30 | 2017-01-04 | 中国石油化工股份有限公司 | Eroded Thickness computational methods |
-
2017
- 2017-01-10 CN CN201710015647.5A patent/CN106526679B/en not_active Expired - Fee Related
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103424773A (en) * | 2012-05-25 | 2013-12-04 | 中国石油化工股份有限公司 | S Stratification position flattening method-based ancient landform restoration method |
| WO2014056541A1 (en) * | 2012-10-11 | 2014-04-17 | Omnisens S.A. | Method for monitoring soil erosion around buried devices, instrumented device and system implementing the method |
| 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 |
| CN106291695A (en) * | 2015-05-30 | 2017-01-04 | 中国石油化工股份有限公司 | Eroded Thickness computational methods |
| CN105093311A (en) * | 2015-06-29 | 2015-11-25 | 成都理工大学 | Stratum denudation thickness measurement method for superposition basin multi-phase difference upheaval region |
| CN105044797A (en) * | 2015-08-06 | 2015-11-11 | 中国石油天然气股份有限公司 | Quantitative recovery method for carbonate rock stratum denudation amount |
| CN105242328A (en) * | 2015-10-28 | 2016-01-13 | 中国石油大学(北京) | Paleo-thermal lithosphere thickness determination method and device |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108957555A (en) * | 2017-05-19 | 2018-12-07 | 中国石油化工股份有限公司 | Sedimentation process planar periodic analysis method |
| CN108957555B (en) * | 2017-05-19 | 2019-11-08 | 中国石油化工股份有限公司 | Sedimentation process planar periodic analysis method |
| CN110263763A (en) * | 2019-07-09 | 2019-09-20 | 兰州大学 | Method based on lithology midpoint Sine Interpolation model analysis High-resolution Cyclostratigraphy |
| CN110263763B (en) * | 2019-07-09 | 2021-06-29 | 兰州大学 | Method for analyzing gyrus stratum based on lithology midpoint sine interpolation model |
| CN113970796A (en) * | 2020-07-23 | 2022-01-25 | 中国石油化工股份有限公司 | Method for accurately recovering ancient water depth of sedimentary basin |
| CN113970796B (en) * | 2020-07-23 | 2023-06-23 | 中国石油化工股份有限公司 | Method for accurately recovering paleo-water depth of sedimentary basin |
| CN113985476A (en) * | 2020-07-27 | 2022-01-28 | 中国石油天然气股份有限公司 | Stratum ancient structure recovery method, device and equipment and readable storage medium |
| CN113985476B (en) * | 2020-07-27 | 2024-04-30 | 中国石油天然气股份有限公司 | Stratum ancient structure recovery method, device, equipment and readable storage medium |
| CN113777657A (en) * | 2021-09-07 | 2021-12-10 | 中国石油大学(华东) | Slope area well-connected profile seismic facies characterization industrialization process based on differential compaction recovery |
| CN113777657B (en) * | 2021-09-07 | 2023-10-17 | 中国石油大学(华东) | Slope zone well-connected profile seismic phase characterization industrialization flow based on differential compaction recovery |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106526679B (en) | 2018-05-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106526679B (en) | A kind of method for recovering unconformity erosion sediment thickness using Cyclic Analysis method | |
| CN101251604B (en) | Method for analyzing and NMO correcting two parameters transformation wave speed | |
| CN103454685B (en) | Method and device for predicting sand body thickness by utilizing logging constrained wave impedance inversion | |
| Chen et al. | Efficient ensemble-based closed-loop production optimization | |
| CN104879103B (en) | Layered water injection effect analysis method | |
| US8271243B2 (en) | System and method of integrating subterranean computer models for oil and gas exploration | |
| CN108802812A (en) | Well-seismic fusion stratum lithology inversion method | |
| CN107329171A (en) | Depth domain reservoir stratum seismic inversion method and device | |
| CN103857876A (en) | System and method for performing wellbore fracture operations | |
| CN104459782B (en) | Seismic velocity modeling method and modeling unit using thin layer chromatography inversion | |
| Liu et al. | Integration of time-lapse seismic data using the onset time approach: The impact of seismic survey frequency | |
| EP3033482A2 (en) | System and method of determining and optimizing waterflood performance | |
| CN103592680B (en) | A kind of log data based on FORWARD AND INVERSE PROBLEMS and Depth Domain seismic profile synthetic method | |
| CN105388518A (en) | Centroid frequency and spectral ratio integrated borehole seismic quality factor inversion method | |
| CN107907910A (en) | A kind of different lithology oil reservoir s-wave logging determines method | |
| CN105089615A (en) | Log data historical retrogression treatment method based on oil reservoir model | |
| Patterson et al. | Integration of microseismic data, completion data, and production data to characterize fracture geometry in the Permian Basin | |
| Alhuthali et al. | Optimal rate control under geologic uncertainty | |
| Iverson | Closure stress calculations in anisotropic formations | |
| CN103852789A (en) | Nonlinear chromatography method and device for seismic data | |
| Wallace et al. | Understanding Completion Performance in Niobrara-Codell Reservoirs Through the Use of Innovative Software-Guided Workflows and Models | |
| Sharma et al. | Optimizing polymerflood via rate control | |
| Gong et al. | Application of multi-level and high-resolution fracture modeling in field-scale reservoir simulation study | |
| Cheng et al. | Compressible streamlines and three-phase history matching | |
| Omatsone et al. | Redevelopment of the cardium formation using fractured horizontal wells: Reservoir engineering perspectives and early case histories |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | 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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180525 Termination date: 20220110 |