CN106908840A - Seismic data Hz noise automatic identification and drawing method based on principal component analysis - Google Patents
Seismic data Hz noise automatic identification and drawing method based on principal component analysis Download PDFInfo
- Publication number
- CN106908840A CN106908840A CN201710320396.1A CN201710320396A CN106908840A CN 106908840 A CN106908840 A CN 106908840A CN 201710320396 A CN201710320396 A CN 201710320396A CN 106908840 A CN106908840 A CN 106908840A
- Authority
- CN
- China
- Prior art keywords
- noise
- matrix
- value
- seismic data
- average
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000000513 principal component analysis Methods 0.000 title claims abstract description 13
- 238000012545 processing Methods 0.000 claims abstract description 7
- 239000013598 vector Substances 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 36
- 238000000605 extraction Methods 0.000 claims description 6
- 238000012512 characterization method Methods 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 239000000284 extract Substances 0.000 abstract description 3
- 238000001914 filtration Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012880 independent component analysis Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
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 present invention relates to a kind of seismic data Hz noise automatic identification based on principal component analysis and drawing method, by calculating each road the average energy value of earthquake record, according to the average energy value otherness automatic identification Hz noise road, extract some neighboring track composition seismic interference road collection comprising interference way, according to the strong correlation of Hz noise noise, characteristic vector is reconstructed using principal component analysis technology reach the purpose of compacting Hz noise.The energy of useful signal is weakened while denoising due to introducing principal component analysis, is that this has carried out relative amplitude preserved processing again to reconstruct data, it is ensured that the trace consistency of earthquake record.Empirical tests, a kind of seismic data Hz noise automatic identification based on principal component analysis disclosed by the invention can realize the automatic identification to Hz noise in seismic data and compacting with drawing method, desert area seismic data quality is effectively increased, it is significant especially for the bad track for repairing the original seismic data caused by Hz noise.
Description
Technical field:
The present invention relates to a kind of method of seismic data process denoising, the especially serious earthquake record denoising of Hz noise
Method, using the average energy value identification interference way in each road of earthquake record and using principal component analysis technology compacting Hz noise
Method.
Background technology:
The influence of Hz noise is suffered from during acquiring seismic exploration data, is often mixed with the earthquake record for collecting
Industrial frequency noise, when Hz noise is very strong, and through one geological data all the time when, the seismic channel often shows as bad track, draws
Seismic data quality is played to decline.By the analysis to a large amount of bad track data, many seismic channels do not stem from the strong of Near Ground
Interference source, often from geophone problem in itself, the time depth with earthquake record is unrelated, currently to the treatment of bad track
It is to reject bad track, this treatment reduces earthquake data quality, and the technology that this patent is provided, can effectively recognize due to power frequency
The bad track for producing is disturbed, this provides possibility further to suppress Hz noise.Do not have very much also currently for Hz noise
The solution of effect.Drawing method main at present has wave trap method, power frequency recurrence to subtract each other suppression method, the class of self-adaptive routing three.
Article " automatic identification of Hz noise and compacting in seismic data " proposes that wave trap method filters industry disturbance, but the method is simultaneously
The seismic signal of 50Hz frequencies can be filtered;Article " extraction of disturbance of industry frequency ripple and removing method " is proposed in time-domain with altogether
Yoke gradient algorithm extracts mono-tone interference, but the party have ignored the changeability of work frequency in actual seismic data, and research shows
Different zones, different time, Hz noise frequency is change;Article " force in compacting seismic data by the cosine of industrial noise
The improvement and application of nearly method " proposes to filter industry disturbance based on cosine approach, but the method cannot realize that automatic identification is disturbed
Road;Article " separating 50Hz industry disturbances in industry disturbance in geophone domain " proposes that the conversion of big gun collection filters industry disturbance, but the method efficiency is low simultaneously
Cannot automatic identification interference way;Article " 50Hz industry disturbances denoising method and application based on Wiener filtering " proposes Wiener filtering
Single-frequency noise is predicted in matching of the device to adjust the amplitude and phase of library track, but the method is not suitable for the power frequency of stationarity difference
Interference;Article " the Hz noise technology for eliminating based on independent component analysis " proposition is isolated separate from mixed signal
Each component of signal reaches the purpose for eliminating man-made noise, but the method lacks strict theoretical foundation by empirical;
" a kind of method of automatic identification based on dual factors with industry disturbance is removed " is proposed disclosed in CN104570118A
Sine and cosine Weighted approximation method processes industry disturbance, but the method have ignored the frequency of power frequency, the not stationarity of phase and amplitude;
" method of alternating current disturbance signal in removal geological data " proposes that frequency domain is extracted disclosed in CN103630935A
The method of work frequency, but the method does not account for the presence of industrial frequency harmonic;
" a kind of automatic identification and the method for eliminating seismic exploration industry electrical interference " proposes base disclosed in CN101907726A
Determine Hz noise in self-correlation theory, but seismic data before first arrival time is depended on during the method automatic identification industrial frequency noise
Quality.
The content of the invention:
The purpose of the present invention is that for above-mentioned the deficiencies in the prior art, there is provided a kind of earthquake based on principal component analysis
Data Hz noise automatic identification and drawing method.
Main idea is that:Seismic data is often influenceed by Hz noise during earthquake-capturing, from
And a large amount of bad tracks or strong jamming road in earthquake record are caused, and the quality of seismic exploration data is had a strong impact on, the present invention is to pass through
Then automatic identification is isolated interference by the seismic channel of Hz noise by principal component analysis technology, then to the data after reconstruct
Guarantor's width is carried out, automatic identification and the compacting of seismic data Hz noise is realized.
The present invention is achieved by the following technical solutions:
Seismic data Hz noise automatic identification and drawing method based on principal component analysis, comprise the following steps:
A, for pending original seismic data, choose the when window T that the record time is located at middle hypomereW, calculate the earthquake
Record is in TWNei Ge roads the average energy value, such as formula
Wherein EiFor when window in the i-th road of earthquake record the average energy value, xi(tj) for when window in the road of earthquake record i-th jth
Individual sampling point value, i=1,2 ... n, n be the total road number of earthquake record, j=1,2 ... k, when k is in window earthquake record sampling number,
When window TWLength suggestion for summary journal lengthArriveIn time;
B, definition l roads characterization factorAccording to λlValue identification interference way, ElFor earthquake record l roads exist
TWInterior the average energy value, q is TWThe interior earthquake Taoist monastic name with minimum average B configuration energy value, wherein l=1,2 ... n recognize interference way
Comprise the following steps that:
If b1, there is λ to all roadsl< 5, shows the earthquake record without notable Hz noise, without follow-up denoising
Journey, otherwise performs step b2;
If b2, λl>=5, show that earthquake record has notable Hz noise, note l roads are interference way, are designated as Ul;
C, extraction and UlIt is adjacent and including UlContinuous 5 trace record constitute sub- seismic interference road collection, X is designated as, per together in X
It is a row, obtains the average of X each columns, and the average of the row is subtracted with each element in X, obtains the sub- seismic channel after X treatment
Collection, is designated as X*;
D, calculating X*Covariance matrix Γ, such as formula
Wherein b is X*Columns, X*TIt is X*Transposed matrix, " " representing matrix multiplication;
E, eigenvalue matrix Λ and eigenvectors matrix R that covariance matrix Γ is calculated using singular value decomposition method, then
There is formula
Γ=R Λ RT (3)
Wherein Λ is characterized the diagonal matrix that value is arranged from big to small, each characteristic vectors for being classified as character pair value of R, RT
It is transposed matrix, meets RTR=RRT=E, E are unit matrix;
F, X are through RTLinear Mapping, obtains principal component matrix, such as formula
Φ=RT·X (4)
Wherein Φ is main component matrix;
G, Φ the first rows are set to 0, obtain Φ ', made
X1 *=R Φ ' (5)
Wherein X1 *It is restructuring matrix;
H, extraction X1 *The data of middle correspondence interference way, are designated as Ul', to Ul' carry out relative amplitude preserved processing, such as formula
Wherein Ul *It is interference way signal, A after relative amplitude preserved processing1It is the average value of interference way signal amplitude absolute value after denoising,
A2It is and UlThe absolute value of amplitude average value sum of adjacent two track data it is average.
Beneficial effect:Through experiment, a kind of seismic data Hz noise based on principal component analysis disclosed by the invention from
Dynamic identification can be realized suppressing the noise that industrial frequency interference source brings in seismic prospecting with drawing method, and being capable of automatic identification
Interference way, the algorithm can repair the bad track that Hz noise causes, and effectively increase desert area seismic data quality, reduce and gather into
This.
Brief description of the drawings:
Fig. 1 real seismic records
Earthquake record after Fig. 2 compacting Hz noises
Specific embodiment:
The present invention is described in further detail with reference to the accompanying drawings and examples:
Seismic data Hz noise automatic identification and drawing method based on principal component analysis, comprise the following steps:
A, using actual single shot record, this example TW=700ms~800ms, calculates TWInterior Ei, such as formula
Wherein EiIt is the i-th road of earthquake record the average energy value, xi(tj) it is j-th sampled point unit in the road of earthquake record i-th
Element, i=1 in this example, 2 ... n, n=48, j=1,2 ... k, k=100;
B, definition l roads characterization factorAccording to λlValue identification interference way, ElFor earthquake record l roads exist
TWInterior the average energy value, q is T in this exampleWInterior first, E1=0.0042, wherein l=1,2 ... n, identification interference stage property
Body step is as follows:
If b1, there is λ to all roadsl< 5, shows the earthquake record without notable Hz noise, without follow-up denoising
Journey, meets λ in this examplelThe road of < 5 is the 1st~7 road, the 9th~42 road, otherwise the 44th~48 road, execution step b2;
If b2, λl>=5, show that earthquake record has notable Hz noise, note l roads are interference way, are designated as Ul, it is full in this example
Sufficient λl>=5 road is the 8th and the 43rd road, is designated as U8And U43, E8=0.078, E43=0.23, λ8=18.57, λ43=54.76;
C, with U8As a example by, extract the 6th, 7,8,9,10 road and constitute sub- seismic interference road collection, X is designated as, every in X is together one
Row, obtain the average of X each columns, and the average of the row is subtracted with each element in X, obtain the sub- seismic channel set after X treatment, note
It is X*;
D, calculating X*Covariance matrix Γ, such as formula
Wherein b is X*Columns, X*TIt is X*Transposed matrix, " " representing matrix multiplication, b=5 in formula;
E, eigenvalue matrix Λ and eigenvectors matrix R that covariance matrix Γ is calculated using singular value decomposition method, then
There is formula
Γ=R Λ RT (3)
Wherein Λ is characterized the diagonal matrix that value is arranged from big to small, each characteristic vectors for being classified as character pair value of R, RT
It is transposed matrix, meets RTR=RRT=E, E are unit matrix;
F, X are through RTLinear Mapping, obtains principal component matrix, such as formula
Φ=RT·X (4)
Wherein Φ is main component matrix;
G, Φ the first rows are set to 0, obtain Φ ', made
X1 *=R Φ ' (5)
Wherein X1 *It is restructuring matrix;
H, extraction X1 *The data of middle correspondence interference way, are designated as Ul', to Ul' carry out relative amplitude preserved processing, such as formula
Wherein Ul *It is interference way signal, A after relative amplitude preserved processing1It is the average value of interference way signal amplitude absolute value after denoising,
A2It is and UlAverage, the A in this example of the absolute value of amplitude average value sum of adjacent two track data2=0.0689, A1=0.0057.
Claims (1)
1. a kind of seismic data Hz noise automatic identification and drawing method based on principal component analysis, comprise the following steps:
A, for pending original seismic data, choose the when window T that the record time is located at middle hypomereW, calculate the earthquake record and exist
TWNei Ge roads the average energy value, such as formula
Wherein EiFor when window in the i-th road of earthquake record the average energy value, xi(tj) for when window in j-th of the road of earthquake record i-th adopt
Sample value, i=1,2 ... n, n be the total road number of earthquake record, j=1,2 ... k, when k is in window earthquake record sampling number, when window
TWLength suggestion for summary journal lengthArriveIn time;
B, definition l roads characterization factorAccording to λlValue identification interference way, ElIt is earthquake record l roads in TWIt is interior
The average energy value, q is TWThe interior earthquake Taoist monastic name with minimum average B configuration energy value, wherein l=1,2 ... n, identification interference stage property
Body step is as follows:
If b1, there is λ to all roadsl< 5, shows the earthquake record without notable Hz noise, no without follow-up denoising process
Then perform step b2;
If b2, λl>=5, show that earthquake record has notable Hz noise, note l roads are interference way, are designated as Ul;
C, extraction and UlIt is adjacent and including UlContinuous 5 trace record constitute sub- seismic interference road collection, X is designated as, per being together one in X
Row, obtain the average of X each columns, and the average of the row is subtracted with each element in X, obtain the sub- seismic channel set after X treatment, note
It is X*;
D, calculating X*Covariance matrix Γ, such as formula
Wherein b is X*Columns, X*TIt is X*Transposed matrix, " " representing matrix multiplication;
E, eigenvalue matrix Λ and eigenvectors matrix R that covariance matrix Γ is calculated using singular value decomposition method, then be present
Formula
Γ=R Λ RT (3)
Wherein Λ is characterized the diagonal matrix that value is arranged from big to small, each characteristic vectors for being classified as character pair value of R, RTTo turn
Matrix is put, R is metTR=RRT=E, E are unit matrix;
F, X are through RTLinear Mapping, obtains principal component matrix, such as formula
Φ=RT·X (4)
Wherein Φ is main component matrix;
G, Φ the first rows are set to 0, obtain Φ ', made
X1 *=R Φ ' (5)
Wherein X1 *It is restructuring matrix;
H, extraction X1 *The data of middle correspondence interference way, are designated as Ul', to Ul' carry out relative amplitude preserved processing, such as formula
Wherein Ul *It is interference way signal, A after relative amplitude preserved processing1It is the average value of interference way signal amplitude absolute value after denoising, A2Be with
UlThe absolute value of amplitude average value sum of adjacent two track data it is average.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710320396.1A CN106908840A (en) | 2017-05-09 | 2017-05-09 | Seismic data Hz noise automatic identification and drawing method based on principal component analysis |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710320396.1A CN106908840A (en) | 2017-05-09 | 2017-05-09 | Seismic data Hz noise automatic identification and drawing method based on principal component analysis |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106908840A true CN106908840A (en) | 2017-06-30 |
Family
ID=59210031
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710320396.1A Pending CN106908840A (en) | 2017-05-09 | 2017-05-09 | Seismic data Hz noise automatic identification and drawing method based on principal component analysis |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106908840A (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108732624A (en) * | 2018-05-29 | 2018-11-02 | 吉林大学 | A kind of parallel focus seismic data stochastic noise suppression method based on PCA-EMD |
CN108957552A (en) * | 2018-07-17 | 2018-12-07 | 吉林大学 | Seismic data wave noise drawing method based on SS-PCA |
CN108957550A (en) * | 2018-06-28 | 2018-12-07 | 吉林大学 | The strong industrial noise drawing method of TSP based on SVD-ICA |
CN112036234A (en) * | 2020-07-16 | 2020-12-04 | 成都飞机工业(集团)有限责任公司 | PCA-based aircraft conduit vibration signal power frequency noise suppression method |
CN112180447A (en) * | 2019-07-04 | 2021-01-05 | 中国石油天然气集团有限公司 | Method and system for eliminating strong reflection shielding of reservoir |
CN113257268A (en) * | 2021-07-02 | 2021-08-13 | 成都启英泰伦科技有限公司 | Noise reduction and single-frequency interference suppression method combining frequency tracking and frequency spectrum correction |
-
2017
- 2017-05-09 CN CN201710320396.1A patent/CN106908840A/en active Pending
Non-Patent Citations (2)
Title |
---|
杨重洋: "地震资料中工频干扰的自动识别与压制", 《石油工业计算机应用》 * |
王权海: "子空间分析方法在地震勘探等信号处理中的初步应用研究", 《中国博士学位论文全文数据库 基础科学辑》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108732624A (en) * | 2018-05-29 | 2018-11-02 | 吉林大学 | A kind of parallel focus seismic data stochastic noise suppression method based on PCA-EMD |
CN108957550A (en) * | 2018-06-28 | 2018-12-07 | 吉林大学 | The strong industrial noise drawing method of TSP based on SVD-ICA |
CN108957552A (en) * | 2018-07-17 | 2018-12-07 | 吉林大学 | Seismic data wave noise drawing method based on SS-PCA |
CN112180447A (en) * | 2019-07-04 | 2021-01-05 | 中国石油天然气集团有限公司 | Method and system for eliminating strong reflection shielding of reservoir |
CN112036234A (en) * | 2020-07-16 | 2020-12-04 | 成都飞机工业(集团)有限责任公司 | PCA-based aircraft conduit vibration signal power frequency noise suppression method |
CN113257268A (en) * | 2021-07-02 | 2021-08-13 | 成都启英泰伦科技有限公司 | Noise reduction and single-frequency interference suppression method combining frequency tracking and frequency spectrum correction |
CN113257268B (en) * | 2021-07-02 | 2021-09-17 | 成都启英泰伦科技有限公司 | Noise reduction and single-frequency interference suppression method combining frequency tracking and frequency spectrum correction |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106908840A (en) | Seismic data Hz noise automatic identification and drawing method based on principal component analysis | |
Phillips | Designing matched bandpass and azimuthal filters for the separation of potential-field anomalies by source region and source type | |
CN105510976B (en) | A kind of many subwaves combination adaptive attenuation method | |
CN107144879B (en) | A kind of seismic wave noise-reduction method based on adaptive-filtering in conjunction with wavelet transformation | |
CN107589454A (en) | One kind is based on VMD TFPF compacting seismic prospecting random noise methods | |
CN104345341A (en) | Region constraint-based frequency band division energy seismic surface wave processing method | |
CN107219555B (en) | The strong industrial frequency noise drawing method of parallel focus seismic prospecting data based on principal component analysis | |
CN105445801B (en) | A kind of processing method for eliminating 2-d seismic data random noise | |
CN103399348A (en) | Denoising method for seismic signal based on Shearlet transform | |
CN109633761B (en) | Magnetic resonance signal power frequency noise reduction method based on wavelet transformation modulus maximum value method | |
CN110031899A (en) | Compressed sensing based weak signal extraction algorithm | |
Alsdorf | Noise reduction in seismic data using Fourier correction coefficient filtering | |
Lu et al. | Instantaneous polarization filtering focused on suppression of surface waves | |
CN103792574A (en) | Method for detecting frequency-variable gas in storage layer. | |
Zhou et al. | Unsupervised machine learning for waveform extraction in microseismic denoising | |
Li et al. | Distributed acoustic sensing vertical seismic profile data denoising based on multistage denoising network | |
DE69632892T2 (en) | Method for filtering elliptical waves propagating in a medium | |
CN103076626A (en) | Wave field purification treatment method | |
CN106019377A (en) | Two-dimensional seismic exploration noise removing method based on time-space-domain frequency reduction model | |
Bing et al. | A robust random noise suppression method for seismic data using sparse low-rank estimation in the time-frequency domain | |
Tian et al. | Noise suppression method for magnetic resonance sounding signals based on double singular value decomposition | |
Yu et al. | Research on the seismic signal denoising with the LMD and EMD method | |
Wang et al. | Multicomponent seismic noise attenuation with multivariate order statistic filters | |
Yuan et al. | Ground roll attenuation based on an empirical curvelet transform | |
CN101907726B (en) | Method for automatically identifying and eliminating industrial electrical interference in earthquake exploration |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170630 |
|
RJ01 | Rejection of invention patent application after publication |