CN107219555A - 基于主成分分析的并行震源地震勘探资料强工频噪声压制方法 - Google Patents
基于主成分分析的并行震源地震勘探资料强工频噪声压制方法 Download PDFInfo
- Publication number
- CN107219555A CN107219555A CN201710395695.1A CN201710395695A CN107219555A CN 107219555 A CN107219555 A CN 107219555A CN 201710395695 A CN201710395695 A CN 201710395695A CN 107219555 A CN107219555 A CN 107219555A
- Authority
- CN
- China
- Prior art keywords
- mrow
- mtd
- msub
- mtr
- frequency 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000000513 principal component analysis Methods 0.000 title claims abstract description 11
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 5
- 239000011159 matrix material Substances 0.000 claims description 39
- 238000013507 mapping Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 238000010183 spectrum analysis Methods 0.000 claims description 3
- 230000005284 excitation Effects 0.000 abstract description 4
- 238000012545 processing Methods 0.000 abstract description 4
- 230000000694 effects Effects 0.000 abstract description 3
- 238000004519 manufacturing process Methods 0.000 abstract description 2
- 238000012360 testing method Methods 0.000 abstract description 2
- BIIBYWQGRFWQKM-JVVROLKMSA-N (2S)-N-[4-(cyclopropylamino)-3,4-dioxo-1-[(3S)-2-oxopyrrolidin-3-yl]butan-2-yl]-2-[[(E)-3-(2,4-dichlorophenyl)prop-2-enoyl]amino]-4,4-dimethylpentanamide Chemical class CC(C)(C)C[C@@H](C(NC(C[C@H](CCN1)C1=O)C(C(NC1CC1)=O)=O)=O)NC(/C=C/C(C=CC(Cl)=C1)=C1Cl)=O BIIBYWQGRFWQKM-JVVROLKMSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000012847 principal component analysis method Methods 0.000 description 1
- 238000011084 recovery Methods 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
-
- 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
本发明公开一种基于主成分分析的并行震源地震勘探资料强工频噪声压制方法,并行震源地震勘探是在勘探区域内不同炮点处布置多个震源,采用近同时激发技术采集资料,可提高地震勘探生产效率和地震资料质量。在人文活动区并行震源激发采集得到的原始地震记录中经常混合了强工频噪声,对于干扰道根据相空间理论构造Hankel矩阵,并对其进行主成分分解与恢复有用信号。经验证,本发明能够实现对并行震源采集得到的地震勘探资料进行有效的强工频噪声压制,与传统的强工频噪声压制方法相比,处理数据快,不需要精确估计强工频噪声的频率,无相位偏移,信号保护性好,能够有效改善并行震源勘探的地震资料质量,降低了并行震源地震勘探资料处理成本。
Description
技术领域:
本发明涉及一种基于主成分分析的并行震源地震勘探资料强工频噪声压制方法,该方法用于并行震源地震勘探时强工频噪声的压制。
背景技术:
并行震源地震勘探在采集过程中往往受到强工频噪声的影响,采集得到并行地震资料中常混有强工频噪声,因此得到的地震资料品质不高,进而影响了后期地震资料解释和偏移成像质量。而关于并行震源地震勘探资料强工频噪声压制方法的文献尚少见,目前其他领域对强工频噪声抑制的方法主要有频域法和时域法两大类,频域法如陷波法,小波滤波法等都在频率域进行工频压制,虽然该类方法简单,方便,但当有用信号频率和工频的频率互相混叠时,这类方法会对有用信号造成一定程度的损害,同时也会产生一定的相移。时域法如正余弦逼近法,工频回归相减法,自适应滤波法等是将工频噪声表示为以振幅,频率,相位为变量的函数,再通过正余弦函数逼近法来估计强工频噪声,但这类方法需要对强工频噪声频率进行精确估计,同时这类方法对多个工频谐波存在时的情况需要多次处理,从而算法效率较低。可见,上述两类方法都不适应于并行震源地震勘探资料强工频噪声的压制。
发明内容:
本发明的目的就在于针对上述现有技术的不足,提供一种基于主成分分析的并行震源地震勘探资料强工频噪声压制方法。
本发明的主要思想是:
并行震源地震勘探技术提高了地震资料的质量和工作效率,但采集得到的地震资料往往受到强工频噪声的影响,这样获得的地震资料经常影响了后期地震资料解释和偏移成像质量,本发明是通过将干扰道信号根据相空间理论构造Hankel矩阵,并对其进行主成分分解与有效信号的恢复,从而实现了并行震源地震勘探资料强工频噪声的压制。
本发明是通过以下技术方案实现的:
基于主成分分析的并行震源地震勘探资料强工频噪声压制方法,包括以下步骤:
a、对并行地震勘探资料进行频谱分析,根据强工频噪声在频谱中幅值异常突出和频带窄,且频带在50Hz附近的特点,识别含有强工频噪声的干扰道。若干扰道只有一道,记干扰道信号为xi(l),其中i为并行地震勘探资料单炮地震记录的道号,l为采样点序号,l=1,2,···,N,N为干扰道信号的采样点数,执行一次步骤b~g;若干绕道为多个,对每一干扰道信号重复执行步骤b~g;
b、根据相空间重构理论,对xi(l)构造Hankel矩阵
该矩阵的行数记为m,列数记为n,若N为偶数,则令m=N/2,n=N/2,若N为奇数,则令m=(N/2)+1,n=(N/2)+1;
c、计算Hi的协方差矩阵Γi,如公式
Hi T为Hi的转置矩阵,“·”表示矩阵乘法;
d、利用奇异值分解法,计算协方差矩阵Γi的特征值矩阵Λi和特征向量矩阵Ri,则存在公式
Γi=Ri·Λi·Ri T (3)
其中Λi为由大到小排列的特征值矩阵,Ri为各个特征值对应的特征向量矩阵,Ri T为Ri的转置矩阵,且满足Ri T·Ri=Ri·Ri T=Ei,Ei为单位矩阵;
e、Hi经线性映射,得到主成分矩阵
Φi=Ri T·Hi (4)
f、设强工频基波和其主要谐波的数目为r,将Φi第一至第2r行所有元素置0,则得到重构的主成分矩阵Φ'i,令
Yi=Ri·Φi' (5)
则Yi为压制工频干扰后的Hankel矩阵,其具体形式可记为
g、定义xi'(l)=[yi(1),yi(2),···,yi(N)],则xi'(l)即为对应xi(l)的压制强工频噪声后的信号。
有益效果:
本发明涉及基于主成分分析的并行震源地震勘探资料强工频噪声压制方法,并行震源地震勘探是在勘探区域内不同炮点处布置多个震源,采用近同时激发技术采集资料,可提高地震勘探生产效率和地震资料质量。在人文活动区并行震源激发采集得到的原始地震记录中经常混合了强工频噪声。本发明所述的强工频噪声表现为噪声幅值较大,在地震资料上可见强烈的周期性干扰,有用信号被部分或全部淹没,造成地震资料特别是深部地震资料质量的大幅度下降,严重影响了后期地震资料解释和成像质量。针对该问题,本发明提出了基于主成分分析的并行震源地震勘探资料强工频噪声压制方法,对于干扰道根据相空间理论构造Hankel矩阵,并对其进行主成分分解与恢复有用信号。经验证,本发明公开的基于主成分分析的并行震源地震勘探资料强工频噪声压制方法能够实现对并行震源采集得到的地震勘探资料进行有效的强工频噪声压制,与传统的强工频噪声压制方法相比,该算法处理数据快,不需要精确估计强工频噪声的频率,无相位偏移,信号保护性好,能够有效改善并行震源勘探的地震资料质量,降低了并行震源地震勘探资料处理成本,
附图说明:
图1原始单道记录
图2压制强工频噪声后单道记录
具体实施方式:
下面结合附图和实施例做进一步的详细说明:
在本例中使用两个震源为一组的方法进行激发,雷克子波主频为50Hz,记录时间为1.2S,采样率为1000,在单炮记录的任意单道含频率为50Hz,振幅为1.2,相位为1.23的余弦信号,构成干扰道信号。
基于主成分分析的并行震源地震勘探资料强工频噪声压制方法,包括以下步骤:
a、对并行地震勘探资料进行频谱分析,根据强工频噪声在频谱中幅值异常突出和频带窄,且频带在50Hz附近的特点,识别含有强工频噪声的干扰道。若干扰道只有一道,记干扰道信号为xi(l),其中i为并行地震勘探资料单炮地震记录的道号,l为采样点序号,l=1,2,···,N,N为干扰道信号的采样点数,执行一次步骤b~g;若干绕道为多个,对每一干扰道信号重复执行步骤b~g,本例中干扰道为x10(l),l=1,2,···,1201,N=1201;
b、根据相空间重构理论,对x10(l)构造Hankel矩阵
该矩阵的行数记为m,列数记为n,若N为偶数,则令m=N/2,n=N/2,若N为奇数,则令m=(N/2)+1,n=(N/2)+1,本例中m=601,n=601;
c、计算Hi的协方差矩阵Γi,如公式
Hi T为Hi的转置矩阵,“·”表示矩阵乘法,本例中n=601;
d、利用奇异值分解法,计算协方差矩阵Γ10的特征值矩阵Λ10和特征向量矩阵R10,则存在公式
Γ10=R10·Λ10·R10 T (3)
其中Λ10为由大到小排列的特征值矩阵,R10为各个特征值对应的特征向量矩阵,R10 T为R10的转置矩阵,且满足R10 T·R10=R10·R10 T=E10,E10为单位矩阵;
e、H10经线性映射,得到主成分矩阵
Φ10=R10 T·H10 (4)
f、设强工频基波和其主要谐波的数目为r,将Φi第一至第2r行所有元素置0,则得到重构的主成分矩阵Φ'i,本例中r=1,将Φi第一至第2行所有元素置0,令
Y10=R10·Φ10' (5)
则Y10为压制工频干扰后的Hankel矩阵,其具体形式可记为
g、定义x10'(l)=[y10(1),y10(2),···,y10(N)],则x10'(l)即为对应x10(l)的压制强工频噪声后的信号。
Claims (2)
1.基于主成分分析的并行震源地震勘探资料强工频噪声压制方法,其特征在于,包括以下步骤:
a、对并行地震勘探资料进行频谱分析,根据强工频噪声在频谱中幅值异常突出和频带窄,且频带在50Hz附近的特点,识别含有强工频噪声的干扰道;
b、根据相空间重构理论,对xi(l)构造Hankel矩阵
<mrow>
<msub>
<mi>H</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mi>L</mi>
</mtd>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mi>L</mi>
</mtd>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>M</mi>
</mtd>
<mtd>
<mi>M</mi>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mi>M</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mi>n</mi>
<mo>+</mo>
<mn>2</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mi>L</mi>
</mtd>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
该矩阵的行数记为m,列数记为n,若N为偶数,则令m=N/2,n=N/2,若N为奇数,则令m=(N/2)+1,n=(N/2)+1;
c、计算Hi的协方差矩阵Γi,如公式
<mrow>
<msub>
<mi>&Gamma;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>H</mi>
<mi>i</mi>
</msub>
<mo>&CenterDot;</mo>
<msup>
<msub>
<mi>H</mi>
<mi>i</mi>
</msub>
<mi>T</mi>
</msup>
</mrow>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
Hi T为Hi的转置矩阵,“·”表示矩阵乘法;
d、利用奇异值分解法,计算协方差矩阵Γi的特征值矩阵Λi和特征向量矩阵Ri,则存在公式
Γi=Ri·Λi·Ri T (3)
其中Λi为由大到小排列的特征值矩阵,Ri为各个特征值对应的特征向量矩阵,Ri T为Ri的转置矩阵,且满足Ri T·Ri=Ri·Ri T=Ei,Ei为单位矩阵;
e、Hi经线性映射,得到主成分矩阵
Φi=Ri T·Hi (4)
f、设强工频基波和其主要谐波的数目为r,将Φi第一至第2r行所有元素置0,则得到重构的主成分矩阵Φ'i,令
Yi=Ri·Φi' (5)
则Yi为压制工频干扰后的Hankel矩阵,其具体形式可记为
<mrow>
<msub>
<mi>Y</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mi>L</mi>
</mtd>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mi>L</mi>
</mtd>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>M</mi>
</mtd>
<mtd>
<mi>M</mi>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mi>M</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mi>n</mi>
<mo>+</mo>
<mn>2</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mi>L</mi>
</mtd>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
g、定义xi'(l)=[yi(1),yi(2),···,yi(N)],则xi'(l)即为对应xi(l)的压制强工频噪声后的信号。
2.按照权利要求1所述的基于主成分分析的并行震源地震勘探资料强工频噪声压制方法,其特征在于,a步骤所述的识别含有强工频噪声的干扰道有两种情况:若干扰道只有一道,记干扰道信号为xi(l),其中i为并行地震勘探资料单炮地震记录的道号,l为采样点序号,l=1,2,···,N,N为干扰道信号的采样点数,执行一次步骤b~g;若干绕道为多个,对每一干扰道信号重复执行步骤b~g。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710395695.1A CN107219555B (zh) | 2017-05-31 | 2017-05-31 | 基于主成分分析的并行震源地震勘探资料强工频噪声压制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710395695.1A CN107219555B (zh) | 2017-05-31 | 2017-05-31 | 基于主成分分析的并行震源地震勘探资料强工频噪声压制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107219555A true CN107219555A (zh) | 2017-09-29 |
CN107219555B CN107219555B (zh) | 2018-09-14 |
Family
ID=59947019
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710395695.1A Active CN107219555B (zh) | 2017-05-31 | 2017-05-31 | 基于主成分分析的并行震源地震勘探资料强工频噪声压制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107219555B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108732624A (zh) * | 2018-05-29 | 2018-11-02 | 吉林大学 | 一种基于pca-emd的并行震源地震数据随机噪声压制方法 |
CN108957550A (zh) * | 2018-06-28 | 2018-12-07 | 吉林大学 | 基于svd-ica的tsp强工业电干扰压制方法 |
CN108957552A (zh) * | 2018-07-17 | 2018-12-07 | 吉林大学 | 基于ss-pca的地震数据海浪噪声压制方法 |
CN111290024A (zh) * | 2020-03-05 | 2020-06-16 | 吉林大学 | 一种svd自适应地震数据噪声压制方法 |
CN111854930A (zh) * | 2020-07-21 | 2020-10-30 | 长春理工大学 | 一种基于先验预估的振动信号工频干扰压制方法 |
CN112147692A (zh) * | 2019-06-28 | 2020-12-29 | 中国石油化工股份有限公司 | 用于海洋地震资料的二次震源衰减方法 |
CN112577752A (zh) * | 2020-12-25 | 2021-03-30 | 吉林大学 | 在城市复杂噪声环境中用地震检波器测量汽车转速的方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1999064896A1 (en) * | 1998-06-09 | 1999-12-16 | Geco As | Seismic data interpretation method |
CN102854533A (zh) * | 2011-07-01 | 2013-01-02 | 中国石油化工股份有限公司 | 一种基于波场分离原理提高地震资料信噪比的去噪方法 |
WO2013106720A1 (en) * | 2012-01-12 | 2013-07-18 | Schlumberger Canada Limited | Method for constrained history matching coupled with optimization |
CN103344989A (zh) * | 2013-07-11 | 2013-10-09 | 吉林大学 | 可控震源地震记录中脉冲噪声干扰的分析方法 |
CN104360393A (zh) * | 2014-10-20 | 2015-02-18 | 李闯 | 一种地震数据重建方法 |
US9417349B1 (en) * | 2013-02-13 | 2016-08-16 | Ihs Global Inc. | Picking faults in a seismic volume using a cost function |
CN106154330A (zh) * | 2016-06-20 | 2016-11-23 | 吉林大学 | 电磁式可控震源抛物线Radon变换压制并行震源串扰噪声方法 |
-
2017
- 2017-05-31 CN CN201710395695.1A patent/CN107219555B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1999064896A1 (en) * | 1998-06-09 | 1999-12-16 | Geco As | Seismic data interpretation method |
CN102854533A (zh) * | 2011-07-01 | 2013-01-02 | 中国石油化工股份有限公司 | 一种基于波场分离原理提高地震资料信噪比的去噪方法 |
WO2013106720A1 (en) * | 2012-01-12 | 2013-07-18 | Schlumberger Canada Limited | Method for constrained history matching coupled with optimization |
US9417349B1 (en) * | 2013-02-13 | 2016-08-16 | Ihs Global Inc. | Picking faults in a seismic volume using a cost function |
CN103344989A (zh) * | 2013-07-11 | 2013-10-09 | 吉林大学 | 可控震源地震记录中脉冲噪声干扰的分析方法 |
CN104360393A (zh) * | 2014-10-20 | 2015-02-18 | 李闯 | 一种地震数据重建方法 |
CN106154330A (zh) * | 2016-06-20 | 2016-11-23 | 吉林大学 | 电磁式可控震源抛物线Radon变换压制并行震源串扰噪声方法 |
Non-Patent Citations (3)
Title |
---|
周宪英 等: ""主成分分析法及其在数据降噪中的应用"", 《兵工自动化》 * |
王保利 等: ""工频干扰信号的时频域压制方法研究"", 《中国煤炭地质》 * |
钱广华 等: ""离心式压缩机转子振动信号的PCA降噪方法研究"", 《状态监测与分析》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108732624A (zh) * | 2018-05-29 | 2018-11-02 | 吉林大学 | 一种基于pca-emd的并行震源地震数据随机噪声压制方法 |
CN108957550A (zh) * | 2018-06-28 | 2018-12-07 | 吉林大学 | 基于svd-ica的tsp强工业电干扰压制方法 |
CN108957552A (zh) * | 2018-07-17 | 2018-12-07 | 吉林大学 | 基于ss-pca的地震数据海浪噪声压制方法 |
CN112147692A (zh) * | 2019-06-28 | 2020-12-29 | 中国石油化工股份有限公司 | 用于海洋地震资料的二次震源衰减方法 |
CN111290024A (zh) * | 2020-03-05 | 2020-06-16 | 吉林大学 | 一种svd自适应地震数据噪声压制方法 |
CN111854930A (zh) * | 2020-07-21 | 2020-10-30 | 长春理工大学 | 一种基于先验预估的振动信号工频干扰压制方法 |
CN111854930B (zh) * | 2020-07-21 | 2022-10-14 | 长春理工大学 | 一种基于先验预估的振动信号工频干扰压制方法 |
CN112577752A (zh) * | 2020-12-25 | 2021-03-30 | 吉林大学 | 在城市复杂噪声环境中用地震检波器测量汽车转速的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107219555B (zh) | 2018-09-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107219555B (zh) | 基于主成分分析的并行震源地震勘探资料强工频噪声压制方法 | |
US10677957B2 (en) | Method for random noise reduction from MRS oscillating signal using joint algorithms of EMD and TFPF | |
CN104459809B (zh) | 一种基于独立成分分析的全波核磁共振信号噪声滤除方法 | |
CN102221708B (zh) | 基于分数阶傅里叶变换的随机噪声压制方法 | |
CN102116868B (zh) | 一种地震波分解方法 | |
CN104345340B (zh) | 一种地震数据规则化中抖动稀疏采样方法 | |
CN107817527A (zh) | 基于块稀疏压缩感知的沙漠地震勘探随机噪声压制方法 | |
CN103926622B (zh) | 一种基于l1范数多道匹配滤波压制多次波的方法 | |
Baig et al. | Denoising seismic noise cross correlations | |
CN109164483A (zh) | 多分量地震数据矢量去噪方法及多分量地震数据矢量去噪装置 | |
CN105549076B (zh) | 一种基于交替方向法和全变分理论的地震数据处理方法 | |
CN104849756A (zh) | 一种提高地震数据分辨率增强有效弱信号能量的方法 | |
CN106680874A (zh) | 基于波形形态特征稀疏化建模的谐波噪声压制方法 | |
CN107179550B (zh) | 一种数据驱动的地震信号零相位反褶积方法 | |
CN103048684B (zh) | 一种多分量地震资料面波压制方法 | |
CN107144879A (zh) | 一种基于自适应滤波与小波变换结合的地震波降噪方法 | |
CN102262243B (zh) | 一种滤波法可控震源地震数据谐波干扰压制方法 | |
CN104345341A (zh) | 一种基于区域约束的分频段能量地震面波处理方法 | |
CN105259579A (zh) | 一种基于地震数据瞬时属性的强振幅屏蔽层剔除方法 | |
CN109885903A (zh) | 一种基于模型的地面核磁共振信号尖峰噪声去除方法 | |
CN104597502A (zh) | 一种新的石油地震勘探数据去噪方法 | |
CN110208856A (zh) | 一种基于流形分区2d-vmd的沙漠复杂噪声压制方法 | |
CN107436450A (zh) | 一种基于连续小波变换的地震信号带宽拓展方法 | |
CN104133248B (zh) | 一种高保真声波干扰压制方法 | |
CN103792574A (zh) | 一种储层频变含气性检测方法 |
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 |