CN106597536A - Processing method and device of seismic imaging - Google Patents
Processing method and device of seismic imaging Download PDFInfo
- Publication number
- CN106597536A CN106597536A CN201611116919.2A CN201611116919A CN106597536A CN 106597536 A CN106597536 A CN 106597536A CN 201611116919 A CN201611116919 A CN 201611116919A CN 106597536 A CN106597536 A CN 106597536A
- Authority
- CN
- China
- Prior art keywords
- signal
- cymoscope
- benchmark
- direct
- detector
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
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. analysis, for interpretation, for correction
- G01V1/282—Application of seismic models, synthetic seismograms
Abstract
The invention provides a processing method and device of seismic imaging. The method comprises that 1) a reference wave detector is arranged in an excitation well in a preset distance to an object; 2) a reflection wave signal received by a first wave detector in the excitation well is obtained; 3) a direct wave signal excited by a first seismic center and received by the reference wave detector is obtained, and correlation between the direct wave signal and the reflection wave signal is calculated; and 4) reflection wave signals excited by other seismic centers are received from the first wave detector are obtained one by one, the steps 2) and 3) are repeated, correlation signals between the direct waves received by the reference wave detector and the reflection wave signals are calculated, a sum of the correlation signals is calculated, and a reconstruction signal of the reference wave detector is generated. The method and device can reduce the sensitivity of static correction and a near-land-surface speed model, and improve the imaging precision of a high steep structure.
Description
Technical field
The present invention relates to technical field of geophysical exploration, more particularly to a kind of processing method and processing device of seismic imaging.
Background technology
With the deep development of seismic prospecting, mid-deep strata, ultra deep oil-gas exploration have got growing concern for.It is deep
The imaging tight association that layer oil-gas exploration is constructed with earth formation deep, and the accuracy of infrastructure imaging is by the shadow of many factors
Ring, static correction, near-surface velocity modeling, offset parameter are arranged etc..Wherein, near-surface velocity modeling factors are related to speed
The key technologies such as modeling quality, static correction precision, imaging deep.High-dip structure is a class common in infrastructure, is accurately visited
High-dip structure is surveyed for petroleum industry assessment trap, the scope of Reservoir Body, oil-gas migration situation and nuclear industry process underground radiation
Property rubbish etc. has great significance.Especially in western China, surface relief is larger, and near surface exists below a large amount of high steep
Structure realm, is the static correction related near surface and velocity modeling etc. to the greatest problem that high-dip structure region is imaged
Problem.
Seismic exploration technique of the prior art mainly includes level ground seismic exploration technique and vertical seismic profiling (VSP)
(VSP) seismic exploration technique.In the seismic prospecting of level ground, often focus and cymoscope are placed in into ground, echo
Need just be detected device reception through near-surface reflection twice, so, the reflected energy loss for receiving is larger.Therefore,
Although level ground seismic prospecting can effectively Underground horizontal media information, the detectivity to inclined reflection layer
Decline, therefore be difficult to obtain the imaging in high-precision high-dip structure region.In VSP seismic prospectings, by focus and cymoscope point
From, focus is located at earth's surface, and cymoscope is placed in well, and seismic signal can just be detected device and receive through a near surface, because
This has higher surveying accuracy.However, due to being constrained by observation system, the areas imaging of traditional VSP is not only limited, and
Also it is faced with static correction and the problems such as near-surface velocity is modeled, it is equally not high to the imaging precision of high-dip structure.
The content of the invention
The purpose of the application is to provide a kind for the treatment of method and apparatus of seismic imaging, can be reduced to static correction and near
The sensitivity of earth's surface rate pattern, improves the imaging precision of high-dip structure.
To achieve these goals, the invention provides a kind of processing method and processing device of seismic imaging, methods described and
What device was specifically realized in：
A kind of processing method of seismic imaging, methods described includes：
Step 1：Benchmark cymoscope is set in the excitation well at distance objective body predeterminable range；
Step 2：Obtain the reflection wave signal that first detector in the excitation well is received, the reflection wave signal includes the
One epicenter excitation and reflex to the signal of the first detector from the objective body；
Step 3：The direct-path signal by first epicenter excitation that the benchmark cymoscope is received is obtained, and will be described
Direct-path signal and the reflection wave signal correlation computations, generate benchmark cymoscope letter related to the first detector
Number；
Step 4：Obtain reflection wave signal of the first detector from other epicenter excitations one by one, repeat step 2 and 3,
The coherent signal being calculated between the direct wave of the benchmark cymoscope reception and the reflection wave signal, to multiple phases
OFF signal is sued for peace, and generates the reconstruction signal of the benchmark cymoscope.
Optionally, it is described by the direct-path signal and the reflection wave signal phase in one embodiment of the application
Close and calculate, generate the coherent signal of the benchmark cymoscope and the first detector, including：
The wave field function of the direct-path signal, the reflection wave signal in frequency domain is obtained respectively；
Calculate the reflection wave signal frequency domain wave field function and the direct-path signal frequency domain wave field function
Conjugate function product, using the real part of the product as the benchmark cymoscope it is related to the first detector letter
Number.
Optionally, in one embodiment of the application, the wave field function includes Green's function.
Optionally, in one embodiment of the application, the objective body includes tomography of the gradient more than first threshold
Tectosome.
Optionally, in one embodiment of the application, at least one focus is placed on predetermined position, described default
Position causes the direct-path signal to partially overlap with the reflection wave signal.
A kind of processing meanss of seismic imaging, described device includes：
Benchmark cymoscope setting unit, for arranging benchmark detection in the excitation well at distance objective body predeterminable range
Device；
Echo receiving unit, for obtaining the excitation well in the reflection wave signal that receives of first detector, it is described anti-
Ejected wave signal include the first epicenter excitation and reflex to the signal of the first detector from the objective body；
Correlation calculation unit, believes for obtaining the direct wave by first epicenter excitation that the benchmark cymoscope is received
Number, and by the direct-path signal and the reflection wave signal correlation computations, the benchmark cymoscope is generated with the described first inspection
The coherent signal of ripple device；
Reconfiguration unit, for obtaining reflection wave signal of the first detector from other epicenter excitations one by one, calculates
Coherent signal between the direct wave received to the benchmark cymoscope and the reflection wave signal, to multiple coherent signals
Sued for peace, generated the reconstruction signal of the benchmark cymoscope.
Optionally, in one embodiment of the application, the correlation calculation unit includes：
Wave field function acquiring unit, for obtaining the direct-path signal, the reflection wave signal respectively in frequency domain
Wave field function；
Coherent signal acquiring unit, believes in frequency domain wave field function for calculating the reflection wave signal with the direct wave
Number frequency domain wave field function conjugate function product, using the real part of the product as the benchmark cymoscope with it is described
The coherent signal of first detector.
Optionally, in one embodiment of the application, the wave field function includes Green's function.
Optionally, in one embodiment of the application, the objective body includes tomography of the gradient more than first threshold
Tectosome.
Optionally, in one embodiment of the application, at least one focus is placed on predetermined position, described default
Position causes the direct-path signal to partially overlap with the reflection wave signal.
The application provide seismic imaging processing method and processing device, can obtain benchmark cymoscope acquisition direct wave with
The echo that first detector is obtained carries out correlation computations, builds the wave field with cymoscope as image source.So as to VSP wave fields be turned
Turn to SWP wave fields so that VSP observation systems are converted into SWP observation systems.The reconstructing method of virtual focus is avoided on objective body
Complicated upper covering structure, therefore the rate pattern information between cymoscope and earth's surface focus need not be obtained, on avoiding
The deformation effect of complex dielectricss is covered, and it is insensitive to the static correction value in well.In addition, the SWP observation systems after reconstruct are more
Plus close to objective body, therefore obtain the objective body imaging of higher resolution.
Description of the drawings
In order to be illustrated more clearly that the embodiment of the present application or technical scheme of the prior art, below will be to embodiment or existing
The accompanying drawing to be used needed for having technology description is briefly described, it should be apparent that, drawings in the following description are only this
Some embodiments described in application, for those of ordinary skill in the art, in the premise for not paying creative labor
Under, can be with according to these other accompanying drawings of accompanying drawings acquisition.
Fig. 1 is a kind of method flow diagram of embodiment of the seismic imaging processing method that the application is provided；
Fig. 2 is the cymoscope energy acceptance graph of a relation that the application is provided；
Fig. 3 is a kind of method flow diagram of embodiment of the calculating coherent signal method that the application is provided；
Fig. 4 is the schematic diagram that the focus that the application is provided is steady phase point；
Fig. 5 is the walkaway VSP model that the application is provided；
Fig. 6 (a) is the SWP single shot records obtained by seismic interference process that the application is provided；
Fig. 6 (b) is the single shot record of the true observation system record that the application is provided；
Fig. 7 (a) -7 (d) is four kinds of different near-surface models that the application is provided；
Fig. 8 be the application provide carry out the list after SWP conversions for four kinds in Fig. 7 (a) -7 (d) different near-surface models
Big gun is recorded；
Fig. 9 is a kind of modular structure schematic diagram of embodiment of the seismic imaging processing meanss that the application is provided.
Specific embodiment
In order that those skilled in the art more fully understand the technical scheme in the application, below in conjunction with the application reality
The accompanying drawing in example is applied, the technical scheme in the embodiment of the present application is clearly and completely described, it is clear that described enforcement
Example is only some embodiments of the present application, rather than the embodiment of whole.Based on the embodiment in the application, this area is common
The every other embodiment that technical staff is obtained under the premise of creative work is not made, should all belong to the application protection
Scope.
Fig. 1 is a kind of method flow diagram of embodiment of the seismic imaging processing method that the application is provided, as shown in figure 1,
Methods described can include：
Step 1：Benchmark cymoscope is set in the excitation well at distance objective body predeterminable range.
As described above, in seismic prospecting, high-dip structure body is often difficult to be imaged, the objective body in the present embodiment can
With including the high-dip structure body, for example, the tectosome such as the high-dip structure body tomography, salt dome.At one of the application
In embodiment, the objective body can include fault tectonic body of the gradient more than first threshold, in the present embodiment method, work as institute
When the gradient for stating objective body is less than predetermined threshold value, the cymoscope in excitation well may not received from the objective body and reflected
The echo for coming.When the gradient of the objective body is bigger, it is imaged more accurate (being specifically described below).For example, described
One threshold value can be not less than 60 degree.
In the present embodiment, can excitation well be set at the objective body predeterminable range, and in the excitation well
Benchmark cymoscope is set.The benchmark cymoscope is used to record the wave field at the benchmark cymoscope position, subsequently may be used
Coherent calculation is carried out with the wave field that the wave field for recording in the benchmark cymoscope is obtained with other cymoscopes in the excitation well, it is raw
Into the benchmark cymoscope and the coherent signal of other cymoscopes.
Step 2：Obtain the reflection wave signal that first detector in the excitation well is received, the reflection wave signal includes the
One epicenter excitation and reflex to the signal of the first detector from the objective body.
In the present embodiment, the reflection wave signal that first detector is received in the excitation well, the echo can be obtained
Signal can include the first epicenter excitation and reflex to the signal of the first detector from the objective body.First inspection
Ripple device is the cymoscope of first detector predetermined position described in distance in the excitation well, or appoints in the excitation well
Cymoscope at meaning position, the application here is not limited.In the present embodiment, focus is positioned over the earth's surface around the excitation well
On, the energy that the epicenter excitation goes out can be reflected away by the objective body, and now, the first detector can receive institute
State the energy of objective body reflection.The reflection wave signal that the first detector is received is obtained by the first detector, equally
Ground, the echo is used to calculate the coherent signal for generating the benchmark cymoscope and the first detector.
Step 3：The direct-path signal by first epicenter excitation that the benchmark cymoscope is received is obtained, and will be described
Direct-path signal and the reflection wave signal correlation computations, generate benchmark cymoscope letter related to the first detector
Number.
In the present embodiment, the direct wave letter by first epicenter excitation that the benchmark cymoscope is received can be obtained
Number, and by the direct-path signal and the reflection wave signal correlation computations, the benchmark cymoscope is generated with the described first inspection
The coherent signal of ripple device.The first focus in the present embodiment can be the focus of predetermined position in earth's surface, or earth's surface
The focus of upper any position, the application here is not limited.As described above, the energy that produces of first epicenter excitation can be with
Carried by the target and being received after reflection and by the first detector, certainly, the energy that first epicenter excitation is produced is same
Sample can be received by the benchmark cymoscope.Fig. 2 is the cymoscope energy acceptance graph of a relation that the application is provided, as shown in Fig. 2 base
Quasi- cymoscope is direct-path signal from the signal that first focus is received, and the signal that the first detector is received is echo
Signal, and the reflection wave signal be first epicenter excitation and reflex to the first detector from the objective body
Signal.
After the direct wave and the echo is obtained, the direct wave can be calculated related to the echo
Signal.Specifically, Fig. 3 is a kind of method flow diagram of embodiment of the calculating coherent signal method that the application is provided, such as
It is described by the direct-path signal and the reflection wave signal correlation computations shown in Fig. 3, generate the benchmark cymoscope with it is described
The coherent signal of first detector, including：
S31：The wave field function of the direct-path signal, the reflection wave signal in frequency domain is obtained respectively.
S32：Calculate the reflection wave signal frequency domain wave field function and the direct-path signal frequency domain wave field
The product of the conjugate function of function is related to the first detector as the benchmark cymoscope using the real part of the product
Signal.
In the present embodiment, the wave field letter of the direct-path signal, the reflection wave signal in frequency domain can be respectively calculated
Number.In one embodiment of the application, the wave field function can be Green's function, and the Green's function can be counted exactly
Calculation obtains the wave field around the benchmark cymoscope.Illustrate by Green's function of wave field function, benchmark cymoscope is received
Direct-path signal can be expressed as in frequency domain：
Wherein, cymoscope on the basis of A, S is the first focus, and ω is circular frequency, τSAIt is the first focus to the focus detection
The direct wave hourage of device, ν is medium velocity.
The reflection wave signal that first detector is received can be expressed as in frequency domain：
Wherein, B is first detector, and S is the first focus, and ω is circular frequency, τSRFor the through of the first focus to objective body
Ripple hourage, τRBFor the echo hourage of objective body to the first detector, ν is medium velocity.
Now, the direct-path signal frequency domain wave field function with the reflection wave signal in frequency domain wave field function
The product of conjugate function can be expressed as：
G(B,S,ω)refG*(A,S,ω)dir (3)
Step 4：Obtain reflection wave signal of the first detector from other epicenter excitations one by one, repeat step 2 and 3,
The coherent signal being calculated between the direct wave of the benchmark cymoscope reception and the reflection wave signal, to multiple phases
OFF signal is sued for peace, and generates the reconstruction signal of the benchmark cymoscope.
In the present embodiment, reflection wave signal of the first detector from other epicenter excitations can be one by one obtained, be repeated
Step 2 and 3, the coherent signal being calculated between the direct wave of the benchmark cymoscope reception and the reflection wave signal is right
Multiple coherent signals are sued for peace, and generate the reconstruction signal of the benchmark cymoscope.Based on above-mentioned formula (1)-(3), can
To calculate direct-path signal that other focus excite at the benchmark cymoscope one by one and excite at the first detector
Reflection wave signal coherent signal, the expression formula of the coherent signal can be：
Wherein, B is first detector, cymoscope on the basis of A, and S is all focus in earth's surface,Expression takes G (B, A, ω)
Real part, ω is circular frequency, and ρ is Media density, and ν is medium velocity, τSAFor the direct wave trip of focus to the benchmark cymoscope
Row time, τSRFor the direct wave hourage of focus to objective body, τRBFor the echo trip of objective body to the first detector
The row time.
In one embodiment of the application, steady facies analyses are carried out to the formula (4), Fig. 4 is the shake that the application is provided
Source is the schematic diagram of steady phase point, as shown in figure 4, first focus can cause penetrating for the echo that first detector receives
The ray path of the direct wave that thread path is received with the benchmark cymoscope partially overlaps, that is, meet following formula (5)：
τSR=τSA+τAR (5)
Wherein, τSRFor the direct wave hourage of the first focus to objective body, τSAIt is the first focus to benchmark cymoscope
Direct wave hourage, τAROn the basis of cymoscope to objective body hourage.
Based on formula (5), then formula (4) can be expressed as：
Formula (6) for individual well section (SWP) approximate expression, by formula (2), (4), (5) can be by vertical seismic
Section (VSP) wave field is converted into individual well section (SWP) wave field, and formula (4) is equivalent to the direct wave and the echo
Seismic interference process.SWP wave fields after conversion represent that using the benchmark cymoscope in excitation well as new focus first detector is made
For the seismic response that receptor is received.Now, the position of the first focus shown in Fig. 4 be steady phase point, the steady phase point
Neighborhood be steady phase region.But, it is not that each focus is all located at steady phase region, therefore, ray path differs and establishes a capital coincidence, and
Actual seismic wavelet has certain hour to postpone.But sum formula (4) is that all focus are sued for peace, by being overlapped mutually between wavelet
Interfere, the focus response positioned at steady phase region is strengthened, non-steady phase region focus response weakens, and finally can effectively reconstruct SWP
Wave field.
Based on above-mentioned analysis, in one embodiment of the application, at least one focus is placed on predetermined position, institute
Predeterminated position is stated so that the direct-path signal partially overlaps with the reflection wave signal.At this point it is possible to ensure focus scope
Inside there is steady phase region so that the focus response of steady phase region is strengthened, and can effectively build SWP wave fields.
The present embodiment method is illustrated below by a specific scene, walkaway VSP mould as shown in Figure 5 is set up
Type.The size of the model is 2000m × 2000m, and is respectively the horizontal constant speed layer of 1500m/s, 3000m/s by two speed
Composition.Vertical excitation well is located at horizontal range 1400m, and high steep objective body, the high steep target are provided with the right side of excitation well
The speed of body is 4000m/s.60 big gun focus are uniformly distributed in earth's surface, shot interval 20m, and cymoscope starts to put from the 400m of underground
Put, spacing is 16m, the burnt ripple device totally 101.By forward simulation, 60 big gun VSP records are simulated altogether.With shown in Fig. 1
Method flow, to VSP single shot records seismic interference process is carried out.Fig. 6 (a) is represented
Benchmark cymoscope (i.e. image source), the SWP single shot records obtained by seismic interference process, section direct wave is in the same direction with echo
Axle is clear, and Fig. 6 (b) represents the single shot record that the true observation system blown out at (1400m, 400m) place recorded.Contrast Fig. 6
A () and Fig. 6 (b) understand, direct wave and basically identical from the echo form of high steep objective body.However, Fig. 6 (b) exist it is bright
Aobvious upgoing wave, by analysis, the upgoing wave is from the interface of bottom level course.Although there is also upgoing wave in Fig. 6 (a),
But energy is weaker, it is clear that from the emphasis that the reflection wave signal of high steep objective body is the high steep objective body imaging, therefore, institute
It is a kind of advantage to state imaging of the weak upgoing wave to high steep objective body.To sum up, through simulation contrast, demonstrate at seismic interference
Reason realizes the feasibility that VSP wave fields are converted to SWP wave fields, realizes that each cymoscope is " image source " by processing, and is prevented effectively from
Destruction of the actual focus to well and subsurface structure in well.
Below by another scene explanation above-described embodiment method, Fig. 7 (a) -7 (b) be the application provide four kinds not
Same near-surface model, wherein, Fig. 7 (a) represents that focus is uniformly distributed in level ground, earth's surface uniform dielectric；Fig. 7 (b) is represented
Focus is located at earth's surface, but skewness, and longitudinal buried depth is inconsistent；Fig. 7 (c) represents that focus is uniformly distributed in ground, and underground is shallow
Layer has velocity variations layer；Fig. 7 (d) represents that focus is located at relief surface, and earth's surface velocity variations have low velocity layer.By simulating
To four VSP records, through seismic interference process, realize that SWP wave fields are converted, respectively obtain the four SWP single-shots note shown in Fig. 8
Record.Four SWP single shot records correspond to respectively four kinds of near surface situations and hypocenter distributing situation.By analysis, reconstruct it is straight
There is good matching up to ripple and echo form, waveform and Energy distribution are normal.In the case of different near surfaces, reconstruct it is straight
There are portion disturbances up to ripple, the disturbance is mainly caused by focus boundary effect.However, migration imaging focus on from
The echo at high steep interface, therefore in processing procedure, by direct wave excision.By can be with by the contrast of four single shot records
It was found that, no matter earth's surface is any situation, and the echo form and energy of the SWP single-shots of reconstruct is basically identical, shows VSP wave fields
Focus static correction problem can be automatically eliminated to the conversion of SWP wave fields.Because the focus and cymoscope of the SWP sections of reconstruct are all located at again
Under miscellaneous coating, it is not necessary to understand the speed condition of near surface, therefore can avoid automatically in well more than first cymoscope
Unknown velocity information, skips near-surface velocity modeling, the local structure rate pattern directly set up near well.It is right respectively
Four SWP sections of Fig. 7 (a) -7 (b) reconstruct carry out Kirchhoff depth migrations, obtain respectively as shown in front four width figure in Fig. 8
Four migrated sections.Last width figure of Fig. 8 represents the migration result of traditional VSP imaging methods, and the migration result can be with
The profile of level course is imaged well, but the interface of high steep objective body can not be imaged exactly.Before in Fig. 8
The imaging results of four migrated sections shown in four width figures show that in addition to tiny disturbance, imaging results are basically identical, clearly
The form profile of high steep objective body is imaged clearly, and four migrated section bottoms have a trunnion axis, just with
The reflection of last width figure level course is corresponding in Fig. 8.In sum, using the SWP wave fields of reconstruct, near-earth can effectively be avoided
Table velocity modeling process, can be imaged exactly to the steep objective body of height near excitation well.
The processing method of the seismic imaging that the application is provided, can obtain direct wave and first inspection of the acquisition of benchmark cymoscope
The echo that ripple device is obtained carries out correlation computations, builds the wave field with cymoscope as image source.So as to VSP wave fields be converted into
SWP wave fields so that VSP observation systems are converted into SWP observation systems.The reconstructing method of virtual focus avoids complicated on objective body
Upper covering structure, therefore the rate pattern information between cymoscope and earth's surface focus need not be obtained, so as to avoid upper covering
The deformation effect of complex dielectricss, and it is insensitive to the static correction value in well.In addition, the SWP observation systems more adjunction after reconstruct
Objective body is bordering on, therefore obtains the objective body imaging of higher resolution.
Another aspect of the present invention also provides a kind of processing meanss of seismic imaging, and Fig. 9 is that the embodiment of the present application provides earthquake
A kind of modular structure schematic diagram of embodiment of imaging processing device, with reference to accompanying drawing 9, the device 90 includes：
Benchmark cymoscope setting unit 91, for arranging benchmark detection in the excitation well at distance objective body predeterminable range
Device；
Echo receiving unit 92, for obtaining the excitation well in first detector receive reflection wave signal, it is described
Reflection wave signal include the first epicenter excitation and reflex to the signal of the first detector from the objective body；
Correlation calculation unit 93, for obtaining the direct wave by first epicenter excitation that the benchmark cymoscope is received
Signal, and by the direct-path signal and the reflection wave signal correlation computations, generate the benchmark cymoscope and described first
The coherent signal of cymoscope；
Reconfiguration unit 94, for obtaining reflection wave signal of the first detector from other epicenter excitations one by one, calculates
The coherent signal between the direct wave and the reflection wave signal of the benchmark cymoscope reception is obtained, to multiple related letters
Number sued for peace, generated the reconstruction signal of the benchmark cymoscope.
The processing meanss of the seismic imaging that the application is provided, can obtain direct wave and first inspection of the acquisition of benchmark cymoscope
The echo that ripple device is obtained carries out correlation computations, builds the wave field with cymoscope as image source.So as to VSP wave fields be converted into
SWP wave fields so that VSP observation systems are converted into SWP observation systems.The reconstructing method of virtual focus avoids complicated on objective body
Upper covering structure, therefore the rate pattern information between cymoscope and earth's surface focus need not be obtained, so as to avoid upper covering
The deformation effect of complex dielectricss, and it is insensitive to the static correction value in well.In addition, the SWP observation systems more adjunction after reconstruct
Objective body is bordering on, therefore obtains the objective body imaging of higher resolution.
Optionally, in one embodiment of the application, the correlation calculation unit 93 includes：
Wave field function acquiring unit, for obtaining the direct-path signal, the reflection wave signal respectively in frequency domain
Wave field function；
Coherent signal acquiring unit, believes in frequency domain wave field function for calculating the reflection wave signal with the direct wave
Number frequency domain wave field function conjugate function product, using the real part of the product as the benchmark cymoscope with it is described
The coherent signal of first detector.
Optionally, in one embodiment of the application, the wave field function includes Green's function.
Optionally, in one embodiment of the application, the objective body includes tomography of the gradient more than first threshold
Tectosome.
Optionally, in one embodiment of the application, at least one focus is placed on predetermined position, described default
Position causes the direct-path signal to partially overlap with the reflection wave signal.
Each embodiment in this specification is described by the way of progressive, what each embodiment was stressed be with
The difference of other embodiment, between each embodiment identical similar part mutually referring to.Especially for system reality
For applying example, because it is substantially similar to embodiment of the method, so description is fairly simple, related part is referring to embodiment of the method
Part explanation.
Foregoing description in the application involved by each embodiment is only the application in some embodiments in the application,
Embodiment amended slightly can also carry out the scheme of each embodiment of above-mentioned the application on the basis of some methods.Certainly,
Meeting other deformations without creativeness of the process method step described in the application the various embodiments described above, still can be real
Existing identical application, will not be described here.
Although this application provides the method operating procedure as described in embodiment or flow chart, based on routine or noinvasive
The work of the property made can include more or less operating procedures.The step of enumerating in embodiment order is only numerous steps
A kind of mode in execution sequence, does not represent unique execution sequence.When device or client production in practice is performed, can
So that according to embodiment, either method shown in the drawings order performs either executed in parallel (such as at parallel processor or multithreading
The environment of reason).
Device or module that above-described embodiment is illustrated, specifically can be realized by computer chip or entity, or by having
The product of certain function is realizing.For convenience of description, it is divided into various modules with function when describing apparatus above to describe respectively.
The function of each module can be realized in same or multiple softwares and/or hardware when the application is implemented.It is of course also possible to
The module for realizing certain function is combined into realization by multiple submodule or subelement.
Method described herein, device or module can with computer readable program code mode realize controller by
Any appropriate mode realizes, for example, controller can take such as microprocessor or processor and the storage can be by (micro-)
The computer-readable medium of the computer readable program code (such as software or firmware) of computing device, gate, switch, specially
With integrated circuit (Application Specific Integrated Circuit, ASIC), programmable logic controller (PLC) and embedding
Enter the form of microcontroller, the example of controller includes but is not limited to following microcontroller：ARC 625D、Atmel AT91SAM、
Microchip PIC18F26K20 and Silicone Labs C8051F320, Memory Controller is also implemented as depositing
A part for the control logic of reservoir.It is also known in the art that except with pure computer readable program code mode reality
Beyond existing controller, controller can be caused by the way that method and step is carried out into programming in logic completely with gate, switch, special
The form of integrated circuit, programmable logic controller (PLC) and embedded microcontroller etc. is realizing identical function.Therefore this controller
A kind of hardware component is considered, and the device for realizing various functions included to its inside can also be considered as hardware
Structure in part.Or or even, can by be considered as the device for realizing various functions both can be implementation method software
Module can be again the structure in hardware component.
Part of module in herein described device can be in the general of computer executable instructions
Described in context, such as program module.Usually, program module includes performing particular task or realizes specific abstract data class
Routine, program, object, component, data structure, class of type etc..The application can also be in a distributed computing environment put into practice,
In these distributed computing environment, by the remote processing devices connected by communication network task is performed.Distributed
In computing environment, program module may be located at including in the local and remote computer-readable storage medium including storage device.
As seen through the above description of the embodiments, those skilled in the art can be understood that the application can
Realize by the mode of software plus required hardware.Based on such understanding, the technical scheme of the application is substantially in other words
The part contributed to prior art can be embodied in the form of software product, it is also possible to by the enforcement of Data Migration
During embody.The computer software product can be stored in storage medium, such as ROM/RAM, magnetic disc, CD, bag
Include some instructions to use so that a computer equipment (can be that personal computer, mobile terminal, server, or network set
It is standby etc.) perform method described in some parts of each embodiment of the application or embodiment.
Each embodiment in this specification is described by the way of progressive, same or analogous portion between each embodiment
Divide mutually referring to what each embodiment was stressed is the difference with other embodiment.The whole of the application or
Person part can be used in numerous general or special purpose computing system environments or configuration.For example：Personal computer, server are calculated
Machine, handheld device or portable set, laptop device, mobile communication terminal, multicomputer system, based on microprocessor it is
System, programmable electronic equipment, network PC, minicomputer, mainframe computer, including the distribution of any of the above system or equipment
Formula computing environment etc..
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application have it is many deformation and
Change is without deviating from spirit herein, it is desirable to which appended claim includes these deformations and changes without deviating from the application's
Spirit.
Claims (10)
1. a kind of processing method of seismic imaging, it is characterised in that methods described includes：
Step 1：Benchmark cymoscope is set in the excitation well at distance objective body predeterminable range；
Step 2：The reflection wave signal that first detector is received in the excitation well is obtained, the reflection wave signal includes the first shake
Source excitation and reflex to the signal of the first detector from the objective body；
Step 3：The direct-path signal by first epicenter excitation that the benchmark cymoscope is received is obtained, and will be described through
Ripple signal and the reflection wave signal correlation computations, generate the coherent signal of the benchmark cymoscope and the first detector；
Step 4：Reflection wave signal of the first detector from other epicenter excitations is obtained one by one, and repeat step 2 and 3 is calculated
The coherent signal between the direct wave and the reflection wave signal of the benchmark cymoscope reception is obtained, to multiple related letters
Number sued for peace, generated the reconstruction signal of the benchmark cymoscope.
2. method according to claim 1, it is characterised in that described by the direct-path signal and the reflection wave signal
Correlation computations, generate the coherent signal of the benchmark cymoscope and the first detector, including：
The wave field function of the direct-path signal, the reflection wave signal in frequency domain is obtained respectively；
Calculate the reflection wave signal to be total in the wave field function of frequency domain in frequency domain wave field function and the direct-path signal
The product of yoke function, using the real part of the product as the benchmark cymoscope and the first detector coherent signal.
3. method according to claim 2, it is characterised in that the wave field function includes Green's function.
4. method according to claim 1, it is characterised in that the objective body includes that gradient is disconnected more than first threshold
Layer tectosome.
5. method according to claim 1, it is characterised in that at least one focus is placed on predetermined position, described
Predeterminated position causes the direct-path signal to partially overlap with the reflection wave signal.
6. a kind of processing meanss of seismic imaging, it is characterised in that described device includes：
Benchmark cymoscope setting unit, for arranging benchmark cymoscope in the excitation well at distance objective body predeterminable range；
Echo receiving unit, for obtaining the excitation well in first detector receive reflection wave signal, the echo
Signal include the first epicenter excitation and reflex to the signal of the first detector from the objective body；
Correlation calculation unit, for obtaining the direct-path signal by first epicenter excitation that the benchmark cymoscope is received,
And by the direct-path signal and the reflection wave signal correlation computations, generate the benchmark cymoscope and the first detector
Coherent signal；
Reconfiguration unit, for obtaining reflection wave signal of the first detector from other epicenter excitations one by one, is calculated institute
The coherent signal between the direct wave and the reflection wave signal of benchmark cymoscope reception is stated, multiple coherent signals are carried out
Summation, generates the reconstruction signal of the benchmark cymoscope.
7. device according to claim 6, it is characterised in that the correlation calculation unit includes：
Wave field function acquiring unit, for obtaining the wave field of the direct-path signal, the reflection wave signal in frequency domain respectively
Function；
Coherent signal acquiring unit, exists in frequency domain wave field function for calculating the reflection wave signal with the direct-path signal
The product of the conjugate function of the wave field function of frequency domain, using the real part of the product as the benchmark cymoscope and described first
The coherent signal of cymoscope.
8. device according to claim 7, it is characterised in that the wave field function includes Green's function.
9. device according to claim 6, it is characterised in that the objective body includes that gradient is disconnected more than first threshold
Layer tectosome.
10. device according to claim 6, it is characterised in that at least one focus is placed on predetermined position, described
Predeterminated position causes the direct-path signal to partially overlap with the reflection wave signal.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611116919.2A CN106597536A (en) | 2016-12-07 | 2016-12-07 | Processing method and device of seismic imaging |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611116919.2A CN106597536A (en) | 2016-12-07 | 2016-12-07 | Processing method and device of seismic imaging |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106597536A true CN106597536A (en) | 2017-04-26 |
Family
ID=58596525
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611116919.2A Pending CN106597536A (en) | 2016-12-07 | 2016-12-07 | Processing method and device of seismic imaging |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106597536A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110058299A (en) * | 2018-09-14 | 2019-07-26 | 南方科技大学 | Method of earthquake locating, device and terminal device |
| CN111352151A (en) * | 2020-03-05 | 2020-06-30 | 广州地铁设计研究院股份有限公司 | Well source-deficiency method |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100118653A1 (en) * | 2008-11-08 | 2010-05-13 | Ruiqing He | Vertical seismic profiling velocity estimation method |
| CN104570073A (en) * | 2013-10-16 | 2015-04-29 | 中国石油化工股份有限公司 | Bi-reflection seismic wave imaging method applicable to complex, high and steep structure |
| CN105093292A (en) * | 2014-05-14 | 2015-11-25 | 中国石油天然气股份有限公司 | Data processing method and device for earthquake imaging |
-
2016
- 2016-12-07 CN CN201611116919.2A patent/CN106597536A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100118653A1 (en) * | 2008-11-08 | 2010-05-13 | Ruiqing He | Vertical seismic profiling velocity estimation method |
| CN104570073A (en) * | 2013-10-16 | 2015-04-29 | 中国石油化工股份有限公司 | Bi-reflection seismic wave imaging method applicable to complex, high and steep structure |
| CN105093292A (en) * | 2014-05-14 | 2015-11-25 | 中国石油天然气股份有限公司 | Data processing method and device for earthquake imaging |
Non-Patent Citations (5)
| Title |
|---|
| G. T. SCHUSTER ET AL.: "Interferometric/daylight seismic imaging", 《GEOPHYS.J.INT.》 * |
| 杨文采: "《反射地震学理论纲要》", 29 February 2012 * |
| 梁上林等: "一种基于地震干涉测量的高陡构造成像技术", 《中国地球科学联合学术年会2016》 * |
| 梁上林等: "不同参数对地震干涉测量实现VSP向SWP 数据重构的影响", 《SPG/SEG 北京2016国际地球物理会议》 * |
| 许卓等: "虚拟震源方法用于探测地下陡倾角地质构造", 《世界地质》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110058299A (en) * | 2018-09-14 | 2019-07-26 | 南方科技大学 | Method of earthquake locating, device and terminal device |
| CN111352151A (en) * | 2020-03-05 | 2020-06-30 | 广州地铁设计研究院股份有限公司 | Well source-deficiency method |
| CN111352151B (en) * | 2020-03-05 | 2022-06-03 | 广州地铁设计研究院股份有限公司 | Well source-deficiency method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Decker et al. | Diffraction imaging and time-migration velocity analysis using oriented velocity continuation | |
| Wittkamp et al. | Individual and joint 2-D elastic full-waveform inversion of Rayleigh and Love waves | |
| CN107407736B (en) | Generate the multistage full wave field inversion processing of the data set without multiple wave | |
| CN102841376A (en) | Retrieval method for chromatography speed based on undulating surface | |
| CN107817523B (en) | The analysis method and device of diffracted wave migration velocity | |
| Anquez et al. | Automatic correction and simplification of geological maps and cross-sections for numerical simulations | |
| CN106772593B (en) | The imaging method and device of diffracted wave | |
| Cosma et al. | 3D VSP migration by image point transform | |
| Li et al. | Finite difference seismic forward modeling method for fluid–solid coupled media with irregular seabed interface | |
| CN105137479B (en) | A kind of computational methods and device of bin degree of covering | |
| US20220221604A1 (en) | Method of analysing seismic data | |
| Datta et al. | Full-waveform inversion of salt models using shape optimization and simulated annealing | |
| Yang et al. | Using image warping for time-lapse image domain wavefield tomography | |
| Xue et al. | Predictive painting across faults | |
| CN105467445A (en) | Method and device for establishing three-dimensional seismic observation system | |
| CN105353406B (en) | A kind of method and apparatus for generating angle gathers | |
| CN106597536A (en) | Processing method and device of seismic imaging | |
| EP3612863B1 (en) | Post-stack kirchhoff depth de-migration method for tilted transverse isotropic (tti) and heterogeneous media based on ray tracing on migrated data | |
| CN110118993A (en) | Diffraction wave imaging method and device | |
| CN108594300B (en) | Landform imaging method and device and computer storage medium | |
| CN105572735B (en) | A kind of method and device for improving pre-stack depth migration imaging precision | |
| Wang et al. | Velocity model estimation of karstic fault reservoirs using full-waveform inversion accelerated on graphics processing unit | |
| CN104267434A (en) | Three-dimensional multi-component earthquake observation system geophone offset distribution obtaining method and device | |
| CN106772616A (en) | A kind of processing method and processing device of seismic imaging | |
| US20180059275A1 (en) | System and method for mapping horizons in seismic images |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170426 |