发明名称 一种提高矿用微震检波器信噪比的方法
摘要 本发明公开了一种提高矿用微震检波器信噪比的方法,该方法包括:在矿用微震检波器信中的检波器敏感单元最大灵敏度方向上安装n+1个检波器敏感单元;分析同一方向上的n+1个通道信号的频谱相关性,提取有效信号,实现随机噪声的抑制;并通过设置单个检波器敏感单元的过采样率,利用随机噪声的不相关,对同一检波器敏感单元过采样数据进行均值降噪,从而提高信噪比。通过采用本发明公开的方法,提高了对能量、频率等参数信息均未知震源信号的拾取能力,且有效的抑制了随机噪声的干扰,提高了信噪比。
申请公布号 CN103852784B 申请公布日期 2016.12.07
申请号 CN201410090687.2 申请日期 2014.03.12
申请人 北京矿冶研究总院 发明人 杨小聪;冀虎;张达;刘冠洲
分类号 G01V1/18(2006.01)I 主分类号 G01V1/18(2006.01)I
代理机构 北京凯特来知识产权代理有限公司 11260 代理人 郑立明;赵镇勇
主权项 一种提高矿用微震检波器信噪比的方法,其特征在于,该方法包括:在矿用微震检波器中的检波器敏感单元最大灵敏度方向上安装n+1个检波器敏感单元;分析同一方向上的n+1个通道信号的频谱相关性,提取有效信号,实现随机噪声的抑制;并通过设置单个检波器敏感单元的过采样率,利用随机噪声的不相关,对同一检波器敏感单元过采样数据进行均值降噪,从而提高信噪比;其中,同一方向上的n+1个检波器敏感单元拾取振动信号在时域下叠加后的表达式为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>V</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>v</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>v</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>&Delta;</mi><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mo>...</mo><mo>+</mo><mi>v</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>n</mi><mi>&Delta;</mi><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mi>n</mi></munderover><mi>v</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>i</mi><mi>&Delta;</mi><mi>t</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0001009270130000011.GIF" wi="1179" he="134" /></maths>其中,Δt表示同一方向上两个相邻检波器敏感单元拾取振动信号的时差,v(t)表示输入信号;对表达式进行傅里叶变换,获得频域下检波器敏感单元拾取振动信号的方程:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>V</mi><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mi>n</mi></munderover><mi>v</mi><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><msup><mi>e</mi><mrow><mi>j</mi><mi>&omega;</mi><mi>i</mi><mi>&Delta;</mi><mi>t</mi></mrow></msup><mo>=</mo><mi>v</mi><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mi>n</mi></munderover><msup><mi>e</mi><mrow><mi>j</mi><mi>&omega;</mi><mi>i</mi><mi>&Delta;</mi><mi>t</mi></mrow></msup><mo>=</mo><mi>v</mi><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><msup><mi>e</mi><mrow><mo>-</mo><mi>j</mi><mi>&omega;</mi><mfrac><mi>n</mi><mn>2</mn></mfrac><mi>&Delta;</mi><mi>t</mi></mrow></msup><mfrac><mrow><mi>sin</mi><mrow><mo>(</mo><mrow><mi>&omega;</mi><mfrac><mrow><mrow><mo>(</mo><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow><mo>)</mo></mrow><mi>&Delta;</mi><mi>t</mi></mrow><mn>2</mn></mfrac></mrow><mo>)</mo></mrow></mrow><mrow><mi>sin</mi><mrow><mo>(</mo><mrow><mi>&omega;</mi><mfrac><mrow><mi>&Delta;</mi><mi>t</mi></mrow><mn>2</mn></mfrac></mrow><mo>)</mo></mrow></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0001009270130000012.GIF" wi="1437" he="158" /></maths>其中,v(ω)为v(t)在频域下的表达式,ω为傅里叶变换后的频域变量,j为时域转换到频域后的虚数单位;记<img file="FDA0001009270130000013.GIF" wi="661" he="151" />则:V(ω)=v(ω)·H(ω);表示n+1个检波器敏感单元其拾取振动信号的总振幅与信号的到达时间无关,仅与信号的频率和到达各检波器敏感单元的相对时差有关;所述分析同一方向上的n+1个通道信号的频谱相关性,提取有效信号,实现随机噪声的抑制包括:当同一方向上的n+1个检波器敏感单元进行信号叠加后,输出信号为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>V</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>v</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>v</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>&Delta;</mi><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mo>...</mo><mo>+</mo><mi>v</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>n</mi><mi>&Delta;</mi><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>i</mi><mi>&Delta;</mi><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><mi>n</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mi>s</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msqrt><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msqrt><mover><mrow><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>&OverBar;</mo></mover><mo>;</mo></mrow>]]></math><img file="FDA0001009270130000014.GIF" wi="1718" he="135" /></maths>其中,s(t)表示有效信号,u(t)表示随机噪声;输出信号取平均后,得:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mover><mrow><mi>V</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mrow><mi>V</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>=</mo><mi>s</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mfrac><mover><mrow><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>&OverBar;</mo></mover><msqrt><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msqrt></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0001009270130000021.GIF" wi="558" he="135" /></maths>所述通过设置单个检波器敏感单元的过采样率,利用随机噪声的不相关,对同一检波器敏感单元过采样数据进行均值降噪,从而提高信噪比包括:将所述单个检波器敏感单元的过采样率设为大于被测微震信号频率L倍,当同一个检波器敏感单元实现过采样时,由于过采样频率远大于被采集振动信号频率,其过采样输出信号为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>V</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>v</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>v</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>&Delta;</mi><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mo>...</mo><mo>+</mo><mi>v</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>m</mi><mi>&Delta;</mi><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mi>m</mi></munderover><mi>v</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>i</mi><mi>&Delta;</mi><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><mi>m</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mi>s</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msqrt><mrow><mi>m</mi><mo>+</mo><mn>1</mn></mrow></msqrt><mover><mrow><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>&OverBar;</mo></mover><mo>;</mo></mrow>]]></math><img file="FDA0001009270130000022.GIF" wi="1574" he="125" /></maths>其中,m+1表示单点过采样次数;输出信号取平均后,得:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mover><mrow><mi>V</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mrow><mi>V</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mrow><mi>m</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>=</mo><mi>s</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mfrac><mover><mrow><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>&OverBar;</mo></mover><msqrt><mrow><mi>m</mi><mo>+</mo><mn>1</mn></mrow></msqrt></mfrac><mo>.</mo></mrow>]]></math><img file="FDA0001009270130000023.GIF" wi="586" he="135" /></maths>
地址 100160 北京市丰台区南四环西路188号总部基地十八区23号楼