CN106680875A - Analysis method and apparatus of wide-azimuth seismic migration speed - Google Patents
Analysis method and apparatus of wide-azimuth seismic migration speed Download PDFInfo
- Publication number
- CN106680875A CN106680875A CN201710006684.XA CN201710006684A CN106680875A CN 106680875 A CN106680875 A CN 106680875A CN 201710006684 A CN201710006684 A CN 201710006684A CN 106680875 A CN106680875 A CN 106680875A
- Authority
- CN
- China
- Prior art keywords
- correction
- result
- dynamic
- velocity
- angle gathers
- 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.)
- Pending
Links
- 238000013508 migration Methods 0.000 title claims abstract description 72
- 230000005012 migration Effects 0.000 title claims abstract description 72
- 238000004458 analytical method Methods 0.000 title claims abstract description 22
- 238000012545 processing Methods 0.000 claims abstract description 43
- 238000012937 correction Methods 0.000 claims description 138
- 238000000034 method Methods 0.000 claims description 62
- 238000006243 chemical reaction Methods 0.000 claims description 31
- 230000010354 integration Effects 0.000 claims description 19
- 238000000605 extraction Methods 0.000 claims description 18
- 238000009499 grossing Methods 0.000 claims description 11
- 230000009977 dual effect Effects 0.000 claims description 10
- 230000011514 reflex Effects 0.000 claims description 7
- 235000013399 edible fruits Nutrition 0.000 claims description 5
- 238000010030 laminating Methods 0.000 claims description 4
- 238000005516 engineering process Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 3
- 230000008030 elimination Effects 0.000 description 3
- 238000003379 elimination reaction Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000000717 retained effect Effects 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000005610 quantum mechanics 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
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- 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/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/51—Migration
- G01V2210/514—Post-stack
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/52—Move-out correction
-
- 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/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
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 invention provides an analysis method and apparatus of wide-azimuth seismic migration speed, and relates to the technical field of geological prospecting. The analysis method of wide-azimuth seismic migration speed includes the steps: obtaining trace gather data after performing offset processing on an offset vector sheet; classifying the trace gather data after performing offset processing on the offset vector sheet according to an azimuth angle, and obtaining the azimuth angle trace gather data of the preset number through classification; performing residual speed scanning processing on each azimuth trace gather data, and obtaining the speed scanning result; and performing dual-path integral stacking processing on the speed scanning result so as to extract the speed parameters of each azimuth trace gather data and obtain the speed model of each azimuth trace gather data. Therefore, the analysis method and apparatus of wide-azimuth seismic migration speed relieves the technical problem that in the prior art, the efficiency of eliminating azimuthal anisotropism is relatively lower.
Description
Technical field
The present invention relates to the technical field of geological prospecting, more particularly, to a kind of analysis side of wide-azimuth seismic migration velocity
Method and device.
Background technology
With the fast development of geological exploration field, and the complexity in the face of geological prospecting target and to seismic prospecting essence
The raising that degree is required, has become main flow and the direction of geological exploration techniques development using the geological exploration techniques of wide-azimuth.Width side
The purpose of position earthquake data acquisition is to obtain observed bearing, geophone offset and degree of covering distribution three-dimensional data as uniform as possible
Body.
But, due to the anisotropic property of underground medium, for example, anisotropism and the presence in crack is aligned,
Cause seismic wave different in ground lower edge different azimuth spread speed, it is right that this speed difference causes to be difficult to using single image taking speed
The data in all orientation accurately playback, and will generally carry out high-precision point of azimuthal velocity analysis, and using the optimal of each orientation
Velocity shifts are imaged.
The method that azimuthal anisotropy is eliminated in prior art mainly has two kinds：First is ellipse based on fitting orientation NMO
The remaining orientation dynamic(al) correction method that circle is carried out；Second is HTI anisotropy migration processings.
But, the method for eliminating azimuthal anisotropy at present is less efficient, and cost is generally higher.For example, based on fitting side
Road collection after geophone offset slice of vector (OVT) skew is typically sorted into by the remaining orientation dynamic(al) correction method that position NMO ellipses are carried out
Several orientation angle gathers, as the orientation divided is more, signal to noise ratio is generally relatively low, often will complete each side by manual pickup
To velocity pick, cost is high, and the error and workload due to picking up, it is difficult to take into account the precision in each orientation；HTI is each
Anisotropy skew is often also required to more accurate anisotropic parameters, and the acquisition of the parameter generally needs inverting or speed mould
Type is obtained, and is not very universal method.
The content of the invention
It is an object of the invention to provide the analysis method and device of a kind of wide-azimuth seismic migration velocity, existing to alleviate
Less efficient technical problem in technology to elimination azimuthal anisotropy.
A kind of one side according to embodiments of the present invention, there is provided analysis method of wide-azimuth seismic migration velocity, bag
Include：Obtain the road collection data after geophone offset slice of vector migration processing；To the road collection number after the geophone offset slice of vector migration processing
According to being classified according to azimuth, classification obtains predetermined number orientation angle gathers data；To angle gathers number in orientation each described
According to residual velocity scan process is carried out, velocity scanning result is obtained；The velocity scanning result is carried out at two-way integration superposition
Reason, to extract the speed parameter of each orientation angle gathers data, obtains the speed mould of each orientation angle gathers data
Type.
Further, residual velocity scan process is carried out to angle gathers data in orientation each described, obtains velocity scanning knot
Fruit includes：Inverse dynamic correction is carried out to angle gathers data in orientation each described by the migration velocity of angle gathers data in orientation each described
Just, obtain the first dynamic(al) correction result；Dynamic(al) correction scanning is carried out to the first dynamic(al) correction result, the second dynamic(al) correction result is obtained,
And using the second dynamic(al) correction result as the velocity scanning result.
Further, by the migration velocity of angle gathers data in orientation each described to angle gathers data in orientation each described
Carrying out inverse dynamic correction just includes：By formulaThe reaction correction is carried out to the orientation angle gathers data,
The first dynamic(al) correction result is obtained, wherein, t0It it is zero-offset reflex time, x is offset distance, v is NMO velocity, t1It is
The time difference of the reaction correction；Dynamic(al) correction is carried out to the first dynamic(al) correction result, obtaining the second dynamic(al) correction result includes：To institute
The first dynamic(al) correction result is stated according to formulaDynamic(al) correction scanning is carried out, the second dynamic(al) correction result is obtained,
Wherein, t2It is the time difference of the dynamic(al) correction.
Further, two-way integration overlap-add procedure is carried out to the velocity scanning result, to extract each described azimuth
The speed parameter of road collection data, the rate pattern for obtaining each orientation angle gathers data include：The velocity scanning is tied
Fruit carries out dual path integration overlap-add procedure, obtains stack result；Each described orientation angle gathers is extracted according to the stack result
The speed parameter of data；Interpolation smoothing process is carried out to the speed parameter, the speed of each orientation angle gathers data is obtained
Degree model.
Further, according to the velocity scanning result, dual path product is carried out to angle gathers data in orientation each described
Point overlap-add procedure includes：By formulaTo the velocity scanning result
First time superposition is carried out, the first stack result is obtained, wherein,Represent the velocity scanning
As a result, U (t, h) represents the first dynamic(al) correction result, and h represents offset distance, and v represents stack velocity, τ (t0,x0,h；V) represent base
In the superposition time curve of v；By formulaIt is simultaneously right
The scanning result of the speed carries out second superposition, wherein, w (v) represents weighting function, and the numerical value of the weighting function depends on
In the degree of evening up of road collection, or, the focus level of road collection.
Further, the speed parameter for extracting each orientation angle gathers data according to the stack result includes：Root
According to formulaThe speed parameter is extracted, wherein, α0For the speed parameter.
Other side according to embodiments of the present invention, additionally provides a kind of analysis dress of wide-azimuth seismic migration velocity
Put, including：Acquiring unit, for obtaining the road collection data after geophone offset slice of vector migration processing；Taxon, for described
Road collection data after geophone offset slice of vector migration processing are classified according to azimuth, and classification obtains predetermined number azimuth road
Collection data；Scanning element, for carrying out residual velocity scan process to angle gathers data in orientation each described, obtains velocity scanning
As a result；Extraction unit, for two-way integration overlap-add procedure is carried out to the velocity scanning result, to extract each described azimuth
The speed parameter of road collection data, obtains the rate pattern of each orientation angle gathers data.
Further, the scanning element includes：First correction module, for angle gathers data in orientation each described
Migration velocity carries out reaction correction, obtains the first dynamic(al) correction result；Second correction module, for the first dynamic(al) correction result
Dynamic(al) correction scanning is carried out, the second dynamic(al) correction result is obtained, and using the second dynamic(al) correction result as the velocity scanning result.
Further, first correction module is used for by formulaTo the orientation angle gathers
Data carry out the reaction correction, obtain the first dynamic(al) correction result, wherein, t0It it is zero-offset reflex time, x is skew
It is NMO velocity away from, v, t1It is the time difference of the reaction correction；Second correction module is for first dynamic(al) correction
As a result according to formulaDynamic(al) correction scanning is carried out, the second dynamic(al) correction result is obtained, wherein, t2It is institute
State the time difference of dynamic(al) correction.
Further, the extraction unit includes：Laminating module, for carrying out dual path product to the velocity scanning result
Divide overlap-add procedure, obtain stack result；Extraction module, for extracting each described orientation angle gathers number according to the stack result
According to speed parameter；Smoothing module, for carrying out interpolation smoothing process to the speed parameter, obtains each described orientation
The rate pattern of angle gathers data.
In embodiments of the present invention, road collection data after geophone offset slice of vector migration processing are obtained first, then, to geophone offset
After slice of vector migration processing, data are classified according to azimuth, obtain predetermined number orientation angle gathers data, next, right
Each orientation angle gathers data carries out residual velocity scanning, obtains scanning result, finally, extracts azimuth road according to scanning result
The speed parameter of collection data, obtains rate pattern, to complete the extraction to velocity anisotropy.In embodiments of the present invention, by
Speed is scanned in two-way integration superposition algorithm is employed, the precision of the speed parameter obtained by algorithm scanning and close
Degree is higher, and relative to the anisotropic scheme of removal of prior art, the embodiment of the present invention is offseted in alleviating prior art
Except the less efficient technical problem of azimuthal anisotropy, it is achieved thereby that improving the technology of the efficiency for eliminating azimuthal anisotropy
Effect.
Description of the drawings
In order to be illustrated more clearly that the specific embodiment of the invention or technical scheme of the prior art, below will be to concrete
Needed for embodiment or description of the prior art, accompanying drawing to be used is briefly described, it should be apparent that, in describing below
Accompanying drawing is some embodiments of the present invention, for those of ordinary skill in the art, before creative work is not paid
Put, can be with according to these other accompanying drawings of accompanying drawings acquisition.
Fig. 1 is the flow chart of the analysis method of wide-azimuth seismic migration velocity；
Fig. 2 is a kind of flow chart of residual velocity scan process according to embodiments of the present invention；
Fig. 3 is one kind according to embodiments of the present invention alternatively speed superimposed curves；
Fig. 4 is optional a kind of flow chart of the analysis method of wide-azimuth seismic migration velocity according to embodiments of the present invention；
Fig. 5 is a kind of schematic diagram of the analytical equipment of wide-azimuth seismic migration velocity according to embodiments of the present invention.
Specific embodiment
Technical scheme is clearly and completely described below in conjunction with accompanying drawing, it is clear that described enforcement
Example is a part of embodiment of the invention, rather than the embodiment of whole.Based on the embodiment in the present invention, ordinary skill
The every other embodiment obtained under the premise of creative work is not made by personnel, belongs to the scope of protection of the invention.
In describing the invention, it should be noted that term " " center ", " on ", D score, "left", "right", " vertical ",
The orientation or position relationship of the instruction such as " level ", " interior ", " outward " be based on orientation shown in the drawings or position relationship, merely to
Be easy to description the present invention and simplify description, rather than indicate or imply indication device or element must have specific orientation,
With specific azimuth configuration and operation, therefore it is not considered as limiting the invention.Additionally, term " first ", " second ",
" the 3rd " is only used for describing purpose, and it is not intended that indicating or implying relative importance.
In describing the invention, it should be noted that unless otherwise clearly defined and limited, term " installation ", " phase
Company ", " connection " should be interpreted broadly, for example, it may be being fixedly connected, or being detachably connected, or be integrally connected；Can
Being to be mechanically connected, or electrically connect；Can be joined directly together, it is also possible to be indirectly connected to by intermediary, Ke Yishi
The connection of two element internals.For the ordinary skill in the art, above-mentioned term can be understood at this with concrete condition
Concrete meaning in invention.
According to embodiments of the present invention, there is provided a kind of embodiment of the analysis method of wide-azimuth seismic migration velocity, need
Illustrate, can be in the such as computer system of one group of computer executable instructions the step of the flow process of accompanying drawing is illustrated
Perform, and, although show logical order in flow charts, but in some cases, can be with suitable different from herein
Sequence performs shown or described step.
Fig. 1 is a kind of flow chart of the analysis method of wide-azimuth seismic migration velocity according to embodiments of the present invention, such as Fig. 1
Shown, the method comprises the steps：
Step S102, obtains the road collection data after geophone offset slice of vector migration processing.
Geophone offset represents shot point the distance between to geophone station, in original shot gather data, including the number of multiple geophone offsets
According to, therefore, in embodiments of the present invention, need to extract geophone offset identical shot gather data in original shot gather data, group
Into geophone offset vector sheet data.Then, geophone offset vector sheet data is carried out after migration processing, obtains the skew of geophone offset slice of vector
Road collection data after process, wherein, each azimuthal letter can be retained in the road collection data after geophone offset slice of vector migration processing
Breath.
Road collection data after geophone offset slice of vector migration processing are classified according to azimuth, are classified by step S104
To predetermined number orientation angle gathers data.
In embodiments of the present invention, after the road collection data after geophone offset slice of vector migration processing is got, can be by
The road collection data after geophone offset slice of vector migration processing are classified according to the azimuthal scope for pre-setting, obtain present count
Amount orientation angle gathers data.
For example, azimuth is divided into into following interval：[0,30 °], [30 °, 60 °], [60 °, 90 °], [90 °, 120 °],
[120 °, 150 °], [150 °, 180 °]；Then, by the road collection data after geophone offset slice of vector migration processing according to above-mentioned 6 areas
Between classified, obtain 6 prescription parallactic angle road collection data.
Each orientation angle gathers data is carried out residual velocity scan process, obtains velocity scanning result by step S106.
In embodiments of the present invention, after predetermined number orientation angle gathers data are obtained, can be to every prescription parallactic angle
Road collection data carry out velocity scanning process, obtain velocity scanning result.Specifically, residual velocity scan method will be in following enforcements
Describe in detail in example.
Step S108, carries out two-way integration overlap-add procedure, to extract each orientation angle gathers data to velocity scanning result
Speed parameter, obtain the rate pattern of each orientation angle gathers data.
In embodiments of the present invention, the road collection data after geophone offset slice of vector migration processing are obtained first, then, big gun are examined
Data after slice of vector migration processing are classified according to azimuth, obtain predetermined number orientation angle gathers data, connect down
Come, residual velocity scanning is carried out to each orientation angle gathers data, velocity scanning result is obtained, finally, is carried according to scanning result
The speed parameter of orientation angle gathers data is taken, rate pattern is obtained, to complete the extraction to velocity anisotropy.Of the invention real
Apply in example, speed is scanned as a result of two-way integration superposition algorithm, the speed parameter obtained by the algorithm scanning
Precision and density it is higher, relative to anisotropic scheme is removed in prior art, the embodiment of the present invention alleviates existing
Less efficient technical problem in technology to elimination azimuthal anisotropy, it is achieved thereby that improve eliminating azimuthal anisotropy
The technique effect of efficiency.
Fig. 2 is a kind of flow chart of residual velocity scan process according to embodiments of the present invention, as shown in Fig. 2 to each
Orientation angle gathers data carry out residual velocity scan process, obtain velocity scanning result and comprise the steps：
Step S201, the migration velocity of the road collection data after geophone offset slice of vector migration processing is to each azimuth road
Collection data carry out reaction correction, obtain the first dynamic(al) correction result；
Step S202, carries out dynamic(al) correction scanning to the first dynamic(al) correction result, obtains the second dynamic(al) correction result, and second is moved
Correction result is used as velocity scanning result.
Further, by the migration velocity of angle gathers data in orientation each described to angle gathers data in orientation each described
Carrying out inverse dynamic correction just includes：By formulaAzimuthal road collection data carry out reaction correction, obtain first
Dynamic(al) correction result, wherein, t0It it is zero-offset reflex time, x is offset distance, v is NMO velocity, t1When being that reaction is corrected
Difference；
Dynamic(al) correction scanning is carried out to the first dynamic(al) correction result, obtaining the second dynamic(al) correction result includes：First dynamic(al) correction is tied
Fruit is according to formulaDynamic(al) correction scanning is carried out, the second dynamic(al) correction result is obtained, wherein, t2It is dynamic(al) correction
The time difference.
In embodiments of the present invention, by the migration velocity of orientation angle gathers data to the azimuth road collection data according to anti-
Dynamic(al) correction formulaReaction correction is carried out, the first dynamic(al) correction result is obtained, wherein, t0It is zero-offset reflection
Time, x are offset distances, and v is NMO velocity, and t is reaction timing difference correction result.
After the first dynamic(al) correction result (that is, reaction correction result) is obtained, action can be entered to reaction correction result
Correction rate is scanned.Specifically, by dynamic(al) correction formulaDynamic(al) correction is carried out to reaction correction result.
It should be noted that in embodiments of the present invention, above-mentioned residual velocity scanning process is (that is, reaction is corrected and dynamic
Correction) process carried out on the basis of the migration velocity of the road collection data after geophone offset slice of vector migration processing, for example, presses
Velocity scanning is carried out according to the 90% to 110% of original offset speed, and 1% interval, and in velocity scanning, by moving school
Positive processing method modification orientation road collection data, obtain dynamic(al) correction result (that is, the second dynamic(al) correction result), and dynamic(al) correction result are entered
Row storage.
After by above-mentioned process, the data inhomogeneities after geophone offset slice of vector migration processing and geophone offset are changed
Impact alleviate, become be defined by the datum level for specifying, (common reflection point) CRP road collection data of offset midpoint self excitation and self receiving, connect
Get off, it is possible to which process is overlapped to CRP road collection data.
In embodiments of the present invention, two-way integration overlap-add procedure is carried out to above-mentioned velocity scanning result, to extract each side
The speed parameter of parallactic angle road collection data, the rate pattern for obtaining each orientation angle gathers data comprise the steps：
Step S1081, carries out dual path integration overlap-add procedure, obtains stack result to velocity scanning result；
Step S1082, extracts the speed parameter of each orientation angle gathers data according to stack result；
Step S1083, carries out interpolation smoothing process to speed parameter, obtains the speed mould of each orientation angle gathers data
Type.
Specifically, in embodiments of the present invention, overlap-add procedure algorithm is integrated simultaneously to above-mentioned speed by dual path first
Scanning result is processed, and obtains stack result twice, wherein, dual path integration is superposed to while carrying out to velocity scanning result
It is superimposed twice.After stack result twice is obtained, it is possible to extract each orientation angle gathers data according to stack result twice
Speed parameter, and then, according to speed parameter build rate pattern.In embodiments of the present invention, can be to each azimuth road
The speed parameter of collection data carries out interpolation smoothing process, to obtain the rate pattern of each orientation angle gathers data.
Alternatively, in embodiments of the present invention, it is according to velocity scanning result, double to the carrying out of each orientation angle gathers data
Path integral overlap-add procedure comprises the steps：
Step S1, by formulaVelocity scanning result is carried out
It is superimposed for the first time, obtains the first stack result, wherein,Expression velocity scanning result, U (t,
H) the first dynamic(al) correction result is represented, h represents offset distance, and v represents stack velocity, τ (t0,x0,h；When v) representing superposition based on v away from
Curve；
Specifically, in embodiments of the present invention, above-mentioned formula can be passed through
First time path integral superposition is carried out to velocity scanning result, the first stack result is obtained
Step S2, by formulaSimultaneously to speed
Scanning result carries out second superposition, wherein, w (v) represents weighting function, and what the numerical value of weighting function depended on road collection evens up journey
Degree, or, the focus level of road collection；
Specifically, in embodiments of the present invention, can be according to above-mentioned formula
To the first stack resultSecond path integral superposition is carried out, second stack result is obtained
It should be noted that in the embodiment of the present invention, w (v) typically obeys the probability distribution in quantum mechanics, weighting parameters
By the focus level for evening up degree or road collection of road collection, the attributes such as the slope of road collection, curvature, coherent value ginseng can be can also be
Number.As shown in Figure 3 is a kind of alternatively speed superimposed curves.As shown in figure 3, when speed is low, being superimposed time curve
Unreasonable, now, w (v) is less；When speed is higher, superposition time curve is equally unreasonable, and now, w (v) is equally less.Cause
This, in embodiments of the present invention, w (v) depends primarily on the focus level for evening up degree and road collection of road collection.
Obtaining the first stack resultWith the second stack resultAfterwards, it is possible to according to first
Stack resultWith the second stack resultEach orientation angle gathers data speed parameter is extracted, to build
Divide the rate pattern in orientation.
Alternatively, the speed parameter for extracting each orientation angle gathers data according to stack result includes：According to formulaExtraction rate parameter, wherein, α0For speed parameter.
Fig. 4 is optional a kind of flow chart of the analysis method of wide-azimuth seismic migration velocity according to embodiments of the present invention,
As shown in figure 4, the method comprises the steps：
Step S401, obtains geophone offset vector sheet data；
Step S402, the geophone offset vector sheet data to getting carry out migration processing, and road collection data are obtained after process；
Road collection data after geophone offset slice of vector migration processing are classified according to azimuth, are classified by step S403
To predetermined number orientation angle gathers data；In embodiments of the present invention, angle gathers packet in orientation includes orientation angle gathers data
1, orientation angle gathers data 2 ..., orientation angle gathers data n, specifically, in embodiments of the present invention, to each orientation angle gathers
The processing procedure of data is identical, therefore, illustrate by taking the processing procedure of orientation angle gathers data 1 as an example；
In embodiments of the present invention, after the road collection data after geophone offset slice of vector migration processing is got, can be by
The road collection data after geophone offset slice of vector migration processing are classified according to the azimuthal scope for pre-setting, obtain present count
Amount orientation angle gathers data.
For example, azimuth is divided into into following interval：[0,30 °], [30 °, 60 °], [60 °, 90 °], [90 °, 120 °],
[120 °, 150 °], [150 °, 180 °].Then, by the road collection data after geophone offset slice of vector migration processing according to above-mentioned 6 areas
Between classified, obtain 6 prescription parallactic angle road collection data.
Step S404, azimuthal road collection data 1 carry out reaction correction, obtain reaction correction result；
In embodiments of the present invention, by formulaReaction correction is carried out to migration velocity, is obtained
One dynamic(al) correction result, wherein, t0It it is zero-offset reflex time, x is offset distance, v is NMO velocity, t1It is reaction correction
The time difference.
Step S405, carries out residual NMO correction velocity scanning to reaction correction result, obtains dynamic(al) correction result；
After the first dynamic(al) correction result (that is, reaction correction result) is obtained, action can be entered to reaction correction result
Correction rate is scanned.Specifically, by dynamic(al) correction formulaDynamic(al) correction is carried out to reaction correction result.
Step S406, carries out path integral superposition calculation twice to dynamic(al) correction result, obtains stack result；
Specifically, by formulaVelocity scanning result is carried out
It is superimposed for the first time, obtains the first stack result, wherein,Expression velocity scanning result, U (t,
H) the first dynamic(al) correction result is represented, h represents offset distance, and v represents stack velocity, τ (t0,x0,h；When v) representing superposition based on v away from
Curve；By formulaIs carried out to the first stack result
Secondary superposition, wherein, w (v) represents weighting function, and the numerical value of weighting function depends on the degree of evening up of road collection, or, road collection
Focus level.
Step S407, according to stack result extraction rate parameter；
In embodiments of the present invention, can be according to formulaExtraction rate parameter, wherein, α0For speed parameter.
Step S408, builds the rate pattern in point orientation according to speed parameter；
In embodiments of the present invention, interpolation smoothing process can be carried out to the speed parameter of each orientation angle gathers data,
To obtain the rate pattern of each orientation angle gathers data.
The embodiment of the present invention additionally provides a kind of analytical equipment of wide-azimuth seismic migration velocity, the wide-azimuth seismic migration
The analytical equipment of speed is mainly used in dividing for the wide-azimuth seismic migration velocity provided by execution embodiment of the present invention the above
Analysis method, does concrete introduction below to the analytical equipment of wide-azimuth seismic migration velocity provided in an embodiment of the present invention.
Fig. 5 is a kind of schematic diagram of the analytical equipment of wide-azimuth seismic migration velocity according to embodiments of the present invention, such as Fig. 5
Shown, the analytical equipment of the wide-azimuth seismic migration velocity mainly includes：Acquiring unit 51, taxon 52, scanning element 53
With superposition extraction unit 54, wherein,
Acquiring unit 51, for obtaining the road collection data after geophone offset slice of vector migration processing；
Geophone offset represents shot point the distance between to geophone station, in original shot gather data, including the number of multiple geophone offsets
According to, therefore, in embodiments of the present invention, need to extract geophone offset identical shot gather data in original shot gather data, group
Into geophone offset vector sheet data.Then, geophone offset vector sheet data is carried out after migration processing, obtains the skew of geophone offset slice of vector
Road collection data after process, wherein, each azimuthal letter can be retained in the road collection data after geophone offset slice of vector migration processing
Breath.
Taxon 52, for classifying according to azimuth to the road collection data after geophone offset slice of vector migration processing,
Classification obtains predetermined number orientation angle gathers data；
In embodiments of the present invention, after the road collection data after geophone offset slice of vector migration processing is got, can be by
The road collection data after geophone offset slice of vector migration processing are classified according to the azimuthal scope for pre-setting, obtain present count
Amount orientation angle gathers data.
For example, azimuth is divided into into following interval：[0,30 °], [30 °, 60 °], [60 °, 90 °], [90 °, 120 °],
[120 °, 150 °], [150 °, 180 °]；Then, by the road collection data after geophone offset slice of vector migration processing according to above-mentioned 6 areas
Between classified, obtain 6 prescription parallactic angle road collection data.
Scanning element 53, for carrying out residual velocity scan process to each orientation angle gathers data, obtains velocity scanning
As a result；
In embodiments of the present invention, after predetermined number orientation angle gathers data are obtained, can be to every prescription parallactic angle
Road collection data carry out velocity scanning process, obtain velocity scanning result.Specifically, residual velocity scan method will be in following enforcements
Describe in detail in example.
Superposition extraction unit 54, for two-way integration overlap-add procedure is carried out to velocity scanning result, to extract each orientation
The speed parameter of angle gathers data, obtains the rate pattern of each orientation angle gathers data.
In embodiments of the present invention, the road collection data after geophone offset slice of vector migration processing are obtained first, then, big gun are examined
Data after slice of vector migration processing are classified according to azimuth, obtain predetermined number orientation angle gathers data, connect down
Come, residual velocity scanning is carried out to each orientation angle gathers data, velocity scanning result is obtained, finally, is carried according to scanning result
The speed parameter of orientation angle gathers data is taken, rate pattern is obtained, to complete the extraction to velocity anisotropy.Of the invention real
Apply in example, speed is scanned as a result of two-way integration superposition algorithm, the speed parameter obtained by the algorithm scanning
Precision and density it is higher, relative to anisotropic scheme is removed in prior art, the embodiment of the present invention alleviates existing
Less efficient technical problem in technology to elimination azimuthal anisotropy, it is achieved thereby that improve eliminating azimuthal anisotropy
The technique effect of efficiency.
Alternatively, scanning element includes：First correction module, for entering to the migration velocity of each orientation angle gathers data
Row reaction is corrected, and obtains the first dynamic(al) correction result；Second correction module, sweeps for carrying out dynamic(al) correction to the first dynamic(al) correction result
Retouch, obtain the second dynamic(al) correction result, and using the second dynamic(al) correction result as velocity scanning result.
Alternatively, the first correction module is used for by formulaAzimuthal road collection data carry out reaction
Correction, obtains the first dynamic(al) correction result, wherein, t0It it is zero-offset reflex time, x is offset distance, v is NMO velocity, t1It is
The time difference of reaction correction；Second correction module for the first dynamic(al) correction result according to formulaEnter action
Correct scan, obtains the second dynamic(al) correction result, wherein, t2It is the time difference of dynamic(al) correction.
Alternatively, being superimposed extraction unit includes：Laminating module, for dual path integration superposition is carried out to velocity scanning result
Process, obtain stack result；Extraction module, for the speed parameter of each orientation angle gathers data is extracted according to stack result；
Smoothing module, for carrying out interpolation smoothing process to speed parameter, obtains the rate pattern of each orientation angle gathers data.
Alternatively, laminating module includes：First superposition submodule, for by formula
First time superposition is carried out to the velocity scanning result simultaneously, the first stack result is obtained, wherein,Velocity scanning result is represented, U (t, h) represents the first dynamic(al) correction result, and h represents offset distance,
V represents stack velocity, τ (t0,x0,h；V) represent the superposition time curve based on v；First superposition submodule, for by formulaSecond superposition is carried out to the first stack result, its
In, w (v) represents weighting function, and the numerical value of weighting function depends on the degree of evening up of road collection, or, the focus level of road collection.
Alternatively, extraction module includes：Extracting sub-module, for according to formulaExtraction rate parameter, wherein, α0
For speed parameter.
Finally it should be noted that：Various embodiments above only to illustrate technical scheme, rather than a limitation；To the greatest extent
Pipe has been described in detail to the present invention with reference to foregoing embodiments, it will be understood by those within the art that：Its according to
So the technical scheme described in foregoing embodiments can be modified, or which part or all technical characteristic are entered
Row equivalent；And these modifications or replacement, do not make the essence of appropriate technical solution depart from various embodiments of the present invention technology
The scope of scheme.
Claims (10)
1. a kind of analysis method of wide-azimuth seismic migration velocity, it is characterised in that include：
Obtain the road collection data after geophone offset slice of vector migration processing；
Road collection data after the geophone offset slice of vector migration processing are classified according to azimuth, classification obtains predetermined number
Individual orientation angle gathers data；
Residual velocity scan process is carried out to angle gathers data in orientation each described, velocity scanning result is obtained；
Two-way integration overlap-add procedure is carried out to the velocity scanning result, to extract the speed of each orientation angle gathers data
Parameter, obtains the rate pattern of each orientation angle gathers data.
2. analysis method according to claim 1, it is characterised in that residue is carried out to angle gathers data in orientation each described
Velocity scanning is processed, and obtaining velocity scanning result includes：
Reaction correction is carried out to angle gathers data in orientation each described by the migration velocity of angle gathers data in orientation each described,
Obtain the first dynamic(al) correction result；
Dynamic(al) correction scanning is carried out to the first dynamic(al) correction result, the second dynamic(al) correction result is obtained, and by second dynamic(al) correction
As a result as the velocity scanning result.
3. analysis method according to claim 2, it is characterised in that
Reaction correction is carried out to angle gathers data in orientation each described by the migration velocity of angle gathers data in orientation each described
Including：By formulaThe reaction correction is carried out to the orientation angle gathers data, described first is obtained
Dynamic(al) correction result, wherein, t0It it is zero-offset reflex time, x is offset distance, v is NMO velocity, t1It is the reaction correction
The time difference；
Dynamic(al) correction is carried out to the first dynamic(al) correction result, obtaining the second dynamic(al) correction result includes：First dynamic(al) correction is tied
Fruit is according to formulaDynamic(al) correction scanning is carried out, the second dynamic(al) correction result is obtained, wherein, t2It is described dynamic
The time difference of correction.
4. analysis method according to claim 2, it is characterised in that two-way integration is carried out to the velocity scanning result folded
Plus process, to extract the speed parameter of each orientation angle gathers data, obtain the speed of each orientation angle gathers data
Degree model includes：
Dual path integration overlap-add procedure is carried out to the velocity scanning result, stack result is obtained；
The speed parameter of each orientation angle gathers data is extracted according to the stack result；
Interpolation smoothing process is carried out to the speed parameter, the rate pattern of each orientation angle gathers data is obtained.
5. analysis method according to claim 4, it is characterised in that dual path integration is carried out to the velocity scanning result
Overlap-add procedure, obtaining stack result includes：
By formulaThe velocity scanning result is carried out for the first time
Superposition, obtains the first stack result, wherein,Represent the velocity scanning result, U (t, h)
The first dynamic(al) correction result is represented, h represents offset distance, and v represents stack velocity, τ (t0,x0,h；When v) representing the superposition based on v
Away from curve；
By formulaSimultaneously the velocity scanning is tied
Fruit carries out second superposition, wherein, w (v) represents weighting function, and what the numerical value of the weighting function depended on road collection evens up journey
Degree, or, the focus level of road collection.
6. analysis method according to claim 5, it is characterised in that each described orientation is extracted according to the stack result
The speed parameter of angle gathers data includes：
According to formulaThe speed parameter is extracted, wherein, α0For the speed parameter.
7. a kind of analytical equipment of wide-azimuth seismic migration velocity, it is characterised in that include：
Acquiring unit, for obtaining the road collection data after geophone offset slice of vector migration processing；
Taxon, for classifying according to azimuth to the road collection data after the geophone offset slice of vector migration processing, point
Class obtains predetermined number orientation angle gathers data；
Scanning element, for carrying out residual velocity scan process to angle gathers data in orientation each described, obtains velocity scanning knot
Really；
Superposition extraction unit, for two-way integration overlap-add procedure is carried out to the velocity scanning result, to extract each described side
The speed parameter of parallactic angle road collection data, obtains the rate pattern of each orientation angle gathers data.
8. analytical equipment according to claim 7, it is characterised in that the scanning element includes：
First correction module, for carrying out reaction correction to the migration velocity of angle gathers data in orientation each described, obtains first
Dynamic(al) correction result；
Second correction module, for carrying out dynamic(al) correction scanning to the first dynamic(al) correction result, obtains the second dynamic(al) correction result, and
Using the second dynamic(al) correction result as the velocity scanning result.
9. analytical equipment according to claim 8, it is characterised in that
First correction module is used for by formulaThe orientation angle gathers data are carried out described anti-
Dynamic(al) correction, obtains the first dynamic(al) correction result, wherein, t0It it is zero-offset reflex time, x is offset distance, and v is dynamic(al) correction speed
Degree, t1It is the time difference of the reaction correction；
Second correction module for the first dynamic(al) correction result according to formulaCarry out dynamic(al) correction
Scanning, obtains the second dynamic(al) correction result, wherein, t2It is the time difference of the dynamic(al) correction.
10. analytical equipment according to claim 8, it is characterised in that the extraction unit includes：
Laminating module, for dual path integration overlap-add procedure is carried out to the velocity scanning result, obtains stack result；
Extraction module, for the speed parameter of each orientation angle gathers data is extracted according to the stack result；
Smoothing module, for carrying out interpolation smoothing process to the speed parameter, obtains each described orientation angle gathers number
According to rate pattern.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710006684.XA CN106680875A (en) | 2017-01-05 | 2017-01-05 | Analysis method and apparatus of wide-azimuth seismic migration speed |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710006684.XA CN106680875A (en) | 2017-01-05 | 2017-01-05 | Analysis method and apparatus of wide-azimuth seismic migration speed |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106680875A true CN106680875A (en) | 2017-05-17 |
Family
ID=58849953
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710006684.XA Pending CN106680875A (en) | 2017-01-05 | 2017-01-05 | Analysis method and apparatus of wide-azimuth seismic migration speed |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106680875A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110737019A (en) * | 2018-07-19 | 2020-01-31 | 中国石油化工股份有限公司 | Processing method and system for improving velocity spectrum precision |
| CN112666606A (en) * | 2019-10-15 | 2021-04-16 | 中国石油天然气集团有限公司 | Method and device for determining anisotropic speed of near-surface azimuth |
| CN113049202A (en) * | 2021-03-08 | 2021-06-29 | 中国地震局工程力学研究所 | Local weighted regression correction method and system for acceleration integral displacement |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104166161A (en) * | 2014-08-19 | 2014-11-26 | 成都理工大学 | Method and device for predicating fractures based on elliptical velocity inversion of anisotropism |
| CN105093299A (en) * | 2015-07-24 | 2015-11-25 | 中国石油天然气集团公司 | Observation system optimization method based on offset vector tile technology and apparatus thereof |
| CN105242308A (en) * | 2015-10-14 | 2016-01-13 | 中国石油天然气集团公司 | Superposition method and apparatus of wide-azimuth earthquake data |
| CN105301641A (en) * | 2015-10-30 | 2016-02-03 | 中国石油天然气集团公司 | Azimuthal anisotropy speed inversion method and apparatus |
| CN106094029A (en) * | 2016-08-24 | 2016-11-09 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | The method utilizing offset distance vector sheet geological data Predicating Reservoir Fractures |
| CN106199717A (en) * | 2016-07-08 | 2016-12-07 | 中国石油天然气集团公司 | Azimuthal anisotropy Residual moveout correction method and device |
-
2017
- 2017-01-05 CN CN201710006684.XA patent/CN106680875A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104166161A (en) * | 2014-08-19 | 2014-11-26 | 成都理工大学 | Method and device for predicating fractures based on elliptical velocity inversion of anisotropism |
| CN105093299A (en) * | 2015-07-24 | 2015-11-25 | 中国石油天然气集团公司 | Observation system optimization method based on offset vector tile technology and apparatus thereof |
| CN105242308A (en) * | 2015-10-14 | 2016-01-13 | 中国石油天然气集团公司 | Superposition method and apparatus of wide-azimuth earthquake data |
| CN105301641A (en) * | 2015-10-30 | 2016-02-03 | 中国石油天然气集团公司 | Azimuthal anisotropy speed inversion method and apparatus |
| CN106199717A (en) * | 2016-07-08 | 2016-12-07 | 中国石油天然气集团公司 | Azimuthal anisotropy Residual moveout correction method and device |
| CN106094029A (en) * | 2016-08-24 | 2016-11-09 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | The method utilizing offset distance vector sheet geological data Predicating Reservoir Fractures |
Non-Patent Citations (3)
| Title |
|---|
| E.LANDA, ETC: "Path-integral seismic imaging", 《GEOPHYSICAL PROSPECTING》 * |
| 张保庆等: "宽方位地震资料处理技术及应用效果", 《石油地球物理勘探》 * |
| 田梦等: "宽方位角地震勘探与常规地震勘探对比研究", 《大庆石油地质与开发》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110737019A (en) * | 2018-07-19 | 2020-01-31 | 中国石油化工股份有限公司 | Processing method and system for improving velocity spectrum precision |
| CN110737019B (en) * | 2018-07-19 | 2021-06-29 | 中国石油化工股份有限公司 | Processing method and system for improving velocity spectrum precision |
| CN112666606A (en) * | 2019-10-15 | 2021-04-16 | 中国石油天然气集团有限公司 | Method and device for determining anisotropic speed of near-surface azimuth |
| CN113049202A (en) * | 2021-03-08 | 2021-06-29 | 中国地震局工程力学研究所 | Local weighted regression correction method and system for acceleration integral displacement |
| CN113049202B (en) * | 2021-03-08 | 2022-07-12 | 中国地震局工程力学研究所 | Local weighted regression correction method and system for acceleration integral displacement |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104142514B (en) | Three-dimensional earthquake observing system quantitative designing method | |
| CN103927731B (en) | Low-altitude remote sensing image rapid and automatic splicing method without POS assisting | |
| CN102914791B (en) | Kirchhoff prestack time migration method for processing seismic data of undulating surface | |
| CN102221709B (en) | Velocity analysis and dynamic correction method based on formation parameter information | |
| CN104932010B (en) | A kind of diffracted wave separation method based on the sparse Radon transformation of shortcut fringing | |
| CN101598803B (en) | Method for directly obtaining stacked section of converted wave | |
| CN105954797B (en) | The fault recognition method and apparatus of seismic data | |
| CN104570110B (en) | A kind of multi component signal based on the matching of ripple in length and breadth combines velocity analysis method | |
| EP2920615A1 (en) | Process for creating image gathers | |
| CN102401908B (en) | Method for suppressing multiple waves by the aid of different-mode weighting sparse parabola Radon transform | |
| CN107561589A (en) | A kind of near surface S-wave interval velocity method for establishing model | |
| CN105259570A (en) | Seaborne time shifting earthquake receiving point displacement correction method | |
| CN104199103B (en) | Static correction method and static correction device for converted wave | |
| CN102901985A (en) | Depth domain layer speed correcting method suitable for undulating surface | |
| CN105866833A (en) | VSP-CDP stacking method and three-dimensional VSP coverage number calculation method | |
| CN105931216A (en) | Image processing technology-based airport pavement distress detection method | |
| CN101609167B (en) | Cross-hole seismic wave equation pre-stack depth migration imaging method based on relief surface | |
| CN105510975B (en) | Improve the method and device of geological data signal to noise ratio | |
| CN106680875A (en) | Analysis method and apparatus of wide-azimuth seismic migration speed | |
| CN104360384B (en) | Micro-seismic event localization method based on wave energy automatically scanning in length and breadth and device | |
| CN107861156A (en) | The extracting method and device of diffracted wave | |
| CN105301648A (en) | Method of acquiring common reflection surface stacking parameters | |
| CN109884693A (en) | Adaptively move towards normal-moveout spectrum acquiring method and system | |
| CN107942391A (en) | A kind of seabed geophone station localization method and device | |
| CN106338766A (en) | Pre-stack time migration method based on split-step Fourier algorithm |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170517 |