CN105954798B - A kind of method and apparatus of the pre-stack time migration speed of determining relief surface - Google Patents
A kind of method and apparatus of the pre-stack time migration speed of determining relief surface Download PDFInfo
- Publication number
- CN105954798B CN105954798B CN201610265242.2A CN201610265242A CN105954798B CN 105954798 B CN105954798 B CN 105954798B CN 201610265242 A CN201610265242 A CN 201610265242A CN 105954798 B CN105954798 B CN 105954798B
- Authority
- CN
- China
- Prior art keywords
- time migration
- stack time
- stack
- velocity analysis
- result
- 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
- 230000005012 migration Effects 0.000 title claims abstract description 213
- 238000013508 migration Methods 0.000 title claims abstract description 213
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000004458 analytical method Methods 0.000 claims abstract description 121
- 238000012545 processing Methods 0.000 claims abstract description 19
- 238000004364 calculation method Methods 0.000 claims abstract description 16
- 238000003384 imaging method Methods 0.000 claims description 22
- 238000012937 correction Methods 0.000 claims description 20
- 238000001228 spectrum Methods 0.000 claims description 17
- 230000003247 decreasing effect Effects 0.000 claims description 4
- 239000011148 porous material Substances 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 10
- 230000005284 excitation Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 235000013399 edible fruits Nutrition 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 230000007423 decrease Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000011017 operating method Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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/30—Analysis
-
- 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
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 embodiment of the present application provides a kind of method and apparatus of the pre-stack time migration speed of determining relief surface.This method comprises: obtaining the round trip time difference of the common point plane of reference at multiple pre-stack time migration velocity analysis position to datum level；Respectively at velocity analysis position, one group of scanning speed function changed with seismic wave reflection interval is formed；To each seismic channel in multiple input seismic channels, at velocity analysis position, using scanning speed function and the corresponding round trip time difference, calculation formula is corrected when travelling by relief surface and carries out pre-stack time migration calculating, obtains the pre-stack time migration result of each seismic channel in multiple input seismic channels；Pre-stack time migration result is merged into processing, the pre-stack time migration result after merging is changed into pre-stack time migration velocity analysis result.The pre-stack time migration velocity analysis result of the complicated earth surfaces such as relief surface can be more accurately determined using technical solution provided by the embodiments of the present application.
Description
Technical field
The present invention relates to geophysical exploration technical field of data processing more particularly to a kind of prestacks of determining relief surface
The method and apparatus of time migration speed.
Background technique
In geophysical exploration technical field of data processing, pre-stack time migration is common imaging method, is had become
Conventional treatment means need corresponding pre-stack time migration velocity field in the application process of pre-stack time migration imaging method
As the parameter that must be used.Therefore, pre-stack time migration speed is accurately determined out, corresponding pre-stack time migration speed is established
It is very crucial to spend field, is the deciding factor that pre-stack time migration can popularize.
In the prior art, the application of Patent No. 200910080622.9 proposes a kind of determining pre-stack time migration most
The method of good velocity profile, pre-stack time migration normal-moveout spectrum are one of the methods of true pre-stack time migration speed, and its advantages are
Application method is simple, essentially identical with the usage mode of stacking velocity analysis.But the technical solution of the patent is for water
Level land table, for the exploration work area of non-horizontal earth's surface, existing way is the method using static correction seismic exploration data
It is corrected to datum level from the common point plane of reference, as shown in Fig. 1, is allowed to meet the hypothesis of horizontal earth's surface.S is excitation in Fig. 1
Point, R are receiving point, and I is imaging point, and V indicates velocity analysis position.This correction is that handle is excited in S point in the received earthquake of R point
Data postpone the corresponding time, put received seismic data in R ' as in S ' point excitation, doing so can make subsequent processing generate one
Fixed error and is considered because path is S '-S-I-R-R ' when the travelling of seismic data after correction in subsequent processing
S '-I-R ', the two is not identical, and difference increases as correcting value increases.As exploration work area is to complicated earth surface region
Extend, the survey area very big for some surface relief amplitudes, the mode of existing static correction to datum level can make subsequent place
Reason generates unacceptable error, is also in this way, influencing relief surface pre-stack time migration on pre-stack time migration velocity analysis
Application effect.
Therefore, the method that the prior art needs a kind of pre-stack time migration speed of determining relief surface, can be more acurrate
The pre-stack time migration velocity analysis of the complicated earth surfaces such as relief surface is determined as a result, expanding pre-stack time migration speed point in ground
The application range of analysis.
Summary of the invention
The purpose of the application is to provide a kind of method and apparatus of the pre-stack time migration speed of determining relief surface, can be with
More accurately determine the pre-stack time migration velocity analysis of the complicated earth surfaces such as relief surface as a result, expanding pre-stack time migration
The application range of velocity analysis.
The method and apparatus of the pre-stack time migration speed of determining relief surface provided by the present application are achieved in that
A kind of method of the pre-stack time migration speed of determining relief surface, which comprises
Obtain the round trip of the common point plane of reference at multiple pre-stack time migration velocity analysis position to datum level
The time difference；
Respectively at the multiple pre-stack time migration velocity analysis position, changed with given with seismic wave reflection interval
Speed centered on speed, using the given speed interval with the variation of seismic wave reflection interval as step-length, with the center
Based on speed, 1 is increased or decreased to N number of step-length, forms one group of scanning speed function changed with seismic wave reflection interval,
In, N is preset parameter；
To each seismic channel in multiple input seismic channels, at the multiple pre-stack time migration velocity analysis position,
Corrected using the scanning speed function and the corresponding round trip time difference, when being travelled by relief surface calculation formula into
Row pre-stack time migration calculates, and obtains the pre-stack time migration result of each seismic channel in the multiple input seismic channel；
The pre-stack time migration result of each seismic channel is merged into processing, by the pre-stack time migration after merging
As a result it is changed into pre-stack time migration velocity analysis result.
In a preferred embodiment, correction calculation formula includes following formula when the relief surface is travelled:
In above formula, t indicates the hourage on pre stack data road, tiIndicate the time of imaging point, x1Indicate shot point to current
Pre-stack time migration velocity analysis position horizontal distance, x2Indicate geophone station to current pre-stack time migration velocity analysis
The horizontal distance of position, v (ti) indicate imaging point scanning speed, Δ t indicate input seismic channel the common point plane of reference extremely
The round trip time difference of datum level and present speed analyze the TWT of the common point plane of reference at position to datum level
Difference between difference, i indicate imaging point.
In a preferred embodiment, the value of N is 5 to 100.
In a preferred embodiment, the pre-stack time migration velocity analysis result includes:
Migration velocity spectra, CRP gather and migrated section section.
In a preferred embodiment, the pre-stack time migration result by after merging is changed into pre-stack time migration
Velocity analysis result includes:
The form that the pre-stack time migration result that different scanning speed function is generated is changed into Migration velocity spectra is defeated
Out, the pre-stack time migration velocity analysis result of the Migration velocity spectra form includes that preset quantity scanning speed function is corresponding
CRP gather and migrated section section.
A kind of device of the pre-stack time migration speed of determining relief surface, described device include:
The time difference obtains module, for obtaining the common point plane of reference at multiple pre-stack time migration velocity analysis position extremely
The round trip time difference of datum level；
Scanning speed function determination module is used for respectively at the multiple pre-stack time migration velocity analysis position, with
Speed centered on the given speed with the variation of seismic wave reflection interval, with the given speed changed with seismic wave reflection interval
Degree interval is used as step-length, based on the central speed, increases or decreases 1 to N number of step-length, forms one group with seismic reflection
The scanning speed function of time change, wherein N is preset parameter；
Data computation module, for it is multiple input seismic channels in each seismic channels, it is inclined in the multiple prestack time
It moves at velocity analysis position, using the scanning speed function and the corresponding round trip time difference, passes through relief surface trip
Calculation formula is corrected when row and carries out pre-stack time migration calculating, obtains the prestack of each seismic channel in the multiple input seismic channel
Time migration result；
Data processing module will be closed for the pre-stack time migration result of each seismic channel to be merged processing
Pre-stack time migration result after and is changed into pre-stack time migration velocity analysis result.
In a preferred embodiment, correction calculation formula includes following formula when the relief surface is travelled:
In above formula, t indicates the hourage on pre stack data road, tiIndicate the time of imaging point, x1Indicate shot point to current
Pre-stack time migration velocity analysis position horizontal distance, x2Indicate geophone station to current pre-stack time migration velocity analysis
The horizontal distance of position, v (ti) indicate imaging point scanning speed, Δ t indicate input seismic channel the common point plane of reference extremely
The round trip time difference of datum level and present speed analyze the TWT of the common point plane of reference at position to datum level
Difference between difference, i indicate imaging point.
In a preferred embodiment, the value of N is 5 to 100.
In a preferred embodiment, the pre-stack time migration velocity analysis result includes:
Migration velocity spectra, CRP gather and migrated section section.
In a preferred embodiment, the pre-stack time migration result by after merging is changed into pre-stack time migration
Velocity analysis result includes:
The form that the pre-stack time migration result that different scanning speed function is generated is changed into Migration velocity spectra is defeated
Out, the pre-stack time migration velocity analysis result of the Migration velocity spectra form includes that preset quantity scanning speed function is corresponding
CRP gather and migrated section section.
The application in the multiple pre-stack time migration speed by dividing to each seismic channel in multiple input seismic channels
It analyses at position, using the multiple pre-stack time migration speed of the acquisition of the scanning speed function changed with seismic wave reflection interval and acquisition
The common point plane of reference at degree analysis position was corrected when being travelled by relief surface and is calculated to the round trip time difference of datum level
Formula carries out pre-stack time migration calculating, obtains the pre-stack time migration knot of each seismic channel in the multiple input seismic channel
Fruit；Finally, the pre-stack time migration result of each seismic channel is merged processing, by the pre-stack time migration after merging
As a result it is changed into pre-stack time migration velocity analysis result.It compared with prior art, can using technical solution provided by the present application
More accurately to determine the pre-stack time migration velocity analysis of the complicated earth surfaces such as relief surface as a result, to expand prestack time inclined
Move the application range of velocity analysis.
Detailed description of the invention
In order to illustrate the technical solutions in the embodiments of the present application or in the prior art more clearly, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
The some embodiments recorded in application, for those of ordinary skill in the art, in the premise of not making the creative labor property
Under, it is also possible to obtain other drawings based on these drawings.
Fig. 1 is the schematic diagram for seismic data being corrected to from the common point plane of reference in the prior art datum level；
Fig. 2 is a kind of stream of embodiment of the method for the pre-stack time migration speed of determining relief surface provided by the present application
Cheng Tu；
Fig. 3 is a kind of school of embodiment of the method for the pre-stack time migration speed of determining relief surface provided by the present application
Positive schematic diagram；
Fig. 4 is that the pre-stack time migration velocity analysis result obtained using the embodiment of the present application and existing prestack time are inclined
Move a kind of embodiment schematic diagram that velocity analysis result is compared；
Fig. 5 is in a kind of embodiment of the device of the pre-stack time migration speed of determining relief surface provided by the present application
Structural schematic diagram.
Specific embodiment
In order to make those skilled in the art better understand the technical solutions in the application, below in conjunction with the application reality
The attached drawing in example is applied, the technical scheme in the embodiment of the application is clearly and completely described, it is clear that described implementation
Example is only a part of the embodiment of the application, instead of all the embodiments.Based on the embodiment in the application, this field is general
The application guarantor all should belong in logical technical staff every other embodiment obtained without making creative work
The range of shield.
The specific implementation of the embodiment of the present application is described in detail with several specific examples below.
Introduce a kind of a kind of implementation of the method for the pre-stack time migration speed of determining relief surface of the application first below
Example.Fig. 2 is a kind of flow chart of embodiment of the method for the pre-stack time migration speed of determining relief surface provided by the present application,
This application provides the method operating procedures as described in embodiment or flow chart, but based on conventional or can without creative labor
To include more or less operating procedure.The step of enumerating in embodiment sequence is only in numerous step execution sequences
A kind of mode does not represent and unique executes sequence.It, can be according to embodiment when system or client production in practice executes
Either method sequence shown in the drawings executes or executes (such as environment of parallel processor or multiple threads) parallel.Tool
Body as shown in Fig. 2, the method may include:
S110: the round trip of the common point plane of reference at multiple pre-stack time migration velocity analysis position to datum level is obtained
Travel-time difference.
The common point plane of reference in the embodiment of the present application, at available multiple pre-stack time migration velocity analysis position
To the round trip time difference of datum level.
Specifically, the quantity of pre-stack time migration velocity analysis position may include the practical exploration previously according to target area
Situation setting.The common point plane of reference to datum level the round trip time difference can from storage survey data database,
Seismic data trace header or corresponding T-V obtain function.
S120: respectively at the multiple pre-stack time migration velocity analysis position, when giving with seismic reflection
Between speed centered on the speed that changes, using the given speed interval with the variation of seismic wave reflection interval as step-length, with institute
Based on stating central speed, 1 is increased or decreased to N number of step-length, forms one group of scanning speed changed with seismic wave reflection interval
Function, wherein N is preset parameter.
Specifically, in the embodiment of the present application, with it is given with the variation of seismic wave reflection interval speed (such as: at some
Position is 500ms 1800m/s, 2000ms 2650m/s, 3000ms 2900m/s, 5000ms 3800m/s etc.) centered on speed
Degree, then again with change with the seismic reflection time speed interval (such as: 500ms 15m/s, 2000ms 18m/s,
3000ms 21m/s, 5000ms 25m/s etc.) it is used as step-length.After central speed and step-length has been determined, so that it may in this
Based on heart speed, one, two, three ... speed intervals are subtracted, or add one, two, three ... speed intervals,
To form one group of scanning speed function changed with seismic wave reflection interval.For example, the value of above-mentioned N can be 5 to 100
In any one positive integer, it is preferred that any one positive integer in 10 to 50 can be can be with the value of N, specifically taken certainly
Value can in the light of actual conditions be chosen with demand, and the application is it is not limited here.
S130: to each seismic channel in multiple input seismic channels, in the multiple pre-stack time migration velocity analysis position
Place is set, using the scanning speed function and the corresponding round trip time difference, corrects and calculates when being travelled by relief surface
Formula carries out pre-stack time migration calculating, obtains the pre-stack time migration knot of each seismic channel in the multiple input seismic channel
Fruit.
In the embodiment of the present application, step S120 obtain at the velocity analysis position corresponding scanning speed function it
Afterwards, each seismic channel in multiple input seismic channels can be adopted at the multiple pre-stack time migration velocity analysis position
With the scanning speed function and the corresponding round trip time difference, calculation formula is corrected when travelling by relief surface and is carried out
Pre-stack time migration calculates, and obtains the pre-stack time migration result of each seismic channel in the multiple input seismic channel.
Specifically, correction calculation formula may include following formula when the relief surface is travelled:
In above formula, t indicates the hourage on pre stack data road, tiIndicate the time of imaging point, x1Indicate shot point to current
Pre-stack time migration velocity analysis position horizontal distance, x2Indicate geophone station to current pre-stack time migration velocity analysis
The horizontal distance of position, v (ti) indicate imaging point scanning speed, Δ t indicate input seismic channel the common point plane of reference extremely
The round trip time difference of datum level and present speed analyze the TWT of the common point plane of reference at position to datum level
Difference between difference, i indicate imaging point.
In above-mentioned the embodiment of the present application, in order to avoid seismic data is corrected to datum level to rear from the common point plane of reference
Continuous processing generates too big error, at the same meet pre-stack time migration algorithm requirement horizontal earth's surface it is assumed that the application in prestack
In the pore diameter range of time migration, at the common point plane of reference seismic data corrections to velocity analysis position, such as Fig. 3 institute
Show, Fig. 3 is that a kind of correction of embodiment of the method for the pre-stack time migration speed of determining relief surface provided by the present application is shown
It is intended to.In figure, S is excitation point, and R is receiving point, and I is imaging point, and V indicates velocity analysis position, and S ' is the excitation point after correction,
R ' is the receiving point after correction.Because the aperture of pre-stack time migration is less than exploration work area, the correcting value done so is less than
Correcting value from common point plane of reference the to datum level, to reduce the error of subsequent processing caused by static correction.This side
Formula make each velocity analysis position pre-stack time migration velocity analysis calculate meet horizontal earth's surface it is assumed that although each speed
The horizontal ground level of degree analysis position may be different.Therefore, the pre-stack time migration that can be obtained by the way of the application
Velocity analysis result is better than seismic data being corrected to from the common point plane of reference pre-stack time migration speed calculated after datum level
Degree analysis result.
The embodiment of the present application is by combining the corresponding common point plane of reference extremely in pre-stack time migration velocity analysis position
Round trip time difference of datum level carries out corresponding correction process, by the calculating of pre-stack time migration velocity analysis according to
The common point plane of reference to the time difference of datum level and present speed for inputting seismic channel analyzes the common point plane of reference of position extremely
The time difference of datum level is corrected, and so that pre-stack time migration velocity analysis result is can be used in relief surface, when expanding prestack
Between migration velocity analysis application range.
S140: merging processing for the pre-stack time migration result of each seismic channel, when by prestack after merging
Between migration result be changed into pre-stack time migration velocity analysis result.
It, can will be described after step S130 obtains the pre-stack time migration result data in the embodiment of the present application
The pre-stack time migration result of each seismic channel merges processing, and the pre-stack time migration result after merging is changed into prestack
Time migration velocity analysis result.Specifically, the pre-stack time migration velocity analysis result can include at least in following
It is a kind of:
Migration velocity spectra, CRP gather and migrated section section.
Specifically, the pre-stack time migration result by after merging be changed into pre-stack time migration velocity analysis result can
To include:
The form that the pre-stack time migration result that different scanning speed function is generated is changed into Migration velocity spectra is defeated
Out, the pre-stack time migration velocity analysis result of the Migration velocity spectra form includes that preset quantity scanning speed function is corresponding
CRP gather and migrated section section.
In a specific embodiment, the complicated earth surfaces such as relief surface are chosen corresponding at a certain velocity analysis position
Pre-stack time migration velocity analysis result be compared with existing pre-stack time migration velocity analysis result.Fig. 4 is to utilize
The embodiment of the present application obtain pre-stack time migration velocity analysis result and existing pre-stack time migration velocity analysis result into
A kind of embodiment schematic diagram that row compares.As shown in figure 4, the pre-stack time migration result of 410 expression relief surfaces carries out in figure
(result of the renewal speed analysis can obtain better pre-stack time migration knot to the schematic diagram of the result of renewal speed analysis
Fruit, in general, the result that the renewal speed can be analyzed is as pre-stack time migration velocity analysis result reference standard),
420 indicate the schematic diagram of the pre-stack time migration velocity analysis result obtained using the embodiment of the present application, and 430 indicate the prior art
The middle schematic diagram the pre-stack time migration velocity analysis result after seismic data corrections to datum level, 440 indicate the prior art
Pre-stack time migration speed after moving on to the common point plane of reference when the middle pre-stack time migration velocity analysis result by datum level
The schematic diagram of degree analysis result.It can be seen from the figure that the pre-stack time migration velocity analysis result and renewal speed of the application
Analysis result essentially coincides, and only small difference illustrates that the pre-stack time migration velocity analysis result of the application is substantially accurate.
On the contrary, the result difference of other two kinds of pre-stack time migration velocity analysis results in the prior art and renewal speed analysis
Larger, the result difference of pre-stack time migration velocity analysis result and renewal speed analysis on datum level is more than 30%, even if
When result difference also with renewal speed analysis of the result that moves on to after the common point plane of reference be more than 15%.Therefore, this Shen is utilized
Please the pre-stack time migration velocity analysis result accuracy rate that obtains of the technical solution that passes through of embodiment it is higher.
By a kind of embodiment of the method for the pre-stack time migration speed of determining relief surface of above the application as it can be seen that this Shen
Please by being adopted at the multiple pre-stack time migration velocity analysis position to each seismic channel in multiple input seismic channels
At multiple pre-stack time migration velocity analysis position with the scanning speed function and acquisition that change with seismic wave reflection interval
The common point plane of reference to datum level the round trip time difference, when being travelled by relief surface correct calculation formula carry out prestack when
Between calculations of offset, obtain it is the multiple input seismic channel in each seismic channel pre-stack time migration result；Finally, by described each
The pre-stack time migration result of a seismic channel merges processing, when the pre-stack time migration result after merging is changed into prestack
Between migration velocity analysis result.Compared with prior art, it can more accurately be determined using technical solution provided by the present application
The pre-stack time migration velocity analysis of the complicated earth surfaces such as relief surface is as a result, expand the application of pre-stack time migration velocity analysis
Range.
Based on the same inventive concept, a kind of pre-stack time migration of determining relief surface is additionally provided in the embodiment of the present invention
The device of speed, as described in the following examples.Due to determining that the device solution of the pre-stack time migration speed of relief surface is asked
The principle of topic is similar to the method for pre-stack time migration speed of relief surface is determined, it is thus determined that the prestack time of relief surface
The implementation of the device of migration velocity may refer to determine the implementation of the method for the pre-stack time migration speed of relief surface, repetition
Place repeats no more.Used below, term " unit " or " module " can be the software and/or hardware for realizing predetermined function
Combination.Although device described in following embodiment is preferably realized with software, hardware or software and hardware
The realization of combination is also that may and be contemplated.Fig. 5 is the pre-stack time migration speed of determining relief surface provided by the present application
Device a kind of embodiment in structural schematic diagram；As shown in figure 5, described device 500 can also include:
The time difference obtains module 510, can be used for obtaining the common point at multiple pre-stack time migration velocity analysis position
The plane of reference to datum level the round trip time difference；
Scanning speed function determination module 520 can be used for respectively in the multiple pre-stack time migration velocity analysis position
Place is set, the speed centered on the given speed with the variation of seismic wave reflection interval, with given with seismic wave reflection interval
The speed interval of variation is as step-length, based on the central speed, increases or decreases 1 to N number of step-length, forms one group everywhere
The scanning speed function of seismic wave reflection interval variation, wherein N is preset parameter；
Data computation module 530 can be used for each seismic channel in multiple input seismic channels, in the multiple prestack
At time migration velocity analysis position, using the scanning speed function and the corresponding round trip time difference, pass through fluctuating
Correction calculation formula carries out pre-stack time migration calculating when earth's surface is travelled, and obtains each seismic channel in the multiple input seismic channel
Pre-stack time migration result；
Data processing module 540 can be used for the pre-stack time migration result of each seismic channel merging place
Reason, is changed into pre-stack time migration velocity analysis result for the pre-stack time migration result after merging.
In a preferred embodiment, correction calculation formula may include following formula when the relief surface is travelled:
In above formula, t indicates the hourage on pre stack data road, tiIndicate the time of imaging point, x1Indicate shot point to current
Pre-stack time migration velocity analysis position horizontal distance, x2Indicate geophone station to current pre-stack time migration velocity analysis
The horizontal distance of position, v (ti) indicate imaging point scanning speed, Δ t indicate input seismic channel the common point plane of reference extremely
The round trip time difference of datum level and present speed analyze the TWT of the common point plane of reference at position to datum level
Difference between difference, i indicate imaging point.
In a preferred embodiment, the value of N can be 5 to 100.
In a preferred embodiment, the pre-stack time migration velocity analysis result may include:
Migration velocity spectra, CRP gather and migrated section section.
In a preferred embodiment, the pre-stack time migration result by after merging is changed into pre-stack time migration
Velocity analysis result may include:
The form that the pre-stack time migration result that different scanning speed function is generated is changed into Migration velocity spectra is defeated
Out, the pre-stack time migration velocity analysis result of the Migration velocity spectra form includes that preset quantity scanning speed function is corresponding
CRP gather and migrated section section.
It can by a kind of embodiment of the method and apparatus of the pre-stack time migration speed of determining relief surface of above the application
See, the application by it is multiple input seismic channels in each seismic channels, in the multiple pre-stack time migration velocity analysis position
Place is set, using the multiple pre-stack time migration speed of the acquisition of the scanning speed function changed with seismic wave reflection interval and acquisition point
The common point plane of reference at position is analysed to the round trip time difference of datum level, calculation formula is corrected when travelling by relief surface
Pre-stack time migration calculating is carried out, the pre-stack time migration result of each seismic channel in the multiple input seismic channel is obtained；Most
Afterwards, the pre-stack time migration result of each seismic channel is merged into processing, by the pre-stack time migration result after merging
It is changed into pre-stack time migration velocity analysis result.It compared with prior art, can be more using technical solution provided by the present application
The pre-stack time migration velocity analysis of the complicated earth surfaces such as relief surface is accurately determined out as a result, expanding pre-stack time migration speed
Spend the application range of analysis.
Obviously, those skilled in the art should be understood that each module of the above-mentioned embodiment of the present invention or each step can be with
It is realized with general computing device, they can be concentrated on a single computing device, or be distributed in multiple computing devices
On composed network, optionally, they can be realized with the program code that computing device can perform, it is thus possible to by it
Store and be performed by computing device in the storage device, and in some cases, can be held with the sequence for being different from herein
The shown or described step of row, perhaps they are fabricated to each integrated circuit modules or will be multiple in them
Module or step are fabricated to single integrated circuit module to realize.In this way, the embodiment of the present invention be not limited to it is any specific hard
Part and software combine.
The foregoing is only a preferred embodiment of the present invention, is not intended to restrict the invention, for the skill of this field
For art personnel, the embodiment of the present invention can have various modifications and variations.All within the spirits and principles of the present invention, made
Any modification, equivalent substitution, improvement and etc. should all be included in the protection scope of the present invention.
Claims (2)
1. a kind of method of the pre-stack time migration speed of determining relief surface, which is characterized in that the described method includes:
In the pore diameter range of pre-stack time migration, the common point plane of reference of seismic data corrections to velocity analysis position
Place；
Obtain the round trip time difference of the common point plane of reference at multiple pre-stack time migration velocity analysis position to datum level；
Respectively at the multiple pre-stack time migration velocity analysis position, with the given speed changed with seismic wave reflection interval
Speed centered on degree, using the given speed interval with the variation of seismic wave reflection interval as step-length, with the central speed
Based on, 1 is increased or decreased to N number of step-length, forms one group of scanning speed function changed with seismic wave reflection interval, wherein N
For preset parameter, the value of N is any one positive integer in 10 to 50；
Each seismic channel in multiple input seismic channels is used at the multiple pre-stack time migration velocity analysis position
The scanning speed function and the corresponding round trip time difference correct calculation formula and are folded when travelling by relief surface
Preceding time migration calculates, and obtains the pre-stack time migration result of each seismic channel in the multiple input seismic channel；
The pre-stack time migration result of each seismic channel is merged into processing, by the pre-stack time migration result after merging
It is changed into pre-stack time migration velocity analysis result；
Wherein, correction calculation formula includes following formula when the relief surface is travelled:
In above formula, t indicates the hourage on pre stack data road, tiIndicate the time of imaging point, x1Indicate that shot point to current is folded
The horizontal distance of preceding time migration velocity analysis position, x2Indicate geophone station to current pre-stack time migration velocity analysis position
Horizontal distance, v (ti) indicating the scanning speed of imaging point, Δ t indicates the common point plane of reference of input seismic channel to benchmark
The common point plane of reference at round trip time difference in face and present speed analysis position to datum level the round trip time difference it
Between difference, i indicate imaging point；
Wherein, the pre-stack time migration result by after merging is changed into pre-stack time migration velocity analysis result and includes:
The pre-stack time migration result that different scanning speed function is generated is changed into the form output of Migration velocity spectra, institute
The pre-stack time migration velocity analysis result for stating Migration velocity spectra form includes that preset quantity scanning speed function is corresponding anti-altogether
Exit point trace gather and migrated section section.
2. a kind of device of the pre-stack time migration speed of determining relief surface, which is characterized in that described device includes:
The time difference obtains module, for obtaining the common point plane of reference at multiple pre-stack time migration velocity analysis position to benchmark
The round trip time difference in face；
Scanning speed function determination module is used for respectively at the multiple pre-stack time migration velocity analysis position, with given
The speed with the variation of seismic wave reflection interval centered on speed, between the given speed with the variation of seismic wave reflection interval
Every as step-length, based on the central speed, 1 is increased or decreased to N number of step-length, forms one group with seismic wave reflection interval
The scanning speed function of variation, wherein N is preset parameter, and the value of N is any one positive integer in 10 to 50；
Data computation module, for each seismic channels in multiple input seismic channels, in the multiple pre-stack time migration speed
At degree analysis position, using the scanning speed function and the corresponding round trip time difference, when being travelled by relief surface
It corrects calculation formula and carries out pre-stack time migration calculating, obtain the prestack time of each seismic channel in the multiple input seismic channel
Migration result；
Data processing module, in the pore diameter range of pre-stack time migration, seismic data corrections to velocity analysis position
The common point plane of reference at, be also used to the pre-stack time migration result of each seismic channel merging processing, will close
Pre-stack time migration result after and is changed into pre-stack time migration velocity analysis result；Wherein, the prestack by after merging
Time migration result is changed into pre-stack time migration velocity analysis result
The pre-stack time migration result that different scanning speed function is generated is changed into the form output of Migration velocity spectra, institute
The pre-stack time migration velocity analysis result for stating Migration velocity spectra form includes that preset quantity scanning speed function is corresponding anti-altogether
Exit point trace gather and migrated section section；
Wherein, correction calculation formula includes following formula when the relief surface is travelled:
In above formula, t indicates the hourage on pre stack data road, tiIndicate the time of imaging point, x1Indicate that shot point to current is folded
The horizontal distance of preceding time migration velocity analysis position, x2Indicate geophone station to current pre-stack time migration velocity analysis position
Horizontal distance, v (ti) indicating the scanning speed of imaging point, Δ t indicates the common point plane of reference of input seismic channel to benchmark
The common point plane of reference at round trip time difference in face and present speed analysis position to datum level the round trip time difference it
Between difference, i indicate imaging point.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610265242.2A CN105954798B (en) | 2016-04-26 | 2016-04-26 | A kind of method and apparatus of the pre-stack time migration speed of determining relief surface |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610265242.2A CN105954798B (en) | 2016-04-26 | 2016-04-26 | A kind of method and apparatus of the pre-stack time migration speed of determining relief surface |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105954798A CN105954798A (en) | 2016-09-21 |
| CN105954798B true CN105954798B (en) | 2019-05-07 |
Family
ID=56916249
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610265242.2A Active CN105954798B (en) | 2016-04-26 | 2016-04-26 | A kind of method and apparatus of the pre-stack time migration speed of determining relief surface |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105954798B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112230285B (en) * | 2019-07-15 | 2023-09-26 | 中国石油天然气集团有限公司 | Correction method and correction device for earthquake data excitation delay time |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101036069A (en) * | 2004-10-04 | 2007-09-12 | Pgs美洲公司 | Amplitude preserving prestack migration method |
| EP2020609A1 (en) * | 2006-09-04 | 2009-02-04 | Geosystem S.r.l. | Methods and apparatus for geophysical exploration via joint inversion |
| CN102914791A (en) * | 2011-08-05 | 2013-02-06 | 中国石油天然气集团公司 | Kirchhoff prestack time migration method for processing seismic data of undulating surface |
| CN105137480A (en) * | 2015-08-20 | 2015-12-09 | 中国石油集团东方地球物理勘探有限责任公司 | Method and apparatus for determining pre-stack time migration velocity of VTI medium |
-
2016
- 2016-04-26 CN CN201610265242.2A patent/CN105954798B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101036069A (en) * | 2004-10-04 | 2007-09-12 | Pgs美洲公司 | Amplitude preserving prestack migration method |
| EP2020609A1 (en) * | 2006-09-04 | 2009-02-04 | Geosystem S.r.l. | Methods and apparatus for geophysical exploration via joint inversion |
| CN102914791A (en) * | 2011-08-05 | 2013-02-06 | 中国石油天然气集团公司 | Kirchhoff prestack time migration method for processing seismic data of undulating surface |
| CN105137480A (en) * | 2015-08-20 | 2015-12-09 | 中国石油集团东方地球物理勘探有限责任公司 | Method and apparatus for determining pre-stack time migration velocity of VTI medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105954798A (en) | 2016-09-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| RU2694621C1 (en) | Method and apparatus for processing seismic data | |
| CN107765298B (en) | A kind of method and device of determining near-surface velocity model | |
| CN105425289B (en) | The method and apparatus for determining low frequency wave impedance | |
| CN104422963B (en) | Variable-offset VSP data time difference correction method | |
| Cameron et al. | Seismic velocity estimation from time migration | |
| CN103760603B (en) | The prestack time migration method of converted wave geological data and device | |
| CN104199103B (en) | Static correction method and static correction device for converted wave | |
| CN103399300A (en) | Wave packet superposition microseism ground location method | |
| CN107817525B (en) | A kind of near surface constant speed fill method and device | |
| CN107229071B (en) | A kind of subsurface structure inversion imaging method | |
| CN104570076A (en) | Automatic seismic wave first-arrival picking method based on dichotomy | |
| CN114994754B (en) | Seismic source mechanism joint inversion method based on direct wave and depth seismic phase initial motion polarity | |
| CN106054252B (en) | A kind of method and device of pre-stack time migration | |
| CN106842302B (en) | Method and device for editing first arrivals in batches | |
| CN105954798B (en) | A kind of method and apparatus of the pre-stack time migration speed of determining relief surface | |
| CN104914469B (en) | Static correcting method and device in a kind of converted shear wave | |
| CN108957553B (en) | Method and device for dynamic correction without stretching distortion by recursion correction of dynamic correction value | |
| CN105445793B (en) | A kind of method and device for determining bad track data | |
| CN106707345B (en) | Method and device for identifying lithology of angle elastic parameter | |
| CN105353406A (en) | Method and device for generating angle gather | |
| CN109581494B (en) | Pre-stack migration method and device | |
| CN105425290B (en) | A kind of method and device of pre-stack time migration | |
| CN107255832B (en) | A kind of inversion method of subsurface structure | |
| CN105137484B (en) | The method and apparatus for estimating the center of geological data trace gather energy value | |
| CN111208558A (en) | Method and device for establishing ultra-deep low-amplitude three-dimensional geological structure |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |