CN103412332A - Method for quantitative calculation of thickness of thin reservoir layer - Google Patents
Method for quantitative calculation of thickness of thin reservoir layer Download PDFInfo
- Publication number
- CN103412332A CN103412332A CN2013100324415A CN201310032441A CN103412332A CN 103412332 A CN103412332 A CN 103412332A CN 2013100324415 A CN2013100324415 A CN 2013100324415A CN 201310032441 A CN201310032441 A CN 201310032441A CN 103412332 A CN103412332 A CN 103412332A
- Authority
- CN
- China
- Prior art keywords
- seismic
- target
- data
- reservoir
- thickness
- 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
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a method for quantitative calculation of thickness of a thin reservoir layer. The method comprises the steps of carrying out three-dimensional multiple coverage data acquisition on the earth's surface to obtain seismic data; carrying out "three-high" processing on the seismic data to obtain three-dimensional post-stack seismic data; conducting logging-seismic calibration by use of logging data, and making clear of the position of a target reservoir layer; carrying out three-dimensional seismic horizon interpretation on the target layer to obtain a time window extraction guide layer; under the guidance of the time window extraction guide layer, giving an appropriate time window according to estimated thickness of the target reservoir layer and seismic data dominant frequency; extracting relative spectral peak frequency increment seismic attributes, and calculating two-way travel time thickness of the target thin reservoir layer; and with P-wave velocity data, converting the time thickness of the target thin reservoir layer into strata thickness. Compared with a current quantitative prediction method for a thin reservoir layer, the method is not subject to an absolute reflection coefficient of top and bottom sides of the thin reservoir layer and does not rely on well data, and the prediction accuracy is high.
Description
Technical field
The present invention relates to the method for seismic prospecting in applied geophysics, particularly about thin reservoirs exploration technical field, is that a kind of relative spectrum peak frequency increment that utilizes quantitatively calculates the method for thin reservoir thickness specifically.
Background technology
Method of seismic prospecting in applied geophysics refers on earth's surface with the manual method earthquake-wave-exciting, to underground propagation the time, meet the different rock stratum interphase of medium character, reflection and transmission will occur in seismic event, on earth's surface, receive seismic event with wave detector, obtain seismologic record.Feature on seismologic record is relevant with the character of subterranean strata and structure, by seismologic record being processed and being explained, can infer the character and form of subterranean strata.Aspect precision and accuracy reconnoitring, seismic prospecting all is better than other geophysical exploration method, thereby in oil-gas exploration, occupies very important status.
But nearly decades are along with oil-gas exploration deepens continuously, the increase of difficulties in exploration, structural trapping hydrocarbon-bearing pool sweeping, that be easy to find is excavated totally, and our exploration targets has to turn to thin reservoir, lithology and complex oil and gas reservoir.Thin reservoir is important one of type of preserving of China's Continental Petroliferous Basins, is the main goal in research of present stage petroleum exploration in China.Be limited by earthquake-capturing, process the impact of links, conventional reservoir prediction means resolution is inadequate, does not reach the effect of thin reservoir being carried out to detailed predicting, with actual production require gap larger.Thereby new thin (mutually) layer EARTHQUAKE QUANTIFICATION Forecasting Methodology of research and development has important theoretical and practical significance.
For thin RESERVOIR RECOGNITION and thickness quantitative examination, many scholars did discussion both at home and abroad.1973, Widess points out when thickness of thin layer is less than λ/8, the maximal value of thickness of thin layer and reflection configuration amplitude is linear, and provided thickness equations, but he has only discussed reflection coefficient equal and opposite in direction at the bottom of top, opposite polarity thin layer, and computing formula needs the absolute size of known reflection coefficient.Koefoed and Voogd (1981), Robertson and Nogami (1984) etc. also did similar research, had drawn similar conclusion.Chung and Lawton (1995) improve the approximate formula of Widess (1973), have improved the computational accuracy of superthin layer.Sun Lu equality (2009,2010) crest frequency and thin layer time thickness relationship are discussed, proposed a kind of thinking of utilizing crest frequency quantitatively to calculate thin layer time thickness, do not obtained velocity of longitudinal wave and thin layer time thickness is converted to the method for zone thickness but provide.Partyka and Gridley etc. (1999) propose the spectral factorization method, utilize sunken cycle of frequency and the thin layer time thickness of thin layer frequency spectrum to exist reciprocal relation to carry out thickness of thin layer quantitative estimation (Partyka (1999,2005); Marfurt and Kirlin (2001); Huang Xude (2002) etc.), but the method is subjected to the impact of seismic band width, and the frequency of the inadequate time-frequency spectrum of the earthquake frequency range cycle of falling into is difficult to determine.On the basis of spectral factorization method, the research team of Castagna (2006) has proposed the spectrum inversion technique, adopts prior imformation and Spectral Decomposition Technique to improve thin layer imaging effect (Chopra and Castagna etc., 2007 that are less than tuning thickness; Puryear and Castagna etc., 2008), the method precision is higher, but requires very high to the signal to noise ratio (S/N ratio) of wavelet accuracy and data.In addition, Zeng Hongliu (2005) proposes 90 ° of phase conversion technique, more is conducive to the explanation to thin layer; Dou Yisheng (1995) proposes spectral amplitude duplicate ratio standard measure and explains thickness of thin layer.Huang Zhenping etc. (1997), Zhang Enke etc. (2006) adopt the non-linear inversion algorithm predicts thickness of thin layer of optimizing, and these class methods need to have the thin-layer sample participation in learning of known thickness, can only be used in the area of well.In sum, current thin reservoir quantitative forecasting technique or need in advance the absolute size of known thin reservoir top, the reflection coefficient at the end or the acquisition process of seismic data is had to harsh requirement, or strictly depend on drilling well, thereby these methods are very restricted when practical application.
Summary of the invention
The objective of the invention is, the method for the thin reservoir thickness of a kind of quantitative calculating is provided.A kind of new seismic properties has been proposed---relative spectrum peak frequency increment, by the thin reservoir formation thickness of this property calculation, the method is not limited by the restriction of the reflection coefficient absolute size of bottom surface, thin reservoir top, applied widely and computational accuracy is higher.
Content of the present invention comprises:
Gather three dimensional seismic data, described seismic data is carried out to " three height " (high-fidelity, high s/n ratio, high resolving power) and process the three-dimensional post-stack seismic data of acquisition;
Utilize well logging and 3D seismic data to carry out well-shake composite traces and demarcate, dark relation while determining, hard objectives reservoir position, utilize statistical method to extract the statistics seismic wavelet on the geological data within comprising the target reservoir section;
In the target interval, carrying out 3-D seismics layer position explains, during acquisition, window extracts guide layer, the time window extract under the guidance of guide layer, extract the geological data daughter that comprises the target reservoir section, the time geological data in window that comprises the target reservoir section is carried out to the extraction of relative spectrum peak frequency increment seismic properties;
Utilize relative spectrum peak frequency increment seismic properties quantitatively to calculate the two-way travel time thickness of the thin reservoir of target;
Well-log information is carried out to statistical study, obtain the two relation of p-wave impedance and velocity of longitudinal wave in the target interval;
3D seismic data is carried out to the inverting of poststack p-wave impedance, obtains three-dimensional p-wave impedance data volume, so the time window extract under the guidance of guide layer, obtain the p-wave impedance data daughter that comprises the target reservoir section;
Extract target interval p-wave impedance plane seismic properties, utilize the two relation of target interval p-wave impedance and velocity of longitudinal wave, target interval p-wave impedance attribute is converted to the velocity of longitudinal wave attribute, obtain the velocity of longitudinal wave plane attribute of target interval;
Utilize the velocity of longitudinal wave plane attribute of target interval, the two-way travel time thickness of the thin reservoir of target is converted to zone thickness.
Beneficial effect of the present invention is, the method of the thin reservoir thickness of calculating that the present invention proposes, be not subjected to the restriction of thin layer top end reflection coefficient absolute size, also without the value of known reflection coefficient, computation process is not controlled by well simultaneously, in the area that the well data is few, also can apply, and have very high computational accuracy.
The accompanying drawing explanation
In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, below will the accompanying drawing of required use in embodiment or description of the Prior Art be briefly described, obviously, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skills, under the prerequisite of not paying creative work, can also obtain according to these accompanying drawings other accompanying drawing.
Fig. 1 is that the relative spectrum peak frequency increment that utilizes of the embodiment of the present invention quantitatively calculates the method flow diagram of thin reservoir thickness;
Fig. 2 is the relative spectrum peak frequency increment seismic properties figure on certain Jurassic systerm top, work area of the embodiment of the present invention;
Fig. 3 is the thickness planimetric map of the sand body time by relative spectrum peak frequency increment property calculation of the embodiment of the present invention;
Fig. 4 is the velocity of longitudinal wave planimetric map on certain Jurassic systerm top, work area of the embodiment of the present invention;
Fig. 5 is the real thickness chart of calculating of certain Jurassic systerm top, work area sand body of the embodiment of the present invention;
Fig. 6 is the wedge-like thin film model schematic diagram of the embodiment of the present invention;
Fig. 7 is the wedge-like thin film model theogram schematic diagram of the embodiment of the present invention;
Fig. 8 is the spectral amplitude schematic diagram in each road of wedge-like thin film model theogram of the embodiment of the present invention;
Fig. 9 is the relative spectrum peak frequency increment of the embodiment of the present invention and the matched curve schematic diagram of thickness of thin layer;
Figure 10 is the computational accuracy of the embodiment of the present invention and forefathers' computational accuracy comparative analysis schematic diagram;
Embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is clearly and completely described, obviously, described embodiment is only the present invention's part embodiment, rather than whole embodiment.Based on the embodiment in the present invention, those of ordinary skills, not making under the creative work prerequisite the every other embodiment obtained, belong to the scope of protection of the invention.
Fig. 1 is that the relative spectrum peak frequency increment attribute that utilizes of the embodiment of the present invention quantitatively calculates the method flow diagram of thin reservoir thickness.As shown in Figure 1, the method comprises the steps:
S101, on earth's surface, gather the seismic prospecting of multi-fold high-resolution three-dimension and obtain three dimensional seismic data.Particularly, in step S101, comprise and carry out the design of multi-fold stereo observing system, according to the recording geometry design, construct, in earth's surface specific location (focal point), carry out exciting of man-made explosion (as explosive), earth's surface specific location (acceptance point) utilization, receive instrument (as three-component seismometer) and record seismic data.
S102, the three dimensional seismic data that collection is obtained carry out " three height " (high-fidelity, high s/n ratio, high resolving power) processing, to obtain three-dimensional post-stack seismic data.Particularly, in step S102, comprise denoising, static correction, energy compensating, deconvolution and pre-stack time migration.
S103, utilization well logging and 3D seismic data are logged well-the synthetic seismic record demarcation, dark relation while determining, hard objectives reservoir present position.Particularly, comprise the log data of collecting in work area, utilize interval transit time and density logging Curves compilation p-wave impedance curve; On seismologic record, extract seismic wavelet W (t), by seismic wavelet and p-wave impedance curve plotting theogram road; Synthetic traces and seismic trace near well are carried out to optimum matching, and dark relation while determining, clearly mark the reservoir present position.
S104, carry out 3-D seismics layer position and explain, window extracts guide layer t when obtaining one.The time window extract guide layer t and can be end face, the bottom surface of target interval or all can in the target interval.In a preferred embodiment of the present invention, when described window extract guide layer t preferably want can the thin reservoir of indicating target in the temporal information at diverse location place.
S105, the time window extract under the guidance of guide layer t, extract the geological data daughter S (t) that comprises the target reservoir section.Particularly, according to given window width T when suitable of the dominant frequency of estimating thickness and seismic data of target reservoir, if the time window extracts that guide layer t explains is the end face of target reservoir, the time window that extracts the geological data daughter is [t, t+T]; The time window extracts that guide layer t explains is the bottom surface of target reservoir, the time window that extracts the geological data daughter is [t-T, t]; If the time window extracts that guide layer t explains is the centre of target reservoir, the time window that extracts the geological data daughter is [t-T/2, t+T/2].
In a preferred embodiment of the present invention, we find that there is following rule: when target reservoir estimate large or dominant frequency seismic data of thickness when on the low side, the time window T should select bigger; Otherwise, the time window T should be less.
S106, upper at the geological data daughter S (t) that comprises the target reservoir section, carry out relative spectrum peak frequency increment seismic properties K to each seismic trace
RPFI, jExtraction, K
RPFI, jComputing formula as follows:
Wherein: f
WpMeet
F
Sp, jMeet
K
RPFI, jBe the relative spectrum peak frequency increment of j road seismic trace, n comprises the number of seismic trace, S in geological data daughter S (t)
j(t) be earthquake record in j road in geological data daughter S (t), A
w(f) be the spectral amplitude of seismic wavelet W (t), A
s(f) be that S is recorded in the earthquake of j road
j(t) spectral amplitude, f
WpFor the spectrum peak frequency of seismic wavelet W (t), f
Sp, jBe that S is recorded in the earthquake of j road
j(t) spectrum peak frequency.
The relative spectrum peak frequency increment seismic properties K extracted
RPFI, jAs shown in Figure 2, in figure, the arrow indication is typical fluvial facies deposit feature.
S107, utilize relative spectrum peak frequency increment seismic properties K
RPFI, jQuantitatively calculate the two-way travel time thickness d of the thin reservoir of target
j.When
The time, adopt formula
Calculate; When
The time, adopt formula
Calculate.Wherein, d
jFor the two-way travel time thickness of the thin reservoir of each seismic trace position target that calculates, K
RPFI, jBe the relative spectrum peak frequency increment of j road seismic trace, R
1For reflection coefficient corresponding to thin reservoir interface, top, R
2For reflection coefficient corresponding to thin reservoir bottom boundary.
Calculate the two-way travel time thickness d of the thin reservoir of target
jAs shown in Figure 3.The channel deposit feature is clear, and can reflect obviously that center, river course stratum time thickness is large, the feature reduced gradually to both sides.
So far (step S101-S107), the two-way travel time thickness d of the thin reservoir of target
jCalculated acquisition, if expect the actual (real) thickness of the thin reservoir of target, needed the velocity information of target reservoir section, the embodiment of the present invention also provides a kind of method (step S108-S112) that obtains velocity information.
S109, well-log information is carried out to statistical study, to obtain the two relation of target interval p-wave impedance and velocity of longitudinal wave.Concrete operation step comprises: 1, collect the log data that comprises objective interval in work area, mainly comprise each mouthful well: densimetric curve Den, interval transit time curve Dt, lithology curve (carrying out the achievement of lithologic interpretation by logging trace).2, utilize interval transit time curve Dt to calculate velocity of longitudinal wave curve Vp, Vp=1/Dt, by the synthetic p-wave impedance curve Zp of densimetric curve Den, velocity of longitudinal wave curve Vp, Zp=Den*Vp.3, divide the interior p-wave impedance Zp of drafting objective interval of lithology and the X plot of velocity of longitudinal wave Vp, the X plot transverse axis is p-wave impedance Zp, and the longitudinal axis is velocity of longitudinal wave Vp.4, the two quantitative relationship of the p-wave impedance Zp of match sandstone and velocity of longitudinal wave Vp on X plot, obtain being calculated by p-wave impedance Zp the formula of velocity of longitudinal wave Vp, is designated as Vp=f (Zp).
S110, the inverting of three-dimensional poststack p-wave impedance, obtain three-dimensional p-wave impedance data volume.According to the work area actual conditions, inversion method adopts constraint Sparse Pulse Inversion, recurrence inversion or Application of Logging-constrained Inversion all can.In a preferred embodiment of the present invention, adopted the constraint Sparse Pulse Inversion.
S111, the time window extract under the guidance of guide layer, obtain the p-wave impedance data daughter that comprises the target reservoir section, and extract target interval p-wave impedance plane attribute Z
j.The time window determine that mode is identical with the invention process step S105.
S112, calculating target interval velocity of longitudinal wave plane attribute.By the quantitative relationship Vp=f (Zp) of the determined p-wave impedance Zp of the invention process step S109 and velocity of longitudinal wave Vp and the target interval p-wave impedance plane attribute Z of the invention process step S111 extraction
jCalculate target interval velocity of longitudinal wave plane attribute V
j, obtain each seismic trace position velocity of longitudinal wave V
j.The target interval velocity of longitudinal wave V calculated
jThe plane attribute as shown in Figure 4.
The two-way travel time thickness d of S113, the thin reservoir of target
jBe converted to zone thickness H
j.The two-way travel time thickness d of the thin reservoir of target obtained by the invention process step S107
jAnd the velocity of longitudinal wave V provided by the invention process step S112
jInformation, calculate the zone thickness H of the thin reservoir of target
j, computing formula is: H
j=(d
j* V
j)/2.Finally, as shown in Figure 5, the thickness of sand body, between 4-30m, meets the existing geological knowledge in this area to the zone thickness result of the thin reservoir of target, well portrayed the planar distribution feature of river channel sand, the favo(u)rable target of finding next step for this area provides good foundation.
In order to further illustrate in the invention process step S107 by relative spectrum peak frequency increment K
RPFI, jCalculate the two-way travel time thickness d
jThe origin of computing formula, provide following explanation (Fig. 6-Figure 10 describes by wedge shape thin film model and seismologic record thereof).
Fig. 6 a, b be respectively the rhythm type of law (
) and the alternation type (
) two kinds of sphenoid thin film models, more directly perceived for the size that makes reflection coefficient, do not consider the variation of density here.Fig. 7 is the theogram that wedge shape thin film model and wavelet convolution obtain, and seismic wavelet has adopted spectrum peak frequency f
WpThe Ricker wavelet of=40Hz.The position of crest-trough that in Fig. 7, dotted line representative is picked up on composite traces, as seen because thickness of thin layer is less than tuning thickness, can not from synthetic waveform peak-top two, end reflecting interface is distinguished in the paddy position.Fig. 8 is the spectral amplitude of wedge shape thin layer seismologic record, and on spectral amplitude, solid vertical line means the position at wavelet spectrum peak; With the line of arrow dotted line for each channel amplitude spectrum spectrum peak, the direction that direction of arrow indication thickness of thin layer reduces.Visible, along with reducing of thickness of thin layer, in two kinds of situations of the rhythm type of law and alternation type, all along with the trend that reduces to be increase of thickness of thin layer, the present invention has utilized this feature to carry out the quantitative calculating of thickness of thin layer to the spectrum peak frequency in each road of composite traces just.
In reality, the wavelet spectrum peak frequency is more difficultly accurately to ask for out, in order to remove wavelet spectrum peak frequency f
WpOn the impact of computation process, example of the present invention is not used the absolute value of seismic trace spectrum peak frequency, but proposes to use relative spectrum peak frequency increment K
RPFICarry out thickness of thin layer calculating.Fig. 9 be to the rhythm type of law (
) and the alternation type (
) in two kinds of situations, relative spectrum peak frequency increment K
RPFIWith thickness of thin layer, concern fitting result, parabolic mathematical form f (x)=ax is adopted in match
2+ c carries out, and the longitudinal axis is relative spectrum peak frequency increment, and transverse axis is thickness of thin layer.In Fig. 9, loose point is actual data point, and solid line is matched curve.Visible, the fitting precision in two kinds of situations is all very high, and fitting precision R-square all reaches more than 99.5%.By Fig. 9, obtained, the rhythm type of law (
) fitting coefficient a=-2.33, c=0.224 in situation; The alternation type (
) fitting coefficient a=-5.12, c=-0.00522 in situation, obtain thus using relative spectrum peak frequency increment K
RPFICarry out the formula of thickness of thin layer calculating, as follows:
When
The time,
When
The time,
Finally, Figure 10 provides the comparison diagram of actuarial precision, and in figure, four lines represent respectively actual value, (1973) method that proposes of utilizing Widess, utilize the method that Chung (1995) proposes and the result of calculation of utilizing method provided by the invention.Visible, result of calculation of the present invention and actual value are more approaching, and precision is higher, and the method for Widess (1973) and Chung (1995) all needs known thin layer top end reflection R
1, R
2Absolute value, and method of the present invention does not rely on R
1, R
2Size, the scope of application is wider.
In sum, useful achievement of the present invention and advantage are: proposed a kind of relative spectrum peak frequency increment seismic properties of utilizing and quantitatively calculated the method for thin reservoir formation thickness.With current existing thin reservoir quantitative forecasting technique, compare, the method is not limited by the reflection coefficient absolute size of bottom surface, thin reservoir top, does not rely on the well data, and applied widely, and precision of prediction is high.
Claims (10)
1. method of quantitatively calculating thin reservoir thickness, characteristics are to comprise the steps:
(1) on earth's surface, gather multi-fold high-resolution three-dimension seismic exploration data data;
(2) seismic data is carried out to " three height " (high-fidelity, high s/n ratio, high resolving power) and process, to obtain three-dimensional post-stack seismic data;
(3) utilize well logging and 3D seismic data to log well-the synthetic seismic record demarcation, dark relation while determining, hard objectives reservoir present position;
(4) well-log information is carried out to statistical study, to obtain the two relation of target interval p-wave impedance and velocity of longitudinal wave;
(5) in the target interval, carry out 3-D seismics layer position and explain, window extracts guide layer when obtaining one;
(6) the time window extract under the guidance of guide layer, according to the given suitable time window of the dominant frequency of estimating thickness and seismic data of target reservoir, to extract the geological data daughter that comprises the target reservoir section;
(7) the time geological data in window that comprises the target reservoir section is carried out to the extraction of relative spectrum peak frequency increment seismic properties;
(8) utilize relative spectrum peak frequency increment seismic properties quantitatively to calculate the two-way travel time thickness of the thin reservoir of target;
(9) 3D seismic data step (2) obtained carries out the inverting of poststack p-wave impedance, to obtain three-dimensional p-wave impedance data volume;
(10) the time window that utilizes step (6) to obtain, the three-dimensional p-wave impedance data volume that step (9) is obtained extracts, and with acquisition, comprises the p-wave impedance data daughter of target reservoir section;
(11) seismic attributes analysis is carried out in the p-wave impedance data daughter that comprises target reservoir section step (10) obtained, and extracts target interval p-wave impedance plane seismic properties;
(12) utilize the two relation of target interval p-wave impedance that step (5) obtains and velocity of longitudinal wave, and the target interval p-wave impedance plane seismic properties of step (11) acquisition, target interval p-wave impedance attribute is converted to the velocity of longitudinal wave attribute, to obtain the velocity of longitudinal wave plane attribute of target interval;
(13) utilize the velocity of longitudinal wave plane attribute of the target interval of step (12) acquisition, the two-way travel time thickness of the thin reservoir of target of step (8) calculating is converted to zone thickness.
2. method according to claim 1, characteristics are the described multi-fold high-resolution three-dimension of step (1) seismic exploration data data, when gathering in the wild, the design view examining system is the acquisition system of multi-fold, utilize the 3-component earthquake detector collection signal, utilize the magnetic recording original earthquake data.
3. method according to claim 1, characteristics are that step (2) described " three height " (high-fidelity, high s/n ratio, high resolving power) is processed, and " three height " requires to relate to the links that data are processed.
4. method according to claim 1, to be that step (4) is described carry out statistical study to well-log information and comprise characteristics: at objective interval, p-wave impedance and velocity of longitudinal wave are carried out to histogram, crossplot analysis and its, and in the two mathematical relation of X plot upper returning target interval p-wave impedance and velocity of longitudinal wave, if there are lithology curve or lithology information in work area, to mainly utilize reservoir sandstone lithology section returning when p-wave impedance and velocity of longitudinal wave concern.
5. method according to claim 1, characteristics be step (5) when described window extract guide layer t and can be end face, the bottom surface of target interval or all can in the target interval, if the time window extracts that guide layer t explains is the end face of target reservoir, the time window that extracts the geological data daughter is [t, t+T]; If the time window extracts that guide layer t explains is the bottom surface of target reservoir, the time window that extracts the geological data daughter is [t-T, t]; If the time window extracts that guide layer t explains is the centre of target reservoir, the time window that extracts the geological data daughter is [t-T/2, t+T/2].
6. method according to claim 1, characteristics are given window T when suitable of the described dominant frequency of estimating thickness and seismic data according to target reservoir of step (6), when target reservoir, estimate large or dominant frequency seismic data of thickness when on the low side, window when while selecting, window T is large as far as possible, otherwise during selection, window T is a hour window.
7. method according to claim 1, characteristics are the described relative spectrum of step (7) peak frequency increment seismic properties, it is defined as:
Wherein: f
WpMeet
F
Sp, jMeet
In formula (1):
K
RPF, jBe the relative spectrum peak frequency increment of j road seismic trace, i.e. j road seismic trace S
j(t) spectrum peak frequency that spectrum peak frequency deducts seismic wavelet W (t) is again divided by the spectrum peak frequency of seismic wavelet W (t);
N is Seismic Traces S
j(t) number;
A
w(f) be the spectral amplitude of seismic wavelet W (t);
A
w(f) be j road seismic trace S
j(t) spectral amplitude;
F
WpSpectrum peak frequency for seismic wavelet W (t);
F
Sp, jBe j road seismic trace S
j(t) spectrum peak frequency.
8. method according to claim 8, the preparation method that is characterized in seismic wavelet W (t) is from the geological data daughter that comprises the target reservoir section, being extracted by statistical method; A
w(f) be the spectral amplitude of seismic wavelet W (t), obtain by time domain seismic wavelet W (t) is carried out to Short Time Fourier Transform; A
s(f) be Seismic Traces S
j(t) spectral amplitude, by each the seismic trace S in the geological data daughter that comprises the target reservoir section that step (6) is obtained
j(t) carrying out Short Time Fourier Transform obtains.
9. method according to claim 1, characteristics are the described two-way travel time thickness that utilizes relative spectrum peak frequency increment quantitatively to calculate the thin reservoir of target of step (8), are undertaken by following computing formula:
When
The time,
Formula (2-1)
When
The time,
Formula (2-2)
Formula (2-1) and (2-2) in:
D
jTwo-way travel time thickness for the thin reservoir of j road seismic trace position target that calculates;
K
RPFI, jBe the relative spectrum peak frequency increment of j road seismic trace, its account form as claimed in claim 8;
R
1For reflection coefficient corresponding to thin reservoir interface, top;
R
2For reflection coefficient corresponding to thin reservoir bottom boundary.
10. method according to claim 10 is characterized in: when
The time, representing the thin reservoir of target top, bottom boundary reflection coefficient equal and opposite in direction, opposite in sign, be the thin reservoir of the constant amplitude rhythm type of law; When
The time, representing the thin reservoir of target top, bottom boundary reflection coefficient equal and opposite in direction, symbol is identical, is the thin reservoir of constant amplitude alternation type.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310032441.5A CN103412332B (en) | 2013-01-22 | 2013-01-22 | A kind of method of determining thin reservoir thickness |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310032441.5A CN103412332B (en) | 2013-01-22 | 2013-01-22 | A kind of method of determining thin reservoir thickness |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103412332A true CN103412332A (en) | 2013-11-27 |
| CN103412332B CN103412332B (en) | 2016-05-25 |
Family
ID=49605359
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310032441.5A Expired - Fee Related CN103412332B (en) | 2013-01-22 | 2013-01-22 | A kind of method of determining thin reservoir thickness |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103412332B (en) |
Cited By (26)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103941286A (en) * | 2014-04-16 | 2014-07-23 | 中国海洋石油总公司 | Rapid three-dimensional horizon interpretation method based on horizontal navigation |
| CN104375177A (en) * | 2014-11-03 | 2015-02-25 | 中国石油天然气股份有限公司 | Reservoir top surface identification method and device by utilizing inversion data volume |
| CN105093306A (en) * | 2014-05-15 | 2015-11-25 | 中国石油化工股份有限公司 | Method for automatic interpretation and thickness calculation of reservoir in geophysical exploration |
| CN105204066A (en) * | 2015-10-26 | 2015-12-30 | 中国矿业大学 | Spectral decomposition-based method for directly indicating intrusion position of igneous rock into coal seam |
| CN105467449A (en) * | 2014-09-04 | 2016-04-06 | 中国石油化工股份有限公司 | Deep thin interbedded reservoir quantitative characterization method based on seismic grading sensitive attribute fusion |
| CN106255902A (en) * | 2014-05-08 | 2016-12-21 | 雪佛龙美国公司 | Use the system and method for seismic data analysis geologic feature |
| CN106338767A (en) * | 2015-07-07 | 2017-01-18 | 中国石油化工股份有限公司 | Multi-parameter seismic attribute analysis method and system |
| CN106405647A (en) * | 2016-12-22 | 2017-02-15 | 成都晶石石油科技有限公司 | Tune inversion method about sedimentary stratum thickness |
| CN106483562A (en) * | 2016-09-28 | 2017-03-08 | 中国海洋石油总公司 | A kind of method based on seismic amplitude predictably layer thickness variation |
| CN106707338A (en) * | 2016-11-18 | 2017-05-24 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Reservoir high-precision prediction method under strong shielding |
| CN106772598A (en) * | 2016-12-12 | 2017-05-31 | 中国石油大学(华东) | Using the method for receiver function periodic measurement sedimentary formation time thickness |
| CN106842314A (en) * | 2016-12-09 | 2017-06-13 | 中国石油天然气股份有限公司 | Method for determining formation thickness |
| CN106873036A (en) * | 2017-04-28 | 2017-06-20 | 中国石油集团川庆钻探工程有限公司 | Denoising method based on well-seismic combination |
| CN107450100A (en) * | 2016-06-01 | 2017-12-08 | 中国石油化工股份有限公司 | Thin reservoir thickness prediction method and device |
| CN107688197A (en) * | 2017-07-12 | 2018-02-13 | 中国石油天然气集团公司 | Thin bed prediction method and apparatus |
| CN109188520A (en) * | 2018-09-17 | 2019-01-11 | 中国石油天然气股份有限公司 | Thin reservoir thickness prediction method and device |
| CN109283577A (en) * | 2017-07-20 | 2019-01-29 | 中国石油化工股份有限公司 | A kind of seismic layer labeling method |
| CN109470187A (en) * | 2018-11-12 | 2019-03-15 | 中国海洋石油集团有限公司 | Reservoir thickness prediction method based on three attribute of earthquake |
| CN109633750A (en) * | 2019-01-18 | 2019-04-16 | 中国地质大学（北京） | Nonlinear mapping relationship inversion method based on well logging phase wave impedance and seismic waveform |
| CN110456415A (en) * | 2019-07-17 | 2019-11-15 | 中国石油大港油田勘探开发研究院 | A kind of thin RESERVOIR INTERPRETATION method and system of terrestrial facies based on local peaking's frequency |
| CN111158048A (en) * | 2020-01-04 | 2020-05-15 | 杨林海 | Analysis method for improving reservoir prediction precision through seismic waveform envelope interpretation |
| CN111221038A (en) * | 2018-11-26 | 2020-06-02 | 中国石油天然气股份有限公司 | Method and device for quantitatively predicting thickness of thin reservoir |
| CN112946754A (en) * | 2019-12-10 | 2021-06-11 | 中国石油天然气集团有限公司 | Reservoir porosity prediction method and device |
| CN114002743A (en) * | 2021-10-20 | 2022-02-01 | 中国海洋石油集团有限公司 | Method and device for calculating reservoir thickness of natural gas hydrate |
| CN114152985A (en) * | 2021-12-14 | 2022-03-08 | 中国地质大学（北京） | Method for determining boundary of underground ancient river channel and thickness of thin sand body in boundary |
| CN114428322A (en) * | 2020-10-29 | 2022-05-03 | 中国石油天然气股份有限公司 | Method and device for predicting thickness of thin reservoir based on frequency attribute |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6243650B1 (en) * | 1998-09-11 | 2001-06-05 | Diamond Geoscience Research Corporation | Method of determining net reservoir thickness |
| US20080294345A1 (en) * | 2007-05-22 | 2008-11-27 | Chevron U.S.A. Inc. | Method for determining attributes associated with net-sand thickness |
| CN101408624A (en) * | 2007-10-08 | 2009-04-15 | 陶庆学 | Forecasting and evaluating technologies of three-dimensional earthquake optimum time window river course sand body storage layer |
| US20100165791A1 (en) * | 2008-09-05 | 2010-07-01 | Statoilhydro Asa | Method for quantitatively making a thickness estimate of thin geological layers based on seismic reflection signals in the frequency domain |
| CN102109611A (en) * | 2009-12-23 | 2011-06-29 | 中国石油天然气集团公司 | Fast and convenient method for predicting high-quality petroleum reservoir in virtue of seism attributes |
| CN102109613A (en) * | 2009-12-23 | 2011-06-29 | 中国石油天然气股份有限公司 | Method for determining effective thickness of target reservoir under complex geological conditions |
| CN102736107A (en) * | 2011-04-07 | 2012-10-17 | 中国石油天然气股份有限公司 | Energy constraint heterogeneous reservoir thickness identification system |
-
2013
- 2013-01-22 CN CN201310032441.5A patent/CN103412332B/en not_active Expired - Fee Related
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6243650B1 (en) * | 1998-09-11 | 2001-06-05 | Diamond Geoscience Research Corporation | Method of determining net reservoir thickness |
| US20080294345A1 (en) * | 2007-05-22 | 2008-11-27 | Chevron U.S.A. Inc. | Method for determining attributes associated with net-sand thickness |
| CN101408624A (en) * | 2007-10-08 | 2009-04-15 | 陶庆学 | Forecasting and evaluating technologies of three-dimensional earthquake optimum time window river course sand body storage layer |
| US20100165791A1 (en) * | 2008-09-05 | 2010-07-01 | Statoilhydro Asa | Method for quantitatively making a thickness estimate of thin geological layers based on seismic reflection signals in the frequency domain |
| CN102109611A (en) * | 2009-12-23 | 2011-06-29 | 中国石油天然气集团公司 | Fast and convenient method for predicting high-quality petroleum reservoir in virtue of seism attributes |
| CN102109613A (en) * | 2009-12-23 | 2011-06-29 | 中国石油天然气股份有限公司 | Method for determining effective thickness of target reservoir under complex geological conditions |
| CN102736107A (en) * | 2011-04-07 | 2012-10-17 | 中国石油天然气股份有限公司 | Energy constraint heterogeneous reservoir thickness identification system |
Cited By (37)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103941286B (en) * | 2014-04-16 | 2016-08-17 | 中国海洋石油总公司 | A kind of quick three-dimensional layer position based on horizontal navigation means of interpretation |
| CN103941286A (en) * | 2014-04-16 | 2014-07-23 | 中国海洋石油总公司 | Rapid three-dimensional horizon interpretation method based on horizontal navigation |
| CN106255902B (en) * | 2014-05-08 | 2019-01-04 | 雪佛龙美国公司 | Use the system and method for seismic data analysis geologic feature |
| CN106255902A (en) * | 2014-05-08 | 2016-12-21 | 雪佛龙美国公司 | Use the system and method for seismic data analysis geologic feature |
| CN105093306A (en) * | 2014-05-15 | 2015-11-25 | 中国石油化工股份有限公司 | Method for automatic interpretation and thickness calculation of reservoir in geophysical exploration |
| CN105467449B (en) * | 2014-09-04 | 2018-01-05 | 中国石油化工股份有限公司 | Deep layer Thin oil sandwich quantitatively characterizing method based on earthquake classification Sensitive Attributes fusion |
| CN105467449A (en) * | 2014-09-04 | 2016-04-06 | 中国石油化工股份有限公司 | Deep thin interbedded reservoir quantitative characterization method based on seismic grading sensitive attribute fusion |
| CN104375177A (en) * | 2014-11-03 | 2015-02-25 | 中国石油天然气股份有限公司 | Reservoir top surface identification method and device by utilizing inversion data volume |
| CN106338767A (en) * | 2015-07-07 | 2017-01-18 | 中国石油化工股份有限公司 | Multi-parameter seismic attribute analysis method and system |
| CN105204066A (en) * | 2015-10-26 | 2015-12-30 | 中国矿业大学 | Spectral decomposition-based method for directly indicating intrusion position of igneous rock into coal seam |
| CN105204066B (en) * | 2015-10-26 | 2017-07-21 | 中国矿业大学 | A kind of direct indicating means in coal seam Igneous rock invasion position based on spectral factorization |
| CN107450100A (en) * | 2016-06-01 | 2017-12-08 | 中国石油化工股份有限公司 | Thin reservoir thickness prediction method and device |
| CN107450100B (en) * | 2016-06-01 | 2019-04-02 | 中国石油化工股份有限公司 | Thin reservoir thickness prediction method and device |
| CN106483562A (en) * | 2016-09-28 | 2017-03-08 | 中国海洋石油总公司 | A kind of method based on seismic amplitude predictably layer thickness variation |
| CN106707338A (en) * | 2016-11-18 | 2017-05-24 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Reservoir high-precision prediction method under strong shielding |
| CN106842314B (en) * | 2016-12-09 | 2019-09-03 | 中国石油天然气股份有限公司 | Method for determining formation thickness |
| CN106842314A (en) * | 2016-12-09 | 2017-06-13 | 中国石油天然气股份有限公司 | Method for determining formation thickness |
| CN106772598A (en) * | 2016-12-12 | 2017-05-31 | 中国石油大学(华东) | Using the method for receiver function periodic measurement sedimentary formation time thickness |
| CN106405647B (en) * | 2016-12-22 | 2018-04-13 | 成都晶石石油科技有限公司 | A kind of tuning inversion method on sedimentary formation thickness |
| CN106405647A (en) * | 2016-12-22 | 2017-02-15 | 成都晶石石油科技有限公司 | Tune inversion method about sedimentary stratum thickness |
| CN106873036A (en) * | 2017-04-28 | 2017-06-20 | 中国石油集团川庆钻探工程有限公司 | Denoising method based on well-seismic combination |
| CN107688197A (en) * | 2017-07-12 | 2018-02-13 | 中国石油天然气集团公司 | Thin bed prediction method and apparatus |
| CN109283577B (en) * | 2017-07-20 | 2023-02-28 | 中国石油化工股份有限公司 | Seismic horizon calibration method |
| CN109283577A (en) * | 2017-07-20 | 2019-01-29 | 中国石油化工股份有限公司 | A kind of seismic layer labeling method |
| CN109188520B (en) * | 2018-09-17 | 2020-05-08 | 中国石油天然气股份有限公司 | Thin reservoir thickness prediction method and device |
| CN109188520A (en) * | 2018-09-17 | 2019-01-11 | 中国石油天然气股份有限公司 | Thin reservoir thickness prediction method and device |
| CN109470187A (en) * | 2018-11-12 | 2019-03-15 | 中国海洋石油集团有限公司 | Reservoir thickness prediction method based on three attribute of earthquake |
| CN111221038A (en) * | 2018-11-26 | 2020-06-02 | 中国石油天然气股份有限公司 | Method and device for quantitatively predicting thickness of thin reservoir |
| CN111221038B (en) * | 2018-11-26 | 2022-02-01 | 中国石油天然气股份有限公司 | Method and device for quantitatively predicting thickness of thin reservoir |
| CN109633750A (en) * | 2019-01-18 | 2019-04-16 | 中国地质大学（北京） | Nonlinear mapping relationship inversion method based on well logging phase wave impedance and seismic waveform |
| CN110456415A (en) * | 2019-07-17 | 2019-11-15 | 中国石油大港油田勘探开发研究院 | A kind of thin RESERVOIR INTERPRETATION method and system of terrestrial facies based on local peaking's frequency |
| CN112946754A (en) * | 2019-12-10 | 2021-06-11 | 中国石油天然气集团有限公司 | Reservoir porosity prediction method and device |
| CN112946754B (en) * | 2019-12-10 | 2024-03-01 | 中国石油天然气集团有限公司 | Reservoir porosity prediction method and device |
| CN111158048A (en) * | 2020-01-04 | 2020-05-15 | 杨林海 | Analysis method for improving reservoir prediction precision through seismic waveform envelope interpretation |
| CN114428322A (en) * | 2020-10-29 | 2022-05-03 | 中国石油天然气股份有限公司 | Method and device for predicting thickness of thin reservoir based on frequency attribute |
| CN114002743A (en) * | 2021-10-20 | 2022-02-01 | 中国海洋石油集团有限公司 | Method and device for calculating reservoir thickness of natural gas hydrate |
| CN114152985A (en) * | 2021-12-14 | 2022-03-08 | 中国地质大学（北京） | Method for determining boundary of underground ancient river channel and thickness of thin sand body in boundary |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103412332B (en) | 2016-05-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103412332B (en) | A kind of method of determining thin reservoir thickness | |
| Mi et al. | Estimating near-surface shear-wave-velocity structures via multichannel analysis of Rayleigh and love waves: An experiment at the Boise hydrogeophysical research site | |
| EP3488267B1 (en) | Seismic spectral balancing | |
| CN102109616B (en) | Wave impedance inversion method for deposition cycle restriction | |
| CN101329407B (en) | Method for quick switching wave direct simulation to determine formation lithology and lithofacies change | |
| CN104237945B (en) | A kind of seismic data self adaptation high resolution processing method | |
| CN104570067B (en) | Phased seismic inversion method in a kind of geophysical exploration | |
| CA2818790C (en) | Seismic trace attribute | |
| CN109541685B (en) | River channel sand body identification method | |
| CN104516018A (en) | Porosity inversion method under lithological constraint in geophysical exploration | |
| CN104155693A (en) | Angle gather seismic response numerical computation method of reservoir fluid fluidity | |
| CN111722284B (en) | Method for establishing speed depth model based on gather data | |
| CN103257363A (en) | Method for detecting inclination angle of fissure in underground fissure type reservoir stratum | |
| Khan et al. | 3D structural modeling integrated with seismic attribute and petrophysical evaluation for hydrocarbon prospecting at the Dhulian Oilfield, Pakistan | |
| Khalil et al. | Engineering applications of seismic refraction method: A field example at Wadi Wardan, Northeast Gulf of Suez, Sinai, Egypt | |
| CN102385066B (en) | Pre-stack earthquake quantitative imaging method | |
| CN102565852B (en) | Angle domain pre-stack offset data processing method aiming to detect oil-gas-bearing property of reservoir | |
| Yan et al. | Building subsurface models with horizon-guided interpolation and deep learning: Application to the Volve field data | |
| AlHakeem | 3D seismic attribute analysis and machine learning for reservoir characterization in Taranaki Basin, New Zealand | |
| CN104991271A (en) | Prestack depth migration method based on VTI medium bidirectional wave field superposition | |
| CN109613615B (en) | Geologic body scale quantitative estimation method based on pre-stack seismic response analysis | |
| CN116068663A (en) | Igneous rock wave impedance inversion method based on magnetic vibration combined low-frequency modeling | |
| CN106226814A (en) | Utilize converted shear wave seismic data inversion reservoir S-wave impedance and the method for density | |
| Ivanov et al. | Refraction nonuniqueness studies at levee sites using the refraction-tomography and JARS methods | |
| Hammed et al. | THREE DIMENSIONAL SEISMIC, WELL LOG AND STRUCTURAL ANALYSIS OF" IGBOBI" FIELD, OFFSHORE NIGER DELTA. |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160525 Termination date: 20180122 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |