CN102967878B - Method and device for predicting oil and gas reservoirs through passive seismic wave field - Google Patents

Method and device for predicting oil and gas reservoirs through passive seismic wave field Download PDF

Info

Publication number
CN102967878B
CN102967878B CN201210443023.0A CN201210443023A CN102967878B CN 102967878 B CN102967878 B CN 102967878B CN 201210443023 A CN201210443023 A CN 201210443023A CN 102967878 B CN102967878 B CN 102967878B
Authority
CN
China
Prior art keywords
information
passiveseismic
frequency
frequency field
field
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201210443023.0A
Other languages
Chinese (zh)
Other versions
CN102967878A (en
Inventor
王孝
王述江
高建虎
王宇超
王小卫
胡自多
曾华会
寇龙江
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Natural Gas Co Ltd
Original Assignee
China Petroleum and Natural Gas Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Petroleum and Natural Gas Co Ltd filed Critical China Petroleum and Natural Gas Co Ltd
Priority to CN201210443023.0A priority Critical patent/CN102967878B/en
Publication of CN102967878A publication Critical patent/CN102967878A/en
Application granted granted Critical
Publication of CN102967878B publication Critical patent/CN102967878B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

An embodiment of the invention relates to the field of oil and gas reservoir prediction and particularly relates to a method and device for predicting oil and gas reservoirs through a passive seismic wave field. The method includes: time domain passive seismic information is obtained, and the time domain passive seismic information is converted into frequency domain passive seismic information; the frequency domain passive seismic information is sorted to a common detector gather to perform overlying treatment on frequency domains, and frequency domain overlying information is obtained; the frequency domain overlying information is displayed as a frequency domain overlying section according to the sequence of detection point number; and abnormal energy range is recognized from the frequency domain overlying section so as to predict the oil and gas reservoirs. Abnormal information of the passive seismic wave field is used for determining positions of the oil and gas reservoirs, no artificial stimulation is needed, a large amount of cost and environmental cost can be saved, and the method and device can further be used for reflecting complex exploration areas with failure seismic exploration.

Description

A kind of method of passiveseismic wave field predicting hydrocarbon reservoirs and device
Technical field
The present invention relates to prediction of oil-gas reserve field, particularly relate to a kind of method and device of passiveseismic wave field predicting hydrocarbon reservoirs.
Background technology
Common seismic exploration engineering, through long-run development, has had comparatively proven technique, from initial two-dimensional seismic survey, three-dimensional exploration, to multi-wave and multi-component, the exploration of high density wide-azimuth three-dimensional of newly-developed.
But, the focus of common seismic needs artificial excitation, the exploration of a block at least needs 1 to two years from design, field acquisition to data process, reservoir prediction, oil and gas detection, well site deployment whole cycle, used device is many, its human resources and cost consumption are quite huge, bring huge pressure to environment simultaneously.And the focus of common seismic is generally on earth's surface, source location is known, its exploration idea on the basis of subsurface formations accurate imaging, carries out the analysis of many attributes multidisciplinary synthesis come look for oil and gas Tibetan, therefore the development of conventional surveying techniques, but exploration cycle is more and more faced long, cost is high, the challenge that environmental pressure is large.
To sum up, find the oil-gas exploration technical method that a kind of efficiency is higher, more economical, environment is more friendly to be imminent.
Summary of the invention
The embodiment of the present invention provides a kind of method and device of passiveseismic wave field predicting hydrocarbon reservoirs, artificial excitation is needed for solving common seismic exploration engineering in prior art, and it is long to face exploration cycle, cost is high, the problem that its human resources and cost consumption are quite huge.
The method of a kind of passiveseismic wave field predicting hydrocarbon reservoirs for solving the problems of the technologies described above that the embodiment of the present invention provides, wherein, comprising:
Acquisition time territory passiveseismic information;
Convert described time domain passiveseismic information to frequency field passiveseismic information;
Described frequency field passiveseismic information is sorted into the overlap-add procedure that common detector gather does frequency field, obtains frequency field overlapped information;
Described frequency field overlapped information is shown as frequency field stacked section according to the order of detection period;
From the scope of described frequency field stacked section identification energy exception, and then predicting hydrocarbon reservoirs.
The method of above-mentioned passiveseismic wave field predicting hydrocarbon reservoirs, wherein, described convert described time domain passiveseismic information to frequency field passiveseismic information after, described method also comprises:
Described frequency field passiveseismic information is calculated the root mean square energy of point frequency segment information in described frequency field, and draw the root mean square energy trace of point frequency segment information;
From the root mean square energy trace of described point of frequency segment information, identify energy abnormal ranges, and then verify described hydrocarbon-bearing pool.
The method of above-mentioned passiveseismic wave field predicting hydrocarbon reservoirs, wherein, described acquisition time territory passiveseismic information comprises: from initial passiveseismic information, filtering is lower than very first time length value and the information exceeding first frequency value; Wherein: described very first time length is not higher than 500ms, and described first frequency value is not less than 40Hz.
The method of above-mentioned passiveseismic wave field predicting hydrocarbon reservoirs, wherein, according to formula:
S stk ( ω ) = S nor m 1 ( ω ) + S norm 2 ( ω ) + S norm 3 ( ω ) + . . . + S norm n ( ω ) n
The overlap-add procedure of frequency field is done at described common detector gather; Wherein, S stk(ω) be described frequency field overlapped information at described common detector gather, for the seismic trace of common detector gather, 1,2,3 ... n is the number of channels of common detector gather.
The method of above-mentioned passiveseismic wave field predicting hydrocarbon reservoirs, wherein, according to formula:
A rms = 1 L ( Σ S stk ( ω i ) 2 ) 1 2 Calculate the root mean square energy of described point frequency segment information; Wherein, A rmsfor the root mean square energy value of point frequency segment information, L is number of samples in point frequency band, S stki) be the energy value of frequency field passiveseismic information at i sampling point place.
The embodiment of the present invention also provides a kind of device of passiveseismic wave field predicting hydrocarbon reservoirs, comprising:
Information acquisition unit, for acquisition time territory passiveseismic information;
Converting unit, for converting described time domain passiveseismic information to frequency field passiveseismic information;
Processing unit, for described frequency field passiveseismic information is sorted into the overlap-add procedure that common detector gather does frequency field, obtains frequency field overlapped information;
Section display unit, for being shown as frequency field stacked section by described frequency field overlapped information according to the order of detection period;
Predicting unit, for identifying the scope of energy exception from described frequency field stacked section, and then predicting hydrocarbon reservoirs.
The device of above-mentioned passiveseismic wave field predicting hydrocarbon reservoirs, wherein, described device also comprises: storage unit, for described frequency field passiveseismic information to be calculated the root mean square energy of point frequency segment information in described frequency field, and draws the root mean square energy trace of point frequency segment information;
Described predicting unit also for identifying energy abnormal ranges from the root mean square energy trace of described point of frequency segment information, and then verifies described hydrocarbon-bearing pool.
The device of above-mentioned passiveseismic wave field predicting hydrocarbon reservoirs, wherein, described device also comprises: filter element, for filtering from initial passiveseismic information lower than very first time length value and the information exceeding first frequency value; Wherein: described very first time length is not higher than 500ms, and described first frequency value is not less than 40Hz.
The device of above-mentioned passiveseismic wave field predicting hydrocarbon reservoirs, wherein, described processing unit is according to formula:
S stk ( ω ) = S nor m 1 ( ω ) + S norm 2 ( ω ) + S norm 3 ( ω ) + . . . + S norm n ( ω ) n
The overlap-add procedure doing frequency field at described common detector gather obtains frequency field overlapped information, wherein, and S stk(ω) be the information after the superposition of common detector gather frequency field, for the seismic trace of common detector gather, 1,2,3 ... n is the number of channels of common detector gather.
The device of above-mentioned passiveseismic wave field predicting hydrocarbon reservoirs, wherein, described device also comprises: computing unit, for according to formula:
A rms = 1 L ( Σ S stk ( ω i ) 2 ) 1 2
Calculate the root mean square energy of described point frequency segment information; Wherein, A rmsfor the root mean square energy of point frequency segment information, L is number of samples in point frequency band, S stki) be the energy value of frequency field passiveseismic information at i sampling point place.
By method and the device of the passiveseismic wave field predicting hydrocarbon reservoirs of the embodiment of the present invention, do not need artificial excitation, and passiveseismic is the position utilizing various abnormal information directly to determine hydrocarbon-bearing pool, therefore can save a large amount of manpowers, time, financial cost and environmental cost.
Accompanying drawing explanation
Accompanying drawing described herein is used to provide a further understanding of the present invention, forms a application's part, does not form limitation of the invention.In the accompanying drawings:
The method flow diagram of a kind of passiveseismic wave field predicting hydrocarbon reservoirs that Fig. 1 provides for the embodiment of the present invention;
The source book information schematic diagram that in the method for a kind of passiveseismic wave field predicting hydrocarbon reservoirs that Fig. 2 provides for the embodiment of the present invention, common seismic gathers;
Initial passiveseismic information schematic diagram in the method for a kind of passiveseismic wave field predicting hydrocarbon reservoirs that Fig. 3 provides for the embodiment of the present invention;
Time domain passiveseismic information schematic diagram in the method for a kind of passiveseismic wave field predicting hydrocarbon reservoirs that Fig. 4 provides for the embodiment of the present invention;
The method medium frequency territory passiveseismic information of a kind of passiveseismic wave field predicting hydrocarbon reservoirs that Fig. 5 provides for the embodiment of the present invention is in the frequency field diagrammatic cross-section of common detector gather single-coverage;
The frequency field stacked section schematic diagram of method medium frequency territory passiveseismic information after common detector gather superposition of a kind of passiveseismic wave field predicting hydrocarbon reservoirs that Fig. 6 provides for the embodiment of the present invention;
The method medium frequency territory passiveseismic information of a kind of passiveseismic wave field predicting hydrocarbon reservoirs that Fig. 7 provides for the embodiment of the present invention is at the root mean square energy trace of point frequency segment information of frequency field;
The device schematic diagram of a kind of passiveseismic wave field predicting hydrocarbon reservoirs that Fig. 8 provides for the embodiment of the present invention.
Embodiment
For making the object, technical solutions and advantages of the present invention clearly, below in conjunction with accompanying drawing, embodiment of the present invention is described further in detail.At this, schematic description and description of the present invention is for explaining the present invention, but not as a limitation of the invention.
The method flow diagram of a kind of passiveseismic wave field predicting hydrocarbon reservoirs that Fig. 1 provides for the embodiment of the present invention, comprising:
Step 101, acquisition time territory passiveseismic information; Concrete, from initial passiveseismic information, filtering is lower than very first time length value and the information exceeding first frequency value; Wherein: described very first time length is not higher than 500ms, described first frequency value is not less than 40Hz, as shown in Figure 3, be the initial passiveseismic information in the present embodiment, shown in Fig. 4 for the initial passiveseismic information shown in Fig. 3 is carried out the time domain passiveseismic information after filtering.Preferably, wherein very first time length is 500ms, first frequency value is 40Hz, because time span there will be low-frequency information distortion lower than the short passiveseismic information of time domain of 500ms when being transformed into frequency field, equally, frequency is serious more than the usual interfering noise of the high-frequency information of 40Hz, and they all may produce illusion, therefore first filtering is needed, with the impact making the time domain passiveseismic information obtained eliminate distortion and noise.
In one embodiment of the invention, initial passiveseismic information can obtain from the source book information of common seismic collection; As shown in Figure 2, be the source book information schematic diagram that common seismic gathers, the region wherein marking A place is initial passiveseismic information region, and the region of mark B is artificial earthquake information area; In the source book information that the common seismic shown in Fig. 2 gathers, filter out artificial earthquake information, only retain man-made explosion primary wave time passiveseismic information in the past as initial passiveseismic information; Subsequently, lower than the information of very first time length and the information exceeding first frequency value in the initial passiveseismic information of filtering, time domain passiveseismic information can be obtained.
In another embodiment of the invention, initial passive information directly can also receive passiveseismic source book from wave detector and read, because the passiveseismic information utilizing wave detector directly to receive only comprises initial passiveseismic information, do not comprise artificial earthquake information, can facilitate, directly obtain initial passiveseismic information.
Step 102, converts described time domain passiveseismic information to frequency field passiveseismic information; Thus carry out follow-up frequency field overlap-add procedure and at frequency field curve plotting.Preferably, can according to formula:
S ( ω , τ ) = ∫ - ∞ ∞ s ( t ) g ( t - τ ) e - jωτ dt
Convert described time domain passiveseismic information to frequency field passiveseismic information, wherein, S (ω, τ) is frequency field passiveseismic information, and s (t) is time domain passiveseismic information, and g (t) is window function.Concrete, if the limited passiveseismic information of time domain is s (t), s (t) comprises two parts content usually: time domain passiveseismic wave field and noise information, frequency field passiveseismic information is S (ω, τ), utilize Short Time Fourier Transform formula, passiveseismic wave field limited for time domain and noise information s (t) are transformed to frequency field passiveseismic information S (ω, τ); Wherein g (t) is window function, the effective length of cutting passiveseismic s (t) is gone in time domain window function g (t), Fourier transform is made to cutting down the local signal come, because s (t) message length and window function are determined, S (ω, τ) can be write a Chinese character in simplified form into S (ω).The above-mentioned frequency field that time domain converted to is because identical wave field passiveseismic information frequency characteristic of field is similar, random noise does not then have the feature of similar frequency field, therefore in follow-up process, just passiveseismic information signal can be strengthened by superposition, compacting noise.
Step 103, is sorted into the frequency field passiveseismic information obtained in step 102 overlap-add procedure that common detector gather does frequency field, obtains frequency field overlapped information; Need to obtain the passiveseismic accurate delay time owing to doing superposition in time domain, and in the passiveseismic information usually extracted, because the signal of passiveseismic information is extremely faint, want to obtain passiveseismic accurately very difficult for time delay.Therefore utilize frequency field passiveseismic information at common detector gather frequency field stacking method in the embodiment of the present invention, without the need to obtaining the passiveseismic accurate delay time, and utilize passive wave field information effectively can also suppress the random noise of frequency field in the superposition of common detector gather frequency field, strengthen the signal of stable passiveseismic information.
Step 104, is shown as frequency field stacked section by described frequency field overlapped information according to the order of detection period, so that the distribution on sectional view of observed frequency territory overlapped information more intuitively;
Step 105, identifies the scope of energy exception from described frequency field stacked section, and then whether prediction is containing hydrocarbon-bearing pool.
By the method in the above embodiment of the present invention, achieve and utilize the passiveseismic information obtained in initial passiveseismic information to carry out predicting hydrocarbon reservoirs distribution range, if the passive wave field received comprises the reservoir information of underground, then can occur that energy is abnormal at the band segment of oily scope, and without the need to artificial excitation, both the time of later stage seism processing, reservoir prediction can have been saved, also the precision of prediction of oil-gas reserve can be improved, and the emphasis of later stage seism processing, reservoir prediction can be focused on, further look for oil and gas hides the evidence existed.
In the method for the passiveseismic wave field predicting hydrocarbon reservoirs that the embodiment of the present invention provides, preferably, after converting described time domain passiveseismic information to frequency field passiveseismic information, described method also comprises: the root mean square energy described frequency field passiveseismic information being calculated point frequency segment information in described frequency field, and draw the root mean square energy trace of point frequency segment information, analyze with this and verify the change of energy in frequency and space.
Energy abnormal ranges is identified from the root mean square energy trace of described point of frequency band, in conjunction with the method common identification hydrocarbon-bearing pool of the above-mentioned scope in frequency field stacked section identification energy exception, can improve the reliability of identification, reach and make the efficiency of prediction higher, accuracy is more accurate.
Wherein, in the method for the passiveseismic wave field predicting hydrocarbon reservoirs that the embodiment of the present invention provides, preferably, according to formula: S stk ( ω ) = S nor m 1 ( ω ) + S norm 2 ( ω ) + S norm 3 ( ω ) + . . . + S norm n ( ω ) n , The overlap-add procedure of frequency field is done at described common detector gather; Wherein, S stk(ω) be the information after the superposition of common detector gather frequency field, for the seismic trace of common detector gather, 1,2,3 ... n is the number of channels of common detector gather.Be illustrated in figure 5 frequency field passiveseismic information in the frequency field diagrammatic cross-section of common detector gather through single-coverage, by frequency field passiveseismic information at common detector gather through multi-fold, namely do the overlap-add procedure of frequency field, frequency field stacked section as shown in Figure 6 can be obtained; Wherein marking C position energy occurs abnormal, namely can regard as drilling well position, in above-described embodiment, common detector gather frequency field be done overlap-add procedure and not only suppressed noise later, also makes energy abnormal more obvious.
In the method for the passiveseismic wave field predicting hydrocarbon reservoirs that the embodiment of the present invention provides, preferably, according to formula:
A rms = 1 L ( Σ S stk ( ω i ) 2 ) 1 2
Calculate the root mean square energy of described point frequency segment information; Wherein, A rmswhat divide frequency segment information is root mean square energy, and L is number of samples in point frequency band, S stki) be the energy value of frequency field passiveseismic information at i sampling point place.As shown in Figure 7, after described frequency field passiveseismic information is divided into different frequency bands, the root mean square energy trace drawn by the root mean square energy value obtained, wherein marking D position is energy exception, therefore can regard as drilling well position.Known from Fig. 6 and Fig. 7, frequency field stacked section is consistent with the out-of-the way position of the root mean square energy trace of point frequency segment information, more demonstrates drilling well result thus.
By said method, do not need artificial excitation, and passiveseismic information is the position utilizing various abnormal information directly to determine hydrocarbon-bearing pool, therefore can save a large amount of manpowers, time, financial cost and environmental cost.
The embodiment of the present invention also provides a kind of device of passiveseismic wave field predicting hydrocarbon reservoirs, and as shown in Figure 8, it comprises:
Information acquisition unit 801, for filtering acquisition time territory passiveseismic information;
Converting unit 802, for converting described time domain passiveseismic information to frequency field passiveseismic information; Preferably, described converting unit 802 utilizes formula: S ( ω , τ ) = ∫ - ∞ ∞ s ( t ) g ( t - τ ) e - jωτ dt , Convert described time domain passiveseismic information to frequency field passiveseismic information, wherein, S (ω, τ) is frequency field passiveseismic information, and s (t) is time domain passiveseismic information, and g (t) is window function.
The above-mentioned frequency field that time domain converted to is because identical wave field passiveseismic information frequency characteristic of field is similar, random noise does not then have the feature of similar frequency field, therefore in follow-up process, just passiveseismic information signal can be strengthened by superposition, compacting noise.
Processing unit 803, the overlap-add procedure doing frequency field for described frequency field passiveseismic information being sorted into common detector gather obtains frequency field overlapped information; Need to obtain the passiveseismic accurate delay time owing to doing superposition in time domain, and in the passiveseismic information usually extracted, because the signal of passiveseismic information is extremely faint, want to obtain passiveseismic accurately very difficult for time delay.Therefore utilize frequency field passiveseismic information at common detector gather frequency field stacking method in the embodiment of the present invention, without the need to obtaining the passiveseismic accurate delay time, and utilize passive wave field information effectively can also suppress the random noise of frequency field in the superposition of common detector gather frequency field, strengthen the signal of stable passiveseismic information.
Section display unit 804, for described frequency field overlapped information is shown as frequency field stacked section according to the order of detection period, so that the distribution on sectional view of observed frequency territory overlapped information more intuitively;
Predicting unit 806, for identifying the scope of energy exception from described frequency field stacked section, and then whether prediction is containing hydrocarbon-bearing pool.
In the device of the passiveseismic wave field predicting hydrocarbon reservoirs provided in the embodiment of the present invention, preferably, described device also comprises:
Storage unit 805, for described frequency field passiveseismic information to be calculated the root mean square energy of point frequency segment information in described frequency field, and draws the root mean square energy trace of point frequency segment information; With the change of this analysing energy in frequency and space.
Predicting unit 806 also for identifying energy abnormal ranges from the root mean square energy trace of described point of frequency segment information, in conjunction with method common identification and the checking hydrocarbon-bearing pool of the above-mentioned scope in frequency field stacked section identification energy exception, the reliability of identification can be improved, reach and make the efficiency of prediction higher, accuracy is more accurate.
In the device of the passiveseismic wave field predicting hydrocarbon reservoirs provided in the embodiment of the present invention, preferably, described device also comprises: filter element 807, for filtering from initial passiveseismic information lower than very first time length value and the information exceeding first frequency value; Wherein: described very first time length is not higher than 500ms, and described first frequency value is not less than 40Hz.Preferably, wherein very first time length is 500ms, first frequency value is 40Hz, because time span there will be low-frequency information distortion lower than the short passiveseismic information of time domain of 500ms when being transformed into frequency field, equally, frequency is serious more than the usual interfering noise of the high-frequency information of 40Hz, and they all may produce illusion, therefore first filtering is needed, with the impact making the time domain passiveseismic information obtained eliminate distortion and noise.
Wherein, in one embodiment of the invention, initial passiveseismic information can obtain from the source book information of common seismic collection; In the source book information that common seismic gathers, filter out artificial earthquake information, only retain man-made explosion primary wave time passiveseismic information in the past as initial passiveseismic information; Subsequently, lower than the information of very first time length and the information exceeding first frequency value in the initial passiveseismic information of filtering, time domain passiveseismic information can be obtained.
In another embodiment of the invention, initial passive information directly can also receive passiveseismic source book from wave detector and read, because the passiveseismic information utilizing wave detector directly to receive only comprises initial passiveseismic information, do not comprise artificial earthquake information, can facilitate, directly obtain initial passiveseismic information.
In the device of the passiveseismic wave field predicting hydrocarbon reservoirs provided in the embodiment of the present invention, preferably, described processing unit 803 is according to formula:
S stk ( ω ) = S nor m 1 ( ω ) + S norm 2 ( ω ) + S norm 3 ( ω ) + . . . + S norm n ( ω ) n
The overlap-add procedure doing frequency field at described common detector gather obtains frequency field overlapped information, wherein, and S stk(ω) be the information after the superposition of common detector gather frequency field, for the seismic trace of common detector gather, 1,2,3 ... n is the number of channels of common detector gather, and the superposition of common detector gather frequency field has not only suppressed noise later, also makes energy abnormal more obvious.
In the device of the passiveseismic wave field predicting hydrocarbon reservoirs that the embodiment of the present invention provides, preferably, also comprise: computing unit 808, for according to formula:
A rms = 1 L ( Σ S stk ( ω i ) 2 ) 1 2
Calculate the root mean square energy of described point frequency segment information; Wherein, A rmsfor the root mean square energy of point frequency segment information, L is number of samples in point frequency band, S stki) be the energy value of frequency field passiveseismic information at i sampling point place, by the root mean square energy trace of point frequency segment information that obtains in conjunction with the common predicting hydrocarbon reservoirs of frequency field stacked section, the efficiency of prediction can be made higher, and accuracy is more accurate.
Whereby, the common predicting hydrocarbon reservoirs of root mean square energy trace of the frequency field stacked section that the embodiment of the present invention utilizes passiveseismic information to obtain and point frequency segment information, do not need artificial excitation, and a large amount of manpowers, time, financial cost and environmental cost can be saved.
Above-described specific embodiment; object of the present invention, technical scheme and beneficial effect are further described; be understood that; the foregoing is only specific embodiments of the invention; the protection domain be not intended to limit the present invention; within the spirit and principles in the present invention all, any amendment made, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.

Claims (8)

1. a method for passiveseismic wave field predicting hydrocarbon reservoirs, is characterized in that, described method comprises:
Acquisition time territory passiveseismic information;
Convert described time domain passiveseismic information to frequency field passiveseismic information;
Described frequency field passiveseismic information is sorted into the overlap-add procedure that common detector gather does frequency field, obtains frequency field overlapped information;
Described frequency field overlapped information is shown as frequency field stacked section according to the order of detection period;
From the scope of described frequency field stacked section identification energy exception, and then predicting hydrocarbon reservoirs;
Wherein, according to formula:
S stk ( ω ) = S norm 1 ( ω ) + S norm 2 ( ω ) + S norm 3 ( ω ) + . . . + S norm n ( ω ) n
The overlap-add procedure of frequency field is done at described common detector gather; Wherein, S stk(ω) be described frequency field overlapped information at described common detector gather, for the seismic trace of common detector gather, 1,2,3 ... n is the number of channels of common detector gather.
2. the method for passiveseismic wave field predicting hydrocarbon reservoirs according to claim 1, is characterized in that, described convert described time domain passiveseismic information to frequency field passiveseismic information after, described method also comprises:
Described frequency field passiveseismic information is calculated the root mean square energy of point frequency segment information in described frequency field, and draw the root mean square energy trace of point frequency segment information;
From the root mean square energy trace of described point of frequency segment information, identify the scope of energy exception, and then verify described hydrocarbon-bearing pool.
3. the method for passiveseismic wave field predicting hydrocarbon reservoirs according to claim 1, it is characterized in that, described acquisition time territory passiveseismic information comprises: from initial passiveseismic information, filtering is lower than very first time length value and the information exceeding first frequency value; Wherein: described very first time length is not higher than 500ms, and described first frequency value is not less than 40Hz.
4. the method for passiveseismic wave field predicting hydrocarbon reservoirs according to claim 2, is characterized in that, according to formula:
A rms = 1 L ( Σ S stk ( ω i ) 2 ) 1 2
Calculate the root mean square energy of described point frequency segment information; Wherein, A rmsfor the root mean square energy value of point frequency segment information, L is number of samples in point frequency band, S stki) be the energy value of frequency field passiveseismic information at i sampling point place.
5. a device for passiveseismic wave field predicting hydrocarbon reservoirs, is characterized in that, comprising:
Information acquisition unit, for acquisition time territory passiveseismic information;
Converting unit, for converting described time domain passiveseismic information to frequency field passiveseismic information;
Processing unit, for described frequency field passiveseismic information is sorted into the overlap-add procedure that common detector gather does frequency field, obtains frequency field overlapped information;
Section display unit, for being shown as frequency field stacked section by described frequency field overlapped information according to the order of detection period;
Predicting unit, for identifying the scope of energy exception from described frequency field stacked section, and then predicting hydrocarbon reservoirs;
Wherein, described processing unit is according to formula:
S stk ( ω ) = S norm 1 ( ω ) + S norm 2 ( ω ) + S norm 3 ( ω ) + . . . + S norm n ( ω ) n
The overlap-add procedure doing frequency field at described common detector gather obtains frequency field overlapped information, wherein, and S stk(ω) be the information after the superposition of common detector gather frequency field, for the seismic trace of common detector gather, 1,2,3 ... n is the number of channels of common detector gather.
6. the device of passiveseismic wave field predicting hydrocarbon reservoirs according to claim 5, it is characterized in that, described device also comprises:
Storage unit, for described frequency field passiveseismic information being calculated in described frequency field the root mean square energy of point frequency segment information, and draws the root mean square energy trace of point frequency segment information;
Described predicting unit also for identifying energy abnormal ranges from the root mean square energy trace of described point of frequency segment information, and then verifies described hydrocarbon-bearing pool.
7. the device of passiveseismic wave field predicting hydrocarbon reservoirs according to claim 5, it is characterized in that, described device also comprises: filter element, for filtering from initial passiveseismic information lower than very first time length value and the information exceeding first frequency value; Wherein: described very first time length is not higher than 500ms, and described first frequency value is not less than 40Hz.
8. the device of passiveseismic wave field predicting hydrocarbon reservoirs according to claim 6, it is characterized in that, described device also comprises: computing unit, for according to formula:
A rms = 1 L ( Σ S stk ( ω i ) 2 ) 1 2
Calculate the root mean square energy of described point frequency segment information; Wherein, A rmsfor the root mean square energy of point frequency segment information, L is number of samples in point frequency band, S stki) be the energy value of frequency field passiveseismic information at i sampling point place.
CN201210443023.0A 2012-11-08 2012-11-08 Method and device for predicting oil and gas reservoirs through passive seismic wave field Active CN102967878B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210443023.0A CN102967878B (en) 2012-11-08 2012-11-08 Method and device for predicting oil and gas reservoirs through passive seismic wave field

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210443023.0A CN102967878B (en) 2012-11-08 2012-11-08 Method and device for predicting oil and gas reservoirs through passive seismic wave field

Publications (2)

Publication Number Publication Date
CN102967878A CN102967878A (en) 2013-03-13
CN102967878B true CN102967878B (en) 2015-07-08

Family

ID=47798146

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210443023.0A Active CN102967878B (en) 2012-11-08 2012-11-08 Method and device for predicting oil and gas reservoirs through passive seismic wave field

Country Status (1)

Country Link
CN (1) CN102967878B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104765063B (en) * 2015-03-13 2018-08-14 中国石油天然气股份有限公司 The gas-oil detecting method and device of attenuation by absorption attribute are calculated based on frequency spectrum

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101923177B (en) * 2010-07-09 2012-05-23 中国石油天然气股份有限公司 Velocity analysis method and device of multi-wave seismic data

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8441891B2 (en) * 2009-12-30 2013-05-14 Westerngeco L.L.C. Seismic data acquisition and processing quality control

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101923177B (en) * 2010-07-09 2012-05-23 中国石油天然气股份有限公司 Velocity analysis method and device of multi-wave seismic data

Also Published As

Publication number Publication date
CN102967878A (en) 2013-03-13

Similar Documents

Publication Publication Date Title
Lecocq et al. Monitoring ground water storage at mesoscale using seismic noise: 30 years of continuous observation and thermo-elastic and hydrological modeling
Hsiao et al. Development of earthquake early warning system in Taiwan
CN104142519B (en) Mud rock crack oil deposit predicting method
CN102373923B (en) Reservoir stratum identification method
CN102809762B (en) Reservoir imaging technique based on full-frequency-band seismic information mining
CN103293562B (en) Method and device for determining carbonate reservoir geological reserves
CN103235339A (en) Time-frequency decomposition earthquake-fluid recognition method
Martin-Short et al. Mantle flow geometry from ridge to trench beneath the Gorda–Juan de Fuca plate system
Koçkar et al. Evaluation of the site effects of the Ankara basin, Turkey
CN104329079A (en) Method and system for recognizing gas logging oil and gas reservoir
Maxwell et al. Microseismic geomechanics of hydraulic-fracture networks: Insights into mechanisms of microseismic sources
CN104737035A (en) System and method for detecting a fracture in a rock formation using an electromagnetic source
CN102914799B (en) Forward modeling method and device for nonequivalent wave field
CN107255831A (en) A kind of extracting method of prestack frequency dispersion attribute
Hanssen et al. Pitfalls in the analysis of low frequency passive seismic data
Ouenes et al. Integrated characterization and simulation of the fractured Tensleep Reservoir at Teapot Dome for CO2 injection design
CN102305940B (en) Method for extracting fluid factor
CN103758511A (en) Method and device for recognizing hidden reservoirs by underground reverse-time migration imaging
CN104007462A (en) Fracture prediction method based on attenuation anisotropy
CN107356668A (en) A kind of detection method of safety of dam body identification
CN103792575A (en) Densified sandstone air layer identification evaluation method free from influence of frame
Shan et al. Spectral decomposition and a waveform cluster to characterize strongly heterogeneous paleokarst reservoirs in the Tarim basin, China
CN102914797A (en) Method and device for acquiring anisotropy coefficient of stratum
CN102967878B (en) Method and device for predicting oil and gas reservoirs through passive seismic wave field
CN112711068B (en) Method and device for predicting effective reservoir of oil gas in sandstone

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