CN104990774A - Seismic data interpolation method - Google Patents
Seismic data interpolation method Download PDFInfo
- Publication number
- CN104990774A CN104990774A CN201510354493.3A CN201510354493A CN104990774A CN 104990774 A CN104990774 A CN 104990774A CN 201510354493 A CN201510354493 A CN 201510354493A CN 104990774 A CN104990774 A CN 104990774A
- Authority
- CN
- China
- Prior art keywords
- geological data
- interpolation
- matrix
- geological
- linear combination
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Abstract
The embodiment of the invention discloses a seismic data interpolation method. The method comprises the steps that first seismic data are obtained; based on the first seismic data, a distance matrix M is generated through a formula Mi,j=[B+d2(xi,xj)]a, and an inverse matrix M-1 of the distance matrix M is obtained; Fourier transforming is carried out on the first seismic data in the time direction, and second seismic data are obtained; based on the second seismic data and the inverse matrix M-1, a liner combination coefficient matrix C is obtained through a formula C=yM-1; based on the liner combination coefficient matrix C, an interpolation of the first seismic data in a frequency domain is obtained; Fourier inversion is carried out on the interpolation of the first seismic data in the frequency domain, and an interpolation of the first seismic data is obtained. According to the seismic data interpolation method, the stability of a seismic data interpolation result can be improved.
Description
Technical field
The application relates to technical field of geophysical exploration, particularly a kind of geological data interpolation method.
Background technology
In the gatherer process of geological data, due to the restriction of the objective condition such as construction environment and workload, usually there is the problem that on some direction, track pitch is larger in the seismologic record collected, makes space sampling frequency wretched insufficiency, can produce have a strong impact on migration imaging process etc.Address this problem when the most direct method is naturally and gathers in the wild and reduce spatial sampling interval by increasing survey line quantity; But the thing followed is increasing substantially of acquisition cost, and some work area is due to the restriction of actual terrestrial reference condition, is difficult to the space sampling frequency improving field work.
Geological data interpolation is under the prerequisite not increasing acquisition cost, reduces the Perfected process of track pitch.In prior art, conventional geological data difference approach, based on Fourier transform.Its detailed process is as follows:
1) geological data (is supposed to have N number of, is designated as u
1, u
2..., u
n) do Fourier transform, obtain one group of Fourier coefficient;
2) calculate the absolute value of each Fourier coefficient, and therefrom choose the Fourier coefficient of maximum absolute value, Fourier inversion is done to the Fourier coefficient chosen, obtain one group of fitting data and (be designated as a
1, a
2..., a
n);
3) formula is passed through
digital simulation coefficient c, wherein a
irepresent a fitting data, u
irepresent a geological data, the size of subscript i is from 0 to N, a
i *represent plural a
iconjugate complex number;
4) fitting coefficient c is utilized, by formula u
j=u
i– c a
iobtain the interpolation u of geological data
j;
5) by formula E=u
j *u
jcalculate the residual error E of each geological data interpolation, and described residual error E and predetermined threshold value are compared.When residual error E is less than described predetermined threshold value, terminate whole Interpolation Process.When residual error E is greater than described predetermined threshold value, by step 4) in the geological data interpolation that obtains and step 1) in geological data combine, generate new geological data, and step 1 is performed to new geological data) to step 5), till residual error E is less than described predetermined threshold value.
Realizing in the application's process, inventor finds that in prior art, at least there are the following problems:
Usually, Fourier transform has Gibbs' effect (Gibbs Phenomenon), and namely the fourier coefficient of one group of matching input data can generate the fitting function of a concuss.Because above-mentioned prior art calculates the interpolation of geological data based on fourier coefficient, therefore, the geological data interpolation that above-mentioned prior art obtains normally shake with instability.
Summary of the invention
The object of the embodiment of the present application is to provide a kind of geological data interpolation method, to improve the stability of geological data interpolation result.
For solving the problems of the technologies described above, the embodiment of the present application provides a kind of geological data interpolation method to be achieved in that
A kind of geological data interpolation method, comprising:
Obtain the first geological data;
Based on described first geological data, by formula M
i,j=[β+d
2(x
i, x
j)]
αgenerate distance matrix M, and ask for the inverse matrix M of described distance matrix M
-1, wherein, α and β is constant, x
iand x
jbe respectively two the first geological datas, d (x
i, x
j) be distance function, M
i,jthe numerical value of the capable j row of the i for distance matrix M;
Along time orientation, Fourier transform is done to described first geological data, obtains the second geological data;
Based on described second geological data and inverse matrix M
-1, by formula C=y M
-1ask for linear combination coefficient Matrix C, wherein, y is the second geological data;
Based on described linear combination coefficient Matrix C, obtain the interpolation of described first geological data in frequency field;
In the difference of frequency field, Fourier inversion is carried out to described first geological data, obtains the interpolation of the first geological data.
The technical scheme provided from above the embodiment of the present application, the geological data interpolation method of the embodiment of the present application, ask for linear combination coefficient Matrix C by distance matrix M, then by linear combination coefficient Matrix C, interpolation is carried out to geological data, obtain the interpolation result of geological data.Compared with prior art, the geological data interpolation method of the embodiment of the present application does not rely on fourier coefficient, thus not by the impact of Gibbs' effect, can ensure the stability of geological data interpolation result.
Accompanying drawing explanation
In order to be illustrated more clearly in the embodiment of the present application or technical scheme of the prior art, be briefly described to the accompanying drawing used required in embodiment or description of the prior art below, apparently, the accompanying drawing that the following describes is only some embodiments recorded in the application, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.
Fig. 1 is the process flow diagram of the embodiment of the present application geological data interpolation method.
Embodiment
Technical scheme in the application is understood better in order to make those skilled in the art person, below in conjunction with the accompanying drawing in the embodiment of the present application, technical scheme in the embodiment of the present application is clearly and completely described, obviously, described embodiment is only some embodiments of the present application, instead of whole embodiments.Based on the embodiment in the application, those of ordinary skill in the art are not making the every other embodiment obtained under creative work prerequisite, all should belong to the scope of the application's protection.
Above-mentioned interpolation method of the prior art carries out based on iteration.The stop condition of iteration is residual error little to a certain extent (being less than predetermined threshold value).Like this, often there is certain residual error in the interpolation result obtained by this interpolation method.Therefore, this interpolation method often can not ensure the correctness of geological data interpolation result.
Introduce an embodiment of the application's geological data interpolation method below.As shown in Figure 1, described embodiment comprises:
S101: obtain the first geological data.
Described first geological data is five dimension geological datas.Particularly, field excites Artificial Seismic Wave, is received and record seismic wave field by wave detector, obtains five dimension geological data x=x (a
m, a
n, b
r, b
s, t).Described five dimensions are generally a
m, a
n, b
r, b
sand t.Wherein, a
mand a
nrepresent the locus of two focal point m and n respectively, b
rand b
srepresentative and a respectively
mand a
nthe locus of two corresponding geophone station r and s, t represents observation time.The quantity of described first geological data is generally multiple.
Further, in described first geological data, observation time t is equally distributed.Namely at interval of the identical time, the first geological data is obtained.
S102: based on described first geological data, by formula M
i,j=[β+d
2(x
i, x
j)]
αgenerate distance matrix M, and ask for the inverse matrix M of described distance matrix M
-1.
Particularly, based on multiple first geological datas got in step S101, by formula M
i,j=[β+d
2(x
i, x
j)]
αdistance matrix M can be generated.Wherein, α and β is constant (being commonly defined as β=1, α=0.5), x
iand x
jrepresent two the first geological datas respectively, i and j represents the numbering of the first geological data, M
i,jrepresent the numerical value of the capable j row of the i of distance matrix M.Because the first geological data is focal point a
m, a
nand geophone station b
r, b
sfunction, therefore, can each first geological data be regarded as four dimensional vectors.So, function d (x
i, x
j) can be distance function, represent and ask for the first geological data x
iand x
jbetween distance (namely ask for two four-dimensional vector x
iand x
jbetween distance).
When after generation distance matrix M, matrix M of can adjusting the distance does low-rank decomposition, thus obtains inverse matrix M
-1.
S103: along time orientation, Fourier transform is done to described first geological data, obtains the second geological data.
Particularly, certain first geological data x can be chosen from the first geological data
i, then to this first geological data x
iafter doing Fourier transform along time orientation, can obtain and the first geological data x
ithe second corresponding geological data y
i.After all Fourier transform is done along time orientation to each first geological data, second geological data corresponding with the first geological data can be obtained.Described second geological data is the geological data of frequency field.
S104: based on described second geological data and inverse matrix M
-1, by formula C=y M
-1ask for linear combination coefficient Matrix C.
Particularly, certain second geological data y can be chosen from the second geological data
i, then according to formula c
i=y
im
-1obtain and this second geological data y
icorresponding linear combination coefficient c
i.When after the linear combination coefficient getting each second geological data, linear combination coefficient Matrix C can be generated.
S105: based on described linear combination coefficient Matrix C, obtains the interpolation of described first geological data in frequency field.
Particularly, the first geological data x can be chosen from the first geological data
k, then choose from linear combination coefficient Matrix C and the first geological data x
kcorresponding linear combination coefficient c
k, then pass through formula
obtain the first geological data x about the first geological data x
kat the interpolation L (x) of frequency field.
S106: in the difference of frequency field, Fourier inversion is carried out to described first geological data, obtains the interpolation of the first geological data.
Particularly, Fourier inversion can be carried out to described first geological data at the difference L (x) of frequency field, thus obtain the first geological data x about the first geological data x
kin the interpolation of time domain.
The geological data interpolation method of the embodiment of the present application, asks for linear combination coefficient Matrix C by distance matrix M, then carries out interpolation by linear combination coefficient Matrix C to geological data, obtains the interpolation result of geological data.Compared with prior art, on the one hand, the geological data interpolation method of the embodiment of the present application does not rely on fourier coefficient, thus not by the impact of Gibbs' effect, can ensure the stability of geological data interpolation result.On the other hand, the geological data interpolation method of the embodiment of the present application does not need iterative, thus does not have information dropout in Interpolation Process, can ensure the correctness of geological data interpolation result.
Although depict the application by embodiment, those of ordinary skill in the art know, the application has many distortion and change and do not depart from the spirit of the application, and the claim appended by wishing comprises these distortion and change and do not depart from the spirit of the application.
Claims (5)
1. a geological data interpolation method, is characterized in that, comprising:
Obtain the first geological data;
Based on described first geological data, by formula M
i,j=[β+d
2(x
i, x
j)]
αgenerate distance matrix M, and ask for the inverse matrix M of described distance matrix M
-1, wherein, α and β is constant, x
iand x
jbe respectively two the first geological datas, d (x
i, x
j) be distance function, M
i,jthe numerical value of the capable j row of the i for distance matrix M;
Along time orientation, Fourier transform is done to described first geological data, obtains the second geological data;
Based on described second geological data and inverse matrix M
-1, by formula C=y M
-1ask for linear combination coefficient Matrix C, wherein, y is the second geological data;
Based on described linear combination coefficient Matrix C, obtain the interpolation of described first geological data in frequency field;
In the difference of frequency field, Fourier inversion is carried out to described first geological data, obtains the interpolation of the first geological data.
2. the method for claim 1, is characterized in that, described based on described linear combination coefficient Matrix C, obtains described first geological data in the interpolation of frequency field, specifically comprises:
The first geological data x is chosen from the first geological data
k, obtain and described first geological data x from linear combination coefficient Matrix C
kcorresponding linear combination coefficient c
k, then pass through formula
obtain the first geological data x about the first geological data x
kat the interpolation L (x) of frequency field.
3. the method for claim 1, is characterized in that, described based on described second geological data and inverse matrix M
-1, by formula C=y M
-1ask for linear combination coefficient Matrix C, specifically comprise:
The second geological data y is chosen from the second geological data
i, then according to formula c
i=y
im
-1obtain and the second geological data y
icorresponding linear combination coefficient c
i.
4. the method for claim 1, is characterized in that, does Fourier transform, obtain the second geological data, specifically comprise described first geological data along time orientation:
The first geological data x is chosen from the first geological data
i, to this first geological data x
iafter doing Fourier transform along time orientation, obtain and described first geological data x
ithe second corresponding geological data y
i.
5. the method for claim 1, is characterized in that, described first geological data comprises the locus of two focal points, the locus of two geophone stations corresponding with these two focal points and observation time.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510354493.3A CN104990774B (en) | 2015-06-24 | 2015-06-24 | A kind of geological data interpolation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510354493.3A CN104990774B (en) | 2015-06-24 | 2015-06-24 | A kind of geological data interpolation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104990774A true CN104990774A (en) | 2015-10-21 |
| CN104990774B CN104990774B (en) | 2018-01-05 |
Family
ID=54302615
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510354493.3A Active CN104990774B (en) | 2015-06-24 | 2015-06-24 | A kind of geological data interpolation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104990774B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4594662A (en) * | 1982-11-12 | 1986-06-10 | Schlumberger Technology Corporation | Diffraction tomography systems and methods with fixed detector arrays |
| CN103901466A (en) * | 2012-12-28 | 2014-07-02 | 中国石油天然气集团公司 | Three-dimensional seismic data interpolation method |
| CN104360393A (en) * | 2014-10-20 | 2015-02-18 | 李闯 | Seismic data reconstruction method |
| CN104536044A (en) * | 2015-01-16 | 2015-04-22 | 中国石油大学(北京) | Interpolation and denoising method and system for seismic data |
-
2015
- 2015-06-24 CN CN201510354493.3A patent/CN104990774B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4594662A (en) * | 1982-11-12 | 1986-06-10 | Schlumberger Technology Corporation | Diffraction tomography systems and methods with fixed detector arrays |
| CN103901466A (en) * | 2012-12-28 | 2014-07-02 | 中国石油天然气集团公司 | Three-dimensional seismic data interpolation method |
| CN104360393A (en) * | 2014-10-20 | 2015-02-18 | 李闯 | Seismic data reconstruction method |
| CN104536044A (en) * | 2015-01-16 | 2015-04-22 | 中国石油大学(北京) | Interpolation and denoising method and system for seismic data |
Non-Patent Citations (2)
| Title |
|---|
| SIDE JIN: "5D seismic data regularization by a damped least-norm fourier inversion", 《GEOPHYSICS》 * |
| 苏欣等: "一种数据插值算法——在地震数据解释成果中的应用", 《承德石油高等专科学校学报》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104990774B (en) | 2018-01-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Trad | Five-dimensional interpolation: Recovering from acquisition constraints | |
| CN105277985A (en) | OVT-domain seismic data regularization method based on image processing | |
| Xue et al. | High-order sparse Radon transform for AVO-preserving data reconstruction | |
| Bakulin et al. | Application of supergrouping to enhance 3D prestack seismic data from a desert environment | |
| Zhang et al. | 2D and 3D prestack seismic data regularization using an accelerated sparse time-invariant Radon transform | |
| CN104597499A (en) | Adjacent shot interference suppressing method and device for independent simultaneous sweeping seismic data of seismic vibrators | |
| CN105676291A (en) | Multiple wave matching attenuation method based on optimized phase axis tracking | |
| CN107861156B (en) | The extracting method and device of diffracted wave | |
| Kong et al. | Deep prior-based seismic data interpolation via multi-res u-net | |
| CN110244353B (en) | Seismic data regularization method based on sparse norm optimization algorithm | |
| Liu et al. | Blind-source separation of seismic signals based on information maximization | |
| CN105510975A (en) | Method and device for improving signal-to-noise ratio of seismic data | |
| Xiang‐Bo et al. | Anisotropic Radon transform and its application to demultiple | |
| CN105301648A (en) | Method of acquiring common reflection surface stacking parameters | |
| Huang et al. | Least-squares seismic inversion with stochastic conjugate gradient method | |
| CN102903104B (en) | Subtractive clustering based rapid image segmentation method | |
| CN104990774A (en) | Seismic data interpolation method | |
| CN103513288B (en) | A kind of compensation direction filtering method of two-dimensional grid data | |
| CN111337973B (en) | Seismic data reconstruction method and system | |
| CN105301655B (en) | A kind of common imaging gather line noise minimizing technology and device | |
| CN110161565B (en) | Seismic data reconstruction method | |
| CN112200069A (en) | Tunnel filtering method and system combining time-frequency domain spectral subtraction and empirical mode decomposition | |
| WO2002065159A1 (en) | Method of seismic imaging using direct travel time computing | |
| Dmitriev et al. | Efficient four-dimensional supergrouping algorithm for enhancement of high-channel count seismic data | |
| Guruacharya et al. | Off-grid compressed sensing for radars with unknown transmit waveform |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |