CN107589455A - A kind of method and device of land and water detector seismic data merging treatment - Google Patents
A kind of method and device of land and water detector seismic data merging treatment Download PDFInfo
- Publication number
- CN107589455A CN107589455A CN201710565184.XA CN201710565184A CN107589455A CN 107589455 A CN107589455 A CN 107589455A CN 201710565184 A CN201710565184 A CN 201710565184A CN 107589455 A CN107589455 A CN 107589455A
- Authority
- CN
- China
- Prior art keywords
- mrow
- land
- seismic data
- water
- data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Abstract
The embodiment of the present application discloses a kind of method and device of land and water detector seismic data merging treatment.Methods described includes：Phase Processing is carried out to land detector seismic data；The amplitude processing is carried out to detector seismic data in water and to the land detector seismic data after Phase Processing；Fourier transformation is carried out respectively to detector seismic data in the water after the amplitude processing and land detector seismic data；Calibration factor is calculated according to detector seismic data in the water after the amplitude processing and land detector seismic data, and according to detector seismic data in the water after Fourier transformation and land detector seismic data, calculates sea water advanced and bottom reflection coefficient；Based on sea water advanced, mirror image connector is determined；Mirror image combined treatment is carried out to detector seismic data in the water of target area and land detector seismic data, obtains land and water wave detector merging data.The embodiment that the application provides, the degree of accuracy of land and water detector seismic data merging treatment can be improved.
Description
Technical field
The application is related to oil exploration, exploitation and development technique field, more particularly to a kind of land and water detector seismic data
The method and device of merging treatment.
Background technology
With the development of seismic exploration technique, the difficulty and depth of exploratory engineering of off-shore petroleum/gas reservoir are also increasing, to seismic data
Signal to noise ratio and resolution requirement also more and more higher.OBC (Ocean Bottom Cable, submarine cable), it is a kind of joint sea
Wave detector, is fixed on seabed, can obtain high-resolution three-dimension geological data by upper and land seismic data acquisition.
In OBC data acquisitions, at least three ships：One energy source boat, air gun source arrangement is only pulled, carry out seismic wave and excite；One
Taken-over vessel, is fixed, is connected to submarine cable, receives seismic wave；One ship or several ships, build a submarine cable
With recovery submarine cable.In OBC data acquisitions, because seabed and sea are all stronger reflecting interfaces.With epicenter excitation
A seismic wavelet reach seabed from hypocentral location, or a reflection seismic wavelet is issued to seabed, submarine cable from ground
In wave detector, sense and record this reflection seismic wavelet.This reflection wavelet continues up advance and reaches sea, by
The reflection on sea, then change direction and propagate downwards, reach seabed.Wave detector in submarine cable, senses and records again
This lower seismic wavelet.This seismic wavelet is reflected by seabed simultaneously, then changes direction and upwardly propagates, reaches sea,
Reflected by sea, then change direction and propagate downwards, reach seabed.This cyclic process repeats.And these are original
Reflection seismic wavelet is undesirable secondary and subsequently reaches, and is exactly the more subwaves of seawater singing (reverberation).Seawater singing is multiple
Ripple is noise jamming maximum in offshore seismic exploration data.The more subwave noise jammings of seawater singing are eliminated, are offshore earthquake numbers
According to step mostly important in processing.
The data of OBC collections, there is provided detector seismic data and land detector seismic data two in same position water
Kind data.Both data use wave detector and land wave detector in water to record respectively.Wave detector is a kind of pressure detection in water
Device, record is pressure change caused by seismic wave；Land wave detector is a kind of particle velocity wave detector, and record is particle speed
Degree change.Because the recording mechanism of both wave detectors is different, for the multiple wave interference of seawater singing at same position, show
Different characteristic.Compared with the multiple wave interference of seawater singing that wave detector in water records, the seawater singing of land wave detector record is more
Subwave interference shows polarity and amplitude characteristic difference.The multiple wave interference of seawater singing of two kinds of wave detector records, its polarity are
Opposite, amplitude is different, and differs a constant proportional to bottom reflection coefficient, this constant value be exactly demarcate because
Son.Therefore this amplitude and polar character difference are utilized, can effectively eliminate the multiple wave interference of seawater singing.So eliminate seawater
The step of singing multiple wave interference, includes：(1) in each receiving point opening position, detector seismic data and land in water are recorded
Two kinds of data of detector seismic data；(2) using the transducer sensitivity (transmission constant) of two kinds of wave detectors, land detection is adjusted
The amplitude of device geological data, to match the amplitude of detector seismic data in water；(3) calculate determine calibration factor, it is sea water advanced,
Bottom reflection coefficient parameter；(4) calibration factor is utilized, demarcation adjusts the land detector seismic data after amplitude；(5) demarcation
Land detector seismic data afterwards, it is added with detector seismic data in corresponding water, obtains water-land geophone calibration number
According to；(6) sea water advanced, bottom reflection coefficient parameter is utilized, eliminates the multiple wave interference of seawater singing.So, calculate calibration factor,
The parameters such as sea water advanced, bottom reflection coefficient, constitute indoor marine seismic data processing and eliminate the multiple wave interference of seawater singing
The basic method and committed step of method.
Conventional treatment method, calculated using scan method and determine calibration factor.Using presetting a calibration factor model
Value and scanning step are enclosed, a series of calibration factor value is provided using scan method, then calculates detector seismic data in water
With the data of land detector seismic data and, then to data and calculate auto-correlation function, calculated by auto-correlation function most generous
Differential mode, finally by the maximum variance modulus value of maximum, determine calibration factor numerical value.This method need substantial amounts of autocorrelation calculation and
Maximum variance module calculates, therefore calculates very time-consuming.Conventional treatment method, using upstream wave field and down-going wave fields data cross-correlation
It is sea water advanced that method calculates determination.Using detector seismic data in water and land detector seismic data, upstream wave field is calculated
With down-going wave fields data；Then upstream wave field and down-going wave fields data cross-correlation are calculated, by cross-correlation function maximum, is determined
Sea water advanced value.Due in detector seismic data and land detector seismic data in actual water, including various noises
Interference.Especially since land and water inspection data effective band scope is different, low-frequency noise (such as face ripple etc.) and useless high frequency division
Cloth is also different, particularly land inspection data, includes stronger surface wave interference.In addition, also comprising strong in prestack land and water inspection data
Amplitude energy disturbs.The cross-correlation function that detector seismic data and land detector seismic data calculate in water is so used,
Also various noise contributions are included.Therefore the cross-correlation function comprising noise is used, that determines is sea water advanced, exists very big
Error, it is difficult to meet real data processing requirement.Conventional treatment method, calculated using scan method and determine bottom reflection coefficient.
Using a bottom reflection coefficient value range and scanning step is preset, a series of reflectance factor is provided using scan method
Value, to a series of bottom reflection coefficients, calculates a series of more subwave inverse filter operators of Backus seawater singings, to frequency domain water
Land inspection upstream wave field data (that is, land and water detector seismic data sum), which are multiplied by the more subwave inverse filters of Backus seawater singings, to be calculated
Son, the upstream wave field data after filtering process are obtained, then using inverse Fourier transform, upstream wave field data are transformed into the time
Domain, a series of wavefield data energy, the wherein bottom reflection coefficient corresponding to least energy are calculated in time-domain, be exactly it is required most
Good bottom reflection coefficient.This method needs substantial amounts of autocorrelation calculation and maximum variance module to calculate, therefore calculates very time-consuming.
Therefore, how accurately, timesaving calculate and determine calibration factor, sea water advanced and bottom reflection coefficient, so as to it is quick,
Accurately realize that wave detector and the merging treatment of land detector seismic data are that seawater singing is more in elimination geological data in water
Subwave interference effect, improve geological data signal to noise ratio and the problem of resolution ratio urgent need to resolve.
The content of the invention
The purpose of the embodiment of the present application is to provide a kind of method and device of land and water detector seismic data merging treatment, with
Quickly and accurately enter the merging treatment of detector seismic data and land detector seismic data in water-filling, reach elimination earthquake
The influence of the multiple wave interference of seawater singing in data, improve the purpose of geological data signal to noise ratio and resolution ratio.
In order to solve the above technical problems, the embodiment of the present application provides a kind of side of land and water detector seismic data merging treatment
What method and device were realized in：
A kind of method of land and water detector seismic data merging treatment, there is provided have wave detector earthquake number in the water of target area
Include according to land detector seismic data, methods described：
Phase Processing is carried out to the land detector seismic data, obtains the land wave detector earthquake number after Phase Processing
According to；
First the amplitude processing is carried out to detector seismic data in the water, with obtaining in the water after the amplitude processing wave detector
Data are shaken, and the second the amplitude processing is carried out to the land detector seismic data after the Phase Processing, obtain the amplitude processing
Land detector seismic data afterwards；
Detector seismic data in water after the amplitude processing and land detector seismic data are carried out in Fu respectively
Leaf transformation, obtain the land wave detector earthquake number after detector seismic data and Fourier transformation in the water after Fourier transformation
According to；
According to detector seismic data in the water after the amplitude processing and the land wave detector after the amplitude processing
Shake data and calculate calibration factor, and according to detector seismic data in the water after the Fourier transformation and the Fourier transformation
Land detector seismic data afterwards, calculate the sea water advanced and bottom reflection coefficient of the target area；
Based on described sea water advanced, mirror image connector is determined；
Based on the calibration factor, described sea water advanced and described bottom reflection coefficient, and the mirror image connector,
Mirror image is carried out to the land detector seismic data of detector seismic data in the water of the target area and the target area
Combined treatment, obtain the land and water wave detector merging data of the target area.
In preferred scheme, the mirror image combined treatment is carried out using following formula：
S=OH[n]*H+OG[n]*G
Wherein, S represents the land and water wave detector merging data, and H represents wave detector earthquake number in the water of the target area
According to G represents the land detector seismic data of the target area；* convolution operation symbol is represented；OH[n] represents water inspection data
Combination operators, OG[n] expression land inspection data combination operators, the n expression worthwhile child-sequence-numbers of mirror set, n=-NOp ,-NOp+1 ,-
NOp+2 ..., -1,0,1 ..., NOp-2, NOp-1, NOp, (2NOp+1) represent mirror image combination operators length；Δ t represents earthquake
The data time sampling interval；ω represents angular frequency；α represents optimal calibration factor；A represents that data it is expected that absolute value shakes after merging
Width；AHRepresent wave detector trace gather data average absolute value amplitude in water, AGRepresent that land wave detector trace gather data average root-mean-square is shaken
Width；RsRepresent water-surface reflection coefficient；R represents optimal bottom reflection coefficient；Z represents delay operator, and (1-Z) and (1+Z) represents institute
State mirror image connector；L is imaginary unit, and l2=-1, ω represents angular frequency, and D represents sea water advanced, and V represents seawater speed,
Φ represents the phase correction factor obtained when carrying out Phase Processing to land detector seismic data.
In preferred scheme, the first the amplitude processing is carried out to detector seismic data in the water using following formula：
Wherein,Represent detector seismic data in the water after the amplitude processing, Hi,jRepresent wave detector earthquake in water
Data, AHRepresent that whole road catchments middle detector seismic data average absolute value amplitude, the serial number of window library track, i when i is represented
=1,2 ..., L, L window parameter total road number when representing, the serial number of window parameter temporal sampling point when j is represented, j=1,2 ..., N, N
Window parameter temporal number of samples during expression, symbol | | | | for the operator that takes absolute value.
In preferred scheme, second is carried out to the land detector seismic data after the Phase Processing using following formula and is shaken
Width processing：
Wherein,Represent the land detector seismic data after the amplitude processing, Gi,jAfter representing the Phase Processing
Land detector seismic data, AGRepresent whole trace gather land detector seismic data average absolute value amplitude, window when i is represented
The serial number of library track, i=1,2 ..., L, L represent when the total road number of window parameter, j represent when window parameter temporal sampling point order
Number, j=1,2 ..., N, N window parameter temporal number of samples when representing, symbol | | | | to take absolute value operator.
In preferred scheme, in the water according to after the amplitude processing after detector seismic data and the amplitude processing
Land detector seismic data calculate calibration factor, including：
The auto-correlation function of detector seismic data and the amplitude processing in the water after the amplitude processing are calculated respectively
The auto-correlation function of land detector seismic data afterwards, and with detector seismic data in the water after the amplitude processing and
The cross-correlation function that land detector seismic data after the amplitude processing is associated；
According to detector seismic data in the water after the amplitude processing and the land wave detector after the amplitude processing
The auto-correlation function of data, and the cross-correlation function are shaken, calculates maximum variance module equation coefficient；
Calibration factor characteristic equation coefficient is calculated according to the maximum variance module equation coefficient；
According to the calibration factor characteristic equation coefficients to construct calibration factor characteristic equation, and it is special to solve the calibration factor
Equation is levied, obtains calibration factor characteristic equation root；
Maximum variance module, and root are calculated according to the maximum variance module equation coefficient and the calibration factor characteristic equation root
Optimal calibration factor is determined according to the maximum variance module.
It is described to be calculated according to the maximum variance module equation coefficient and the calibration factor characteristic equation root in preferred scheme
Maximum variance module, and optimal calibration factor is determined according to the maximum variance module, including：
Maximum variance module is calculated using following formula：
Wherein, Q (αn) represent the maximum variance module, PmRepresent the maximum variance module equation coefficient, m=1,2 ...,
8, αnRepresent the calibration factor characteristic equation root, n=1,2 ..., 7；
The maximum of the maximum variance module is determined using following formula, and optimal mark is determined according to the maximum variance module
Determine the factor：
Wherein, Q (αbest) represent the maximum of the maximum variance module, αbestRepresent the maximum with the maximum variance module
Optimal calibration factor corresponding to value.
In preferred scheme, detector seismic data and the Fourier become in the water according to after the Fourier transformation
Land detector seismic data after changing, the sea water advanced of the target area is calculated, including：
According to detector seismic data in the water after the Fourier transformation and the land detection after the Fourier transformation
Device geological data, calculate the average coherence spectra of land and water detector seismic data；
According to the average coherence spectra of the land and water detector seismic data, it is averagely mutual to calculate land and water detector seismic data
Close function；
According to the average cross-correlation function of the land and water detector seismic data, the sea water advanced of the target area is calculated.
In preferred scheme, detector seismic data and the Fourier become in the water according to after the Fourier transformation
Land detector seismic data after changing, the bottom reflection coefficient of the target area is calculated, including：
According to detector seismic data in the water after the Fourier transformation and the land detection after the Fourier transformation
Device geological data, calculate the first wave field data splitting, the second wave field data splitting and the 3rd wave field data splitting；
According to the first wave field data splitting, the second wave field data splitting and the 3rd wave field data splitting, seabed is calculated
Reflection coefficient characteristic equation coefficient；
According to the bottom reflection coefficient characteristic equation coefficients to construct bottom reflection coefficient characteristic equation, and solve the sea
Bottom reflection coefficient characteristic equation, obtain multiple bottom reflection coefficient values；
Calculated respectively according to each bottom reflection coefficient value corresponding respectively with each bottom reflection coefficient value
Wavefield data energy, and the optimal sub-bottom reflection in the multiple bottom reflection coefficient value is determined according to the wavefield data energy
Coefficient.
In preferred scheme, detector seismic data and the Fourier become in the water according to after the Fourier transformation
Land detector seismic data after changing, calculate the first wave field data splitting, the second wave field data splitting and the combination of the 3rd wave field
Data, including：
According to detector seismic data in the water after the Fourier transformation and the land detection after the Fourier transformation
Device geological data, calculate upstream wave field data, first-order lag upstream wave field data, first-order lag down-going wave fields data, second order and prolong
Slow upstream wave field data, scond-order lag down-going wave fields data and three ranks delay down-going wave fields data；
Based on the upstream wave field data and the first-order lag down-going wave fields data, the first wave field number of combinations is calculated
According to；
Based on the first-order lag upstream wave field data and the scond-order lag down-going wave fields data, second ripple is calculated
Field data splitting；
Based on the scond-order lag upstream wave field data and three rank delay down-going wave fields data, the 3rd ripple is calculated
Field data splitting.
A kind of device of land and water detector seismic data merging treatment, described device provide target area water in wave detector
Geological data and land detector seismic data, described device include：Phase processing module, the amplitude processing module, Fourier become
Change the mold block, parameter calculating module, mirror image connector determining module and land and water inspection data combiners block；Wherein,
The phase processing module, for carrying out Phase Processing to the land detector seismic data, obtain at phase
Land detector seismic data after reason；
The amplitude processing module, for carrying out the first the amplitude processing to detector seismic data in the water, shaken
Detector seismic data in water after width processing, and second is carried out to the land detector seismic data after the Phase Processing
The amplitude processing, obtain the land detector seismic data after the amplitude processing；
The fourier transformation module, for detector seismic data in the water after the amplitude processing and land detection
Device geological data carries out Fourier transformation respectively, obtains detector seismic data and Fourier transformation in the water after Fourier transformation
Land detector seismic data afterwards；
The parameter calculating module, for according to detector seismic data in the water after the amplitude processing and the amplitude
Land detector seismic data after processing calculates calibration factor, and according to wave detector earthquake in the water after the Fourier transformation
Land detector seismic data after data and the Fourier transformation, the sea water advanced and seabed for calculating the target area are anti-
Penetrate coefficient；
The mirror image connector determining module, for based on described sea water advanced, determining mirror image connector；
Data combiners block is examined in the land and water, for anti-based on the calibration factor, described sea water advanced and described seabed
Coefficient, and the mirror image connector are penetrated, to detector seismic data and the target area in the water of the target area
Land detector seismic data carry out mirror image combined treatment, obtain the land and water wave detector merging data of the target area.
The embodiment of the present application provides a kind of method and device of land and water detector seismic data merging treatment, it is not necessary to pre-
First set relevant parameter value range and scanning step, but optimal calibration factor and optimal is directly calculated using correlation function algorithm
Bottom reflection coefficient, amount of calculation is small, calculating speed is fast；Simultaneously by advance to the Phase Processing of land detector seismic data with
And the amplitude processing of land and water detector seismic data, reduce low-frequency noise and reduce amplitude energy interference, so as to improve
Sea water advanced computational accuracy.Calibration factor, sea water advanced and bottom reflection coefficient of the embodiment of the present application according to calculating, draw simultaneously
Enter mirror image connector, fast and accurately realize the merging of detector seismic data and land detector seismic data in water
Processing, influenceed so as to effectively eliminate the multiple wave interference of seawater singing in geological data, improve geological data signal to noise ratio and
Resolution ratio.
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 land and water detector seismic data merging treatment of the application；
Fig. 2 (a) is that the inspection of the embodiment land and water of the application one merges common detector gather water inspection schematic diagram data；
Fig. 2 (b) is that the inspection of the embodiment land and water of the application one merges common detector gather land inspection schematic diagram data；
Fig. 2 (c) is that existing land and water inspection merges common detector gather land and water inspection merging data schematic diagram；
Fig. 2 (d) is that the inspection of the embodiment land and water of the application one merges common detector gather land and water inspection merging data schematic diagram；
Fig. 3 (a) is that the inspection of the embodiment land and water of the application one merges common detector gather water inspection data spectrum schematic diagram；
Fig. 3 (b) is that the inspection of the embodiment land and water of the application one merges common detector gather land inspection data spectrum schematic diagram；
Fig. 3 (c) is that existing land and water inspection merges common detector gather land and water inspection merging data spectrum diagram；
Fig. 3 (d) is that the inspection of the embodiment land and water of the application one merges common detector gather land and water inspection merging data spectrum diagram；
Fig. 4 (a) is that the inspection of the embodiment land and water of the application one merges common-shot-gather water inspection schematic diagram data；
Fig. 4 (b) is that the inspection of the embodiment land and water of the application one merges common-shot-gather land inspection schematic diagram data；
Fig. 4 (c) is that existing land and water inspection merges common-shot-gather land and water inspection merging data schematic diagram；
Fig. 4 (d) is that the inspection of the embodiment land and water of the application one merges common-shot-gather land and water inspection merging data schematic diagram；
Fig. 5 (a) is that the inspection of the embodiment land and water of the application one merges common-shot-gather water inspection data spectrum schematic diagram；
Fig. 5 (b) is that the inspection of the embodiment land and water of the application one merges common-shot-gather land inspection data spectrum schematic diagram；
Fig. 5 (c) is that existing land and water inspection merges common-shot-gather land and water inspection merging data spectrum diagram；
Fig. 5 (d) is that the inspection of the embodiment land and water of the application one merges common-shot-gather land and water inspection merging data spectrum diagram；
Fig. 6 (a) is that the inspection of the embodiment land and water of the application one merges common midpoint gather superposition water inspection schematic diagram data；
Fig. 6 (b) is that the inspection of the embodiment land and water of the application one merges common midpoint gather superposition land inspection schematic diagram data；
Fig. 6 (c) is that existing land and water inspection merges common midpoint gather superposition land and water inspection merging data schematic diagram；
Fig. 6 (d) is that the inspection of the embodiment land and water of the application one merges common midpoint gather superposition land and water inspection merging data schematic diagram；
Fig. 7 (a) is that the inspection of the embodiment land and water of the application one merges common midpoint gather superposition water inspection data spectrum schematic diagram；
Fig. 7 (b) is that the inspection of the embodiment land and water of the application one merges common midpoint gather superposition land inspection data spectrum schematic diagram；
Fig. 7 (c) is that existing land and water inspection merges common midpoint gather superposition land and water inspection merging data spectrum diagram；
Fig. 7 (d) is that the embodiment land and water of the application one inspection merging common midpoint gather superposition land and water inspection merging data frequency spectrum shows
It is intended to；
Fig. 8 is the composition structure chart of the device embodiment of the application land and water detector seismic data merging treatment.
Embodiment
The embodiment of the present application provides a kind of method and device of land and water detector seismic data merging treatment.
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 land and water detector seismic data merging treatment.Methods described is provided with mesh
Mark detector seismic data and land detector seismic data in the water in region.
Fig. 1 is a kind of flow chart of the embodiment of the method for land and water detector seismic data merging treatment of the application.Such as Fig. 1 institutes
Show, the method for the land and water detector seismic data merging treatment, comprise the following steps.
S101：Phase Processing is carried out to the land detector seismic data, obtains the land wave detector after Phase Processing
Geological data.
Because the reception mechanism of wave detector and land wave detector in water is different, two kinds of wave detectors are caused to receive between data
Phase is different, therefore, before two kinds of data merge, for active balance land and water phase, it is necessary to land wave detector earthquake number
According to progress Phase Processing.Obtained by carrying out phase analysis to detector seismic data in actual water and land detector seismic data
To phase correction factor.
In the present embodiment, described Phase Processing, can specifically include：
Triumphant damp window Hilbert factor sequence h is calculated according to following formula：
Wherein, I0(x) it is first kind zeroth order modified Bessel function, β represents triumphant damp window function parameter, and M represents triumphant damp window
Function length, α=M/2.Triumphant damp window Hilbert factor sequence is described by two parameter betas and M.
Land detector seismic data Hilbert transform G is calculated according to following formulaH：
GH=G0* h (formula 2)
Wherein, symbol " * " represents convolution operation, and h is triumphant damp window Hilbert factor sequence.
Phase Processing is carried out to land detector seismic data according to following formula：
G=G0cosΦ-GHSin Φ (formula 3)
Wherein, G be Phase Processing after land detector seismic data, G0For land detector seismic data, GHFor land
The Hilbert transform of detector seismic data, Φ are phase correction factors.
S102：First the amplitude processing is carried out to detector seismic data in the water, obtains examining in the water after the amplitude processing
Ripple device geological data, and the second the amplitude processing is carried out to the land detector seismic data after the Phase Processing, shaken
Land detector seismic data after width processing.
In water during detector seismic data and land detector seismic data merging treatment, due to wave detector in water
It is different with the reception mechanism of land wave detector, cause two kinds of wave detectors to receive the energy between data different with effective width.Cause
This is easy to follow-up land and water detector seismic data to merge, it is necessary to carry out the amplitude processing to land and water detector seismic data.
In the present embodiment, following formula can be used to carry out at the first amplitude detector seismic data in the water
Reason：
Wherein,Represent detector seismic data in the water after the amplitude processing, Hi,jRepresent wave detector earthquake in water
Data, AHRepresent that whole road catchments middle detector seismic data average absolute value amplitude, the serial number of window library track, i when i is represented
=1,2 ..., L, L window parameter total road number when representing, the serial number of window parameter temporal sampling point when j is represented, j=1,2 ..., N, N
Window parameter temporal number of samples during expression, symbol | | | | for the operator that takes absolute value.
In the present embodiment, the land detector seismic data after the Phase Processing can be entered using following formula
The amplitude processing of row second：
Wherein,Represent the land detector seismic data after the amplitude processing, Gi,jAfter representing the Phase Processing
Land detector seismic data, AGRepresent whole trace gather land detector seismic data average absolute value amplitude, window when i is represented
The serial number of library track, i=1,2 ..., L, L represent when the total road number of window parameter, j represent when window parameter temporal sampling point order
Number, j=1,2 ..., N, N window parameter temporal number of samples when representing, symbol | | | | to take absolute value operator.
S103：Detector seismic data in water after the amplitude processing and land detector seismic data are carried out respectively
Fourier transformation, obtain the land wave detector earthquake after detector seismic data and Fourier transformation in the water after Fourier transformation
Data.
In the present embodiment, detector seismic data H in the water after Fourier transformation is calculated according to following formulai[k]
With the land detector seismic data G after Fourier transformationi[k]：
NL=2m≥N
Wherein, Hi,jFor detector seismic data in the water after the amplitude processing, Gi,jFor the land wave detector after the amplitude processing
Geological data；I represent when window library track serial number, i=1,2 ..., L, L represent when the total road number of window parameter, j represent when window join
Number the time sampling points serial number, j=1,2 ..., N, N represent when window parameter temporal number of samples.L is imaginary unit, and l2=-1,
K=0,1,2 ..., NL-1, NL be data Fourier transformation number of samples, m is an appropriate positive integer.
S104：According to detector seismic data in the water after the amplitude processing and the land detection after the amplitude processing
Device geological data calculates calibration factor, and according to detector seismic data in the water after the Fourier transformation and the Fourier
Land detector seismic data after conversion, calculate the sea water advanced and bottom reflection coefficient of the target area.
In the OBC data comprising the more subwaves of seawater singing, due to the energy of data be by subsurface reflective ripple energy and
The energy of the more subwaves of seawater singing is formed.The auto-correlation of data includes the auto-correlation and the more subwaves of seawater singing for the lower back wave that lands
Auto-correlation.Because the auto-correlation function energy of subsurface reflective ripple concentrates on auto-correlation function zero crossings, and seawater singing is more
The auto-correlation energy of subwave is located remotely from auto-correlation function zero point.After eliminating the more subwaves of seawater singing, due in geological data,
Having eliminated multiple wave energy, be left the energy of subsurface reflective ripple, i.e. zero crossings auto-correlation function value increases, and away from zero
Auto-correlation function value reduces at point.Auto-correlation function variance mould can reflect this relation, thus in order to effectively estimate demarcation because
Son, it is minimum as determination calibration factor criterion usually using auto-correlation function variance mould.In addition, eliminating the more subwaves of seawater singing
Afterwards, due in geological data, having eliminated multiple wave energy, it is left the energy of subsurface reflective ripple.Therefore in order to effectively estimate
Bottom reflection coefficient, it is minimum as the criterion for determining bottom reflection coefficient usually using energy.
In the present embodiment, after according to detector seismic data in the water after the amplitude processing and the amplitude processing
Land detector seismic data calculate calibration factor, can specifically include：
(1) calculate respectively in the water after the amplitude processing at the auto-correlation function of detector seismic data and the amplitude
The auto-correlation function of land detector seismic data after reason, and with detector seismic data in the water after the amplitude processing
The associated cross-correlation function with the land detector seismic data after the amplitude processing；
Specifically, land and water inspection combination is carried out according to following formula：
Wherein, Si,jLand and water detector array data are represented,Represent wave detector earthquake in the water after the amplitude processing
Data,The land detector seismic data after the amplitude processing is represented, α represents calibration factor；Window library track when i is represented
Serial number, i=1,2 ..., L, L represent when the total road number of window parameter, j represent when window parameter temporal sampling point serial number, j=1,
2 ..., N, N represent when window parameter temporal number of samples；
Every one of land and water detector array auto-correlation function A is calculated according to following formulai,k：
Wherein, k represents correlation function order of delay number, and M represents the branch length of correlation function half；
Formula 7 is substituted into formula 8, and deployed, then is had：
The both sides of formula 9 are summed to road serial number i, then had：
Order：
Wherein, akRepresent detector seismic data average autocorrelation function in water, bkRepresent that land and water detector seismic data is put down
Equal cross-correlation function, ckRepresent land detector seismic data average autocorrelation function；
(2) according to detector seismic data in the water after the amplitude processing and the land wave detector after the amplitude processing
The auto-correlation function of geological data, and the cross-correlation function, calculate maximum variance module equation coefficient；
Specifically, formula 11 is substituted into formula 10, then had：
Ak=ak+bkα+ckα2(formula 12)
By the both sides square of formula 12, then have：
By the both sides square of formula 13, then have：
The both sides of formula 14 are summed to auto-correlation function serial number k, then are had：
In formula,
Auto-correlation function maximum variance module calculation formula is
Wherein, PmExpression maximum variance module equation coefficient, m=0,1,2 ..., 8；
(3) calibration factor characteristic equation coefficient is calculated according to the maximum variance module equation coefficient；
Specifically, the both sides of formula 17 then have to calibration factor α derivations：
And it is 0 to make its derivative, obtains calibration factor α 7 rank characteristic equations
T0+T1α+T2α2+…+T7α7=0 (formula 19)
Wherein, calibration factor characteristic equation coefficient TkFor
Tk=(k+1) Pk+1, k=0,1,2 ..., 7 (formula 20)
(4) according to the calibration factor characteristic equation coefficients to construct calibration factor characteristic equation, and solve the demarcation because
Subcharacter equation, obtain calibration factor characteristic equation root；
Specifically, according to calibration factor characteristic equation coefficient, calibration factor characteristic equation, Ran Houqiu are constructed according to formula 19
Calibration factor characteristic equation is solved, obtains 7 characteristic root α corresponding to calibration factor characteristic equationn, n=1,2,3 ..., 7；
(5) maximum variance module is calculated according to the maximum variance module equation coefficient and the calibration factor characteristic equation root,
And optimal calibration factor is determined according to the maximum variance module.
Specifically, following formula can be used to calculate maximum variance module：
Wherein, Q (αn) represent the maximum variance module, PmRepresent the maximum variance module equation coefficient, m=1,2 ...,
8, αnRepresent the calibration factor characteristic equation root, n=1,2 ..., 7；
The maximum of the maximum variance module is determined using following formula, and optimal mark is determined according to the maximum variance module
Determine the factor：
Wherein, Q (αbest) represent the maximum of the maximum variance module, αbestRepresent the maximum with the maximum variance module
Optimal calibration factor corresponding to value.
In the present embodiment, become according to detector seismic data in the water after the Fourier transformation and the Fourier
Land detector seismic data after changing, calculates the sea water advanced of the target area, can specifically include：
(1) examined according to detector seismic data in the water after the Fourier transformation and the land after the Fourier transformation
Ripple device geological data, calculate the average coherence spectra of land and water detector seismic data；
Specifically, following formula can be used to calculate the average coherence spectra of the land and water detector seismic data：
Wherein, R [k] represents the average coherence spectra of the land and water detector seismic data, Hi[k] represents that the Fourier becomes
Detector seismic data in water after changing, Gi[k] represents the land detector seismic data after the Fourier transformation,WithRespectively with Hi[k] and Gi[k] is complex conjugate relationship；
(2) according to the average coherence spectra of the land and water detector seismic data, calculate land and water detector seismic data and be averaged
Cross-correlation function；
Specifically, the average cross-correlation function r [n] of land and water detector seismic data is calculated according to following formula：
Wherein, n=1,2,3 ..., NL；
(3) according to the average cross-correlation function of the land and water detector seismic data, the seawater for calculating the target area is deep
Degree.
Specifically, the position corresponding to the average cross-correlation function maximum of land and water detector seismic data is exactly seawater round trip
Hourage sample value m：
Sea water advanced D is calculated according to the seawater two-way travel time sample value m：
Wherein, Δ t represents geological data time sampling interval, and unit be the second (s), and V represents seawater speed, unit for rice/
Second (m/s).
In the present embodiment, become according to detector seismic data in the water after the Fourier transformation and the Fourier
Land detector seismic data after changing, the bottom reflection coefficient of the target area is calculated, can specifically be included：
(1) examined according to detector seismic data in the water after the Fourier transformation and the land after the Fourier transformation
Ripple device geological data, calculate the first wave field data splitting, the second wave field data splitting and the 3rd wave field data splitting；
Specifically, according to detector seismic data in the water after the Fourier transformation and the land after the Fourier transformation
Ground detector seismic data, calculate upstream wave field data, first-order lag upstream wave field data, first-order lag down-going wave fields data,
Scond-order lag upstream wave field data, scond-order lag down-going wave fields data and three ranks delay down-going wave fields data；
The more subwave inverse filter operator expression formulas of Backus seawater singings are：
In formula, R is bottom reflection coefficient, RsFor the reflection coefficient of sea surface, Z is delay operator, and its expression formula is
Z=elωτ(formula 28)
In formula, l is imaginary unit, and l2=-1, ω is angular frequency, and unit is radian per second (Ω/s), and τ is water layer round trip
Hourage, unit are the second (s), are with the relation between sea water advanced
In formula, V is seawater speed, and unit is meter per second (m/s), and D is sea water advanced, and unit is rice (m)；
According to focus singing and micro- multiple wave pattern of flexion reflex singing, in order to effectively eliminate focus singing and the ring of micro- flexion reflex
More subwaves are shaken, using the more subwave inverse filter operators of Backus seawater singings, are multiplied by upstream wave field and down-going wave fields data splitting.
I.e.
Si(ω)=[(1-Z) Hi(ω)+(1+Z)Gi(ω)] B (ω) (formula 30)
In formula, Si(ω) is to eliminate the multiple wave interference wave field Fourier transformation data of seawater singing, Hi(ω) becomes for Fourier
Detector seismic data in water after changing, Gi(ω) is the land detector seismic data after Fourier transformation, and i represents detection altogether
Point trace gather data track serial number；I=1,2,3 ..., L, L represent the total road number of common detector gather data；B (ω) is Backus seas
The more subwave inverse filter operators of water singing；
Formula 27 is substituted into formula 30, then had：
Its time-domain expression formula is
Wherein, " * " represents convolution operation, Si,jIt is to eliminate the multiple wave interference wavefield data of seawater singing,It is at amplitude
Detector seismic data in water after reason,It is the land detector seismic data after the amplitude processing；O1[j] represents single order
The more subwave inverse filter operators of Backus seawater singings, also referred to as first-order lag operator；O2[j] represents the ring of second order Backus seawater
Shake more subwave inverse filter operators, also referred to as scond-order lag operator；O3[j] represents the inverse filter of the three more subwaves of rank Backus seawater singings
Ripple device operator, also referred to as three rank delay operators；And
O1[j]=δ [j- τ/Δ t]
O2[j]=δ [j-2 τ/Δ t]
O3[j]=δ [j-3 τ/Δ t] (formula 33)
In formula, δ [j] represents unit pulse operator, and
Based on the upstream wave field data and the first-order lag down-going wave fields data, the first wave field number of combinations is calculated
According to；Based on the first-order lag upstream wave field data and the scond-order lag down-going wave fields data, the second wave field group is calculated
Close data；Based on the scond-order lag upstream wave field data and three rank delay down-going wave fields data, the 3rd ripple is calculated
Field data splitting；
Order
ui,j=(Hi,j+Gi,j)-(Hi,j-Gi,j)*O1[j]
vi,j=-2Rs(Hi,j+Gi,j)*O1[j]+2Rs(Hi,j-Gi,j)*O2[j]
In formula, ui,jRepresent the first wave field data splitting, vi,jRepresent the second wave field data splitting, wi,jRepresent the 3rd wave field
Data splitting；
(2) according to the first wave field data splitting, the second wave field data splitting and the 3rd wave field data splitting, sea is calculated
Bottom reflection coefficient characteristic equation coefficient；
Specifically, formula 35 is substituted into formula 32, then had：
Si,j=ui,j+vi,jR+wi,jR2(formula 36)
The both sides square of formula 36, then have：
The both sides of formula 37 are summed to i, then are had：
Order
Wherein, E represents data capacity, d0、d1、d2、d3And d4Represent bottom reflection coefficient characteristic equation coefficient；
(3) according to the bottom reflection coefficient characteristic equation coefficients to construct bottom reflection coefficient characteristic equation, and institute is solved
Bottom reflection coefficient characteristic equation is stated, obtains multiple bottom reflection coefficient values；
Specifically, formula 39 is substituted into formula 38, then had：
E=d0+d1R+d2R2+d3R3+d4R4(formula 40)
The both sides of formula 40 are differentiated to bottom reflection coefficient R, then are had：
It is 0 to make energy derivative, obtains bottom reflection coefficient R characteristic equation：
d1+2d2R+3d3R2+4d4R3=0 (formula 42)
Characteristic formula 42 is solved, three bottom reflection coefficient numerical value R can be obtained1、R2And R3；
(4) calculated respectively according to each bottom reflection coefficient value corresponding respectively with each bottom reflection coefficient value
Wavefield data energy, and determine that the optimal seabed in the multiple bottom reflection coefficient value is anti-according to the wavefield data energy
Penetrate coefficient.
Specifically, these three bottom reflection coefficient numerical value are substituted into formula 40, obtains three energy value E1、E2And E3.I.e.
The optimal bottom reflection coefficient R in the multiple bottom reflection coefficient value is determined according to following formulabest：
S105：Based on described sea water advanced, mirror image connector is determined.
In the present embodiment, the mirror image connector can be determined using following formula：
Wherein, Z represents delay operator, (1-Z) and (1+Z) expression mirror image connector, and l is imaginary unit, and l2
=-1, ω represents angular frequency, and D represents sea water advanced, and V represents seawater speed.
S106：Combined based on the calibration factor, described sea water advanced and described bottom reflection coefficient, and the mirror image
The factor, the land detector seismic data of detector seismic data in the water of the target area and the target area is carried out
Mirror image combined treatment, obtain the land and water wave detector merging data of the target area.
In the present embodiment, following formula can be used to carry out the mirror image combined treatment：
S=OH[n]*H+OG[n]*G
Wherein, S represents the land and water wave detector merging data, and H represents wave detector earthquake number in the water of the target area
According to G represents the land detector seismic data of the target area；* convolution operation symbol is represented；OH[n] represents water inspection data
Combination operators, OG[n] expression land inspection data combination operators, the n expression worthwhile child-sequence-numbers of mirror set, n=-NOp ,-NOp+1 ,-
NOp+2 ..., -1,0,1 ..., NOp-2, NOp-1, NOp, (2NOp+1) represent mirror image combination operators length；Δ t represents earthquake
The data time sampling interval；ω represents angular frequency；α represents optimal calibration factor；A represents that data it is expected that absolute value shakes after merging
Width；AHRepresent wave detector trace gather data average absolute value amplitude in water, AGRepresent that land wave detector trace gather data average root-mean-square is shaken
Width；RsRepresent water-surface reflection coefficient；R represents optimal bottom reflection coefficient；Z represents delay operator, and (1-Z) and (1+Z) represents institute
State mirror image connector；L is imaginary unit, and l2=-1, ω represents angular frequency, and D represents sea water advanced, and V represents seawater speed,
Φ represents the phase correction factor obtained when carrying out Phase Processing to land detector seismic data.
Before step S101, the embodiment of the present application can also include to detector seismic data in the water of the target area
Pre-processed with land detector seismic data；It is corresponding, wave detector earthquake in the water of target area described in step S102
Data are detector seismic data in pretreated water, and the land detector seismic data of the target area is after pre-processing
Land detector seismic data.Pretreatment described in the embodiment of the present application specifically includes：To wave detector earthquake in the water
Data and land detector seismic data put label, define observation system, detector seismic data in the water and land are examined
The operations such as ripple device geological data is separated, denoising, filtering, velocity analysis and overlap-add procedure, it is important to note, however, that above-mentioned
The pretreatment operation enumerated can also be handled using other pretreatment modes merely to the present invention is better described,
Itself please this is not construed as limiting.
After step S106, the embodiment of the present application can also include drawing land and water wave detector merging data section and storage water
Land wave detector merging data.
The embodiment of the method for the land and water detector seismic data merging treatment, it is not necessary to preset relevant parameter scope
Value and scanning step, but optimal calibration factor and optimal bottom reflection coefficient are directly calculated using correlation function algorithm, calculate
Measure that small, calculating speed is fast；Pass through the Phase Processing to land detector seismic data in advance and land and water wave detector earthquake simultaneously
The amplitude processing of data, reduce low-frequency noise and reduce amplitude energy interference, so as to improve sea water advanced computational accuracy.
Calibration factor, sea water advanced and bottom reflection coefficient of the embodiment of the present application according to calculating, while mirror image connector is introduced, soon
Speed, the merging treatment of detector seismic data and land detector seismic data in water accurately is realized, so as to effectively disappear
Except the multiple wave interference of seawater singing influences in geological data, geological data signal to noise ratio and resolution ratio are improved.
For the beneficial effect of explanation the embodiment of the present application of removing, illustrate below in conjunction with the accompanying drawings：
Common detector gather data comparison, Fig. 4 (a) to Fig. 4 (d) land and water are merged by Fig. 2 (a) to Fig. 2 (d) land and water inspection
Inspection, which merges common-shot-gather data comparison, Fig. 6 (a) to Fig. 6 (d) land and water inspection merges common midpoint gather data comparison to see
Go out, the water calculated using the embodiment of the present application examines data average root-mean-square amplitude, data average root-mean-square amplitude is examined in land, demarcation because
Son, bottom reflection coefficient, seawater TWT and the parameter such as sea water advanced, effectively eliminate more subwaves caused by water layer and do
Disturb, effectively widened cable data effective band, improve OBC data SNRs.
Abscissa and ordinate into Fig. 7 (d) are respectively to Fig. 5 (d) and Fig. 7 (a) by Fig. 3 (a) to Fig. 3 (d), Fig. 5 (a)
Frequency and amplitude.The contrast of common detector gather data spectrum, Fig. 5 (a) to Fig. 5 are merged by Fig. 3 (a) to Fig. 3 (d) land and water inspection
(d) land and water inspection merges the contrast of common-shot-gather data spectrum and Fig. 7 (a) to Fig. 7 (d) land and water inspection merges common midpoint gather
Data spectrum contrast is as can be seen that particularly Fig. 3 (a) to Fig. 3 (d), the inspection of the embodiment land and water of Fig. 3 (a) the application one merge inspection altogether
Wave point road, which catchments, examines data, and effective band low frequency Fw1=22.5 hertz (Hz), high frequency Fw2=44.0Hz, centre frequency is (depending on master
Frequently Fpw=33.25Hz), absolute frequency range Bw=21.5Hz, temporal resolution Trw=13.1ms；The embodiment of Fig. 3 (b) the application one
Land and water inspection merges common detector gather land inspection data, and effective band low frequency Fl1=18.0Hz, high frequency Fl2=32.0Hz, center is frequently
Rate (depending on dominant frequency) Fpl=25.0Hz, absolute frequency range Bl=14.0Hz, temporal resolution Trl=17.4ms；The existing land and waters of Fig. 3 (c)
Inspection merges common detector gather land and water inspection merging data, effective band low frequency Fc1=18.0Hz, high frequency Fc2=32.0Hz, center
Frequency (depending on dominant frequency) Fpc=25.0Hz, absolute frequency range Bc=14.0Hz, temporal resolution Trw=17.4ms；Fig. 3 (d) the application
The inspection of one embodiment land and water merges common detector gather land and water inspection merging data, effective band low frequency Fn1=6.5Hz, high frequency Fn2=
62.5Hz, centre frequency (depending on dominant frequency) Fpn=34.5Hz, absolute frequency range Bn=56.0Hz, temporal resolution Trn=12.6ms；
Centre frequency (depending on dominant frequency) Fp, absolute frequency range B and temporal resolution Tr these three parameters are to determine the pass of seismic data resolution
Bond parameter.Compared with examining with existing land and water and merge common detector gather land and water inspection merging data, the inspection of the application land and water merges detection altogether
Point trace gather land and water inspection merging data centre frequency (depending on dominant frequency) is bigger (big 9.5Hz), and absolute frequency range is wider (wide 42.0Hz), wherein
Low frequency is lower (low 11.5Hz), and high frequency is higher (high 30.5Hz), and temporal resolution is higher (high 4.8ms).Therefore, the application is utilized
The water that embodiment calculates examines data average root-mean-square amplitude, land inspection data average root-mean-square amplitude, calibration factor, sub-bottom reflection system
Number, seawater TWT and the parameter such as sea water advanced, have effectively widened cable data effective band, have improved OBC numbers
According to resolution ratio.
Fig. 8 is the composition structure chart of the device embodiment of the application land and water detector seismic data merging treatment.The dress
Put detector seismic data and land detector seismic data in the water that target area can be provided.As shown in figure 8, described device
It can include：Phase processing module 100, the amplitude processing module 200, fourier transformation module 300, parameter calculating module 400, mirror
As connector determining module 500 and land and water inspection data combiners block 600.
The phase processing module 100, it can be used for carrying out Phase Processing to the land detector seismic data, obtain
Land detector seismic data after Phase Processing.
The amplitude processing module 200, it can be used for carrying out the first the amplitude processing to detector seismic data in the water,
Detector seismic data in the water after the amplitude processing is obtained, and the land detector seismic data after the Phase Processing is entered
The amplitude processing of row second, obtain the land detector seismic data after the amplitude processing；
The fourier transformation module 300, can be used for detector seismic data in the water after the amplitude processing and
Land detector seismic data carries out Fourier transformation respectively, obtains detector seismic data and Fu in the water after Fourier transformation
In land detector seismic data after leaf transformation.
The parameter calculating module 400, can be used for according to detector seismic data in the water after the amplitude processing and
Land detector seismic data after the amplitude processing calculates calibration factor, and is examined according in the water after the Fourier transformation
Land detector seismic data after ripple device geological data and the Fourier transformation, calculates the sea water advanced of the target area
And bottom reflection coefficient.
The mirror image connector determining module 500, can be used for based on described sea water advanced, determine mirror image combination because
Son.
Data combiners block 600 is examined in the land and water, can be used for based on the calibration factor, described sea water advanced and described
Bottom reflection coefficient, and the mirror image connector, to detector seismic data and the mesh in the water of the target area
The land detector seismic data for marking region carries out mirror image combined treatment, and the land and water wave detector for obtaining the target area merges number
According to.
The device embodiment of the land and water detector seismic data merging treatment is closed with the land and water detector seismic data
And the embodiment of the method handled is corresponding, it is possible to achieve the embodiment of the method for the land and water detector seismic data merging treatment,
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 land and water detector seismic data merging treatment, it is characterised in that be provided with the water of target area and examine Ripple device geological data and land detector seismic data, methods described include：Phase Processing is carried out to the land detector seismic data, obtains the land detector seismic data after Phase Processing；First the amplitude processing is carried out to detector seismic data in the water, obtains wave detector earthquake number in the water after the amplitude processing According to, and the second the amplitude processing is carried out to the land detector seismic data after the Phase Processing, after obtaining the amplitude processing Land detector seismic data；Fourier's change is carried out respectively to detector seismic data in the water after the amplitude processing and land detector seismic data Change, obtain the land detector seismic data after detector seismic data and Fourier transformation in the water after Fourier transformation；According to the land wave detector earthquake number after detector seismic data in the water after the amplitude processing and the amplitude processing According to calculate calibration factor, and according to detector seismic data in the water after the Fourier transformation and the Fourier transformation after Land detector seismic data, calculate the sea water advanced and bottom reflection coefficient of the target area；Based on described sea water advanced, mirror image connector is determined；Based on the calibration factor, described sea water advanced and described bottom reflection coefficient, and the mirror image connector, to institute State detector seismic data and the land detector seismic data of the target area in the water of target area and carry out mirror image combination Processing, obtains the land and water wave detector merging data of the target area.
- 2. the method for a kind of land and water detector seismic data merging treatment according to claim 1, it is characterised in that use Following formula carry out the mirror image combined treatment：S=OH[n]*H+OG[n]*G<mrow> <msub> <mi>O</mi> <mi>H</mi> </msub> <mo>&lsqb;</mo> <mi>n</mi> <mo>&rsqb;</mo> <mo>=</mo> <mfrac> <mi>A</mi> <mrow> <msub> <mi>A</mi> <mi>H</mi> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msub> <mi>R</mi> <mi>s</mi> </msub> <mi>R</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mfrac> <msubsup> <mo>&Integral;</mo> <mrow> <mo>-</mo> <mi>&infin;</mi> </mrow> <mi>&infin;</mi> </msubsup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>Z</mi> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msub> <mi>R</mi> <mi>s</mi> </msub> <mi>R</mi> <mi>Z</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <msup> <mi>e</mi> <mrow> <mi>l</mi> <mi>&omega;</mi> <mi>n</mi> <mi>&Delta;</mi> <mi>t</mi> </mrow> </msup> <mi>d</mi> <mi>&omega;</mi> </mrow><mrow> <msub> <mi>O</mi> <mi>G</mi> </msub> <mo>&lsqb;</mo> <mi>n</mi> <mo>&rsqb;</mo> <mo>=</mo> <mfrac> <mrow> <mi>A</mi> <mi>&alpha;</mi> </mrow> <mrow> <msub> <mi>A</mi> <mi>G</mi> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msub> <mi>R</mi> <mi>s</mi> </msub> <mi>R</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mfrac> <msubsup> <mo>&Integral;</mo> <mrow> <mo>-</mo> <mi>&infin;</mi> </mrow> <mi>&infin;</mi> </msubsup> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>Z</mi> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msub> <mi>R</mi> <mi>s</mi> </msub> <mi>R</mi> <mi>Z</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <msup> <mi>e</mi> <mrow> <mi>l</mi> <mi>&Phi;</mi> </mrow> </msup> <msup> <mi>e</mi> <mrow> <mi>l</mi> <mi>&omega;</mi> <mi>n</mi> <mi>&Delta;</mi> <mi>t</mi> </mrow> </msup> <mi>d</mi> <mi>&omega;</mi> </mrow><mrow> <mi>Z</mi> <mo>=</mo> <msup> <mi>e</mi> <mrow> <mi>l</mi> <mi>&omega;</mi> <mfrac> <mrow> <mn>2</mn> <mi>D</mi> </mrow> <mi>V</mi> </mfrac> </mrow> </msup> </mrow>Wherein, S represents the land and water wave detector merging data, and H represents detector seismic data in the water of the target area, G Represent the land detector seismic data of the target area；* convolution operation symbol is represented；OH[n] represents water inspection data combination Operator, OG[n] represents land inspection data combination operators, and n represents the worthwhile child-sequence-number of mirror set, n=-NOp ,-NOp+1 ,-NOp+ 2 ..., -1,0,1 ..., NOp-2, NOp-1, NOp, (2NOp+1) represent mirror image combination operators length；Δ t represents geological data Time sampling interval；ω represents angular frequency；α represents optimal calibration factor；A represents that data it is expected absolute amplitude after merging；AH Represent wave detector trace gather data average absolute value amplitude in water, AGRepresent land wave detector trace gather data average root-mean-square amplitude；Rs Represent water-surface reflection coefficient；R represents optimal bottom reflection coefficient；Z represents delay operator, and (1-Z) and (1+Z) represents the mirror image Connector；L is imaginary unit, and l2=-1, ω represents angular frequency, and D represents sea water advanced, and V represents seawater speed, and Φ is represented The phase correction factor obtained when carrying out Phase Processing to land detector seismic data.
- 3. the method for a kind of land and water detector seismic data merging treatment according to claim 1, it is characterised in that use Following formula carry out the first the amplitude processing to detector seismic data in the water：<mrow> <msub> <mover> <mi>H</mi> <mo>^</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>=</mo> <mfrac> <msub> <mi>H</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <msub> <mi>A</mi> <mi>H</mi> </msub> </mfrac> </mrow><mrow> <msub> <mi>A</mi> <mi>H</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>N</mi> <mi>L</mi> </mrow> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>L</mi> </munderover> <munderover> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mo>|</mo> <mo>|</mo> <msub> <mi>H</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>|</mo> <mo>|</mo> </mrow>Wherein,Represent detector seismic data in the water after the amplitude processing, Hi,jDetector seismic data in water is represented, AHRepresent that whole road catchments middle detector seismic data average absolute value amplitude, the serial number of window library track when i is represented, i=1, 2 ..., L, L represent when the total road number of window parameter, j represent when window parameter temporal sampling point serial number, j=1,2 ..., N, N represent When window parameter temporal number of samples, symbol | | | | to take absolute value operator.
- 4. the method for a kind of land and water detector seismic data merging treatment according to claim 1, it is characterised in that use Following formula carry out the second the amplitude processing to the land detector seismic data after the Phase Processing：<mrow> <msub> <mover> <mi>G</mi> <mo>^</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>=</mo> <mfrac> <msub> <mi>G</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <msub> <mi>A</mi> <mi>G</mi> </msub> </mfrac> </mrow><mrow> <msub> <mi>A</mi> <mi>G</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>N</mi> <mi>L</mi> </mrow> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>L</mi> </munderover> <munderover> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mo>|</mo> <mo>|</mo> <msub> <mi>G</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>|</mo> <mo>|</mo> </mrow>Wherein,Represent the land detector seismic data after the amplitude processing, Gi,jRepresent the land after the Phase Processing Ground detector seismic data, AGWhole trace gather land detector seismic data average absolute value amplitude is represented, window refers to when i is represented The serial number in road, i=1,2 ..., L, L represent when the total road number of window parameter, j represent when window parameter temporal sampling point serial number, j= 1,2 ..., N, N window parameter temporal number of samples when representing, symbol | | | | to take absolute value operator.
- 5. the method for a kind of land and water detector seismic data merging treatment according to claim 1, it is characterised in that described According to the land detector seismic data meter after detector seismic data in the water after the amplitude processing and the amplitude processing Calibration factor is calculated, including：Calculate respectively in the water after the amplitude processing after the auto-correlation function of detector seismic data and the amplitude processing The auto-correlation function of land detector seismic data, and with detector seismic data in the water after the amplitude processing and described The cross-correlation function that land detector seismic data after the amplitude processing is associated；According to the land wave detector earthquake number after detector seismic data in the water after the amplitude processing and the amplitude processing According to auto-correlation function, and the cross-correlation function, calculate maximum variance module equation coefficient；Calibration factor characteristic equation coefficient is calculated according to the maximum variance module equation coefficient；According to the calibration factor characteristic equation coefficients to construct calibration factor characteristic equation, and solve the calibration factor feature side Journey, obtain calibration factor characteristic equation root；Maximum variance module is calculated according to the maximum variance module equation coefficient and the calibration factor characteristic equation root, and according to institute State maximum variance module and determine optimal calibration factor.
- 6. the method for a kind of land and water detector seismic data merging treatment according to claim 5, it is characterised in that described Maximum variance module is calculated according to the maximum variance module equation coefficient and the calibration factor characteristic equation root, and according to it is described most Generous differential mode determines optimal calibration factor, including：Maximum variance module is calculated using following formula：<mrow> <mi>Q</mi> <mrow> <mo>(</mo> <msub> <mi>&alpha;</mi> <mi>n</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>m</mi> <mo>=</mo> <mn>0</mn> </mrow> <mn>8</mn> </munderover> <msub> <mi>P</mi> <mi>m</mi> </msub> <msubsup> <mi>&alpha;</mi> <mi>n</mi> <mi>m</mi> </msubsup> </mrow>Wherein, Q (αn) represent the maximum variance module, PmRepresent the maximum variance module equation coefficient, m=1,2 ..., 8, αn Represent the calibration factor characteristic equation root, n=1,2 ..., 7；Determine the maximum of the maximum variance module using following formula, and according to the maximum variance module determine optimal demarcation because Son：<mrow> <mi>Q</mi> <mrow> <mo>(</mo> <msub> <mi>&alpha;</mi> <mi>best</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <munder> <mrow> <mi>Max</mi> <mo>{</mo> <mi>Q</mi> <mrow> <mo>(</mo> <msub> <mi>&alpha;</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>,</mo> <mi>Q</mi> <mrow> <mo>(</mo> <msub> <mi>&alpha;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>,</mo> <mo>.</mo> <mo>.</mo> <mo>.</mo> <mo>,</mo> <mi>Q</mi> <mrow> <mo>(</mo> <msub> <mi>&alpha;</mi> <mi>n</mi> </msub> <mo>)</mo> </mrow> <mo>}</mo> </mrow> <mrow> <msub> <mi>&alpha;</mi> <mi>best</mi> </msub> <mo>&Element;</mo> <mrow> <mo>(</mo> <msub> <mi>&alpha;</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>&alpha;</mi> <mn>2</mn> </msub> <mo>.</mo> <mo>.</mo> <mo>.</mo> <msub> <mi>&alpha;</mi> <mi>n</mi> </msub> <mo>)</mo> </mrow> </mrow> </munder> </mrow>Wherein, Q (αbest) represent the maximum of the maximum variance module, αbestRepresent the maximum pair with the maximum variance module The optimal calibration factor answered.
- 7. the method for a kind of land and water detector seismic data merging treatment according to claim 1, it is characterised in that described According to the land wave detector earthquake number after detector seismic data in the water after the Fourier transformation and the Fourier transformation According to, the sea water advanced of the target area is calculated, including：According to detector seismic data in the water after the Fourier transformation and the land wave detector after the Fourier transformation Data are shaken, calculate the average coherence spectra of land and water detector seismic data；According to the average coherence spectra of the land and water detector seismic data, the average cross-correlation letter of land and water detector seismic data is calculated Number；According to the average cross-correlation function of the land and water detector seismic data, the sea water advanced of the target area is calculated.
- 8. the method for a kind of land and water detector seismic data merging treatment according to claim 1, it is characterised in that described According to the land wave detector earthquake number after detector seismic data in the water after the Fourier transformation and the Fourier transformation According to, the bottom reflection coefficient of the target area is calculated, including：According to detector seismic data in the water after the Fourier transformation and the land wave detector after the Fourier transformation Data are shaken, calculate the first wave field data splitting, the second wave field data splitting and the 3rd wave field data splitting；According to the first wave field data splitting, the second wave field data splitting and the 3rd wave field data splitting, sub-bottom reflection is calculated Coefficient characteristics equation coefficient；According to the bottom reflection coefficient characteristic equation coefficients to construct bottom reflection coefficient characteristic equation, and it is anti-to solve the seabed Coefficient characteristics equation is penetrated, obtains multiple bottom reflection coefficient values；Calculated respectively according to each bottom reflection coefficient value and distinguish corresponding wave field with each bottom reflection coefficient value Data capacity, and determine according to the wavefield data energy the optimal sub-bottom reflection system in the multiple bottom reflection coefficient value Number.
- 9. the method for a kind of land and water detector seismic data merging treatment according to claim 8, it is characterised in that described According to the land wave detector earthquake number after detector seismic data in the water after the Fourier transformation and the Fourier transformation According to, the first wave field data splitting, the second wave field data splitting and the 3rd wave field data splitting are calculated, including：According to detector seismic data in the water after the Fourier transformation and the land wave detector after the Fourier transformation Shake data, calculate upstream wave field data, first-order lag upstream wave field data, first-order lag down-going wave fields data, on scond-order lag Traveling wave field data, scond-order lag down-going wave fields data and three ranks delay down-going wave fields data；Based on the upstream wave field data and the first-order lag down-going wave fields data, the first wave field data splitting is calculated；Based on the first-order lag upstream wave field data and the scond-order lag down-going wave fields data, the second wave field group is calculated Close data；Based on the scond-order lag upstream wave field data and three rank delay down-going wave fields data, the 3rd wave field group is calculated Close data.
- 10. a kind of device of land and water detector seismic data merging treatment, it is characterised in that described device provides target area Detector seismic data and land detector seismic data, described device include in water：Phase processing module, the amplitude processing mould Block, fourier transformation module, parameter calculating module, mirror image connector determining module and land and water inspection data combiners block；Wherein,The phase processing module, for carrying out Phase Processing to the land detector seismic data, after obtaining Phase Processing Land detector seismic data；The amplitude processing module, for carrying out the first the amplitude processing to detector seismic data in the water, obtain at amplitude Detector seismic data in water after reason, and the second amplitude is carried out to the land detector seismic data after the Phase Processing Processing, obtains the land detector seismic data after the amplitude processing；The fourier transformation module, for detector seismic data in the water after the amplitude processing and land wave detector Shake data carry out Fourier transformation respectively, obtain in the water after Fourier transformation after detector seismic data and Fourier transformation Land detector seismic data；The parameter calculating module, for according to detector seismic data in the water after the amplitude processing and the amplitude processing Land detector seismic data afterwards calculates calibration factor, and according to detector seismic data in the water after the Fourier transformation With the land detector seismic data after the Fourier transformation, the sea water advanced and sub-bottom reflection system of the target area is calculated Number；The mirror image connector determining module, for based on described sea water advanced, determining mirror image connector；Data combiners block is examined in the land and water, for based on the calibration factor, described sea water advanced and described sub-bottom reflection system Number, and the mirror image connector, to the land of detector seismic data and the target area in the water of the target area Ground detector seismic data carries out mirror image combined treatment, obtains the land and water wave detector merging data of the target area.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710565184.XA CN107589455B (en) | 2017-07-12 | 2017-07-12 | A kind of method and device of land and water detector seismic data merging treatment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710565184.XA CN107589455B (en) | 2017-07-12 | 2017-07-12 | A kind of method and device of land and water detector seismic data merging treatment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107589455A true CN107589455A (en) | 2018-01-16 |
| CN107589455B CN107589455B (en) | 2019-03-15 |
Family
ID=61042251
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710565184.XA Active CN107589455B (en) | 2017-07-12 | 2017-07-12 | A kind of method and device of land and water detector seismic data merging treatment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107589455B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108363096A (en) * | 2018-02-01 | 2018-08-03 | 中国石油天然气集团有限公司 | A kind of land and water detector seismic data uplink and downlink wave field separation method and device |
| CN111624660A (en) * | 2019-02-27 | 2020-09-04 | 中国石油天然气集团有限公司 | Method and device for separating up-and-down traveling wave fields of seismic data of submarine nodes |
| CN112147716A (en) * | 2020-09-15 | 2020-12-29 | 中海油田服务股份有限公司 | Method and device for merging land and water detection signals |
| CN112327360A (en) * | 2020-10-29 | 2021-02-05 | 中海油田服务股份有限公司 | Method, device, computer storage medium and terminal for realizing noise processing |
| CN112558181A (en) * | 2019-09-26 | 2021-03-26 | 中国石油天然气集团有限公司 | Sensitivity calibration method and device for near-field detector of marine air gun |
| CN116680630A (en) * | 2023-07-27 | 2023-09-01 | 成都雨航创科科技有限公司 | Human-vehicle motion detection method and device based on vibration and image |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8208342B2 (en) * | 2009-09-14 | 2012-06-26 | Pgs Geophysical As | Method for combining signals of pressure and particle motion sensors in marine seismic streamers |
| CN103675910A (en) * | 2013-11-29 | 2014-03-26 | 中国石油天然气集团公司 | Amphibious detector seismic data scaling factor retrieval method |
| CN104181586A (en) * | 2014-08-04 | 2014-12-03 | 中国石油集团东方地球物理勘探有限责任公司 | Inversion method of waterland detector data seabed reflection coefficient |
| CN105116445A (en) * | 2015-09-02 | 2015-12-02 | 中国石油集团东方地球物理勘探有限责任公司 | Method and apparatus of seismic data combination processing of amphibious detector |
-
2017
- 2017-07-12 CN CN201710565184.XA patent/CN107589455B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8208342B2 (en) * | 2009-09-14 | 2012-06-26 | Pgs Geophysical As | Method for combining signals of pressure and particle motion sensors in marine seismic streamers |
| CN103675910A (en) * | 2013-11-29 | 2014-03-26 | 中国石油天然气集团公司 | Amphibious detector seismic data scaling factor retrieval method |
| CN104181586A (en) * | 2014-08-04 | 2014-12-03 | 中国石油集团东方地球物理勘探有限责任公司 | Inversion method of waterland detector data seabed reflection coefficient |
| CN105116445A (en) * | 2015-09-02 | 2015-12-02 | 中国石油集团东方地球物理勘探有限责任公司 | Method and apparatus of seismic data combination processing of amphibious detector |
Non-Patent Citations (2)
| Title |
|---|
| P. HUGONNET ET AL.: "PZ Summation of 3D WAZ OBS Receiver Gathers", 《73RD EAGE CONFERENCE & EXHIBITION INCORPORATING SPE EUROPEC 2011》 * |
| 张绍东等: "《国内外石油技术进展（十一五）-石油物探》", 31 July 2012, 中国石化出版社 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108363096A (en) * | 2018-02-01 | 2018-08-03 | 中国石油天然气集团有限公司 | A kind of land and water detector seismic data uplink and downlink wave field separation method and device |
| CN108363096B (en) * | 2018-02-01 | 2020-02-14 | 中国石油天然气集团有限公司 | Method and device for separating up-going and down-going wave fields of seismic data of land and water detector |
| CN111624660A (en) * | 2019-02-27 | 2020-09-04 | 中国石油天然气集团有限公司 | Method and device for separating up-and-down traveling wave fields of seismic data of submarine nodes |
| CN111624660B (en) * | 2019-02-27 | 2023-05-26 | 中国石油天然气集团有限公司 | Method and device for separating up-and-down wave fields of submarine node seismic data |
| CN112558181A (en) * | 2019-09-26 | 2021-03-26 | 中国石油天然气集团有限公司 | Sensitivity calibration method and device for near-field detector of marine air gun |
| CN112147716A (en) * | 2020-09-15 | 2020-12-29 | 中海油田服务股份有限公司 | Method and device for merging land and water detection signals |
| CN112327360A (en) * | 2020-10-29 | 2021-02-05 | 中海油田服务股份有限公司 | Method, device, computer storage medium and terminal for realizing noise processing |
| CN112327360B (en) * | 2020-10-29 | 2022-12-09 | 中海油田服务股份有限公司 | Method, device, computer storage medium and terminal for realizing noise processing |
| CN116680630A (en) * | 2023-07-27 | 2023-09-01 | 成都雨航创科科技有限公司 | Human-vehicle motion detection method and device based on vibration and image |
| CN116680630B (en) * | 2023-07-27 | 2023-10-13 | 成都雨航创科科技有限公司 | Human-vehicle motion detection method and device based on vibration and image |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107589455B (en) | 2019-03-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107589455B (en) | A kind of method and device of land and water detector seismic data merging treatment | |
| CN105116445B (en) | A kind of method and device of land and water detector seismic data merging treatment | |
| CN103675910B (en) | A kind of land and water detector seismic data calibration factor inversion method | |
| CN106896409B (en) | A kind of varying depth cable ghost reflection drawing method based on wave equation boundary values inverting | |
| AU2008270922A1 (en) | Seismic inversion of data containing surface related multiples | |
| CN108363096B (en) | Method and device for separating up-going and down-going wave fields of seismic data of land and water detector | |
| AU2015238881B2 (en) | Methods and Systems that Attenuate Noise in Seismic Data | |
| Robinson et al. | Digital imaging and deconvolution: The ABCs of seismic exploration and processing | |
| CN107884829A (en) | A kind of method for combining compacting shallow sea OBC Multiple Attenuation in Seismic Data | |
| CN107193045A (en) | A kind of seismic data processing technique and device | |
| US10393901B2 (en) | Wavefield interpolation and regularization in imaging of multiple reflection energy | |
| CN104199087B (en) | Method and device for inverting sea water depth by use of data of underwater detector and land detector | |
| CN101576621B (en) | Method and device for processing data of submarine cable double-detection seismic exploration | |
| Toxopeus et al. | Simulating migrated and inverted seismic data by filtering a geologic model | |
| Métivier et al. | A review of the use of optimal transport distances for high resolution seismic imaging based on the full waveform | |
| CN111999766B (en) | Multi-wave-type wave field separation method and reflection and transmission coefficient acquisition method | |
| Gasperini et al. | Pseudo-3D techniques for analysis and interpretation of high-resolution marine seismic reflection data | |
| Sollid et al. | Scattering-angle migration of ocean-bottom seismic data in weakly anisotropic media | |
| Li et al. | Wavelet-based detection of singularities in acoustic impedances from surface seismic reflection data | |
| CN105572742B (en) | It is a kind of to determine sea water advanced method and apparatus | |
| Nilsen et al. | Inversion of reflection data | |
| US20150185348A1 (en) | Methods and systems for attenuating residual acoustic energy in seismic data | |
| CN104597500B (en) | Land geophone and hydrophone seismic data matching method | |
| Alerini et al. | Adaptive focusing window for seismic angle migration | |
| US20150003197A1 (en) | Inversion techniques using streamers at different depths |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |