主权项 |
1、一种压制低信噪比地震记录中随机噪声的方法,包括以下处理步骤:(1)采集地震数据,采用傅利叶变换将多道地震数据从时间域变换到频率域,形成空间域的矢量数据;(2)在频率域的单个频率点采用常规F-X域预测滤波方法求取预测算子,(3)在频率域,用单点频率预测算子对不同地震道相同频率点的数据进行噪声衰减,消除无规律分布的随机噪音;(4)把消除随机噪音的数据作为有效信号,选定频带的宽度,在整个频率域内分别统计去噪前后不同地震道相同频带宽度的振幅能量,利用去噪前后相同频带宽度的振幅能量比值求取信噪比,信噪比值最大的频带即较高信噪比频带被确认为优势频带;(5)通过对优势频带相应的预测算子分量进行统计,根据优势频带与非优势频带预测算子分量的关联性,推算出非优势频带的预测算子分量,进而求得非优势频带的预测算子;优势频带与非优势频带预测算子分量关联系数为Q,则有:<maths num="001"><![CDATA[ <math> <mrow> <msub> <mi>F</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>m</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>kΔf</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>Σ</mi> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>C</mi> <mrow> <mi>L</mi> <mo>-</mo> <mi>j</mi> </mrow> </msub> </munderover> <msub> <mi>Q</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>m</mi> </mrow> </msub> <msub> <mi>F</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>m</mi> </mrow> </msub> <mo></mo> <mrow> <mo>[</mo> <mrow> <mo>(</mo> <mi>k</mi> <mo>±</mo> <mi>m</mi> <mo>)</mo> </mrow> <mi>Δf</mi> <mo>]</mo> </mrow> <mo></mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow> </math> ]]></maths>(j=1,2,…,L)式(4)中:F<sub>j,m</sub>为某一频率的预测算子分量,用通常的维纳滤波方法求得,C为组合符号,<maths num="002"><![CDATA[ <math> <mrow> <msub> <mi>C</mi> <mrow> <mi>L</mi> <mo>-</mo> <mi>j</mi> </mrow> </msub> <mo>=</mo> <msubsup> <mi>C</mi> <mi>L</mi> <mrow> <mi>L</mi> <mo>-</mo> <mi>j</mi> </mrow> </msubsup> <mo>;</mo> </mrow> </math> ]]></maths>Δf为频率采样间隔;k为优势频段的带宽,为算子外推长度;L为算子长度;式中+表示向低频外推,-表示向高频外推,关联系数为Q从下式求得:<maths num="003"><![CDATA[ <math> <mrow> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>Q</mi> <mrow> <mi>j</mi> <mo>,</mo> <mn>1</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>Q</mi> <mrow> <mi>j</mi> <mo>,</mo> <mn>2</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>Q</mi> <mrow> <mi>j</mi> <mo>,</mo> <mn>3</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <mo>·</mo> <mo>·</mo> <mo>·</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>Q</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <msup> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>F</mi> <mrow> <mi>j</mi> <mo>,</mo> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>F</mi> <mrow> <mi>j</mi> <mo>,</mo> <mn>2</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>F</mi> <mrow> <mi>j</mi> <mo>,</mo> <mn>3</mn> </mrow> </msub> </mtd> <mtd> <mo>·</mo> <mo>·</mo> <mo>·</mo> </mtd> <mtd> <msub> <mi>F</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>F</mi> <mrow> <mi>j</mi> <mo>,</mo> <mn>2</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>F</mi> <mrow> <mi>j</mi> <mo>,</mo> <mn>3</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>F</mi> <mrow> <mi>j</mi> <mo>,</mo> <mn>4</mn> </mrow> </msub> </mtd> <mtd> <mo>·</mo> <mo>·</mo> <mo>·</mo> </mtd> <mtd> <msub> <mi>F</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>k</mi> <mo>+</mo> <mn>2</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>F</mi> <mrow> <mi>j</mi> <mo>,</mo> <mn>3</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>F</mi> <mrow> <mi>j</mi> <mo>,</mo> <mn>4</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>F</mi> <mrow> <mi>j</mi> <mo>,</mo> <mn>5</mn> </mrow> </msub> </mtd> <mtd> <mo>·</mo> <mo>·</mo> <mo>·</mo> </mtd> <mtd> <msub> <mi>F</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>k</mi> <mo>+</mo> <mn>3</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <mo>·</mo> <mo>·</mo> <mo>·</mo> </mtd> <mtd> <mo>·</mo> <mo>·</mo> <mo>·</mo> </mtd> <mtd> <mo>·</mo> <mo>·</mo> <mo>·</mo> </mtd> <mtd> <mo>·</mo> <mo>·</mo> <mo>·</mo> </mtd> <mtd> <mo>·</mo> <mo>·</mo> <mo>·</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>F</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>m</mi> </mrow> </msub> </mtd> <mtd> <mo>·</mo> <mo>·</mo> <mo>·</mo> </mtd> <mtd> <mo>·</mo> <mo>·</mo> <mo>·</mo> </mtd> <mtd> <mo>·</mo> <mo>·</mo> <mo>·</mo> </mtd> <mtd> <msub> <mi>F</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>m</mi> <mo>+</mo> <mi>k</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>·</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>F</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>F</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>k</mi> <mo>+</mo> <mn>2</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>F</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>k</mi> <mo>+</mo> <mn>3</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <mo>·</mo> <mo>·</mo> <mo>·</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>F</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>k</mi> <mo>+</mo> <mi>m</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow> </math> ]]></maths>(6)完成优势频带外推随机噪音压制处理;(7)将地震数据从频率域转换到时间域,得到消除随机噪声的有效信号,生成地震剖面图。 |