CN106125138B - Pre-stack CRP gather frequency broadening processing method - Google Patents
Pre-stack CRP gather frequency broadening processing method Download PDFInfo
- Publication number
- CN106125138B CN106125138B CN201610428193.XA CN201610428193A CN106125138B CN 106125138 B CN106125138 B CN 106125138B CN 201610428193 A CN201610428193 A CN 201610428193A CN 106125138 B CN106125138 B CN 106125138B
- Authority
- CN
- China
- Prior art keywords
- frequency
- crp
- seismic
- energy
- original
- 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
Links
- 238000003672 processing method Methods 0.000 title description 4
- 238000012545 processing Methods 0.000 claims abstract description 21
- 238000000034 method Methods 0.000 claims abstract description 18
- 238000001228 spectrum Methods 0.000 abstract description 5
- 230000000694 effects Effects 0.000 abstract description 4
- 239000011159 matrix material Substances 0.000 description 14
- 238000005516 engineering process Methods 0.000 description 8
- 238000011160 research Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 241001062009 Indigofera Species 0.000 description 1
- 241000219000 Populus Species 0.000 description 1
- 238000007630 basic procedure Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
- 230000002087 whitening effect Effects 0.000 description 1
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. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
- G01V2210/322—Trace stacking
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
A method for frequency extension processing of pre-stack CRP gather. According to the frequency spectrum characteristics of seismic waves, the resolution ratio on a seismic time section is improved, and meanwhile, a good signal-to-noise ratio is kept, so that the seismic time section and the signal-to-noise ratio are combined. The method comprises the following steps: the method comprises the following steps: inverting the target interval by using the logging data and the well-side seismic channel; step two: calculating seismic wavelets according to the number of the destination layers in the step one; step three: calculating the energy and frequency distribution range of the original CRP gather; step four: obtaining a noise type according to the energy and frequency distribution range of the original CRP gather; step five: selecting a typical channel and determining a main energy distribution space of the effective wave; step six: determining a CRP frequency extension operator according to the energy, the frequency distribution range and the main energy distribution space of the active wave of the original CRP gather; step seven: and performing frequency extension processing on the CRP gather according to the seismic wavelets, the noise types and the CRP frequency extension operator. The signal-to-noise ratio with good signal-to-noise ratio can be kept; high reliability, good stability and very obvious effect.
Description
Technical field
The invention belongs to oil field prospecting field, is related to a kind of trace gather processing method, especially CRP trace gathers and opens up frequency processing side
Method, the method that specifically a kind of prestack CRP trace gathers open up frequency processing.
Background technology
With the raising of Periphery Oilfield exploration and development degree, seism processing faces object and become increasingly complex, especially
FY oil-layer, oil water relation is complicated, sand-body distribution is scattered, thickness of thin, is not responded significantly on seismic data, especially T2
Reflecting interface screen " covering effect " should be especially prominent and be also difficult to the spatial feature for portraying thin reservoir in terms of seismic properties,
How to improve the earthquake prediction precision of this kind of reservoir is always " bottleneck " problem that In Oil Field Exploration And Development faces, for this problem,
Many association areas have carried out research of technique, the tackling key problem particularly in terms of seismic data processing technology and reservoir prediction is opened in succession
Exhibition, therefore the research work in terms of the quality of pre stack data quality directly affects latter earthquake attribute and reservoir prediction;
New tool at home and abroad oil field popularization and application of the prestack inversion as reservoir prediction, the reason of numerous prestack inversions
It is mainly the various different approximate expressions to Zoeppritz equation groups by method, its basic procedure is to utilize CRP trace gather numbers
According to, the CRP trace gather data of high quality can not be obtained using conventional processing technique, how to improve the resolution ratio of Prestack seismic data into
The problem of being keen to discuss for current people, the resolution ratio of seismic data is low, and the insufficient or scarce of radio-frequency component is shown as in frequency domain
Lose;Although " spectral whitening " processing and " correction of indigo plant spectrum " processing method for having grown up can preferably solve seismic data
In " effective frequency range " the problem of some radio-frequency component deficiencies, but those frequency contents beyond the frequency range have to pass through certain
Method can just compensate;
Traditional CRP trace gather treatment technologies be can not solve at all on time section of earthquake is kept good signal to noise ratio and
Resolution ratio, the processing of traditional CRP trace gathers are the single signal to noise ratio improved on time section of earthquake, and this way is to have lost ground
Many effective radio-frequency components of seismic wave, the result obtained by traditional method is that signal to noise ratio increases, but signal to noise ratio is big
Amplitude reduction.
The content of the invention
In order to solve the above problems, the present invention, which proposes, a kind of can improve earthquake according to the spectral characteristic of seismic wave
Good signal to noise ratio is kept on time section while resolution ratio so that the two processing method be combineding with each other, its particular technique
Scheme is as follows:
A kind of method that prestack CRP trace gathers open up frequency processing, comprises the following steps:
Step 1:Utilize well-log information and seismic trace near well, inverting interval of interest;
Step 2:The purpose number of plies in step 1 calculates seismic wavelet;
Step 3:Calculate energy, the frequency distribution scope of original CRP trace gathers;
Step 4:According to the energy of original CRP trace gathers, frequency distribution scope, noise types are drawn;
Step 5:Typical track is selected, determines significant wave main energetic distribution space;
Step 6:According to the energy of original CRP trace gathers, frequency distribution scope and significant wave main energetic distribution space, really
Determine CRP and open up frequency operator;
Step 7:Frequency operator is opened up according to seismic wavelet, noise types and CRP, carries out opening up frequency processing in CRP trace gathers.
Beneficial effects of the present invention:Carry out opening up frequency processing present invention is generally directed to prestack CRP trace gathers, by using this method
The CRP trace gathers resolution ratio and signal to noise ratio obtained is all increased substantially, and the various seismic properties for prestack inversion and poststack are ground
Study carefully and true and reliable data volume is provided;Solve that both the upper signal to noise ratio of traditional CRP trace gathers processing and resolution ratio mutually restricts asks
Topic;The good signal to noise ratio of signal to noise ratio can be kept;Reliability is high, stability is good and effect is very notable.
Brief description of the drawings
Fig. 1 is the original CRP trace gathers of example;
Fig. 2 is the CRP trace gathers after opening up frequency of example;
Fig. 3 is the original CRP trace gathers stack result of example;
Fig. 4 is the CRP trace gather stack results after opening up frequency of example.
Embodiment
In order to make it easy to understand, the present invention is further detailed below:
Embodiment 1:A kind of method that prestack CRP trace gathers open up frequency processing, comprises the following steps:
Step 1:Utilize well-log information and seismic trace near well, inverting interval of interest;
Step 2:The purpose number of plies in step 1 calculates seismic wavelet;
Step 3:Calculate energy, the frequency distribution scope of original CRP trace gathers;
Step 4:According to the energy of original CRP trace gathers, frequency distribution scope, noise types are drawn;
Step 5:Typical track is selected, determines significant wave main energetic distribution space;
Step 6:According to the energy of original CRP trace gathers, frequency distribution scope and significant wave main energetic distribution space, really
Determine CRP and open up frequency operator;
Step 7:Frequency operator is opened up according to seismic wavelet, noise types and CRP, carries out opening up frequency processing in CRP trace gathers.
It is the difference of foundation useful signal and noise in frequency spectrum that CRP trace gathers, which open up frequency, and the difference on spatial coherence is carried out
SNR estimation and compensation, the purpose that CRP trace gathers open up frequency is exactly that useful signal is strengthened in the range of effective frequency band distribution, simultaneously
Abate the noise, high quality CRP trace gather data are provided for prestack inversion and poststack attributive analysis;Influence the weight of CRP trace gather quality
It is noise to want factor, therefore need to carry out CRP trace gather spectrum sigtral responses, and noise and useful signal have significantly on spectrum signature
Difference, the distribution of noise is found out by analysis;Original CRP trace gathers analysis is mainly from noise, energy, frequency etc.
Start with, judge whether original CRP trace gathers meet the requirement of reservoir prediction and explanation;
If earthquake record has M road, there is N number of sampled point in each road, and the data of this earthquake record can use a matrix X tables
Show:
X=[xij], i=1,2 ..., M;J=1,2 ..., N
Wherein i represents Taoist monastic name, and j represents sampled point sequence number;According to matrix theory, any one M × N matrix can be made just
Hand over and decompose, i.e.,:
X=U ∑s VT (1-1)
U and V is respectively M × M and N × N rank orthogonal matrixes, and ∑ is M × N matrix;If all row of matrix U are by covariance matrix
XXTEach characteristic vector (M dimensions) formed, V is by XTX each characteristic vector (N-dimensional) composition, then ∑ is M × N diagonal matrixs, main diagonal
Element by X singular value σi(XXTOr XTX ith feature value λiNon-negative square root) composition, off-diagonal element is all zero, this
When (2-1) formula be referred to as matrix X singular value decomposition;
(1-1) formula can also
Here r be matrix X order, uiFor XXTCorrespond to σiCharacteristic vector, viFor XTX corresponds to σiCharacteristic vector;A M × N matrix, referred to as X ith feature as;Obviously, matrix X singular value decomposition formula is not unique, in order to
Using upper convenience, the diagonal element of our regulation diagonal matrix ∑s is arranged by the order successively decreased, i.e.,:
σ1>σ2>…>σr>0
When M roads are separate, r=min (M, N), generally, r≤min (M, N);Matrix X Frobenius models
Number is defined as:
F is taken to be equal to 2, above formula is changed into:
Known by (1-4) formula, the F norms of matrix actually reflect the gross energy of M roads earthquake record (as F=2);Can be with
Prove:
This shows the gross energy of M roads earthquake record, the covariance matrix characteristic value (or square of X singular value) equal to X
Sum;
If M<N, when M Seismic Traces linear independence, σj(j=1,2 ..., M) is not zero, so X complete weight
Structure need all feature as;When X M Seismic Traces linear correlation, then only σ1It is not zero, X Perfect Reconstruction only needs
Want first feature picture;Seismic channel correlation is better, reconstructs the feature required for the record as fewer;
If only with preceding P feature as approximate reconstruction X, the earthquake record of reconstruction is XP, then:
XPGross energy be:
Reconstruction error is:
Wherein (1-6) is that CRP trace gathers open up frequency operator.
With going deep into for degree of prospecting, microstructures, lithological reservoir exploration turn into theme, help poplar river channel sand to identify
Problem is day by day highlighted, and the resolution ratio of common seismic processing is more inadequate, and earthquake is opened up frequency technology and is thus lifted to highly is carried out very much
Exploratory development, the present invention propose CRP trace gathers and open up frequency treatment technology, it is therefore an objective to for middle-shallow strata reservoir sand body scope is small, thing
Property change is fast, resource potential greatly, the characteristics of oil reservoir identification is difficult, clear the difficult point of High Resolution Seismic Technology, explore and protect panel height and divide
Resolution opens up frequency treatment technology, the identification of research thin interbed and description technique, improves recognition capability of the seismic technology to thin reservoir, is formed
Counter-techniques scheme, fine granularing scalability developmental research is serviced, the ground of high quality is provided for prestack AVO invertings, poststack attributes research
Shake data.
Specific embodiment:With reference to Fig. 1, in certain work area trace gather data, its Reservoir Section be mainly distributed on 1000ms-1200ms it
Between, in terms of original road set analysis, its vertical, horizontal energy coincidence is preferable, but the lineups layering near main purpose layer compared with
Difference, especially signal to noise ratio and resolution ratio have some random noises all than relatively low, these be also generally existing in trace gather phenomenon because
This, that improves trace gather quality is critical that lineups are as horizontal as possible and suppress noise as far as possible;It is special according to the trace gather of this area
Sign, window size and maximum time shift amount up and down when definition contrasts, it is that the horizontal key processing of trace gather can that suitable parameter adjusts repeatedly
To improve the resolution ratio of CRP trace gathers and signal to noise ratio to a certain extent.It is with the result CRP roads after present invention processing with reference to Fig. 2
Collection, and the contrast of original trace gather are differentiated, it is apparent that the feature of horizontal distribution is substantially presented in the reflected energy of prestack trace gather
Rate and signal to noise ratio are all greatly improved.
Beneficial effect:The good signal to noise ratio of signal to noise ratio can be kept, reliability is high, stability is good and effect is very notable.
It can keep good signal to noise ratio according to the spectral characteristic of seismic wave while resolution ratio on improving time section of earthquake, make
The two is obtained to be combined with each other.
Claims (1)
1. a kind of method that prestack CRP trace gathers open up frequency processing, it is characterised in that comprise the following steps:
Step 1:Utilize well-log information and seismic trace near well, the inverting purpose number of plies;
Step 2:The purpose number of plies in step 1 calculates seismic wavelet;
Step 3:Calculate energy, the frequency distribution scope of original CRP trace gathers;
Step 4:According to the energy of original CRP trace gathers, frequency distribution scope, noise types are drawn;
Step 5:Typical track is selected, determines significant wave main energetic distribution space;
Step 6:According to the energy of original CRP trace gathers, frequency distribution scope and significant wave main energetic distribution space, CRP is determined
Open up frequency operator;
Step 7:Frequency operator is opened up according to seismic wavelet, noise types and CRP, carries out opening up frequency processing in CRP trace gathers.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610428193.XA CN106125138B (en) | 2016-06-16 | 2016-06-16 | Pre-stack CRP gather frequency broadening processing method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610428193.XA CN106125138B (en) | 2016-06-16 | 2016-06-16 | Pre-stack CRP gather frequency broadening processing method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106125138A CN106125138A (en) | 2016-11-16 |
CN106125138B true CN106125138B (en) | 2018-03-16 |
Family
ID=57469603
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610428193.XA Active CN106125138B (en) | 2016-06-16 | 2016-06-16 | Pre-stack CRP gather frequency broadening processing method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106125138B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113325471B (en) * | 2021-05-21 | 2022-08-23 | 成都理工大学 | Seismic wave field subcomponent extraction method based on singular value decomposition |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU137397U1 (en) * | 2013-08-14 | 2014-02-10 | Джемма Павловна Земцова | SEISMIC EXPLORATION COMPLEX FOR SEARCHING HYDROCARBONS WITH THE MODULE FOR ANALYSIS OF LOW FREQUENCY RESONANCE EMISSION OF GEODYNAMIC NOISE |
CN104849756B (en) * | 2015-03-31 | 2018-04-27 | 中国地质大学(北京) | A kind of seismic data resolution that improves strengthens the method for effective weak signal energy |
-
2016
- 2016-06-16 CN CN201610428193.XA patent/CN106125138B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN106125138A (en) | 2016-11-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105334535B (en) | A kind of method for identifying the thin hidden lithologic deposit of reservoir | |
CN104516018B (en) | Porosity inversion method under lithological constraint in geophysical exploration | |
CN104570067B (en) | Phased seismic inversion method in a kind of geophysical exploration | |
CN107102356B (en) | Seismic signal high resolution data processing methods based on CEEMD | |
CN108802812A (en) | Well-seismic fusion stratum lithology inversion method | |
Du et al. | Seismic facies analysis based on self-organizing map and empirical mode decomposition | |
CN109633750A (en) | Nonlinear mapping relationship inversion method based on well logging phase wave impedance and seismic waveform | |
CN107356965B (en) | Reflection coefficient inverting method for predicting reservoir based on weighted superposition Noise Elimination strategy | |
Chopra et al. | Multispectral, multiazimuth, and multioffset coherence attribute applications | |
Li et al. | Seismic coherence for discontinuity interpretation | |
Xue et al. | Q-factor estimation by compensation of amplitude spectra in synchrosqueezed wavelet domain | |
Chopra et al. | Evolution of seismic interpretation during the last three decades | |
Cai et al. | Automated spectral recomposition with application in stratigraphic interpretation | |
Naeem et al. | Integrated study of seismic and well data for porosity estimation using multi-attribute transforms: a case study of Boonsville Field, Fort Worth Basin, Texas, USA | |
Shen et al. | Separation of diffracted waves via SVD filter | |
Lin et al. | Time-frequency mixed domain multi-trace simultaneous inversion method | |
CN106125138B (en) | Pre-stack CRP gather frequency broadening processing method | |
CN107678065B (en) | The guarantor for improving seismic resolution constructs well control space the Method of Deconvolution and device | |
Li et al. | Joint inversion of quality factor based on upgoing and downgoing waves | |
Vohs | 3D seismic attributes analysis in reservoir characterization: the Morrison NE field & Morrison field, Clark County Kansas | |
Mutlu et al. | Improving seismic resolution of prestack time-migrated data | |
Jiang et al. | Fast least-squares prestack time migration via accelerating the explicit calculation of Hessian matrix with dip-angle Fresnel zone | |
Tschache et al. | Estimation of net-to-gross ratio and net pay from seismic amplitude-variation-with-offset using Bayesian inversion | |
Li et al. | An efficient deep learning method for VSP wavefield separation: A DAS-VSP case | |
CN113311482B (en) | High-resolution medium-deep reservoir prediction method based on prestack spectrum inversion optimization |
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 |