发明名称 提高灰岩出露区地震数据分辨率的分频匹配方法
摘要 本发明公开了一种提高灰岩出露区地震数据分辨率的分频匹配方法。它是考虑利用高灵敏度检波器和常规检波器按同样的条件进行观测。然后,将同等条件下采集的地震数据以20Hz为界分别进行分频处理。之后,将高灵敏度检波器对应的大于20Hz的地震数据作为期望输出,求出匹配滤波器,进而对所有常规检波器采集的大于20Hz的地震数据进行匹配滤波。最后,将常规检波器采集的小于等于20Hz的地震数据和大于20Hz的匹配后的地震数据进行数据重构,就获得了分频匹配的地震数据。本发明使得常规检波器对应地震数据高频成分得到了加强,并且低频成分得到了保持,提高了灰岩出露区地震数据的分辨率。
申请公布号 CN102323616B 申请公布日期 2013.03.13
申请号 CN201110152185.4 申请日期 2011.06.08
申请人 浙江大学 发明人 石战结;田钢;沈洪垒
分类号 G01V1/28(2006.01)I 主分类号 G01V1/28(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 张法高
主权项 1.一种提高灰岩出露区地震数据分辨率的分频匹配方法,其特征在于它的步骤如下:1)利用高灵敏度检波器和常规检波器进行地震数据的对比接收:将野外地震排列某一个接收点的常规检波器换成高灵敏度检波器接到大线上,并将高灵敏度检波器埋置在相邻的一个接收点的位置上,进行数据采集;2)分频处理:设灰岩地表某点实际振动信号为<img file="2011101521854100001DEST_PATH_IMAGE001.GIF" wi="36" he="25" />,常规检波器、高灵敏度检波器的系统响应分别为<img file="2011101521854100001DEST_PATH_IMAGE002.GIF" wi="29" he="22" />、<img file="2011101521854100001DEST_PATH_IMAGE003.GIF" wi="32" he="22" />,若忽略耦合系统差异的影响,则以上检波器实际接收的信号分别为,<img file="2011101521854100001DEST_PATH_IMAGE004.GIF" wi="121" he="25" />(1)<img file="2011101521854100001DEST_PATH_IMAGE005.GIF" wi="125" he="26" />(2)采用余弦函数在频率域构造滤波函数,以20Hz为界,对同等条件下高灵敏度检波器和常规检波器采集的地震数据进行分频处理;若地震数据的采样率为1ms,那么将会获得0-20Hz和20-500Hz两个频段的地震数据,分频后的地震数据表示为,<img file="2011101521854100001DEST_PATH_IMAGE006.GIF" wi="202" he="25" />(3)<img file="2011101521854100001DEST_PATH_IMAGE007.GIF" wi="206" he="26" />(4)其中<img file="2011101521854100001DEST_PATH_IMAGE008.GIF" wi="20" he="25" />为第<img file="2011101521854100001DEST_PATH_IMAGE009.GIF" wi="14" he="20" />个频带(<img file="2011101521854100001DEST_PATH_IMAGE010.GIF" wi="49" he="22" />),<img file="2011101521854100001DEST_PATH_IMAGE011.GIF" wi="18" he="24" />对应0-20Hz,<img file="2011101521854100001DEST_PATH_IMAGE012.GIF" wi="20" he="24" />对应20-500Hz,<img file="2011101521854100001DEST_PATH_IMAGE013.GIF" wi="54" he="25" />为余弦分频滤波函数,<img file="2011101521854100001DEST_PATH_IMAGE014.GIF" wi="58" he="25" />、<img file="2011101521854100001DEST_PATH_IMAGE015.GIF" wi="60" he="26" />分别为常规检波器、高灵敏度检波器对应的分频后的地震数据;3)匹配处理:在20-500Hz频率范围内,利用算子<img file="2011101521854100001DEST_PATH_IMAGE016.GIF" wi="30" he="22" />表示两种检波器的系统响应差异,得出,<img file="2011101521854100001DEST_PATH_IMAGE017.GIF" wi="166" he="26" />(5)利用最小二乘原理,将式(5)构建为维纳方程,<img file="2011101521854100001DEST_PATH_IMAGE018.GIF" wi="50" he="20" />(6)其中,<img file="2011101521854100001DEST_PATH_IMAGE019.GIF" wi="17" he="18" />是一个矩阵,为<img file="2011101521854100001DEST_PATH_IMAGE020.GIF" wi="58" he="25" />的自相关函数; h=(h(t<sub>1</sub>),h(t<sub>2</sub>),…h(t<sub>n</sub>))是系统差异算子;<img file="2011101521854100001DEST_PATH_IMAGE021.GIF" wi="110" he="25" />是<img file="811596DEST_PATH_IMAGE020.GIF" wi="58" he="25" />和<img file="2011101521854100001DEST_PATH_IMAGE022.GIF" wi="60" he="26" />的互相关函数;利用式(6)求出20-500Hz频率范围内的检波器系统响应差异算子,将其称为匹配滤波器,然后,利用<img file="604103DEST_PATH_IMAGE016.GIF" wi="30" he="22" />对同频段的其它常规检波器采集的地震数据进行匹配处理,即可得到匹配处理后的地震数据,表示为,<img file="DEST_PATH_IMAGE023.GIF" wi="168" he="28" />(7)其中,<img file="2011101521854100001DEST_PATH_IMAGE024.GIF" wi="61" he="28" />为利用常规检波器采集的地震数据,<img file="DEST_PATH_IMAGE025.GIF" wi="58" he="25" />匹配后的地震数据;4)数据重构:对0-20Hz频率范围内的常规检波器采集的地震数据与20-500Hz频率范围内常规检波器对应的匹配后的地震数据进行重构,即得到分频匹配处理后的地震数据<img file="2011101521854100001DEST_PATH_IMAGE026.GIF" wi="30" he="22" />,表示为,<img file="DEST_PATH_IMAGE027.GIF" wi="206" he="29" />(8)。
地址 310027 浙江省杭州市西湖区浙大路38号