发明名称 地震弱信号增强处理方法
摘要 本发明提供一种地震弱信号增强处理方法,该地震弱信号增强处理方法包括:步骤a,由目标地震道建立初始围岩反射<img file="14450dest_path_image001.GIF" wi="38" he="25" />;步骤b,根据该初始围岩反射<img file="242169dest_path_image001.GIF" wi="38" he="25" />计算新的围岩反射<img file="808280dest_path_image002.GIF" wi="32" he="22" />;步骤c,判断该新的围岩反射<img file="906686dest_path_image002.GIF" wi="32" he="22" />是否稳定,当该新的围岩反射<img file="329577dest_path_image002.GIF" wi="32" he="22" />稳定时,输出该新的围岩反射<img file="451116dest_path_image002.GIF" wi="32" he="22" />作为围岩信号;以及步骤d,从地震道的地震数据中减去该围岩信号,得到弱信号。该地震弱信号增强处理方法解决了现有方法难以识别与解释储层目标弱信号的问题,可以对检测到的地震信号进行盲分离,对地震叠合干涉的复杂反射信号,根据输入源信号的统计特性,提取所需要特征信号或恢复出源信号各个独立成分。
申请公布号 CN103543467B 申请公布日期 2016.04.27
申请号 CN201210247077.X 申请日期 2012.07.17
申请人 中国石油化工股份有限公司;中国石油化工股份有限公司胜利油田分公司地质科学研究院 发明人 穆星;刘显太;王永诗;王长江;徐希坤;段鸿杰;赵铭海;杨培杰;王庆华;董立生;李敏;贾玉茹;王学军;张景涛;李秀国
分类号 G01V1/28(2006.01)I 主分类号 G01V1/28(2006.01)I
代理机构 济南日新专利代理事务所 37224 代理人 谢省法
主权项 地震弱信号增强处理方法,其特征在于,该地震弱信号增强处理方法包括:步骤a,建立盲源信号混叠模型:<maths num="0001"><math><![CDATA[<mrow><mi>x</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>A</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><msub><mi>F</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>-</mo><msub><mi>t</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>s</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000928387580000011.GIF" wi="656" he="135" /></maths>i=1,2,…,m式中x(i,t)代表某一地震道的地震数据;F<sub>k</sub>代表第k个围岩信号源;A<sub>ik</sub>代表F<sub>k</sub>的振幅大小;t<sub>ik</sub>代表时间延迟;s<sub>i</sub>(t)代表该地震道中的弱信号;m表示地震道数,n表示地震信号的长度;步骤b,根据该盲源信号混叠模型,以中间地震道为研究目标,地震弱信号分离问题的数学表达式为:<maths num="0002"><math><![CDATA[<mrow><msup><mi>y</mi><mo>*</mo></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>y</mi><mrow><mo>(</mo><mfrac><msub><mi>m</mi><mi>e</mi></msub><mn>2</mn></mfrac><mo>+</mo><mn>1</mn><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>A</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><msub><mi>F</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>-</mo><msub><mi>t</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000928387580000012.GIF" wi="796" he="143" /></maths>式中:y<sup>*</sup>(t)为待求的围岩信号,m<sub>e</sub>为相邻地震道的个数,取偶数;步骤c,由目标地震道建立初始围岩反射y<sub>0</sub>(t);步骤d,根据该初始围岩反射y<sub>0</sub>(t)计算新的围岩反射y(t);步骤e,判断该新的围岩反射y(t)是否稳定,当该新的围岩反射y(t)稳定时,输出该新的围岩反射y(t)作为围岩信号,当该新的围岩反射y(t)不稳定时,将该新的围岩反射y(t)作为新的初始围岩反射y<sub>0</sub>(t),并返回步骤d;以及步骤f,从地震道的地震数据中减去该围岩信号,得到弱信号;所述步骤d包括:分别计算各地震道与该初始围岩反射y<sub>0</sub>(t)的互相关矩阵R<sub>xy</sub>,根据互相关计算由各地震道预测y<sup>*</sup>(t)的预测因子;将预测因子作用于地震道,获得各地震道计算出的围岩反射y<sub>i</sub>(t);以及将多个y<sub>i</sub>(t)平均,获得该新的围岩反射y(t)。
地址 100728 北京市朝阳区朝阳门北大街22号