CN108363096A

CN108363096A - A kind of land and water detector seismic data uplink and downlink wave field separation method and device

Info

Publication number: CN108363096A
Application number: CN201810101188.7A
Authority: CN
Inventors: 高少武; 张少华; 钱忠平; 王成祥; 方云峰; 王爱佳
Original assignee: BGP Inc; China Petroleum and Natural Gas Co Ltd
Current assignee: BGP Inc; China Petroleum and Natural Gas Co Ltd
Priority date: 2018-02-01
Filing date: 2018-02-01
Publication date: 2018-08-03
Anticipated expiration: 2038-02-01
Also published as: CN108363096B

Abstract

The embodiment of the present application discloses a kind of land and water detector seismic data uplink and downlink wave field separation method and device.The method is provided with detector seismic data and land detector seismic data in the water of target area, the method includes：Based on detector seismic data in the water and the land detector seismic data, initial uplink wavefield data and initial downlink wavefield data are determined；Based on the initial uplink wavefield data and the initial downlink wavefield data, traveling-wave field correlation function up and down is determined；According to the traveling-wave field correlation function up and down, object function is established；Based on the object function, target wave field separation coefficient is determined；Based on the target wave field separation coefficient, the initial uplink wavefield data and the initial downlink wavefield data, target upstream wave field data and target downlink wavefield data are determined.Technical solution provided by the embodiments of the present application can efficiently separate upstream wave field and down-going wave fields.

Description

A kind of land and water detector seismic data uplink and downlink wave field separation method and device

Technical field

This application involves oil exploration, exploitation and development technique field, more particularly to a kind of land and water detector seismic data Uplink and downlink wave field separation method and device.

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 it is also higher and higher.3-D seismics can accurately describing reservoir.Marine double inspection data acquire skill Art, including OBC (Ocean Bottom Cable, submarine cable), OBN (Ocean Bottom Node, seabed node) and double inspections Towing cable etc. is fixed together wave detector in water and land wave detector, receives pressure wavefield and speed that air gun source generates respectively Wave field is spent, to obtain high-resolution three-dimension seismic data.They all use a fixed arrangement, and wave detector is fixed on seabed Or on the towing cable of sea, seismic wave is received, and a boat only pulls focus and carries out seismic wave excitation.It is acquired in OBC data In, at least three ships：One energy source boat only pulls air gun source arrangement, carries out seismic wave excitation；One taken-over vessel is It is fixed, it is connected to submarine cable, receives seismic wave；One ship or several ships, build a submarine cable and recycle seabed Cable.In some special areas, in order to facilitate oil reservoir monitoring test, reception cable is retained in seabed by design experiment system For many years, to facilitate multiple data to acquire.In the acquisition of OBN data, at least two ships：One energy source boat, pulls air gun source Arrangement carries out seismic wave excitation；One fishline applicating boat, when operation at stationary arrangement position folding and unfolding wave detector (in water wave detector and Land wave detector), extract data.OBN acquisitions are adapted to more deep-sea domain.In the acquisition of double inspection streamer datas, a ship pulls Air gun source arranges, and carries out seismic wave excitation；Several towing cables are pulled simultaneously (wave detector in water is placed at stationary arrangement position With land wave detector), extract data.Underground medium high accuracy three-dimensional image in order to obtain, it is necessary to be accurately determined all big guns The position of point and receiving point.

At sea when double inspection data acquisitions, since seabed and sea are all stronger reflecting interfaces.With epicenter excitation One seismic wavelet reaches seabed or a reflection seismic wavelet from hypocentral location and reaches seabed from underground, in submarine cable Wave detector, incude and record this reflection seismic wavelet.This reflection wavelet continues up advance and reaches sea, by sea Then the reflection in face changes direction and propagates downwards, reaches seabed.Wave detector in submarine cable, incudes and records again This seismic wavelet.This seismic wavelet is reflected by seabed simultaneously, is then changed direction and is upwardly propagated, reaches sea, by To the reflection on sea, then changes direction and propagate downwards, reach seabed.This cyclic process repeats.And these are original anti- It is undesirable secondary and subsequently reach to penetrate seismic wavelet, is exactly seawater singing multiple wave (reverberation).What source wavelet generated Seawater singing multiple wave, referred to as focus singing multiple wave, the seawater singing multiple wave that subsurface reflective wave generates, referred to as micro- flexion reflex Singing multiple wave.Seawater singing multiple wave is maximum noise jamming in offshore seismic exploration data.It is multiple to eliminate seawater singing Wave noise jamming is mostly important step in marine seismic data processing.

The multiple wave interference of seawater singing is can effectively eliminate, using the Method of Deconvolution to restore primary event in phytal zone Seismic wavelet.But for the sea water advanced area more than 10 meters, primary reflection seismic wavelet and follow-up seawater singing multiple wave The time difference between interference is very big so that deconvolution algorithms cannot be removed effectively the lineups that the multiple wave interference of seawater singing generates. If the multiple wave interference of seawater singing cannot be removed effectively from seismic data, several wavelets represent a reflection circle Face can thus obscure geological fault interface.The data of OBC acquisitions, provide detector data and land in same position water Two kinds of data of detector data.Both data use wave detector and land wave detector in water to record respectively.Wave detector is remembered in water What is recorded is the pressure change that seismic wave generates；Land wave detector record is particle velocity variation.Due to both wave detectors Recording mechanism is different, for the multiple wave interference of seawater singing at same position, shows different characteristic.It is recorded with wave detector in water The multiple wave interference of seawater singing compare, it is special that the seawater singing multiple wave interference of land wave detector record shows polarity and amplitude Levy difference.The multiple wave interference of seawater singing of two kinds of wave detectors record, polarity is on the contrary, amplitude is different, and differs one and sea The proportional constant of bottom reflectance factor, this constant value are exactly calibration factor.Therefore, poor using this amplitude and polar character It is different, it can effectively eliminate the multiple wave interference of seawater singing.By calibrated land detector data and wave detector in corresponding water Data are added, and are obtained water-land geophone calibration data (upstream wave field data), are eliminated the multiple wave interference of seawater singing；It will mark Land detector data after fixed is subtracted each other with detector data in corresponding water, has been obtained and sea level related multiple data (down-going wave fields data) combine the processing such as deconvolution, denoising and offset with down-going wave fields, with further for subsequent upstream wave field Improve resolution ratio and signal-to-noise ratio.Noise when due to field acquisition so that also deposited in the upstream wave field data separated in this way There is also a small amount of remaining upstream wave field data in a small amount of remaining down-going wave fields data, down-going wave fields data.Therefore, in order to Upstream wave field and down-going wave fields are efficiently separated, there is an urgent need for a kind of methods of land and water detector seismic data uplink and downlink wave field separation.

Invention content

The purpose of the embodiment of the present application is to provide a kind of land and water detector seismic data uplink and downlink wave field separation method and dress It sets, to efficiently separate upstream wave field and down-going wave fields.

In order to solve the above technical problems, traveling-wave field divides a kind of land and water detector seismic data of the embodiment of the present application offer up and down It is realized in from method and device：

A kind of land and water detector seismic data uplink and downlink wave field separation method, with being provided in the water of target area wave detector Data and land detector seismic data are shaken, including：

Based on detector seismic data in the water and the land detector seismic data, initial uplink wave field number is determined According to initial downlink wavefield data；

Based on the initial uplink wavefield data and the initial downlink wavefield data, traveling-wave field correlation letter up and down is determined Number；

According to the traveling-wave field correlation function up and down, object function is established；

Based on the object function, target wave field separation coefficient is determined；

Based on the target wave field separation coefficient, the initial uplink wavefield data and the initial downlink wavefield data, Determine target upstream wave field data and target downlink wavefield data.

It is described based on detector seismic data and the land detector seismic data in the water in preferred embodiment, really Determine initial uplink wavefield data and initial downlink wavefield data, including：

By the sum of detector seismic data in the water and the land detector seismic data, as the initial uplink Wavefield data；

Detector seismic data in the water is subtracted into the land detector seismic data, and using result of calculation as institute State initial downlink wavefield data.

In preferred embodiment, the traveling-wave field correlation function up and down includes：Upstream wave field auto-correlation function, down-going wave fields are from phase Close function, upstream wave field and down-going wave fields cross-correlation function and down-going wave fields and upstream wave field cross-correlation function；Using following Formula determines the upstream wave field auto-correlation function：

Wherein, a (k) indicates the upstream wave field auto-correlation function, U_i,jIndicate that the initial uplink wavefield data, i indicate When window parameter road serial number, i=1,2 ..., M, M indicate when the total road number of window parameter, j indicate when window parameter temporal sampling point it is suitable Serial number, j=1, window parameter temporal number of samples when 2 ..., N, N are indicated, k expression correlation function delay sampling point serial numbers, k=0,1, 2 ..., L-1, L indicate correlation function length；

The down-going wave fields auto-correlation function is determined using following formula：

Wherein, b (k) indicates the down-going wave fields auto-correlation function, D_i,jIndicate the initial downlink wavefield data；

The upstream wave field and down-going wave fields cross-correlation function are determined using following formula：

Wherein, c (k) indicates the upstream wave field and down-going wave fields cross-correlation function；

The down-going wave fields and upstream wave field cross-correlation function are determined using following formula：

Wherein, d (k) indicates the down-going wave fields and upstream wave field cross-correlation function.

In preferred embodiment, the object function is established using following formula：

Wherein, Q₂Under indicating that the object function, a (k) indicate that the upstream wave field auto-correlation function, b (k) indicate described Traveling-wave field auto-correlation function, c (k) indicate that the upstream wave field and down-going wave fields cross-correlation function, d (k) indicate the down going wave Field and upstream wave field cross-correlation function, k indicate that correlation function postpones sampling point serial number, k=0, and 1,2 ..., L-1, L indicate correlation Function length；β indicates that the energy coefficient for including down-going wave fields in the initial uplink wavefield data, γ indicate the initial downlink Include the energy coefficient of upstream wave field in wavefield data.

In preferred embodiment, the target wave field separation coefficient includes：Target downlink wave field separation coefficient and target upgoing wave Field separation；Wherein, the target downlink wave field separation coefficient is under characterizing and including in the initial uplink wavefield data The target energy coefficient of traveling-wave field, the target upstream wave field separation are wrapped for characterizing in the initial downlink wavefield data Target energy coefficient containing upstream wave field；It is described to be based on the object function, determine target wave field separation coefficient, including：

The object function is pre-processed, pretreated object function is obtained；Wherein, the pretreated mesh Scalar functions are associated with down-going wave fields separation, upstream wave field separation；

The pretreated object function is calculated respectively to down-going wave fields separation and upstream wave field separation Single order local derviation；

Single order local derviation based on calculating is as a result, establish being associated with for down-going wave fields separation and upstream wave field separation System；

First upstream wave field separation is set, and is associated with described according to the first upstream wave field separation System, determines the target downlink wave field separation coefficient and the target upstream wave field separation.

It is described that the object function is pre-processed in preferred embodiment, including：

The object function is multiplied by specified parameter, and using result of calculation as the pretreated object function；Its In, characterize the specified parameter using following formula：

A=(1- β γ)⁴

Wherein, A indicates that the specified parameter, β indicate the energy for including down-going wave fields in the initial uplink wavefield data Coefficient, γ indicate the energy coefficient for including upstream wave field in the initial downlink wavefield data.

In preferred embodiment, the traveling-wave field correlation function up and down includes：Upstream wave field auto-correlation function, down-going wave fields are from phase Close function, upstream wave field and down-going wave fields cross-correlation function and down-going wave fields and upstream wave field cross-correlation function；Using following Formula establishes the incidence relation：

Wherein, β indicates down-going wave fields separation, includes down-going wave fields for characterizing in the initial uplink wavefield data Energy coefficient, γ indicate upstream wave field separation, for characterize in the initial downlink wavefield data include upstream wave field Energy coefficient；A (k) indicates that the upstream wave field auto-correlation function, b (k) indicate the down-going wave fields auto-correlation function, c (k) Indicate that the upstream wave field and down-going wave fields cross-correlation function, d (k) indicate the down-going wave fields and upstream wave field cross-correlation letter Number, k indicate that correlation function postpones sampling point serial number, k=0, and 1,2 ..., L-1, L indicate correlation function length.

It is described according to the first upstream wave field separation and the incidence relation in preferred embodiment, determine the mesh Down-going wave fields separation and the target upstream wave field separation are marked, including：

According to the first upstream wave field separation and the incidence relation, the first down-going wave fields separation is determined, And according to the first down-going wave fields separation and the incidence relation, the second upstream wave field separation is determined, and according to The second upstream wave field separation and the incidence relation determine the second down-going wave fields separation；

When the absolute value of the difference of the first upstream wave field separation and the second upstream wave field separation is less than The absolute value of the difference of predetermined coefficient threshold value and the first down-going wave fields separation and the second down-going wave fields separation When less than the predetermined coefficient threshold value, by the first upstream wave field separation and the second upstream wave field separation it The half of sum, as the target upstream wave field separation, and by the first down-going wave fields separation and described the The half of the sum of two down-going wave fields separations, as the target downlink wave field separation coefficient；Otherwise, by second uplink Wave field separation coefficient is as the first new upstream wave field separation, and using the second down-going wave fields separation as newly First down-going wave fields separation, and generate the second new upstream wave field separation and the second new down-going wave fields segregative line Number, until the absolute value of the difference of the first new upstream wave field separation and the second new upstream wave field separation Less than predetermined coefficient threshold value and the first new down-going wave fields separation and the second new down-going wave fields separation Absolute value of the difference be less than the predetermined coefficient threshold value, and by the first new upstream wave field separation and described new the The half of the sum of two upstream wave field separations, as the target upstream wave field separation, and by described new first The half of down-going wave fields separation and the sum of the second new down-going wave fields separation, as the target downlink wave field Separation.

In preferred embodiment, the target wave field separation coefficient includes：Target downlink wave field separation coefficient and target upgoing wave Field separation；Wherein, the target downlink wave field separation coefficient is under characterizing and including in the initial uplink wavefield data The target energy coefficient of traveling-wave field, the target upstream wave field separation are wrapped for characterizing in the initial downlink wavefield data Target energy coefficient containing upstream wave field；The target upstream wave field data are determined using following formula：

Wherein,Indicate the target upstream wave field data, U_i,jIndicate the initial uplink wavefield data, D_i,jIt indicates The initial downlink wavefield data, the serial number in window parameter road when i is indicated, i=1, window parameter total road when 2 ..., M, M are indicated Number, the serial number of window parameter temporal sampling point when j is indicated, j=1, window parameter temporal number of samples when 2 ..., N, N are indicated, k indicate phase It closes function and postpones sampling point serial number, k=0,1,2 ..., L-1, L indicate correlation function length；β ' indicates the target downlink wave Field separation, γ ' indicate the target upstream wave field separation；

The target downlink wavefield data is determined using following formula：

Wherein,Indicate the target downlink wavefield data.

A kind of land and water detector seismic data uplink and downlink wave field separation device, described device provide target area water in examine Wave device seismic data and land detector seismic data, described device include：Initial traveling-wave field determining module, correlation function up and down Traveling-wave field determining module above and below determining module, object function determining module, wave field separation coefficient determination module and target；Wherein,

The initial traveling-wave field determining module up and down, for based on detector seismic data in the water and land inspection Wave device seismic data, determines initial uplink wavefield data and initial downlink wavefield data；

The correlation function determining module, for being based on the initial uplink wavefield data and the initial downlink wave field number According to determining traveling-wave field correlation function up and down；

The object function determining module, for according to the traveling-wave field correlation function up and down, establishing object function；

The wave field separation coefficient determination module determines target wave field separation coefficient for being based on the object function；

Traveling-wave field determining module above and below the target, for based on the target wave field separation coefficient, the initial uplink Wavefield data and the initial downlink wavefield data, determine target upstream wave field data and target downlink wavefield data.

By the above technical solution provided by the embodiments of the present application as it can be seen that a kind of land and water wave detector provided by the embodiments of the present application Seismic data uplink and downlink wave field separation method and device utilizes the identified target letter that traveling-wave field correlation function is established up and down Number, quickly and accurately inverting can obtain target wave field separation coefficient, i.e. target downlink wave field separation coefficient and target upgoing wave Field separation；Wherein, target downlink wave field separation coefficient is used to characterize in initial uplink wavefield data comprising down-going wave fields Target energy coefficient, target upstream wave field separation is for characterizing the target for including upstream wave field in initial downlink wavefield data Energy coefficient；Finally, it is based on target wave field separation coefficient, detector seismic data and land wave detector earthquake in water may be implemented Upstream wave field and down-going wave fields efficiently separate in data, reach and eliminate seawater multiple wave interference effect in seismic data, to Seismic data signal-to-noise ratio and resolution ratio are improved, and upstream wave field and down-going wave fields data are provided for latter earthquake data processing.

Description of the drawings

In order to illustrate the technical solutions in the embodiments of the present application or in the prior art more clearly, to embodiment or will show below There is attached drawing needed in technology description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this Some embodiments described in application, for those of ordinary skill in the art, in the premise of not making the creative labor property Under, other drawings may also be obtained based on these drawings.

Fig. 1 is a kind of flow chart of land and water detector seismic data uplink and downlink wave field separation embodiment of the method for the application；

Fig. 2 is the schematic diagram of horizontal layer speed and formation thickness model in the embodiment of the present application；

Fig. 3 (a) is the diagrammatic cross-section of land detector seismic data in the embodiment of the present application；

Fig. 3 (b) is the diagrammatic cross-section of detector seismic data in water in the embodiment of the present application；

Fig. 4 (a) is the down going wave isolated from Fig. 3 (a) and Fig. 3 (b) using the application method in the embodiment of the present application The diagrammatic cross-section of field data；

Fig. 4 (b) is the upgoing wave isolated from Fig. 3 (a) and Fig. 3 (b) using the application method in the embodiment of the present application The diagrammatic cross-section of field data；

Fig. 5 (a) is the diagrammatic cross-section of wave detector common-shot-gather data in water in the embodiment of the present application；

Fig. 5 (b) is the diagrammatic cross-section of wave detector common-shot-gather data in land in the embodiment of the present application；

Fig. 5 (c) is the uplink separated from Fig. 5 (a) and Fig. 5 (b) using the application method in the embodiment of the present application The diagrammatic cross-section of wavefield data；

Fig. 5 (d) is the downlink separated from Fig. 5 (a) and Fig. 5 (b) using the application method in the embodiment of the present application The diagrammatic cross-section of wavefield data；

Fig. 6 (a) is the diagrammatic cross-section of wave detector common midpoint stack data in water in the embodiment of the present application；

Fig. 6 (b) is the diagrammatic cross-section of wave detector common midpoint stack data in land in the embodiment of the present application；

Fig. 6 (c) is the uplink separated from Fig. 6 (a) and Fig. 6 (b) using the application method in the embodiment of the present application The diagrammatic cross-section of wavefield data；

Fig. 6 (d) is the downlink separated from Fig. 6 (a) and Fig. 6 (b) using the application method in the embodiment of the present application The diagrammatic cross-section of wavefield data；

Fig. 7 is the composed structure signal of the application land and water detector seismic data uplink and downlink wave field separation device embodiment Figure.

Specific implementation mode

The embodiment of the present application provides a kind of land and water detector seismic data uplink and downlink wave field separation method and device.

In order to make those skilled in the art better understand the technical solutions in the application, below in conjunction with the application reality The attached drawing in example is applied, technical solutions in the embodiments of the present application is clearly and completely described, it is clear that described implementation Example is merely a part but not all of the embodiments of the present application.Based on the embodiment in the application, this field is common The every other embodiment that technical staff is obtained without creative efforts should all belong to the application protection Range.

The embodiment of the present application provides a kind of land and water detector seismic data uplink and downlink wave field separation method.The land and water inspection Wave device seismic data uplink and downlink wave field separation method is provided with detector seismic data and land wave detector in the water of target area Seismic data.

In the present embodiment, the double inspection towing cables in sea, submarine cable (Ocean Bottom Cable, OBC) can be passed through With marine double inspection data acquisition technologys such as seabed node (Ocean Bottom Node, OBN), examined in same position records water Wave device seismic data and land detector seismic data.

Fig. 1 is a kind of flow chart of land and water detector seismic data uplink and downlink wave field separation embodiment of the method for the application.Such as Shown in Fig. 1, the land and water detector seismic data uplink and downlink wave field separation method includes the following steps.

Step S101：Based on detector seismic data in the water and the land detector seismic data, determine initial Upstream wave field data and initial downlink wavefield data.

In the present embodiment, it is based on detector seismic data and the land detector seismic data in the water, really Determine initial uplink wavefield data and initial downlink wavefield data, can specifically include, it can be by wave detector earthquake number in the water According to the sum of with the land detector seismic data, as the initial uplink wavefield data.It can be by wave detector in the water Seismic data subtracts the land detector seismic data, and using result of calculation as the initial downlink wavefield data.For example, Following formula may be used and calculate the initial uplink wavefield data U_i,j：

U_i,j=H_i,j+G_i,j(formula 1)

Wherein, H_i,jIndicate detector seismic data in the water, G_i,jIndicate the land detector seismic data, i tables The serial number in window parameter road when showing, i=1, window parameter total road number when 2 ..., M, M are indicated, window parameter temporal sampling point when j is indicated Serial number, j=1, window parameter temporal number of samples when 2 ..., N, N are indicated.

Following formula may be used and calculate the initial downlink wavefield data D_i,j：

D_i,j=H_i,j-G_i,j(formula 2)

Wherein, H_i,jIndicate detector seismic data in the water, G_i,jIndicate that the land detector seismic data i is indicated When window parameter road serial number, i=1,2 ..., M, M indicate when the total road number of window parameter, j indicate when window parameter temporal sampling point it is suitable Serial number, j=1, window parameter temporal number of samples when 2 ..., N, N are indicated.

Due to the initial uplink wavefield data U_i,jIn include a small amount of down-going wave fields and the initial downlink wave Field data D_i,jIn include a small amount of upstream wave field, in this way, following formula may be used characterizes the initial uplink wave field respectively Data U_i,jWith the initial downlink wavefield data D_i,j：

Wherein,WithIndicate that pure upstream wave field data and pure down-going wave fields data, β indicate the initial uplink respectively Include the energy coefficient of down-going wave fields in wavefield data, γ indicates the energy for including upstream wave field in the initial downlink wavefield data Coefficient of discharge.Wherein, down-going wave fields are not included in the pure upstream wave field data, do not include uplink in the pure down-going wave fields data Wave field.

Based on formula 3 and formula 4, following formula may be used and characterize the pure upstream wave field dataWith it is described it is pure under Traveling wave field data

Wherein,WithIndicate that pure upstream wave field data and pure down-going wave fields data, β indicate the initial uplink respectively Include the energy coefficient of down-going wave fields in wavefield data, γ indicates the energy for including upstream wave field in the initial downlink wavefield data Coefficient of discharge.In this way, obtaining the energy coefficient comprising down-going wave fields and institute in the initial uplink wavefield data in subsequent step State in initial downlink wavefield data comprising upstream wave field energy coefficient after, pure upstream wave field data and pure downlink can be obtained Wavefield data.

Step S102：Based on the initial uplink wavefield data and the initial downlink wavefield data, traveling wave up and down is determined Field correlation function.

In the present embodiment, following formula may be used and establish the pure upstream wave field dataWith the pure downlink Wavefield dataBetween cross-correlation function P (k)：

Wherein, k indicates that correlation function postpones sampling point serial number, k=0, and 1,2 ..., L-1, L indicate correlation function length.

Formula 5 and formula 6 are substituted into formula 7, then had：

The traveling-wave field correlation function up and down can specifically include：Upstream wave field auto-correlation function, down-going wave fields auto-correlation Function, upstream wave field and down-going wave fields cross-correlation function and down-going wave fields and upstream wave field cross-correlation function.It may be used down It states formula and determines the upstream wave field auto-correlation function：

Following formula may be used and determine the down-going wave fields auto-correlation function：

Following formula may be used and determine the upstream wave field and down-going wave fields cross-correlation function：

Following formula may be used and determine the down-going wave fields and upstream wave field cross-correlation function：

In this way, can formula 8 be substituted into for formula 10, formula 11, formula 12 and formula 13, then have：

Step S103：According to the traveling-wave field correlation function up and down, object function is established.

In the present embodiment, following formula may be used and establish the object function：

Wherein, Q₂Indicate that the object function, a (k) indicate that the upstream wave field auto-correlation function, b (k) indicate the down-going wave fields auto-correlation Function, c (k) indicate that the upstream wave field and down-going wave fields cross-correlation function, d (k) indicate the down-going wave fields and upstream wave field Cross-correlation function, k indicate that correlation function postpones sampling point serial number, k=0, and 1,2 ..., L-1, L indicate correlation function length；β tables Show that the energy coefficient for including down-going wave fields in the initial uplink wavefield data, γ indicate to wrap in the initial downlink wavefield data Energy coefficient containing upstream wave field.

Step S104：Based on the object function, target wave field separation coefficient is determined.

In the present embodiment, the target wave field separation coefficient can specifically include：Target downlink wave field separation coefficient With target upstream wave field separation.Wherein, the target downlink wave field separation coefficient is for characterizing the initial uplink wave field Include the target energy coefficient of down-going wave fields in data, the target upstream wave field separation is for characterizing the initial downlink Include the target energy coefficient of upstream wave field in wavefield data.

In the present embodiment, it is based on the object function, target wave field separation coefficient is determined, can specifically include：

(1) object function is pre-processed, obtains pretreated object function；Wherein, after the pretreatment Object function it is associated with down-going wave fields separation, upstream wave field separation；It specifically, can be by the object function It is multiplied by specified parameter, and using result of calculation as the pretreated object function；Wherein, using described in following formula characterization Specified parameter：

A=(1- β γ)⁴

For example, enabling

Formula 13 and formula 15 are substituted into formula 14, then had：

It enables

Wherein,Indicate the pretreated object function；

(2) the pretreated object function is calculated respectively to down-going wave fields separation and upstream wave field separation Single order local derviation；Specifically, 17 both sides of formula differentiate to down-going wave fields separation β, and it is 0 to enable derivative, then has：

2f₂β+2f₄βγ²+2f₅γ-2f₇γ²-2f₈-4f₉βγ+2f₁₀γ=0 (formula 18)

17 both sides of formula differentiate to upstream wave field separation γ, and it is 0 to enable derivative, then has：

2f₁γ+2f₄β²γ+2f₅β-2f₆-4f₇βγ-2f₉β²+2f₁₀β=0 (formula 19)

Wherein, the down-going wave fields separation is used to characterize in the initial uplink wavefield data comprising down-going wave fields Energy coefficient, the upstream wave field separation is for characterizing the energy for including upstream wave field in the initial downlink wavefield data Coefficient；

(3) the single order local derviation based on calculating is as a result, establish the pass of down-going wave fields separation and upstream wave field separation Connection relationship；Specifically, the incidence relation can be established according to formula 18 and formula 19；For example, establishing institute using following formula State incidence relation：

Wherein, β indicates down-going wave fields separation, includes down-going wave fields for characterizing in the initial uplink wavefield data Energy coefficient, γ indicate upstream wave field separation, for characterize in the initial downlink wavefield data include upstream wave field Energy coefficient；

(4) the first upstream wave field separation is set, and according to the first upstream wave field separation and the association Relationship determines the target downlink wave field separation coefficient and the target upstream wave field separation.

In the present embodiment, the first upstream wave field separation γ can be set₁It is 0.

In the present embodiment, according to the first upstream wave field separation and the incidence relation, the mesh is determined Down-going wave fields separation and the target upstream wave field separation are marked, can specifically include：

1) according to the first upstream wave field separation and the incidence relation, the first down-going wave fields segregative line is determined Number；It specifically, can be by the first upstream wave field separation γ₁Formula 20 is substituted into, first down going wave is calculated Field separation β₁；

2) according to the first down-going wave fields separation and the incidence relation, the second upstream wave field segregative line is determined Number；It specifically, can be by the first down-going wave fields separation β₁Formula 21 is substituted into, second upstream wave field is calculated Separation γ₂；

3) according to the second upstream wave field separation and the incidence relation, the second down-going wave fields segregative line is determined Number；It specifically, can be by the first upstream wave field separation γ₂Formula 20 is substituted into, first down going wave is calculated Field separation β₂；

4) when the absolute value of the difference of the first upstream wave field separation and the second upstream wave field separation is small In the absolute of predetermined coefficient threshold value and the difference of the first down-going wave fields separation and the second down-going wave fields separation When value is less than the predetermined coefficient threshold value, by the first upstream wave field separation and the second upstream wave field separation The sum of half, as the target upstream wave field separation, and by the first down-going wave fields separation with it is described The half of the sum of second down-going wave fields separation, as the target downlink wave field separation coefficient；

5) otherwise, using the second upstream wave field separation as the first new upstream wave field separation, and by institute The second down-going wave fields separation is stated as the first new down-going wave fields separation, and generates new the second upstream wave field separation Coefficient and the second new down-going wave fields separation, until the first new upstream wave field separation and described new second The absolute value of the difference of upstream wave field separation be less than predetermined coefficient threshold value and the first new down-going wave fields separation with The absolute value of the difference of the second new down-going wave fields separation is less than the predetermined coefficient threshold value, and by described new first The half of upstream wave field separation and the sum of the second new upstream wave field separation, as the target upstream wave field Separation, and by the first new down-going wave fields separation and new the sum of the second down-going wave fields separation Half, as the target downlink wave field separation coefficient.

Step S105：Based on the target wave field separation coefficient, the initial uplink wavefield data and the initial downlink Wavefield data determines target upstream wave field data and target downlink wavefield data.

In the present embodiment, the target wave field separation coefficient may include：Target downlink wave field separation coefficient and mesh Put on traveling-wave field separation；Wherein, the target downlink wave field separation coefficient is for characterizing the initial uplink wavefield data In include the target energy coefficients of down-going wave fields, the target upstream wave field separation is for characterizing the initial downlink wave field Include the target energy coefficient of upstream wave field in data.Following formula may be used and determine the target upstream wave field data：

The target downlink wavefield data is determined using following formula：

Wherein,Indicate the target downlink wavefield data.Namely by the target downlink wave field separation factor beta ' and The target upstream wave field separation γ ' substitutes into formula 5 and formula 6, and the target upstream wave field data and institute is calculated State target downlink wavefield data.So as to efficiently separate upstream wave field and down-going wave fields.

Before step S101, the embodiment of the present application can also include detector seismic data in the water to the target area It is pre-processed with land detector seismic data；Corresponding, detector seismic data is pre- place in water described in step S101 Detector seismic data in water after reason, the land detector seismic data are pretreated land wave detector earthquake number According to.Pretreatment described in the embodiment of the present application specifically includes：To detector seismic data in the water and land wave detector Shake data set label, define observation system, divide detector seismic data in the water and land detector seismic data From operations such as, denoising, filtering, velocity analysis and overlap-add procedures, it is important to note, however, that the above-mentioned pretreatment operation enumerated is only It is that in order to better illustrate the present invention, other pretreatment modes can also be used to be handled, itself is please not construed as limiting this.

After step S105, the embodiment of the present application can also include drawing target upstream wave field data and target downlink wave field Data section, and storage target upstream wave field data and target downlink wavefield data.

The land and water detector seismic data uplink and downlink wave field separation embodiment of the method utilizes identified traveling-wave field up and down The object function that correlation function is established, successive iteration inverting obtains target wave field separation coefficient, i.e., in initial uplink wavefield data Including including the target energy coefficient of upstream wave field, inverting in the target energy coefficient and initial downlink wavefield data of down-going wave fields Computational accuracy is higher；Finally, it is based on target wave field separation coefficient, detector seismic data and land wave detector in water may be implemented Upstream wave field and down-going wave fields efficiently separate in seismic data, reach and eliminate seawater multiple wave interference effect in seismic data, To improve seismic data signal-to-noise ratio and resolution ratio, and upstream wave field and down going wave number of fields are provided for latter earthquake data processing According to.Moreover, upper and lower traveling-wave field correlation function is the fourth order cumulant of land and water detector seismic data, and fourth-order cumulant measurer There is stronger noise immunity, so that the separating treatment of upstream wave field and down-going wave fields has stronger anti-noise ability.Meanwhile it using When the method for the embodiment of the present application carries out the separating treatment of upstream wave field and down-going wave fields, calculation amount is small, calculating speed is fast.

For the advantageous effect for illustrating the embodiment of the present application of removing, illustrate below in conjunction with the accompanying drawings：

Fig. 2 is the schematic diagram of horizontal layer speed and formation thickness model in the embodiment of the present application.Fig. 3 (a) is the application The diagrammatic cross-section of land detector seismic data in embodiment.Fig. 3 (b) is wave detector earthquake number in water in the embodiment of the present application According to diagrammatic cross-section.Fig. 4 (a) is isolated from Fig. 3 (a) and Fig. 3 (b) using the application method in the embodiment of the present application The diagrammatic cross-section of down-going wave fields data.Fig. 4 (b) is to use the application method from Fig. 3 (a) and Fig. 3 (b) in the embodiment of the present application In the diagrammatic cross-section of upstream wave field data isolated；Wherein, what Fig. 3 (a) and Fig. 3 (b) was characterized respectively is according in Fig. 2 Detector seismic data in the collected land detector seismic data of model and water.Speed thickness model in Fig. 1, i.e., it is horizontal Stratiform speed and formation thickness model, in Fig. 1, S indicate that shot point, R0 indicate that sea level, R1 indicate the interface of underground first (sea Bottom), R2 indicates that the interface of underground second, R3 indicate underground third interface；Between sea level R0 and the interface R1 of underground first The depth of water be 700 meters, corresponding horizontal layer speed is 1500 meter per seconds (m/s), the interface R1 of underground first and underground second Formation thickness between the R2 of interface is 800 meters, and corresponding horizontal layer speed is 2000m/s, the interface R2 of underground second and ground Formation thickness between lower third interface R3 is 700 meters, and corresponding horizontal layer speed is 3000m/s, underground third boundary The corresponding horizontal layer speed in the stratum below face R3 is 4500m/s, and geophone station sinking is in seabed, the i.e. interface of underground first R1, geophone station spacing are 25 meters, and shot point sinking depth is 7.5 meters, and shot point spacing is 25 meters.In Fig. 3 (a) and Fig. 3 (b), 1. It is interface (seabed) primary event (direct wave) (S-R1) of underground first, is 2. the interface primary event (S-R1- of underground second R2-R1), 3. it is first interface in underground sea primary event (S-R1-R0-R1), is 4. underground third interface primary event (S-R1-R2-R3-R2-R1), 5. it is the interface secondary reflection (S-R1-R2-R1-R2-R1) of underground second, is 6. underground second A interface sea primary event (S-R1-R2-R1-R0-R1) and third interface secondary reflection (S-R1-R2-R3-R2-R1), 7. It is first interface in underground sea secondary reflection (S-R1-R0-R1-R0-R1), is 8. the second contact surface two of underground third interface one time Secondary reflection (S-R1-R2-R3-R2-R1-R2-R1) is 9. underground third interface sea primary event (S-R1-R2-R3-R2- ) and third interface triple reflection (S-R1-R2-R3-R2-R3-R2-R3-R2-R1) R1-R0-R1.Such as Fig. 4 (a) and Fig. 4 (b) Shown, the down-going wave fields data of Fig. 4 (a) characterizations include：1. first interface (seabed) direct wave (S-R1), 3. first boundary Face sea primary event (S-R1-R0-R1), 6. in second interface sea primary event (S-R1-R2-R1-R0-R1), 7. First interface sea secondary reflection, (S-R1-R0-R1-R0-R1), 9. in third interface sea primary event (S-R1- R2-R3-R2-R1-R0-R1), be sea level multiple reflections, i.e. sea level multiple wave.The upstream wave field data of Fig. 4 (b) characterizations Include：1. first interface (seabed) primary event (S-R1), 2. second interface primary event (S-R1-R2-R1), 4. Three interface primary events (S-R1-R2-R3-R2-R1), 5. second interface secondary reflection (S-R1-R2-R1-R2-R1), 6. In third interface secondary reflection (S-R1-R2-R3-R2-R1,8. the second contact surface secondary reflection of third interface one time, (S-R1- R2-R3-R2-R1-R2-R1), 9. in the triple reflection of third interface, (S-R1-R2-R3-R2-R3-R2-R3-R2-R1), It is subsurface interface reflection, as significant wave.In this way, may indicate that the land and water detector data using the application or more traveling-wave field point From method, upstream wave field and down-going wave fields can be separated completely, and upstream wave field does not include remaining down-going wave fields, under Traveling-wave field does not include remaining upstream wave field.

Fig. 5 (a) is the diagrammatic cross-section of wave detector common-shot-gather data in water in the embodiment of the present application.Fig. 5 (b) is this Apply for the diagrammatic cross-section of land wave detector common-shot-gather data in embodiment.Fig. 5 (c) is in the embodiment of the present application using this The diagrammatic cross-section for the upstream wave field data that application method is separated from Fig. 5 (a) and Fig. 5 (b).Fig. 5 (d) is the application reality Apply the diagrammatic cross-section for the down-going wave fields data separated from Fig. 5 (a) and Fig. 5 (b) using the application method in example.Fig. 6 (a) be wave detector common midpoint stack data in water in the embodiment of the present application diagrammatic cross-section.Fig. 6 (b) is that the application is implemented The diagrammatic cross-section of land wave detector common midpoint stack data in example.Fig. 6 (c) is that the application side is used in the embodiment of the present application The diagrammatic cross-section for the upstream wave field data that method is separated from Fig. 6 (a) and Fig. 6 (b).Fig. 6 (d) is in the embodiment of the present application Using the diagrammatic cross-section for the down-going wave fields data that the application method is separated from Fig. 6 (a) and Fig. 6 (b).From Fig. 5 (a)~ As can be seen that using the land and water detector data uplink and downlink wave field separation side of the application in Fig. 5 (d) and Fig. 6 (a)~Fig. 6 (d) Upstream wave field data and down-going wave fields data can be effectively separated in method, and eliminate the multiple wave interference of water layer generation, carry The high signal-to-noise ratio of marine seismic data.

Fig. 7 is the composed structure signal of the application land and water detector seismic data uplink and downlink wave field separation device embodiment Figure.The land and water detector seismic data uplink and downlink wave field separation device provide in the water of target area detector seismic data and Land detector seismic data.As shown in fig. 7, the land and water detector seismic data uplink and downlink wave field separation device can wrap It includes：Initial traveling-wave field determining module 100, correlation function determining module 200, object function determining module 300, wave field separation up and down Traveling-wave field determining module 500 above and below coefficient determination module 400 and target.

The initial traveling-wave field determining module 100 up and down, can be used for based on detector seismic data and institute in the water Land detector seismic data is stated, determines initial uplink wavefield data and initial downlink wavefield data.

The correlation function determining module 200, can be used for based on the initial uplink wavefield data and it is described it is initial under Traveling wave field data determines traveling-wave field correlation function up and down.

The object function determining module 300 can be used for, according to the traveling-wave field correlation function up and down, establishing target letter Number.

The wave field separation coefficient determination module 400 can be used for being based on the object function, determine target wave field separation Coefficient.

Target traveling-wave field determining module 500 up and down, can be used for based on the target wave field separation coefficient, it is described just Beginning upstream wave field data and the initial downlink wavefield data, determine target upstream wave field data and target downlink wavefield data.

In the present embodiment, the initial traveling-wave field determining module 100 up and down can be used for wave detector in the water The sum of seismic data and the land detector seismic data as the initial uplink wavefield data, and will be examined in the water Wave device seismic data subtracts the land detector seismic data, and using result of calculation as the initial downlink wavefield data.

In the present embodiment, the traveling-wave field correlation function up and down may include：Upstream wave field auto-correlation function, downlink Wave field auto-correlation function, upstream wave field and down-going wave fields cross-correlation function and down-going wave fields and upstream wave field cross-correlation function； The correlation function determining module 200 can be used for determining the upstream wave field auto-correlation function using following formula：

In the present embodiment, the object function determining module 300 can be used for establishing the mesh using following formula Scalar functions：

In the present embodiment, the target wave field separation coefficient may include：Target downlink wave field separation coefficient and mesh Put on traveling-wave field separation；Wherein, the target downlink wave field separation coefficient is for characterizing the initial uplink wavefield data In include the target energy coefficients of down-going wave fields, the target upstream wave field separation is for characterizing the initial downlink wave field Include the target energy coefficient of upstream wave field in data；The wave field separation coefficient determination module 400 can be used for the mesh Scalar functions are pre-processed, and pretreated object function is obtained；Wherein, the pretreated object function and down-going wave fields Separation, upstream wave field separation are associated；The pretreated object function is calculated respectively to detach down-going wave fields The single order local derviation of coefficient and upstream wave field separation；Single order local derviation based on calculating is as a result, establish down-going wave fields separation With the incidence relation of upstream wave field separation；First upstream wave field separation is set, and according to first upstream wave field Separation and the incidence relation determine the target downlink wave field separation coefficient and the target upstream wave field segregative line Number.

In the present embodiment, the wave field separation coefficient determination module 400 can be used for according to first upgoing wave Field separation and the incidence relation determine the first down-going wave fields separation, and are detached according to first down-going wave fields Coefficient and the incidence relation, determine the second upstream wave field separation, and according to the second upstream wave field separation and The incidence relation determines the second down-going wave fields separation；When on the first upstream wave field separation and described second The absolute value of the difference of traveling-wave field separation is less than predetermined coefficient threshold value and the first down-going wave fields separation and described the When the absolute value of the difference of two down-going wave fields separations is less than the predetermined coefficient threshold value, by the first upstream wave field segregative line The half of number and the sum of the second upstream wave field separation, as the target upstream wave field separation, and by institute The half for stating the sum of the first down-going wave fields separation and the second down-going wave fields separation, as the target downlink wave Field separation；Otherwise, using the second upstream wave field separation as the first new upstream wave field separation, and by institute The second down-going wave fields separation is stated as the first new down-going wave fields separation, and generates new the second upstream wave field separation Coefficient and the second new down-going wave fields separation, until the first new upstream wave field separation and described new second The absolute value of the difference of upstream wave field separation be less than predetermined coefficient threshold value and the first new down-going wave fields separation with The absolute value of the difference of the second new down-going wave fields separation is less than the predetermined coefficient threshold value, and by described new first The half of upstream wave field separation and the sum of the second new upstream wave field separation, as the target upstream wave field Separation, and by the first new down-going wave fields separation and new the sum of the second down-going wave fields separation Half, as the target downlink wave field separation coefficient.

In the present embodiment, traveling-wave field determining module 500 can be used for determining using following formula the target up and down The target upstream wave field data：

The target downlink wavefield data is determined using following formula：

Wherein,Indicate the target downlink wavefield data.

The land and water detector seismic data uplink and downlink wave field separation device embodiment and the land and water wave detector earthquake number It is corresponding according to uplink and downlink wave field separation embodiment of the method, land and water detector seismic data uplink and downlink wave field separation method may be implemented Embodiment, and the technique effect of adquisitiones embodiment.

In the 1990s, the improvement of a technology can be distinguished clearly 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 with 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 to be 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, programmable logic device (Programmable Logic Device, PLD) (such as field programmable gate array (Field Programmable Gate Array, FPGA)) it is exactly such a integrated circuit, logic function determines device programming by user.By designer Voluntarily programming comes a digital display circuit " integrated " on a piece of PLD, designs and makes without asking chip maker Dedicated 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 there are many kind, 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 needing method flow slightly programming in logic and being programmed into integrated circuit with above-mentioned several hardware description languages In, so that it may to be readily available the hardware circuit for realizing the logical method flow.

It is also known in the art that other than 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 so that controller with logic gate, switch, application-specific integrated circuit, programmable The form of logic controller and embedded microcontroller etc. realizes identical function.Therefore this controller is considered one kind Hardware component, and the structure that the device for realizing various functions for including in it can also be considered as in hardware component.Or Even, can will be considered as realizing the device of various functions either the software module of implementation method can be Hardware Subdivision again Structure in part.

Device that above-described embodiment illustrates, module can specifically realize by computer chip or entity, or by having certain The product of function is planted to realize.

For convenience of description, it is divided into various modules when description apparatus above with function to describe respectively.Certainly, implementing this The function of each module is realized can in the same or multiple software and or hardware when application.

As seen through the above description of the embodiments, those skilled in the art can be understood that the application can It is realized by the mode of software plus required general hardware platform.Based on this understanding, the technical solution essence of the application On in other words the part that contributes to existing technology can be expressed in the form of software products, in a typical configuration In, computing device includes one or more processors (CPU), input/output interface, network interface and memory.The computer is soft Part product may include that some instructions are used so that a computer equipment (can be personal computer, server or network Equipment etc.) execute method described in certain parts of each embodiment of the application or embodiment.The computer software product can To store in memory, memory may include the volatile memory in computer-readable medium, random access memory (RAM) and/or the forms such as Nonvolatile memory, such as read-only memory (ROM) or flash memory (flash RAM).Memory is computer The example of readable medium.Computer-readable medium includes that permanent and non-permanent, removable and non-removable media can be by Any method or technique realizes information storage.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 memory (ROM), electrically erasable programmable read-only memory (EEPROM), fast flash memory bank or other memory techniques, CD-ROM are read-only Memory (CD-ROM), digital versatile disc (DVD) or other optical storages, magnetic tape cassette, tape magnetic disk storage or Other magnetic storage apparatus or any other non-transmission medium can be used for storage and can be accessed by a computing device information.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 in a progressive manner, identical similar portion between each embodiment Point just to refer each other, and each embodiment focuses on the differences from other embodiments.Especially for device reality For applying example, since it is substantially similar to the method embodiment, so description is fairly simple, related place 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 device, laptop device, multicomputer system, microprocessor-based system, set Top box, programmable consumer-elcetronics devices, network PC, minicomputer, mainframe computer including any of the above system or equipment Distributed computing environment etc..

The application can describe in the general context of computer-executable instructions executed by a computer, such as program Module.Usually, program module includes routines performing specific tasks or implementing specific abstract data types, program, object, group Part, data structure etc..The application can also be put into practice in a distributed computing environment, in these distributed computing environments, by Task is executed by the connected remote processing devices of communication network.In a distributed computing environment, program module can be with In the local and remote computer storage media including storage device.

Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application there are many deformation and Variation is without departing from spirit herein, it is desirable to which the attached claims include these deformations and change without departing from the application's Spirit.

Claims

1. a kind of land and water detector seismic data uplink and downlink wave field separation method, which is characterized in that be provided with the water of target area Middle detector seismic data and land detector seismic data, including：

Based on detector seismic data in the water and the land detector seismic data, determine initial uplink wavefield data and Initial downlink wavefield data；

Based on the initial uplink wavefield data and the initial downlink wavefield data, traveling-wave field correlation function up and down is determined；

2. according to the method described in claim 1, it is characterized in that, described based on detector seismic data in the water and described Land detector seismic data determines initial uplink wavefield data and initial downlink wavefield data, including：

By the sum of detector seismic data in the water and the land detector seismic data, as the initial uplink wave field Data；

Detector seismic data in the water is subtracted into the land detector seismic data, and using result of calculation as described first Beginning down-going wave fields data.

3. according to the method described in claim 1, it is characterized in that, the traveling-wave field correlation function up and down includes：Upstream wave field Auto-correlation function, down-going wave fields auto-correlation function, upstream wave field and down-going wave fields cross-correlation function and down-going wave fields and uplink Wave field cross-correlation function；The upstream wave field auto-correlation function is determined using following formula：

Wherein, a (k) indicates the upstream wave field auto-correlation function, U_i,jIndicate the initial uplink wavefield data, window when i is indicated The serial number in parameter road, i=1, window parameter total road number when 2 ..., M, M are indicated, the sequence of window parameter temporal sampling point when j is indicated Number, j=1, window parameter temporal number of samples when 2 ..., N, N are indicated, k expression correlation function delay sampling point serial numbers, k=0,1, 2 ..., L-1, L indicate correlation function length；

4. according to the method described in claim 1, it is characterized in that, establishing the object function using following formula：

Wherein, Q₂Indicate that the object function, a (k) indicate that the upstream wave field auto-correlation function, b (k) indicate the down going wave Auto-correlation function, c (k) indicate the upstream wave field and down-going wave fields cross-correlation function, d (k) indicate the down-going wave fields with Upstream wave field cross-correlation function, k indicate that correlation function postpones sampling point serial number, k=0, and 1,2 ..., L-1, L indicate correlation function Length；β indicates that the energy coefficient for including down-going wave fields in the initial uplink wavefield data, γ indicate the initial downlink wave field Include the energy coefficient of upstream wave field in data.

5. according to the method described in claim 1, it is characterized in that, the target wave field separation coefficient includes：Target downlink wave Field separation and target upstream wave field separation；Wherein, the target downlink wave field separation coefficient is described first for characterizing Include the target energy coefficient of down-going wave fields in beginning upstream wave field data, the target upstream wave field separation is for characterizing institute State the target energy coefficient for including upstream wave field in initial downlink wavefield data；It is described to be based on the object function, determine target Wave field separation coefficient, including：

The object function is pre-processed, pretreated object function is obtained；Wherein, the pretreated target letter Number is associated with down-going wave fields separation, upstream wave field separation；

The pretreated object function is calculated respectively to the single order of down-going wave fields separation and upstream wave field separation Local derviation；

Single order local derviation based on calculating is as a result, establish the incidence relation of down-going wave fields separation and upstream wave field separation；

First upstream wave field separation is set, and according to the first upstream wave field separation and the incidence relation, really The fixed target downlink wave field separation coefficient and the target upstream wave field separation.

6. according to the method described in claim 5, it is characterized in that, described pre-process the object function, including：

The object function is multiplied by specified parameter, and using result of calculation as the pretreated object function；Wherein, it adopts The specified parameter is characterized with following formula：

A=(1- β γ)⁴

Wherein, A indicates that the specified parameter, β indicate the energy coefficient for including down-going wave fields in the initial uplink wavefield data, γ indicates the energy coefficient for including upstream wave field in the initial downlink wavefield data.

7. according to the method described in claim 5, it is characterized in that, the traveling-wave field correlation function up and down includes：Upstream wave field Auto-correlation function, down-going wave fields auto-correlation function, upstream wave field and down-going wave fields cross-correlation function and down-going wave fields and uplink Wave field cross-correlation function；The incidence relation is established using following formula：

Wherein, β indicates down-going wave fields separation, for characterizing the energy for including down-going wave fields in the initial uplink wavefield data Coefficient of discharge, γ indicates upstream wave field separation, for characterizing the energy for including upstream wave field in the initial downlink wavefield data Coefficient of discharge；A (k) indicates that the upstream wave field auto-correlation function, b (k) indicate that the down-going wave fields auto-correlation function, c (k) indicate The upstream wave field and down-going wave fields cross-correlation function, d (k) indicate the down-going wave fields and upstream wave field cross-correlation function, k tables Show that correlation function postpones sampling point serial number, k=0,1,2 ..., L-1, L indicate correlation function length.

8. according to the method described in claim 5, it is characterized in that, described according to the first upstream wave field separation and institute Incidence relation is stated, determines the target downlink wave field separation coefficient and the target upstream wave field separation, including：

According to the first upstream wave field separation and the incidence relation, the first down-going wave fields separation, and root are determined According to the first down-going wave fields separation and the incidence relation, the second upstream wave field separation is determined, and according to described Second upstream wave field separation and the incidence relation determine the second down-going wave fields separation；

It is preset when the absolute value of the difference of the first upstream wave field separation and the second upstream wave field separation is less than The absolute value of the difference of coefficient threshold and the first down-going wave fields separation and the second down-going wave fields separation is less than When the predetermined coefficient threshold value, by the sum of the first upstream wave field separation and the second upstream wave field separation Half, as the target upstream wave field separation, and will be under the first down-going wave fields separation and described second The half of the sum of traveling-wave field separation, as the target downlink wave field separation coefficient；Otherwise, by second upstream wave field Separation is as the first new upstream wave field separation, and using the second down-going wave fields separation as first newly Down-going wave fields separation, and the second new upstream wave field separation and the second new down-going wave fields separation are generated, directly Absolute value of the difference to the first new upstream wave field separation and the second new upstream wave field separation is less than The difference of predetermined coefficient threshold value and the first new down-going wave fields separation and the second new down-going wave fields separation Absolute value be less than the predetermined coefficient threshold value, and by the first new upstream wave field separation on described new second The half of the sum of traveling-wave field separation, as the target upstream wave field separation, and by the first new downlink The half of wave field separation coefficient and the sum of the second new down-going wave fields separation, as the target downlink wave field separation Coefficient.

9. according to the method described in claim 1, it is characterized in that, the target wave field separation coefficient includes：Target downlink wave Field separation and target upstream wave field separation；Wherein, the target downlink wave field separation coefficient is described first for characterizing Include the target energy coefficient of down-going wave fields in beginning upstream wave field data, the target upstream wave field separation is for characterizing institute State the target energy coefficient for including upstream wave field in initial downlink wavefield data；The target upgoing wave is determined using following formula Field data：

Wherein,Indicate the target upstream wave field data, U_i,jIndicate the initial uplink wavefield data, D_i,jDescribed in expression Initial downlink wavefield data, the serial number in window parameter road when i is indicated, i=1, window parameter total road number when 2 ..., M, M are indicated, j tables The serial number of window parameter temporal sampling point when showing, j=1, window parameter temporal number of samples when 2 ..., N, N are indicated, k indicate correlation function Postpone sampling point serial number, k=0,1,2 ..., L-1, L indicate correlation function length；β ' indicates the target downlink wave field separation Coefficient, γ ' indicate the target upstream wave field separation；

The target downlink wavefield data is determined using following formula：

Wherein,Indicate the target downlink wavefield data.

10. a kind of land and water detector seismic data uplink and downlink wave field separation device, which is characterized in that described device provides target area Detector seismic data and land detector seismic data, described device include in the water in domain：Initial traveling-wave field up and down determines mould Traveling-wave field determines up and down for block, correlation function determining module, object function determining module, wave field separation coefficient determination module and target Module；Wherein,

The initial traveling-wave field determining module up and down, for based on detector seismic data and the land wave detector in the water Seismic data determines initial uplink wavefield data and initial downlink wavefield data；

The correlation function determining module, for being based on the initial uplink wavefield data and the initial downlink wavefield data, Determine traveling-wave field correlation function up and down；

Traveling-wave field determining module above and below the target, for based on the target wave field separation coefficient, the initial uplink wave field Data and the initial downlink wavefield data, determine target upstream wave field data and target downlink wavefield data.