CN106680879A

CN106680879A - Method and device for sedimentary cycle division

Info

Publication number: CN106680879A
Application number: CN201611201548.8A
Authority: CN
Inventors: 贺文; 杨平; 岳鹏飞; 吴嘉兴
Original assignee: China National Petroleum Corp; BGP Inc
Current assignee: China National Petroleum Corp; BGP Inc
Priority date: 2016-12-22
Filing date: 2016-12-22
Publication date: 2017-05-17
Anticipated expiration: 2036-12-22
Also published as: CN106680879B

Abstract

The embodiment of the invention provides a method and device for sedimentary cycle division. The method comprises the following steps: obtaining a peak frequency attribute curve according to seismic data; performing frequency division smoothing of the peak frequency property curve, and obtaining a low-frequency attribute band-pass frequency component curve and a high-frequency attribute band-pass frequency component curve; performing broken line processing of the low-frequency attribute band-pass frequency component curve and the high-frequency attribute band-pass frequency component curve, and obtaining a first sedimentary cycle curve and a second sedimentary cycle curve. According to the scheme, the frequency division processing is performed through a filter, the curve after frequency division is subjected to broken line processing to obtain a corresponding sedimentary cycle curve. The technical problems are solved that the error is large, the efficiency is low and the qualification or semiquantitative analysis is only performed in the current sedimentary cycle division, and therefore the technical effects are achieved that the seismic data is fully utilized to rapidly and accurately perform the sedimentary cycle division.

Description

The method and apparatus that the cycle of sedimentation is divided

Technical field

The application is related to oil-gas exploration technical field, the method and apparatus that more particularly to a kind of cycle of sedimentation is divided.

Background technology

During oil-gas exploration is carried out, when cycle of sedimentation division is carried out to formation at target locations, generally require first to heavy Product cycle is analyzed.Wherein, the above-mentioned cycle of sedimentation refers to the result of nature periodicity eustasy.Specifically： Have in time periodically, there is in the origin cause of formation similarity, cyclicity is shown in geology, is shown in geophysics Certain response characteristic.For example, in seismic wave field, the regular of seismic wave frequency spectrum feature changes.Therefore, by the cycle of sedimentation Analysis, can further realize analysis, SEQUENCE STRATIGRAPHIC research or even the prediction to reservoir for seismic facies and sedimentary facies.

At present, in order to divide depositional cycle curve, conventional method is to obtain depositional cycle curve according to well-log information.Or Person obtains depositional cycle curve according to seismic data where well-log information is lacked.Wherein, it is described to be obtained according to seismic data The method of depositional cycle curve, typically according to the attribute for calculating curve, according to the trend of curve, is artificially subjectively rolled over Line, and then obtain depositional cycle curve.Using this depositional cycle curve division methods, when being embodied as, often exist Error is big, efficiency is low, can only qualitative or semi-quantitative analyses technical problem.

For the problems referred to above, effective solution is not yet proposed at present.

The content of the invention

The embodiment of the present application provides the method and apparatus that a kind of cycle of sedimentation is divided, and is divided with solving the existing cycle of sedimentation The technical problem that error is big, efficiency is low existed when dividing depositional cycle curve in method.

The embodiment of the present application provides a kind of method that cycle of sedimentation is divided, including：

Time frequency analysis are carried out to seismic data, time-frequency data volume is obtained；

According to the time-frequency data volume, crest frequency attribute curve is obtained；

Frequency-division filter is carried out to the crest frequency attribute curve, the band logical frequency component curve and height of low frequency attribute is obtained The band logical frequency component curve of frequency attribute；

Broken line process is carried out to the band logical frequency component curve of the low frequency attribute, the first depositional cycle curve is obtained, Broken line process is carried out to the band logical frequency component curve of the high frequency nature, the second depositional cycle curve is obtained；

According to first depositional cycle curve and second depositional cycle curve, cycle of sedimentation division is carried out.

In one embodiment, frequency-division filter is carried out to the crest frequency attribute curve, obtains the band of low frequency attribute The band logical frequency component curve of passband rate component curve and high frequency nature, including：

According to cycle of sedimentation rank to be studied, mid frequency is determined；

According to the mid frequency and pass band width, band filter is designed；

Frequency-division filter is carried out to the crest frequency attribute curve by the band filter, the band of low frequency attribute is obtained The band logical frequency component curve of passband rate component curve and high frequency nature.

In one embodiment, according to the mid frequency and pass band width, band filter is designed, including：

Using the difference between the half of the mid frequency and the pass band width as lowpass frequency, with center frequency The additive value of rate and the half of the pass band width is used as high-pass equipment；

According to the lowpass frequency and the high-pass equipment, the band filter is designed.

In one embodiment, broken line process is carried out to the band logical frequency component curve of the low frequency attribute, is obtained First depositional cycle curve, carries out broken line process to the band logical frequency component curve of the high frequency nature, obtains the second deposition Cycle curve, including：

According to the band logical frequency component curve of the low frequency attribute, low frequency attribute most value sequence is obtained；

Numerical value process is carried out to the low frequency attribute most value sequence, the low frequency attribute most value sequence after being processed；

According to the low frequency attribute most value sequence after the process, first cycle of sedimentation is obtained by linear interpolation bent Line；

According to the band logical frequency component curve of the high frequency nature, high frequency nature most value sequence is obtained；

Numerical value process is carried out to the high frequency nature most value sequence, the most value sequence of the high frequency nature after being processed；

According to the most value sequence of the high frequency nature after the process, second cycle of sedimentation is obtained by linear interpolation bent Line.

In one embodiment, numerical value process is carried out to the low frequency attribute most value sequence, the low frequency after being processed Attribute most value sequence, including：

According to the cycle of sedimentation rank to be studied, the first cycle threshold value is determined；

According to the low frequency attribute most value sequence, the meansigma methodss of low frequency attribute most value sequence are determined；

It is according to putting in order for low frequency attribute most value sequence, adjacent with two in the low frequency attribute most value sequence respectively Most being worth and subelement being processed as low frequency, process subelement to each low frequency in the low frequency attribute most value sequence respectively two are most It is poor that value is carried out, and obtains the absolute value that corresponding low frequency processes the difference of subelement；

Absolute value, the meansigma methodss of the low frequency attribute most value sequence and the institute of the difference of subelement are processed according to the low frequency The first cycle threshold value is stated, numerical value process is carried out to the low frequency attribute most value sequence, the low frequency attribute after being processed most is worth sequence Row.

In one embodiment, the absolute value of difference of subelement, the low frequency attribute are processed most according to the low frequency The meansigma methodss of value sequence and the first cycle threshold value, carry out numerical value process to the low frequency attribute most value sequence, including：

In the case where the absolute value of the difference of low frequency process subelement is more than or equal to the first cycle threshold value, then Retain the low frequency and process two most values and the distinguished corresponding time coordinate of described two most values in subelement；

The absolute value that the difference of subelement is processed in the low frequency is less than the first cycle threshold value, and the low frequency is processed In the case that two in subelement are most worth both less than equal to the meansigma methodss of the low frequency attribute most value sequence, then retain described low Frequency process subelement in two most values in less most value and it is described it is less be most worth corresponding to time coordinate；

The absolute value that the difference of subelement is processed in the low frequency is less than the first cycle threshold value, and the low frequency is processed In the case that two in subelement are most worth the both greater than meansigma methodss of the low frequency attribute most value sequence, then retain at the low frequency In two most values in reason subelement larger most value and it is described it is larger be most worth corresponding to time coordinate；

The absolute value that the difference of subelement is processed in the low frequency is less than the first cycle threshold value, and the low frequency is processed Two in subelement be most worth in a larger meansigma methodss of the most value more than the low frequency attribute most value sequence, one less to be most worth In the case of meansigma methodss less than the low frequency attribute most value sequence, then retain in two most values that the low frequency is processed in sub- list The most value larger with the absolute value of the meansigma methodss difference of the low frequency attribute most value sequence and corresponding time coordinate.

In one embodiment, numerical value process is carried out to the high frequency nature most value sequence, the high frequency after being processed Attribute most value sequence, including：

According to the cycle of sedimentation rank to be studied, the second cycle threshold value is determined；

According to the high frequency nature most value sequence, the meansigma methodss of high frequency nature most value sequence are determined；

It is according to putting in order for high frequency nature most value sequence, adjacent with two in the high frequency nature most value sequence respectively Most it is worth as high-frequency therapeutic treatment subelement, respectively to two of each high-frequency therapeutic treatment subelement in the high frequency nature most value sequence most It is poor that value is carried out, and obtains the absolute value of the difference of corresponding high-frequency therapeutic treatment subelement；

Absolute value, the meansigma methodss of the high frequency nature most value sequence and institute according to the difference of the high-frequency therapeutic treatment subelement The second cycle threshold value is stated, numerical value process is carried out to the high frequency nature most value sequence, the high frequency nature after being processed most is worth sequence Row.

In one embodiment, according to the difference of the high-frequency therapeutic treatment subelement absolute value, the high frequency nature are most The meansigma methodss of value sequence and the second cycle threshold value, carry out numerical value process to the high frequency nature most value sequence, including：

In the case where the absolute value of the difference of the high-frequency therapeutic treatment subelement is more than or equal to the second cycle threshold value, then Retain two most values in the high-frequency therapeutic treatment subelement and described two be most worth corresponding time coordinate；

The second cycle threshold value, and the high-frequency therapeutic treatment are less than in the absolute value of the difference of the high-frequency therapeutic treatment subelement In the case that two in subelement are most worth both less than equal to the meansigma methodss of the high frequency nature most value sequence, then retain the height Frequency process subelement in two most values in less most value and it is described it is less be most worth corresponding to time coordinate；

The second cycle threshold value, and the high-frequency therapeutic treatment are less than in the absolute value of the difference of the high-frequency therapeutic treatment subelement In the case that two in subelement are most worth the both greater than meansigma methodss of the high frequency nature most value sequence, then retain the high frequency treatment In two most values in reason subelement larger most value and it is described it is larger be most worth corresponding to time coordinate；

The second cycle threshold value, and the high-frequency therapeutic treatment are less than in the absolute value of the difference of the high-frequency therapeutic treatment subelement Two in subelement be most worth in a larger meansigma methodss of the most value more than the high frequency nature most value sequence, one less to be most worth In the case of meansigma methodss less than the high frequency nature most value sequence, then retain in two most values in the sub- list of the high-frequency therapeutic treatment The most value larger with the absolute value of the meansigma methodss difference of the high frequency nature most value sequence and corresponding time coordinate.

In one embodiment, broken line process is carried out to the band logical frequency component curve of the low frequency attribute, is obtained First depositional cycle curve, carries out broken line process to the band logical frequency component curve of the high frequency nature, obtains the second deposition After cycle curve, methods described also includes：

Carry out multiple broken lineization to process to first depositional cycle curve, the first cycle curve after process meets Preset requirement, will meet the first depositional cycle curve of preset requirement as first depositional cycle curve,

Carry out multiple broken lineization to process to second depositional cycle curve, the second cycle curve after process meets Preset requirement, will meet the second depositional cycle curve of preset requirement as second depositional cycle curve.

Based on identical inventive concept, the embodiment of the present application additionally provides the device that a kind of cycle of sedimentation is divided, including：

First acquisition module, for carrying out time frequency analysis to seismic data, obtains time-frequency data volume；

Second acquisition module, for according to the time-frequency data volume, obtaining crest frequency attribute curve；

Frequency-division filter module, for carrying out frequency-division filter to the crest frequency attribute curve, obtains the band of low frequency attribute The band logical frequency component curve of passband rate component curve and high frequency nature；

Broken line processing module, for carrying out broken line process to the band logical frequency component curve of the low frequency attribute, obtains To the first depositional cycle curve, broken line process is carried out to the band logical frequency component curve of the high frequency nature, obtain second and sink Product cycle curve；

Cycle of sedimentation division module, for according to first depositional cycle curve and second depositional cycle curve, Carry out cycle of sedimentation division.

In one embodiment, the frequency-division filter module includes：

Determine filter cell, for according to cycle of sedimentation rank to be studied, determining mid frequency, and according in described Frequency of heart and pass band width, determine band filter；

Frequency-division filter unit, for carrying out frequency dividing filter to the crest frequency attribute curve by the band filter Ripple, obtains the band logical frequency component curve of the band logical frequency component curve and high frequency nature of low frequency attribute.

In one embodiment, the broken line processing module includes：

Low frequency most value sequence acquiring unit, for the band logical frequency component curve according to the low frequency attribute, obtains low frequency Attribute most value sequence；

Low frequency most value sequence processing unit, for carrying out numerical value process to the low frequency attribute most value sequence, is processed Low frequency attribute most value sequence afterwards；

Low frequency broken line processing unit, for according to the low frequency attribute most value sequence after the process, by linear interpolation Obtain first depositional cycle curve；

High frequency most value sequence acquiring unit, for the band logical frequency component curve according to the high frequency nature, obtains high frequency Attribute most value sequence；

High frequency most value sequence processing unit, for carrying out numerical value process to the high frequency nature most value sequence, is processed The most value sequence of high frequency nature afterwards；

High frequency broken line processing unit, for according to the most value sequence of the high frequency nature after the process, by linear interpolation Obtain second depositional cycle curve.

In the embodiment of the present application, first passing through wave filter carries out scaling down processing, then the curve after frequency dividing is rolled over respectively Lineization process, obtains corresponding depositional cycle curve.Existing cycle of sedimentation division side is solved by way of broken line is processed Carry out in method that the error that exists when depositional cycle curve is divided is big, efficiency is low, can only qualitative or semi-quantitative analyses technical problem, Reach and made full use of seismic data, technique effect that is quick, accurately carrying out cycle of sedimentation division.

Description of the drawings

In order to be illustrated more clearly that the embodiment of the present application or technical scheme of the prior art, below will be to embodiment or existing Accompanying drawing to be used needed for having technology description is briefly described, it should be apparent that, drawings in the following description are only this Some embodiments described in application, for those of ordinary skill in the art, in the premise for not paying creative labor Under, can be with according to these other accompanying drawings of accompanying drawings acquisition.

Fig. 1 is the process chart of the method divided according to the cycle of sedimentation of the embodiment of the present application；

Fig. 2 is the composition structure chart of the device divided according to the cycle of sedimentation of the embodiment of the present application；

Fig. 3 is the process schematic that the method, device that the cycle of sedimentation is divided is provided using the embodiment of the present application；

Fig. 4 is the trial zone seismic profile that the method, device process that the cycle of sedimentation is divided is provided using the embodiment of the present application Figure；

When Fig. 5 is that single track data, the single track of the method, device acquisition that the cycle of sedimentation is divided is provided using the embodiment of the present application Frequency analysis and crest frequency attribute schematic diagram；

Fig. 6 be using the embodiment of the present application provide the cycle of sedimentation divide method, device obtain frequency-division filter after when Frequency attributed graph；

Fig. 7 is the depositional cycle curve that the method, device acquisition that the cycle of sedimentation is divided is provided using the embodiment of the present application.

Specific embodiment

In order that those skilled in the art more fully understand the technical scheme in the application, below in conjunction with the application reality The accompanying drawing in example is applied, the technical scheme in the embodiment of the present application is clearly and completely described, it is clear that described enforcement Example is only some embodiments of the present application, rather than the embodiment of whole.Based on the embodiment in the application, this area is common The every other embodiment obtained under the premise of creative work is not made by technical staff, should all belong to the application protection Scope.

In view of existing cycle of sedimentation division methods, when being embodied as, it usually needs first divide the corresponding cycle of sedimentation Curve.It is typically all the artificial subjectivity according to the frequency attribute curve for calculating, according to the trend of curve to divide depositional cycle curve Carry out broken line.Therefore, when depositional cycle curve is divided in causing existing depositional cycle curve division methods often there is mistake Difference is big, efficiency is low, can only qualitative or semi-quantitative analyses technical problem.Basic reason for producing above-mentioned technical problem, this Application consideration first can obtain crest frequency attribute curve according to seismic data, then carry out frequency dividing filter to the curve by wave filter Ripple, carries out corresponding broken lineization process respectively to curve is obtained after filtering, obtains corresponding depositional cycle curve.So as to solve When depositional cycle curve is divided in existing depositional cycle curve division methods, the error that exists is big, efficiency is low, can only it is qualitative or The technical problem of semi-quantitative analyses, realizes and makes full use of seismic data, technology effect that is quick, accurately carrying out cycle of sedimentation division Really, the purpose of the quantitative analyses cycle of sedimentation has been reached.

Based on above-mentioned thinking thinking, this application provides a kind of method that cycle of sedimentation is divided.Refer to Fig. 1.The application The method that the cycle of sedimentation of offer is divided, specifically may comprise steps of.

Step 101：Time frequency analysis are carried out to seismic data, time-frequency data volume is obtained.

In one embodiment, in order to obtain time-frequency data according to seismic data, can be specifically：To poststack earthquake number Generalized S-transform is carried out according to by road, time-frequency data volume is obtained.It should be noted that the generalized S-transform being previously mentioned in present embodiment It is intended merely to better illustrate the time frequency analysis in the embodiment of the present application, cited one is schematically illustrate.Except being become by broad sense S Outside swap-in row time frequency analysis, time frequency analysis can also be carried out by methods such as the change of Fourier in short-term, Wavelet transformation and S-transformations, Obtain time-frequency data volume.When being embodied as, flexibly appropriate method can be selected over the ground as the case may be, according to specific requirement Shake data carry out time frequency analysis, obtain corresponding time-frequency data volume.In this regard, the application is not construed as limiting.Involved by the application In embodiment, time frequency analysis are carried out to geological data by generalized S-transform mainly.It should be noted that in present embodiment In, from time domain seismic data cube, then frequency domain data body can be a process for rising dimension.

Step 102：According to the time-frequency data volume, crest frequency attribute curve is obtained.

In the present embodiment, obtain instantaneous spectrum crest frequency attribute volume and can be equal to obtaining the crest frequency Attribute curve.I.e. when instantaneous spectrum crest frequency attribute volume is obtained, corresponding peak value can be obtained by simple method Frequency attribute curve.For example, crest frequency attribute curve can be obtained according to instantaneous spectrum crest frequency attribute data volume drawing. Corresponding peak value frequency can be automatically generated by software for calculation while instantaneous spectrum crest frequency attribute volume is calculated Rate attribute curve.In this regard, the application is not construed as limiting.

In one embodiment, in order to according to the acquisition of time-frequency data volume instantaneous spectrum crest frequency attribute volume, i.e. peak Value frequency attribute curve.Can be with：Instantaneous spectrum crest frequency attribute is obtained according to time-frequency data volume first, wherein, the instantaneous spectral peak Value frequency attribute can be the frequency on time orientation corresponding to single-point amplitude spectrum maximum.Therefore, it can by shaking to single-point Width spectral sequence is compared asks for instantaneous spectrum crest frequency attribute, and then obtains wink according to the instantaneous spectrum crest frequency attribute tried to achieve When spectrum peak frequency attribute data volume, i.e. crest frequency attribute curve.It is it should be noted that in the present embodiment, relative, It is the process of a dimensionality reduction from time-frequency domain data volume to instantaneous spectrum crest frequency data volume.

Step 103：Frequency-division filter is carried out to the crest frequency attribute curve, the band logical frequency component of low frequency attribute is obtained The band logical frequency component curve of curve and high frequency nature.

In one embodiment, in order to obtain the band passband of the band logical frequency component curve and high frequency nature of low frequency attribute Rate component curve, when being embodied as, can perform according to following steps：

S1：According to cycle of sedimentation rank to be studied, mid frequency is determined.

In the present embodiment, it should be noted that depending on the cycle rank can be the demand according to user, Sometimes need to study the first cycle, then set corresponding first cycle rank；Sometimes need to study the second cycle, then set Corresponding second cycle rank.In this regard, the application is not construed as limiting.In addition it is also necessary to explanation, mid frequency can be according to rotation Return rank to determine.For example, the mid frequency that the first cycle is selected is relatively low, and the second cycle selects mid frequency higher.Specifically applying In work, to study cycle, often most multipotency reaches three-level to commonly used seismic data.The mid frequency that such as one-level cycle is selected For 10HZ, the mid frequency that second-order cycle is selected is 15HZ, by that analogy.It is concrete how center frequency to be determined according to cycle rank Rate, can flexibly set according to the situation that is embodied as and specific requirement.In this regard, the application is not construed as limiting.

S2：According to the mid frequency and pass band width, band filter is designed.

In one embodiment, frequency-division filter is carried out in order to design satisfactory band filter, specifically can be with：With Difference between the half of the mid frequency and the pass band width is as lowpass frequency, logical with described with the mid frequency The additive value of the half with bandwidth is used as high-pass equipment；According to the lowpass frequency and the high-pass equipment, the band logical is designed Wave filter.

In one embodiment, the wave filter is Butterworth band filter.Certainly can also be using others Wave filter.When being embodied as, suitable wave filter can be selected as the case may be.In this regard, the application is not construed as limiting.

S3：Frequency-division filter is carried out to the crest frequency attribute curve by the band filter, low frequency attribute is obtained Band logical frequency component curve and high frequency nature band logical frequency component curve.

In the present embodiment, it should be noted that the high frequency and the low frequency here are all relative concepts.Not Refer in particular to certain specific frequency values scope.With reference to present embodiment, can illustrate：When being obtained by wave filter frequency-division filter During to two band passband rate component curves, can be with the of a relatively high band passband rate component curve of wherein mid frequency as height The band logical frequency component curve of frequency attribute；Can be with the relatively low band passband rate component curve of wherein mid frequency as low The band logical frequency component curve of frequency attribute.For specific frequency range, the application is not limited thereto.

It should be noted that in the embodiment of the present application, first cycle refers specifically to rudimentary cycle, described second Cycle refers specifically to senior cycle.Accordingly, first depositional cycle curve refers specifically to rudimentary depositional cycle curve, Second depositional cycle curve refers specifically to senior depositional cycle curve.The first cycle rank refers specifically to rudimentary Cycle rank, the second cycle rank refer specifically to senior cycle rank.The first cycle threshold value refers specifically to low Level cycle threshold value, the second cycle threshold value refer specifically to senior cycle threshold value.

Step 104：Broken line process is carried out to the band logical frequency component curve of the low frequency attribute, the rotation of the first deposition is obtained Curve is returned, broken line process is carried out to the band logical frequency component curve of the high frequency nature, the second depositional cycle curve is obtained.

In one embodiment, it is bent in order to obtain first cycle of sedimentation according to the band logical frequency component curve of low frequency attribute Line, in order to obtain the second depositional cycle curve according to the band logical frequency component curve of high frequency nature, can be according to when being embodied as Following steps are performed：

S1：According to the band logical frequency component curve of the low frequency attribute, low frequency attribute most value sequence is obtained.

In one embodiment, in order to obtain low frequency attribute most value sequence, when being embodied as, Ke Yishi：Along described low The band logical frequency component curve of frequency attribute, every time from the point of the first predetermined number on the band logical frequency component curve of low frequency attribute Take out maximum and minima；Time coordinate corresponding to the maximum and minima of record taking-up；By taken out it is multigroup most Big value and minima are arranged according to the time sequencing that time coordinate is represented, you can obtain low frequency attribute most value sequence.In this enforcement In mode, it should be noted that the setting of the first predetermined number can be multiplied by 1000 institutes with 1 business obtained divided by mid frequency The numerical value for obtaining.For example, following point is sequentially arranged on the band logical frequency component curve of a low frequency attribute：1、3、 4、6、7、9、1、4、7、10、5、6、3、9、8、5、3、6、2、1.Assume calculated first predetermined number numerical value be 4, then often 4 points of secondary taking-up, select maximum and minima from 4 points.For the first time, maximum 6 is selected from 1,3,4,6 this 4 points With minima 1, and 6 and 1 selected corresponding time coordinate is recorded；Second, maximum 9 is selected from 7,9,1,4 this 4 points With minima 1, and 9 and 1 corresponding time coordinate here is recorded；For the third time, maximum is selected from 7,10,5,6 this 4 points 10 and minima 5, and record here 10 and 5 corresponding to time coordinate；4th time, select most from 3,9,8,5 this 4 points Big value 9 and minima 3, and record here 9 and 3 corresponding to time coordinate；5th time, select from 3,6,2,1 this 4 points Maximum 6 and minima 1, and record time coordinate here corresponding to 6 and 1.According to the time coordinate of each point selected, will Each maximum of points or minimum point selected is arranged according to the time sequencing that the corresponding time coordinate of each point is represented, Obtain low frequency attribute most value sequence：1、6、9、1、10、5、3、9、6、1.

S2：Numerical value process is carried out to the low frequency attribute most value sequence, the low frequency attribute most value sequence after being processed.

In one embodiment, in order to obtain the low frequency attribute most value sequence after satisfactory process, need to low Frequency attribute most value sequence carries out numerical value process, is embodied as to perform according to following steps：

S2-1：According to the cycle of sedimentation rank to be studied, the first cycle threshold value is determined.

In one embodiment, the first cycle threshold value is typically relatively large, for example, can be by the first cycle threshold value It is set to 1.Certainly, it is set to 1 here to be intended merely to the embodiment of the present application is better described.When being embodied as, may be arranged as With 1 other numerical value being close to, suitable numerical value can be set as the case may be used as the first cycle threshold value.In this regard, the application It is not construed as limiting.

S2-2：According to the low frequency attribute most value sequence, the meansigma methodss of low frequency attribute most value sequence are determined.

In one embodiment, such as low frequency attribute most value sequence is：1st, 6,9,1,10,5,3,9,6,1, to the sequence In row, each point is averaged as 5.1, then 5.1 be the meansigma methodss of the low frequency attribute most value sequence.

S2-3：According to putting in order for low frequency attribute most value sequence, respectively with two in the low frequency attribute most value sequence Adjacent most value processes subelement as low frequency, processes subelement to each low frequency in the low frequency attribute most value sequence respectively It is poor that two most values are carried out, and obtains the absolute value that corresponding low frequency processes the difference of subelement.

In one embodiment, in order to carry out low frequency attribute most value sequence numerical value process, need first to belong to according to low frequency Property most value sequence put in order, take two adjacent most values successively as a low frequency and process subelement, calculate low frequency respectively Process the absolute value of the difference of subelement.For example, for low frequency attribute most value sequence 1,6,9,1,10,5,3,9,6,1.Can be by According to putting in order for the low frequency attribute most value sequence, 1 and 6 are taken out successively and processes subelement as first low frequency, take out 9 and 1 Subelement is processed as second low frequency, 10 and 5 is taken out and subelement is processed as the 3rd low frequency, take out 3 and 9 as the 4th Low frequency processes subelement, takes out 6 and 1 and processes subelement as the 5th low frequency.Calculate each low frequency respectively to process in subelement The absolute value of the difference of two values, the absolute value for obtaining the difference that first low frequency processes subelement is 5, at second low frequency The absolute value of the difference of reason subelement is 8, and the absolute value of the difference of the 3rd low frequency process subelement is 5, at the 4th low frequency The absolute value of the difference of reason subelement is 6, and the absolute value of the difference of the 5th low frequency process subelement is 5.

S2-4：According to the low frequency process the absolute value of difference of subelement, the low frequency attribute most value sequence it is average Value and the first cycle threshold value, carry out numerical value process to the low frequency attribute most value sequence, the low frequency attribute after being processed Most value sequence.

In one embodiment, for the absolute value, described of subelement is processed according to the low frequency on request difference The meansigma methodss of low frequency attribute most value sequence and the first cycle threshold value, are carried out at numerical value to the low frequency attribute most value sequence Reason, specifically can be performed in accordance with the following methods：

The absolute value that the difference of subelement is processed in the low frequency is less than the first cycle threshold value, and the low frequency is processed Two in subelement be most worth in a larger meansigma methodss of the most value more than the low frequency attribute most value sequence, one less to be most worth In the case of meansigma methodss less than the low frequency attribute most value sequence, then retain in two most values that the low frequency is processed in sub- list The most value larger with the absolute value of the meansigma methodss difference of the high frequency nature most value sequence and corresponding time coordinate.

For example, for low frequency attribute most value sequence 1,6,9,1,10,5,3,9,6,1, through being calculated first It is 5 that low frequency processes the absolute value of the difference of subelement, and it is 8 that second low frequency processes the absolute value of the difference of subelement, the 3rd It is 5 that low frequency processes the absolute value of the difference of subelement, and it is 6 that the 4th low frequency processes the absolute value of the difference of subelement, the 5th It is 5 that low frequency processes the absolute value of the difference of subelement.Additionally, the meansigma methodss of the low frequency attribute most value sequence are 5.1, for convenience Illustrated, the first cycle threshold hypothesis are 6 herein.As the second, the 4th low frequency processes dividing for the difference of subelement 8 and 6 Wei not be all higher than equal to the first cycle threshold value, then each for retaining that the second, the 4th low frequency processed in subelement is most worth, i.e., Retain 9,1 and 3,9, and this 4 are most worth corresponding time coordinate.As the first, the 3rd, the 5th low frequency processes subelement The absolute value of difference is all 5, is most worth less than two in the first cycle threshold value 6, and above-mentioned each low frequency process subelement and accords with Meansigma methodss of one of them larger most value more than the low frequency attribute most value sequence are closed, a less most value belongs to less than the low frequency The situation of the meansigma methodss of property most value sequence.For example, first low frequency processes in subelement 1 and is less than 5.1, and 6 are more than 5.1, equally 3rd low frequency is processed in subelement, and 10 are more than 5.1, and 5 are less than 5.1, and the 5th low frequency is processed in subelement, and 6 are more than 5.1, And 1 is less than 5.1.Therefore, retain two in each sub- list of low frequency process respectively to above-mentioned each low frequency process subelement to be most worth In the most value larger with the absolute value of the meansigma methodss difference of the high frequency nature most value sequence and corresponding time coordinate.For example it is right In first low frequency processes subelement, the absolute value of most value 1 and 5.1 differences is 4.9, the most absolute value of the difference of value 6 and 5.1 For 0.9, it is more than 0.9 due to 4.9, therefore first low frequency processes subelement and retain most value 1 and corresponding time coordinate.In the same manner Three low frequencies process subelement and retain most value 10 and corresponding time coordinate, and the 5th low frequency processes subelement and retain most 1 He of value Corresponding time coordinate.Again being most worth for last reservation is arranged according to the time sequencing represented by corresponding time coordinate, is obtained Low frequency attribute most value sequence to after process：1、9、1、10、3、9、1.

S3：According to the low frequency attribute most value sequence after the process, first cycle of sedimentation is obtained by linear interpolation Curve.

In one embodiment, in order to carry out linear interpolation to the low frequency attribute most value sequence after process, and then obtain First cycle curve, specifically can perform according to following steps：

S3-1：According to putting in order for the low frequency attribute most value sequence after process, the low frequency attribute from after the process is most Choose continuous three to be most worth as a low frequency broken line treatment group in value sequence, obtain multiple low frequency broken line treatment groups, wherein, Two adjacent low frequency broken line treatment groups share one and are most worth；

S3-2：Each low frequency broken line treatment group in the plurality of low frequency broken line treatment group is held according to following steps respectively OK, obtain the low frequency broken line group after each process：

S3-2-1：Determine the maximum in three most values in the low frequency broken line treatment group, and by the maximum assignment For 100, the time coordinate of the maximum keeps constant；

S3-2-2：Determine the minima in three most values in the low frequency broken line treatment group, and the minima is assigned It is worth for 0, the time coordinate of the minima keeps constant；

S3-2-3：Determine the intermediate value in three most values in the low frequency broken line treatment group, and remove the intermediate value；

S3-3：According to the low frequency broken line group after handled, first depositional cycle curve is obtained by linear interpolation.

For example, the low frequency attribute most value sequence after process is：1、4、9、1、3、4、7、9、8.It is low after according to process Putting in order for frequency attribute most value sequence, takes out 3 in the low frequency attribute most value sequence from after process successively and is most worth to one group Low frequency broken line treatment group, takes out 1,4,9 as the first low frequency broken line treatment group, takes out 9,1,3 as the process of the second low frequency broken line Group, takes out 3,4,7 as the 3rd low frequency broken line treatment group, takes out 7,9,8 as the 4th low frequency broken line treatment group.Need explanation It is that two adjacent here low frequency broken line treatment groups share one and are most worth.For example, the first low frequency broken line treatment group and the second low frequency Broken line treatment group shares most value 9, and the second low frequency broken line treatment group shares most value 3 with the 3rd low frequency broken line treatment group.With reference to rule, Above-mentioned four groups of low frequencies broken line treatment group is handled as follows respectively：By the numerical value maximum 9 in the first low frequency broken line treatment group Be entered as 100, but retain 9 time coordinate, will former most value 9 be changed to 100, and 100 time coordinate or original most value 9 Time coordinate.It is similar, the numerical value minima in the first low frequency broken line treatment group is entered as into 1, and retains the time seat of former most value 1 It is denoted as time coordinate that this group of numerical value is 0.Additionally, in the first low frequency broken line treatment group, the numerical values recited of most value 4 neither Maximum in the group, nor the minima in the group, is the intermediate value in the group, therefore most value 4 is directly deleted.Obtain The first low frequency broken line treatment group after process is 0,100, wherein 0 and 100 according to represented by the original corresponding time coordinate being most worth Time sequencing arrangement.With reference to above-mentioned to the first low frequency broken line treatment group processing method, the second low frequency after can be processed Broken line treatment group is 100,0, and the 3rd low frequency broken line treatment group after process is 0,100, and the 4th low frequency broken line after process is processed Group is 0,100.And then can according to process after the first low frequency broken line group, process after the second low frequency broken line group, process after The 4th low frequency broken line treatment group after 3rd low frequency broken line treatment group and process, is processed to the low frequency broken line after above-mentioned each process Value in group after assignment carries out linear interpolation according to time coordinate position, you can obtain corresponding first depositional cycle curve.

S4：According to the band logical frequency component curve of the high frequency nature, high frequency nature most value sequence is obtained.

In one embodiment, in order to obtain high frequency nature most value sequence, when being embodied as, Ke Yishi：Along the height The band logical frequency component curve of frequency attribute, every time from the point of the first predetermined number on the band logical frequency component curve of high frequency nature Take out maximum and minima；Time coordinate corresponding to the maximum and minima of record taking-up；By taken out it is multigroup most Big value and minima are arranged according to the time sequencing that time coordinate is represented, you can obtain high frequency nature most value sequence.In this enforcement In mode, it should be noted that the setting of the first predetermined number can be multiplied by 1000 institutes with 1 business obtained divided by mid frequency The numerical value for obtaining.It is embodied as being referred to the acquisition of low frequency attribute most value sequence.Here, the application is repeated no more.

S5：Numerical value process is carried out to the high frequency nature most value sequence, the most value sequence of the high frequency nature after being processed.

In one embodiment, in order to obtain the most value sequence of the high frequency nature after satisfactory process, need to height Frequency attribute most value sequence carries out numerical value process, is embodied as to perform according to following steps：

S5-1：According to the cycle of sedimentation rank to be studied, the second cycle threshold value is determined.

In one embodiment, the second cycle threshold value is typically relatively small, for example, can be by the second cycle threshold value It is set to 0.5.Certainly, it is set to 0.5 here to be intended merely to the embodiment of the present application is better described.When being embodied as, can also set It is set to and 0.5 other numerical value being close to, suitable numerical value can be set as the case may be used as the second cycle threshold value.In this regard, The application is not construed as limiting.

S5-2：According to the high frequency nature most value sequence, the meansigma methodss of high frequency nature most value sequence are determined.

In the present embodiment, determine the method for the meansigma methodss of high frequency nature most value sequence and determine that low frequency attribute is most worth sequence The method of the meansigma methodss of row is similar to, and is specifically referred to the method for the meansigma methodss of determination low frequency attribute most value sequence.Here, this Shen Please repeat no more.

S5-3：According to putting in order for high frequency nature most value sequence, respectively with two in the high frequency nature most value sequence Adjacent most value as high-frequency therapeutic treatment subelement, respectively to each high-frequency therapeutic treatment subelement in the high frequency nature most value sequence It is poor that two most values are carried out, and obtains the absolute value of the difference of corresponding high-frequency therapeutic treatment subelement.

In one embodiment, in order to carry out high frequency nature most value sequence numerical value process, need first to belong to according to high frequency Putting in order for property most value sequence, takes two adjacent most values successively as a high-frequency therapeutic treatment subelement, calculates high frequency respectively Process the absolute value of the difference of subelement.Be embodied as it is similar with the absolute value for calculating the difference that low frequency processes subelement, can be with With reference to the process of the absolute value of the difference for calculating low frequency process subelement.Here, the application is repeated no more.

S5-4：Absolute value, the high frequency nature most value sequence according to the difference of the high-frequency therapeutic treatment subelement it is average Value and the second cycle threshold value, carry out numerical value process to the high frequency nature most value sequence, the high frequency nature after being processed Most value sequence.

In one embodiment, for absolute value on request according to the difference of the high-frequency therapeutic treatment subelement, described The meansigma methodss of high frequency nature most value sequence and the second cycle threshold value, are carried out at numerical value to the high frequency nature most value sequence Reason, specifically can be performed in accordance with the following methods：

The second cycle threshold value, and the high-frequency therapeutic treatment are less than in the absolute value of the difference of the high-frequency therapeutic treatment subelement Two in subelement be most worth in a larger meansigma methodss of the most value more than the high frequency nature most value sequence, one less to be most worth In the case of meansigma methodss less than the high frequency nature most value sequence, then retain in two most values that the low frequency is processed in sub- list The most value larger with the absolute value of the meansigma methodss difference of the high frequency nature most value sequence and corresponding time coordinate.

When being embodied as, the high frequency nature most value sequence is carried out numerical value process with to the low frequency attribute most value sequence Carry out numerical value and process similar, be referred to carry out the low frequency attribute most value sequence process of numerical value process.Here, the application Repeat no more.

S6：According to the most value sequence of the high frequency nature after the process, second cycle of sedimentation is obtained by linear interpolation Curve.

In one embodiment, in order to carry out linear interpolation to the most value sequence of the high frequency nature after process, and then obtain Second cycle curve, specifically can perform according to following steps：

S6-1：According to putting in order for the most value sequence of the high frequency nature after process, the high frequency nature from after the process is most Choose continuous three to be most worth as a high frequency broken line treatment group in value sequence, obtain multiple high frequency broken line treatment groups, wherein, Two adjacent high frequency broken line treatment groups share one and are most worth；

S6-2：Each high frequency broken line treatment group in the plurality of high frequency broken line treatment group is held according to following steps respectively OK, obtain the high frequency broken line group after each process：

S6-2-1：Determine the maximum in three most values in the high frequency broken line treatment group, and by the maximum assignment For 100, the time coordinate of the maximum keeps constant；

S6-2-2：Determine the minima in three most values in the high frequency broken line treatment group, and the minima is assigned It is worth for 0, the time coordinate of the minima keeps constant；

S6-2-3：Determine the intermediate value in three most values in the high frequency broken line treatment group, and remove the intermediate value；

S6-3：According to the high frequency broken line group after handled, second depositional cycle curve is obtained by linear interpolation.

When being embodied as, according to the most value sequence of the high frequency nature after the process, described second is obtained by linear interpolation The concrete processing procedure of depositional cycle curve, can be with to try to achieve the first depositional cycle curve according to low frequency attribute most value sequence similar With reference to the process for solving the first depositional cycle curve.Here, the application is repeated no more.

Step 105：According to first depositional cycle curve and second depositional cycle curve, carry out the cycle of sedimentation and draw Point.

In one embodiment, the result that can be divided according to the cycle of sedimentation, carries out geological prospecting.Wherein, described According to the result that the cycle of sedimentation is divided, geological prospecting is carried out, can be included：According to the result that the cycle of sedimentation is divided, to seismic facies and Sedimentary facies are analyzed, and analyze SEQUENCE STRATIGRAPHIC or predicting reservoir.Specifically, the result that can be divided according to the cycle of sedimentation, enters Row well site deployment.It should be noted that above-mentioned cited several utilization cycle of sedimentation division results carry out the mode of geological prospecting It is intended merely to the embodiment of the present application is better described, when being embodied as, it is also possible to as the case may be or construction requirement, utilizes The result that the cycle of sedimentation is divided carries out other corresponding geological prospecting work.In this regard, the application is not construed as limiting.

In the embodiment of the present application, scaling down processing is carried out to crest frequency attribute curve due to first passing through wave filter, then it is right Curve after frequency dividing carries out broken line process respectively, obtains corresponding depositional cycle curve, solves the existing cycle of sedimentation and divides When in method, depositional cycle curve is divided, the error that exists is big, efficiency is low, can only qualitative or semi-quantitative analyses technical problem, it is real Show and made full use of seismic data, technique effect that is quick, accurately dividing the cycle of sedimentation has reached the quantitative analyses cycle of sedimentation Purpose, and then deployment foundation can be provided for well site deployment.

In one embodiment, it is bent in order to obtain more accurate first depositional cycle curve and second cycle of sedimentation Line, the low frequency attribute most value sequence after according to the process, obtains first depositional cycle curve, root by linear interpolation It is according to the most value sequence of the high frequency nature after the process, after obtaining second depositional cycle curve by linear interpolation, described Method can also include：Repeat linear difference to process, multiple broken lineization is carried out to first depositional cycle curve and is processed, until The first cycle curve after process meets preset requirement, will meet the first depositional cycle curve of preset requirement as described first Depositional cycle curve, carries out multiple broken lineization and processes to second depositional cycle curve, and the second cycle after process is bent Line meets preset requirement, will meet the second depositional cycle curve of preset requirement as second depositional cycle curve.Further The accuracy of the first depositional cycle curve and the second depositional cycle curve for finally obtaining can be improved, error is further reduced.

In an embodiment, can further expand, can be according to specific requirement or condition of construction, to 3 grades of cycles Rank is analyzed.It is due to being that analysis object has 3 cycle ranks, corresponding to be needed to consider when design wave filter divides To 3 cycle ranks, during design wave filter, one-level cycle mid frequency and second-order cycle mid frequency are set respectively, according to one The corresponding dividing filter of level cycle mid frequency and second-order cycle center frequency designations institute, obtains 3 bands by the wave filter Pass filtered component curve, then carry out corresponding numerical value process and broken lineization process, tool to this 3 band pass filtered component curve respectively Body is referred to the above-mentioned process for obtaining the first depositional cycle curve according to low-frequency band pass filtered component curve and is processed, and then Corresponding 3 depositional cycle curves can be obtained, the analyzing and processing to 3 grades of cycle ranks is realized.Certainly for 1 grade of cycle rank Analyzing and processing, obtain 1 depositional cycle curve, be referred to 2 grades of cycle ranks of above-mentioned process, obtain 2 cycle of sedimentation songs 3 grades of cycle ranks of line and process, obtain the situation of 3 depositional cycle curves.Here, the application is repeated no more.

Based on same inventive concept, the device that a kind of cycle of sedimentation is divided in the embodiment of the present invention, is additionally provided, it is such as following Embodiment described in.Due to the principle of device solve problem it is similar to the method that the cycle of sedimentation is divided, therefore the cycle of sedimentation divide Device enforcement may refer to the cycle of sedimentation division method enforcement, repeat part repeat no more.It is used below, art Language " unit " or " module " can realize the combination of the software and/or hardware of predetermined function.Although described by following examples Device preferably with software realizing, but hardware, or the realization of the combination of software and hardware is also may and to be contemplated 's.Fig. 2 is referred to, is that a kind of the of device of the cycle of sedimentation division of the embodiment of the present invention constitutes structure chart, the device can be wrapped Include：First acquisition module 201, the second acquisition module 202, frequency-division filter module 203, broken line processing module 204 and deposition rotation Division module 205 is returned, below the structure is specifically described.

First acquisition module 201, for carrying out time frequency analysis to seismic data, obtains time-frequency data volume；

Second acquisition module 202, for according to the time-frequency data volume, obtaining crest frequency attribute curve；

Frequency-division filter module 203, for carrying out frequency-division filter to the crest frequency attribute curve, obtains low frequency attribute The band logical frequency component curve of band logical frequency component curve and high frequency nature；

Broken line processing module 204, for carrying out broken line process to the band logical frequency component curve of the low frequency attribute, The first depositional cycle curve is obtained, broken line process is carried out to the band logical frequency component curve of the high frequency nature, second is obtained Depositional cycle curve；

Cycle of sedimentation division module 205, for bent according to first depositional cycle curve and second cycle of sedimentation Line, carries out cycle of sedimentation division.

In one embodiment, in order to according to cycle rank, carry out frequency-division filter, obtain corresponding a plurality of bandpass filtering Component curve, the frequency-division filter module 204 can include：

In one embodiment, in order to carry out corresponding numerical value process and broken line respectively to band logical frequency component curve Process, obtain corresponding depositional cycle curve, the broken line processing module 205 can include：

In one embodiment, it is bent in order to obtain more accurate first depositional cycle curve and second cycle of sedimentation Line, after first depositional cycle curve and second depositional cycle curve are obtained by broken line processing module 205, First depositional cycle curve of broken line processing module 205 pairs can also be passed through and second depositional cycle curve enters respectively The multiple broken lineization of row is processed, and the second cycle curve and the first depositional cycle curve after process meets preset requirement, will Meet second depositional cycle curve and the first depositional cycle curve of preset requirement respectively as second depositional cycle curve With the first depositional cycle curve.

In an embodiment, can further expand, can be according to specific requirement or condition of construction, to 3 grades of cycles Rank is analyzed.Due to being that analysis object has 3 cycle ranks, the determination filter cell in frequency-division filter module 204 exists When performing concrete operations, it is contemplated that 3 cycle ranks, corresponding wave filter is determined：One-level cycle center frequency is set respectively first Rate and second-order cycle mid frequency, are divided accordingly further according to one-level cycle mid frequency and second-order cycle center frequency designations Frequency-division filter unit in wave filter, and then frequency-division filter module 204 can obtain 3 corresponding band logical filters by the wave filter Wave component curve, then this 3 band pass filtered component curve carries out corresponding numerical value process respectively by broken line processing module 205 With broken lineization process, described above is specifically referred to.Here, the application is repeated no more.Dividing certainly for 1 grade of cycle rank Analysis process, obtain 1 depositional cycle curve, be referred to 2 grades of cycle ranks of above-mentioned process, obtain 2 depositional cycle curves and 3 grades of cycle ranks are processed, the situation of 3 depositional cycle curves is obtained.Here, the application is equally also repeated no more.

Each embodiment in this specification is described by the way of progressive, identical similar portion between each embodiment Divide mutually referring to what each embodiment was stressed is the difference with other embodiment.Especially for system reality For applying example, as which is substantially similar to embodiment of the method, so description is fairly simple, related part is referring to embodiment of the method Part explanation.

It should be noted that system, device, module or unit that above-mentioned embodiment is illustrated, specifically can be by computer Chip or entity realize, or by the product with certain function realizing.For convenience of description, in this manual, retouch It is divided into various units with function when stating apparatus above to describe respectively.Certainly, can be the function of each unit when the application is implemented Realize in same or multiple softwares and/or hardware.

Additionally, in this manual, adjective as such as first and second is can be only used for an element or dynamic Work is made a distinction with another element or action, without requiring or implying any actual this relation or order.Permit in environment Perhaps, in the case of, in only element, part or step is should not be interpreted as limited to reference to element or part or step (s) It is individual, and can be one or more in element, part or step etc..

As can be seen from the above description, the embodiment of the present application provide the cycle of sedimentation divide method and apparatus by Scaling down processing is carried out to crest frequency attribute curve in wave filter is first passed through, then the curve after frequency dividing is carried out at broken line respectively Reason, obtains corresponding depositional cycle curve, thus solve the error existed when depositional cycle curve in existing method is divided it is big, Efficiency is low, can only qualitative or semi-quantitative analyses technical problem；As this method and device are to obtain peak value according to seismic data Frequency attribute curve, carries out follow-up process further according to crest frequency attribute curve, therefore is sufficiently used seismic data, solution Determine the technical problem that the area lacked in well-log information cannot be divided according to well-log information depositional cycle curve, reach extension The technique effect of the scope of application, and then can realize that the cycle of sedimentation is divided, foundation is provided for well site deployment；Again due to this method and Device is processed by carrying out corresponding numerical value respectively to the curve that frequency dividing is filtered, then carries out corresponding broken lineization process, is improved The accuracy that depositional cycle curve is divided, so realize it is quick, accurate determine deposition facies cycle curve, and then deposited The technique effect that cycle is divided.

In a specific embodiment, the method, device divided using the cycle of sedimentation that the application is provided enters to somewhere The research of the specific cycle of sedimentation of row.Specifically can perform in accordance with the following methods.Wherein, specific construction process can be refering to figure 3。

Please referring initially to Fig. 4, pending is somewhere domain actual time pure wave seismic data, wherein, sample rate is 4ms.

Step 1：Generalized S-transform time frequency analysis are carried out to the seismic data cube of time domain, is obtained based on generalized S-transform time-frequency The time-frequency domain data volume of analysis, instant frequency data volume.Wherein original seismic profile can be refering to shown in Fig. 4.Selected the earthquake In section, the single-channel seismic data of certain point carry out frequency division during generalized S-transform refering to shown in the left figure in Fig. 5 to the single track data Analysis result can be referred to shown in the middle graph in Fig. 5.

Step 2：To the single track time-frequency analysis result shown in the middle graph in Fig. 5, single-point amplitude spectrum is asked on time orientation most The corresponding frequency of big value, i.e. instantaneous spectrum crest frequency attribute.Concrete outcome can be referred to shown in the right figure in Fig. 5.

Step 3：According to the cycle of sedimentation rank (Ta~Tb) of required research, it is determined that the required depositional cycle curve for taking Mid frequency (F=1/T).Thus design band filter so that lowpass frequency and high-pass equipment are respectively：Fa=F-M and Fb=F+M (M typically takes 5~10)；Bandpass filtering is carried out to crest frequency attribute curve, the band passband of corresponding frequencies composition is obtained Rate component curve：The band logical frequency component curve CF2 of the band logical frequency component curve CF1 and high frequency nature of low frequency attribute.Specifically Fig. 6 can be referred to.

Step 4：Respectively the band logical frequency component curve of the band logical frequency component curve and high frequency nature of low frequency attribute is entered Row broken line, specifically can be bent to the band logical frequency component after bandpass filtering according to designed bandpass filter centre frequency F At interval of asking for maximum Amax and minima Amin respectively on T sampling point and record its corresponding time coordinate C_ on line CF Amax and C_Amin.According to the striked most value sequence of time coordinate order arrangement.Ask for meansigma methodss M_mean of most value sequence. According to cycle of sedimentation rank (Ta~Tb) threshold value to be studied, usually, the first cycle threshold value arranges larger as 1.0 is left The right side, the second cycle threshold value arrange less such as 0.5 or so.Two neighboring most value in most value sequence is made the difference, retains interpolation big In the most value sequence equal to threshold value.For most value of the difference less than threshold value does following process：If adjacent two most values are less than average Value M_mean, then retain less most value and corresponding time coordinate；, whereas if adjacent two most values are more than meansigma methodss M_ Mean, then retain larger most value and corresponding time coordinate；If additionally, adjacent two most values one are more than meansigma methodss M_ Mean, one is less than meansigma methodss M_mean, then retain the most value larger with meansigma methodss M_mean difference and corresponding time coordinate. To the most value sequence after above-mentioned process, continuous 3 are often selected as one group, wherein next group will include upper one group last One value, it is then determined that maximum and minima in these three values and maximum is entered as 100, minima is entered as 0. Linear interpolation is carried out according to coordinate position to the value after assignment, the required depositional cycle curve for taking just has been obtained：First deposition Cycle curve and the second depositional cycle curve.Concrete outcome can be refering to shown in Fig. 7.

The method, device divided by the cycle of sedimentation proposed using the application is regional without well according to the one kind in trial zone Poststack seismic data carries out the determination of depositional cycle curve.Due to being for without well area, so by treating drill site in work area Result of log interpretation divide verifying its cycle of sedimentation around the well point.By comparing, the heavy of the application offer is demonstrated The method and apparatus that product cycle is divided solves present in existing method depositional cycle curve and divides that error is big, efficiency is low, only The technical problem of the qualitative or semi-quantitative analyses of energy, realizes and makes full use of seismic data, quick, accurate to determine deposition facies cycle song The technique effect of line, has reached the purpose of the quantitative analyses cycle of sedimentation.

Although the method or apparatus that the different cycles of sedimentation are divided are mentioned in teachings herein, the application not office It is limited to be industry standard or the situation described by embodiment etc., some industry standards or uses self-defined mode or enforcement Example description practice processes on embodiment amended slightly can also realize above-described embodiment it is identical, equivalent or it is close, or Foreseeable implementation result after deformation.Using these modifications or the data acquisition after deformation, process, output, judgment mode etc. Embodiment, within the scope of still may belong to the optional embodiment of the application.

Although this application provides the method operating procedure as described in embodiment or flow chart, is based on routine or noinvasive The means of the property made can include more or less operating procedures.The step of enumerating in embodiment order is only numerous steps A kind of mode in execution sequence, does not represent unique execution sequence.When device or client production in practice is performed, can To perform or executed in parallel (such as at parallel processor or multithreading according to embodiment or method shown in the drawings order The environment of reason, even distributed data processing environment).Term " including ", "comprising" or its any other variant are intended to contain Lid nonexcludability is included, so that a series of process, method, product or equipment including key elements not only will including those Element, but also including other key elements being not expressly set out, or also include for this process, method, product or equipment Intrinsic key element.In the absence of more restrictions, it is not excluded that including the process of the key element, method, product or Also there are other identical or equivalent elements in person's equipment.

Device that above-described embodiment is illustrated or module etc., specifically can be realized by computer chip or entity, or by having There is the product of certain function to realize.For convenience of description, it is divided into various modules with function when describing apparatus above to retouch respectively State.Certainly, the function of each module can be realized in same or multiple softwares and/or hardware when the application is implemented, The module for realizing same function can be realized by the combination of multiple submodule etc..Device embodiment described above is only Schematically, for example, the division of the module, only a kind of division of logic function can have other drawing when actually realizing Point mode, such as multiple module or components can with reference to or be desirably integrated into another system, or some features can be ignored, Or do not perform.

It is also known in the art that in addition to realizing controller in pure computer readable program code mode, it is complete Entirely can by by method and step carry out programming in logic cause controller with gate, switch, special IC, may be programmed The form of logic controller and embedded microcontroller etc. is realizing identical function.Therefore this controller is considered one kind Hardware component, and the device for realizing various functions included to its inside can also be considered as the structure in hardware component.Or Person even, can be used for realizing that the device of various functions be considered as not only being the software module of implementation method but also being hardware Structure in part.

The application can be described in the general context of computer executable instructions, such as program Module.Usually, program module includes execution particular task or realizes the routine of particular abstract data type, program, object, group Part, data structure, class etc..The application is put into practice in a distributed computing environment can also, in these distributed computing environment, Task is performed by the remote processing devices connected by communication network.In a distributed computing environment, program module can With positioned at including in the local and remote computer-readable storage medium including storage device.

As seen through the above description of the embodiments, those skilled in the art can be understood that the application can By software plus required general hardware platform mode realizing.Based on such understanding, the technical scheme essence of the application On part that in other words prior art is contributed can be embodied in the form of software product, the computer software product Can be stored in storage medium, such as ROM/RAM, magnetic disc, CD etc., use so that a computer equipment including some instructions (can be personal computer, mobile terminal, server, or network equipment etc.) performs each embodiment of the application or enforcement Method described in some parts of example.

Each embodiment in this specification is described by the way of progressive, same or analogous portion between each embodiment Divide mutually referring to what each embodiment was stressed is the difference with other embodiment.The application can be used for crowd In more general or special purpose computing system environments or configuration.For example：Personal computer, server computer, handheld device or Portable set, laptop device, multicomputer system, set based on the system of microprocessor, set top box, programmable electronics Standby, network PC, minicomputer, mainframe computer, including the distributed computing environment etc. of any of the above system or equipment.

Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application have it is many deformation and Change is without deviating from spirit herein, it is desirable to which appended claim includes these deformations and changes without deviating from the application.

Claims

1. a kind of method that cycle of sedimentation is divided, it is characterised in that include：

Frequency-division filter is carried out to the crest frequency attribute curve, the band logical frequency component curve and high frequency category of low frequency attribute is obtained The band logical frequency component curve of property；

Broken line process is carried out to the band logical frequency component curve of the low frequency attribute, the first depositional cycle curve is obtained, to institute The band logical frequency component curve for stating high frequency nature carries out broken line process, obtains the second depositional cycle curve；

2. method according to claim 1, it is characterised in that frequency-division filter is carried out to the crest frequency attribute curve, The band logical frequency component curve of the band logical frequency component curve and high frequency nature of low frequency attribute is obtained, including：

According to the mid frequency and pass band width, band filter is designed；

Frequency-division filter is carried out to the crest frequency attribute curve by the band filter, the band passband of low frequency attribute is obtained The band logical frequency component curve of rate component curve and high frequency nature.

3. method according to claim 2, it is characterised in that according to the mid frequency and pass band width, designs band logical Wave filter, including：

Using the difference between the half of the mid frequency and the pass band width as lowpass frequency, with the mid frequency with The additive value of the half of the pass band width is used as high-pass equipment；

4. method according to claim 1, it is characterised in that the band logical frequency component curve of the low frequency attribute is carried out Broken lineization process, obtains the first depositional cycle curve, the band logical frequency component curve of the high frequency nature is carried out at broken line Reason, obtains the second depositional cycle curve, including：

According to the low frequency attribute most value sequence after the process, first depositional cycle curve is obtained by linear interpolation；

According to the most value sequence of the high frequency nature after the process, second depositional cycle curve is obtained by linear interpolation.

5. method according to claim 4, it is characterised in that numerical value process is carried out to the low frequency attribute most value sequence, Low frequency attribute most value sequence after being processed, including：

According to putting in order for low frequency attribute most value sequence, respectively with two adjacent most values in the low frequency attribute most value sequence Subelement is processed as low frequency, respectively each low frequency in the low frequency attribute most value sequence is processed two of subelement be most worth into Row is poor, obtains the absolute value that corresponding low frequency processes the difference of subelement；

The absolute value of difference of subelement, the meansigma methodss of the low frequency attribute most value sequence and described the are processed according to the low frequency One cycle threshold value, carries out numerical value process to the low frequency attribute most value sequence, the low frequency attribute most value sequence after being processed.

6. method according to claim 5, it is characterised in that processed according to the low frequency subelement difference it is absolute Value, the meansigma methodss of the low frequency attribute most value sequence and the first cycle threshold value, are carried out to the low frequency attribute most value sequence Numerical value process, including：

The absolute value of difference of subelement is processed in the low frequency more than or equal in the case of the first cycle threshold value, retains institute State low frequency and process two most values and the distinguished corresponding time coordinate of described two most values in subelement；

The absolute value that the difference of subelement is processed in the low frequency is less than the first cycle threshold value, and the low frequency processes son list In the case that two in unit are most worth both less than equal to the meansigma methodss of the low frequency attribute most value sequence, then retain at the low frequency Reason subelement in two most values in less most value and it is described it is less be most worth corresponding to time coordinate；

The absolute value that the difference of subelement is processed in the low frequency is less than the first cycle threshold value, and the low frequency processes son list In the case that two in unit are most worth the both greater than meansigma methodss of the low frequency attribute most value sequence, then retain the low frequency and process son Larger most value and described larger most it is worth corresponding time coordinate in two most values in unit；

The absolute value that the difference of subelement is processed in the low frequency is less than the first cycle threshold value, and the low frequency processes son list Two in unit be most worth in a larger meansigma methodss of the most value more than the low frequency attribute most value sequence, a less most value is less than In the case of the meansigma methodss of the low frequency attribute most value sequence, then retain in two most values that the low frequency is processed in sub- list with institute State the meansigma methodss difference of low frequency attribute most value sequence the larger most value of absolute value and this be most worth corresponding time coordinate.

7. method according to claim 4, it is characterised in that numerical value process is carried out to the high frequency nature most value sequence, High frequency nature after being processed most value sequence, including：

According to putting in order for high frequency nature most value sequence, respectively with two adjacent most values in the high frequency nature most value sequence As high-frequency therapeutic treatment subelement, respectively two of each high-frequency therapeutic treatment subelement in the high frequency nature most value sequence are most worth into Row is poor, obtains the absolute value of the difference of corresponding high-frequency therapeutic treatment subelement；

Absolute value, the meansigma methodss of the high frequency nature most value sequence according to the difference of the high-frequency therapeutic treatment subelement and described Two cycle threshold values, carry out numerical value process to the high frequency nature most value sequence, the most value sequence of the high frequency nature after being processed.

8. method according to claim 7, it is characterised in that according to the absolute of the difference of the high-frequency therapeutic treatment subelement Value, the meansigma methodss of the high frequency nature most value sequence and the second cycle threshold value, are carried out to the high frequency nature most value sequence Numerical value process, including：

In the case where the absolute value of the difference of the high-frequency therapeutic treatment subelement is more than or equal to the second cycle threshold value, then retain Two most values in the high-frequency therapeutic treatment subelement and described two most it is worth corresponding time coordinate；

The second cycle threshold value is less than in the absolute value of the difference of the high-frequency therapeutic treatment subelement, and high-frequency therapeutic treatment is single In the case that two in unit are most worth both less than equal to the meansigma methodss of the high frequency nature most value sequence, then retain the high frequency treatment Reason subelement in two most values in less most value and it is described it is less be most worth corresponding to time coordinate；

The second cycle threshold value is less than in the absolute value of the difference of the high-frequency therapeutic treatment subelement, and high-frequency therapeutic treatment is single In the case that two in unit are most worth the both greater than meansigma methodss of the high frequency nature most value sequence, then retain high-frequency therapeutic treatment Larger most value and described larger most it is worth corresponding time coordinate in two most values in unit；

The second cycle threshold value is less than in the absolute value of the difference of the high-frequency therapeutic treatment subelement, and high-frequency therapeutic treatment is single Two in unit be most worth in a larger meansigma methodss of the most value more than the high frequency nature most value sequence, a less most value is less than In the case of the meansigma methodss of the high frequency nature most value sequence, then retain in two most values in the sub- list of the high-frequency therapeutic treatment with institute State the meansigma methodss difference of high frequency nature most value sequence the larger most value of absolute value and this be most worth corresponding time coordinate.

9. method according to claim 1, it is characterised in that the band logical frequency component curve of the low frequency attribute is carried out Broken lineization process, obtains the first depositional cycle curve, the band logical frequency component curve of the high frequency nature is carried out at broken line Reason, after obtaining the second depositional cycle curve, methods described also includes：

Carry out multiple broken lineization to process to first depositional cycle curve, the first cycle after the process of multiple broken lineization is bent Line meets preset requirement, will meet the first depositional cycle curve of preset requirement as first depositional cycle curve；

Carry out multiple broken lineization to process to second depositional cycle curve, the second cycle after the process of multiple broken lineization is bent Line meets preset requirement, will meet the second depositional cycle curve of preset requirement as second depositional cycle curve.

10. the device that a kind of cycle of sedimentation is divided, it is characterised in that include：

Frequency-division filter module, for carrying out frequency-division filter to the crest frequency attribute curve, obtains the band passband of low frequency attribute The band logical frequency component curve of rate component curve and high frequency nature；

Broken line processing module, for carrying out broken line process to the band logical frequency component curve of the low frequency attribute, obtains One depositional cycle curve, carries out broken line process to the band logical frequency component curve of the high frequency nature, obtains the rotation of the second deposition Return curve；

Cycle of sedimentation division module, for according to first depositional cycle curve and second depositional cycle curve, carrying out The cycle of sedimentation is divided.

11. devices according to claim 10, it is characterised in that the frequency-division filter module includes：

Determine filter cell, for according to cycle of sedimentation rank to be studied, determining mid frequency, and according to center frequency Rate and pass band width, determine band filter；

Frequency-division filter unit, for carrying out frequency-division filter to the crest frequency attribute curve by the band filter, obtains To the band logical frequency component curve of the band logical frequency component curve and high frequency nature of low frequency attribute.

12. devices according to claim 10, it is characterised in that the broken line processing module includes：

Low frequency most value sequence processing unit, for carrying out numerical value process to the low frequency attribute most value sequence, after being processed Low frequency attribute most value sequence；

Low frequency broken line processing unit, for according to the low frequency attribute most value sequence after the process, being obtained by linear interpolation First depositional cycle curve；

High frequency most value sequence acquiring unit, for the band logical frequency component curve according to the high frequency nature, obtains high frequency nature Most value sequence；

High frequency most value sequence processing unit, for carrying out numerical value process to the high frequency nature most value sequence, after being processed High frequency nature most value sequence；

High frequency broken line processing unit, for according to the most value sequence of the high frequency nature after the process, being obtained by linear interpolation Second depositional cycle curve.