CN110261912A - The interpolation and denoising method and system of a kind of seismic data - Google Patents

The interpolation and denoising method and system of a kind of seismic data Download PDF

Info

Publication number
CN110261912A
CN110261912A CN201910664806.3A CN201910664806A CN110261912A CN 110261912 A CN110261912 A CN 110261912A CN 201910664806 A CN201910664806 A CN 201910664806A CN 110261912 A CN110261912 A CN 110261912A
Authority
CN
China
Prior art keywords
seismic data
interpolation
denoising
sparse transformation
noise
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
Application number
CN201910664806.3A
Other languages
Chinese (zh)
Other versions
CN110261912B (en
Inventor
曹静杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hebei GEO University
Original Assignee
Hebei GEO University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hebei GEO University filed Critical Hebei GEO University
Priority to CN201910664806.3A priority Critical patent/CN110261912B/en
Publication of CN110261912A publication Critical patent/CN110261912A/en
Application granted granted Critical
Publication of CN110261912B publication Critical patent/CN110261912B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/30Noise handling
    • G01V2210/32Noise reduction
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/40Transforming data representation
    • G01V2210/48Other transforms

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention discloses a kind of interpolation of seismic data and denoising method and system.This method comprises: obtaining the seismic data of Noise;Regular parameter sequence is determined according to the seismic data of the Noise;Regular parameter in the regular parameter sequence is sequentially reduced;The seismic data of interpolation and denoising after calculating sparse transformation according to the regular parameter sequence;Seismic data by sparse transformation method, after interpolation and denoising are calculated according to the seismic data of interpolation and denoising after sparse transformation.The present invention can be directed to the random noise of different-energy, realize interpolation and the denoising simultaneously of full automatic seismic data, avoid manual adjustment parameter, save human cost, computing resource and calculating time.

Description

The interpolation and denoising method and system of a kind of seismic data
Technical field
The present invention relates to seismic data noise attenuation fields, interpolation and denoising method more particularly to a kind of seismic data and are System.
Background technique
In seismic data acquisition process, due to the limitation of surface obstructions object and cost, the presence in river, lake etc. makes The seismic data that must be acquired is not able to satisfy the design of observation system, the case where seismic data missing often occurs;In addition, based on pressure The earthquake data acquisition of contracting perception theory can be random arrangement wave detector, therefore obtain the data of undersampling.Above-mentioned acquisition Imperfect seismic data is required to recover the seismic data of complete sample using seismic data interpolation technique.Meanwhile it making an uproar at random Sound is generally existing in seismic exploration data, therefore often contains random noise in acquiring incomplete seismic data.Containing with Interpolation and denoising can realize two problems of interpolation and denoising in a process while seismic data of machine noise, for The reconstruction efficiency for improving seismic data is of great significance.It can be realized to contain currently based on the interpolation method of sparse transformation and make an uproar at random Interpolation and denoising while sound data.But currently based on interpolation while the seismic data of 1 norm regularization and denoising method The energy for noise is needed, the artificial suitable threshold parameter of selection could obtain reliable effect for different data. Since seismic exploration data amount is big, noise energy is different, the method for artificial selection threshold parameter can expend a large amount of computing resource, Human cost, time cost.
Summary of the invention
The object of the present invention is to provide a kind of interpolation of seismic data and denoising method and system, can be directed to different-energy Random noise, realize full automatic seismic data simultaneously interpolation and denoising, avoid manual adjustment parameter, save manpower at This, computing resource and calculating time.
To achieve the above object, the present invention provides following schemes:
A kind of interpolation and denoising method of seismic data, the interpolation and denoising method include:
Obtain the seismic data of Noise;
Regular parameter sequence is determined according to the seismic data of the Noise;Regular parameter in the regular parameter sequence It is sequentially reduced;
The seismic data of interpolation and denoising after calculating sparse transformation according to the regular parameter sequence;
By sparse transformation method, after calculating interpolation and denoising according to the seismic data of interpolation and denoising after sparse transformation Seismic data.
Optionally, the seismic data according to the Noise determines regular parameter sequence, specifically includes:
Sparse transformation is carried out to the seismic data of the Noise;
Maximum regular parameter in the regular sequence is determined according to the seismic data of the Noise after sparse transformation.
Optionally, the seismic data of the interpolation after the calculating sparse transformation according to the regular parameter sequence and denoising, It specifically includes:
Obtain Optimized model;The optimal solution of the Optimized model be sparse transformation after interpolation and denoising seismic data;
Multiple solutions of the Optimized model are calculated according to the regular parameter sequence;
Multiple control parameters are calculated according to the regular parameter sequence and the solution of the Optimized model;Maximum control parameter The solution of corresponding Optimized model is the optimal solution.
A kind of interpolation and denoising system of seismic data, the interpolation and denoising system include:
Module is obtained, for obtaining the seismic data of Noise;
Determining module, for determining regular parameter sequence according to the seismic data of the Noise;The regular parameter sequence Regular parameter in column is sequentially reduced;
First computing module, for the earthquake according to interpolation and denoising after regular parameter sequence calculating sparse transformation Data;
Second computing module, for passing through sparse transformation method, according to the earthquake number of interpolation and denoising after sparse transformation According to the seismic data after calculating interpolation and denoising.
Optionally, the determining module, specifically includes:
Sparse transformation unit carries out sparse transformation for the seismic data to the Noise;
Determination unit, it is maximum in the regular sequence for being determined according to the seismic data of the Noise after sparse transformation Regular parameter.
Optionally, first computing module, specifically includes:
Acquiring unit, for obtaining Optimized model;The optimal solution of the Optimized model is the interpolation after sparse transformation and goes The seismic data made an uproar;
Computing unit, for calculating multiple solutions of the Optimized model according to the regular parameter sequence;
Control parameter computing unit, it is multiple for being calculated according to the solution of the regular parameter sequence and the Optimized model Control parameter;The solution of the corresponding Optimized model of maximum control parameter is the optimal solution.
Compared with prior art, the present invention has following technical effect that the present invention according to the seismic data of the Noise Determine regular parameter sequence;The seismic data of interpolation and denoising after calculating sparse transformation according to the regular parameter sequence;It is logical Sparse transformation method is crossed, the earthquake number after interpolation and denoising are calculated according to the seismic data of interpolation and denoising after sparse transformation According to.This method automatically determines regular parameter using the thought of inversion theory, therefore can automatically realize the same of seismic data When interpolation and denoising, the present invention can save cost of labor, calculate cost and time cost, so that seismic data is slotting simultaneously Value and denoising become an automatic process.
Detailed description of the invention
It in order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, below will be to institute in embodiment Attached drawing to be used is needed to be briefly described, it should be apparent that, the accompanying drawings in the following description is only some implementations of the invention Example, for those of ordinary skill in the art, without any creative labor, can also be according to these attached drawings Obtain other attached drawings.
Fig. 1 is the flow chart of the interpolation and denoising method of seismic data of the embodiment of the present invention;
Fig. 2 is the original schematic diagram without noise simulation data;
Fig. 3 is one Noise seismic data of embodiment, signal-to-noise ratio=0.5333dB schematic diagram;
Fig. 4 is the seismic data that embodiment one samples, the schematic diagram that sample rate is 70%;
Fig. 5 is one adaptive-interpolation of embodiment and denoising result, signal-to-noise ratio=9.0824db schematic diagram;
Fig. 6 is the schematic diagram of two Noise Data of embodiment;
Fig. 7 is the schematic diagram of the sampled data of embodiment 2 60%;
Fig. 8 is the schematic diagram of embodiment two full-automatic interpolation and denoising simultaneously;
Fig. 9 is the schematic diagram for the noise that two automatic interpolation of embodiment and denoising method remove;
Figure 10 is the interpolation of seismic data and the structural block diagram of denoising system.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other Embodiment shall fall within the protection scope of the present invention.
The object of the present invention is to provide a kind of interpolation of seismic data and denoising method and system, can be directed to different-energy Random noise, realize full automatic seismic data simultaneously interpolation and denoising, avoid manual adjustment parameter, save manpower at This, computing resource and calculating time.
In order to make the foregoing objectives, features and advantages of the present invention clearer and more comprehensible, with reference to the accompanying drawing and specific real Applying mode, the present invention is described in further detail.
As shown in Figure 1, the interpolation and denoising method of a kind of seismic data include:
Step 101: obtaining the seismic data of Noise.
Step 102: regular parameter sequence is determined according to the seismic data of the Noise;In the regular parameter sequence Regular parameter is sequentially reduced.It specifically includes:
Sparse transformation is carried out to the seismic data of the Noise;
Maximum regular parameter in the regular sequence is determined according to the seismic data of the Noise after sparse transformation.
Step 103: the seismic data of interpolation and denoising after calculating sparse transformation according to the regular parameter sequence.Specifically Include:
Obtain Optimized model;The optimal solution of the Optimized model be sparse transformation after interpolation and denoising seismic data;
Multiple solutions of the Optimized model are calculated according to the regular parameter sequence;
Multiple control parameters are calculated according to the regular parameter sequence and the solution of the Optimized model;Maximum control parameter The solution of corresponding Optimized model is the optimal solution.
Step 104: by sparse transformation method, interpolation being calculated according to the seismic data of interpolation and denoising after sparse transformation With the seismic data after denoising.
The principle of this method is as follows:
Assuming that seismic data is d, noise ε, seismic data sample operator is Φ, the noise-containing seismic data of observation For dobs, relationship between three are as follows:
Φ d+ ε=dobs (1)
Interpolation and denoising are exactly by observation Noise Data d to seismic data simultaneouslyobsSeek the d of not Noise.According to sparse The denoising method of transformation, this patent is using warp wavelet as sparse transformation, then expression x=Ψ d of the d after sparse transformation is Sparse, therefore above formula becomes:
ΦΨ*X+ ε=dobs (2)
Wherein Ψ indicates warp wavelet (can be two dimension, be also possible to three-dimensional), Ψ*Indicate the inverse transformation of warp wavelet, Due to the sparsity of x, optimization problem below is established:
Wherein λiReferred to as regular parameter, A=Φ Ψ*, A is compound, the x of sample operator Φ and warp wavelet operator ΨiIt indicates By regular parameter λiDetermining seismic data warp wavelet domain coefficient, | | xi||1Indicate xi1 norm.
Give a regular parameter sequence λ being gradually reduced1> λ2> λ3> ... > λN, decline mode can be index side Formula is also possible to linear mode, λ1Take Ψ dobs(ΨdobsIt indicates to observation data dobsDo warp wavelet) maximum value 1/2, other ratios such as 1/3 or 2/3, λNFor one close to 0 positive real number, N indicate sequence length.
For problemDefinition initial solution is x0=dobs;Using iteration method to its into Primary solve of row obtains x1=soft (x0+AT(dobs-Ax0),λ1), soft (x, λ)=sign (x) .*max (abs (x)-λ, 0) table Show and soft-threshold operation is carried out to each component of x, it can also be using hard -threshold and other threshold operations;For each λi, i= 2 ..., N equally use iteration method pairIt carries out primary solve and obtains xi=soft (xi-1+ AT(dobs-Axi-1),λi), wherein xi-1For with iteration method pairAn iteration is carried out to obtain The approximate solution arrived.
DefinitionWork as Pi(PiIt is a customized parameter, is automatically stopped for control algolithm) it obtains When maximum, result d=Ψ is exported*xi, d is the result of interpolation and denoising simultaneously.
Embodiment one: Fig. 2 is the simulation seismic data of a not Noise, size of data 726*96, and time shaft has 726 A sampled point, time sampling interval 4ms, Fig. 3 are to increase noise to Fig. 2 as a result, signal-to-noise ratio is 0.5333dB, and Fig. 4 is pair The result of Fig. 3 stochastical sampling 70%;Fig. 5 is adaptive-interpolation and denoising result, signal-to-noise ratio 9.0824dB, proposed by the present invention Automatic interpolation and denoising method do not have to adjust any parameter, can obtain reliable interpolation and denoising result.
Embodiment two: for the true seismic data of a 276*400, Temporal sampling 4ms, time sampling number Be 276, spatially there is 400, original Noise Data as shown in fig. 6, its sampling 60% result such as Fig. 7, automatic interpolation and The denoising result of denoising is shown in that the difference of Fig. 8, Fig. 6 and Fig. 8 are shown in Fig. 9, it can be seen that this method can be realized automatic interpolation and denoising.
The specific embodiment provided according to the present invention, the invention discloses following technical effects: the present invention contains according to described The seismic data of noise determines regular parameter sequence;Interpolation and denoising after calculating sparse transformation according to the regular parameter sequence Seismic data;By sparse transformation method, interpolation is calculated according to the seismic data of interpolation and denoising after sparse transformation and is gone Seismic data after making an uproar.This method automatically determines regular parameter using the thought of inversion theory, therefore can automatically realize Interpolation and denoising while seismic data, the present invention can save cost of labor, cost and time cost be calculated, so that ground Interpolation and denoising become an automatic process to shake data simultaneously.
As shown in Figure 10, a kind of interpolation of seismic data and denoising system include:
Module 1001 is obtained, for obtaining the seismic data of Noise.
Determining module 1002, for determining regular parameter sequence according to the seismic data of the Noise;The canonical ginseng Regular parameter in Number Sequence is sequentially reduced.
The determining module 1002, specifically includes:
Sparse transformation unit carries out sparse transformation for the seismic data to the Noise;
Determination unit, it is maximum in the regular sequence for being determined according to the seismic data of the Noise after sparse transformation Regular parameter.
First computing module 1003, for according to the interpolation and denoising after regular parameter sequence calculating sparse transformation Seismic data.
First computing module 1003, specifically includes:
Acquiring unit, for obtaining Optimized model;The optimal solution of the Optimized model is the interpolation after sparse transformation and goes The seismic data made an uproar;
Computing unit, for calculating multiple solutions of the Optimized model according to the regular parameter sequence;
Control parameter computing unit, it is multiple for being calculated according to the solution of the regular parameter sequence and the Optimized model Control parameter;The solution of the corresponding Optimized model of maximum control parameter is the optimal solution.
Second computing module 1004, for passing through sparse transformation method, according to the ground of interpolation and denoising after sparse transformation It shakes data and calculates the seismic data after interpolation and denoising.
Each embodiment in this specification is described in a progressive manner, the highlights of each of the examples are with other The difference of embodiment, the same or similar parts in each embodiment may refer to each other.For system disclosed in embodiment For, since it is corresponded to the methods disclosed in the examples, so being described relatively simple, related place is said referring to method part It is bright.
Used herein a specific example illustrates the principle and implementation of the invention, and above embodiments are said It is bright to be merely used to help understand method and its core concept of the invention;At the same time, for those skilled in the art, foundation Thought of the invention, there will be changes in the specific implementation manner and application range.In conclusion the content of the present specification is not It is interpreted as limitation of the present invention.

Claims (6)

1. the interpolation and denoising method of a kind of seismic data, which is characterized in that the interpolation and denoising method include:
Obtain the seismic data of Noise;
Regular parameter sequence is determined according to the seismic data of the Noise;Regular parameter in the regular parameter sequence is successively Reduce;
The seismic data of interpolation and denoising after calculating sparse transformation according to the regular parameter sequence;
Ground by sparse transformation method, after interpolation and denoising are calculated according to the seismic data of interpolation and denoising after sparse transformation Shake data.
2. the interpolation and denoising method of seismic data according to claim 1, which is characterized in that described according to described noisy The seismic data of sound determines regular parameter sequence, specifically includes:
Sparse transformation is carried out to the seismic data of the Noise;
Maximum regular parameter in the regular sequence is determined according to the seismic data of the Noise after sparse transformation.
3. the interpolation and denoising method of seismic data according to claim 1, which is characterized in that described according to the canonical The seismic data of interpolation and denoising after argument sequence calculating sparse transformation, specifically includes:
Obtain Optimized model;The optimal solution of the Optimized model be sparse transformation after interpolation and denoising seismic data;
Multiple solutions of the Optimized model are calculated according to the regular parameter sequence;
Multiple control parameters are calculated according to the regular parameter sequence and the solution of the Optimized model;Maximum control parameter is corresponding Optimized model solution be the optimal solution.
4. the interpolation and denoising system of a kind of seismic data, which is characterized in that the interpolation and denoising system include:
Module is obtained, for obtaining the seismic data of Noise;
Determining module, for determining regular parameter sequence according to the seismic data of the Noise;In the regular parameter sequence Regular parameter be sequentially reduced;
First computing module, for the earthquake number according to interpolation and denoising after regular parameter sequence calculating sparse transformation According to;
Second computing module, for passing through sparse transformation method, according to the seismic data meter of interpolation and denoising after sparse transformation Seismic data after calculating interpolation and denoising.
5. the interpolation and denoising system of seismic data according to claim 1, which is characterized in that the determining module, tool Body includes:
Sparse transformation unit carries out sparse transformation for the seismic data to the Noise;
Determination unit, for determining maximum canonical in the regular sequence according to the seismic data of the Noise after sparse transformation Parameter.
6. the interpolation and denoising system of seismic data according to claim 4, which is characterized in that described first calculates mould Block specifically includes:
Acquiring unit, for obtaining Optimized model;The optimal solution of the Optimized model is the interpolation and denoising after sparse transformation Seismic data;
Computing unit, for calculating multiple solutions of the Optimized model according to the regular parameter sequence;
Control parameter computing unit, for calculating multiple controls according to the solution of the regular parameter sequence and the Optimized model Parameter;The solution of the corresponding Optimized model of maximum control parameter is the optimal solution.
CN201910664806.3A 2019-07-23 2019-07-23 Interpolation and denoising method and system for seismic data Active CN110261912B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910664806.3A CN110261912B (en) 2019-07-23 2019-07-23 Interpolation and denoising method and system for seismic data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910664806.3A CN110261912B (en) 2019-07-23 2019-07-23 Interpolation and denoising method and system for seismic data

Publications (2)

Publication Number Publication Date
CN110261912A true CN110261912A (en) 2019-09-20
CN110261912B CN110261912B (en) 2020-06-16

Family

ID=67927751

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910664806.3A Active CN110261912B (en) 2019-07-23 2019-07-23 Interpolation and denoising method and system for seismic data

Country Status (1)

Country Link
CN (1) CN110261912B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111007566A (en) * 2019-12-27 2020-04-14 西南石油大学 Curvature-driven diffusion full-convolution network seismic data bad channel reconstruction and denoising method
CN117055107A (en) * 2023-06-29 2023-11-14 成都理工大学 Seismic interpolation method based on interaction of Framelet transformation and Lp pseudo-norms

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111007566A (en) * 2019-12-27 2020-04-14 西南石油大学 Curvature-driven diffusion full-convolution network seismic data bad channel reconstruction and denoising method
CN117055107A (en) * 2023-06-29 2023-11-14 成都理工大学 Seismic interpolation method based on interaction of Framelet transformation and Lp pseudo-norms
CN117055107B (en) * 2023-06-29 2024-04-09 成都理工大学 Seismic interpolation method based on interaction of Framelet transformation and Lp pseudo-norms

Also Published As

Publication number Publication date
CN110261912B (en) 2020-06-16

Similar Documents

Publication Publication Date Title
CN102998706B (en) Method and system for attenuating seismic data random noise
CN104536044B (en) The interpolation denoising method and system of a kind of geological data
CN105974468B (en) A kind of method that can be carried out at the same time five dimension Reconstruction of seismic data and noise compacting
SE515924C2 (en) Methods for time-lapse examination of an oil / gas source
CN106842321B (en) seismic data reconstruction method and device
CN104914466B (en) A kind of method for improving seismic data resolution
CN109633752B (en) Offshore towing cable data self-adaptive ghost wave compression method based on three-dimensional fast Radon transformation
CN107179550B (en) A kind of seismic signal zero phase deconvolution method of data-driven
CN110261912A (en) The interpolation and denoising method and system of a kind of seismic data
CN107144879A (en) A kind of seismic wave noise-reduction method combined based on adaptive-filtering with wavelet transformation
CN110174702A (en) A kind of method and system that marine seismic data low frequency weak signal is restored
CN110687597A (en) Wave impedance inversion method based on joint dictionary
CN107367760B (en) Based on the surface-related multiple and higher-order spectra method and system for accelerating linear Bregman algorithm
Sun et al. Multiple attenuation using λ-f domain high-order and high-resolution Radon transform based on SL0 norm
CN105676292A (en) 3D earthquake data de-noising method based on 2D curvelet transform
CN106950597B (en) Mixing source data separation method based on the filtering of three sides
CN110749925A (en) Broadband reverse time migration imaging processing method
CN116203634A (en) Ghost wave removing method based on low-rank constraint
CN116719088A (en) Aviation transient electromagnetic data noise suppression method
CN116520418A (en) Efficient extraction method for elastic wave angle domain common imaging point gather
CN103901470B (en) Seabed settled layer reflection strength and time delay estimation method
CN109212590B (en) Amplitude-preserving reverse time migration low-frequency noise suppression method and system
Duan et al. Band-limited impedance perturbation inversion using cross-correlative least-squares RTM
Lopez et al. SRME and estimation of primaries by sparse inversion: a hybrid approach
CN113281808B (en) Anti-dispersion seismic wave forward modeling method, system, device and medium

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant