CN104089964A - Dating method based on logging Milankovitch cycle analysis method - Google Patents
Dating method based on logging Milankovitch cycle analysis method Download PDFInfo
- Publication number
- CN104089964A CN104089964A CN201410353207.7A CN201410353207A CN104089964A CN 104089964 A CN104089964 A CN 104089964A CN 201410353207 A CN201410353207 A CN 201410353207A CN 104089964 A CN104089964 A CN 104089964A
- Authority
- CN
- China
- Prior art keywords
- depth
- data
- record data
- sediment
- sample
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention relates to the field of petroleum geological exploration, and particularly relates to a dating method based on a logging Milankovitch cycle analysis method. According to the dating method based on a logging Milankovitch cycle analysis method, based on a conventional dating method of a Milankovitch cycle analysis method, by judging whether sediments in different depths in an interval of interest is eventfully deposited or not, further, corresponding data and corresponding depth values containing eventful deposition are eliminated from deposition recording data and corresponding depth values of sediments in different depths, wherein the sediments are obtained by using the conventional Milankovitch cycle analysis method, and ages of the sediment in the appointed depth is calculated by using the data and the depth value obtained after the corresponding data and corresponding depth values containing eventful deposition are eliminated. Therefore, the influence on dating of the sediments due to eventful deposition is prevented so that the dated age of the sediment is accurate.
Description
Technical field
The present invention relates to petroleum geology exploration field, in particular to the survey year method based on well logging Milankovitch Cycle Hypothesis analytical approach.
Background technology
Along with the propelling of geological research to quantification direction, Ages of Sediments definite had to urgent demand, determining when setting up basin etc. Stratigraphic framework and carrying out oil-gas exploration of Ages of Sediments is all significant.Definite isotope that mainly contains of Ages of Sediments is surveyed a year method, zircon survey year method and a traditional Milankovitch survey year method (traditional Milankovitch Cycle Hypothesis is surveyed a year technology).
The isotope that the precision dating of ancient sediments is used is surveyed year method and Zircon dating method, when specifically surveying year, testing expense is high, sample quality is required strict, namely two kinds of methods are limit by sample and cost all, cannot effectively use on a large scale.Another kind of traditional Milankovitch cycle means are determined the method in year.Cardinal principle is the also generating period variation of radiant quantity that the cyclical variation of astronomical cycle parameter is accepted earth's surface.The radiant quantity variation in ancient times can be deposited thing and record, and its information can be embodied in the parameters such as sediment magnetic susceptibility, remanence, natural gamma.So when using traditional Milankovitch survey year method, take and calculate as example according to natural gamma data, first to obtain objective interval natural gamma (GR) log data (degree of depth and corresponding natural gamma data) of certain well, then the GR log data obtaining is carried out to spectrum analysis, to obtain corresponding GR curve (degree of depth and natural gamma data institute constituent curve), and find out and consult frequency corresponding to astronomical cycle of obtaining by the frequency curve of GR data, to processing with the data filtering of astronomical cycle respective frequencies in GR data, again filtered GR data and curves is matched on the theoretical curve in astronomical cycle, (each point in astronomical cycle has known precise time, the age and the degree of depth that by data consultation, obtain the boundary on stratum are also known, time value by stratigraphic boundary place matched with the astronomical cycle, according to the degree of depth at the degree of depth in GR data and stratigraphic boundary place, match again, and then filtered GR data and curves is matched on the theoretical curve in astronomical cycle).Again the astronomical cycle in GR data and curves and ancient times is compared, determinacy based on the astronomical cycle, just can obtain the sedimental age of respective depth, the astronomical cycle is as shown in Figure 1 mated schematic diagram with ancient sediments respective cycle filtering result, uninterrupted curve represents astronomical a certain frequency filtering curve of cycle, unbroken curve represents the filter curve of sedimentary record and astronomical corresponding frequency of cycle, and arrow represents the matching result of astronomical cycle and sedimentary record.
But, for traditional Milankovitch Cycle Hypothesis, survey a year technology, the sediment forming if there are unexpected incidents in sediment, will cause larger error to result of spectrum analysis and filtering result so.And then caused sediment (sediment comprises event sediment) survey year to cause certain error.
Summary of the invention
The object of the present invention is to provide the survey year method based on well logging Milankovitch Cycle Hypothesis analytical approach, to solve the above problems.
Survey year method based on well logging Milankovitch Cycle Hypothesis analytical approach is provided in an embodiment of the present invention, has comprised:
Obtain the sedimental sedimentary record data of different depth and the corresponding depth value of objective interval;
Obtain the sedimental core sample of different depth and the depth value of described objective interval, and described core sample is analyzed, to determine the degree of depth at event sediment place;
In described sedimentary record data, reject the sedimentary record data of the event sediment place degree of depth, to obtain, treat spectrum analysis data and described depth value;
By described, treat that spectrum analysis data and described depth value mated with the astronomical cycle of obtaining in advance, the time being produced to obtain sediment.
Preferably, described described core sample analysis is comprised: adopt core observation method and/or sizing assay test method to carry out described core sample.
Preferably, in described acquisition after spectrum analysis data, described by described treat that spectrum analysis data and described depth value mated with the astronomical cycle of obtaining in advance before, also comprise:
To described, treat that spectrum analysis data and described depth value carry out compacting recovery, and use described in overcompaction recovery and treat that spectrum analysis data and described depth value mated with the astronomical cycle of obtaining in advance.
Preferably, describedly to described, treat that spectrum analysis data and described depth value carry out compacting and recover to comprise:
According to the spacing of the depth value of different depth and the adjacent degree of depth, calculate the sedimentary record data of the different depth after compacting recovers;
The sedimentary record data of the different depth after compacting is recovered resample, and to obtain described in overcompaction recovery, treat spectrum analysis data.
Preferably, described according to the spacing of the depth value of different depth and the adjacent degree of depth, the sedimentary record data of calculating the different depth after compacting recovers comprise:
According to the depth value of different depth, calculate the sedimental compaction coefficient of different depth;
According to the spacing of described compaction coefficient and the adjacent degree of depth obtained in advance, calculate the sedimentary record data of the different depth after compacting recovers.
Preferably, the formula of the sedimental compaction coefficient of described calculating different depth is: k=b-a*lnH, and wherein k is compaction coefficient, and a and b are constant, and determine by voidage, pore texture and the rock composition of the deposition of designated depth, and H is the degree of depth.
Preferably, described to described after spectrum analysis data are carried out compacting recovery, described, described treating also comprised before spectrum analysis data and described depth value mated with the astronomical cycle of obtaining in advance:
Whether the amplitude of variation of the sedimentary record data of contrast different depth scope conforms to, and described depth range comprises continuous a plurality of depth values;
If not, the sedimentary record data of different depth scope are normalized to computing, to obtain sedimentary record data and the corresponding depth value after normalization computing;
Sedimentary record data after described normalization computing and corresponding depth value are equidistantly resampled, to obtain sedimentary record data and the corresponding depth value at even depth interval;
Use sedimentary record data and the corresponding depth value at described even depth interval to mate with the astronomical cycle of obtaining in advance.
Preferably, described sedimentary record data comprise sediment magnetic susceptibility data, sediment remanence data, natural gamma data.
Preferably, described sizing assay test method comprises:
To described core sample carry out broken sample, remove organic matter, remove cementing matter, pickling, grinding and discrete particles process, to obtain sample to be measured;
By laser particle analyzer preheating 30 minutes, and described laser particle analyzer is carried out adjustment and measures background;
With glass bar, the suspension of described sample to be measured is stirred;
The suspension of described sample to be measured being added in the sample cell of laser particle analyzer, is that the obscurity of the suspension of described sample to be measured reaches 10%-20%;
Suspension in described sample cell is carried out to ultrasonic processing 5min;
Repeatedly measure the suspension in described sample cell, in the hope of the mean value of measurement data;
According to the mean value of described measurement data, judge whether described core sample is event sediment;
If so, the degree of depth of determining described core sample place is the degree of depth at event sediment place.
Preferably, described measurement data comprises: the corresponding percent by volume of particle diameters at different levels and distribution curve, median particle diameter and the corresponding particle diameter of cumulative volume at different levels.
The survey year method based on well logging Milankovitch Cycle Hypothesis analytical approach that the embodiment of the present invention provides, directly objective interval natural gamma (GR) log data (the deposition degree of depth and GR data) that obtains certain well is processed with of the prior art, and the formed curve of GR data after processing is corresponding with the astronomical cycle, thereby obtain the age of the deposition of designated depth, but cannot get rid of the survey year deviation of bringing due to event deposition, it is by analyzing the core sample of the objective interval different depth getting, to determine whether the sedimental core sample of different depth has event sediment to produce, if had, corresponding, the sedimentary record data (as GR data) of the rejecting event sediment place degree of depth, to obtain, treat spectrum analysis data and corresponding depth value, finally will treat that spectrum analysis data and depth value mated with the astronomical cycle of obtaining in advance, to obtain sedimental generation time, thereby avoided, due to what event deposition was brought, deposition is surveyed to the impact in year.
Accompanying drawing explanation
Fig. 1 shows in conventional art the astronomical cycle and mates schematic diagram with ancient sediments respective cycle filtering result;
Fig. 2 shows the survey year method basic flow sheet based on well logging Milankovitch Cycle Hypothesis analytical approach of the embodiment of the present invention;
Fig. 3 a shows the event that the is subject to deposition affects and the contrast oscillogram that is not subject to event deposition affects of the embodiment of the present invention;
Fig. 3 b shows the event that the is subject to deposition affects and the spectrogram that is not subject to event deposition affects of the embodiment of the present invention;
Fig. 3 c shows the event that the is subject to deposition affects and the filtering result figure that is not subject to event deposition affects of the embodiment of the present invention;
Fig. 4 shows the tempestite sreen 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 the spectrogram after event deposition affects, compacting recovery and rate of sedimentation normalization are corrected;
The natural gamma 125kyr filtering that Fig. 7 shows the embodiment of the present invention is mated schematic diagram with astronomical cycle 125kyr filtering;
Fig. 8 shows the tempestite development position of the embodiment of the present invention and determines schematic diagram;
Fig. 9 shows the data spectrum comparison diagram before and after correcting through rejecting event deposition affects, compacting recovery and rate of sedimentation normalization of the embodiment of the present invention.
Embodiment
Below by specific embodiment, also by reference to the accompanying drawings the present invention is described in further detail.The embodiment of the present invention 1 provide based on well logging Milankovitch Cycle Hypothesis analytical approach survey year method basic procedure, as shown in Figure 2, comprise the steps:
S201, obtains the sedimental sedimentary record data of different depth of objective interval and corresponding depth value;
S202, obtains the sedimental core sample of different depth and the depth value of objective interval, and core sample is analyzed, to determine the degree of depth at event sediment place;
S203 rejects the sedimentary record data of the event sediment place degree of depth in sedimentary record data, to obtain, treats spectrum analysis data and depth value;
S204, will treat that spectrum analysis data and depth value mated with the astronomical cycle of obtaining in advance, the time being produced to obtain sediment.
Wherein, step S201, as described in the background art, the formed time of deposition of measuring different depth need to first obtain the data of two aspects, i.e. the depth value of objective interval and the corresponding sedimentary record data of this depth value.It should be noted that, sedimentary record data comprise: natural gamma data (GR data), sediment magnetic susceptibility data and sediment remanence data, these three kinds of data can both form corresponding data and curves with depth value, and then obtain the time that designated depth sediment produces.But the cost factor of considering, the difficulty and the cost that obtain GR data are less than difficulty and the cost that obtains sediment magnetic susceptibility data and sediment remanence data, but precision do not have magnetic susceptibility and remanence data high.Thus, sedimentary record data of the present invention, preferably adopt GR data as sample, but consider that the raw data that different regions obtains has different, therefore in indivedual areas, use sediment magnetic susceptibility data and sediment remanence data to survey the effect in year, effect compared to use GR data may be better, thus when concrete operations, can be according to the difference of region, or obtain the difference of raw data difficulty and cost, carry out choice for use natural gamma data (GR data), one or more in sediment magnetic susceptibility data and sediment remanence data, the sediment of different depth is surveyed to a year operation.
Step S202, as the sediment of being introduced in background technology is surveyed year method, if include the sedimental data of event in sedimentary record data, due to event sediment and the formed speed of general sediment be have certain difference (owing to being subject to events affecting, the speed of event deposition is conventionally faster than the speed of generality deposition, wherein event refers to as storm, tsunamis etc. make environment cause the spontaneous phenomenon of acute variation), the impact bringing in order to weed out event deposition in sedimentary record data, first to determine which partial data is event deposition data, namely to first obtain the core sample of the different depth of the objective interval that obtains sedimentary record data in step S101, and these core samples are analyzed, to determine the degree of depth at event sediment place.Wherein, core sample analysis can adopt as, adopt core observation method and/or sizing assay test method, whether what use these two kinds of methods to determine to obtain core sample is event sediment (storm wind rock, turbidite etc.).Certainly, while obtaining the sedimental sedimentary record data of objective interval different depth, adjacent degree of depth spacing is less, and the sample sequence being namely comprised of depth value is more intensive, is more conducive to the accuracy of final calculation result.Corresponding, obtain the sedimental core sample of different depth and also should guarantee the dense degree of sampling interval as far as possible.
Step S203 after the degree of depth at the event sediment place in having determined sedimentary record data, gets rid of the sedimentary record data of the event deposition place degree of depth in sedimentary record data, has just obtained treating spectrum analysis data and corresponding depth value.
Step S204, middlely treats that spectrum analysis data and described depth value mated with the astronomical cycle of obtaining in advance, the time being produced to obtain sediment by described.Concrete, then treat spectrum analysis data and carry out spectrum analysis, to obtain corresponding spectrum analysis curve (degree of depth and sedimentary record data institute constituent curve), and find out and consult frequency corresponding to astronomical cycle of obtaining by the frequency curve of sedimentary record data, to processing with the data filtering of astronomical cycle respective frequencies in sedimentary record data, again filtered sedimentary record data and curves is matched on the theoretical curve in astronomical cycle, (each point in astronomical cycle has known precise time, the age and the degree of depth that by data consultation, obtain the boundary on stratum are also known, time value by stratigraphic boundary place matched with the astronomical cycle, according to the degree of depth at the degree of depth in sedimentary record data and stratigraphic boundary place, match again, and then filtered sedimentary record data and curves is matched on the theoretical curve in astronomical cycle).The astronomical cycle in sedimentary record data and curves and ancient times is compared, the determinacy based on the astronomical cycle, just can obtain the sedimental age of respective depth again.Just completed thus the operation of surveying year to sedimental.
For traditional Milankovitch Cycle Hypothesis, survey a year technology, if there is paroxysmal event deposit thing in sediment, will cause larger error to result of spectrum analysis and filtering result so.Suppose to exist the astronomical cycle of 125kyr, without event deposit in the situation that, its sedimentary record is as shown in the figure of part on Fig. 3 a, if there is Storm Events effect at 155kyr place, at storm, between the puberty, deposit rapidly the sediment that one section of thickness is equivalent to 30kyr deposition under normal sedimentation speed, as shown in the figure of part under Fig. 3 a.To the sedimentary record without event deposit be subject to the sedimentary record spectrum analysis of event deposition affects to obtain under the figure of part Fig. 3 b on and Fig. 3 b figure partly, can see that outgoing event deposition can make respective frequencies peak value in spectrum analysis diminish, frequency dispersion respectively.Sedimentary record shown in the figure of part under the figure of part on Fig. 3 a and Fig. 3 a is carried out to filter action to be obtained as the figure of part under figure, Fig. 3 c of part on Fig. 3 c, in filtering result matching process, Storm Events in short-term fast deposition is mistaken as the common deposition that has continued 30kyr, therefore filtering result is when second period mates, normal sedimentary record has compressed 80.6%, and this has caused certain error to fixed year of sediment.
Survey year method based on well logging Milankovitch Cycle Hypothesis analytical approach provided by the present invention, by the core sample of the objective interval different depth getting is analyzed, to determine whether the sedimental core sample of different depth has event sediment to produce, if had, corresponding, the sedimentary record data (as GR data) of the rejecting event sediment place degree of depth, to obtain, treat spectrum analysis data and corresponding depth value, finally will treat that spectrum analysis data and depth value mated with the astronomical cycle of obtaining in advance, to obtain sedimental generation time, and avoided, due to what event deposition was brought, deposition is surveyed to the impact in year.
The embodiment of the present invention 2 provide based on well logging Milankovitch Cycle Hypothesis analytical approach survey year method optimization method, on the basis of embodiment 1, concrete, described core sample analysis is comprised: adopt core observation method and/or sizing assay test method to carry out described core sample.
It should be noted that, the different event sediments general sediment that compares has obvious characteristic, the storm wind rock of take herein describes core observation method as example, but it should be noted, the external appearance characteristic of variety of event sediment (as storm wind rock, turbidite etc.) is well-known to those skilled in the art, just or can directly obtain by inspection information.Certainly, rock core described herein is rock core after treatment, can observe directly its external appearance characteristic by naked eyes.Storm wind rock take below as example, core observation method is described.
On the rock core of tempestite, there is significant feature.The feature of tempestite is lower thick upper thin, and rock core is clumpy structure.Tempestite bottom has the vestige of washing away, and granularity exists sudden change, and middle part has diagonal stratification, and top is the fine-grained sediment of horizontal stratification.In storm mud stone, grow and cut off structure, mud stone is torn the phenomenons such as bits, storm mud boulder, breeze layer, accurate contemporaneous deformation structure, biological escape trace, there is to rub wrinkle distortion in the mud stone in rock core, when this is storm generation, the drag interaction of bottom mud stone is formed, and be the characteristic feature of tempestite.
Sizing assay test method is that another kind can play to distinguish whether sediment is the sedimental method of event, compared to direct core observation, sends out, and sizing assay test method is more accurate, but the cost of expense and time is also more, not as core observation method direct., consider when short grained rock core is distinguished meanwhile, adopt core observation method cannot effectively play identification effect, thus, when short grained rock core is analyzed, preferably adopt sizing assay test method, take and distinguish more accurately whether sediment is event deposition.
First sizing assay test method carries out pre-service by sediment sample to be measured, so that sample can carry out sreen analysis.Use afterwards laser particle analyzer to test sample, to obtain size-grade distribution result.The information such as the sedimentary particle percent by volume finally, providing by the corresponding software of laser particle analyzer and distribution curve, median particle diameter, the corresponding particle diameter of cumulative volumes at different levels judge whether sediment is event sediment.
Concrete, sizing assay test method comprises the steps:
To described core sample carry out broken sample, remove organic matter, remove cementing matter, pickling, grinding and discrete particles process, to obtain sample to be measured;
By laser particle analyzer preheating 30 minutes, and described laser particle analyzer is carried out adjustment and measures background;
With glass bar, the suspension of described sample to be measured is stirred;
The suspension of described sample to be measured being added in the sample cell of laser particle analyzer, is that the obscurity of the suspension of described sample to be measured reaches 10%-20%;
Suspension in described sample cell is carried out to ultrasonic processing 5min;
Repeatedly measure the suspension in described sample cell, in the hope of the mean value of measurement data;
According to the mean value of described measurement data, judge whether described core sample is event sediment;
If so, the degree of depth of determining described core sample place is the degree of depth at event sediment place.
Core sample is being prepared in the process of sample to be measured, the process of broken sample and reagent adding processing sample is all carried out in the containers such as beaker, and container is cleaned, and can not cause the loss of sample particle.During Separation of Solid and Liquid, under the High Rotation Speed of hydro-extractor, superfine little sample particle and fluid separation applications can be opened again.So, through the preparation process of above-mentioned measurement sample, can preserve in theory sample to be measured grain fraction completely, reach the sreen analysis to sample full constituent particle to be measured.
Wherein, repeatedly measure the suspension in described sample cell, during in the hope of the mean value of measurement data, circulation is measured and is advisable for 3-4 time conventionally, and each Measuring Time is set to 8s.And instrument automatic analysis, detects data and averages, and draws size-grade distribution result.
Further, as used Ma Erwen 2000 laser particle analyzers, the software that can provide by it provides the information such as the corresponding percent by volume of particle diameters at different levels and distribution curve, median particle diameter, the corresponding particle diameter of cumulative volumes at different levels.By the test data to this experiment, process and obtain Fig. 4.In core sample, suspending components content reaches 90%, and this meets the characteristic feature of tempestite in granularity, so just can determine that rock sample to be measured is tempestite.
Preferably, after step S203, before step S204, also comprise: to described, treat that spectrum analysis data and described depth value carry out compacting recovery, and use described in overcompaction recovery and treat that spectrum analysis data and described depth value mated with the astronomical cycle of obtaining in advance.
Except sediment being surveyed the impact causing in year due to the sedimental data of event, formation compaction effect can be surveyed and impact in year sediment equally.Therefore, in order to overcome this impact, should be after step S203, the sedimentary record data (treating spectrum analysis data) of rejecting event sediment impact have namely been obtained afterwards, also to treat spectrum analysis data and corresponding depth value carries out compacting recovery, to eliminate formation compaction effect, sediment be surveyed to the impact that year brings.
Concrete, to described, treat that spectrum analysis data and described depth value carry out compacting and recover to comprise:
According to the spacing of the depth value of different depth and the adjacent degree of depth, calculate the sedimentary record data of the different depth after compacting recovers;
The sedimentary record data of the different depth after compacting is recovered resample, and to obtain described in overcompaction recovery, treat spectrum analysis data.
Further, according to the spacing of the depth value of different depth and the adjacent degree of depth, the sedimentary record data of calculating the different depth after compacting recovers comprise:
According to the depth value of different depth, calculate the sedimental compaction coefficient of different depth;
According to the spacing of described compaction coefficient and the adjacent degree of depth obtained in advance, calculate the sedimentary record data of the different depth after compacting recovers.
And the formula that calculates the sedimental compaction coefficient of different depth is: k=b-a*lnH, wherein k is compaction coefficient, and a and b are constant, and determine by voidage, pore texture and the rock composition of the deposition of designated depth, and H is the degree of depth.It should be noted that, if rock is mud shale, a can value 0.1, and b can value 1.46, if rock is sandstone, a can value 0.08, and b can value 1.37.A is to obtain by consulting relevant data to the value of b.
Wherein, the sampling interval of consecutive point (spacing of the adjacent degree of depth) just obtains the primary deposit thickness between the adjacent degree of depth divided by compaction coefficient, again the original thickness of different depth is stacked up and obtain a new degree of depth sequence, the new degree of depth sequence being comprised of a plurality of depth values does not become with the sedimentary record data (natural gamma value) of the degree of depth sequence same sequence number being comprised of a plurality of depth values before (from the sequence forming according to the depth value of a plurality of degree of depth of depth size arrangement, the corresponding numbering in sequence of each depth value) correspondence.
New degree of depth sequence is carried out to interpolation with corresponding natural gamma, generate the sample sequence at even depth interval, namely through sediment compacting, recover, and the sequence that forms with the corresponding degree of depth of the sedimentary record data after equidistantly resampling.
The sedimentary record data of obtaining of take describe as natural gamma numerical value as example, and sand shale has significant difference in natural gamma value, and sandstone natural gamma numerical value is less, and mud stone natural gamma numerical value is larger.By the contrast reference to natural gamma numerical value and core data, select suitable natural gamma value as the separatrix of sand shale, as get average 90 as cut off value, for the mud shale that is considered as that is greater than 90, with k=1.46-0.1*lnH formula, calculate compaction coefficient, the sandstone that is considered as for being less than 90, calculates compaction coefficient k with k=1.37-0.08*lnH formula.
Except event deposition and formation compaction effect meeting impact sedimental surveys year, at the uniform velocity to deposit that the bottom age causing cannot accurately measure be also an important reason due to non-.In order to overcome this problem, described to described after spectrum analysis data are carried out compacting recovery, described, described treating also comprised before spectrum analysis data and described depth value mated with the astronomical cycle of obtaining in advance:
Whether the amplitude of variation of the sedimentary record data of contrast different depth scope conforms to, and described depth range comprises continuous a plurality of depth values;
If not, the sedimentary record data of different depth scope are normalized to computing, to obtain sedimentary record data and the corresponding depth value after normalization computing;
Sedimentary record data after described normalization computing and corresponding depth value are equidistantly resampled, to obtain sedimentary record data and the corresponding depth value at even depth interval;
Use sedimentary record data and the corresponding depth value at described even depth interval to mate with the astronomical cycle of obtaining in advance.
Concrete, as shown in Figure 5, can significantly observe, gamma ray curve from the angle of mean value size can be divided into three sections (as corresponding in the identification section of figure middle and upper part, be divided into three sections of left, center, rights), the horizontal ordinate of wherein changing plan is depth value, ordinate is the numerical value of GR data.So, the amplitude of variation that just can judge the sedimentary record data of different depth scope does not conform to, namely uneven, now tackles this kind of data and carries out rate of sedimentation normalization computing.Certainly, in this application, sampled point is more intensive, and the data that obtain are also just more accurate, and so, the analysis result going out by data observation is also just more accurate.
The amplitude of variation that can confirm the sedimentary record data (natural gamma data) of different depth scope by Fig. 5 does not conform to.Now, the sedimentary record data of tackling different depth scope are normalized computing, to obtain sedimentary record data and the corresponding depth value after normalization computing.Concrete, natural gamma is worth size and rate of sedimentation to be inverse relation, adopt following formula v=(150-ave)/100 pair of three sections of rate of sedimentation to be normalized computing, according to mean value, obtain every section of corresponding rate of sedimentation of looking, by the sampling interval of every section divided by the depth interval sequence obtaining depending on rate of sedimentation after rate of sedimentation homogenization, its depth interval is added up and obtains new degree of depth sequence, and wherein ave is the numerical value of natural gamma data.Certainly, if use the sedimentary record data (as sediment magnetic susceptibility data and sediment remanence data) of other kinds as computational data, also the ave correspondence in formula should be adjusted to corresponding numerical value unit.
Afterwards, new degree of depth sequence the preceding paragraph being obtained and sedimentary record data (natural gamma data) are carried out the sample sequence that equidistantly resampling obtains the even depth interval after rate of sedimentation is proofreaied and correct.And use the sample sequence (comprising sedimentary record data and corresponding depth value) at even depth interval to mate with the astronomical cycle of obtaining in advance, namely the sedimentary record data in the sample sequence of reciprocity depth interval (namely through overcompaction recover and rate of sedimentation normalization computing after treat spectrum analysis data) carry out spectrum analysis, to obtain corresponding spectrum analysis curve (degree of depth and sedimentary record data institute constituent curve), and find out and consult frequency corresponding to astronomical cycle of obtaining by the frequency curve of sedimentary record data, to processing with the data filtering of astronomical cycle respective frequencies in sedimentary record data, again filtered sedimentary record data and curves is matched on the theoretical curve in astronomical cycle, (each point in astronomical cycle has known precise time, the age and the degree of depth that by data consultation, obtain the boundary on stratum are also known, time value by stratigraphic boundary place matched with the astronomical cycle, according to the degree of depth at the degree of depth in sedimentary record data and stratigraphic boundary place, match again, and then filtered sedimentary record data and curves is matched on the theoretical curve in astronomical cycle).The astronomical cycle in sedimentary record data and curves and ancient times is compared, the determinacy based on the astronomical cycle, just can obtain the sedimental age of respective depth again.
Concrete, the formed curve of sample sequence of reciprocity depth interval carries out spectrum analysis, and result is as shown in Figure 6, obtain frequency and be five main peaks (5 stains in figure) of 0.00859,0.00982,0.0175,0.0249,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 cycle, with theoretical cycle 125kyr, 96kyr, 54kyr, 40kyr, the proportionate relationship 1,0.768,0.432 of 23kyr, 0.32,0.184 has good corresponding relation.Therefore above-mentioned six cycles are records in sediment of astronomical cycle.Demarcate by known stratum separation age and astronomical cycle, can obtain natural gamma 125kyr filtering and astronomical cycle 125kyr filtering matching result, as Fig. 7, natural gamma 125kyr filtering is mated with astronomical cycle 125kyr filtering, and natural gamma filtering and astronomical periodical filtering have good matching.
On this basis, by natural gamma 125kyr filtering and astronomical cycle 125kyr filtering matching result, can obtain the time that event deposition (as tempestite) is occurred.The astronomical cycle of 125kyr is identical with the time of the corresponding crest of natural gamma filter curve.The time of the corresponding sampled point between crest can obtain by 40kyr filter curve, and method is as follows:
First determine tempestite development position in natural gamma filter curve between any two peak values, as determined the corresponding time of Fig. 8 middle and upper part branch A, point A is between natural gamma 125kyr filter curve second and the 3rd peak value, second and the 3rd peak value of the corresponding red astronomical periodical filtering curve of this two peak value, by second time corresponding to peak value of the known blue curve of second time to peak of red curve, blue curve and yellow curve are synchronous, second peak value of blue curve can correspond to straight down and in yellow curve, obtain a B, again A point is corresponded to lower yellow curve and obtain C, corresponding 40kyr periodicity P between estimation B and C point, the corresponding time of A point is: T+P*40.So just, can draw the time point that the event sediment of rejecting occurs.
By rejecting impact that event deposition causes, eliminate the impact that bottom compaction causes and solved the non-at the uniform velocity impact that deposition causes before and after data spectrum contrast, as shown in Figure 8, can find out:
1, in the data spectrum of the natural gamma data after processing, can identify above-mentioned five frequency main peaks, and in raw data frequency spectrum, be difficult to the peak value of identification and astronomical cycle respective frequencies.
2, five frequency peak that identify are compared with raw data respective frequencies peak value, and its energy all improves more than 30%.
3, in the data spectrum after processing, dispersion phenomenon is effectively suppressed, and frequency dispersion rate reduces more than 50%.
Survey year method based on well logging Milankovitch Cycle Hypothesis analytical approach provided by the present invention is compared with traditional survey year method, the error of surveying year can be dwindled even up to a million years of hundreds of thousands year.
Survey year method based on well logging Milankovitch Cycle Hypothesis analytical approach provided by the present invention, by the core sample of the objective interval different depth getting is analyzed, to determine whether the sedimental core sample of different depth has event sediment to produce, if had, corresponding, the sedimentary record data (as GR data) of the rejecting event sediment place degree of depth, to obtain, treat spectrum analysis data and corresponding depth value, and then will treat that spectrum analysis data are carried out compacting recovery and rate of sedimentation normalization is corrected, finally treat spectrum analysis data and the depth value that through overcompaction, recover and rate of sedimentation normalization is corrected were mated with the astronomical cycle of obtaining in advance, to obtain sedimental generation time, further avoided, due to what event deposition was brought, deposition is surveyed to the impact in year.
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.Within the spirit and principles in the present invention all, any modification of doing, be equal to replacement, improvement etc., within all should being included in protection scope of the present invention.
Claims (10)
1. the survey year method based on well logging Milankovitch Cycle Hypothesis analytical approach, is characterized in that, comprising:
Obtain the sedimental sedimentary record data of different depth and the corresponding depth value of objective interval;
Obtain the sedimental core sample of different depth and the depth value of described objective interval, and described core sample is analyzed, to determine the degree of depth at event sediment place;
In described sedimentary record data, reject the sedimentary record data of the event sediment place degree of depth, to obtain, treat spectrum analysis data and described depth value;
By described, treat that spectrum analysis data and described depth value mated with the astronomical cycle of obtaining in advance, the time being produced to obtain sediment.
2. the survey year method based on well logging Milankovitch Cycle Hypothesis analytical approach according to claim 1, is characterized in that, described described core sample analyzes is comprised: employing core observation method and/or sizing assay test method are carried out described core sample.
3. the survey year method based on well logging Milankovitch Cycle Hypothesis analytical approach according to claim 1, it is characterized in that, in described acquisition after spectrum analysis data, described by described treat that spectrum analysis data and described depth value mated with the astronomical cycle of obtaining in advance before, also comprise:
To described, treat that spectrum analysis data and described depth value carry out compacting recovery, and use described in overcompaction recovery and treat that spectrum analysis data and described depth value mated with the astronomical cycle of obtaining in advance.
4. the survey year method based on well logging Milankovitch Cycle Hypothesis analytical approach according to claim 3, is characterized in that, describedly to described, treats that spectrum analysis data and described depth value carry out compacting recovery and comprise:
According to the spacing of the depth value of different depth and the adjacent degree of depth, calculate the sedimentary record data of the different depth after compacting recovers;
The sedimentary record data of the different depth after compacting is recovered resample, and to obtain described in overcompaction recovery, treat spectrum analysis data.
5. the survey year method based on well logging Milankovitch Cycle Hypothesis analytical approach according to claim 4, is characterized in that, described according to the spacing of the depth value of different depth and the adjacent degree of depth, and the sedimentary record data of calculating the different depth after compacting recovers comprise:
According to the depth value of different depth, calculate the sedimental compaction coefficient of different depth;
According to the spacing of described compaction coefficient and the adjacent degree of depth obtained in advance, calculate the sedimentary record data of the different depth after compacting recovers.
6. the survey year method based on well logging Milankovitch Cycle Hypothesis analytical approach according to claim 5, it is characterized in that, the formula of the sedimental compaction coefficient of described calculating different depth is: k=b-a*lnH, wherein k is compaction coefficient, a and b are constant, and voidage, pore texture and rock composition by the deposition of designated depth determine, H is the degree of depth.
7. the survey year method based on well logging Milankovitch Cycle Hypothesis analytical approach according to claim 3, it is characterized in that, described to described after spectrum analysis data are carried out compacting recovery, described, described treating also comprised before spectrum analysis data and described depth value mated with the astronomical cycle of obtaining in advance:
Whether the amplitude of variation of the sedimentary record data of contrast different depth scope conforms to, and described depth range comprises continuous a plurality of depth values;
If not, the sedimentary record data of different depth scope are normalized to computing, to obtain sedimentary record data and the corresponding depth value after normalization computing;
Sedimentary record data after described normalization computing and corresponding depth value are equidistantly resampled, to obtain sedimentary record data and the corresponding depth value at even depth interval;
Use sedimentary record data and the corresponding depth value at described even depth interval to mate with the astronomical cycle of obtaining in advance.
8. the survey year method based on well logging Milankovitch Cycle Hypothesis analytical approach according to claim 1, is characterized in that, described sedimentary record data comprise sediment magnetic susceptibility data, sediment remanence data, natural gamma data.
9. the survey year method based on well logging Milankovitch Cycle Hypothesis analytical approach according to claim 2, is characterized in that, described sizing assay test method comprises:
To described core sample carry out broken sample, remove organic matter, remove cementing matter, pickling, grinding and discrete particles process, to obtain sample to be measured;
By laser particle analyzer preheating 30 minutes, and described laser particle analyzer is carried out adjustment and measures background;
With glass bar, the suspension of described sample to be measured is stirred;
The suspension of described sample to be measured being added in the sample cell of laser particle analyzer, is that the obscurity of the suspension of described sample to be measured reaches 10%-20%;
Suspension in described sample cell is carried out to ultrasonic processing 5min;
Repeatedly measure the suspension in described sample cell, in the hope of the mean value of measurement data;
According to the mean value of described measurement data, judge whether described core sample is event sediment;
If so, the degree of depth of determining described core sample place is the degree of depth at event sediment place.
10. the survey year method based on well logging Milankovitch Cycle Hypothesis analytical approach according to claim 9, it is characterized in that, described measurement data comprises: the corresponding percent by volume of particle diameters at different levels and distribution curve, median particle diameter and the corresponding particle diameter of 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 true CN104089964A (en) | 2014-10-08 |
| CN104089964B 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 (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108427146A (en) * | 2017-02-13 | 2018-08-21 | 中国石油化工股份有限公司 | A kind of method and system of Milankovitch Cycles in identification sedimentary formation |
| CN108919369A (en) * | 2018-05-16 | 2018-11-30 | 中国海洋石油集团有限公司 | A kind of method of quantitative judge Milankovitch Cycles |
| CN109001840A (en) * | 2018-07-18 | 2018-12-14 | 中国地质大学（北京） | A kind of mud shale genetic analysis method and device |
| CN109031443A (en) * | 2017-06-09 | 2018-12-18 | 中国石油化工股份有限公司 | The erosion restoration method in year is determined based on reservoir rock and caprock |
| CN110231660A (en) * | 2019-05-15 | 2019-09-13 | 防灾科技学院 | The method of the dislocation amount of buried active fault is determined using magnetic susceptibility log method |
| CN110454155A (en) * | 2019-08-23 | 2019-11-15 | 防灾科技学院 | A kind of method determining the Quaternary Stratigraphic age using magnetic susceptibility log method |
| CN112950551A (en) * | 2021-02-05 | 2021-06-11 | 自然资源部第一海洋研究所 | Ocean sediment rock core age rapid calibration method |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107884830B (en) * | 2016-09-29 | 2020-03-24 | 中国石油化工股份有限公司 | Sedimentary formation cycle dividing method and device |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5999485A (en) * | 1996-12-13 | 1999-12-07 | Anstey; Nigel A. | Method of distinguishing geological sequences and their boundaries |
-
2014
- 2014-07-23 CN CN201410353207.7A patent/CN104089964B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5999485A (en) * | 1996-12-13 | 1999-12-07 | Anstey; Nigel A. | Method of distinguishing geological sequences and their boundaries |
Non-Patent Citations (6)
| Title |
|---|
| HELGE W. ARZ，ET AL.: "Correlated Millennial-Scale Changes in Surface Hydrography and Terrigenous Sediment Yield Inferred from Last-Glacial Marine Deposits off Northeastern Brazil", 《QUATERNARY RESEARCH》, vol. 50, 31 December 1998 (1998-12-31), pages 157 - 166 * |
| YUANFU ZHANG，ZAIXING，JIANG: "outcrop characterization of an early miocene slop fan system，chelif basin，algeria", 《ENERGY EXPLORATION&EXPLOITATION》, vol. 29, no. 5, 31 December 2011 (2011-12-31), pages 633 - 646 * |
| ZHANG YUANFU，ET AL.: "comparative analysis of sequence characteristics among different superimposed stages of the chelif basin,algeria", 《ACTA GEOLOGICA SINICA》, vol. 83, no. 6, 31 December 2009 (2009-12-31), pages 1041 - 1051 * |
| 张元福，李洪奇: "基于测井资料的层序地层划分方法研究", 《中国西部油气地质》, vol. 2, no. 4, 31 December 2006 (2006-12-31), pages 427 - 431 * |
| 邱桂强，刘军鄂，帅萍: "米氏旋回基本原理及其在陆相湖盆分析中的应用前景", 《油气地质与采收率》, vol. 8, no. 5, 31 October 2001 (2001-10-31), pages 5 - 9 * |
| 郎旭娟，余继峰等: "济阳坳陷车镇凹陷沙河街组旋回地层学研究", 《煤炭学报》, vol. 37, no. 11, 30 November 2012 (2012-11-30), pages 1901 - 1908 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108427146A (en) * | 2017-02-13 | 2018-08-21 | 中国石油化工股份有限公司 | A kind of method and system of Milankovitch Cycles in identification sedimentary formation |
| CN108427146B (en) * | 2017-02-13 | 2019-08-30 | 中国石油化工股份有限公司 | A kind of method and system of Milankovitch Cycles in identification sedimentary formation |
| CN109031443A (en) * | 2017-06-09 | 2018-12-18 | 中国石油化工股份有限公司 | The erosion restoration method in year is determined based on reservoir rock and caprock |
| CN109031443B (en) * | 2017-06-09 | 2021-05-14 | 中国石油化工股份有限公司 | Method for recovering denudation amount based on fixed years of Milnaciaceae Virginia gyrus |
| CN108919369A (en) * | 2018-05-16 | 2018-11-30 | 中国海洋石油集团有限公司 | A kind of method of quantitative judge Milankovitch Cycles |
| CN109001840A (en) * | 2018-07-18 | 2018-12-14 | 中国地质大学（北京） | A kind of mud shale genetic analysis method and device |
| CN110231660A (en) * | 2019-05-15 | 2019-09-13 | 防灾科技学院 | The method of the dislocation amount of buried active fault is determined using magnetic susceptibility log method |
| CN110454155A (en) * | 2019-08-23 | 2019-11-15 | 防灾科技学院 | A kind of method determining the Quaternary Stratigraphic age using magnetic susceptibility log 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 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104089964B (en) | 2016-08-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104089964A (en) | Dating method based on logging Milankovitch cycle analysis method | |
| Chen et al. | Multiparameter adjoint tomography of the crust and upper mantle beneath East Asia: 1. Model construction and comparisons | |
| CN105675635B (en) | Tight rock component relative amount and brittleness index determine method and apparatus | |
| CN103809217A (en) | Method for calculating shale content of high and low gamma interaction reservoirs | |
| CN105275456B (en) | A kind of method utilizing well-log information identification high-quality mud shale | |
| 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 | |
| CN106285660A (en) | A kind of multilamellar sandstone oil reservoir Low Resistivity Reservoir Identification Methods and device | |
| CN103616731B (en) | Altered Volcanic Rocks at some Place Effective Reservoirs defining method and device in a kind of oil-gas exploration | |
| CN103744109B (en) | Method for identifying a weathering crust structure of clastic rock in a area of covering no well | |
| CN103675907A (en) | AVO inversion hydrocarbon detection method based on petrographic constraints | |
| CN102789003A (en) | Method for inverting stratum parameters by array induction logging data and device thereof | |
| CN104047598A (en) | Method for predicating productivity of nonhomogeneity ancient karst carbonate reservoir | |
| 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 | |
| CN109298465A (en) | Phased sand body main structure body block sandstone quantitatively distinguishes method and device | |
| CN102914797A (en) | Method and device for acquiring anisotropy coefficient of stratum | |
| Venkataramanan et al. | New method to estimate porosity more accuately from NMR data with short relaxation times | |
| CN105590018A (en) | Oil-water layer identification method for sandstone and mudstone thin interbed oil reservoir | |
| CN112034514A (en) | Shale oil-gas system based on seismic sedimentology theory and identification method thereof | |
| Gartmair et al. | One size does not fit all: Refining zircon provenance interpretations via integrated grain shape, geochronology, and Hf isotope analysis | |
| Feng et al. | A new method for logging identification of fluid properties in tight sandstone gas reservoirs based on gray correlation weight analysis—A case study of the Middle Jurassic Shaximiao Formation on the eastern slope of the Western Sichuan Depression, China | |
| CN113960288B (en) | Quantitative evaluation method for source-sink system | |
| CN112878999B (en) | Method and device for calculating water saturation of anisotropic stratum | |
| Liu et al. | Seismic-Geological Integrated Study on Sedimentary Evolution and Peat Accumulation Regularity of the Shanxi Formation in Xinjing Mining Area, Qinshui Basin | |
| Rotimi et al. | Reservoir characterization and modeling of lateral heterogeneity using multivariate analysis |
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 |