CN104089964B - The dating methods of method is analyzed based on well logging Milankovitch Cycle Hypothesis - Google Patents
The dating methods of method is analyzed based on well logging Milankovitch Cycle Hypothesis Download PDFInfo
- Publication number
- CN104089964B CN104089964B CN201410353207.7A CN201410353207A CN104089964B CN 104089964 B CN104089964 B CN 104089964B CN 201410353207 A CN201410353207 A CN 201410353207A CN 104089964 B CN104089964 B CN 104089964B
- Authority
- CN
- China
- Prior art keywords
- depth
- data
- deposit
- depth value
- record data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention relates to petroleum geology exploration field, in particular to the dating methods analyzing method based on well logging Milankovitch Cycle Hypothesis.The dating methods of method should be analyzed based on well logging Milankovitch Cycle Hypothesis, on the basis of traditional Milankovitch Cycle Hypothesis analyzes the dating methods of method, whether it is that event deposits by judging the deposit of interval of interest different depth, and then traditional Milankovitch Cycle Hypothesis will be used to analyze in the sedimentary record data of the different depth deposit obtained and the depth value of correspondence of method, corresponding data and corresponding depth value containing event deposition are rejected, and make the data after spending rejecting and depth value calculate the age of the deposit of designated depth.Thus avoid the impact that deposit carries out surveying year brought due to event deposition, make the age of deposit recorded the most accurate.
Description
Technical field
The present invention relates to petroleum geology exploration field, in particular to based on well logging Michaelis
The dating methods of gyration analysis method.
Background technology
Along with geological research has to the propelling in quantification direction, the determination to Ages of Sediments
Urgent demand, the determination of Ages of Sediments is to setting up basin chronostratigraphic architecture and carrying out
Oil-gas exploration is the most significant.The determination of Ages of Sediments mainly have isotope tracling method method,
Zircon surveys year method and traditional Milankovitch surveys year method (traditional Milankovitch Cycle Hypothesis survey year skill
Art).
Isotope tracling method method that the precision dating of ancient sediments is used and Zircon dating method,
Specifically carrying out surveying year when, testing expense is high, requires sample quality strict, also
It is exactly that two kinds of methods are all limited by sample and cost, it is impossible to effectively use on a large scale.Another
Plant the method that traditional reservoir rock and caprock means determine year.Cardinal principle is astrodisasters science
The cyclically-varying of parameter makes the amount of radiation also generating period change that earth's surface accepts.Ancient times
Amount of radiation change can be deposited thing and recorded, its information can be embodied in deposit magnetization
In the parameters such as rate, remanence, natural gamma.So using traditional Milankovitch to survey
During year method, as a example by calculating according to natural gamma data, first have to obtain the mesh of certain well
Interval natural gamma (GR) log data (the natural gamma data of the degree of depth and correspondence), then
The GR log data obtained is carried out spectrum analysis, (deep to obtain the GR curve of correspondence
Degree and natural gamma data institute constituent curve), and found out by the frequency curve of GR data and
Consult the frequency that the astrodisasters science obtained is corresponding, to frequency corresponding with astrodisasters science in GR data
The data filtering of rate processes, more filtered GR data and curves is matched astrodisasters science
On theoretical curve, (each point of astrodisasters science has known precise time, by money
Expect to consult the age of the boundary obtaining stratum and the degree of depth is also known, at stratigraphic boundary
Time value match with astrodisasters science, further according to the degree of depth in GR data and stratigraphic boundary
The degree of depth at place matches, and then filtered GR data and curves is matched astrodisasters science
On theoretical curve).Again the astrodisasters science of GR data and curves with ancient times is compared, based on
The definitiveness of astrodisasters science, just can get the age of respective depth deposit, as shown in Figure 1
Astrodisasters science mate schematic diagram with ancient sediments respective cycle filter result, uninterruptedly bent
Line represents astrodisasters science a certain frequency filtering curve, and unbroken curve represents that sedimentary record is with astronomical
The filter curve of cycle corresponding frequency, arrow represents mating of astrodisasters science and sedimentary record
Result.
But, for traditional Milankovitch Cycle Hypothesis dating technology, if deposit existing sudden
Event and the deposit that formed, then result of spectrum analysis and filter result will be caused relatively
Big error.And then result in deposit (deposit includes event deposit) survey
Year causes certain error.
Summary of the invention
It is an object of the invention to provide the survey year side analyzing method based on well logging Milankovitch Cycle Hypothesis
Method, to solve above-mentioned problem.
Provide the survey year analyzing method based on well logging Milankovitch Cycle Hypothesis in an embodiment of the present invention
Method, including:
Obtain sedimentary record data and the degree of depth of correspondence of the different depth deposit of interval of interest
Value;
Obtain core sample and the depth value of the different depth deposit of described interval of interest, and
Described core sample is analyzed, to determine the degree of depth at event deposit place;
The sedimentary record of the event deposit place degree of depth is rejected in described sedimentary record data
Data, treat spectrum analysis data and described depth value to obtain;
Treat that spectrum analysis data and described depth value enter with the astrodisasters science obtained in advance by described
Row coupling, to obtain the time produced by deposit.
Preferably, described it is analyzed including to described core sample: use core observation method
And/or described core sample is analyzed by sizing assay test method.
Preferably, in described acquisition after spectrum analysis data, treat frequency spectrum described by described
Before analytical data and described depth value mate with the astrodisasters science obtained in advance, also include:
Treat that spectrum analysis data and described depth value carry out compacting and recover to described, and use warp
Overcompaction treats spectrum analysis data and described depth value and the astronomy obtained in advance described in recovering
Cycle mates.
Preferably, described treat that spectrum analysis data are compacted extensive with described depth value to described
Include again:
Depth value according to different depth and the spacing of the adjacent degree of depth, after calculating compacting recovery
The sedimentary record data of different depth;
The sedimentary record data of the different depth after recovering compacting carry out resampling, to obtain
Spectrum analysis data are treated described in recovering through overcompaction.
Preferably, the described depth value according to different depth and the spacing of the adjacent degree of depth, calculate
The sedimentary record data of the different depth after compacting recovery include:
According to the depth value of different depth, calculate the compacting factor of different depth deposit;
According to described compacting factor and the spacing of the adjacent degree of depth that obtains in advance, calculate compacting extensive
The sedimentary record data of the different depth after Fu.
Preferably, the formula of the compacting factor of described calculating different depth deposit is: k=b-
A*lnH, wherein k is compacting factor, a and b is constant, and by the deposition of designated depth
Voidage, pore structure and rock composition determined, H is the degree of depth.
Preferably, described to described until spectrum analysis data carry out compacting recover after, in institute
State and treat that spectrum analysis data and described depth value are carried out with the astrodisasters science obtained in advance by described
Also include before coupling:
Whether the amplitude of variation of the sedimentary record data of contrast different depth scope is consistent, described
Depth bounds includes the multiple depth value of continuous print;
If it is not, then the sedimentary record data of different depth scope are normalized computing, with
Obtain the sedimentary record data after normalization computing and the depth value of correspondence;
The depth value of the sedimentary record data after described normalization computing and correspondence is carried out
Spacing resampling, to obtain sedimentary record data and the depth value of correspondence at even depth interval;
Use the depth value of sedimentary record data that described even depth is spaced and correspondence with in advance
The astrodisasters science obtained mates.
Preferably, described sedimentary record data include that deposit susceptibility data, deposit are surplus
Residual magnetism data, natural gamma data.
Preferably, described sizing assay test method includes:
Described core sample is carried out broken sample, remove organic, remove cement, pickling,
Grind and discrete particles processes, to obtain sample to be measured;
Laser particle analyzer is preheated 30 minutes, and described laser particle analyzer is carried out adjustment and survey
Amount background;
With Glass rod, the suspension of described sample to be measured is stirred;
The suspension of described sample to be measured is added in the sample cell of laser particle analyzer, be
The obscurity of the suspension of described sample to be measured reaches 10%-20%;
Suspension in described sample cell is carried out supersound process 5min;
Suspension in sample cell described in repetitive measurement, in the hope of the meansigma methods of measurement data;
Meansigma methods according to described measurement data judges whether described core sample is that event is sunk
Long-pending thing;
If, it is determined that the degree of depth at described core sample place is event deposit place
The degree of depth.
Preferably, described measurement data includes: percentage by volume corresponding to particle diameters at different levels and
Distribution curve, median particle diameter and the particle diameter corresponding to cumulative volume at different levels.
The survey year side of the method for analyzing based on well logging Milankovitch Cycle Hypothesis that the embodiment of the present invention provides
In method, with prior art, directly the interval of interest natural gamma (GR) obtaining certain well is surveyed
Well data (deposit depth and GR data) process, and will process after GR data
The curve formed is corresponding with astrodisasters science, thus obtains the age of the deposition of designated depth,
But cannot get rid of owing to event deposits the survey year deviation brought, it will be by getting
The core sample of interval of interest different depth be analyzed, to determine different depth
Whether the core sample of deposit has event deposit to produce, if it has, then corresponding
, reject the sedimentary record data (such as GR data) of the event deposit place degree of depth,
Treat the depth value of spectrum analysis data and correspondence to obtain, finally will treat spectrum analysis data and
Depth value mates with the astrodisasters science obtained in advance, to obtain the generation time of deposit,
Thus avoid the impact that deposition carries out surveying year brought due to event deposition.
Accompanying drawing explanation
Fig. 1 shows astrodisasters science and ancient sediments respective cycle filtering knot in conventional art
Fruit coupling schematic diagram;
Fig. 2 shows the survey year of the method for analyzing based on well logging Milankovitch Cycle Hypothesis of the embodiment of the present invention
Method basic flow sheet;
Fig. 3 a shows being affected and not by thing by event deposition of the embodiment of the present invention
The contrast oscillogram of part deposition impact;
Fig. 3 b shows being affected and not by thing by event deposition of the embodiment of the present invention
The spectrogram of part deposition impact;
Fig. 3 c shows being affected and not by thing by event deposition of the embodiment of the present invention
The filter result figure of part deposition impact;
Fig. 4 shows the tempestite grain size analysis figure of the embodiment of the present invention;
Fig. 5 shows the natural gamma data sectional schematic diagram of the embodiment of the present invention;
Fig. 6 shows through rejecting event deposition impact, compacting recovers and sedimentation rate is returned
One changes the spectrogram after correcting;
Fig. 7 shows natural gamma 125kyr filtering and the astrodisasters science of the embodiment of the present invention
125kyr filtering coupling schematic diagram;
Fig. 8 shows that the tempestite development position of the embodiment of the present invention determines schematic diagram;
Fig. 9 show the embodiment of the present invention through rejecting event deposition impact, be compacted extensive
Data spectrum comparison diagram again and before and after sedimentation rate normalization rectification.
Detailed description of the invention
Below by specific embodiment and combine accompanying drawing and the present invention is done further in detail
Describe.The embodiment of the present invention 1 provides the survey year side analyzing method based on well logging Milankovitch Cycle Hypothesis
The basic procedure of method, as in figure 2 it is shown, comprise the steps:
S201, obtains sedimentary record data and the correspondence of the different depth deposit of interval of interest
Depth value;
S202, obtains core sample and the depth value of the different depth deposit of interval of interest,
And core sample is analyzed, to determine the degree of depth at event deposit place;
S203, rejects the deposition note of the event deposit place degree of depth in sedimentary record data
Record data, treat spectrum analysis data and depth value to obtain;
S204, will treat that spectrum analysis data and depth value are carried out with the astrodisasters science obtained in advance
Coupling, to obtain the time produced by deposit.
Wherein, step S201, as described in the background art, measure the heavy of different depth
The long-pending time formed need first to obtain of both data, i.e. interval of interest depth value and
Sedimentary record data corresponding to this depth value.It should be noted that sedimentary record packet
Include: natural gamma data (GR data), deposit susceptibility data and deposit remanence
Data, these three data can form corresponding data and curves, and then obtain with depth value
Time produced by designated depth deposit.But the cost factor considered, it is thus achieved that GR number
According to difficulty and cost be less than acquisition deposit susceptibility data and deposit remanence data
Difficulty and cost, but precision does not has susceptibility and remanence data high.Thus, the present invention
The sedimentary record data used, it is preferred to use GR data are as sample, but consider not
The initial data that same region is obtained has different, therefore in indivedual areas, uses
Deposit susceptibility data and deposit remanence data carry out surveying the effect in year, compared to making
May be more preferable by the effect of GR data, thus when concrete operations, can be according to region
Difference, or obtain initial data difficulty and the difference of cost, select to use nature gal
In horse data (GR data), deposit susceptibility data and deposit remanence data one
Plant or multiple, carry out surveying year operation to the deposit of different depth.
Step S202, the deposit dating methods as described in background technology, if deposition
Record data include the data of event deposit, owing to event deposit is with general
The speed that property deposit is formed be have certain difference (owing to being affected by event, event
Property deposition speed generally than generality deposition speed faster, wherein event refer to as storm,
Tsunamis etc. make environment cause the natural phenomena of acute variation), in order to pick in sedimentary record data
Event of removing deposits the impact brought, and first has to determine which part data is that event is sunk
Volume data, namely first to obtain the interval of interest obtaining sedimentary record data in step S101
The core sample of different depth, and these core samples are analyzed, to determine thing
The degree of depth at part deposit place.Wherein, core sample is analyzed can using such as, adopts
By core observation method and/or sizing assay test method, both approaches is used to determine
Whether to core sample is event deposit (storm wind rock, turbidite etc.).Certainly, obtain
When taking the sedimentary record data of interval of interest different depth deposit, adjacent degree of depth spacing is more
Little, the sample sequence being namely made up of depth value is the most intensive, is more conducive to final calculating
The accuracy of result.Corresponding, the core sample obtaining different depth deposit also should be tried one's best
Ensure the dense degree in sampling interval.
Step S203, the event deposit place in determining sedimentary record data deep
After degree, by the sedimentary record data of the event deposition place degree of depth in sedimentary record data
Get rid of, just obtained treating spectrum analysis data and corresponding depth value.
Step S204, middle treat spectrum analysis data and described depth value by described and obtain in advance
Astrodisasters science mate, to obtain the time produced by deposit.Concrete, then
Treat spectrum analysis data and carry out spectrum analysis, (deep to obtain the spectrum analysis curve of correspondence
Degree and sedimentary record data institute constituent curve), and looked for by the frequency curve of sedimentary record data
Go out the frequency corresponding with consulting the astrodisasters science obtained, to all with astronomy in sedimentary record data
The data filtering of phase respective frequencies processes, more filtered sedimentary record data and curves is mated
On the theoretical curve of astrodisasters science, (each point of astrodisasters science has known accurate
Time, the age and the degree of depth that are obtained the boundary on stratum by data consultation are also known, logical
Cross the time value at stratigraphic boundary to match with astrodisasters science, further according in sedimentary record data
The degree of depth and stratigraphic boundary at the degree of depth match, and then by filtered sedimentary record data
Curve Matching is on the theoretical curve of astrodisasters science).Again by sedimentary record data and curves and ancient times
Astrodisasters science compare, definitiveness based on astrodisasters science, just can get respective depth
The age of deposit.The most just the operation that deposit is surveyed year is completed.
For traditional Milankovitch Cycle Hypothesis dating technology, if deposit exists paroxysmal thing
Part deposit, then result of spectrum analysis and filter result will be caused bigger error.
Assuming to exist the astrodisasters science of 125kyr, in the case of without event deposit, its sedimentary record is such as
Shown in the figure of Fig. 3 a upper part, if there is Storm Events effect at 155kyr, at storm
During growth, rapidly one section of thickness of deposition is equivalent to the heavy of 30kyr deposition under normal sedimentation speed
Long-pending thing, as shown in the figure of Fig. 3 a lower part.Respectively to without the sedimentary record of event deposit be subject to
The sedimentary record spectrum analysis of event deposit impact obtains under figure and Fig. 3 b of Fig. 3 b upper part
The figure of part, it can be seen that event deposition can make respective frequencies peak value in spectrum analysis diminish,
Frequency dispersion.To the sedimentary record shown in the figure of Fig. 3 a upper part and the figure of Fig. 3 a lower part
It is filtered acting on figure, the figure of Fig. 3 c lower part obtained such as Fig. 3 c upper part, in filtering
In result matching process, Storm Events fast deposition in short-term is mistaken as continue for 30kyr
Common deposition, therefore filter result is when second period mates, normal sedimentary record
Have compressed 80.6%, deposit is determined to cause certain error year by this.
The dating methods analyzing method based on well logging Milankovitch Cycle Hypothesis provided by the present invention, passes through
The core sample of the interval of interest different depth got is analyzed, different deep to determine
Whether the core sample of the deposit of degree has event deposit to produce, if it has, then corresponding
, reject the sedimentary record data (such as GR data) of the event deposit place degree of depth, with
Obtain and treat the depth value of spectrum analysis data and correspondence, finally will treat spectrum analysis data and deeply
Angle value is mated with the astrodisasters science obtained in advance, to obtain the generation time of deposit,
And avoid the impact that deposition carries out surveying year brought due to event deposition.
The embodiment of the present invention 2 provides the dating methods analyzing method based on well logging Milankovitch Cycle Hypothesis
Optimization method, on the basis of embodiment 1, concrete, described core sample is carried out
Analysis includes: use core observation method and/or sizing assay test method to enter described core sample
Row is analyzed.
It should be noted that different event deposits compares, general deposit is to have
Obvious characteristic, as a example by storm wind rock, core observation method is illustrated herein, but should
It is noted that the external appearance characteristic of various event deposits (such as storm wind rock, turbidite etc.)
It is well-known to those skilled in the art, or just can be directly obtained by inspection information.
Certainly, rock core described herein is the rock core after passing through process, it is possible to pass through naked eyes
Observe directly its external appearance characteristic.Below as a example by storm wind rock, core observation method is said
Bright.
Significant feature is had on the rock core of tempestite.The feature of tempestite is lower thick upper thin,
Rock core is massive structure.Tempestite bottom has washes away vestige, and granularity exists sudden change, middle part
Having cross-bedding, top is the fine-grained sediment of horizontal bedding.Storm mud stone is grown and cuts
Cut structure, mud stone tear bits, storm mud boulder, breeze layer, quasi-contemporaneous deformation structure, biology
The phenomenons such as escape trace, the mud stone in rock core there occurs that scrunch deforms, the end of to when this is storm generation
The drag interaction of portion's mud stone is formed, and is the typical characteristic of tempestite.
Sizing assay test method is that another kind can play and distinguishes that whether deposit is that event is sunk
The method of long-pending thing, sends out compared to direct core observation, and sizing assay test method is the most accurate,
But it is the most more, be not as direct as core observation method that expense and time spend.Simultaneously, it is contemplated that
When short grained rock core is distinguished, core observation method is used cannot effectively to play identification
Effect, thus, when being analyzed short grained rock core, it is preferred to use grain size analysis tries
Test method, to distinguish whether deposit is event deposition more accurately.
First sediment sample to be measured is carried out pretreatment by sizing assay test method, in order to sample
Product can carry out grain size analysis.Use laser particle analyzer that sample is tested afterwards, with
To granularity distribution result.Finally, the deposition be given by the software that laser particle analyzer is corresponding
Corresponding to composition granule percentage by volume and distribution curve, median particle diameter, cumulative volume at different levels
The information such as particle diameter judge whether deposit is event deposit.
Concrete, sizing assay test method comprises the steps:
Described core sample is carried out broken sample, remove organic, remove cement, pickling,
Grind and discrete particles processes, to obtain sample to be measured;
Laser particle analyzer is preheated 30 minutes, and described laser particle analyzer is carried out adjustment and survey
Amount background;
With Glass rod, the suspension of described sample to be measured is stirred;
The suspension of described sample to be measured is added in the sample cell of laser particle analyzer, be
The obscurity of the suspension of described sample to be measured reaches 10%-20%;
Suspension in described sample cell is carried out supersound process 5min;
Suspension in sample cell described in repetitive measurement, in the hope of the meansigma methods of measurement data;
Meansigma methods according to described measurement data judges whether described core sample is that event is sunk
Long-pending thing;
If, it is determined that the degree of depth at described core sample place is event deposit place
The degree of depth.
During core sample is prepared as sample to be measured, broken sample and reagent adding process
The process of sample is all carried out in the containers such as beaker, is carried out container, does not results in sample
The loss of product granule.During solid-liquid separation, under the high speed rotating of centrifuge, again can be by superfine
Little sample particle is separated with liquid.So, through the preparation process of above-mentioned measuring samples,
Sample to be measured grain fraction completely can be preserved in theory, reach to treat measuring samples complete
The grain size analysis of component particles.
Wherein, the suspension in sample cell described in repetitive measurement, putting down in the hope of measurement data
During average, generally circulation is measured and is advisable for 3-4 time, and each measurement set of time is 8s.And
And instrument automatically analyzes, detection data are also averaged, and draw granularity distribution result.
Further, as used Malvern 2000 laser particle analyzer, the software that can provide by it
Provide the percentage by volume corresponding to particle diameter at different levels and distribution curve, median particle diameter, at different levels tired
The information such as long-pending particle diameter corresponding to volume.Processed by the test data that this is tested and obtain figure
4.In core sample, suspending components content reaches 90%, and this meets tempestite in granularity
Typical characteristic, then just may determine that rock sample to be measured is tempestite.
Preferably, after step S203, before step S204, also include: to institute
State and treat that spectrum analysis data and described depth value carry out compacting and recover, and use through overcompaction extensive
Treat described in multiple that spectrum analysis data and described depth value are carried out with the astrodisasters science obtained in advance
Coupling.
In addition to deposit being surveyed the impact caused in year due to the data of event deposit, ground
Deposit can be surveyed and impact in year by layer compaction equally.Therefore, in order to overcome this impact,
Should namely obtain rejecting the deposition note of event deposit impact after step S203
After record data (treating spectrum analysis data), spectrum analysis data to be treated and correspondence
Depth value carries out compacting and recovers, the shadow brought in year to eliminate formation compaction effect to survey deposit
Ring.
Concrete, treat that spectrum analysis data and described depth value carry out compacting and recover bag to described
Include:
Depth value according to different depth and the spacing of the adjacent degree of depth, after calculating compacting recovery
The sedimentary record data of different depth;
The sedimentary record data of the different depth after recovering compacting carry out resampling, to obtain
Spectrum analysis data are treated described in recovering through overcompaction.
Further, according to depth value and the spacing of the adjacent degree of depth of different depth, calculate compacting
The sedimentary record data of the different depth after recovery include:
According to the depth value of different depth, calculate the compacting factor of different depth deposit;
According to described compacting factor and the spacing of the adjacent degree of depth that obtains in advance, calculate compacting extensive
The sedimentary record data of the different depth after Fu.
Further, the formula of the compacting factor calculating different depth deposit is: k=b-a*lnH,
Wherein k is compacting factor, a and b is constant, and the space of the deposition by designated depth
Degree, pore structure and rock composition are determined, H is the degree of depth.If it should be noted that
Rock is mud shale, then a can be with value 0.1, and b can be with value 1.46, if rock is sandstone,
Then a can be with value 0.08, and b can be with value 1.37.The value of a with b can be by consulting
Relevant data obtains.
Wherein, the sampling interval (spacing of the adjacent degree of depth) of consecutive points is divided by compacting factor just
Obtain the primary deposit thickness between the adjacent degree of depth, then by the original thickness superposition of different depth
Get up to obtain a new degree of depth sequence, the new degree of depth sequence being made up of multiple depth values with
The degree of depth sequence same sequence number being before made up of multiple depth values (arranges from according to depth size
The sequence that formed of the depth value of multiple degree of depth, corresponding to each depth value in the sequence
Numbering) corresponding sedimentary record data (natural gamma value) do not become.
New degree of depth sequence is carried out with corresponding natural gamma interpolation, generates even depth interval
Sample sequence, namely recover through deposit compacting, and equidistantly after resampling
The sequence that formed with the corresponding degree of depth of sedimentary record data.
By obtain sedimentary record data for illustrating as a example by natural gamma numerical value, sand shale
Having significant difference in natural gamma value, sandstone natural gamma numerical value is less, and mud stone is certainly
So gamma numerical value is bigger.By the contrast reference to natural gamma numerical value Yu core data, choosing
Select the suitable natural gamma value demarcation line as sand shale, as taken average 90 as boundary
Value, for the mud shale that is considered as more than 90, uses k=1.46-0.1*lnH formula to calculate pressure
Real coefficient, for the sandstone that is considered as less than 90, uses k=1.37-0.08*lnH formula to count
Calculate compacting factor k.
Can impacting in survey year deposit except event deposition and formation compaction effect
Outward, due to non-at the uniform velocity deposit the bottom age caused cannot accurately measure also be one important
Reason.In order to overcome this problem, to treat that spectrum analysis data are compacted extensive to described described
After Fu, treat spectrum analysis data and described depth value and the sky obtained in advance described by described
The literary composition cycle also includes before mating:
Whether the amplitude of variation of the sedimentary record data of contrast different depth scope is consistent, described
Depth bounds includes the multiple depth value of continuous print;
If it is not, then the sedimentary record data of different depth scope are normalized computing, with
Obtain the sedimentary record data after normalization computing and the depth value of correspondence;
The depth value of the sedimentary record data after described normalization computing and correspondence is carried out
Spacing resampling, to obtain sedimentary record data and the depth value of correspondence at even depth interval;
Use the depth value of sedimentary record data that described even depth is spaced and correspondence with in advance
The astrodisasters science obtained mates.
Concrete, as shown in Figure 5, it is possible to significantly observe, gamma ray curve is from flat
The angle of average size can be divided into three sections (as corresponding to the identification section of figure middle and upper part, to divide
For three sections, left, center, right), the abscissa wherein changed plan is depth value, and vertical coordinate is GR data
Numerical value.So, just may determine that the change width of the sedimentary record data of different depth scope
Degree does not corresponds, the most uneven, now tackles this kind of data and carries out sedimentation rate
Normalization computing.Certainly, in this application, sampled point is the most intensive, and the data obtained are the most just
It is the most accurate, then, the analysis result observed by data is the most accurate.
Sedimentary record data (the natural gamma number of different depth scope is can confirm that by Fig. 5
According to) amplitude of variation do not correspond.Now, then the deposition of reply different depth scope is remembered
Record data are normalized computing, with obtain the sedimentary record data after normalization computing and
Corresponding depth value.Concrete, natural gamma is worth size and sedimentation rate to be inverse relation,
Use equation below v=(150-ave)/100 that three sections of sedimentation rates are normalized computing, root
According to meansigma methods obtain every section corresponding regarding sedimentation rate, the sampling interval of every period is heavy divided by regarding
Long-pending speed obtains the depth interval sequence after sedimentation rate homogenization, is added up by its depth interval
Obtaining new degree of depth sequence, wherein ave is the numerical value of natural gamma data.Certainly, if
Use other kinds of sedimentary record data (as deposit susceptibility data and deposit remain
Magnetic data) as calculating data, the ave correspondence in formula also should be adjusted to corresponding number
Value unit.
Afterwards, new degree of depth sequence the preceding paragraph obtained and sedimentary record data (natural gal
Horse data) carry out adopting of the even depth interval after equidistant resampling obtains sedimentation rate correction
Sample sequence.And use sample sequence that even depth is spaced (include sedimentary record data and
Corresponding depth value) mate with the astrodisasters science obtained in advance, namely to even depth
Sedimentary record data in the sample sequence at interval (are namely recovered and deposition speed through overcompaction
Spectrum analysis data are treated after rate normalization computing) carry out spectrum analysis, to obtain correspondence
Spectrum analysis curve (degree of depth and sedimentary record data institute constituent curve), and pass through sedimentary record
The frequency corresponding with consulting the astrodisasters science obtained is found out in the frequency curve of data, to deposition note
In record data, the data filtering with astrodisasters science respective frequencies processes, then by filtered deposition
Record data and curves matches (each point of astrodisasters science on the theoretical curve of astrodisasters science
There is known precise time, obtained the age of boundary on stratum by data consultation with deep
Degree is also known, is matched with astrodisasters science by the time value at stratigraphic boundary, then root
Match according to the degree of depth at the degree of depth in sedimentary record data and stratigraphic boundary, and then will filtering
After sedimentary record data and curves match on the theoretical curve of astrodisasters science).Again deposition is remembered
The astrodisasters science in record data and curves and ancient times is compared, definitiveness based on astrodisasters science,
Just can get the age of respective depth deposit.
Concrete, the curve that the sample sequence of equity depth interval is formed carries out spectrum analysis,
Result as shown in Figure 6, obtain frequency be 0.00859,0.00982,0.0175,0.0249,
Five main peaks (5 stains in figure) of 0.0476, the corresponding cycle is 116.41,101.83,
57.14,40.16,21.01, its proportionate relationship is 1,0.875,0.490,0.345,0.180
In the cycle, with theoretical cycle 125kyr, the ratio of 96kyr, 54kyr, 40kyr, 23kyr is closed
Being 1,0.768,0.432,0.32,0.184 has preferable corresponding relation.Therefore
Above-mentioned six cycles are astrodisasters science records in deposit.Demarcated by known stratum
Point age and astrodisasters science are demarcated, and i.e. can get natural gamma 125kyr filtering and astrodisasters science
125kyr filters matching result, such as Fig. 7, natural gamma 125kyr filtering and astrodisasters science
125kyr filters coupling, and natural gamma filtering has good matching with astrodisasters science filtering.
On this basis, natural gamma 125kyr filtering filter with astrodisasters science 125kyr
Matching result can obtain the time that event deposition (such as tempestite) is occurred.125kyr
The time of astrodisasters science crest corresponding with natural gamma filter curve is identical.Between crest
Time of corresponding sampled point can be obtained by 40kyr filter curve, method is as follows:
First determine tempestite development position natural gamma filter curve any two peak values it
Between, as determined the corresponding time of Fig. 8 middle and upper part branch A, some A is positioned at natural gamma 125kyr
Between filter curve second and the 3rd peak value, this two peak values correspondence astrodisasters science filter curve
Second and the 3rd peak value of (astrodisasters science 125kyr curve), by astrodisasters science 125kyr
It is corresponding that second time to peak of curve understands second peak value of Gr125kyr filter curve
Time, Gr125kyr filter curve and Gr40kyr filter curve are to synchronize, Gr125kyr
Second peak value of filter curve can correspond to obtain in Gr40kyr filter curve straight down
Point B, then A point is corresponded to lower Gr40kyr filter curve obtain C, estimation B Yu C point
Between corresponding 40kyr periodicity P, the A point correspondence time be: T+P*40.The most just can
Draw the time point that the event deposit of rejecting is occurred.
The impact that causes by rejecting event to deposit, eliminate the shadow that causes of bottom compaction
Before and after ringing and solving the non-impact at the uniform velocity depositing and causing, data spectrum contrasts, as shown in Figure 8,
It can be seen that
1, the data spectrum of the natural gamma data after process can identify above-mentioned five frequencies
Rate main peak, and initial data frequency spectrum is difficult to and the peak value of astrodisasters science respective frequencies.
2, five frequency peak identified compared with initial data respective frequencies peak value, its energy
Amount all improves more than 30%.
3, in the data spectrum after process, dispersion phenomenon is effectively suppressed, and frequency dispersion rate reduces
More than 50%.
Dating methods and the tradition analyzing method based on well logging Milankovitch Cycle Hypothesis provided by the present invention
Dating methods compare, it is possible to by survey year error reduce year the most up to a million in hundreds of thousands year.
The dating methods analyzing method based on well logging Milankovitch Cycle Hypothesis provided by the present invention, passes through
The core sample of the interval of interest different depth got is analyzed, different deep to determine
Whether the core sample of the deposit of degree has event deposit to produce, if it has, then corresponding
, reject the sedimentary record data (such as GR data) of the event deposit place degree of depth, with
Obtain the depth value treating spectrum analysis data and correspondence, will treat that spectrum analysis data are entered the most again
Row compacting recovers and sedimentation rate normalization is corrected, and finally will recover through overcompaction and deposit speed
What rate normalization was corrected treats that spectrum analysis data and depth value enter with the astrodisasters science obtained in advance
Row coupling, to obtain the generation time of deposit, further avoid owing to event deposits
The impact that deposition carries out surveying year brought.
The foregoing is only the preferred embodiments of the present invention, be not limited to the present invention,
For a person skilled in the art, the present invention can have various modifications and variations.All
Within the spirit and principles in the present invention, any modification, equivalent substitution and improvement etc. made,
Should be included within the scope of the present invention.
Claims (7)
1. analyze the dating methods of method based on well logging Milankovitch Cycle Hypothesis, it is characterised in that bag
Include:
Obtain sedimentary record data and the degree of depth of correspondence of the different depth deposit of interval of interest
Value;
Obtain core sample and the depth value of the different depth deposit of described interval of interest, and
Described core sample is analyzed, to determine the degree of depth at event deposit place;
The sedimentary record of the event deposit place degree of depth is rejected in described sedimentary record data
Data, treat spectrum analysis data and described depth value to obtain;
Treat that spectrum analysis data and described depth value enter with the astrodisasters science obtained in advance by described
Row coupling, to obtain the time produced by deposit;
Described it is analyzed including to described core sample: use core observation method and/or granularity
Analyze test method(s) described core sample is analyzed;
Described sedimentary record data include deposit susceptibility data, deposit remanence data,
Natural gamma data;
Described sizing assay test method includes:
Described core sample is carried out broken sample, remove organic, remove cement, pickling,
Grind and discrete particles processes, to obtain sample to be measured;
Laser particle analyzer is preheated 30 minutes, and described laser particle analyzer is carried out adjustment and survey
Amount background;
With Glass rod, the suspension of described sample to be measured is stirred;
The suspension of described sample to be measured is added in the sample cell of laser particle analyzer, be
The obscurity of the suspension of described sample to be measured reaches 10%-20%;
Suspension in described sample cell is carried out supersound process 5min;
Suspension in sample cell described in repetitive measurement, in the hope of the meansigma methods of measurement data;
Meansigma methods according to described measurement data judges whether described core sample is that event is sunk
Long-pending thing;
If, it is determined that the degree of depth at described core sample place is event deposit place
The degree of depth.
The survey year analyzing method based on well logging Milankovitch Cycle Hypothesis the most according to claim 1
Method, it is characterised in that in described acquisition after spectrum analysis data, described by described
Before treating that spectrum analysis data and described depth value mate with the astrodisasters science obtained in advance,
Also include:
Treat that spectrum analysis data and described depth value carry out compacting and recover to described, and use warp
Overcompaction treats spectrum analysis data and described depth value and the astronomy obtained in advance described in recovering
Cycle mates.
The survey year analyzing method based on well logging Milankovitch Cycle Hypothesis the most according to claim 2
Method, it is characterised in that described treat that spectrum analysis data are carried out with described depth value to described
Compacting recovery includes:
Depth value according to different depth and the spacing of the adjacent degree of depth, after calculating compacting recovery
The sedimentary record data of different depth;
The sedimentary record data of the different depth after recovering compacting carry out resampling, to obtain
Spectrum analysis data are treated described in recovering through overcompaction.
The survey year analyzing method based on well logging Milankovitch Cycle Hypothesis the most according to claim 3
Method, it is characterised in that the described depth value according to different depth and the spacing of the adjacent degree of depth,
The sedimentary record data calculating the different depth after compacting recovers include:
According to the depth value of different depth, calculate the compacting factor of different depth deposit;
According to described compacting factor and the spacing of the adjacent degree of depth that obtains in advance, calculate compacting extensive
The sedimentary record data of the different depth after Fu.
The survey year analyzing method based on well logging Milankovitch Cycle Hypothesis the most according to claim 4
Method, it is characterised in that the formula of the compacting factor of described calculating different depth deposit is:
K=b-a*lnH, wherein k is compacting factor, a and b is constant, and by designated depth
The voidage of deposition, pore structure and rock composition determined, H is the degree of depth.
The survey year analyzing method based on well logging Milankovitch Cycle Hypothesis the most according to claim 2
Method, it is characterised in that described to described until spectrum analysis data carry out compacting recover after,
Spectrum analysis data and described depth value and the astrodisasters science obtained in advance is treated by described described
Also include before mating:
Whether the amplitude of variation of the sedimentary record data of contrast different depth scope is consistent, described
Depth bounds includes the multiple depth value of continuous print;
If it is not, then the sedimentary record data of different depth scope are normalized computing, with
Obtain the sedimentary record data after normalization computing and the depth value of correspondence;
The depth value of the sedimentary record data after described normalization computing and correspondence is carried out
Spacing resampling, to obtain sedimentary record data and the depth value of correspondence at even depth interval;
Use the depth value of sedimentary record data that described even depth is spaced and correspondence with in advance
The astrodisasters science obtained mates.
The survey year analyzing method based on well logging Milankovitch Cycle Hypothesis the most according to claim 1
Method, it is characterised in that described measurement data includes: the volume hundred corresponding to particle diameters at different levels
Mark and distribution curve, median particle diameter and the particle diameter corresponding to cumulative volume at different levels.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410353207.7A CN104089964B (en) | 2014-07-23 | 2014-07-23 | The dating methods of method is analyzed based on well logging Milankovitch Cycle Hypothesis |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410353207.7A CN104089964B (en) | 2014-07-23 | 2014-07-23 | The dating methods of method is analyzed based on well logging Milankovitch Cycle Hypothesis |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104089964A CN104089964A (en) | 2014-10-08 |
CN104089964B true CN104089964B (en) | 2016-08-24 |
Family
ID=51637700
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410353207.7A Active CN104089964B (en) | 2014-07-23 | 2014-07-23 | The dating methods of method is analyzed based on well logging Milankovitch Cycle Hypothesis |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104089964B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107884830A (en) * | 2016-09-29 | 2018-04-06 | 中国石油化工股份有限公司 | Sedimentary formation cycle partitioning method and device |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108427146B (en) * | 2017-02-13 | 2019-08-30 | 中国石油化工股份有限公司 | A kind of method and system of Milankovitch Cycles in identification sedimentary formation |
CN109031443B (en) * | 2017-06-09 | 2021-05-14 | 中国石油化工股份有限公司 | Method for recovering denudation amount based on fixed years of Milnaciaceae Virginia gyrus |
CN108919369B (en) * | 2018-05-16 | 2019-10-11 | 中国海洋石油集团有限公司 | A kind of method of quantitative judge Milankovitch Cycles |
CN109001840B (en) * | 2018-07-18 | 2019-08-09 | 中国地质大学(北京) | A kind of mud shale genetic analysis method and device |
CN110231660B (en) * | 2019-05-15 | 2021-04-27 | 防灾科技学院 | Method for determining dislocation quantity of blind active fault by using susceptibility logging method |
CN110454155B (en) * | 2019-08-23 | 2023-04-11 | 防灾科技学院 | Method for determining age of quaternary stratum by using susceptibility logging method |
CN112950551A (en) * | 2021-02-05 | 2021-06-11 | 自然资源部第一海洋研究所 | Ocean sediment rock core age rapid calibration method |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
IE960881A1 (en) * | 1996-12-13 | 1998-07-01 | Ronan Francis O Doherty | Method of distinguishing geological sequences and their¹boundaries |
-
2014
- 2014-07-23 CN CN201410353207.7A patent/CN104089964B/en active Active
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107884830A (en) * | 2016-09-29 | 2018-04-06 | 中国石油化工股份有限公司 | Sedimentary formation cycle partitioning method and device |
CN107884830B (en) * | 2016-09-29 | 2020-03-24 | 中国石油化工股份有限公司 | Sedimentary formation cycle dividing method and device |
Also Published As
Publication number | Publication date |
---|---|
CN104089964A (en) | 2014-10-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104089964B (en) | The dating methods of method is analyzed based on well logging Milankovitch Cycle Hypothesis | |
Boggs | Petrology of sedimentary rocks | |
CN103122762B (en) | Detection method and device for effective fractured intervals in unconventional shale oil and gas reservoir | |
Davies et al. | Sandstone composition and depositional environment | |
CN105675635B (en) | Tight rock component relative amount and brittleness index determine method and apparatus | |
Yarmohammadi et al. | An integrated approach for heterogeneity analysis of carbonate reservoirs by using image log based porosity distributions, NMR T2 curves, velocity deviation log and petrographic studies: A case study from the South Pars gas field, Persian Gulf Basin | |
CN105044797B (en) | A kind of carbonate formation denudation quantitative reconstruction method | |
CN107748399B (en) | Method for identifying deep tectonic layer of mountain front zone by utilizing gravity interface inversion | |
CN106285660A (en) | A kind of multilamellar sandstone oil reservoir Low Resistivity Reservoir Identification Methods and device | |
CN103513286B (en) | Dam, beach construction unit method of discrimination under a kind of geological model constraint | |
CN103744109B (en) | Method for identifying a weathering crust structure of clastic rock in a area of covering no well | |
CN107688206A (en) | A kind of HIGH-RESOLUTION SEQUENCE STRATIGRAPHIC DIVISION and control methods | |
CN107609253B (en) | Carbonate rock deposition numerical simulation method | |
Cosgrove et al. | Filter or conveyor? Establishing relationships between clinoform rollover trajectory, sedimentary process regime, and grain character within intrashelf clinothems, offshore New Jersey, USA | |
Zhang et al. | Reservoir quality variations within a sinuous deep water channel system in the Niger Delta Basin, offshore West Africa | |
CN105301647B (en) | The method for distinguishing grey matter mud stone and sandstone | |
Burton et al. | Quantitative shale characterization of the tidally influenced Sego Sandstone | |
CN107884831B (en) | A kind of AVO type identification method based on proximal support vector machines | |
CN103412335B (en) | A kind of method utilizing earthquake thing phase body predicting reservoir | |
CN107817260A (en) | A kind of mud shale high frequency sequence recognition methods | |
CN111580183B (en) | Method for carrying out deep quantitative reduction on ancient lake water | |
Gartmair et al. | One size does not fit all: Refining zircon provenance interpretations via integrated grain shape, geochronology, and Hf isotope analysis | |
CN105590018A (en) | Oil-water layer identification method for sandstone and mudstone thin interbed oil reservoir | |
Fitch | Heterogeneity in the petrophysical properties of carbonate reservoirs | |
Jackson et al. | Gamma-Ray Logs and U–Pb Zircon Geochronology—Essential Tools to Constrain Lithofacies Interpretation of Paleoproterozoic Depositional Systems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |