CN105929446A - Data processing method and device in all-waveform inversion - Google Patents

Abstract

The present invention provides a data processing method and device in all-waveform inversion. The method comprises a step of obtaining seismic data to be processed and finding the commutation gun point positions of all corresponding gun points with the gun point of single gun data as a geophone, a step of carrying out forward processing, recording the data received at the commutation gun point position in the forward processing process, and calculating and obtaining a Hessian adjustment factor, a step of carrying out translation related processing on a source wave field to obtain the autocorrelation bandpass matrix of a source wave field, a step of multiplying the autocorrelation bandpass matrix of the source wave field with the Hessian adjustment factor, summing the single gun Hessian matrix data corresponding to all guns in the seismic data to be processed, and obtaining the total Hessian matrix data in the all-waveform inversion. By using each embodiment in the invention, the Hessian matrix in the all-waveform inversion processing can be calculated, the obtained Hessian matrix is used for optimization, and the purposes of rapid convergence and obtaining a global optimization result are achieved. According to the invention, the calculation efficiency of the all-waveform inversion can be greatly improved.

Description

Data processing method in a kind of full waveform inversion and device

Technical field

The present invention relates to seismic data processing technology in oil exploration, particularly relate to the data processing method in a kind of full waveform inversion And device.

Background technology

Velocity modeling is the committed step in seismic data process, the quality of its rate pattern built up, and directly affects geological data The end result processed: the quality in seismic profile face.Velocity modeling method is divided into two classes: the first kind is traditional chasing after based on ray The modeling method of track (ray tracing), such as reflection tomographic inversion method (Reflection Tomography).Its advantage is meter Calculation amount little (even unit can complete), its shortcoming is the model meeting smoother built up, and lacks details, and this is to be chased after by ray The limit decision of track method, ray tracing itself is to assume that (High frequency assumption) is to wave equation at high frequency An approximate solution, therefore rate pattern is had certain seriality requirement, if rate pattern has a lot of details (the most obvious Border), then ray tracing can failure.

Because modeling method based on ray tracing has above restriction, simultaneously as the raising of computing power, the most another kind of Modeling method: modeling method based on wave equation develops rapidly.Being characterized in can be with direct solution wave equation, therefore There are the potentiality building up refined model.Full waveform inversion (Full Waveform Inversion) is exactly the most most representative one, Full waveform inversion is converted into an optimization problem seismic modeling, its object function be analog data and observation data between European Distance:

$\begin{array}{c}E=R^{T*}R\\ R=u-d\end{array}---\left(1\right)$

In formula (1), R is residual error, is forward simulation data u and the difference of observation data d.R, u, d are vectors, T* Represent to measuring conjugate transpose.E is object function, is a numeral.In formula (1), u is the function of model m, and d is normal Number.After setting up such a object function, the task of seismic modeling reforms into an optimization problem, it may be assumed that find a geology Model so that the difference of analog data and observation data is minimum.

The method solving such optimization problem in state of the art conventional is Newton method.The method first calculating target function The E single order local derviation (being designated as g, be a vector) to model m, then calculating target function E is (usual to the second order local derviation of model It is referred to as Hessian matrix, is designated as H), then optimum results should be equal to:

Am=H^-1g

But, in actual production, it being rarely employed the method, reason is that Hessian matrix is difficult to calculate, Inversion Calculation inefficiency. Therefore people have to take the second best, and solve by steepest descent method or conjugate gradient method, and these methods not only restrain relatively slow, and more hold Easily it is absorbed in local minimum, thus obtains error result.Therefore, needing one in prior art badly can quickly, effectively, reliably It is calculated Hession matrix to use modeling method based on wave equation easily.

Summary of the invention

Present invention aim at providing the data processing method in a kind of full waveform inversion and device, can calculate at full waveform inversion Hessian matrix in reason, the Hessian matrix of gained is used for optimizing, and reaches Fast Convergent and obtains the purpose of global optimization result. The present invention can be greatly improved the computational efficiency of full waveform inversion.

Data processing method and device in a kind of full waveform inversion that the application provides are achieved in that

A kind of data processing method in full waveform inversion, described method includes:

Obtain pending geological data, for the single big gun data in described pending geological data, find with described single big gun data Shot point as all shot points corresponding during geophone station to easy shot position；

Described pending geological data is done and just drills process, record described just drilling in processing procedure and connect at easy shot position described The data received, utilize the described data to receiving at easy shot point to be calculated extra large gloomy Dynamic gene；

The source wavefield of described pending geological data does translation relevant treatment obtain the auto-correlation band of source wavefield and lead to matrix；

The auto-correlation band of described source wavefield is led to Matrix Multiplication with the gloomy Dynamic gene in described sea, obtain the gloomy matrix data in single big gun sea, with And the gloomy matrix data in single big gun sea corresponding for all big guns in described pending geological data is added, obtain sea total in full waveform inversion Gloomy matrix data.

In preferred embodiment, described method also includes:

Utilize the gloomy matrix data in calculated described sea to build the seismic velocity model in full waveform inversion, utilize described structure Seismic velocity model processes seismic profile data.

In preferred embodiment, the gloomy Dynamic gene in described sea includes using following manner to calculate:

It is calculated scale factor, the scale factor of all seismic channels in single big gun is added and obtains extra large gloomy Dynamic gene.

In preferred embodiment, the relation computing formula of described scale factor is:

$Factor=\frac{\underset{\mathrm{\ω}}{\Σ}F\left\{{\mathrm{\ω}}^2{u}_i\right\}{\left(t\right)}^2}{\underset{t}{\Σ}f{\left(t\right)}^2}$

In above formula, F{ ω²u_iIt is ω²u_iFourier transform, be a time series, by time domain is recorded time Between sequence u be ω²Filtering obtain.

In preferred embodiment, the described source wavefield to described pending geological data does translation relevant treatment and obtains source wavefield Auto-correlation band lead to matrix and include that using following manner to be calculated auto-correlation band leads to matrix H:

$H={\left(\frac{\∂u}{\∂m}\right)}^{*}\left(\frac{\∂u}{\∂m}\right)$

M is the vector of a length of K,It it is the matrix of N row K row；Employing following manner is calculated

Su=f

$\frac{\∂S}{\∂m}u+S\frac{\∂u}{\∂m}=0$

$\frac{\∂u}{\∂m}=-S^{-1}\frac{\∂S}{\∂m}u$

It is a three-dimensional tensor,It is a matrix,For non-zero item, value isSu=f is fluctuation side Journey.

In preferred embodiment, following manner is used to calculate matrixThe element of the i-th row kth row

ω²u_iIt is placed on i-th geophone station position u as wavelet_iOn do and just drill, take m_kThe response at place is multiplied byObtain matrix The element of the i-th row kth row

A kind of data processing equipment in full waveform inversion, described device includes:

Easy shot point is confirmed module, is used for obtaining pending geological data, for the single big gun data in described pending geological data, Find using the shot point of described single big gun data as all shot points corresponding during geophone station to easy shot position；

The gloomy factor computing module in sea, just drills process for doing described pending geological data, records described just drilling in processing procedure In the described data to receiving at easy shot position, the described data to receiving at easy shot point are utilized to be calculated extra large gloomy adjustment The factor；

Matrix calculus module, obtains source wavefield for the source wavefield of described pending geological data does translation relevant treatment Auto-correlation band leads to matrix；

The gloomy matrix calculus module in single big gun sea, for the auto-correlation band of described source wavefield is led to Matrix Multiplication with the gloomy adjustment in described sea because of Son, obtains the gloomy matrix data in single big gun sea；

Result output module, for the gloomy matrix data in single big gun sea corresponding for all big guns in described pending geological data is added, The gloomy matrix data in sea total in full waveform inversion.

In preferred embodiment, described device also includes:

Model construction module, for utilizing the gloomy matrix data in calculated described sea to build the seismic velocity mould in full waveform inversion Type, utilizes the seismic velocity model of described structure to process seismic profile data.

In preferred embodiment, described sea gloomy factor computing module includes:

Scale factor calculation module, may be used for using following formula to be calculated scale factor:

$Factor=\frac{\underset{\mathrm{\ω}}{\Σ}F\left\{{\mathrm{\ω}}^2{u}_i\right\}{\left(t\right)}^2}{\underset{t}{\Σ}f{\left(t\right)}^2}$

In above formula, F{ ω²u_iIt is ω²u_iFourier transform, be a time series, by time domain is recorded time Between sequence u be ω²Filtering obtain.

In preferred embodiment, the described source wavefield to described pending geological data does translation relevant treatment and obtains source wavefield Auto-correlation band lead to matrix and include that using following manner to be calculated auto-correlation band leads to matrix H:

$H={\left(\frac{\∂u}{\∂m}\right)}^{*}\left(\frac{\∂u}{\∂m}\right)$

M is the vector of a length of K,It it is the matrix of N row K row；Employing following manner is calculated

Su=f

$\frac{\∂S}{\∂m}u+S\frac{\∂u}{\∂m}=0$

$\frac{\∂u}{\∂m}=-S^{-1}\frac{\∂S}{\∂m}u$

It is a three-dimensional tensor,It is a matrix,For non-zero item, value isSu=f is fluctuation side Journey.

In preferred embodiment, following manner is used to calculate matrixThe element of the i-th row kth row

ω²u_iIt is placed on i-th geophone station position u as wavelet_iOn do and just drill, take m_kThe response at place is multiplied byObtain matrix The element of the i-th row kth row

Data processing method in a kind of full waveform inversion that the application provides and device, (the most right based on " to easy shot point " Easily focus) method, for every big gun data, can first find its " to easy shot point " position, and record " to easily during just drilling Shot point " data of position.The data obtained with sampling calculate a Dynamic gene, and the wave field of shot point self does spatial domain simultaneously Staggered relevant, obtain a logical matrix of band.This band leads to Matrix Multiplication with the Dynamic gene obtained before, it is simply that the Hessian of single big gun Matrix.The Hessian matrix of all big guns is added, and just obtains Hessian matrix total in geological data full waveform inversion.Gained Hessian matrix is used for optimizing, and reaches Fast Convergent and obtains the purpose of global optimization result, is greatly improved the calculating of full waveform inversion Efficiency.

Accompanying drawing explanation

In order to be illustrated more clearly that the embodiment of the present application or technical scheme of the prior art, below will be to embodiment or prior art In description, the required accompanying drawing used is briefly described, it should be apparent that, the accompanying drawing in describing below is only to remember in the application Some embodiments carried, for those of ordinary skill in the art, on the premise of not paying creative work, it is also possible to Other accompanying drawing is obtained according to these accompanying drawings.

Fig. 1 is the method flow schematic diagram of a kind of embodiment of data processing method in herein described a kind of full waveform inversion；

Fig. 2 is the method flow schematic diagram of the data processing method another kind embodiment in herein described a kind of full waveform inversion；

Fig. 3 is the modular structure schematic diagram of a kind of embodiment of data processing equipment in herein described a kind of full waveform inversion；

Fig. 4 is the modular structure schematic diagram of the data processing equipment another kind embodiment in herein described a kind of full waveform inversion.

Detailed description of the invention

For the technical scheme making those skilled in the art be more fully understood that in the application, below in conjunction with in the embodiment of the present application Accompanying drawing, the technical scheme in the embodiment of the present application is clearly and completely described, it is clear that described embodiment is only It is some embodiments of the present application rather than whole embodiments.Based on the embodiment in the application, ordinary skill people The every other embodiment that member is obtained under not making creative work premise, all should belong to the scope of the application protection.

Fig. 1 is the method flow diagram of one embodiment of data processing method in herein described a kind of full waveform inversion.Although this Application provides such as following embodiment or method operating procedure shown in the drawings or apparatus structure, but based on conventional or without creating The work of property can include more or less operating procedure or modular structure in described method or apparatus.In logicality not Existing in necessary causal step or structure, it is real that the execution sequence of these steps or the modular structure of device are not limited to the application Execute execution sequence or modular structure that example provides.Described method or the device in practice of modular structure or end product are applied Time, can connect according to embodiment or method shown in the drawings or modular structure carry out order perform or executed in parallel (such as Parallel processor or the environment of multiple threads).

The most as described in Figure 1, the data processing method in described a kind of full waveform inversion may include that

S1: obtain pending geological data, for the single big gun data in described pending geological data, finds with described single big gun number According to shot point as all shot points corresponding during geophone station to easy shot position.

Common geological data can include the data record of multiple shot point.Seismic data recording described in the embodiment of the present invention The data record of multiple shot point can be included, can also be for only including the data record of a shot point in other application scenarios, this Invent without limitation.After obtaining pending geological data, based on to easy principle, each in pending geological data Single big gun data of big gun, can first find using the shot point of described single big gun data as all shot points corresponding during geophone station to easy big gun Point position.Here shot position corresponding during data that the shot point data using single big gun receive as geophone station can be referred to as this shot point To easy shot point.

After obtaining pending geological data, can find for the single big gun data in described pending geological data with described single big gun number According to shot point as all shot points corresponding during geophone station to easy shot position.

S2: do described pending geological data and just drill process, records described just drilling in processing procedure described easy shot position The data that place receives, utilize the described data to receiving at easy shot point to be calculated extra large gloomy Dynamic gene.

Concrete, can use in time domain described pending geological data to be done and just drill process, and record to easy shot position The signal received at place.Described just drilling typically refers in geophysics, it is known that the property information of ball medium, as ( Seismic wave propagation velocities etc.), ask the process of when walking of seismic wave (parameter information such as seismic wave propagation time in the earth). It can be ui to the signal record of easy shot point by i-th that the present invention is just drilling in the implementation process of process.

Then can utilize acquisition and recording the data at easy shot point are calculated this list big gun the gloomy Dynamic gene in sea (or It is referred to as Hessian Dynamic gene).

Concrete, first we indicated that the object function to shape such as (1), its Hessian matrix is:

$H={\left(\frac{\∂u}{\∂m}\right)}^{*}\left(\frac{\∂u}{\∂m}\right)$

M is vector (assuming a length of K),It it is the matrix of N row K row.Problem is how to calculate now It is understood that u meets wave equation:

Su=f

So,

$\frac{\∂S}{\∂m}u+S\frac{\∂u}{\∂m}=0$

$\frac{\∂u}{\∂m}=-S^{-1}\frac{\∂S}{\∂m}u---\left(2\right)$

In above formula,It is a three-dimensional tensor (tensor),It it is a matrix.The non-zero item of only one of which, value is

For calculatingWe can use following manner to calculate matrixThe element of the i-th row kth row

ω²u_iIt is placed on i-th geophone station position u as wavelet_iOn do and just drill, take m_kThe response at place is multiplied byObtain matrix The element of the i-th row kth row

The most once can obtain all k m_kResponse.That is we can be obtained by by performing the most just drilling

ObtainAfter, we just can calculate

Specific practice is rightTranslating and be correlated with (shift and correlate), concrete can be the most rightIt is x, y The translation in direction, then andDoing is correlated with obtains relevant dilution.More than deriving is all in frequency domain, and derivation described above is all Being for single frequency, big gun to first suing for peace frequency, then should be sued for peace by correct result.The i.e. present invention can have two kinds to realize shape Formula, first can realize in frequency domain, now first should sue for peace frequency, then sue for peace big gun.Certainly, the present invention can also be Time domain realizes.

Sequence f to frequency domain, seeks f^*F's should be equal to the dot product in the corresponding sequence of time domain with result.It is to say, If the Fourier's series of f is F, then,

$\underset{\mathrm{\ω}}{\Σ}f{\left(\mathrm{\ω}\right)}^{*}f\left(\mathrm{\ω}\right)=\underset{t}{\Σ}F{\left(t\right)}^2$

That is in order to calculate in time domain:

$\underset{\mathrm{\ω}}{\Σ}{\left(\frac{\∂u_i}{\∂m}\right)}^{*}\left(\frac{\∂u_i}{\∂m}\right)$

We should be at u_iA ω is put on position²u_iWavelet, do and just drill, then calculate and wave field is relevant (the shift and of translation correlate).Acquired results is exactly the single geophone station contribution to Hessian diagonal angle item in single big gun.

Concrete considers how that realizing all geophone stations to all big guns sues for peace, it is contemplated that geophone station can also be shot point simultaneously.This Invention, in time domain full-wave type inverting, will calculate the illumination of shot point to each big gun.The one of the method for the invention is concrete real Executing in example, the gloomy Dynamic gene in described sea can include using following manner to calculate:

S201: be calculated scale factor, is added the scale factor of all seismic channels in single big gun and obtains extra large gloomy Dynamic gene.

Here scale factor is used to calculate Hessian Dynamic gene, is an intermediate object program.Concrete, described ratio The factor can use F{ ω²u_iDot product obtains divided by the dot product of wavelet.Source wavelet rather than ω can be used at this²u_i's Wavelet.Therebetween it is to there is certain proportionate relationship.In one embodiment of the present invention,

The relation computing formula of S202: this described scale factor can be:

$Factor=\frac{\underset{\mathrm{\ω}}{\Σ}F\left\{{\mathrm{\ω}}^2{u}_i\right\}{\left(t\right)}^2}{\underset{t}{\Σ}f{\left(t\right)}^2}---\left(3\right)$

In above formula (3), F{ ω²u_iIt is ω²u_iFourier transform, be a time series, can be by time domain record The u got off is ω²Filtering obtain, concrete implementation method can include first this time series of u being done Fourier transform, Obtain a frequency sequence, to each item in frequency sequence, be multiplied by respective frequencies square, the then frequency sequence to gained Do inverse fourier transform, obtain a time series.It it is the dot product of wavelet.Scale factor is equal to F{ ω²u_iDot product Dot product divided by wavelet.Further with regards to u_iCalculating also, in principle, u_i(claim from other shot points when doing geophone station equal to this shot point For to easy shot point) receive just drill signal, u can be calculated based on this_i.Embodiments of the invention can utilize easily Principle is u_iIt is considered to do from shot point just to drill, at the signal that easy shot point is received.

In concrete implementation process, for each big gun, it can be first found to do all shot positions corresponding during geophone station, I Referred to as " to easy shot point ".Then when doing time domain and just drilling, the record signal at " to easy shot point " place, it is designated as ui.. Simultaneously, that can add up source wavefield square obtains focus illumination.After just drilling end, per pass is accepted at easy shot point To signal be ω²Filtering.Scale factor is calculated with formula (3) before.The scale factor in all roads is added, and is somebody's turn to do The gloomy Dynamic gene in sea (Hessian Dynamic gene) of single big gun data.

Described pending geological data is done and just drills process, record described just drilling in processing procedure and connect at easy shot position described The data received, utilize the described data to receiving at easy shot point to be calculated extra large gloomy Dynamic gene

S3: the source wavefield of described pending geological data is done translation relevant treatment and obtains the auto-correlation band of source wavefield and lead to square Battle array.

As it was previously stated, source wavefield can be done the staggered relevant place self doing spatial domain of translation relevant treatment, i.e. shot point wave field Reason, can obtain a logical matrix of band, and this band leads to the auto-correlation band that matrix is source wavefield and leads to matrix.

S4: the auto-correlation band of described source wavefield leads to Matrix Multiplication with the gloomy Dynamic gene in described sea, obtains the gloomy matrix data in single big gun sea, And the gloomy matrix data in single big gun sea corresponding for all big guns in described pending geological data is added, obtain in full waveform inversion total The gloomy matrix data in sea.

Then can lead to Matrix Multiplication with the gloomy Dynamic gene in described sea with described auto-correlation band, obtain the gloomy square in single big gun sea of these list big gun data Battle array data, the Hessian of i.e. single big gun adjusts matrix.It is calculated the list of all big guns in pending geological data according to the method described above The gloomy matrix data in big gun sea, is added the gloomy matrix data in described all of single big gun sea, can obtain the gloomy square in sea total in full waveform inversion Battle array data.

Data processing method in full waveform inversion of the present invention provides a kind of full waveform inversion and calculates Hessian matrix Method, the method utilizing each embodiment to provide can be calculated Hessian matrix, the Hessian of gained effectively, reliably Matrix can accelerate the convergence rate that full waveform inversion optimizes.Utilize the Hessian matrix of gained, can be excellent with Gauss-Newton Change method carries out full waveform inversion, restrains fast than conventional steepest descent method and conjugate gradient method.

Certainly, after being calculated Hessian matrix, full waveform inversion method can be used to build seismic velocity model, such structure The model built out has higher precision, establishes good basis for seismic profile analysis.Fig. 2 is herein described a kind of Full wave shape The method flow schematic diagram of the data processing method another kind embodiment in inverting, as in figure 2 it is shown, the method for the invention is another In a kind of embodiment, described method can also include:

S5: utilize the gloomy matrix data in calculated described sea to build the seismic velocity model in full waveform inversion, utilize described structure The seismic velocity model built processes seismic profile data.

Based on method described above, the application also provides for the data processing equipment in a kind of full waveform inversion.Fig. 3 is the application institute State the modular structure schematic diagram of a kind of embodiment of the data processing equipment in a kind of full waveform inversion, concrete, as it is shown on figure 3, Described device may include that

Easy shot point is confirmed module 101, may be used for obtaining pending geological data, in described pending geological data Single big gun data, find using the shot point of described single big gun data as all shot points corresponding during geophone station to easy shot position；

The gloomy factor computing module 102 in sea, may be used for doing described pending geological data just drilling process, records and described is just drilling place In the described data to receiving at easy shot position during reason, the described data to receiving at easy shot point are utilized to be calculated The gloomy Dynamic gene in sea；

Matrix calculus module 103, may be used for that the source wavefield of described pending geological data does translation relevant treatment and is shaken The auto-correlation band of source wave field leads to matrix；

The gloomy matrix calculus module 104 in single big gun sea, may be used for the auto-correlation band of described source wavefield is led to Matrix Multiplication with described Hai Sen Dynamic gene, obtains the gloomy matrix data in single big gun sea；

Result output module 105, may be used for the gloomy matrix data in single big gun sea corresponding for all big guns in described pending geological data It is added, obtains the gloomy matrix data in sea total in full waveform inversion.

Fig. 4 is the modular structure schematic diagram of a kind of embodiment of data processing equipment in herein described a kind of full waveform inversion, separately In a kind of embodiment, described device can also include:

Model construction module 106, may be used for the ground utilizing the gloomy matrix data in calculated described sea to build in full waveform inversion Shake rate pattern, utilizes the seismic velocity model of described structure to process seismic profile data.

In another kind of embodiment, described sea gloomy factor computing module includes:

Scale factor calculation module, may be used for using following formula to be calculated scale factor:

$Factor=\frac{\underset{\mathrm{\ω}}{\Σ}F\left\{{\mathrm{\ω}}^2{u}_i\right\}{\left(t\right)}^2}{\underset{t}{\Σ}f{\left(t\right)}^2}$

In above formula, F{ ω²u_iIt is ω²u_iFourier transform, be a time series, by time domain is recorded time Between sequence u be ω²Filtering obtain.

In another kind of embodiment, the described source wavefield to described pending geological data does translation relevant treatment and obtains source wavefield Auto-correlation band lead to matrix and include that using following manner to be calculated auto-correlation band leads to matrix H:

$H={\left(\frac{\∂u}{\∂m}\right)}^{*}\left(\frac{\∂u}{\∂m}\right)$

M is the vector of a length of K,It it is the matrix of N row K row；Employing following manner is calculated

Su=f

$\frac{\∂S}{\∂m}u+S\frac{\∂u}{\∂m}=0$

$\frac{\∂u}{\∂m}=-S^{-1}\frac{\∂S}{\∂m}u$

It is a three-dimensional tensor,It is a matrix,For non-zero item, value isSu=f is fluctuation side Journey.

In a kind of embodiment, following manner can be used to calculate matrixThe element of the i-th row kth row

ω²u_iIt is placed on i-th geophone station position u as wavelet_iOn do and just drill, take m_kThe response at place is multiplied byObtain matrix The element of the i-th row kth row

Data processing method in a kind of full waveform inversion that the application provides and device, (the most right based on " to easy shot point " Easily focus) method, for every big gun data, can first find its " to easy shot point " position, and record " to easily during just drilling Shot point " data of position.The data obtained with sampling calculate a Dynamic gene, and the wave field of shot point self does spatial domain simultaneously Staggered relevant, obtain a logical matrix of band.This band leads to Matrix Multiplication with the Dynamic gene obtained before, it is simply that the Hessian of single big gun Matrix.The Hessian matrix of all big guns is added, and just obtains Hessian matrix total in geological data full waveform inversion.Gained Hessian matrix is used for optimizing, and reaches Fast Convergent and obtains the purpose of global optimization result.Utilize the Hessian matrix of gained, Full waveform inversion can be carried out with Gauss-Newton optimization method, restrain fast than conventional steepest descent method and conjugate gradient method, It is greatly improved the computational efficiency of full waveform inversion.

Although teachings herein being mentioned extra large gloomy matrix, just drilling, translate relevant or the like description, but, the application not office It is limited to must be to comply fully with the situation described by canonical algorithm or embodiment.On the basis of some algorithm or embodiment description slightly Amended embodiment can also carry out above-described embodiment identical, equivalent or close or deformation after foreseeable implementation result. Certainly, even if not using aforesaid way, as long as the data meeting the application the various embodiments described above process, information is mutual and information is sentenced Disconnected feedback system, still can realize identical application, not repeat them here.

Although this application provides the method operating procedure as described in embodiment or flow chart, but based on routine or without creativeness Means can include more or less operating procedure.The sequence of steps enumerated in embodiment is only numerous step execution sequences In a kind of mode, do not represent unique execution sequence.When device in practice or client production perform, can be according to reality Execute example or method shown in the drawings order performs or executed in parallel (environment of such as parallel processor or multiple threads).

Device that above-described embodiment illustrates or module, specifically can be realized by computer chip or entity, or by having certain merit The product of energy realizes.For convenience of description, it is divided into various module to be respectively described with function when describing apparatus above.Certainly, The function of each module can be realized in same or multiple softwares and/or hardware when implementing the application, it is also possible to will realize same The module of one function is realized by the combination of multiple submodules or subelement.

It is also known in the art that in addition to realizing controller in pure computer readable program code mode, the most permissible Make controller with gate, switch, special IC, FPGA control by method step carries out programming in logic The form of device processed and embedding microcontroller etc. realizes identical function.The most this controller is considered a kind of Hardware Subdivision Part, and its inside is included can also be considered as the structure in hardware component for the device realizing various function.Or even, In can being considered as the device being used for realizing various function not only can being the software module of implementation method but also can being hardware component Structure.

The application can be described in the general context of computer executable instructions, such as program module. Usually, program module include perform particular task or realize the routine of particular abstract data type, program, object, assembly, Data structure, class etc..The application can also be put into practice in a distributed computing environment, in these distributed computing environment, by The remote processing devices connected by communication network performs task.In a distributed computing environment, program module can position In the local and remote computer-readable storage medium including storage device.

As seen through the above description of the embodiments, those skilled in the art it can be understood that to the application can be by soft Part adds the mode of required general hardware platform and realizes.Based on such understanding, the technical scheme of the application is the most in other words The part contributing prior art can embody with the form of software product, and this computer software product can be stored in In storage medium, such as ROM/RAM, magnetic disc, CD etc., use so that a computer equipment is (permissible including some instructions Be personal computer, mobile terminal, server, or the network equipment etc.) perform each embodiment of the application or embodiment Method described in some part.

Each embodiment in this specification uses the mode gone forward one by one to describe, and between each embodiment, same or analogous part is mutual Seeing, what each embodiment stressed is the difference with other embodiments.The application can be used for numerous general or In special computing system environments or configuration.Such as: personal computer, server computer, handheld device or portable set Standby, laptop device, multicomputer system, system based on microprocessor, set top box, programmable electronic equipment, network PC, minicomputer, mainframe computer, the distributed computing environment including any of the above system or equipment etc..

Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application have many deformation and a change and Without departing from spirit herein, it is desirable to appended claim includes that these deformation and change are without deviating from spirit herein.

Claims (11)

1. the data processing method in a full waveform inversion, it is characterised in that described method includes:

Obtain pending geological data, for the single big gun data in described pending geological data, find with described single big gun data Shot point as all shot points corresponding during geophone station to easy shot position；

Described pending geological data is done and just drills process, record described just drilling in processing procedure and connect at easy shot position described The data received, utilize the described data to receiving at easy shot point to be calculated extra large gloomy Dynamic gene；

The source wavefield of described pending geological data does translation relevant treatment obtain the auto-correlation band of source wavefield and lead to matrix；

The auto-correlation band of described source wavefield is led to Matrix Multiplication with the gloomy Dynamic gene in described sea, obtain the gloomy matrix data in single big gun sea, with And the gloomy matrix data in single big gun sea corresponding for all big guns in described pending geological data is added, obtain sea total in full waveform inversion Gloomy matrix data.

Data processing method in a kind of full waveform inversion the most as claimed in claim 1, it is characterised in that described method is also Including:

Utilize the gloomy matrix data in calculated described sea to build the seismic velocity model in full waveform inversion, utilize described structure Seismic velocity model processes seismic profile data.

Data processing method in a kind of full waveform inversion the most as claimed in claim 1 or 2, it is characterised in that described sea Gloomy Dynamic gene includes using following manner to calculate:

It is calculated scale factor, the scale factor of all seismic channels in single big gun is added and obtains extra large gloomy Dynamic gene.

Data processing method in a kind of full waveform inversion the most as claimed in claim 3, it is characterised in that described ratio because of The relation computing formula of son is:

$Factor=\frac{\underset{\mathrm{\ω}}{\Σ}F\left\{{\mathrm{\ω}}^2{u}_i\right\}{\left(t\right)}^2}{\underset{t}{\Σ}f{\left(t\right)}^2}$

In above formula, F{ ω²u_iIt is ω²u_iFourier transform, be a time series, by time domain is recorded time Between sequence u be ω²Filtering obtain.

Data processing method in a kind of full waveform inversion the most as claimed in claim 3, it is characterised in that described to described The source wavefield of pending geological data does translation relevant treatment and obtains the auto-correlation band of source wavefield and lead to matrix and include using following Mode is calculated auto-correlation band and leads to matrix H:

$H={\left(\frac{\∂u}{\∂m}\right)}^{*}\left(\frac{\∂u}{\∂m}\right)$

M is the vector of a length of K,It it is the matrix of N row K row；Employing following manner is calculated

Su=f

$\frac{\∂S}{\∂m}u+S\frac{\∂u}{\∂m}=0$

$\frac{\∂u}{\∂m}=-S^{-1}\frac{\∂S}{\∂m}u$

It is a three-dimensional tensor,It is a matrix,For non-zero item, value isSu=f is wave equation.

Data processing method in a kind of full waveform inversion the most as claimed in claim 5, it is characterised in that use following side Formula calculates matrixThe element of the i-th row kth row

ω²u_iIt is placed on i-th geophone station position u as wavelet_iOn do and just drill, take m_kThe response at place is multiplied byObtain matrix The element of the i-th row kth row

7. the data processing equipment in a full waveform inversion, it is characterised in that described device includes:

Easy shot point is confirmed module, is used for obtaining pending geological data, for the single big gun data in described pending geological data, Find using the shot point of described single big gun data as all shot points corresponding during geophone station to easy shot position；

The gloomy factor computing module in sea, just drills process for doing described pending geological data, records described just drilling in processing procedure In the described data to receiving at easy shot position, the described data to receiving at easy shot point are utilized to be calculated extra large gloomy adjustment The factor；

Matrix calculus module, obtains source wavefield for the source wavefield of described pending geological data does translation relevant treatment Auto-correlation band leads to matrix；

The gloomy matrix calculus module in single big gun sea, for the auto-correlation band of described source wavefield is led to Matrix Multiplication with the gloomy adjustment in described sea because of Son, obtains the gloomy matrix data in single big gun sea；

Result output module, for the gloomy matrix data in single big gun sea corresponding for all big guns in described pending geological data is added, The gloomy matrix data in sea total in full waveform inversion.

Data processing equipment in a kind of full waveform inversion the most as claimed in claim 7, it is characterised in that described device is also Including:

Model construction module, for utilizing the gloomy matrix data in calculated described sea to build the seismic velocity mould in full waveform inversion Type, utilizes the seismic velocity model of described structure to process seismic profile data.

Data processing equipment in a kind of full waveform inversion the most as claimed in claim 7 or 8, it is characterised in that described sea Gloomy factor computing module includes:

Scale factor calculation module, is used for using following formula to be calculated scale factor:

$Factor=\frac{\underset{\mathrm{\ω}}{\Σ}F\left\{{\mathrm{\ω}}^2{u}_i\right\}{\left(t\right)}^2}{\underset{t}{\Σ}f{\left(t\right)}^2}$

In above formula, F{ ω²u_iIt is ω²u_iFourier transform, be a time series, by time domain is recorded time Between sequence u be ω²Filtering obtain.

Data processing equipment in a kind of full waveform inversion the most as claimed in claim 7 or 8, it is characterised in that described right The source wavefield of described pending geological data does translation relevant treatment and obtains the auto-correlation band of source wavefield and lead to matrix and include using Following manner is calculated auto-correlation band and leads to matrix H:

$H={\left(\frac{\∂u}{\∂m}\right)}^{*}\left(\frac{\∂u}{\∂m}\right)$

M is the vector of a length of K,It it is the matrix of N row K row；Employing following manner is calculated Su=f

$\frac{\∂S}{\∂m}u+S\frac{\∂u}{\∂m}=0$

$\frac{\∂u}{\∂m}=-S^{-1}\frac{\∂S}{\∂m}u$

It is a three-dimensional tensor,It is a matrix,For non-zero item, value isSu=f is fluctuation side Journey.

Data processing equipment in 11. a kind of full waveform inversion as claimed in claim 10, it is characterised in that use following Mode calculates matrixThe element of the i-th row kth row

ω²u_iIt is placed on i-th geophone station position u as wavelet_iOn do and just drill, take m_kThe response at place is multiplied byObtain matrix The element of the i-th row kth row