CN107817517A - Trace gather preprocess method and device for azimuthal AVO inverting - Google Patents
Trace gather preprocess method and device for azimuthal AVO inverting Download PDFInfo
- Publication number
- CN107817517A CN107817517A CN201610819769.5A CN201610819769A CN107817517A CN 107817517 A CN107817517 A CN 107817517A CN 201610819769 A CN201610819769 A CN 201610819769A CN 107817517 A CN107817517 A CN 107817517A
- Authority
- CN
- China
- Prior art keywords
- trace gather
- angle
- interpolation
- neighborhood
- small
- 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
-
- 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
- G01V1/32—Transforming one recording into another or one representation into another
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/626—Physical property of subsurface with anisotropy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/646—Fractures
Abstract
The invention provides a kind of trace gather preprocess method for azimuthal AVO inverting, it is characterised in that the described method comprises the following steps:Input reference azimuth angle gathers (A1,A2,…,Am);By the reference azimuth angle gathers (A1,A2,…,Am) be converted to neighborhood trace gather (G1,G2,…,Gn);In the neighborhood trace gather (G1,G2,…,Gn) enterprising row interpolation;Neighborhood trace gather after interpolation is converted into orientation angle gathers, obtains the orientation angle gathers (C of interpolated pretreatment1,C2,…,Cm).The trace gather preprocess method and device for AVO invertings of the present invention, solve the problems, such as to have airway in reference azimuth trace gather, lack road, the input trace gather of high quality is provided for azimuthal AVO inverting.
Description
Technical field
The present invention relates to geophysical prospecting for oil technical field, in particular it relates to a kind of road for azimuthal AVO inverting
Collect preprocess method and device.
Background technology
With the application of wide-azimuth (or comprehensive) seismic prospecting, how to be excavated from Omnibearing earthquake auto data richer
Rich underground medium information turns into an important study hotspot.Because fracture-type reservoir is commonly present obvious HTI anisotropy
Influence, be mainly shown as that amplitude, speed, reflection configuration and phase change with survey line Orientation differences, therefore fold in wide-azimuth
Crack can be predicted in preceding logging data application with the change in orientation using compressional wave information.Such as when being travelled using orientation, orientation declines
Subtract, azimuthal velocity (VVAZ, speed is with azimuthal variation) and orientation amplitude variations (AVAZ, amplitude is with azimuthal variation) etc. come it is pre-
Survey crack.In terms of using compressional wave azimuthal anisotropy feature prediction crack, Chinese scholar has been carried out more application to this and ground
Study carefully, such as Cheng Bingjie is based on wide-azimuth conversion wave seismic data using converted wave AVAZ crack detection methods to the western deeper clefts in river
Property gas reservoir has carried out Crack Detection;Sun Wei for the seismic properties of compressional wave azimuthal anisotropy FRACTURE PREDICTION to preferably being ground
Study carefully, and utilize the reliability and accuracy of prediction result of the crack Information Authentication on imaging logging and rock core information.It is wherein square
Position AVO inversion techniques make use of the azimuth characteristic of Omnibearing earthquake auto data, the inverting crack information from Omnibearing earthquake auto data, be
Crack elimination prediction provides important analysis technological means.
Inventor has found that the original trace gather for azimuthal AVO inverting generally has the problem of airway, the quality of trace gather compared with
Difference, carry out needing azimuthal trace gather to be pre-processed to avoid depositing under some azimuths or incidence angle before azimuthal AVO analysis
In the phenomenon of airway.Occur that the data in whole survey grid under a certain orientation and incidence angle are sky in real data pretreatment
Road, lacked data can not be also obtained even if by OVT domains interpolation, thus have impact on the application effect of AVO inversion techniques.Cause
This, it is necessary to propose a kind of method that high quality input trace gather is provided for azimuthal AVO inverting.
The information for being disclosed in background of invention part is merely intended to deepen the reason of the general background technology to the present invention
Solution, and be not construed as recognizing or imply known to those skilled in the art existing of the information structure in any form
Technology.
The content of the invention
The purpose of the present invention is that the input trace gather of azimuthal AVO inverting is pre-processed, and solves sky in original trace gather be present
Road, road problem is lacked, providing high quality for azimuthal AVO inverting inputs trace gather.
According to an aspect of the invention, it is proposed that a kind of trace gather preprocess method for AVO invertings.This method can wrap
Include following steps:Input reference azimuth angle gathers (A1,A2,…,Am);By the reference azimuth angle gathers (A1,A2,…,Am) turn
It is changed to neighborhood trace gather (G1,G2,…,Gn);In the neighborhood trace gather (G1,G2,…,Gn) enterprising row interpolation;By the neighborhood after interpolation
Trace gather is converted to orientation angle gathers, obtains the orientation angle gathers (C of interpolated pretreatment1,C2,…,Cm), wherein, m represents orientation
Angle number, n represent incidence angle number.
Preferably, A1,A2,…,AmRespectively A1(φ1;θ1,θ2,…,θn),A2(φ2;θ1,θ2,…,θn),…,Am(φm;
θ1,θ2,…,θn), φ1,φ2,…,φmFor the azimuth arranged from small to large, θ1,θ2,…,θnEnter for what is arranged from small to large
Firing angle (A1,A2,…,Am)。
Preferably, by the orientation angle gathers (A1,A2,…,Am) be converted to neighborhood trace gather (G1,G2,…,Gn) include:
Incidence angle θ1,θ2,…,θnPlace obtains G according to the small seismic channel for arriving big order greatly and extracting different orientations successively1(θ1;φ1,
φ2,…,φm),G2(θ2;φ1,φ2,…,φm),…,Gn(θn;φ1,φ2,…,φm);By the seismic channel of all extractions according to
The order of incidence angle from small to large is arranged, and obtains the neighborhood trace gather (G1,G2,…,Gn)。
Preferably, in the neighborhood trace gather (G1,G2,…,Gn) interpolation method includes MINIMUM WEIGHT used by enterprising row interpolation
Norm interpolation method, the interpolation method based on Fourier transformation.
Preferably, the neighborhood trace gather after interpolation is converted into orientation angle gathers includes:In method angle φ1,φ2,…,φmPlace
Extract the seismic channel of different incidence angles successively according to order from small to large, obtain C1(φ1;θ1,θ2,…,θn),C2(φ2;θ1,
θ2,…,θn),…,Cm(φm;θ1,θ2,…,θn);The seismic channel of all extractions is carried out according to the order of azimuth from small to large
Arrangement, obtains the orientation angle gathers (C of interpolated pretreatment1,C2,…,Cm)。
According to another aspect of the invention, it is proposed that a kind of trace gather pretreatment unit for AVO invertings.The device can wrap
Include:For inputting reference azimuth angle gathers (A1,A2,…,Am) unit;For by the reference azimuth angle gathers (A1,A2,…,
Am) be converted to neighborhood trace gather (G1,G2,…,Gn) unit;For in the neighborhood trace gather (G1,G2,…,Gn) enterprising row interpolation
Unit;For the neighborhood trace gather after interpolation to be converted into orientation angle gathers, the orientation angle gathers of interpolated pretreatment are obtained
(C1,C2,…,Cm) unit, wherein, m represent azimuth number, n represent incidence angle number.
Preferably, A1,A2,…,AmRespectively A1(φ1;θ1,θ2,…,θn),A2(φ2;θ1,θ2,…,θn),…,Am(φm;
θ1,θ2,…,θn), φ1,φ2,…,φmFor the azimuth arranged from small to large, θ1,θ2,…,θnEnter for what is arranged from small to large
Firing angle (A1,A2,…,Am)。
Preferably, by the orientation angle gathers (A1,A2,…,Am) be converted to neighborhood trace gather (G1,G2,…,Gn) include:
Incidence angle θ1,θ2,…,θnPlace obtains G according to the small seismic channel for arriving big order greatly and extracting different orientations successively1(θ1;φ1,
φ2,…,φm),G2(θ2;φ1,φ2,…,φm),…,Gn(θn;φ1,φ2,…,φm);By the seismic channel of all extractions according to
The order of incidence angle from small to large is arranged, and obtains the neighborhood trace gather (G1,G2,…,Gn)。
Preferably, in the neighborhood trace gather (G1,G2,…,Gn) interpolation method includes MINIMUM WEIGHT used by enterprising row interpolation
Norm interpolation method, the interpolation method based on Fourier transformation.
Preferably, the neighborhood trace gather after interpolation is converted into orientation angle gathers includes:In method angle φ1,φ2,…,φmPlace
Extract the seismic channel of different incidence angles successively according to order from small to large, obtain C1(φ1;θ1,θ2,…,θn),C2(φ2;θ1,
θ2,…,θn),…,Cm(φm;θ1,θ2,…,θn);The seismic channel of all extractions is carried out according to the order of azimuth from small to large
Arrangement, obtains the orientation angle gathers (C of interpolated pretreatment1,C2,…,Cm)。
The reference azimuth angle gathers of input are transformed into neighborhood trace gather by the present invention, in the enterprising row interpolation of neighborhood trace gather, then
Neighborhood trace gather after interpolation is converted back into orientation angle gathers.Because neighborhood trace gather is higher compared with the lateral likeness of orientation angle gathers,
The more physical significance of Trace Interpolation is carried out in neighborhood trace gather.
Methods and apparatus of the present invention has other characteristics and advantage, and these characteristics and advantage are attached from what is be incorporated herein
It will be apparent in figure and subsequent specific embodiment, or by the accompanying drawing and subsequent specific implementation being incorporated herein
Stated in detail in example, these the drawings and specific embodiments are provided commonly for explaining the certain principles of the present invention.
Brief description of the drawings
Disclosure illustrative embodiments are described in more detail in conjunction with the accompanying drawings, the disclosure above-mentioned and its
Its purpose, feature and advantage will be apparent, wherein, in disclosure illustrative embodiments, identical reference number
Typically represent same parts.
Fig. 1 shows the stream of the trace gather preprocess method according to an embodiment of the invention for azimuthal AVO inverting
Cheng Tu.
Fig. 2 is the example orientations angle gathers as input data.
Fig. 3 is that field trace gather extracts schematic diagram.
Fig. 4 is exemplary fields trace gather.
Fig. 5 a and Fig. 5 b, which are that field trace gather interpolation is front and rear, to be contrasted.
Fig. 6 a and Fig. 6 b, which are that angle gathers pretreatment in orientation is front and rear, to be contrasted.
Embodiment
The preferred embodiment of the disclosure is more fully described below with reference to accompanying drawings.Although the disclosure is shown in accompanying drawing
Preferred embodiment, however, it is to be appreciated that may be realized in various forms the disclosure without the embodiment party that should be illustrated here
Formula is limited.On the contrary, these embodiments are provided so that the disclosure is more thorough and complete, and can be by the disclosure
Scope is intactly communicated to those skilled in the art.
Embodiment 1
Fig. 1 shows the stream of the trace gather preprocess method according to an embodiment of the invention for azimuthal AVO inverting
Cheng Tu.
In the present embodiment, this method specifically includes following steps:
1) reference azimuth angle gathers (A is inputted1,A2,…,Am)。
It is the characteristics of (inline, crossline) puts azimuthal trace gather:A corresponding Zu Jiao roads under each azimuth
Collection.Therefore, can be by A1,A2,…,AmIt is embodied as:A1(φ1;θ1,θ2,…,θn),A2(φ2;θ1,θ2,…,θn),…,Am
(φm;θ1,θ2,…,θn).Wherein, m represents azimuth number, φ1,φ2,…,φmFor the azimuth arranged from small to large;N tables
It is shown into firing angle number, θ1,θ2,…,θnFor the incidence angle arranged from small to large.
It will be appreciated by those skilled in the art that φ1,φ2,…,φmAnd θ1,θ2,…,θnCan also be according to from big to small
Order arranges.
2) by reference azimuth angle gathers (A1,A2,…,Am) be converted to neighborhood trace gather (G1,G2,…,Gn)。
In one example, reference azimuth angle gathers are converted to the method for neighborhood trace gather to be included:In incidence angle θ1,
θ2,…,θnPlace obtains G according to the small seismic channel for arriving big order greatly and extracting different orientations successively1(θ1;φ1,φ2,…,
φm),G2(θ2;φ1,φ2,…,φm),…,Gn(θn;φ1,φ2,…,φm);By the seismic channel of all extractions according to incidence angle
Order from small to large is arranged, and obtains the neighborhood trace gather (G1,G2,…,Gn)。
That is, first in minimum incidence angle θ1Azimuth φ is extracted successively in place1,φ2,…,φmSeismic channel, obtain G1(θ1;
φ1,φ2,…,φm);Then in next incidence angle θ2Azimuth φ is extracted successively in place1,φ2,…,φmSeismic channel, obtain
G2(θ2;φ1,φ2,…,φm);By that analogy, finally in maximum incident angle θnAzimuth φ is extracted successively in place1,φ2,…,φm
Seismic channel, obtain Gn(θn;φ1,φ2,…,φm).The seismic channel of all extractions is entered according to incidence angle is small to big order
Row arrangement, that is, obtained field trace gather (G1,G2,…,Gn)。
In reference azimuth angle gathers the arrangement of seismic channel be according under specified azimuth by the ascending arrangement of incidence angle, according to
Angular-trace gather corresponding to secondary all azimuths of arrangement.On space plane it is not between reference azimuth angular-trace gather Zhong Dao and road
Adjacent, similitude is poor, and the interpolation method of prior art requires that signal has the continuity on time or locus in addition, because
This is not suitable for carrying out on reference azimuth angle gathers.And after reference azimuth angle gathers are converted into neighborhood trace gather, in neighborhood trace gather
On space plane position it is adjacent between seismic channel, the similitude between Dao Yu roads is compared with reference azimuth angle gathers Zhong Dao and road
Between similitude be obviously improved, the enterprising row interpolation of neighborhood trace gather make use of seismic channel space plane be distributed position phase
Adjacent feature so that signal has continuity in space plane distribution, therefore has more physical significance compared with reference azimuth angle gathers interpolation,
It is also more suitable for being handled using the interpolation method of prior art.
3) in neighborhood trace gather (G1,G2,…,Gn) enterprising row interpolation.
In one example, MINIMUM WEIGHT norm interpolation method can be used or based on the interpolation method of Fourier transformation in neighborhood
Trace gather (G1,G2,…,Gn) enterprising row interpolation.It will be appreciated by those skilled in the art that any appropriate interpolation method can be used
In the enterprising row interpolation of neighborhood trace gather.
4) the neighborhood trace gather after interpolation is converted into orientation angle gathers, obtains the orientation angle gathers (C of interpolated pretreatment1,
C2,…,Cm)。
In one example, the neighborhood trace gather after interpolation being converted into orientation angle gathers can include:In method angle φ1,
φ2,…,φmPlace extracts the seismic channel of different incidence angles according to order from small to large successively, obtains C1(φ1;θ1,θ2,…,
θn),C2(φ2;θ1,θ2,…,θn),…,Cm(φm;θ1,θ2,…,θn);By the seismic channel of all extractions according to azimuth from it is small to
Big order is arranged, and obtains the orientation angle gathers (C of interpolated pretreatment1,C2,…,Cm)。
That is, first in minimum azimuth φ1Incidence angle θ is extracted successively in place1,θ2,…,θnSeismic channel, obtain C1(φ1;
θ1,θ2,…,θn);Then in next azimuth φ2Incidence angle θ is extracted successively in place1,θ2,…,θnSeismic channel, obtain C2
(φ2;θ1,θ2,…,θn);By that analogy, finally in maximum azimuth φmIncidence angle θ is extracted successively in place1,θ2,…,θnEarthquake
Road, obtain Cm(φm;θ1,θ2,…,θn).The seismic channel of all extractions is arranged according to azimuth is small to big order,
Neighborhood trace gather (C is obtained1,C2,…,Cm)。
Embodiment 2
According to another embodiment of the present invention, there is provided a kind of trace gather pretreatment unit for azimuthal AVO inverting.The dress
Putting to include:For inputting reference azimuth angle gathers (A1,A2,…,Am) unit;For by the reference azimuth angle gathers
(A1,A2,…,Am) be converted to neighborhood trace gather (G1,G2,…,Gn) unit;For in the neighborhood trace gather (G1,G2,…,Gn) on
Enter the unit of row interpolation;For the neighborhood trace gather after interpolation to be converted into orientation angle gathers, the orientation of interpolated pretreatment is obtained
Angle gathers (C1,C2,…,Cm) unit.Wherein, m represents azimuth number, and n represents incidence angle number.
In one example, A1,A2,…,AmCan be respectively A1(φ1;θ1,θ2,…,θn),A2(φ2;θ1,θ2,…,
θn),…,Am(φm;θ1,θ2,…,θn), φ1,φ2,…,φmFor the azimuth arranged from small to large, θ1,θ2,…,θnFor from small
To the incidence angle (A of longer spread1,A2,…,Am)。
In one example, by the orientation angle gathers (A1,A2,…,Am) be converted to neighborhood trace gather (G1,G2,…,Gn) can
With including:In incidence angle θ1,θ2,…,θnPlace obtains according to the small seismic channel for arriving big order greatly and extracting different orientations successively
G1(θ1;φ1,φ2,…,φm),G2(θ2;φ1,φ2,…,φm),…,Gn(θn;φ1,φ2,…,φm);By all extractions
Seismic channel is arranged according to the order of incidence angle from small to large, obtains the neighborhood trace gather (G1,G2,…,Gn)。
In one example, in the neighborhood trace gather (G1,G2,…,Gn) interpolation method can be with used by enterprising row interpolation
Including MINIMUM WEIGHT norm interpolation method, the interpolation method based on Fourier transformation.
In one example, the neighborhood trace gather after interpolation being converted into orientation angle gathers can include:In method angle φ1,
φ2,…,φmPlace extracts the seismic channel of different incidence angles according to order from small to large successively, obtains C1(φ1;θ1,θ2,…,
θn),C2(φ2;θ1,θ2,…,θn),…,Cm(φm;θ1,θ2,…,θn);By the seismic channel of all extractions according to azimuth from it is small to
Big order is arranged, and obtains the orientation angle gathers (C of interpolated pretreatment1,C2,…,Cm)。
Using example
For ease of understanding the scheme of the embodiment of the present invention and its effect, a concrete application example given below.This area
It should be understood to the one skilled in the art that the example, only for the purposes of understanding the present invention, its any detail is not intended to be limited in any way
The system present invention.
Illustrate the effect of the present invention with the actual seismic data instance in somewhere below.
Fig. 2 is the exemplary reference azimuth angle gathers as input data.As shown in Fig. 2 azimuth (0 °, 36 °, 72 °,
36 °, 108 °, 144 °) in the corresponding one group of angle gathers being made up of different incidence angles in each azimuth, incidence angle is ascending
Arrangement, the reference azimuth angle gathers inputted are the set of this five groups of angle gathers.As seen from the figure, in the reference azimuth angle gathers
Airway be present under some incidence angles.
Fig. 3 is that neighborhood trace gather extracts schematic diagram.As shown in figure 3, on reference azimuth angle gathers, opened at minimum incidence angle
Begin, all azimuthal seismic channels corresponding to the incidence angle are extracted successively according to order from small to large, in next incidence
All azimuthal seismic channels corresponding to it are gathered successively also according to order from small to large at angle, until having extracted maximum
All azimuthal seismic channels corresponding to incidence angle.Due to be according to incidence angle from small to large, azimuth from small to large suitable
Sequence is extracted successively, thus is adjacent between seismic channel.
Fig. 4 is exemplary neighborhood trace gather, is taken out for Fig. 2 reference azimuth angle gathers according to the extraction mode shown in Fig. 3
The result taken.From fig. 4, it can be seen that in the neighborhood trace gather, existing airway is more concentrated compared with the airway of reference azimuth trace gather.
Fig. 5 a and Fig. 5 b, which are that neighborhood trace gather interpolation is front and rear, to be contrasted.Wherein, Fig. 5 a are the neighborhood trace gather before interpolation, wherein free
Road;Fig. 5 b are the neighborhood trace gather after interpolation.
Fig. 6 a and Fig. 6 b, which are that angle gathers pretreatment in orientation is front and rear, to be contrasted.Wherein, Fig. 6 a are the original-party parallactic angle before pretreatment
Trace gather, many airways be present under some of which incidence angle;Fig. 6 b are the orientation angle gathers by pretreatment.It is it can be seen that of the invention
For the trace gather preprocess method and device of AVO invertings, solve the problems, such as to have airway in reference azimuth trace gather, lack road, be
Azimuthal AVO inverting provides the input trace gather of high quality.
It is described above the presently disclosed embodiments, described above is exemplary, and non-exclusive, and
It is not limited to disclosed each embodiment.In the case of without departing from the scope and spirit of illustrated each embodiment, for this skill
Many modifications and changes will be apparent from for the those of ordinary skill in art field.The selection of term used herein, purport
The principle of each embodiment, practical application or improvement to the technology in market are best being explained, or is making the art
Other those of ordinary skill are understood that each embodiment disclosed herein.
Claims (10)
1. a kind of trace gather preprocess method for azimuthal AVO inverting, it is characterised in that the described method comprises the following steps:
Input reference azimuth angle gathers (A1,A2,…,Am);
By the reference azimuth angle gathers (A1,A2,…,Am) be converted to neighborhood trace gather (G1,G2,…,Gn);
In the neighborhood trace gather (G1,G2,…,Gn) enterprising row interpolation;
Neighborhood trace gather after interpolation is converted into orientation angle gathers, obtains the orientation angle gathers (C of interpolated pretreatment1,C2,…,
Cm),
Wherein, m represents azimuth number, and n represents incidence angle number.
2. the trace gather preprocess method according to claim 1 for azimuthal AVO inverting, it is characterised in that A1,A2,…,
AmRespectively A1(φ1;θ1,θ2,…,θn),A2(φ2;θ1,θ2,…,θn),…,Am(φm;θ1,θ2,…,θn), φ1,φ2,…,
φmFor the azimuth arranged from small to large, θ1,θ2,…,θnFor the incidence angle (A arranged from small to large1,A2,…,Am)。
3. the trace gather preprocess method according to claim 1 for azimuthal AVO inverting, it is characterised in that by the side
Parallactic angle trace gather (A1,A2,…,Am) be converted to neighborhood trace gather (G1,G2,…,Gn) include:
In incidence angle θ1,θ2,…,θnPlace obtains G according to the small seismic channel for arriving big order greatly and extracting different orientations successively1
(θ1;φ1,φ2,…,φm),G2(θ2;φ1,φ2,…,φm),…,Gn(θn;φ1,φ2,…,φm);
The seismic channel of all extractions is arranged according to the order of incidence angle from small to large, obtains the neighborhood trace gather (G1,
G2,…,Gn)。
4. the trace gather preprocess method according to claim 1 for azimuthal AVO inverting, it is characterised in that in the neighbour
Domain trace gather (G1,G2,…,Gn) interpolation method is included MINIMUM WEIGHT norm interpolation method, become based on Fourier used by enterprising row interpolation
The interpolation method changed.
5. the trace gather preprocess method according to claim 1 for azimuthal AVO inverting, it is characterised in that after interpolation
Neighborhood trace gather be converted to orientation angle gathers and include:
In method angle φ1,φ2,…,φmPlace extracts the seismic channel of different incidence angles according to order from small to large successively, obtains C1
(φ1;θ1,θ2,…,θn),C2(φ2;θ1,θ2,…,θn),…,Cm(φm;θ1,θ2,…,θn);
The seismic channel of all extractions is arranged according to the order of azimuth from small to large, obtains the orientation of interpolated pretreatment
Angle gathers (C1,C2,…,Cm)。
6. a kind of trace gather pretreatment unit for azimuthal AVO inverting, it is characterised in that described device includes:
For inputting reference azimuth angle gathers (A1,A2,…,Am) unit;
For by the reference azimuth angle gathers (A1,A2,…,Am) be converted to neighborhood trace gather (G1,G2,…,Gn) unit;
For in the neighborhood trace gather (G1,G2,…,Gn) enterprising row interpolation unit;
For the neighborhood trace gather after interpolation to be converted into orientation angle gathers, the orientation angle gathers (C of interpolated pretreatment is obtained1,
C2,…,Cm) unit,
Wherein, m represents azimuth number, and n represents incidence angle number.
7. the trace gather pretreatment unit according to claim 6 for azimuthal AVO inverting, it is characterised in that A1,A2,…,
AmRespectively A1(φ1;θ1,θ2,…,θn),A2(φ2;θ1,θ2,…,θn),…,Am(φm;θ1,θ2,…,θn), φ1,φ2,…,
φmFor the azimuth arranged from small to large, θ1,θ2,…,θnFor the incidence angle (A arranged from small to large1,A2,…,Am)。
8. the trace gather pretreatment unit according to claim 6 for azimuthal AVO inverting, it is characterised in that by the side
Parallactic angle trace gather (A1,A2,…,Am) be converted to neighborhood trace gather (G1,G2,…,Gn) include:
In incidence angle θ1,θ2,…,θnPlace obtains G according to the small seismic channel for arriving big order greatly and extracting different orientations successively1
(θ1;φ1,φ2,…,φm),G2(θ2;φ1,φ2,…,φm),…,Gn(θn;φ1,φ2,…,φm);
The seismic channel of all extractions is arranged according to the order of incidence angle from small to large, obtains the neighborhood trace gather (G1,
G2,…,Gn)。
9. the trace gather pretreatment unit according to claim 6 for azimuthal AVO inverting, it is characterised in that in the neighbour
Domain trace gather (G1,G2,…,Gn) interpolation method is included MINIMUM WEIGHT norm interpolation method, become based on Fourier used by enterprising row interpolation
The interpolation method changed.
10. the trace gather pretreatment unit according to claim 6 for azimuthal AVO inverting, it is characterised in that after interpolation
Neighborhood trace gather be converted to orientation angle gathers, obtaining the orientation angle gathers of interpolated pretreatment includes:
In method angle φ1,φ2,…,φmPlace extracts the seismic channel of different incidence angles according to order from small to large successively, obtains C1
(φ1;θ1,θ2,…,θn),C2(φ2;θ1,θ2,…,θn),…,Cm(φm;θ1,θ2,…,θn);
The seismic channel of all extractions is arranged according to the order of azimuth from small to large, obtains the orientation of interpolated pretreatment
Angle gathers (C1,C2,…,Cm)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610819769.5A CN107817517B (en) | 2016-09-12 | 2016-09-12 | Trace gather preprocess method and device for azimuthal AVO inverting |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610819769.5A CN107817517B (en) | 2016-09-12 | 2016-09-12 | Trace gather preprocess method and device for azimuthal AVO inverting |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107817517A true CN107817517A (en) | 2018-03-20 |
CN107817517B CN107817517B (en) | 2019-10-08 |
Family
ID=61600749
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610819769.5A Active CN107817517B (en) | 2016-09-12 | 2016-09-12 | Trace gather preprocess method and device for azimuthal AVO inverting |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107817517B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113917533A (en) * | 2020-07-10 | 2022-01-11 | 中国石油化工股份有限公司 | Systematic implementation method of double-linkage omnibearing imaging of TI medium |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104297800A (en) * | 2014-10-09 | 2015-01-21 | 中国海洋石油总公司 | Self-phase-control prestack inversion method |
US20150098627A1 (en) * | 2013-10-08 | 2015-04-09 | Shin-Ju Chu Ye | Automatic Dip Picking From Wellbore Azimuthal Image Logs |
CN104678434A (en) * | 2013-12-02 | 2015-06-03 | 中国石油化工股份有限公司 | Method for predicting storage layer crack development parameters |
WO2015145257A2 (en) * | 2014-03-28 | 2015-10-01 | Cgg Services Sa | Method and apparatus for characterizing azimuthal anisotropy using statistical moments |
CN105487111A (en) * | 2014-09-18 | 2016-04-13 | 中国石油化工股份有限公司 | Automatic azimuth gather partial stacking method |
-
2016
- 2016-09-12 CN CN201610819769.5A patent/CN107817517B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150098627A1 (en) * | 2013-10-08 | 2015-04-09 | Shin-Ju Chu Ye | Automatic Dip Picking From Wellbore Azimuthal Image Logs |
CN104678434A (en) * | 2013-12-02 | 2015-06-03 | 中国石油化工股份有限公司 | Method for predicting storage layer crack development parameters |
WO2015145257A2 (en) * | 2014-03-28 | 2015-10-01 | Cgg Services Sa | Method and apparatus for characterizing azimuthal anisotropy using statistical moments |
CN105487111A (en) * | 2014-09-18 | 2016-04-13 | 中国石油化工股份有限公司 | Automatic azimuth gather partial stacking method |
CN104297800A (en) * | 2014-10-09 | 2015-01-21 | 中国海洋石油总公司 | Self-phase-control prestack inversion method |
Non-Patent Citations (1)
Title |
---|
刘依谋 等: ""宽方位地震勘探技术新进展"", 《石油地球物理勘探》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113917533A (en) * | 2020-07-10 | 2022-01-11 | 中国石油化工股份有限公司 | Systematic implementation method of double-linkage omnibearing imaging of TI medium |
CN113917533B (en) * | 2020-07-10 | 2023-04-28 | 中国石油化工股份有限公司 | TI medium double-linkage omnibearing imaging systematic realization method |
Also Published As
Publication number | Publication date |
---|---|
CN107817517B (en) | 2019-10-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wu et al. | Seismic tomography of Taiwan: Improved constraints from a dense network of strong motion stations | |
CN102221709B (en) | Velocity analysis and dynamic correction method based on formation parameter information | |
CN105301641B (en) | A kind of method and device of azimuthal anisotropy velocity inversion | |
CN102053261B (en) | Method for processing seismic data | |
CN103713323B (en) | Omnibearing aeolotropy amplitude-preservation imaging and gather extracting method | |
Li et al. | Crustal thickness variation in the northern Appalachian Mountains: Implications for the geometry of 3‐D tectonic boundaries within the crust | |
CN106556861B (en) | A kind of azimuthal AVO inversion method based on Omnibearing earthquake auto data | |
CN102401908B (en) | Method for suppressing multiple waves by the aid of different-mode weighting sparse parabola Radon transform | |
CN105093292A (en) | Data processing method and device for earthquake imaging | |
CN103675897A (en) | Seismic diffracted wave separating and imaging method | |
CN106154330B (en) | Electromagnetic type controlled source parabolic Radon transform suppresses parallel focus crosstalk noise method | |
Long et al. | Seismic imaging of mid-crustal structure beneath central and eastern North America: Possibly the elusive Grenville deformation? | |
Amoroso et al. | S‐wave identification by polarization filtering and waveform coherence analyses | |
CN104932010A (en) | Diffraction wave separating method based on shortcut edging sparse Radon transformation | |
CN103454681B (en) | Evaluate the method and apparatus of 3 D seismic observation system imaging effect | |
CN104391324A (en) | Seismic trace set dynamic correction stretching correction pre-processing technology before AVO inversion depending on frequency | |
CN107807393A (en) | Separate unit station collection preliminary wave Enhancement Method based on seismic interference method | |
CN107861156A (en) | The extracting method and device of diffracted wave | |
CN102866425A (en) | Blind source seismic signal stable-superposition model based blind separation method | |
Bulut et al. | Spatiotemporal earthquake clusters along the North Anatolian fault zone offshore Istanbul | |
CN107817517B (en) | Trace gather preprocess method and device for azimuthal AVO inverting | |
CN107340537A (en) | A kind of method of P-SV converted waves prestack reverse-time depth migration | |
CN105866839A (en) | Static correction method and static correction device on basis of common-attitude trace gathers | |
CN104422955A (en) | Method for using travel-time variation quantity to extract anisotropism parameters | |
CN106526664A (en) | Time-space domain geometric transform horizontal stacking seismic imaging method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |