主权项 |
一种基于波束定向原理的单炮随机噪声压制方法,其特征在于,包括以下步骤:a、输入一炮地震数据X,道数为n,将X数据复制得到X<sub>1</sub>、X<sub>2</sub>、X<sub>3</sub>、……、X<sub>q‑2</sub>、X<sub>q‑1</sub>、X<sub>q</sub>,q≥2且为整数,q取奇数;b、利用延时基数τ,对1~q炮数据做延时,延时参数分别为τ<sub>1</sub>、τ<sub>2</sub>、τ<sub>3</sub>、…、τ<sub>k</sub>、…τ<sub>q</sub>,k=1,2,3,...,q;令<img file="FDA0000922087450000011.GIF" wi="423" he="149" />τ<sub>(k+1)</sub>=τ<sub>k</sub>+τ,则τ<sub>(q+1)/2</sub>=0,<math><![CDATA[<mrow><msub><mi>τ</mi><mi>q</mi></msub><mo>=</mo><mfrac><mrow><mo>(</mo><mi>q</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mn>2</mn></mfrac><mo>·</mo><mi>τ</mi><mo>;</mo></mrow>]]></math><img file="FDA00009220874500000112.GIF" wi="381" he="150" /></maths>c、将X<sub>1</sub>数据延时τ<sub>1</sub>,得到Y<sub>1</sub>=X<sub>1</sub>(t+τ<sub>1</sub>);将X<sub>2</sub>数据延时τ<sub>2</sub>,得到Y<sub>2</sub>=X<sub>2</sub>(t+τ<sub>2</sub>);以此类推,将第X<sub>k</sub>炮数据延时τ<sub>k</sub>,得到Y<sub>k</sub>=X<sub>k</sub>(t+τ<sub>k</sub>),最后,将第X<sub>q</sub>炮数据延时τ<sub>q</sub>,得到Y<sub>q</sub>=X<sub>q</sub>(t+τ<sub>q</sub>);d、对延时后的数据Y<sub>1</sub>、Y<sub>2</sub>、Y<sub>3</sub>、…、Y<sub>q</sub>进行有效道分选,分选道基数为λ,λ取正整数,<img file="FDA0000922087450000013.GIF" wi="534" he="87" />1~q炮延时后的数据分选道数为λ<sub>1</sub>、λ<sub>2</sub>、λ<sub>3</sub>、…、λ<sub>k</sub>、…λ<sub>q</sub>,k=1,2,3,…,q;令<img file="FDA0000922087450000014.GIF" wi="430" he="151" />λ<sub>k+1</sub>=λ<sub>k</sub>+λ,则λ<sub>(q+1)/2</sub>=0,<img file="FDA00009220874500000113.GIF" wi="389" he="143" />e、Y<sub>(q+1)/2</sub>为中心炮数据,λ<sub>(q+1)/2</sub>=0,保持Y<sub>(q+1)/2</sub>数据不变;<math><![CDATA[<mrow><mi>i</mi><mo>∈</mo><mo>[</mo><mn>1</mn><mo>,</mo><mrow><mo>(</mo><mfrac><mrow><mo>(</mo><mi>q</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mn>2</mn></mfrac><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>]</mo><mo>,</mo></mrow>]]></math><img file="FDA0000922087450000016.GIF" wi="510" he="172" /></maths>令<math><![CDATA[<mrow><msub><mi>Y</mi><mi>i</mi></msub><mo>=</mo><mrow><mo>(</mo><msubsup><mi>y</mi><mn>1</mn><mi>i</mi></msubsup><mo>,</mo><msubsup><mi>y</mi><mn>2</mn><mi>i</mi></msubsup><mo>,</mo><mo>...</mo><mo>,</mo><msubsup><mi>y</mi><mi>n</mi><mi>i</mi></msubsup><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000922087450000017.GIF" wi="518" he="95" /></maths>分选后得到数据记为<img file="FDA0000922087450000018.GIF" wi="649" he="159" /><math><![CDATA[<mrow><mi>j</mi><mo>∈</mo><mo>[</mo><mrow><mo>(</mo><mfrac><mrow><mo>(</mo><mi>q</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mn>2</mn></mfrac><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><mi>q</mi><mo>]</mo><mo>,</mo></mrow>]]></math><img file="FDA0000922087450000019.GIF" wi="531" he="173" /></maths>令<math><![CDATA[<mrow><msub><mi>Y</mi><mi>j</mi></msub><mo>=</mo><mrow><mo>(</mo><msubsup><mi>y</mi><mn>1</mn><mi>j</mi></msubsup><mo>,</mo><msubsup><mi>y</mi><mn>2</mn><mi>j</mi></msubsup><mo>,</mo><mo>...</mo><mo>,</mo><msubsup><mi>y</mi><mi>n</mi><mi>j</mi></msubsup><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA00009220874500000110.GIF" wi="526" he="95" /></maths>分选后得到数据记为<img file="FDA00009220874500000111.GIF" wi="671" he="165" />f、以λ<sub>1</sub>分选Y<sub>1</sub>数据,得到Z<sub>1</sub>=Y<sub>1</sub>(x+λ<sub>1</sub>),即以λ<sub>2</sub>分选Y<sub>2</sub>数据,得到Z<sub>2</sub>=Y<sub>2</sub>(x+λ<sub>2</sub>);以此类推,以λ<sub>k</sub>分选Y<sub>k</sub>数据,得到Z<sub>k</sub>=Y<sub>k</sub>(x+λ<sub>k</sub>);最后,以λ<sub>q</sub>分选Y<sub>q</sub>数据,得到Z<sub>q</sub>=Y<sub>q</sub>(x+λ<sub>q</sub>);g、根据Z<sub>1</sub>~Z<sub>q</sub>的结果,将波束定向后畸变的数据置0,得到Z<sub>1</sub>'~Z<sub>q</sub>';h、将Z<sub>1</sub>'~Z<sub>q</sub>'的数据做线性叠加,合成针对延时基数τ和分选道基数为λ的数据W<sub>τ,λ</sub>;i、采用新的延时基数τ'和分选道基数为λ',按照b~h的步骤进行处理得到针对延时基数τ'和分选道基数为λ'的数据W<sub>τ',λ'</sub>;j、最后,将所有对应延时基数和分选道基数的结果做线性叠加,得到最终随机噪声压制后的数据。 |