CN107561587A - A kind of method and device of phasing - Google Patents

A kind of method and device of phasing Download PDF

Info

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
Application number
CN201710827572.0A
Other languages
Chinese (zh)
Inventor
王万里
杨午阳
魏新建
禄娟
何欣
陈德武
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Natural Gas Co Ltd
Original Assignee
China Petroleum and Natural Gas Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Petroleum and Natural Gas Co Ltd filed Critical China Petroleum and Natural Gas Co Ltd
Priority to CN201710827572.0A priority Critical patent/CN107561587A/en
Publication of CN107561587A publication Critical patent/CN107561587A/en
Pending legal-status Critical Current

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

A kind of method and device of phasing
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)

  1. 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. 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. 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. 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. 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. 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 θs
    I'(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. 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>&amp;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>&amp;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>&amp;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. 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. 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. 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.
CN201710827572.0A 2017-09-14 2017-09-14 A kind of method and device of phasing Pending CN107561587A (en)

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)

* Cited by examiner, † Cited by third party
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

Patent Citations (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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