CN106291686A - The method and device of spectral balancing between single-frequency data volume - Google Patents
The method and device of spectral balancing between single-frequency data volume Download PDFInfo
- Publication number
- CN106291686A CN106291686A CN201610576825.7A CN201610576825A CN106291686A CN 106291686 A CN106291686 A CN 106291686A CN 201610576825 A CN201610576825 A CN 201610576825A CN 106291686 A CN106291686 A CN 106291686A
- Authority
- CN
- China
- Prior art keywords
- frequency
- data volume
- spectral
- frequency data
- spectral balancing
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
Abstract
The embodiment of the present application provides the method and device of spectral balancing between a kind of single-frequency data volume, and the method comprises determining that the analysis window of target zone;Time domain seismic data cube in the range of analysis window is done time frequency analysis, obtains crest frequency and effective band;The discrete frequency of predetermined quantity is selected in effective band scope;Time domain seismic data cube is done spectral decomposition, obtains the single-frequency data volume that each discrete frequency is corresponding;By the absolute value sum of the amplitude of all sampled points in the range of analysis window of all seismic channels in single-frequency data volume corresponding for crest frequency, all seismic channels ratio of the absolute value sum of the amplitude of all sampled points in the range of analysis window in the single-frequency data volume corresponding with each discrete frequency, as the weighting function of each discrete frequency correspondence single-frequency data volume;Each single-frequency data volume is multiplied with its weighting function, obtains the single-frequency data volume after the spectral balancing of correspondence.The embodiment of the present application has suppressed abnormal strong energy, it is ensured that the energy relativeness between seismic channel.
Description
Technical field
The application relates to seismic data processing technology field, especially relates to a kind of method of spectral balancing between single-frequency data volume
And device.
Background technology
The Energy distribution of seismic wavelet is unbalanced, and at dominant frequency, energy is the strongest, in dominant frequency both sides along with frequency change by
The least.And the amplitude information of earthquake record is the convolution of seismic wavelet and reflection coefficient, thus amplitude spectrum can be folded by wavelet
The impact of print (wavelet overprint), i.e. energy is unbalanced in each frequency distribution, is concentrated mainly near main band,
Therefore generally after time frequency analysis, (its image is commonly referred to as frequency division section, the most different frequencies to different single-frequency data volumes
The corresponding different single-frequency section of rate) carry out spectral balancing by weighting function.Spectral balancing can be follow-up inverting and layer description carries
For convenient.
The most conventional spectral balancing method is: choose the dominant frequency of source wavelet as reference frequency f0, reference frequency pair
Answer the ratio of maximum of the maximum single-frequency data volume corresponding with other frequency of single-frequency data volume as weighting function, this weighting
Function acquiring method can represent by equation below:
In formula, t is sampled point sequence number, and n is seismic channel, fmFor the frequency values of m-th assigned frequency, f0For reference frequency.
max(A(n,t,fm)) represent assigned frequency fmPeak swing value in the corresponding all sampled points in single-frequency data volume the n-th road, max (A
(n,t,f0)) represent reference frequency f0Peak swing value in the corresponding all sampled points in single-frequency data volume the n-th road, wt(n,fm) generation
Table assigned frequency fmThe weighting function that n-th road of corresponding single-frequency data volume is corresponding.Utilize the weighting function asked for, obtain spectrum all
Single-frequency data volume after weighing apparatus:
M (n, t, f)=S (n, t, f) w (n, f)
As can be seen here, above weighting function is both for often carrying out spectral balancing together, due to the power letter of every application together
Number difference, therefore causes the energy of each single-frequency section after equilibrium to change along with the change of frequency, and this is by for follow-up reservoir
Inverting brings difficulty.Additionally, the method is to ask for the maximum of every one single-frequency data, utilize per pass single-frequency data peak swing
Zhi Zhu road carries out spectral balancing, does not accounts for other sampling points in addition to the maximum correspondence sampling point of single-frequency data of per pass.Once
Run into abnormal strong energy in seismic data, again with abnormal strong energy after this road will be made to be equalized, i.e. cannot suppress exception
Strong energy.
Summary of the invention
The purpose of the embodiment of the present application is to provide the method and device of spectral balancing between a kind of single-frequency data volume, to realize
While the abnormal strong energy of compacting, it is to avoid the energy of each single-frequency section after equilibrium changes along with the change of frequency, is beneficial to
Improve the inversion result of follow-up reservoir inversion.
For reaching above-mentioned purpose, on the one hand, the embodiment of the present application provides a kind of method of spectral balancing between single-frequency data volume,
Comprise the following steps:
Determine the analysis window of target zone;
It is pointed to the time domain seismic data cube in the range of described analysis window and carries out time frequency analysis, obtain crest frequency and have
Effect frequency band;
Predetermined quantity is selected to be used for the discrete frequency of spectral balancing in the range of described effective band, described for spectral balancing
Discrete frequency include described crest frequency;
Described time domain seismic data cube is carried out spectral decomposition, obtains the described predetermined quantity the discrete frequency for spectral balancing
The single-frequency data volume that in point, each discrete frequency is corresponding;
Adopt in single-frequency data volume corresponding for described crest frequency, all seismic channels are all in the range of described analysis window
The absolute value sum of the amplitude of sampling point, in the single-frequency data volume corresponding with described each discrete frequency, all seismic channels are at described point
The ratio of the absolute value sum of the amplitude of all sampled points in the range of window during analysis, as each discrete frequency correspondence single-frequency data volume
Weighting function;
Each single-frequency data volume is multiplied with its weighting function, obtains the single-frequency data volume after the spectral balancing of correspondence.
The method of spectral balancing between the single-frequency data volume of the embodiment of the present application, the scope foundation of the analysis window of described target zone
The scope of described target zone determines.
The method of spectral balancing between the single-frequency data volume of the embodiment of the present application, described time frequency analysis uses Fourier transformation real
Existing.
The method of spectral balancing between the single-frequency data volume of the embodiment of the present application, described spectral decomposition is become by Fourier in short-term
Change, S-transformation or add uncle conversion realize.
The method of spectral balancing between the single-frequency data volume of the embodiment of the present application, the quantity of the described discrete frequency for spectral balancing
For odd number, and wherein the various discrete frequency in addition to described crest frequency is uniformly distributed in described crest frequency both sides.
On the other hand, the embodiment of the present application additionally provides the device of spectral balancing between a kind of single-frequency data volume, time window determine mould
Block, for determining the analysis window of target zone;
Time frequency analysis module, frequency division when the time domain seismic data cube being pointed in the range of described analysis window is carried out
Analysis, obtains crest frequency and effective band;
Frequently point selection module, for selecting predetermined quantity for the discrete frequency of spectral balancing in the range of described effective band
Point, the described discrete frequency for spectral balancing includes described crest frequency;
Single-frequency data volume acquisition module, for described time domain seismic data cube is carried out spectral decomposition, obtains described predetermined
Quantity the single-frequency data volume that each discrete frequency is corresponding in the discrete frequency of spectral balancing;
Weighting function determines module, is used for all seismic channels in single-frequency data volume corresponding for described crest frequency described
The absolute value sum of the amplitude of all sampled points in the range of analysis window, the single-frequency data volume corresponding with described each discrete frequency
In all seismic channels ratio of the absolute value sum of the amplitude of all sampled points in the range of described analysis window, as respectively from
Dissipate the weighting function of frequency correspondence single-frequency data volume;
Spectral balancing module, for being multiplied each single-frequency data volume with its weighting function, after obtaining the spectral balancing of correspondence
Single-frequency data volume.
The device of spectral balancing between the single-frequency data volume of the embodiment of the present application, the scope foundation of the analysis window of described target zone
The scope of described target zone determines.
The device of spectral balancing between the single-frequency data volume of the embodiment of the present application, described time frequency analysis uses Fourier transformation real
Existing.
The device of spectral balancing between the single-frequency data volume of the embodiment of the present application, described spectral decomposition is become by Fourier in short-term
Change, S-transformation or add uncle conversion realize.
The device of spectral balancing between the single-frequency data volume of the embodiment of the present application, the quantity of the described discrete frequency for spectral balancing
For odd number, and wherein the various discrete frequency in addition to described crest frequency is uniformly distributed in described crest frequency both sides.
First the embodiment of the present application determines the analysis window of target zone;Next is pointed in the range of analysis window in time domain
Shake data volume carries out time frequency analysis, obtains crest frequency and effective band;Then in the range of effective band, select predetermined quantity
The individual discrete frequency for spectral balancing;Then time domain seismic data cube is carried out spectral decomposition, obtain predetermined quantity for composing
The single-frequency data volume that in the discrete frequency of equilibrium, each discrete frequency is corresponding;Secondly by single-frequency data volume corresponding for crest frequency
All seismic channels are the absolute value sum of the amplitude of all sampled points in the range of analysis window, the list corresponding with each discrete frequency
Frequency according to internal all seismic channels ratio of the absolute value sum of the amplitude of all sampled points in the range of analysis window, as
The weighting function of each discrete frequency correspondence single-frequency data volume;Finally each single-frequency data volume is multiplied with its weighting function, obtains
The corresponding single-frequency data volume after spectral balancing.Due to the embodiment of the present application construct weighting function time, using single-frequency data volume as
Entirety, all lists its all seismic channels amplitude of all sampled points in the range of analysis window in computer capacity, it is ensured that ground
Energy relativeness between shake road, has broken conventional spectral balancing and has built the thinking of weighting function by road, it is to avoid Liao Zhu road has built
Use maximum acquiring method to be easily mixed into incorrect energy during weighting function thus affect weight function and ask for the problem of accuracy, from
And achieving while the abnormal strong energy of compacting, it is to avoid the energy of each single-frequency section after equilibrium becomes along with the change of frequency
Change.Thus be conducive to improving the inversion result of follow-up reservoir inversion.
Accompanying drawing explanation
Accompanying drawing described herein is used for providing being further appreciated by the embodiment of the present application, constitutes the embodiment of the present application
A part, is not intended that the restriction to the embodiment of the present application.In the accompanying drawings:
Fig. 1 be the embodiment of the present application single-frequency data volume between the flow chart of method of spectral balancing;
Fig. 2 a-2d is respectively in the application one embodiment, and spectral balancing processes the list that front 10Hz, 20Hz, 30Hz, 40Hz are corresponding
Frequently section;
Fig. 3 a-3d is respectively in the application one embodiment, obtains after the spectral balancing scheme of the embodiment of the present application processes
The single-frequency section that 10Hz, 20Hz, 30Hz, 40Hz are corresponding;
Fig. 4 a-4d is respectively in the application one embodiment, the 10Hz that obtains after conventional spectral balancing scheme processes, 20Hz,
The single-frequency section that 30Hz, 40Hz are corresponding;
Fig. 5 a is the inversion result of the single-frequency section that spectral balancing scheme based on the embodiment of the present application obtains;
Fig. 5 b is the inversion result of the single-frequency section obtained based on conventional spectral balancing scheme;
Fig. 6 be the embodiment of the present application single-frequency data volume between the structured flowchart of device of spectral balancing.
Detailed description of the invention
For making the purpose of the embodiment of the present application, technical scheme and advantage clearer, below in conjunction with embodiment and attached
Figure, is described in further details the embodiment of the present application.Here, the schematic description and description of the embodiment of the present application is used for
Explain the embodiment of the present application, but be not intended as the restriction to the embodiment of the present application.
Below in conjunction with the accompanying drawings, the detailed description of the invention of the embodiment of the present application is described in further detail.
With reference to shown in Fig. 1, between the single-frequency data volume of the embodiment of the present application, the method for spectral balancing comprises the following steps:
Step S101, determine the analysis window of target zone.
General, the analysis window of target zone is chosen and is wanted suitably, and analysis window is chosen too small, and destination layer can be made to exceed process
Scope, and too small time window also can affect the effectiveness of some parameter;And time window choose excessive, the feature of target zone can be made
Extremely it is submerged in Anomalies of Backgrounds.Therefore, the scope of the analysis window of target zone needs to determine according to the scope of target zone.The most right
One there is N road, per pass has the time domain seismic data cube a of T sampled point (n, t), wherein n is seismic channel sequence number, and scope is 1
~N, t are sampled point sequence number, scope is 1~T, determines analysis window scope t according to the scope of target zone1~t2, 1≤t1< t2
≤T。
In a specific embodiment of the application, the most known seismic channel serial number 1~70, sampled point serial number 1~
651, then analysis window scope can be chosen as between the 450th~550 sampled point.
Step S102, it is pointed to the time domain seismic data cube in the range of described analysis window and carries out time frequency analysis, obtain peak
Value frequency and effective band.
In a specific embodiment of the application, such as to time domain seismic data cube a, (n, t), at analysis window model
Carry out spectrum analysis in enclosing (such as the 450th~550 sampled point), obtain the crest frequency f of seismic amplitude spectrum0And the most frequently
Band scope fmin~fmax。
The most described time frequency analysis can use Fourier transformation to realize.
Step S103, select in the range of described effective band predetermined quantity for the discrete frequency of spectral balancing, described
Discrete frequency for spectral balancing includes described crest frequency.
In the embodiment of the present application, the described quantity for the discrete frequency of spectral balancing is odd number, and wherein except described
Various discrete frequency outside crest frequency is uniformly distributed in described crest frequency both sides.A specific embodiment in the application
In, such as in effective band scope fmin~fmaxIn can select 5 discrete frequencies f for spectral balancing0、f1、f2、f3、f4,
fmin≤f1< f2< f0< f3< f4≤fmax.Assume f0For 28Hz, then can select f1=10Hz, f2=20Hz, f3=30Hz,
f4=40Hz.
Step S104, described time domain seismic data cube is carried out spectral decomposition, obtain described predetermined quantity for composing all
The single-frequency data volume that in the discrete frequency of weighing apparatus, each discrete frequency is corresponding.
In the embodiment of the present application, described spectral decomposition by Short Time Fourier Transform, S-transformation or adds uncle conversion (Gabor
Conversion) etc. realize.To time domain seismic data cube a, (n, after t) carrying out spectral decomposition, can obtain each discrete frequency fmCorresponding
Single-frequency data volume Am(n,t,fm), wherein m=0,1,2,3,4.
In a specific embodiment of the application, select to add uncle's conversion and described time domain seismic data cube is carried out frequency spectrum divide
Solve, available A1(n,t,f1), A2(n,t,f2), A3(n,t,f3), A4(n,t,f4), i.e. composing as shown in Fig. 2 a~Fig. 2 d
The single-frequency section of equilibrium.Wherein, Fig. 2 a, Fig. 2 b, single-frequency corresponding for Fig. 2 c with 2d respectively 10Hz, 20Hz, 30Hz, 40Hz are cutd open
Face.From Fig. 2 a~Fig. 2 d, can be seen that the single-frequency section energy that Fig. 2 b, i.e. 20Hz are corresponding is the strongest, Fig. 2 a, Fig. 2 c and Fig. 2 d couple
The single-frequency section energy answered is relatively weak, and, by comparison diagram 2b, Fig. 2 c, Fig. 2 d it can be seen that frequency is more remote from peak value
Frequency f0, the energy of the single-frequency section of its correspondence is the most weak.
Step S105, by all seismic channels in single-frequency data volume corresponding for described crest frequency in described analysis window scope
The absolute value sum of the amplitude of interior all sampled points, all seismic channels in the single-frequency data volume corresponding with described each discrete frequency
The ratio of the absolute value sum of the amplitude of all sampled points in the range of described analysis window, as each discrete frequency correspondence list
Frequently the weighting function of data volume.
In the embodiment of the present application, to each single-frequency data volume A (n, t, fm) all roads exist in (m=0,1,2,3,4)
Described analysis window scope (t1~t2The ratio of the absolute value sum of the amplitude of all sampled points in), available each single-frequency number
Amplitude and S according to bodym, it is expressed as with formulaBy described crest frequency correspondence single-frequency data volume
Amplitude and S0, with discrete frequency fmThe amplitude of corresponding single-frequency data volume and Sm(m=0,1,2,3,4) ratio is discrete to should be used as
Frequency fmCorresponding single-frequency data volume A (n, t, fm) weighting function w (fm), i.e.
Step S106, each single-frequency data volume is multiplied with its weighting function, obtains the single-frequency number after the spectral balancing of correspondence
According to body.
In the embodiment of the present application, by each single-frequency data volume A (n, t, fm) (m=0,1,2,3,4) corresponding weighting
Function w (fm) be multiplied, single-frequency data volume B (n, t, the f after available spectral balancingm), i.e. B (n, t, fm)=A (n, t, fm)w(fm),
Thus realize spectral balancing.
Fig. 3 a~Fig. 3 d is that each discrete frequency of 10Hz, 20Hz, 30Hz, 40Hz is at the spectral balancing scheme of the embodiment of the present application
The corresponding single-frequency section obtained after reason.Contrast 3a~Fig. 3 d it is found that through the embodiment of the present application spectral balancing scheme process after,
The energy of the single-frequency section that each discrete frequency is corresponding no longer changes along with the change of frequency, thus is conducive to improving follow-up reservoir
The inversion result of inverting.
For the ease of contrast, as shown in Fig. 4 a~Fig. 4 d, the embodiment of the present application also show 10Hz, 20Hz, 30Hz, 40Hz
Each discrete frequency is the corresponding single-frequency section obtained after conventional spectral balancing scheme processes.Comparison diagram 4a~Fig. 4 d is it is found that warp
After conventional spectral balancing scheme processes, although each single-frequency section energy no longer changes along with the change of frequency, but at Fig. 4 b, figure
Incorrect energy is occurred in that near 28th road of 4c and Fig. 4 d, when this is owing to using tradition spectral balancing algorithm, former geological data the
Being originally used near 28 roads calculates the continuously strong energy along layer of weighting function, is mixed into the abnormal strong of other layer of position of superficial part
Energy, and the wavelet of two positions has bigger difference, causes the weighting function near 28 roads exception also occur, thus at this
Single-frequency section there will be exception.If the single-frequency data volume that conventional for these application spectral balancing scheme obtains is used for reservoir inversion
After, the inversion result obtained is often with such incorrect energy.
Fig. 5 a and Fig. 5 b is the single-frequency that spectral balancing scheme based on the embodiment of the present application, conventional spectral balancing scheme obtain respectively
The inversion result of section.Comparison diagram 5a and Fig. 5 b is it is found that create incorrect energy near the 28th road in Fig. 5 b, and its reason is
The single-frequency section obtained through conventional spectral balancing scheme occurs in that incorrect energy in this position.
Therefore, energy between the single-frequency data volume Dao Yu road after the spectral balancing scheme of the embodiment of the present application ensure that spectral balancing
Relativeness, it is to avoid the generation of incorrect energy, to improving follow-up inversion result and to the effect of layer description the most meaningful.
Although procedures described above flow process includes the multiple operations occurred with particular order, it should however be appreciated that understand,
These processes can include more or less of operation.
With reference to shown in Fig. 6, between the single-frequency data volume of the embodiment of the present application, the device of spectral balancing includes:
Time window determine module 61, for determining the analysis window of target zone.
General, the analysis window of target zone is chosen and is wanted suitably, and analysis window is chosen too small, and destination layer can be made to exceed process
Scope, and too small time window also can affect the effectiveness of some parameter;And time window choose excessive, the feature of target zone can be made
Extremely it is submerged in Anomalies of Backgrounds.Therefore, the scope of the analysis window of target zone needs to determine according to the scope of target zone.The most right
One there is N road, per pass has the time domain seismic data cube a of T sampled point (n, t), wherein n is seismic channel sequence number, and scope is 1
~N, t are sampled point sequence number, scope is 1~T, determines analysis window scope t according to the scope of target zone1~t2, 1≤t1< t2
≤T。
In a specific embodiment of the application, the most known seismic channel serial number 1~70, sampled point serial number 1~
651, then analysis window scope can be chosen as between the 450th~550 sampled point.
Time frequency analysis module 62, frequency division when the time domain seismic data cube being pointed in the range of described analysis window is carried out
Analysis, obtains crest frequency and effective band.
In a specific embodiment of the application, such as to time domain seismic data cube a, (n, t), at analysis window model
Carry out spectrum analysis in enclosing (such as the 450th~550 sampled point), obtain the crest frequency f of seismic amplitude spectrum0And the most frequently
Band scope fmin~fmax。
The most described time frequency analysis can use Fourier transformation to realize.
Frequently point selection module 63, for selecting predetermined quantity to be used for the discrete of spectral balancing in the range of described effective band
Frequency, the described discrete frequency for spectral balancing includes described crest frequency.
In the embodiment of the present application, the described quantity for the discrete frequency of spectral balancing is odd number, and wherein except described
Various discrete frequency outside crest frequency is uniformly distributed in described crest frequency both sides.A specific embodiment in the application
In, such as in effective band scope fmin~fmaxIn can select 5 discrete frequencies f for spectral balancing0、f1、f2、f3、f4,
fmin≤f1< f2< f0< f3< f4≤fmax.Assume f0For 28Hz, then can select f1=10Hz, f2=20Hz, f3=30Hz,
f4=40Hz.
Single-frequency data volume acquisition module 64, for described time domain seismic data cube is carried out spectral decomposition, obtains described pre-
Determined number the single-frequency data volume that each discrete frequency is corresponding in the discrete frequency of spectral balancing.
In the embodiment of the present application, described spectral decomposition by Short Time Fourier Transform, S-transformation or adds uncle conversion (Gabor
Conversion) etc. realize.To time domain seismic data cube a, (n, after t) carrying out spectral decomposition, can obtain each discrete frequency fmCorresponding
Single-frequency data volume Am(n,t,fm), wherein m=0,1,2,3,4.
In a specific embodiment of the application, select to add uncle's conversion and described time domain seismic data cube is carried out frequency spectrum divide
Solve, available A1(n,t,f1), A2(n,t,f2), A3(n,t,f3), A4(n,t,f4), i.e. composing as shown in Fig. 2 a~Fig. 2 d
The single-frequency section of equilibrium.Wherein, Fig. 2 a, Fig. 2 b, single-frequency corresponding for Fig. 2 c with 2d respectively 10Hz, 20Hz, 30Hz, 40Hz are cutd open
Face.From Fig. 2 a~Fig. 2 d, can be seen that the single-frequency section energy that Fig. 2 b, i.e. 20Hz are corresponding is the strongest, Fig. 2 a, Fig. 2 c and Fig. 2 d couple
The single-frequency section energy answered is relatively weak, and, by comparison diagram 2b, Fig. 2 c, Fig. 2 d it can be seen that frequency is more remote from peak value
Frequency f0, the energy of the single-frequency section of its correspondence is the most weak.
Weighting function determines module 65, is used for all seismic channels in single-frequency data volume corresponding for described crest frequency in institute
The absolute value sum of the amplitude of all sampled points in the range of stating analysis window, the single-frequency data corresponding with described each discrete frequency
Internal all seismic channels are the ratio of the absolute value sum of the amplitude of all sampled points in the range of described analysis window, as from
Dissipate frequency fmThe weighting function of corresponding single-frequency data volume.
In the embodiment of the present application, to each single-frequency data volume A (n, t, fm) all roads exist in (m=0,1,2,3,4)
Described analysis window scope (t1~t2The ratio of the absolute value sum of the amplitude of all sampled points in), available each single-frequency number
Amplitude and S according to bodym, it is expressed as with formulaBy described crest frequency correspondence single-frequency data volume
Amplitude and S0, with discrete frequency fmThe amplitude of corresponding single-frequency data volume and Sm(m=0,1,2,3,4) ratio is discrete to should be used as
Frequency fmCorresponding single-frequency data volume A (n, t, fm) weighting function w (fm), i.e.Spectral balancing module 66, being used for will
Each single-frequency data volume is multiplied with its weighting function, obtains the single-frequency data volume after the spectral balancing of correspondence.
In the embodiment of the present application, by each single-frequency data volume A (n, t, fm) (m=0,1,2,3,4) corresponding weighting
Function w (fm) be multiplied, single-frequency data volume B (n, t, the f after available spectral balancingm), i.e. B (n, t, fm)=A (n, t, fm)w(fm),
Thus realize spectral balancing.
Fig. 3 a~Fig. 3 d is that each discrete frequency of 10Hz, 20Hz, 30Hz, 40Hz is at the spectral balancing scheme of the embodiment of the present application
The corresponding single-frequency section obtained after reason.Contrast 3a~Fig. 3 d it is found that through the embodiment of the present application spectral balancing scheme process after,
The energy of the single-frequency section that each discrete frequency is corresponding no longer changes along with the change of frequency, thus is conducive to improving follow-up reservoir
The inversion result of inverting.
For the ease of contrast, as shown in Fig. 4 a~Fig. 4 d, the embodiment of the present application also show 10Hz, 20Hz, 30Hz, 40Hz
Each discrete frequency is the corresponding single-frequency section obtained after conventional spectral balancing scheme processes.Comparison diagram 4a~Fig. 4 d is it is found that warp
After conventional spectral balancing scheme processes, although each single-frequency section energy no longer changes along with the change of frequency, but at Fig. 4 b, figure
Incorrect energy is occurred in that near 28th road of 4c and Fig. 4 d, when this is owing to using tradition spectral balancing algorithm, former geological data the
Being originally used near 28 roads calculates the continuously strong energy along layer of weighting function, is mixed into the abnormal strong of other layer of position of superficial part
Energy, and the wavelet of two positions has bigger difference, causes the weighting function near 28 roads exception also occur, thus at this
Single-frequency section there will be exception.If the single-frequency data volume that conventional for these application spectral balancing scheme obtains is used for reservoir inversion
After, the inversion result obtained is often with such incorrect energy.
Fig. 5 a and Fig. 5 b is the single-frequency that spectral balancing scheme based on the embodiment of the present application, conventional spectral balancing scheme obtain respectively
The inversion result of section.Comparison diagram 5a and Fig. 5 b is it is found that create incorrect energy near the 28th road in Fig. 5 b, and its reason is
The single-frequency section obtained through conventional spectral balancing scheme occurs in that incorrect energy in this position.
Therefore, energy between the single-frequency data volume Dao Yu road after the spectral balancing scheme of the embodiment of the present application ensure that spectral balancing
Relativeness, it is to avoid the generation of incorrect energy, to improving follow-up inversion result and to the effect of layer description the most meaningful.
Although procedures described above flow process includes the multiple operations occurred with particular order, it should however be appreciated that understand,
These processes can include more or less of operation.
Particular embodiments described above, has been carried out the most in detail purpose, technical scheme and the beneficial effect of the application
Describe in detail bright, be it should be understood that the specific embodiment that the foregoing is only the embodiment of the present application, be not used to limit this Shen
Protection domain please, all within spirit herein and principle, any modification, equivalent substitution and improvement etc. done, all should wrap
Within being contained in the protection domain of the application.
Claims (10)
1. the method for spectral balancing between a single-frequency data volume, it is characterised in that comprise the following steps:
Determine the analysis window of target zone;
It is pointed to the time domain seismic data cube in the range of described analysis window and carries out time frequency analysis, obtain crest frequency and effective frequency
Band;
Select in the range of described effective band predetermined quantity for the discrete frequency of spectral balancing, described for spectral balancing from
Scattered frequency includes described crest frequency;
Described time domain seismic data cube is carried out spectral decomposition, obtains described predetermined quantity in the discrete frequency of spectral balancing
The single-frequency data volume that each discrete frequency is corresponding;
By all seismic channels all sampled points in the range of described analysis window in single-frequency data volume corresponding for described crest frequency
The absolute value sum of amplitude, in the single-frequency data volume corresponding with described each discrete frequency, all seismic channels are when described analysis
The ratio of the absolute value sum of the amplitude of all sampled points in the range of window, as adding of each discrete frequency correspondence single-frequency data volume
Weight function;
Each single-frequency data volume is multiplied with its weighting function, obtains the single-frequency data volume after the spectral balancing of correspondence.
The method of spectral balancing between single-frequency data volume the most according to claim 1, it is characterised in that the analysis of described target zone
Time window scope determine according to the scope of described target zone.
The method of spectral balancing between single-frequency data volume the most according to claim 1, it is characterised in that described time frequency analysis uses
Fourier transformation realizes.
The method of spectral balancing between single-frequency data volume the most according to claim 1, it is characterised in that described spectral decomposition is passed through
Short Time Fourier Transform, S-transformation or add uncle conversion realize.
The method of spectral balancing between single-frequency data volume the most according to claim 1, it is characterised in that described for spectral balancing
The quantity of discrete frequency is odd number, and wherein the various discrete frequency in addition to described crest frequency is uniformly distributed in described peak value
Frequency both sides.
6. the device of spectral balancing between a single-frequency data volume, it is characterised in that including:
Time window determine module, for determining the analysis window of target zone;
Time frequency analysis module, the time domain seismic data cube in the range of being pointed to described analysis window carries out time frequency analysis, obtains
Take crest frequency and effective band;
Frequently point selection module, for selecting predetermined quantity to be used for the discrete frequency of spectral balancing in the range of described effective band,
The described discrete frequency for spectral balancing includes described crest frequency;
Single-frequency data volume acquisition module, for described time domain seismic data cube is carried out spectral decomposition, obtains described predetermined quantity
The individual single-frequency data volume that each discrete frequency is corresponding in the discrete frequency of spectral balancing;
Weighting function determines module, is used for all seismic channels in single-frequency data volume corresponding for described crest frequency in described analysis
Time window in the range of the absolute value sum of amplitude of all sampled points, institute in the single-frequency data volume corresponding with described each discrete frequency
There is seismic channel ratio of the absolute value sum of the amplitude of all sampled points in the range of described analysis window, as each discrete frequency
The weighting function of the corresponding single-frequency data volume of point;
Spectral balancing module, for being multiplied with its weighting function by each single-frequency data volume, obtains the single-frequency after the spectral balancing of correspondence
Data volume.
The device of spectral balancing between single-frequency data volume the most according to claim 6, it is characterised in that the analysis of described target zone
Time window scope determine according to the scope of described target zone.
The device of spectral balancing between single-frequency data volume the most according to claim 6, it is characterised in that described time frequency analysis uses
Fourier transformation realizes.
The device of spectral balancing between single-frequency data volume the most according to claim 6, it is characterised in that described spectral decomposition is passed through
Short Time Fourier Transform, S-transformation or add uncle conversion realize.
The device of spectral balancing between single-frequency data volume the most according to claim 6, it is characterised in that described for spectral balancing
The quantity of discrete frequency be odd number, and wherein the various discrete frequency in addition to described crest frequency is uniformly distributed in described peak
Value frequency both sides.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610576825.7A CN106291686B (en) | 2016-07-20 | 2016-07-20 | The method and device of spectral balancing between single-frequency data volume |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610576825.7A CN106291686B (en) | 2016-07-20 | 2016-07-20 | The method and device of spectral balancing between single-frequency data volume |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106291686A true CN106291686A (en) | 2017-01-04 |
CN106291686B CN106291686B (en) | 2019-02-15 |
Family
ID=57651773
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610576825.7A Active CN106291686B (en) | 2016-07-20 | 2016-07-20 | The method and device of spectral balancing between single-frequency data volume |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106291686B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107505665A (en) * | 2017-07-07 | 2017-12-22 | 陕西省煤田地质集团有限公司 | A kind of stratigraphic anormaly earthquake detection method based on window Fourier transform |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080232193A1 (en) * | 2007-03-20 | 2008-09-25 | Geocyber Solutions, Inc. | Methods for Noise Removal and/or Attenuation from Seismic Data by Wavelet Selection |
CN102478667A (en) * | 2010-11-30 | 2012-05-30 | 中国石油天然气集团公司 | Inversion method of frequency dispersion amplitude versus offset (AVO) |
CN103207408A (en) * | 2013-03-01 | 2013-07-17 | 中煤科工集团西安研究院 | Passive earthquake monitoring data compression method and control system |
CN103713321A (en) * | 2014-01-08 | 2014-04-09 | 孙赞东 | Crack fluid type identifying method based on longitudinal wave frequency depending on amplitude versus offset (AVO) and azimuth |
CN104216008A (en) * | 2013-06-05 | 2014-12-17 | 中国石油天然气集团公司 | Downhole fracturing microseismic event identification method |
-
2016
- 2016-07-20 CN CN201610576825.7A patent/CN106291686B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080232193A1 (en) * | 2007-03-20 | 2008-09-25 | Geocyber Solutions, Inc. | Methods for Noise Removal and/or Attenuation from Seismic Data by Wavelet Selection |
CN102478667A (en) * | 2010-11-30 | 2012-05-30 | 中国石油天然气集团公司 | Inversion method of frequency dispersion amplitude versus offset (AVO) |
CN103207408A (en) * | 2013-03-01 | 2013-07-17 | 中煤科工集团西安研究院 | Passive earthquake monitoring data compression method and control system |
CN104216008A (en) * | 2013-06-05 | 2014-12-17 | 中国石油天然气集团公司 | Downhole fracturing microseismic event identification method |
CN103713321A (en) * | 2014-01-08 | 2014-04-09 | 孙赞东 | Crack fluid type identifying method based on longitudinal wave frequency depending on amplitude versus offset (AVO) and azimuth |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107505665A (en) * | 2017-07-07 | 2017-12-22 | 陕西省煤田地质集团有限公司 | A kind of stratigraphic anormaly earthquake detection method based on window Fourier transform |
Also Published As
Publication number | Publication date |
---|---|
CN106291686B (en) | 2019-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20150058388A1 (en) | Fractional scaling digital filters and the generation of standardized noise and synthetic data series | |
CN106483563B (en) | Seismic energy compensation method based on complementary set empirical mode decomposition | |
Li et al. | Scalar seismic-wave equation modeling by a multisymplectic discrete singular convolution differentiator method | |
CN104849757B (en) | Eliminate random noise system and method in seismic signal | |
CN106646615A (en) | Data processing method and data processing device for surface wave dispersion curve | |
CN104502965A (en) | Retrieving method for amplitude compensation factor | |
CN106526665B (en) | A kind of method and device of random noise attenuation | |
CN113341455A (en) | Viscous anisotropic medium seismic wave numerical simulation method, device and equipment | |
CN106680874A (en) | Harmonic noise suppression method based on waveform morphology sparse modeling | |
CN105425298B (en) | The method and apparatus of numerical solidification during a kind of elimination finite difference forward modeling | |
CN105301572A (en) | Multi-peak low sidelobe phase encode signal and mismatched filter joint design method | |
CN105044772A (en) | Method and device for realizing time varying spectrum simulation deconvolution | |
CN105550998B (en) | Image enchancing method and Image Intensified System based on two generation small echo integer transforms | |
Shahjouei et al. | Synthetic seismograms using a hybrid broadband ground‐motion simulation approach: Application to central and eastern United States | |
Banjade et al. | Earthquake accelerogram denoising by wavelet-based variational mode decomposition | |
CN104635264B (en) | The processing method of earthquake data before superposition and equipment | |
CN106291686A (en) | The method and device of spectral balancing between single-frequency data volume | |
CN106125134A (en) | Based on the geological data signal-noise ratio computation method of window during hyperbolic | |
CN114441897A (en) | Method for identifying arrival time difference of partial discharge pulse of distribution cable line | |
CN104267430B (en) | Determine the method and device of the earthquake fluid sensitive factor | |
Kassakian | Convex approximation and optimization with applications in magnitude filter design and radiation pattern synthesis | |
CN107219551A (en) | Widen the method and device of geological data frequency band | |
CN106094020B (en) | A kind of seismic inversion method and device | |
CN106249282A (en) | A kind of frequency domain seismic channel set creation method being applicable to AVAF inverting | |
CN109975867B (en) | Frequency extension method for seismic data with frequency domain signal aliasing |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |