CN101609163A - Multi-scale seismic data joint imaging method based on wave theory - Google Patents
Multi-scale seismic data joint imaging method based on wave theory Download PDFInfo
- Publication number
- CN101609163A CN101609163A CNA2009101600153A CN200910160015A CN101609163A CN 101609163 A CN101609163 A CN 101609163A CN A2009101600153 A CNA2009101600153 A CN A2009101600153A CN 200910160015 A CN200910160015 A CN 200910160015A CN 101609163 A CN101609163 A CN 101609163A
- Authority
- CN
- China
- Prior art keywords
- data
- continuation
- wave
- wave field
- seismic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
Multi-scale seismic data joint imaging method based on wave theory, comprise: to surface-seismic data, cross-well seismic data and/or the vertical seismic profiling (VSP) data of a block is that at least two kind seismic datas of VSP method in all adopt the direct downward method of wave field to measure respectively, obtains the pre-stack depth migration data; The pre-stack depth migration data that surface seismic, crosshole seismic and/or VSP obtained described at least two kinds of seismic datas are respectively carried out normalized respectively; Pre-stack depth migration data after the normalized are superposeed, obtain the multi-scale data combine imaging results of this block.Method of the present invention can combine the advantage of surface-seismic data, VSP data and cross-well seismic data, the horizontal change of speed-adaptive field preferably, and keep the dynamic characteristic of certain degree.
Description
Technical field
The present invention relates to field of seismic exploration, relate in particular to a kind of multi-scale seismic data joint imaging method based on wave theory.
Background technology
Seismic prospecting is an important means of surveying oil and natural gas source before the probing.At aspects such as coalfield and engineering geological investigation, Study on Regional Geology and earth's crust researchs, seismic prospecting also is used widely.The plurality of seismic exploration method is arranged at present, for example surface seismic method, seismic crosshole method and vertical seismic profiling (VSP) (VSP, Vertical Seismic Profiling) measuring method.
But, there is following problem in the above-mentioned prior art at least:
Surface-seismic data is because the earthquake wave frequency is lower, and the resolution of imaging is restricted;
Though VSP and cross-well seismic data resolution height are merely able to provide the image-forming information of two mouthfuls of subranges between the well;
Ununified at present effective method carries out joint imaging to the multi-scale data of surface-seismic data, cross-well seismic data and VSP data.
Summary of the invention
The objective of the invention is to, a kind of method is provided, the multi-scale data of surface-seismic data, cross-well seismic data and VSP data can be carried out joint imaging, thereby the advantage of surface-seismic data, VSP data and cross-well seismic data is combined, improve image quality, realize multi-scale seismic data joint imaging based on wave theory.
For this reason, the invention provides a kind of multi-scale seismic data joint imaging method, comprising based on wave theory:
Adopt the direct downward method of wave field to measure respectively at least two kinds of seismic datas in surface-seismic data, cross-well seismic data and the vertical seismic profiling (VSP) data of a block, obtain the pre-stack depth migration data;
Respectively the pre-stack depth migration data of described at least two kinds of seismic datas are carried out normalized;
Pre-stack depth migration data after the normalized are superposeed, obtain the multi-scale data combine imaging results of this block.
As preferably, the direct downward method of wherein said wave field comprises:
With the wave field in the seismic data from a selected level reference downward continuation;
On each predetermined depth step of continuation, the wave field value on the landform face that is intercepted is added in the continuation wave field.
As preferably, described at least two kinds of seismic datas comprise cross-well seismic data or vertical seismic profiling (VSP) data, and the direct downward method of then described wave field comprises:
Cross-well seismic data or vertical seismic profiling (VSP) data are begun from the face of land to carry out wave field extrapolation, and each runs on the depth layer that the wave field value record is arranged in continuation, and the wave field value and continuation to the wave field value of this layer of this depth layer place record are carried out addition.
As preferably, adopt the direct downward method of wave field to be specially to the vertical seismic profiling (VSP) data:
From the face of land, the wave field that source wavefield and wave detector are received is downward continuation simultaneously;
Continuation is when having new wave detector position, and wave field value that continuation is calculated and wave field value addition that should new wave detector place record obtain the wave field value on the horizontal level of this new wave detector correspondence, utilizes the dependent imaging condition to be extracted into then to look like value to carry out imaging.
As preferably, adopt the direct downward method of wave field to comprise the shot record migration pre-stack depth migration imaging to cross-well seismic data, comprising:
The selected reference face, the initialization wave field value, and with the wave field downward continuation;
In each step size, detect and whether have wave detector, if detect wave detector, the then wave field value addition of the wave field value that continuation is calculated and this new wave detector place record, obtain the wave field value on the horizontal level of this new wave detector correspondence, thereby finish the continuation of the pairing depth layer of this wave detector, continue the depth layer of downward continuation then up to all wave detector places of the intact described cross-well seismic data of continuation;
Utilize the dependent imaging condition to be extracted into the picture value and to be added to corresponding lateral attitude at every turn after the continuation, the depth layer of all wave detectors is all carried out continuation and is extracted into the picture value, thereby realized the pre-stack depth migration imaging of a common-shot-point records;
All shot points to cross-well seismic data all carry out above-mentioned continuation imaging process, and the imaging results with all common big gun records superposes then, realizes pre-stack depth migration imaging.
As preferably, adopt the direct downward method of wave field also to comprise the common receiver pre-stack depth migration imaging to cross-well seismic data, comprising:
The selected reference face, the initialization wave field value, and with the wave field downward continuation;
In each step size, detect and whether have shot point, if detect shot point, the then wave field value addition of the wave field value that continuation is calculated and this new shot point place record, obtain the wave field value on the horizontal level of this new shot point correspondence, thereby finish the continuation of the pairing depth layer of this shot point, continue the depth layer of downward continuation then up to all shot point places of the intact described cross-well seismic data of continuation;
Utilize the dependent imaging condition to be extracted into the picture value and to be added to corresponding lateral attitude at every turn after the continuation, the depth layer of all shot points is all carried out continuation and is extracted into the picture value, thereby realized a pre-stack depth migration imaging that common receiver writes down;
All geophone stations to cross-well seismic data all carry out above-mentioned continuation imaging process, and the imaging results with all common receiver records superposes then, realizes pre-stack depth migration imaging.
As preferably, described pre-stack depth migration data after the normalized are superposeed is specially: the pre-stack depth migration data after the normalized are carried out the adaptive optimization stack.
As preferably, described adaptive optimization stack comprises:
Superpose by different scale seismic trace, obtain standard track certain point on the same depth layer of the pre-stack depth migration data after the described normalized;
Obtain weighting coefficient and then obtain the stacking image value of this point on this same depth layer according to above-mentioned standard track;
Each point on this same depth layer is carried out iteration, obtain the multi-scale data combine imaging results of this same depth layer;
Each layer through the each point iteration carried out iteration, obtain the multi-scale data combine imaging results of this block.
Method of the present invention can combine the advantage of surface-seismic data, VSP data and cross-well seismic data, the horizontal change of speed-adaptive field preferably, and keep the dynamic characteristic of certain degree.
Description of drawings
Fig. 1 is the process flow diagram of an embodiment of the multi-scale seismic data joint imaging method based on wave theory provided by the invention;
Fig. 2 is the process flow diagram of another embodiment of the multi-scale seismic data joint imaging method based on wave theory provided by the invention;
Fig. 3 is the process flow diagram of an embodiment of the multi-scale seismic data joint imaging method based on wave theory provided by the invention;
Fig. 4 is the process flow diagram of another embodiment of the multi-scale seismic data joint imaging method based on wave theory provided by the invention;
Fig. 5 is the process flow diagram of another embodiment of the multi-scale seismic data joint imaging method based on wave theory provided by the invention;
Fig. 6 is the process flow diagram of an embodiment of the multi-scale seismic data joint imaging method based on wave theory provided by the invention;
Fig. 7 is the seismogeology illustraton of model;
Fig. 8, Fig. 9, Figure 10 are respectively surface seismic, VSP, crosshole seismic pre-stack depth migration imaging result;
Figure 11 is the imaging results based on the multi-scale seismic data joint imaging method of wave theory.
Embodiment
Below in conjunction with the accompanying drawing in the embodiment of the invention, the technical scheme in the embodiment of the invention is carried out clear, complete description, obviously, described embodiment only is the present invention's part embodiment, rather than whole embodiment.Based on the embodiment among the present invention, those of ordinary skills belong to the scope of protection of the invention not making the every other embodiment that is obtained under the creative work prerequisite.
Fig. 1 is the process flow diagram of an embodiment of the multi-scale seismic data joint imaging method based on wave theory provided by the invention, Fig. 2 is the process flow diagram of another embodiment of the multi-scale seismic data joint imaging method based on wave theory provided by the invention, as shown in Figure 1, the method for this embodiment comprises:
S01: adopt the direct downward method of wave field to measure respectively at least two kinds of seismic datas in surface-seismic data, cross-well seismic data and the vertical seismic profiling (VSP) data of a block, obtain the pre-stack depth migration data;
S02: respectively the pre-stack depth migration data of described at least two kinds of seismic datas are carried out normalized;
S03: the pre-stack depth migration data after the normalized are superposeed, obtain the multi-scale data combine imaging results of this block.
Wherein, at first adopt the direct downward method of wave field to measure respectively at least two in surface-seismic data, cross-well seismic data and/or the vertical seismic profiling (VSP) data (being the VSP data), obtain the pre-stack depth migration data, for example, adopt the direct downward method of wave field to measure respectively to surface-seismic data, cross-well seismic data and vertical seismic profiling (VSP) data, obtain the pre-stack depth migration data.As shown in Figure 2, the direct downward method of wave field is " progressively-add up " the direct downward method of wave field, be about to initial focus or detection wave field (being the wave field that wave detector receives) from a level reference downward continuation, and on each depth step, focus on the existing landform face or detection wave field value are added in the wave field (initial wave field and the new total wave field that adds wave field) of continuation, so just can carry out migration before stack to the record on the rugged topography.And, crosshole seismic, VSP can regard a kind of observed pattern of similar relief surface as, begin to carry out wave field extrapolation from the face of land, run on the depth layer that wave field value is arranged at each, the wave field value and continuation to the wave field value of this layer at this depth layer place have been carried out addition, thereby measurement seismic migrations such as surface seismic, crosshole seismic, VSP are unified on method.
Then, respectively the pre-stack depth migration data of described at least two kinds of seismic datas are carried out normalized, for example, the pre-stack depth migration data that surface seismic, crosshole seismic and/or VSP obtain are respectively carried out normalized.Since surface seismic and pit shaft earthquake (comprising VSP and crosshole seismic etc.) when gathering the pre-stack depth migration data because of exciting the influence with factor such as detecting action to have shot point and the inconsistent and unmatched problem of acceptance point energy, when distinct methods joint imagings such as realization surface-seismic data, VSP data, RVSP (anti-vertical seismic profiling (VSP)) data, cross-well seismic data, the information that certainly exists different frequency has different effective energies under the different acquisition condition.If directly superpose joint imaging to certainly exist the inconsistent and unmatched problem of energy the data imaging result under these different recording geometrys, therefore different wave fields normalized be need be carried out, joint imaging result's correctness and result's stability guaranteed.At each method obtain pre-stack depth migration data (comprising surface seismic, crosshole seismic and VSP etc.), for the order of magnitude that makes various results is unified, adopted normalized, for the pre-stack depth migration section under every kind of observed pattern, for each point (x z) locates, and has:
X wherein
X, izBe the imaging value at this some place after the normalization, u
Ix, izBe the imaging value behind the depth shift, ix is the lateral attitude, and iz is a depth location transversely.
At last, the pre-stack depth migration data after the normalized are superposeed, obtain the multi-scale data combine imaging results of this block.A kind of adaptive optimization that the pre-stack depth migration data of the surface seismic after the normalization soon, VSP, cross-well seismic data are carried out superposes, and then obtains the multi-scale data combine imaging results of this block.
Fig. 3 is the process flow diagram of an embodiment of the multi-scale seismic data joint imaging method based on wave theory provided by the invention, as shown in Figure 3, in said method, vertical seismic profiling (VSP) data VSP is adopted the direct downward method of wave field, comprising:
S011a: from the face of land, the wave field that source wavefield and wave detector are received is downward continuation simultaneously;
S012a: continuation is when having new wave detector position, wave field value that continuation is calculated and the wave field value addition of should new wave detector place writing down, obtain the wave field value on the horizontal level of this new wave detector correspondence, utilize the dependent imaging condition to be extracted into the picture value then and carry out imaging;
S02: respectively the pre-stack depth migration data of described at least two kinds of seismic datas are carried out normalized;
S03: the pre-stack depth migration data after the normalized are superposeed, obtain the multi-scale data combine imaging results of this block.
Among this embodiment, concrete narration is adopted the direct downward method of wave field to vertical seismic profiling (VSP) data VSP.For single big gun skew, the method that progressively adds up based on wave field, from the face of land, the wave field that source wavefield and wave detector receive is downward continuation simultaneously, when wave field extrapolation when new wave detector position is arranged, wave field value that continuation is calculated and the wave field value addition of should new wave detector place writing down obtain the wave field value on the horizontal level of this new wave detector correspondence, utilize the dependent imaging condition to be extracted into the picture value then and carry out dependent imaging, finally obtain the result of the single big gun skew of VSP; Carry out the migration before stack imaging of single big gun successively, thereby obtain the result of the pre-stack depth migration of VSP.
The pre-stack depth migration data that obtain with surface-seismic data and/or cross-well seismic data are carried out normalized then, and (step S02, S03 are as previously mentioned) again superposes the pre-stack depth migration data after the normalized.
Fig. 4 is the process flow diagram of another embodiment of the multi-scale seismic data joint imaging method based on wave theory provided by the invention, as shown in Figure 4, in the said method, adopt the direct downward method of wave field to comprise the shot record migration pre-stack depth migration imaging, comprising cross-well seismic data:
S011b: the selected reference face, the initialization wave field value, and with the wave field downward continuation;
S012b: in each step size, detect whether there is wave detector, if detect wave detector, the then wave field value addition of the wave field value that continuation is calculated and this new wave detector place record, obtain the wave field value on the horizontal level of this new wave detector correspondence, thereby finish the continuation of the pairing depth layer of this wave detector, continue the depth layer of downward continuation then up to all wave detector places of the intact described cross-well seismic data of continuation;
S013b: utilize the dependent imaging condition to be extracted into the picture value and to be added to corresponding lateral attitude at every turn after the continuation, the depth layer of all wave detectors is all carried out continuation and is extracted into the picture value, thereby realized the pre-stack depth migration imaging of a common-shot-point records;
S014b: all shot points to cross-well seismic data all carry out above-mentioned continuation imaging process, and the imaging results with all common big gun records superposes then, realizes pre-stack depth migration imaging;
S02: respectively the pre-stack depth migration data of described at least two kinds of seismic datas are carried out normalized;
S03: the pre-stack depth migration data after the normalized are superposeed, obtain the multi-scale data combine imaging results of this block.
Adopt the direct downward method of wave field to comprise shot record migration pre-stack depth migration imaging and common receiver pre-stack depth migration imaging to cross-well seismic data, among this embodiment, concrete narration adopts the direct downward method of wave field to comprise the shot record migration pre-stack depth migration imaging to cross-well seismic data.At first, the selected reference face, the initialization wave field value is initialize 0 for example, begin degree of depth continuation then, in each step size, all detect whether wave detector is arranged, the wave field value addition that the wave field value that then continuation calculated is arranged and should new wave detector place write down, obtain the wave field value on the horizontal level of this new wave detector correspondence, thereby finish the continuation of the depth layer of this new wave detector, then do not carry out next step wave field extrapolation, continue the depth layer at downward continuation all wave detector places of seismic data up between the continuation completion then; After each continuation, use the dependent imaging condition (for example in frequency field, source wavefield and detection wave field are multiplied each other, and all frequency content additions) be extracted into the picture value be the different frequency composition add with after amplitude, and being put into ground level position in the place rate pattern that corresponding lateral attitude is the imaging point correspondence, the depth layer at all wave detector places of seismic data obtains the pre-stack depth migration imaging of a common-shot-point records between the continuation completion.Every shot point is carried out top process, and the imaging results with all common-shot-point records superposes then, finally obtains well earthquake reflected wave pre-stack depth migration imaging result.
Then with surface-seismic data and/or cross-well seismic data to the pre-stack depth migration data carry out normalized, (step S02, S03 are as previously mentioned) again superposes the pre-stack depth migration data after the normalized.
Fig. 5 is the process flow diagram of another embodiment of the multi-scale seismic data joint imaging method based on wave theory provided by the invention, as shown in Figure 5, in the said method, adopt the direct downward method of wave field to comprise the common receiver pre-stack depth migration imaging, comprising cross-well seismic data:
S011c: the selected reference face, the initialization wave field value, and with the wave field downward continuation;
S012c: in each step size, detect whether there is shot point, if detect shot point, the then wave field value addition of the wave field value that continuation is calculated and this new shot point place record, obtain the wave field value on the horizontal level of this new shot point correspondence, thereby finish the continuation of the pairing depth layer of this shot point, continue the depth layer of downward continuation then up to all shot point places of the intact described cross-well seismic data of continuation;
S013c: utilize the dependent imaging condition to be extracted into the picture value and to be added to corresponding lateral attitude at every turn after the continuation, the depth layer of all shot points is all carried out continuation and is extracted into the picture value, thereby realized a pre-stack depth migration imaging that common receiver writes down;
S014c: all geophone stations to cross-well seismic data all carry out above-mentioned continuation imaging process, and the imaging results with all common receiver records superposes then, realizes pre-stack depth migration imaging;
S02: respectively the pre-stack depth migration data of described at least two kinds of seismic datas are carried out normalized;
S03: the pre-stack depth migration data after the normalized are superposeed, obtain the multi-scale data combine imaging results of this block.
Among this embodiment, concrete narration adopts the direct downward method of wave field to comprise the common receiver pre-stack depth migration imaging to cross-well seismic data.At first, the selected reference face, the initialization wave field value is initialize 0 for example, begin degree of depth continuation then, in each step size, detect whether there is shot point, if detect shot point, the then wave field value addition of the wave field value that continuation is calculated and this new shot point place record, obtain the wave field value on the horizontal level of this new shot point correspondence, thereby finish the continuation of the pairing depth layer of this shot point, continue the depth layer of downward continuation then up to all shot point places of the intact described cross-well seismic data of continuation; After each continuation, use the dependent imaging condition (for example in frequency field, source wavefield and detection wave field are multiplied each other, and all frequency content additions) be extracted into the picture value be the different frequency composition add with after amplitude, and being put into ground level position in the place rate pattern that corresponding lateral attitude is the imaging point correspondence, the depth layer at all shot point places of seismic data obtains the pre-stack depth migration imaging of a common receiver record between the continuation completion.Each geophone station is carried out top process, and the imaging results with all common receiver records superposes then, finally obtains well earthquake reflected wave pre-stack depth migration imaging result.
Then with surface-seismic data and/or cross-well seismic data to the pre-stack depth migration data carry out normalized, (step S02, S03 are as previously mentioned) again superposes the pre-stack depth migration data after the normalized.
By as can be seen above-mentioned, the skew principle of cross-well seismic data and VSP data is similar, all carry out with the different wave field extrapolation cumulative process of two kinds of data that are based on the prestack depth migration method of relief surface of routine, and surface seismic carries out at Different Ground position different depth at same ground location different depth.Both considered the kinematics character (temporal information) of crosshole seismic based on the prestack depth migration method of the VSP data of wave theory and cross-well seismic data, considered its dynamic characteristic (amplitude information) again, thereby it is more suitable for complex structure and lateral velocity variation medium, higher with respect to its imaging precision of conventional method, the relative maintenance of seismic amplitude is better.
Fig. 6 is the process flow diagram of an embodiment of the multi-scale seismic data joint imaging method based on wave theory provided by the invention, as shown in Figure 6, comprises for above-mentioned:
S01: adopt the direct downward method of wave field to measure respectively at least two kinds of seismic datas in surface-seismic data, cross-well seismic data and the vertical seismic profiling (VSP) data of a block, obtain the pre-stack depth migration data;
S02: respectively the pre-stack depth migration data of described at least two kinds of seismic datas are carried out normalized;
S031: superpose by different scale seismic trace, obtain standard track with certain point on the same depth layer of the pre-stack depth migration data after the described normalized;
S032: obtain weighting coefficient and then obtain the stacking image value of this point on this same depth layer according to above-mentioned standard track;
S033: the each point on this same depth layer is carried out iteration, obtain the multi-scale data combine imaging results of this same depth layer;
S034: each layer through the each point iteration carried out iteration, obtain the multi-scale data combine imaging results of this block.
At first adopt the direct downward method of wave field to measure respectively at least two kinds of seismic datas in surface-seismic data, cross-well seismic data and the vertical seismic profiling (VSP) data of a block, obtain the pre-stack depth migration data, respectively the pre-stack depth migration data of described two kinds of seismic datas are carried out normalized (step S01, S02 such as above-mentioned) then at least.
Superpose by different scale seismic trace then, seek out a standard track X certain point on the same depth layer of the pre-stack depth migration data after the described normalized
∑ jz, promptly
X
∑jz＝X
surf ix，iz+X
vsp ix，iz+X
rvsp ix，iz+X
crosswell ix，iz(ix＝1，2，…，nx)
Although wherein standard track is to consider that surface seismic, VSP, RVSP and crosshole seismic have the characteristics of himself data acquisition and the otherness on the local wave field, but for characterizing same underground geologic bodies, the seismic wave field that various acquisition methods collect has its many-sided general character on its kinematics and dynamic characteristic, so the quality that the seismic data sum-average arithmetic that adopts various acquisition methods to obtain to earthquake signal to noise ratio (S/N ratio) block with higher in a certain exploratory area is weighed the seismic data of each method and technology collection as the earthquake standards road of this block is rationally with reliably.The purpose in production standard road is to be recorded in and to participate in the effectively composition of stack on time and the space in order to control seismic trace that various different acquisition methods obtain, make the composition of the measured road of matter participation stack many, and the composition of ropy road participation stack is few, and the road of poor quality does not participate in stack.
Utilize this standard track to ask for weighting coefficient then,
The stacking image value of certain point can be expressed as on a certain depth layer position:
The rest may be inferred, carries out iteration with certain layer each point, just obtained the multi-scale data combine imaging results of this layer; Further, each layer carried out iteration, just obtained the multi-scale data combine imaging results of this block.
Method of the present invention, the horizontal change of speed-adaptive field preferably, the data volume formation method of different scale is unified, compare with existing method and to have kept dynamic characteristic to a certain degree, and, the self-adaptation stack of the pre-stack depth migration data after the normalization has solved the problem of multi-scale data combine imaging, makes the advantage of different scale data imaging be able to complementation.
Fig. 7 is the seismogeology illustraton of model, and Fig. 8, Fig. 9, Figure 10 are respectively surface seismic, VSP, crosshole seismic pre-stack depth migration imaging result, and Figure 11 is the imaging results based on the multi-scale seismic data joint imaging method of wave theory.
As Fig. 7-shown in Figure 11, not only keep preferably unitary construction based on the imaging results of the multi-scale seismic data joint imaging method of wave theory but also fully shown high-frequency information, improved the resolution of imaging, obtain geological effect preferably.
In a word, the above is embodiments of the invention only, is not to be used to limit protection scope of the present invention, but is used to illustrate the present invention.Within the spirit and principles in the present invention all, any modification of being done, be equal to replacement, improvement etc., all should be included within protection scope of the present invention.
Claims (8)
1, based on the multi-scale seismic data joint imaging method of wave theory, comprising:
Adopt the direct downward method of wave field to measure respectively at least two kinds of seismic datas in surface-seismic data, cross-well seismic data and the vertical seismic profiling (VSP) data of a block, obtain the pre-stack depth migration data;
Respectively the pre-stack depth migration data of described at least two kinds of seismic datas are carried out normalized;
Pre-stack depth migration data after the normalized are superposeed, obtain the multi-scale data combine imaging results of this block.
2, the multi-scale seismic data joint imaging method based on wave theory as claimed in claim 1 is characterized in that, the direct downward method of wherein said wave field comprises:
With the wave field in the seismic data from a selected level reference downward continuation;
On each predetermined depth step of continuation, the wave field value on the landform face that is intercepted is added in the continuation wave field.
3, the multi-scale seismic data joint imaging method based on wave theory as claimed in claim 1 is characterized in that, described at least two kinds of seismic datas comprise cross-well seismic data or vertical seismic profiling (VSP) data, and the direct downward method of then described wave field comprises:
Cross-well seismic data or vertical seismic profiling (VSP) data are begun from the face of land to carry out wave field extrapolation, and each runs on the depth layer that the wave field value record is arranged in continuation, and the wave field value and continuation to the wave field value of this layer of this depth layer place record are carried out addition.
4, the multi-scale seismic data joint imaging method based on wave theory as claimed in claim 3 is characterized in that, adopts the direct downward method of wave field to be specially to the vertical seismic profiling (VSP) data:
From the face of land, the wave field that source wavefield and wave detector are received is downward continuation simultaneously;
Continuation is when having new wave detector position, and wave field value that continuation is calculated and wave field value addition that should new wave detector place record obtain the wave field value on the horizontal level of this new wave detector correspondence, utilizes the dependent imaging condition to be extracted into then to look like value to carry out imaging.
5, the multi-scale seismic data joint imaging method based on wave theory as claimed in claim 3 is characterized in that, adopts the direct downward method of wave field to comprise the shot record migration pre-stack depth migration imaging to cross-well seismic data, comprising:
The selected reference face, the initialization wave field value, and with the wave field downward continuation;
In each step size, detect and whether have wave detector, if detect wave detector, the then wave field value addition of the wave field value that continuation is calculated and this new wave detector place record, obtain the wave field value on the horizontal level of this new wave detector correspondence, thereby finish the continuation of the pairing depth layer of this wave detector, continue the depth layer of downward continuation then up to all wave detector places of the intact described cross-well seismic data of continuation;
Utilize the dependent imaging condition to be extracted into the picture value and to be added to corresponding lateral attitude at every turn after the continuation, the depth layer of all wave detectors is all carried out continuation and is extracted into the picture value, thereby realized the pre-stack depth migration imaging of a common-shot-point records;
All shot points to cross-well seismic data all carry out above-mentioned continuation imaging process, and the imaging results with all common big gun records superposes then, realizes pre-stack depth migration imaging.
6, the multi-scale seismic data joint imaging method based on wave theory as claimed in claim 3 is characterized in that, adopts the direct downward method of wave field also to comprise the common receiver pre-stack depth migration imaging to cross-well seismic data, comprising:
The selected reference face, the initialization wave field value, and with the wave field downward continuation;
In each step size, detect and whether have shot point, if detect shot point, the then wave field value addition of the wave field value that continuation is calculated and this new shot point place record, obtain the wave field value on the horizontal level of this new shot point correspondence, thereby finish the continuation of the pairing depth layer of this shot point, continue the depth layer of downward continuation then up to all shot point places of the intact described cross-well seismic data of continuation;
Utilize the dependent imaging condition to be extracted into the picture value and to be added to corresponding lateral attitude at every turn after the continuation, the depth layer of all shot points is all carried out continuation and is extracted into the picture value, thereby realized a pre-stack depth migration imaging that common receiver writes down;
All geophone stations to cross-well seismic data all carry out above-mentioned continuation imaging process, and the imaging results with all common receiver records superposes then, realizes pre-stack depth migration imaging.
7, the multi-scale seismic data joint imaging method based on wave theory as claimed in claim 1, it is characterized in that described pre-stack depth migration data after the normalized are superposeed is specially: the pre-stack depth migration data after the normalized are carried out the adaptive optimization stack.
8, the multi-scale seismic data joint imaging method based on wave theory as claimed in claim 7 is characterized in that, described adaptive optimization stack comprises:
Superpose by different scale seismic trace, obtain standard track certain point on the same depth layer of the pre-stack depth migration data after the described normalized;
Obtain weighting coefficient and then obtain the stacking image value of this point on this same depth layer according to above-mentioned standard track;
Each point on this same depth layer is carried out iteration, obtain the multi-scale data combine imaging results of this same depth layer;
Each layer through the each point iteration carried out iteration, obtain the multi-scale data combine imaging results of this block.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2009101600153A CN101609163B (en) | 2009-07-17 | 2009-07-17 | Multi-scale seismic data joint imaging method based on fluctuation theory |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2009101600153A CN101609163B (en) | 2009-07-17 | 2009-07-17 | Multi-scale seismic data joint imaging method based on fluctuation theory |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101609163A true CN101609163A (en) | 2009-12-23 |
| CN101609163B CN101609163B (en) | 2012-05-09 |
Family
ID=41482988
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2009101600153A Active CN101609163B (en) | 2009-07-17 | 2009-07-17 | Multi-scale seismic data joint imaging method based on fluctuation theory |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101609163B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102012517A (en) * | 2010-09-29 | 2011-04-13 | 北京吉星吉达科技有限公司 | Underground medium imaging method and device |
| CN103907032A (en) * | 2011-11-07 | 2014-07-02 | 界标制图有限公司 | Seismic imaging systems and methods employing correlation-based stacking |
| CN104570124A (en) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | Continuation imaging method suitable for cross-well seismic large-angle reflection conditions |
| CN104620132B (en) * | 2012-03-09 | 2017-03-15 | 国际壳牌研究有限公司 | For inverse time migration shoots the signal enhancing (DeSSeRT) that piles up using diversity |
| CN109100784A (en) * | 2018-06-08 | 2018-12-28 | 恒泰艾普（北京）能源科技研究院有限公司 | All-wave field imaging method is exchanged in the inspection of three-dimensional VSP source |
| CN112558160A (en) * | 2020-11-06 | 2021-03-26 | 中国海洋石油集团有限公司 | Azimuth difference three-dimensional seismic prestack fusion processing method and system |
-
2009
- 2009-07-17 CN CN2009101600153A patent/CN101609163B/en active Active
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102012517A (en) * | 2010-09-29 | 2011-04-13 | 北京吉星吉达科技有限公司 | Underground medium imaging method and device |
| CN103907032A (en) * | 2011-11-07 | 2014-07-02 | 界标制图有限公司 | Seismic imaging systems and methods employing correlation-based stacking |
| CN104620132B (en) * | 2012-03-09 | 2017-03-15 | 国际壳牌研究有限公司 | For inverse time migration shoots the signal enhancing (DeSSeRT) that piles up using diversity |
| US9939542B2 (en) | 2012-03-09 | 2018-04-10 | Shell Oil Company | Signal enhancement using diversity shot stacking for reverse time migrations (DeSSeRT) |
| CN104570124A (en) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | Continuation imaging method suitable for cross-well seismic large-angle reflection conditions |
| CN109100784A (en) * | 2018-06-08 | 2018-12-28 | 恒泰艾普（北京）能源科技研究院有限公司 | All-wave field imaging method is exchanged in the inspection of three-dimensional VSP source |
| CN109100784B (en) * | 2018-06-08 | 2020-04-28 | 恒泰艾普（北京）能源科技研究院有限公司 | Three-dimensional VSP source detection interchange full-wave-field imaging method |
| CN112558160A (en) * | 2020-11-06 | 2021-03-26 | 中国海洋石油集团有限公司 | Azimuth difference three-dimensional seismic prestack fusion processing method and system |
| CN112558160B (en) * | 2020-11-06 | 2024-05-28 | 中国海洋石油集团有限公司 | Direction difference three-dimensional earthquake prestack fusion processing method and system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101609163B (en) | 2012-05-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107526101B (en) | A kind of acquisition and processing method obtaining earthquake reflected wave | |
| CN104570125B (en) | A kind of method utilizing well data to improve image taking speed model accuracy | |
| AU612552B2 (en) | Model-based depth processing of seismic data | |
| CN101551466B (en) | Method for improving prediction precision of oil and gas reservoir by using seismic attribute related to offset distance | |
| US4779237A (en) | Method of geophysical exploration including processing and displaying seismic data to obtain a measure of subterranean formation rock properties | |
| CN101609163B (en) | Multi-scale seismic data joint imaging method based on fluctuation theory | |
| CN104216014A (en) | Seismic signal frequency division processing method | |
| CN103809216B (en) | A kind of resistivity data is combined speed with geological data and is built field method | |
| AU1336800A (en) | Hydrocarbon edge detection using seismic amplitude | |
| Zaslavsky et al. | Amplification effects from earthquakes and ambient noise in the Dead Sea rift (Israel) | |
| CN103645503A (en) | Three-dimensional time domain illumination analysis and amplitude compensation method | |
| CN103576200A (en) | Low signal-to-noise ratio zone shallow wave impedance interface static correction method | |
| CN100349006C (en) | Method for determining distribution of lithologic character and liquid by using inversion technique of wave impedance | |
| Brodic et al. | Three-component seismic land streamer study of an esker architecture through S-and surface-wave imaging | |
| CN104570116A (en) | Geological marker bed-based time difference analyzing and correcting method | |
| CN107656308B (en) | A kind of common scattering point pre-stack time migration imaging method based on time depth scanning | |
| CN102565852B (en) | Angle domain pre-stack offset data processing method aiming to detect oil-gas-bearing property of reservoir | |
| Farrell et al. | Refraction statics | |
| CN111308558B (en) | Shale gas horizontal well longitudinal wave time difference correction method | |
| Nunziata | Seismic ground motion in Napoli for the 1980 Irpinia earthquake | |
| Ivanov et al. | Near-surface evaluation of Ball Mountain Dam, Vermont, using multi-channel analysis of surface waves (MASW) and refraction tomography seismic methods on land-streamer data | |
| Mazanec et al. | SURFACE WAVES AS A COST-EFFECTIVE TOOL FOR ENHANCING THE INTERPRETATION OF SHALLOW REFRACTION SEISMIC DATA. | |
| Schermerhorn | Fundamental Resonant Frequencies Derived from Shallow Sediment Properties for the Charleston, South Carolina Area | |
| Moran | Three-dimensional P-wave velocity structure in the greater Mount Rainier area from local earthquake tomography | |
| Cartmell et al. | An example of the use of crosshole tomography in dam wall foundation studies |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |