CN107561587A - A kind of method and device of phasing - Google Patents
A kind of method and device of phasing Download PDFInfo
- Publication number
- CN107561587A CN107561587A CN201710827572.0A CN201710827572A CN107561587A CN 107561587 A CN107561587 A CN 107561587A CN 201710827572 A CN201710827572 A CN 201710827572A CN 107561587 A CN107561587 A CN 107561587A
- Authority
- CN
- China
- Prior art keywords
- seismic channel
- channel data
- data
- specified
- phasing
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The embodiment of the present application discloses a kind of method and device of phasing.Methods described is provided with common midpoint gather data, and methods described includes:Model seismic channel data is determined from the common midpoint gather data;Based on the model seismic channel data, to specifying seismic channel data to carry out phasescan processing in the common midpoint gather data, target phase shift angle corresponding to the specified seismic channel data is obtained;Based on the target phase shift angle, phasing processing is carried out to the specified seismic channel data, obtains the specified seismic channel data after phasing.The technical scheme that the embodiment of the present application provides, can improve the degree of accuracy to the phasing in geological data.
Description
Technical field
The application is related to oil exploration Seismic Data Processing Technique field, the method and dress of more particularly to a kind of phasing
Put.
Background technology
As seismic exploration technique is to the development that more becomes more meticulous, requirement more and more higher of the high-resolution seismic survey to precision,
Improving the signal to noise ratio of seismic data and resolution ratio turns into the key issue of seism processing.In seism processing, generally
The signal to noise ratio of seismic data is improved using overlap-add procedure, but the assumed condition being superimposed is that sub-wave phase is identical, if wavelet
Phase is different, and overlap-add procedure incurs loss the energy for making back wave to some extent.Therefore, in seism processing, phase
Position difference problem is to influence one of an important factor for reflected signal realizes in-phase stacking.Phase correction techniques can eliminate the ripple of wavelet
Deformation, reach the purpose for the quality for improving stacked section.
The method for correcting phase of generally use is constant phase correction method in the prior art, and this method is mainly based upon son
The assumed condition that the excess phase of ripple changes independent of frequency, one is carried out independent of the constant of frequency to the phase spectrum of record
Adjustment.However, in actual seismic exploration process, excess phase is differentiated, not constant in each frequency range.Therefore,
The degree of accuracy using phase shift angle determined by the method for correcting phase of prior art is relatively low, can not completely effectively eliminate surplus
The influence of remaining phase.
The content of the invention
The purpose of the embodiment of the present application is to provide a kind of method and device of phasing, to improve in geological data
The degree of accuracy of phasing.
It is realized in order to solve the above technical problems, the embodiment of the present application provides a kind of method and device of phasing
's:
A kind of method of phasing, there is provided have common midpoint gather data, methods described includes:
Model seismic channel data is determined from the common midpoint gather data;
Based on the model seismic channel data, to specifying seismic channel data to carry out phase in the common midpoint gather data
Scan process, obtain target phase shift angle corresponding to the specified seismic channel data;
Based on the target phase shift angle, phasing processing is carried out to the specified seismic channel data, obtains phase school
Specified seismic channel data after just.
In preferred scheme, the determination model seismic channel data from the common midpoint gather data, including:
Processing is overlapped to the seismic channel data in the common midpoint gather data, obtains initially being superimposed seismic channel number
According to;
Determine that the first seismic channel data in the common midpoint gather data is initially superimposed seismic channel data with described
First coefficient correlation;
When first coefficient correlation is less than preset correlation coefficient number threshold value, to removing institute in the common midpoint gather data
State the seismic channel data beyond the first seismic channel data and be overlapped processing, obtain target superposition seismic channel data;
The target is superimposed seismic channel data as the model seismic channel data.
In preferred scheme, the value of the preset correlation coefficient number threshold value is 0.2.
In preferred scheme, the determination model seismic channel data from the common midpoint gather data, in addition to:
When the signal to noise ratio of the second seismic channel data in the midpoint trace gather data altogether is less than default snr threshold, to institute
State the seismic channel data in common midpoint gather data in addition to second seismic channel data and be overlapped processing, obtain target
It is superimposed seismic channel data;
The target is superimposed seismic channel data as the model seismic channel data.
It is described to be based on the model seismic channel data in preferred scheme, to specifying ground in the common midpoint gather data
Shake track data and carry out phasescan processing, obtain target phase shift angle corresponding to the specified seismic channel data, including:
Fourier transformation processing is carried out to the specified seismic channel data, frequency domain is obtained and specifies seismic channel data;
Based on phase shift angle is specified in preset phase shift angular range, seismic channel data is specified to carry out phase shift the frequency domain
Processing, obtain dephased frequency domain and specify seismic channel data;
Specify seismic channel data to carry out inversefouriertransform processing the dephased frequency domain, obtain dephased finger
Determine seismic channel data;
Determine to specify envelope corresponding to the dephased specified seismic channel data;
Determine the specified coefficient correlation of the specified envelope and the model seismic channel data;Wherein, the specified correlation
Coefficient corresponds with the specific phase shift angle;
Using the maximum phase shift angle of coefficient correlation in the preset phase shift angular range as the specified seismic channel data
Corresponding target phase shift angle.
In preferred scheme, seismic channel data is specified to carry out phase-shift processing the frequency domain using following formula:
R'(f)=R (f) cos θ-I (f) sin θs
I'(f)=R (f) sin θ+I (f) cos θ
Wherein, R'(f) and I'(f) represent that the dephased frequency domain specifies the real part and void of seismic channel data respectively
Portion;R (f) and I (f) represents that the frequency domain specifies the real and imaginary parts of seismic channel data respectively;θ represents the specific phase shift angle
Degree.
In preferred scheme, determine that the specified envelope is specified related to the model seismic channel data using following formula
Coefficient:
Wherein, corr represents the specified phase relation;xiRepresent to shake corresponding to ith sample point in the specified envelope
Width, AiRepresent amplitude corresponding to ith sample point in the model seismic channel data;N represents the specified envelope or the mould
The quantity of sampled point in type seismic channel data.
It is described to be based on the target phase shift angle in preferred scheme, phasing is carried out to the specified seismic channel data
Processing, obtains the specified seismic channel data after phasing, including:
Fourier transformation processing is carried out to the specified seismic channel data, frequency domain is obtained and specifies seismic channel data;
Based on the target phase shift angle, specify seismic channel data to carry out phase-shift processing the frequency domain, obtain phase shift
Specify seismic channel data in target frequency domain afterwards;
Specify seismic channel data to carry out inversefouriertransform processing in the dephased target frequency domain, obtain the phase
Specified seismic channel data after bit correction.
In preferred scheme, methods described also includes:
It is described from the common midpoint gather data determine model seismic channel data before, to the CMP road
Collection data enter correction process of taking action, and obtain the common midpoint gather data after dynamic correction;
Correspondingly, the determination model seismic channel data from the common midpoint gather data, including:From the dynamic school
Model seismic channel data is determined after just in common midpoint gather data.
A kind of device of phasing, described device provide common midpoint gather data, and described device includes:Model trace is true
Cover half block, target phase shift angle determining module and phase correction module;Wherein,
The model trace determining module, for determining model seismic channel data from the common midpoint gather data;
The target phase shift angle determining module, for based on the model seismic channel data, to the CMP road
Collecting in data specifies seismic channel data to carry out phasescan processing, obtains target angle of phase displacement corresponding to the specified seismic channel data
Degree;
The phase correction module, for based on the target phase shift angle, phase to be carried out to the specified seismic channel data
Bit correction processing, obtains the specified seismic channel data after phasing.
The embodiment of the present application provides a kind of method and device of phasing, can be from the common midpoint gather data
Middle determination model seismic channel data;The model seismic channel data can be based on, to being specified in the common midpoint gather data
Seismic channel data carries out phasescan processing, obtains target phase shift angle corresponding to the specified seismic channel data;It can be based on
The target phase shift angle, phasing processing is carried out to the specified seismic channel data, obtain after phasing specifiedly
Shake track data.In this way, the phase shift angle of each seismic channel data in the common midpoint gather data can be accurately determined, from
And improve the degree of accuracy of phasing in geological data.
Brief description of the drawings
, below will be to embodiment or existing in order to illustrate more clearly of the embodiment of the present application or technical scheme of the prior art
There is the required accompanying drawing used in technology description to be 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, do not paying the premise of creative labor
Under, other accompanying drawings can also be obtained according to these accompanying drawings.
Fig. 1 is a kind of flow chart of the embodiment of the method for phasing of the application;
Fig. 2 is the diagrammatic cross-section of the common midpoint gather data before phasing in the embodiment of the present application;
Fig. 3 is the diagrammatic cross-section of the common midpoint gather data after phasing in the embodiment of the present application;
Fig. 4 is the composition structure chart of the device embodiment of the application phasing.
Embodiment
The embodiment of the present application provides a kind of method and device of phasing.
In order that those skilled in the art more fully understand the technical scheme in the application, it is real below in conjunction with the application
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 implementation
Example only some embodiments of the present application, rather than whole embodiments.It is common based on the embodiment in the application, this area
The every other embodiment that technical staff is obtained under the premise of creative work is not made, it should all belong to the application protection
Scope.
The embodiment of the present application provides a kind of method of phasing.The method of the phasing is provided with CMP
Trace gather data.
In the present embodiment, can by way of seismic prospecting, gather specified location in purpose work area altogether in
Heart point trace gather data.
Fig. 1 is a kind of flow chart of the embodiment of the method for phasing of the application.As shown in figure 1, the phasing
Method, comprise the following steps.
Step S101:Model seismic channel data is determined from the common midpoint gather data.
In the present embodiment, can be with before model seismic channel data is determined from the common midpoint gather data
Correction process of taking action is entered to the common midpoint gather data, obtains the common midpoint gather data after dynamic correction.Correspondingly, from
Model seismic channel data is determined in the common midpoint gather data, can be included:From common midpoint gather after the dynamic correction
Model seismic channel data is determined in data.
In the present embodiment, model seismic channel data is determined from the common midpoint gather data, can specifically be wrapped
Include, processing can be overlapped to the seismic channel data in the common midpoint gather data, obtain initially being superimposed seismic channel number
According to.It can determine that the first seismic channel data in the common midpoint gather data initial is superimposed the of seismic channel data with described
One coefficient correlation.First seismic channel data can represent any seismic channel data in the common midpoint gather data.When
, can be to removing described first in the common midpoint gather data when first coefficient correlation is less than preset correlation coefficient number threshold value
Seismic channel data beyond seismic channel data is overlapped processing, obtains target superposition seismic channel data.Can be by the target
Seismic channel data is superimposed as the model seismic channel data.Wherein, the value of the preset correlation coefficient number threshold value can be
0.2。
In the present embodiment, model seismic channel data is determined from the common midpoint gather data, specifically can be with
Including, when it is described the signal to noise ratio of the second seismic channel data is less than default snr threshold in the trace gather data of midpoint altogether when, can be right
Seismic channel data in the common midpoint gather data in addition to second seismic channel data is overlapped processing, obtains mesh
Mark superposition seismic channel data.Wherein, the span of the default snr threshold can be 50~60 decibels.Can be by described in
Target is superimposed seismic channel data as the model seismic channel data.
Step S102:Based on the model seismic channel data, to specifying seismic channel number in the common midpoint gather data
According to phasescan processing is carried out, target phase shift angle corresponding to the specified seismic channel data is obtained.
In the present embodiment, based on the model seismic channel data, to specifying ground in the common midpoint gather data
Shake track data carry out phasescan processing, obtain the specified seismic channel data corresponding to target phase shift angle, can specifically wrap
Include, Fourier transformation processing can be carried out to the specified seismic channel data, obtain frequency domain and specify seismic channel data.Based on pre-
If specifying phase shift angle in phase shift angle scope, seismic channel data can be specified to carry out phase-shift processing the frequency domain, can be with
Obtain dephased frequency domain and specify seismic channel data.Seismic channel data can be specified to carry out the dephased frequency domain anti-
Fourier transformation processing, obtains dephased specified seismic channel data.The dephased specified seismic channel data can be determined
Corresponding specified envelope.The specified coefficient correlation of the specified envelope and the model seismic channel data can be determined.Wherein, institute
Specified coefficient correlation is stated to correspond with the specific phase shift angle.Can be by coefficient correlation in the preset phase shift angular range
Maximum phase shift angle is as target phase shift angle corresponding to the specified seismic channel data.
In the present embodiment, following formula can be used to specify seismic channel data to carry out at phase shift the frequency domain
Reason:
R'(f)=R (f) cos θ-I (f) sin θs
I'(f)=R (f) sin θ+I (f) cos θ
Wherein, R'(f) and I'(f) represent that the dephased frequency domain specifies the real part and void of seismic channel data respectively
Portion;R (f) and I (f) represents that the frequency domain specifies the real and imaginary parts of seismic channel data respectively;θ represents the specific phase shift angle
Degree.
In the present embodiment, Hilbert transform processing can be carried out to the dephased specified seismic channel data,
Determine to specify envelope corresponding to the dephased specified seismic channel data.
In the present embodiment, following formula can be used to determine the specified envelope and the model seismic channel data
Specify coefficient correlation:
Wherein, corr represents the specified phase relation;xiRepresent to shake corresponding to ith sample point in the specified envelope
Width, AiRepresent amplitude corresponding to ith sample point in the model seismic channel data;N represents the specified envelope or the mould
The quantity of sampled point in type seismic channel data.
Step S103:Based on the target phase shift angle, phasing processing is carried out to the specified seismic channel data, obtained
Specified seismic channel data after to phasing.
In the present embodiment, based on the target phase shift angle, phasing is carried out to the specified seismic channel data
Processing, obtains the specified seismic channel data after phasing, can specifically include, and the specified seismic channel data can be carried out
Fourier transformation processing, obtain frequency domain and specify seismic channel data., can be to the frequency domain based on the target phase shift angle
Specify seismic channel data to carry out phase-shift processing, obtain dephased target frequency domain and specify seismic channel data.Can be to the phase
Target frequency domain after shifting specifies seismic channel data to carry out inversefouriertransform processing, obtain after the phasing specifiedly
Shake track data.
For example, Fig. 2 is the diagrammatic cross-section of the common midpoint gather data before phasing in the embodiment of the present application.Fig. 3
It is the diagrammatic cross-section of the common midpoint gather data in the embodiment of the present application after phasing.Abscissa in Fig. 2 and Fig. 3 and
Ordinate represents road number and the sampling time of seismic channel respectively.As shown in Figures 2 and 3, using the present processes in Fig. 2
After seismic channel data in common midpoint gather data carries out phasing, the influence of excess phase can be effectively eliminated.
The embodiment of the method for the phasing, model seismic channel number can be determined from the common midpoint gather data
According to;The model seismic channel data can be based on, to specifying seismic channel data to carry out phase in the common midpoint gather data
Scan process, obtain target phase shift angle corresponding to the specified seismic channel data;The target phase shift angle can be based on, it is right
The specified seismic channel data carries out phasing processing, obtains the specified seismic channel data after phasing.In this way, can be accurate
The phase shift angle of each seismic channel data in the common midpoint gather data is determined, so as to improve phase in geological data
The degree of accuracy of correction.
Fig. 4 is the composition structure chart of the device embodiment of the application phasing.The device of the phasing provides altogether
Central point trace gather data.As shown in figure 4, the device of the phasing can include:Model trace determining module 100, target phase
Move angle-determining module 200 and phase correction module 300.
The model trace determining module 100, it can be used for determining model seismic channel from the common midpoint gather data
Data.
The target phase shift angle determining module 200, it can be used for being based on the model seismic channel data, in described be total to
Specify seismic channel data to carry out phasescan processing in heart point trace gather data, obtain target corresponding to the specified seismic channel data
Phase shift angle.
The phase correction module 300, it can be used for being based on the target phase shift angle, to the specified seismic channel data
Phasing processing is carried out, obtains the specified seismic channel data after phasing.
The device embodiment of the phasing is corresponding with the embodiment of the method for the phasing, it is possible to achieve phase
The technical scheme of the embodiment of the method for correction, and the technique effect of adquisitiones embodiment.
In the 1990s, the improvement for a technology can clearly distinguish be on hardware improvement (for example,
Improvement to circuit structures such as diode, transistor, switches) or software on improvement (improvement for method flow).So
And as the development of technology, the improvement of current many method flows can be considered as directly improving for hardware circuit.
Designer nearly all obtains corresponding hardware circuit by the way that improved method flow is programmed into hardware circuit.Cause
This, it cannot be said that the improvement of a method flow cannot be realized with hardware entities module.For example, PLD
(Programmable Logic Device, PLD) (such as field programmable gate array (Field Programmable Gate
Array, FPGA)) it is exactly such a integrated circuit, its logic function is determined by user to device programming.By designer
Voluntarily programming comes a digital display circuit " integrated " on a piece of PLD, without asking chip maker to design and make
Special IC chip.Moreover, nowadays, substitution manually makes IC chip, this programming is also used instead mostly " patrols
Volume compiler (logic compiler) " software realizes that software compiler used is similar when it writes with program development,
And the source code before compiling also write by handy specific programming language, this is referred to as hardware description language
(Hardware Description Language, HDL), and HDL is also not only a kind of, but have many kinds, such as ABEL
(Advanced Boolean Expression Language)、AHDL(Altera Hardware Description
Language)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL
(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(Ruby
Hardware Description Language) etc., VHDL (Very-High-Speed are most generally used at present
Integrated Circuit Hardware Description Language) and Verilog2.Those skilled in the art
It will be apparent to the skilled artisan that only need method flow slightly programming in logic and being programmed into integrated circuit with above-mentioned several hardware description languages
In, it is possible to it is readily available the hardware circuit for realizing the logical method flow.
It is also known in the art that in addition to realizing controller in a manner of pure computer readable program code, it is complete
Entirely can by by method and step carry out programming in logic come controller with gate, switch, application specific integrated circuit, may be programmed
The form of logic controller and embedded microcontroller etc. realizes identical function.Therefore this controller is considered one kind
Hardware component, and it is used to realize that the device of various functions can also to be considered as the structure in hardware component to what is included in it.Or
Even, it not only can be able to will be the software module of implementation method for realizing that the device of various functions is considered as but also can be Hardware Subdivision
Structure in part.
Device that above-described embodiment illustrates, module, it can specifically be realized by computer chip or entity, or by with certain
The product of kind of function is realized.
For convenience of description, it is divided into various modules during description apparatus above with function to describe respectively.Certainly, this is being implemented
The function of each module can be realized in same or multiple softwares and/or hardware during application.
As seen through the above description of the embodiments, those skilled in the art can be understood that the application can
Realized by the mode of software plus required general hardware platform.Based on such understanding, the technical scheme essence of the application
On the part that is contributed in other words to prior art can be embodied in the form of software product, in a typical configuration
In, computing device includes one or more processors (CPU), input/output interface, network interface and internal memory.The computer is soft
Part product can include some instructions make it that a computer equipment (can be personal computer, server, or network
Equipment etc.) perform method described in some parts of each embodiment of the application or embodiment.The computer software product can
To be stored in internal memory, internal memory may include the volatile memory in computer-readable medium, random access memory
(RAM) and/or the form such as Nonvolatile memory, such as read-only storage (ROM) or flash memory (flash RAM).Internal memory is computer
The example of computer-readable recording medium.Computer-readable medium includes permanent and non-permanent, removable and non-removable media can be by
Any method or technique come realize information store.Information can be computer-readable instruction, data structure, the module of program or its
His data.The example of the storage medium of computer includes, but are not limited to phase transition internal memory (PRAM), static RAM
(SRAM), dynamic random access memory (DRAM), other kinds of random access memory (RAM), read-only storage
(ROM), Electrically Erasable Read Only Memory (EEPROM), fast flash memory bank or other memory techniques, read-only optical disc are read-only
Memory (CD-ROM), digital versatile disc (DVD) or other optical storages, magnetic cassette tape, tape magnetic rigid disk storage or
Other magnetic storage apparatus or any other non-transmission medium, the information that can be accessed by a computing device available for storage.According to
Herein defines, and computer-readable medium does not include of short duration computer readable media (transitory media), such as modulation
Data-signal and carrier wave.
Each embodiment in this specification is described by the way of progressive, identical similar portion between each embodiment
Divide mutually referring to what each embodiment stressed is the difference with other embodiment.It is real especially for device
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.
The application can be used in numerous general or special purpose computing system environments or configuration.Such as:Personal computer, clothes
Business device computer, handheld device or portable set, laptop device, multicomputer system, the system based on microprocessor, put
Top box, programmable consumer-elcetronics devices, network PC, minicom, mainframe computer including any of the above system or equipment
DCE etc..
The application can be described in the general context of computer executable instructions, such as program
Module.Usually, program module includes performing particular task or realizes routine, program, object, the group of particular abstract data type
Part, data structure etc..The application can also be put into practice in a distributed computing environment, in these DCEs, by
Task is performed and connected remote processing devices by communication network.In a distributed computing environment, program module can be with
In the local and remote computer-readable storage medium including storage device.
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 departing from spirit herein, it is desirable to which appended claim includes these deformations and changed without departing from the application's
Spirit.
Claims (10)
- A kind of 1. method of phasing, it is characterised in that be provided with common midpoint gather data, methods described includes:Model seismic channel data is determined from the common midpoint gather data;Based on the model seismic channel data, to specifying seismic channel data to carry out phasescan in the common midpoint gather data Processing, obtains target phase shift angle corresponding to the specified seismic channel data;Based on the target phase shift angle, phasing processing is carried out to the specified seismic channel data, after obtaining phasing Specified seismic channel data.
- 2. the method for a kind of phasing according to claim 1, it is characterised in that described from the common midpoint gather Model seismic channel data is determined in data, including:Processing is overlapped to the seismic channel data in the common midpoint gather data, obtains initially being superimposed seismic channel data;Determine that the first seismic channel data in the common midpoint gather data is initially superimposed the first of seismic channel data with described Coefficient correlation;When first coefficient correlation is less than preset correlation coefficient number threshold value, in the common midpoint gather data except described the Seismic channel data beyond one seismic channel data is overlapped processing, obtains target superposition seismic channel data;The target is superimposed seismic channel data as the model seismic channel data.
- 3. the method for a kind of phasing according to claim 2, it is characterised in that the preset correlation coefficient number threshold value Value is 0.2.
- 4. the method for a kind of phasing according to claim 1, it is characterised in that described from the common midpoint gather Model seismic channel data is determined in data, in addition to:When the signal to noise ratio of the second seismic channel data in the midpoint trace gather data altogether is less than default snr threshold, to described common Seismic channel data in central point trace gather data in addition to second seismic channel data is overlapped processing, obtains target superposition Seismic channel data;The target is superimposed seismic channel data as the model seismic channel data.
- 5. the method for a kind of phasing according to claim 1, it is characterised in that described to be based on the model seismic channel Data, to specifying seismic channel data to carry out phasescan processing in the common midpoint gather data, obtain the specified earthquake Target phase shift angle corresponding to track data, including:Fourier transformation processing is carried out to the specified seismic channel data, frequency domain is obtained and specifies seismic channel data;Based on phase shift angle is specified in preset phase shift angular range, seismic channel data is specified to carry out at phase shift the frequency domain Reason, obtain dephased frequency domain and specify seismic channel data;To the dephased frequency domain specify seismic channel data carry out inversefouriertransform processing, obtain it is dephased specifiedly Shake track data;Determine to specify envelope corresponding to the dephased specified seismic channel data;Determine the specified coefficient correlation of the specified envelope and the model seismic channel data;Wherein, the specified coefficient correlation Corresponded with the specific phase shift angle;The maximum phase shift angle of coefficient correlation in the preset phase shift angular range is corresponding as the specified seismic channel data Target phase shift angle.
- 6. the method for a kind of phasing according to claim 5, it is characterised in that using following formula to the frequency Seismic channel data is specified to carry out phase-shift processing in domain:R'(f)=R (f) cos θ-I (f) sin θsI'(f)=R (f) sin θ+I (f) cos θWherein, R'(f) and I'(f) represent that the dephased frequency domain specifies the real and imaginary parts of seismic channel data respectively;R (f) and I (f) represents that the frequency domain specifies the real and imaginary parts of seismic channel data respectively;θ represents the specific phase shift angle.
- 7. the method for a kind of phasing according to claim 5, it is characterised in that determine the finger using following formula Determine the specified coefficient correlation of envelope and the model seismic channel data:<mrow> <mi>c</mi> <mi>o</mi> <mi>r</mi> <mi>r</mi> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>x</mi> <mi>i</mi> </msub> <msub> <mi>A</mi> <mi>i</mi> </msub> </mrow> <msqrt> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msubsup> <mi>x</mi> <mi>i</mi> <mn>2</mn> </msubsup> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msubsup> <mi>A</mi> <mi>i</mi> <mn>2</mn> </msubsup> </mrow> </msqrt> </mfrac> </mrow>Wherein, corr represents the specified phase relation;xiRepresent amplitude, A corresponding to ith sample point in the specified envelopeiTable Show amplitude corresponding to ith sample point in the model seismic channel data;N represents the specified envelope or the model seismic channel The quantity of sampled point in data.
- 8. the method for a kind of phasing according to claim 1, it is characterised in that described to be based on the target angle of phase displacement Degree, phasing processing is carried out to the specified seismic channel data, obtains the specified seismic channel data after phasing, including:Fourier transformation processing is carried out to the specified seismic channel data, frequency domain is obtained and specifies seismic channel data;Based on the target phase shift angle, specify seismic channel data to carry out phase-shift processing the frequency domain, obtain dephased Specify seismic channel data in target frequency domain;Specify seismic channel data to carry out inversefouriertransform processing in the dephased target frequency domain, obtain the phase school Specified seismic channel data after just.
- 9. the method for a kind of phasing according to claim 1, it is characterised in that methods described also includes:Before model seismic channel data is determined from the common midpoint gather data, the common midpoint gather data are entered Action correction process, obtain the common midpoint gather data after dynamic correction;Correspondingly, model seismic channel data is determined from the common midpoint gather data, including:In being total to after the dynamic correction Model seismic channel data is determined in heart point trace gather data.
- 10. a kind of device of phasing, it is characterised in that described device provides common midpoint gather data, described device bag Include:Model trace determining module, target phase shift angle determining module and phase correction module;Wherein,The model trace determining module, for determining model seismic channel data from the common midpoint gather data;The target phase shift angle determining module, for based on the model seismic channel data, to the common midpoint gather number Specify seismic channel data to carry out phasescan processing in, obtain target phase shift angle corresponding to the specified seismic channel data;The phase correction module, for based on the target phase shift angle, phase school to be carried out to the specified seismic channel data Positive processing, obtains the specified seismic channel data after phasing.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710827572.0A CN107561587A (en) | 2017-09-14 | 2017-09-14 | A kind of method and device of phasing |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710827572.0A CN107561587A (en) | 2017-09-14 | 2017-09-14 | A kind of method and device of phasing |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107561587A true CN107561587A (en) | 2018-01-09 |
Family
ID=60980905
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710827572.0A Pending CN107561587A (en) | 2017-09-14 | 2017-09-14 | A kind of method and device of phasing |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107561587A (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4348749A (en) * | 1975-06-30 | 1982-09-07 | Mobil Oil Corporation | Phase correcting seismic traces |
CN101201409A (en) * | 2006-12-14 | 2008-06-18 | 中国石油天然气集团公司 | Method for revising temblor data phase-conversion bit |
CN103217716A (en) * | 2012-01-18 | 2013-07-24 | 中国石油天然气集团公司 | Method of fast correction processing of phase of seismic data |
CN103645500A (en) * | 2013-11-08 | 2014-03-19 | 中国石油大学(北京) | Method for estimating mixed-phase seismic wavelets of frequency domain |
CN106842297A (en) * | 2015-12-04 | 2017-06-13 | 中国石油化工股份有限公司 | Borehole restraint unstable state method for correcting phase |
-
2017
- 2017-09-14 CN CN201710827572.0A patent/CN107561587A/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4348749A (en) * | 1975-06-30 | 1982-09-07 | Mobil Oil Corporation | Phase correcting seismic traces |
CN101201409A (en) * | 2006-12-14 | 2008-06-18 | 中国石油天然气集团公司 | Method for revising temblor data phase-conversion bit |
CN103217716A (en) * | 2012-01-18 | 2013-07-24 | 中国石油天然气集团公司 | Method of fast correction processing of phase of seismic data |
CN103645500A (en) * | 2013-11-08 | 2014-03-19 | 中国石油大学(北京) | Method for estimating mixed-phase seismic wavelets of frequency domain |
CN106842297A (en) * | 2015-12-04 | 2017-06-13 | 中国石油化工股份有限公司 | Borehole restraint unstable state method for correcting phase |
Non-Patent Citations (3)
Title |
---|
刘玉金 等: "基于局部相似度的叠前非稳态相位校正方法", 《石油地球物理勘探》 * |
张玉芬: "常相位校正与分频相位校正方法在储层预测中的应用", 《1996年中国地球物理学会第十二届学术年会》 * |
李正文: "《高分辨率地震勘探》", 30 September 1993, 成都科技大学出版社 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106249299B (en) | The lower weak reflected seismic energy restoration methods of thin layer of strong reflection shielding and device | |
CN106646611B (en) | A kind of method and device of longitudinal wave and converted wave seismic data time match | |
CN106772589B (en) | A kind of prestack seismic inversion method and device | |
CN106842304B (en) | A kind of prestack depth migration method and device | |
CN107765298A (en) | A kind of method and device for determining near-surface velocity model | |
CN107679358A (en) | A kind of method and device for determining reservoir permeability | |
CN107193045A (en) | A kind of seismic data processing technique and device | |
CN106353818A (en) | Seismic data time difference correcting method and system | |
CN106054251A (en) | First break picking method and device | |
CN108387926A (en) | A kind of method and device of determining air-gun array far-field wavelet | |
CN108181649A (en) | A kind of dielectric interface vibration amplitude compensation method and device | |
CN107092031A (en) | A kind of sound wave measuring well curve bearing calibration and device | |
CN108008444A (en) | A kind of method and device of definite petrophysical model | |
CN108181644A (en) | A kind of method of sampling and device | |
CN107179546B (en) | A kind of method and device of determining velocity field | |
CN106291691A (en) | A kind of seismic migration imaging method and device | |
CN108897036A (en) | A kind of seismic data processing technique and device | |
CN104459774B (en) | A kind of geologic lithology difference identification method and system | |
CN108035709A (en) | A kind of definite method and device of shale reservoir quality | |
CN106353816A (en) | Seismic acquisition footprint noise suppression method and system | |
CN106772599B (en) | A kind of method and device calculating formation shear speed | |
CN107561588A (en) | A kind of geological data noise drawing method and device | |
CN107561587A (en) | A kind of method and device of phasing | |
CN108107470B (en) | A kind of seismic data processing technique and device | |
CN108196300B (en) | A kind of seismic data processing technique and device |
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 |
Application publication date: 20180109 |
|
RJ01 | Rejection of invention patent application after publication |