CN105607117B - Time-domain offset imaging method and device - Google Patents
Time-domain offset imaging method and device Download PDFInfo
- Publication number
- CN105607117B CN105607117B CN201510577384.8A CN201510577384A CN105607117B CN 105607117 B CN105607117 B CN 105607117B CN 201510577384 A CN201510577384 A CN 201510577384A CN 105607117 B CN105607117 B CN 105607117B
- Authority
- CN
- China
- Prior art keywords
- mrow
- imaging
- imaging point
- coarse grid
- point
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention provides a kind of time-domain offset imaging method and device, wherein, this method includes：Read one of seismic data；By double flat root formula, ask for the focus in the road to the first of imaging road coarse grid imaging point walk when, and coarse grid imaging point to the second of receiving point walk when, and as the two way travel time of coarse grid imaging point the sum of when described first is walked and when described second walks；By cubic convolution interpolation formula, according to the two way travel time of the coarse grid imaging point, interpolation calculation goes out the two way travel time of refined net imaging point；According to the two way travel time of the refined net imaging point, the interrogation energy on input seismic channel, as the imaging results of the aqueduct, final migration result is the superposition of all input seismic channel imaging results.The present invention solves the technical problem that the high-frequency information in migration stack section in the prior art is seriously damaged, and has reached the technique effect for effectively improving imaging results precision and accuracy.
Description
Technical field
The present invention relates to seismic exploration technique field, more particularly to a kind of time-domain offset imaging method and device.
Background technology
Due to having higher computational efficiency and flexibility concurrently, Kirchhoff offsets are still at present the industrial quarters general time
Domain offset imaging method.The essence of Kirchhoff offsets is superposition of the seismic data along subsurface imaging point diffraction travel time curve
Summation, i.e. calculate each imaging point in underground to the focus corresponding to the road and receiving point firstly for each Seismic Traces
When walking, interrogation energy and be overlapped on this seismic channel when then being walked according to this, final stack result be the point into
Picture value.
However, during the realization of Kirchhoff offsets, in order to improve computational efficiency, usually in sparse imaging grid
When upper (coarse grid imaging point) calculates and amplitude, and coarse grid imaging point is interpolated into by refined net by linear interpolation when walking
On point, the calculating of picture value is then carried out into again.
However, since linear interpolation is a kind of interpolation method of low order, it can not ensure Kirchhoff migration operators institute
The continuity of the single order needed, so as to cause the high-frequency energy that abnormal flat bed shape is produced at interpolation knot.In migration stack
Afterwards, the relevant high-frequency energy meeting coherent enhancement in this same interpolation knot position, causes the high frequency in final migration stack section
Information is seriously damaged, the strong influence precision of follow-up inversion procedure and attributes extraction.
In view of the above-mentioned problems, not yet propose effective solution at present.
The content of the invention
An embodiment of the present invention provides a kind of time-domain offset imaging method, to solve migration stack section in the prior art
In high-frequency information be seriously damaged and have impact on the technical problem of the precision of follow-up inversion procedure and attributes extraction, this method bag
Include：
Read the seismic data of a seismic channel；
By double flat root formula, ask for the focus of the seismic channel to the first of imaging road coarse grid imaging point walk when, and
Coarse grid imaging point to the second of receiving point walk when, and the sum of when when described first is walked, is walked with described second as coarse grid into
The two way travel time of picture point；
By cubic convolution interpolation formula, according to the two way travel time of the coarse grid imaging point, interpolation calculation goes out refined net
The two way travel time of imaging point；
According to the two way travel time of the refined net imaging point, the interrogation energy on the seismic channel, and by the energy of pickup into
Row superposition, the imaging results using stack result as the seismic channel；
Above-mentioned steps are repeated to all seismic channels, obtain the migration imaging result in an imaging road.
In one embodiment, the seismic data of a seismic channel is read, including：
Obtain the CMP trace gather seismic datas in region to be analyzed；
One of seismic data is read from the CMP trace gathers seismic data.
In one embodiment, the double flat root formula is：
Wherein, t represents the two way travel time of coarse grid imaging point, tsRepresent focus to the first of imaging road coarse grid imaging point
When walking, tgRepresent coarse grid imaging point to the second of receiving point walk when, t0Represent the coarse grid imaging point it is corresponding it is vertical walk when
Depth, vrmsRepresent the root mean sequare velocity at the coarse grid imaging point, x represents that the coarse grid imaging point arrives common point
Horizontal distance, h represent half geophone offset.
In one embodiment, the formula of the cubic convolution interpolation is：
tn=Tk-1(-s3+2s2-s)/2+Tk(3s3-5s2+2)/2
+Tk+1(-3s3+4s2+s)/2+Tk+2(s3-s2)/2
Wherein, tnRepresent the two way travel time of refined net imaging point obtained after interpolation calculation, Tk-1、Tk、Tk+1、Tk+2Represent thin
When walking of four coarse grid imaging points near grid imaging point, s represent refined net step-length.
In one embodiment, the coarse grid imaging point be according to predetermined Grid Sampling interval, to imaging road into
Row sparse sampling obtains.
The embodiment of the present invention additionally provides a kind of time-domain migration imaging device, is cutd open with solving migration stack in the prior art
High-frequency information in face is seriously damaged and have impact on the technical problem of the precision of follow-up inversion procedure and attributes extraction, the device
Including：
Read module, for reading the seismic data of a seismic channel；
First two way travel time computing module, for by double flat root formula, asking for the focus of the seismic channel to imaging road
The first of coarse grid imaging point when walking, and coarse grid imaging point to the second of receiving point walk when, and when described first is walked and institute
State the second two way travel time the sum of when walking as coarse grid imaging point；
Second two way travel time computing module, for by cubic convolution interpolation formula, according to the coarse grid imaging point
Two way travel time, interpolation calculation go out the two way travel time of refined net imaging point；
Image-forming module, for the two way travel time according to the refined net imaging point, the interrogation energy on the seismic channel, and will
The energy of pickup is overlapped, the imaging results using stack result as the seismic channel.
In one embodiment, the read module includes：
Acquiring unit, for obtaining the CMP trace gather seismic datas in region to be analyzed；
Reading unit, for reading one of seismic data from the CMP trace gathers seismic data.
In one embodiment, the double flat root formula is：
Wherein, t represents the two way travel time of coarse grid imaging point, tsRepresent focus to the first of imaging road coarse grid imaging point
When walking, tgRepresent coarse grid imaging point to the second of receiving point walk when, t0Represent the coarse grid imaging point it is corresponding it is vertical walk when
Depth, vrmsRepresent the root mean sequare velocity at the coarse grid imaging point, x represents that the coarse grid imaging point arrives common point
Horizontal distance, h represent half geophone offset.
In one embodiment, the formula of the cubic convolution interpolation is：
tn=Tk-1(-s3+2s2-s)/2+Tk(3s3-5s2+2)/2
+Tk+1(-3s3+4s2+s)/2+Tk+2(s3-s2)/2
Wherein, tnRepresent the two way travel time of refined net imaging point obtained after interpolation calculation, Tk-1、Tk、Tk+1、Tk+2Represent thin
When walking of four coarse grid imaging points near grid imaging point, s represent refined net step-length.
In one embodiment, the coarse grid imaging point be according to predetermined Grid Sampling interval, to imaging road into
Row sparse sampling obtains.
In embodiments of the present invention, the two way travel time of coarse grid imaging point, Ran Houtong have been tried to achieve by double flat root formula
The two way travel time that the two way travel time of coarse grid imaging point is changed into refined net imaging point by cubic convolution interpolation is crossed, to meet to deviate
The continuity of operator single order, so that the high-frequency information solved in migration stack section in the prior art is seriously damaged and influences
The problem of subsequent treatment and Explanation Accuracy, the technique effect for effectively improving imaging results precision and accuracy is reached.
Brief description of the drawings
Attached drawing described herein is used for providing a further understanding of the present invention, forms the part of the application, not
Form limitation of the invention.In the accompanying drawings：
Fig. 1 is the method flow diagram of time-domain offset imaging method according to embodiments of the present invention；
Fig. 2 is the principle schematic of cubic convolution interpolation according to embodiments of the present invention；
Fig. 3 is according to the single track migrated section schematic diagram deviated by traditional Kirchhoff；
Fig. 4 is the obtained single track migrated section schematic diagram of the method provided according to embodiments of the present invention；
Fig. 5 be according to by traditional Kirchhoff deviate using the single track migrated section schematic diagram after high-pass filtering；
Fig. 6 is that the obtained single track migrated section using after high-pass filtering of the method provided according to embodiments of the present invention shows
It is intended to；
Fig. 7 is according to the multiple tracks migrated section schematic diagram deviated by traditional Kirchhoff；
Fig. 8 is the obtained multiple tracks migrated section schematic diagram of the method provided according to embodiments of the present invention；
Fig. 9 be according to by traditional Kirchhoff deviate using the multiple tracks migrated section schematic diagram after high-pass filtering；
Figure 10 is the obtained multiple tracks migrated section using after high-pass filtering of the method provided according to embodiments of the present invention
Schematic diagram；
Figure 11 is the structure diagram of time-domain migration imaging device according to embodiments of the present invention.
Embodiment
It is right with reference to embodiment and attached drawing for the object, technical solutions and advantages of the present invention are more clearly understood
The present invention is described in further details.Here, the exemplary embodiment and its explanation of the present invention are used to explain the present invention, but simultaneously
It is not as a limitation of the invention.
A kind of time-domain offset imaging method is provided in the present invention, when can eliminate Kirchhoff by this method
Between deviate in high-frequency noise interference information, in the method using cubic convolution interpolation come into walking when interpolation calculation, with
Meet the continuity of migration operator single order, so as to effectively eliminate the high-frequency noise energy of flat bed shape, ensure migration imaging knot
The correctness of fruit medium-high frequency information.
Specifically, as shown in Figure 1, the time-domain offset imaging method may comprise steps of：
Step 101：Read the seismic data of a seismic channel；
From earthquake record choose seismic channel when, can with sequence together a genuine selection, can also
It is to randomly select a seismic channel, it is specific to choose mode the application as long as all seismic channels may finally be traveled through
It is not construed as limiting.
Step 102：By double flat root formula, the focus of the seismic channel is asked for the first of imaging road coarse grid imaging point
When walking, and coarse grid imaging point to the second of receiving point walk when, and the sum of when described first is walked and when described second walks as
The two way travel time of coarse grid imaging point；
That is, into during line displacement, the seismic data of a seismic channel can be read in first, for any in migrated section
One imaging road, then can utilize double flat root formula calculate the corresponding focus of input seismic to subsurface imaging point with
And subsurface imaging point is to the two way travel time of receiving point.
Wherein, double flat root formula can be expressed as：
Wherein, t represents the two way travel time of coarse grid imaging point, tsRepresent focus to the first of imaging road coarse grid imaging point
When walking, tgRepresent coarse grid imaging point to the second of receiving point walk when, t0Represent the coarse grid imaging point it is corresponding it is vertical walk when
Depth, vrmsRepresent the root mean sequare velocity at the coarse grid imaging point, x represents the coarse grid imaging point to common point (institute
The common point of meaning is exactly imaging point to the center of receiving point line) horizontal distance, h represents half geophone offset.
This therefore is calculated in the enterprising walking of coarse grid imaging point mainly for improving computational efficiency, wherein, coarse grid into
Picture point can be according to certain Grid Sampling interval, carry out what sparse sampling obtained to imaging road (being imaged seismic channel).
Assuming that there is N number of mesh point in seismic imaging road, corresponding refined net walks sequential and is classified as tn, 0≤n≤N, if according to one
Fixed grid interval Δ k, chooses K mesh point therein, then can obtain corresponding coarse grid imaging point walks time series Tk,
0≤k≤K, and there are following relation Tk=tk×Δk。
Step 103：By cubic convolution interpolation formula, according to the two way travel time of the coarse grid imaging point, interpolation calculation
Go out the two way travel time of refined net imaging point；
The two way travel time for the correspondence coarse grid imaging point tried to achieve can be calculated corresponding using cubic convolution interpolation
The two way travel time of refined net imaging point, the formula of the cubic convolution interpolation can be expressed as：
tn=Tk-1(-s3+2s2-s)/2+Tk(3s3-5s2+2)/2
+Tk+1(-3s3+4s2+s)/2+Tk+2(s3-s2)/2
Wherein, which is a kind of interpolation method of high-order, it is ensured that the company of interpolation result at least single order
Continuous property, as shown in Fig. 2, tnRepresent the two way travel time of refined net imaging point obtained after interpolation calculation, Tk-1、Tk、Tk+1、Tk+2Represent
When walking of four coarse grid imaging points near the refined net imaging point, s represent refined net step-length.
Further, in order to improve the computational efficiency of interpolation, it is corresponding that each coarse grid imaging point can be precomputed
Interpolation coefficient：
Cof1=(- s3+2s2-s)/2
Cof2=(3s3-5s2+2)/2
Cof3=(- 3s3+4s2+s)/2
Cof4=(s3-s2)/2
The formula of so above-mentioned cubic convolution interpolation can be reduced to：
tn=Tk-1cof1+Tkcof2+Tk+1cof3+Tk+2cof4
Step 104：According to the two way travel time of the refined net imaging point, the interrogation energy on the seismic channel, and will pickup
Energy be overlapped, the imaging results using stack result as the seismic channel.
Two way travel time i.e. based on the refined net imaging point being calculated interrogation energy, Ran Houtong in earthquake record is inputted
After crossing amplitude weight, assignment on refined net imaging point so that obtain the road seismic data of a corresponding seismic channel currently into
As the imaging results on road.
Then, 101 are repeated the above steps to all seismic channels to step 104, you can obtain the imaging in an imaging road
As a result.Further, above-mentioned steps are carried out to all imaging roads, final imaging results are added up, you can obtain most
Whole migration imaging result.
In the above-described embodiments, the two way travel time of coarse grid imaging point has been tried to achieve by double flat root formula, has then been passed through
The two way travel time of coarse grid imaging point is changed into the two way travel time of refined net imaging point by cubic convolution interpolation, to meet that offset is calculated
The continuity of sub- single order, so that the high-frequency information solved in migration stack section in the prior art is seriously damaged and have impact on
Processing and the technical problem of Explanation Accuracy, have reached the technique effect for effectively improving imaging results precision and accuracy.
Above-mentioned time-domain offset imaging method is specifically described with reference to Fig. 3 to Figure 10.
Time-domain velocity analysis is carried out to the seismic data in certain exploratory area, obtains mean-square-root velocity field, then, it is genuine to read in one
Data are shaken, ask for the two way travel time of the corresponding focus in the road and receiving point to coarse grid imaging point using double flat root formula, so
The two way travel time of obtained coarse grid imaging point is interpolated into refined net imaging point using above-mentioned cubic convolution interpolation formula afterwards
On.Wherein, using cubic convolution interpolation calculate refined net imaging point walk constantly, it is necessary to 4 coarse grids around the point into
When walking of picture point.
After the two way travel time of refined net imaging point is got, energy can be picked up in earthquake record according to two way travel time
Amount, is imaged on mesh point, it is only necessary to repeat the above steps to all imaging roads, you can obtain one using being put into after amplitude weight
The imaging results of a input seismic channel.
It is inclined with reference to figure 3 and Fig. 4, the respectively single track using traditional quadrature method offset and method provided in an embodiment of the present invention
Move as a result, it can be seen from Fig. 3 and Fig. 4 when retaining whole band energy, both are almost without any difference.With reference to 5 He of figure
Fig. 6, the progress that result and method provided in an embodiment of the present invention after high-pass filtering are respectively carried out using traditional quadrature method are high
Single track migration result after pass filter, passes through the carry out high pass of method provided in an embodiment of the present invention it can be seen from Fig. 5 and Fig. 6
The high-frequency noise of flat bed shape can be effectively eliminated after filtering.
All seismic datas are inputted, are repeated the above steps, and imaging results are added up, you can are obtained final inclined
Move stack result.
With reference to figure 7 and Fig. 8, respectively deviated using the multiple tracks of traditional quadrature method and method provided in an embodiment of the present invention folded
Add section, it can be seen from Fig. 7 and Fig. 8 when without filtering process, both differences can not equally differentiate.With reference to 9 He of figure
Figure 10, the progress that result and method provided in an embodiment of the present invention after high-pass filtering are respectively carried out using traditional quadrature method are high
Multiple tracks migration result after pass filter, it can be seen from Fig. 9 and Figure 10 after high-pass filtering is carried out, due to being put down after migration stack
The high-frequency energy of stratiform obtains coherent enhancement, and the superposition of the obtained multiple tracks migration result of method provided in an embodiment of the present invention is cutd open
Face completely eliminates the high-frequency noise of flat bed shape, and can differentiate the high-frequency energy of some imaging layer positions.
In this example, be overlapped velocity analysis obtain root mean sequare velocity on the basis of, using double flat root formula into
The two way travel time of row coarse grid imaging point calculates, and then when walking and shakes using cubic convolution interpolation calculating refined net imaging point
Width, then carries out the calculating of imaging energy using imaging formula, after above-mentioned processing is carried out to all earthquake records, you can completes
Whole migration process.
Based on same inventive concept, a kind of time-domain migration imaging device is additionally provided in the embodiment of the present invention, it is such as following
Embodiment described in.Since the principle that time-domain migration imaging device solves the problems, such as is similar to time-domain offset imaging method, because
The implementation of this time domain migration imaging device may refer to the implementation of time-domain offset imaging method, and overlaps will not be repeated.
Used below, term " unit " or " module " can realize the combination of the software and/or hardware of predetermined function.Although with
The lower described device of embodiment is preferably realized with software, but hardware, or the realization of the combination of software and hardware
May and it be contemplated.Figure 11 is a kind of structure diagram of the time-domain migration imaging device of the embodiment of the present invention, such as Figure 11
It is shown, including：Read module 1101, the first two way travel time computing module 1102, the second two way travel time computing module 1103 and into
As module 1104, the structure is illustrated below.
Read module 1101, for reading the seismic data of a seismic channel；
First two way travel time computing module 1102, for by double flat root formula, ask for the focus of the seismic channel into
When being walked as the first of road coarse grid imaging point, and coarse grid imaging point to the second of receiving point walk when, and when described first is walked
Two way travel time the sum of when being walked with described second as coarse grid imaging point；
Second two way travel time computing module 1103, for by cubic convolution interpolation formula, being imaged according to the coarse grid
The two way travel time of point, interpolation calculation go out the two way travel time of refined net imaging point；
Image-forming module 1104, for the two way travel time according to the refined net imaging point, the interrogation energy on the seismic channel,
And the energy of pickup is overlapped, the imaging results using stack result as the seismic channel.
In one embodiment, read module 1101 can include：Acquiring unit, for obtaining region to be analyzed
CMP trace gather seismic datas；Velocity analysis unit, for reading one of seismic data in the CMP trace gathers seismic data.
In one embodiment, above-mentioned double flat root formula is：
Wherein, t represents the two way travel time of coarse grid imaging point, tsRepresent focus to the first of imaging road coarse grid imaging point
When walking, tgRepresent coarse grid imaging point to the second of receiving point walk when, t0Represent the coarse grid imaging point it is corresponding it is vertical walk when
Depth, vrmsRepresent the root mean sequare velocity at the coarse grid imaging point, x represents that the coarse grid imaging point arrives common point
Horizontal distance, h represent half geophone offset.
In one embodiment, the formula of above-mentioned cubic convolution interpolation is：
tn=Tk-1(-s3+2s2-s)/2+Tk(3s3-5s2+2)/2
+Tk+1(-3s3+4s2+s)/2+Tk+2(s3-s2)/2
Wherein, tnRepresent the two way travel time of refined net imaging point obtained after interpolation calculation, Tk-1、Tk、Tk+1、Tk+2Represent institute
When the walking of four coarse grid imaging points near refined net imaging point is stated, s represents refined net step-length.
In one embodiment, above-mentioned coarse grid imaging point be according to predetermined Grid Sampling interval, to imaging road into
Row sparse sampling obtains.
In another embodiment, a kind of software is additionally provided, which is used to perform above-described embodiment and preferred reality
Apply the technical solution described in mode.
In another embodiment, a kind of storage medium is additionally provided, above-mentioned software is stored with the storage medium, should
Storage medium includes but not limited to：CD, floppy disk, hard disk, scratch pad memory etc..
It can be seen from the above description that the embodiment of the present invention realizes following technique effect：It is public by double flat root
Formula has tried to achieve the two way travel time of coarse grid imaging point, is then turned the two way travel time of coarse grid imaging point by cubic convolution interpolation
It is changed into the two way travel time of refined net imaging point, to meet the continuity of migration operator single order, so as to solve inclined in the prior art
The high-frequency information moved in stacked section is seriously damaged and influences subsequent treatment and the technical problem of inversion accuracy, has reached effective
Improve the technique effect of imaging results precision and accuracy.
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
Realized with general computing device, they can be concentrated on single computing device, or are distributed in multiple computing devices
On the network formed, alternatively, they can be realized with the program code that computing device can perform, it is thus possible to by it
Store and performed in the storage device by computing device, and in some cases, can be to be held different from order herein
They, are either fabricated to each integrated circuit modules or will be multiple in them by the shown or described step of row respectively
Module or step are fabricated to single integrated circuit module to realize.In this way, the embodiment of the present invention be not restricted 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 limit the invention, for the skill of this area
For art personnel, the embodiment of the present invention can have various modifications and variations.Within the spirit and principles of the invention, made
Any modification, equivalent substitution, improvement and etc., should all be included in the protection scope of the present invention.
Claims (8)
- A kind of 1. time-domain offset imaging method, it is characterised in that including：Read the seismic data of a seismic channel；By double flat root formula, ask for the focus of the seismic channel to the first of imaging road coarse grid imaging point walk when, and coarse net Lattice imaging point to the second of receiving point walk when, and as coarse grid imaging point the sum of when described first is walked and when described second walks Two way travel time, the double flat root formula is：<mrow> <mi>t</mi> <mo>=</mo> <msub> <mi>t</mi> <mi>s</mi> </msub> <mo>+</mo> <msub> <mi>t</mi> <mi>g</mi> </msub> <mo>=</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <mfrac> <msub> <mi>t</mi> <mn>0</mn> </msub> <mn>2</mn> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mfrac> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>+</mo> <mi>h</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <msubsup> <mi>v</mi> <mrow> <mi>r</mi> <mi>m</mi> <mi>s</mi> </mrow> <mn>2</mn> </msubsup> </mfrac> </mrow> </msqrt> <mo>+</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <mfrac> <msub> <mi>t</mi> <mn>0</mn> </msub> <mn>2</mn> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mfrac> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <mi>h</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <msubsup> <mi>v</mi> <mrow> <mi>r</mi> <mi>m</mi> <mi>s</mi> </mrow> <mn>2</mn> </msubsup> </mfrac> </mrow> </msqrt> </mrow>Wherein, t represents the two way travel time of coarse grid imaging point, tsRepresent focus to be imaged road coarse grid imaging point first walk when, tgRepresent coarse grid imaging point to the second of receiving point walk when, t0Represent the coarse grid imaging point it is corresponding it is vertical walk when depth, vrmsRepresent the root mean sequare velocity at the coarse grid imaging point, x represent the coarse grid imaging point to common point level away from From h represents half geophone offset；By cubic convolution interpolation formula, according to the two way travel time of the coarse grid imaging point, interpolation calculation goes out refined net imaging The two way travel time of point；According to the two way travel time of the refined net imaging point, the interrogation energy on the seismic channel, and the energy of pickup is folded Add, the imaging results using stack result as the seismic channel；Above-mentioned steps are repeated to all seismic channels, obtain the migration imaging result in an imaging road.
- 2. the method as described in claim 1, it is characterised in that the seismic data of a seismic channel is read, including：Obtain the CMP trace gather seismic datas in region to be analyzed；One of seismic data is read from the CMP trace gathers seismic data.
- 3. the method as described in claim 1, it is characterised in that the formula of the cubic convolution interpolation is：tn=Tk-1(-s3+2s2-s)/2+Tk(3s3-5s2+2)/2+Tk+1(-3s3+4s2+s)/2+Tk+2(s3-s2)/2Wherein, tnRepresent the two way travel time of refined net imaging point obtained after interpolation calculation, Tk-1、Tk、Tk+1、Tk+2Represent refined net When walking of four coarse grid imaging points near imaging point, s represent refined net step-length.
- 4. method as claimed any one in claims 1 to 3, it is characterised in that the coarse grid imaging point is according to predetermined Grid Sampling interval, sparse sampling is carried out to imaging road and is obtained.
- A kind of 5. time-domain migration imaging device, it is characterised in that including：Read module, for reading the seismic data of a seismic channel；First two way travel time computing module, for by double flat root formula, asking for the focus of the seismic channel to imaging road coarse net The first of lattice imaging point when walking, and coarse grid imaging point to the second of receiving point walk when, and when described first is walked and described Two two way travel times the sum of when walking as coarse grid imaging point, the double flat root formula are：<mrow> <mi>t</mi> <mo>=</mo> <msub> <mi>t</mi> <mi>s</mi> </msub> <mo>+</mo> <msub> <mi>t</mi> <mi>g</mi> </msub> <mo>=</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <mfrac> <msub> <mi>t</mi> <mn>0</mn> </msub> <mn>2</mn> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mfrac> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>+</mo> <mi>h</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <msubsup> <mi>v</mi> <mrow> <mi>r</mi> <mi>m</mi> <mi>s</mi> </mrow> <mn>2</mn> </msubsup> </mfrac> </mrow> </msqrt> <mo>+</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <mfrac> <msub> <mi>t</mi> <mn>0</mn> </msub> <mn>2</mn> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mfrac> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <mi>h</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <msubsup> <mi>v</mi> <mrow> <mi>r</mi> <mi>m</mi> <mi>s</mi> </mrow> <mn>2</mn> </msubsup> </mfrac> </mrow> </msqrt> </mrow>Wherein, t represents the two way travel time of coarse grid imaging point, tsRepresent focus to be imaged road coarse grid imaging point first walk when, tgRepresent coarse grid imaging point to the second of receiving point walk when, t0Represent the coarse grid imaging point it is corresponding it is vertical walk when depth, vrmsRepresent the root mean sequare velocity at the coarse grid imaging point, x represent the coarse grid imaging point to common point level away from From h represents half geophone offset；Second two way travel time computing module, for by cubic convolution interpolation formula, according to the round trip of the coarse grid imaging point When walking, interpolation calculation goes out the two way travel time of refined net imaging point；Image-forming module, for the two way travel time according to the refined net imaging point, the interrogation energy on the seismic channel, and will pickup Energy be overlapped, the imaging results using stack result as the seismic channel.
- 6. device as claimed in claim 5, it is characterised in that the read module includes：Acquiring unit, for obtaining the CMP trace gather seismic datas in region to be analyzed；Reading unit, for reading one of seismic data from the CMP trace gathers seismic data.
- 7. device as claimed in claim 5, it is characterised in that the formula of the cubic convolution interpolation is：tn=Tk-1(-s3+2s2-s)/2+Tk(3s3-5s2+2)/2+Tk+1(-3s3+4s2+s)/2+Tk+2(s3-s2)/2Wherein, tnRepresent the two way travel time of refined net imaging point obtained after interpolation calculation, Tk-1、Tk、Tk+1、Tk+2Represent refined net When walking of four coarse grid imaging points near imaging point, s represent refined net step-length.
- 8. the device as any one of claim 5 to 7, it is characterised in that the coarse grid imaging point is according to predetermined Grid Sampling interval, sparse sampling is carried out to imaging road and is obtained.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510577384.8A CN105607117B (en) | 2015-09-11 | 2015-09-11 | Time-domain offset imaging method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510577384.8A CN105607117B (en) | 2015-09-11 | 2015-09-11 | Time-domain offset imaging method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105607117A CN105607117A (en) | 2016-05-25 |
| CN105607117B true CN105607117B (en) | 2018-05-15 |
Family
ID=55987193
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510577384.8A Active CN105607117B (en) | 2015-09-11 | 2015-09-11 | Time-domain offset imaging method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105607117B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109143339B (en) * | 2018-08-14 | 2020-06-09 | 中国石油天然气集团有限公司 | Elastic reverse time migration imaging method and device based on transverse wave stress invariant |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2379271B (en) * | 2001-07-31 | 2004-10-27 | Pgs Americas Inc | True relative amplitude migration |
| CN102841379A (en) * | 2012-09-06 | 2012-12-26 | 中国石油大学(华东) | Method for analyzing pre-stack time migration and speed based on common scatter point channel set |
| CN102914791A (en) * | 2011-08-05 | 2013-02-06 | 中国石油天然气集团公司 | Kirchhoff prestack time migration method for processing seismic data of undulating surface |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7447113B2 (en) * | 2003-10-23 | 2008-11-04 | Pgs Americas, Inc. | Kirchhoff prestack time migration method for PS waves |
-
2015
- 2015-09-11 CN CN201510577384.8A patent/CN105607117B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2379271B (en) * | 2001-07-31 | 2004-10-27 | Pgs Americas Inc | True relative amplitude migration |
| CN102914791A (en) * | 2011-08-05 | 2013-02-06 | 中国石油天然气集团公司 | Kirchhoff prestack time migration method for processing seismic data of undulating surface |
| CN102841379A (en) * | 2012-09-06 | 2012-12-26 | 中国石油大学(华东) | Method for analyzing pre-stack time migration and speed based on common scatter point channel set |
Non-Patent Citations (4)
| Title |
|---|
| Cubic Convolution Interpolation for Digital Image Processing;ROBERT G.KEYS;《IEEE TRANSACTIONS ON ACOUSTICS,SPEECH, AND SIGNAL PROCESSING》;19811231;第ASSP-29卷(第6期);第1154页第2栏第38行-第1155页第2栏第35行 * |
| 克希霍夫叠前时间偏移走时算法应用效果评价;许建国;《中外能源》;20150531;第20卷(第5期);第36-40页 * |
| 基于CPU+GPU联合计算真地表叠前时间偏移实用化研究;马召贵等;《石油地球物理勘探》;20140630;第49卷(第3期);第531页第1栏第3-11行及第531页第2栏第34行-第532页第1栏第10行 * |
| 波前构建法研究现状;韩复兴等;《地球物理学进展》;20110630;第26卷(第3期);摘要及第1046页第1栏第2行-第1048页第1栏第3行 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105607117A (en) | 2016-05-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104237940B (en) | A kind of diffraction wave imaging method based on dynamic characteristic and device | |
| CN102914791B (en) | Kirchhoff prestack time migration method for processing seismic data of undulating surface | |
| CN105549080B (en) | A kind of relief surface waveform inversion method based on auxiliary coordinates | |
| CN107462924B (en) | A kind of absolute wave impedance inversion method independent of well-log information | |
| CN108387933A (en) | A kind of method, apparatus and system of definitely interval quality factors | |
| CN109001813A (en) | A kind of method, apparatus and system of multiple suppression | |
| CN112083487B (en) | Method and device for extracting broadband frequency dispersion curve | |
| CN103926623A (en) | Method for suppressing reverse time migration low frequency noise | |
| CN110471113A (en) | Bearing calibration, device and storage medium are moved in inverting based on unstable state seismic data | |
| CN106291681A (en) | Seismic horizon method for automatic tracking and device | |
| CN105911585A (en) | Method and device for extracting seismic record regular interference waves | |
| CN103744114B (en) | Method and device for estimating quality factor on basis of zero offset VSP (vertical seismic profile) data | |
| CN104330826A (en) | A method for removing various noises under the condition of complex surface | |
| CN108415073A (en) | Angle domain back scattering offset imaging method and device | |
| CN105607117B (en) | Time-domain offset imaging method and device | |
| CN106970417A (en) | Ellipse expansion converted-wave velocity analysis method and system | |
| CN108897036A (en) | A kind of seismic data processing technique and device | |
| CN105353406B (en) | A kind of method and apparatus for generating angle gathers | |
| CN107179546B (en) | A kind of method and device of determining velocity field | |
| CN105044769B (en) | The method for improving the resolution ratio of seismic signal | |
| CN106291691A (en) | A kind of seismic migration imaging method and device | |
| CN107918153B (en) | A kind of seismic signal coherence high-precision detecting method | |
| De Martin et al. | Impact of geometric effects on near‐surface Green’s functions | |
| CN104181600A (en) | Seismic data linear noise attenuation method and device | |
| CN107255833B (en) | The determination method and apparatus of Seismic Stacked Section |
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 |