CN107817517B - 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
- CN107817517B CN107817517B CN201610819769.5A CN201610819769A CN107817517B CN 107817517 B CN107817517 B CN 107817517B CN 201610819769 A CN201610819769 A CN 201610819769A CN 107817517 B CN107817517 B CN 107817517B
- Authority
- CN
- China
- Prior art keywords
- trace gather
- interpolation
- azimuth
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000000605 extraction Methods 0.000 claims description 15
- 239000000284 extract Substances 0.000 claims description 12
- 230000009466 transformation Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000010304 firing Methods 0.000 description 3
- 230000001174 ascending effect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000000151 deposition Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 239000011435 rock Substances 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
-
- 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/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
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
The present invention provides a kind of trace gather preprocess methods for azimuthal AVO inverting, which is characterized 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 interpolated pretreated orientation angle gathers (C1,C2,…,Cm).Trace gather preprocess method and device for AVO inverting of the invention solves the problems, such as to provide the input trace gather of high quality there are airway, scarce road in reference azimuth trace gather for azimuthal AVO inverting.
Description
Technical field
The present invention relates to geophysical prospecting for oil technical fields, and in particular, to a kind of road for azimuthal AVO inverting
Collect preprocess method and device.
Background technique
With the application of wide-azimuth (or comprehensive) seismic prospecting, how to be excavated from Omnibearing earthquake auto data richer
Rich underground medium information becomes an important research hotspot.Since fracture-type reservoir is commonly present apparent HTI anisotropy
Influence, be mainly shown as that amplitude, speed, reflection configuration and phase change with survey line Orientation differences, therefore folded in wide-azimuth
Longitudinal wave information be can use in preceding logging data application with the variation in orientation to predict crack.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 longitudinal wave azimuthal anisotropy feature prediction crack, Chinese scholar has been carried out more application to this and has been ground
Study carefully, such as Cheng Bingjie is based on wide-azimuth conversion wave seismic data using converted wave AVAZ crack detection method to the western deeper clefts in river
Property gas reservoir has carried out Crack Detection;Sun Wei preferably grinds the seismic properties for longitudinal wave azimuthal anisotropy FRACTURE PREDICTION
Study carefully, and utilizes the reliability and accuracy of prediction result of the crack Information Authentication on imaging logging and rock core information.Wherein side
The azimuth characteristic of Omnibearing earthquake auto data is utilized in position AVO inversion technique, and the inverting crack information from Omnibearing earthquake auto data is
Crack elimination prediction provides important analysis technological means.
Inventors have found that usually there is airway in the original trace gather for azimuthal AVO inverting, the quality of trace gather compared with
Difference carries out needing azimuthal trace gather to be pre-processed to avoid depositing under certain azimuths or incidence angle before azimuthal AVO analysis
The airway the phenomenon that.The data that will appear in entire survey grid under a certain orientation and incidence angle in real data pretreatment are sky
Road can not obtain lacked data passing through the domain OVT interpolation, thus affect the application effect of AVO inversion technique.Cause
This, it is necessary to propose it is a kind of for azimuthal AVO inverting provide high quality input trace gather method.
The information for being disclosed in background of invention part is merely intended to deepen the reason to general background technique of the invention
Solution, and it is known to those skilled in the art existing to be not construed as recognizing or imply that the information is constituted in any form
Technology.
Summary of the invention
The purpose of the present invention is pre-processing to the input trace gather of azimuthal AVO inverting, solve the presence of sky in original trace gather
Road lacks road problem, provides high quality for azimuthal AVO inverting and inputs trace gather.
According to an aspect of the invention, it is proposed that a kind of trace gather preprocess method for AVO inverting.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 interpolated pretreated orientation angle gathers (C1,C2,…,Cm), wherein m indicates orientation
Angle number, n indicate 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 sequence greatly and successively extracting different orientations1(θ1;φ1,
φ2,…,φm),G2(θ2;φ1,φ2,…,φm),…,Gn(θn;φ1,φ2,…,φm);By the seismic channel of all extractions according to
The sequence of incidence angle from small to large is arranged, and the neighborhood trace gather (G is obtained1,G2,…,Gn)。
Preferably, in the neighborhood trace gather (G1,G2,…,Gn) interpolation method used by enterprising row interpolation includes MINIMUM WEIGHT
Norm interpolation method, the interpolation method based on Fourier transformation.
Preferably, the neighborhood trace gather after interpolation is converted to orientation angle gathers includes: in azimuth φ1,φ2,…,φmPlace
The seismic channel that different incidence angles are successively extracted according to sequence from small to large, obtains 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 sequence of azimuth from small to large
Arrangement, obtains interpolated pretreated orientation angle gathers (C1,C2,…,Cm)。
According to another aspect of the invention, it is proposed that a kind of trace gather pretreatment unit for AVO inverting.The device can wrap
It 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 to orientation angle gathers, interpolated pretreated orientation angle gathers are obtained
(C1,C2,…,Cm) unit, wherein m indicate azimuth number, n indicate 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 sequence greatly and successively extracting different orientations1(θ1;φ1,
φ2,…,φm),G2(θ2;φ1,φ2,…,φm),…,Gn(θn;φ1,φ2,…,φm);By the seismic channel of all extractions according to
The sequence of incidence angle from small to large is arranged, and the neighborhood trace gather (G is obtained1,G2,…,Gn)。
Preferably, in the neighborhood trace gather (G1,G2,…,Gn) interpolation method used by enterprising row interpolation includes MINIMUM WEIGHT
Norm interpolation method, the interpolation method based on Fourier transformation.
Preferably, the neighborhood trace gather after interpolation is converted to orientation angle gathers includes: in azimuth φ1,φ2,…,φmPlace
The seismic channel that different incidence angles are successively extracted according to sequence from small to large, obtains 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 sequence of azimuth from small to large
Arrangement, obtains interpolated pretreated orientation angle gathers (C1,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.Since 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 advantages, these characteristics and advantages are attached from what is be incorporated herein
It will be apparent in figure and subsequent specific embodiment, or will be in the attached drawing and subsequent specific implementation being incorporated herein
It is stated in detail in example, these the drawings and specific embodiments are used together to explain specific principle of the invention.
Detailed description of the invention
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 label
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 are comparison before and after field trace gather interpolation.
Fig. 6 a and Fig. 6 b are orientation angle gathers pretreatment front and back comparison.
Specific embodiment
The preferred embodiment of the disclosure is more fully described below with reference to accompanying drawings.Although showing the disclosure in attached 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 this disclosure will be more thorough and complete, and can be by the disclosure
Range is completely 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 the following steps:
1) reference azimuth angle gathers (A is inputted1,A2,…,Am)。
The characteristics of (inline, crossline) puts azimuthal trace gather are as follows: under each azimuth a corresponding road Zu Jiao
Collection.It 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 indicates azimuth number, φ1,φ2,…,φmFor the azimuth arranged from small to large;N table
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,…,θnIt can also be according to from big to small
Sequence arranges.
2) by reference azimuth angle gathers (A1,A2,…,Am) be converted to neighborhood trace gather (G1,G2,…,Gn)。
It in one example, may include: in incidence angle θ by the method that reference azimuth angle gathers are converted to neighborhood trace gather1,
θ2,…,θnPlace obtains G according to the small seismic channel for arriving big sequence greatly and successively extracting different orientations1(θ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
Sequence from small to large is arranged, and the neighborhood trace gather (G is obtained1,G2,…,Gn)。
That is, first in minimum incidence angle θ1Successively extract azimuth φ in place1,φ2,…,φmSeismic channel, obtain G1(θ1;
φ1,φ2,…,φm);Then in next incidence angle θ2Successively extract azimuth φ in place1,φ2,…,φmSeismic channel, obtain
G2(θ2;φ1,φ2,…,φm);And so on, finally in maximum incident angle θnSuccessively extract azimuth φ in place1,φ2,…,φm
Seismic channel, obtain Gn(θn;φ1,φ2,…,φm).By the seismic channel of all extractions according to incidence angle it is small to big sequence into
Row arrangement is to get field trace gather (G has been arrived1,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.It is not on space plane between reference azimuth angular-trace gather Zhong Dao and road
Adjacent, similitude is poor, and in addition the interpolation method of the prior art requires signal to have the continuity on time or spatial position, because
This is not suitable for carrying out on reference azimuth angle gathers.And after reference azimuth angle gathers are converted to neighborhood trace gather, in neighborhood trace gather
On space plane position be between seismic channel it is adjacent, the similitude between the road Dao Yu is compared with reference azimuth angle gathers Zhong Dao and road
Between similitude be obviously improved, the enterprising row interpolation of neighborhood trace gather be utilized seismic channel space plane distribution 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 the prior art.
3) in neighborhood trace gather (G1,G2,…,Gn) enterprising row interpolation.
It in one example, can be using MINIMUM WEIGHT norm interpolation method 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 suitable 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 interpolated pretreated orientation angle gathers (C1,
C2,…,Cm)。
In one example, the neighborhood trace gather after interpolation is converted to orientation angle gathers may include: in azimuth φ1,
φ2,…,φmPlace successively extracts the seismic channel of different incidence angles according to sequence from small to large, 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 sequence is arranged, and interpolated pretreated orientation angle gathers (C is obtained1,C2,…,Cm)。
That is, first in minimum azimuth φ1Successively extract incidence angle θ in place1,θ2,…,θnSeismic channel, obtain C1(φ1;
θ1,θ2,…,θn);Then in next azimuth φ2Successively extract incidence angle θ in place1,θ2,…,θnSeismic channel, obtain C2
(φ2;θ1,θ2,…,θn);And so on, finally in maximum azimuth φmSuccessively extract incidence angle θ in place1,θ2,…,θnEarthquake
Road obtains Cm(φm;θ1,θ2,…,θn).The seismic channel of all extractions is arranged according to azimuth is small to big sequence,
Neighborhood trace gather (C is obtained1,C2,…,Cm)。
Embodiment 2
According to another embodiment of the present invention, a kind of trace gather pretreatment unit for azimuthal AVO inverting is provided.The dress
Set may 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
Carry out the unit of interpolation;For the neighborhood trace gather after interpolation to be converted to orientation angle gathers, interpolated pretreated orientation is obtained
Angle gathers (C1,C2,…,Cm) unit.Wherein, m indicates azimuth number, and n indicates incidence angle number.
In one example, A1,A2,…,AmIt can 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
To include: in incidence angle θ1,θ2,…,θnPlace is obtained according to the small seismic channel for arriving big sequence greatly and successively extracting different orientations
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 sequence 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 used by enterprising row interpolation can be with
Including MINIMUM WEIGHT norm interpolation method, based on the interpolation method of Fourier transformation.
In one example, the neighborhood trace gather after interpolation is converted to orientation angle gathers may include: in azimuth φ1,
φ2,…,φmPlace successively extracts the seismic channel of different incidence angles according to sequence from small to large, 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 sequence is arranged, and interpolated pretreated orientation angle gathers (C is obtained1,C2,…,Cm)。
Using example
A concrete application example is given below in the scheme and its effect of the embodiment of the present invention for ease of understanding.This field
It should be understood to the one skilled in the art that the example is only for the purposes of understanding the present invention, any detail is not intended to be limited in any way
The system present invention.
Effect of the invention is illustrated 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 °) each of the corresponding one group of angle gathers being made of different incidence angles in azimuth, incidence angle it 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
There are airways under certain incidence angles.
Fig. 3 is that neighborhood trace gather extracts schematic diagram.As shown in figure 3, being opened from minimum incidence angle on reference azimuth angle gathers
Begin, all azimuthal seismic channels corresponding to the incidence angle is successively extracted according to sequence from small to large, in next incidence
All azimuthal seismic channels corresponding to it are successively acquired also according to sequence 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 successively extracted, thus is adjacent between seismic channel.
Fig. 4 is exemplary neighborhood trace gather, is taken out for the reference azimuth angle gathers of Fig. 2 according to extraction mode shown in Fig. 3
The result taken.From fig. 4, it can be seen that existing airway is more concentrated compared with the airway of reference azimuth trace gather in the neighborhood trace gather.
Fig. 5 a and Fig. 5 b are comparison before and after neighborhood trace gather interpolation.Wherein, Fig. 5 a is the neighborhood trace gather before interpolation, wherein having time
Road;Fig. 5 b is the neighborhood trace gather after interpolation.
Fig. 6 a and Fig. 6 b are orientation angle gathers pretreatment front and back comparison.Wherein, Fig. 6 a is the original-party parallactic angle before pretreatment
Trace gather, there are many airways under some of them incidence angle;Fig. 6 b is to pass through pretreated orientation angle gathers.As it can be seen that of the invention
For the trace gather preprocess method and device of AVO inverting, solve the problems, such as in reference azimuth trace gather there are airway, lack road, be
Azimuthal AVO inverting provides the input trace gather of high quality.
The presently disclosed embodiments is described above, above description is exemplary, and non-exclusive, and
It is not limited to disclosed each embodiment.Without departing from the scope and spirit of illustrated each embodiment, for this skill
Many modifications and changes are obvious for the those of ordinary skill in art field.The selection of term used herein, purport
In the principle, practical application or improvement to the technology in market for best explaining each embodiment, or make the art
Other those of ordinary skill can understand each embodiment disclosed herein.
Claims (8)
1. a kind of trace gather preprocess method for azimuthal AVO inverting, which is characterized 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 interpolated pretreated orientation angle gathers (C1,C2,…,
Cm),
Wherein, m indicates azimuth number, and n indicates incidence angle number;
Wherein, by the orientation angle gathers (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 sequence greatly and successively extracting different orientations1
(θ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 sequence of incidence angle from small to large, obtains the neighborhood trace gather (G1,
G2,…,Gn)。
2. the trace gather preprocess method according to claim 1 for azimuthal AVO inverting, which is characterized 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, which is characterized in that in the neighbour
Domain trace gather (G1,G2,…,Gn) interpolation method used by enterprising row interpolation included MINIMUM WEIGHT norm interpolation method, become based on Fourier
The interpolation method changed.
4. the trace gather preprocess method according to claim 1 for azimuthal AVO inverting, which is characterized in that after interpolation
Neighborhood trace gather be converted to orientation angle gathers and include:
In azimuth φ1,φ2,…,φmPlace successively extracts the seismic channel of different incidence angles according to sequence from small to large, 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 sequence of azimuth from small to large, obtains interpolated pretreated orientation
Angle gathers (C1,C2,…,Cm)。
5. a kind of trace gather pretreatment unit for azimuthal AVO inverting, which is characterized 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 to orientation angle gathers, interpolated pretreated orientation angle gathers (C is obtained1,
C2,…,Cm) unit,
Wherein, m indicates azimuth number, and n indicates incidence angle number;
Wherein, by the orientation angle gathers (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 sequence greatly and successively extracting different orientations1
(θ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 sequence of incidence angle from small to large, obtains the neighborhood trace gather (G1,
G2,…,Gn)。
6. the trace gather pretreatment unit according to claim 5 for azimuthal AVO inverting, which is characterized 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)。
7. the trace gather pretreatment unit according to claim 5 for azimuthal AVO inverting, which is characterized in that in the neighbour
Domain trace gather (G1,G2,…,Gn) interpolation method used by enterprising row interpolation included MINIMUM WEIGHT norm interpolation method, become based on Fourier
The interpolation method changed.
8. the trace gather pretreatment unit according to claim 5 for azimuthal AVO inverting, which is characterized in that after interpolation
Neighborhood trace gather be converted to orientation angle gathers, obtaining interpolated pretreated orientation angle gathers includes:
In azimuth φ1,φ2,…,φmPlace successively extracts the seismic channel of different incidence angles according to sequence from small to large, 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 sequence of azimuth from small to large, obtains interpolated pretreated orientation
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 CN107817517A (en) | 2018-03-20 |
CN107817517B true 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) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113917533B (en) * | 2020-07-10 | 2023-04-28 | 中国石油化工股份有限公司 | TI medium double-linkage omnibearing imaging systematic realization method |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104297800A (en) * | 2014-10-09 | 2015-01-21 | 中国海洋石油总公司 | Self-phase-control prestack inversion method |
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 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015053876A1 (en) * | 2013-10-08 | 2015-04-16 | Exxonmobil Upstream Research Company | Automatic dip picking from wellbore azimuthal image logs |
-
2016
- 2016-09-12 CN CN201610819769.5A patent/CN107817517B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 |
---|
"宽方位地震勘探技术新进展";刘依谋 等;《石油地球物理勘探》;20140630;第49卷(第3期);第596-610页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107817517A (en) | 2018-03-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102221709B (en) | Velocity analysis and dynamic correction method based on formation parameter information | |
Liu et al. | Making reliable shear‐wave splitting measurements | |
CN105301641B (en) | A kind of method and device of azimuthal anisotropy velocity inversion | |
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 | |
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 | |
Barmin et al. | Epicentral location based on Rayleigh wave Empirical Green's Functions from ambient seismic noise | |
CN103713323A (en) | Omnibearing aeolotropy amplitude-preservation imaging and gather extracting 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 | |
CN103901465A (en) | Design method of holographic three-dimensional seismic prospecting and observing system | |
CN103149592A (en) | Method for separating variable offset vertical seismic profile (VSP) wave fields | |
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 | |
CN107817517B (en) | Trace gather preprocess method and device for azimuthal AVO inverting | |
Bulut et al. | Spatiotemporal earthquake clusters along the North Anatolian fault zone offshore Istanbul | |
CN105319587A (en) | Vibroseis shaping filtering method based on first arrival wavelet extraction | |
Cheng et al. | Time-lapse VSP data processing for monitoring CO 2 injection | |
CN106526664A (en) | Time-space domain geometric transform horizontal stacking seismic imaging method | |
Kumar et al. | Seismic anisotropy of the Indian tectonic plate: Deciphering continental scale mantle deformation patterns | |
CN107229069A (en) | A kind of method that velocity analysis is carried out to common receiver data trace gather |
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 |