本发明专利技术涉及地球物理勘探提高地震勘探资料叠前处理信噪比的三维高密度地震勘探记录叠前面波压制方法，采集三维高密度地震数据，解编后按照炮线和检波线进行分选，得到多个正交子集记录，根据Nyquist定理得到空间频率的范围，从第1个正交子集开始，用傅氏变换的公式得到该正交子集记录的三维频率波数谱、FKK域压制面波的滤波器响应，将三维频率波数谱利用傅氏反变换的公式进行变换，得到压制面波后的时间域子集数据，进而完成整个三维高密度空间采样数据的叠前面波压制。本发明专利技术可以很好地在叠前压制面波，提高处理的效率和精度，有助于后续处理和成图。
【技术实现步骤摘要】
本专利技术涉及地球物理勘探技术，属于提高地震勘探资料叠前处理信噪比的三维高 密度地震勘探记录叠前面波压制方法。
技术介绍
在常规的地震勘探记录中，由于大的时间和空间采样信号，经过傅氏变换后，在频 率波数域会出现假频，这种现象导致在频率波数域内，信号和干扰互相干涉很难分离，对于 二维数据来说，这种干涉是面的干涉，对三维数据而言，则更加复杂，是体的干涉。这种情况 下的采样三维傅氏变换压制面波的去噪效果差强人意。高密度空间采样地震勘探是近几年来为了提高地震数据品质而发展起来的采集 和处理技术，它采用了单点接收、可控震源交替扫描、数字检波器等一系列技术。高密度勘 探数据具有时间、空间采样间隔小、覆盖次数高的特点，避免了空间假频的产生。三维高密 度地震的正交子集记录中，面波具有在某一时间切片上呈圆形分布的特征，而且面波在正 交子集的三维频率波数谱内呈圆锥形分布。目前对于三维高密度地震记录压制面波的处理，通常还是采用常规处理中的二维 滤波等方式按逐条测线进行，没有充分利用三维高密度地震的正交子集记录中，面波三维频率波数谱呈圆锥形分布的特征，所以，而提高处理的效率和精度都不高。
技术实现思路
本法明目的在于提供一种提高处理的效率和精度的三维高密度地震勘探记录叠 前面波压制方法。本专利技术采用以下技术步骤实现1)激发地震波，野外采集三维高密度地震数据，解编后按照炮线和检波线进行分 选，得到多个正交子集记录；2)根据子集数据的炮点间的距离Δχ和道间距Ay，根据Nyquist定理得到空间频率的范围权利要求1. 一种，其特征是采用以下技术步骤实现1)激发地震波，野外采集三维高密度地震数据，解编后按照炮线和检波线进行分选，得 到多个正交子集记录；2)根据子集数据的炮点间的距离ΔΧ和道间距Δy，根据Nyquist定理得到空间频率的范围全文摘要本专利技术涉及地球物理勘探提高地震勘探资料叠前处理信噪比的，采集三维高密度地震数据，解编后按照炮线和检波线进行分选，得到多个正交子集记录，根据Nyquist定理得到空间频率的范围，从第1个正交子集开始，用傅氏变换的公式得到该正交子集记录的三维频率波数谱、FKK域压制面波的滤波器响应，将三维频率波数谱利用傅氏反变换的公式进行变换，得到压制面波后的时间域子集数据，进而完成整个三维高密度空间采样数据的叠前面波压制。本专利技术可以很好地在叠前压制面波，提高处理的效率和精度，有助于后续处理和成图。文档编号G01V1/28GK102103215SQ20091024273公开日2011年6月22日 申请日期2009年12月16日 优先权日2009年12月16日专利技术者詹毅, 赵波 申请人:中国石油天然气集团公司, 中国石油集团东方地球物理勘探有限责任公司本文档来自技高网...
【技术保护点】
１．一种三维高密度地震勘探记录叠前面波压制方法，其特征是采用以下技术步骤实现：１）激发地震波，野外采集三维高密度地震数据，解编后按照炮线和检波线进行分选，得到多个正交子集记录；２）根据子集数据的炮点间的距离Δｘ和道间距Δｙ，根据Ｎｙｑｕｉｓｔ定理得到空间频率的范围（ｍａｔｈ）??（ｍｒｏｗ）?（ｍｉ）ｋｘ（／ｍｉ）?（ｍｏ）＝（／ｍｏ）?（ｍｎ）０（／ｍｎ）?（ｍｏ）～（／ｍｏ）?（ｍｆｒａｃ）?（ｍｎ）１（／ｍｎ）?（ｍｒｏｗ）?（ｍｎ）２（／ｍｎ）?（ｍｏ）＊（／ｍｏ）?（ｍｉ）＆Ｄｅｌｔａ；ｘ（／ｍｉ）?（／ｍｒｏｗ）?（／ｍｆｒａｃ）?（ｍｏ），（／ｍｏ）?（／ｍｒｏｗ）?（／ｍａｔｈ）（ｍａｔｈ）??（ｍｒｏｗ）?（ｍｉ）ｋｙ（／ｍｉ）?（ｍｏ）＝（／ｍｏ）?（ｍｎ）０（／ｍｎ）?（ｍｏ）～（／ｍｏ）?（ｍｆｒａｃ）?（ｍｎ）１（／ｍｎ）?（ｍｒｏｗ）?（ｍｎ）２（／ｍｎ）?（ｍｏ）＊（／ｍｏ）?（ｍｉ）＆Ｄｅｌｔａ；ｙ（／ｍｉ）?（／ｍｒｏｗ）?（／ｍｆｒａｃ）?（ｍｏ），（／ｍｏ）?（ｍｏ）；（／ｍｏ）?（／ｍｒｏｗ）?（／ｍａｔｈ）３）根据子集数据的时间采样率Δｔ，根据Ｎｙｑｕｉｓｔ定理得到时间方向的频率范围（ｍａｔｈ）??（ｍｒｏｗ）?（ｍｉ）ｗ（／ｍｉ）?（ｍｏ）＝（／ｍｏ）?（ｍｆｒａｃ）?（ｍｎ）１（／ｍｎ）?（ｍｒｏｗ）?（ｍｎ）２（／ｍｎ）?（ｍｏ）＊（／ｍｏ）?（ｍｉ）＆Ｄｅｌｔａ；ｔ（／ｍｉ）?（／ｍｒｏｗ）?（／ｍｆｒａｃ）?（ｍｏ）；（／ｍｏ）?（／ｍｒｏｗ）?（／ｍａｔｈ）４）然后根据炮点数量ｍ、检波点数量ｎ和每一道的样点数ｐ，得到子集数据中每一个点的空间波数值和频率值ｉ为该点所在的炮数，ｊ为该点所在的道数，ｋ为该点所在的样点数；５）通过叠前记录中面波传播的距离Ｘ和传播该距离所需时间ｔ的比值得到面波的视速度ＶＬ；６）从第１个正交子集开始，利用以下傅氏变换的公式进行三维傅氏变换，得到该正交子集记录的三维频率波数谱；（ｍａｔｈ）??（ｍｒｏｗ）?（ｍｉ）Ｆ（／ｍｉ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｉ）ｗ（／ｍｉ）?（ｍｏ），（／ｍｏ）?（ｍｓｕｂ）?（ｍｉ）ｋ（／ｍｉ）?（ｍｉ）ｘ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｏ），（／ｍｏ）?（ｍｓｕｂ）?（ｍｉ）ｋ（／ｍｉ）?（ｍｉ）ｙ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（ｍｏ）＝（／ｍｏ）?（ｍｕｎｄｅｒｏｖｅｒ）?（ｍｏ）＆Ｉｎｔｅｇｒａｌ；（／ｍｏ）?（ｍｒｏｗ）?（ｍｏ）－（／ｍｏ）?（ｍｏ）＆ｉｎｆｉｎ；（／ｍｏ）?（／ｍｒｏｗ）?（ｍｏ）＆ｉｎｆｉｎ；（／ｍｏ）?（／ｍｕｎｄｅｒｏｖｅｒ）?（ｍｕｎｄｅｒｏｖｅｒ）?（ｍｏ）＆Ｉｎｔｅｇｒａｌ；（／ｍｏ）?（ｍｒｏｗ）?（ｍｏ）－（／ｍｏ）?（ｍｏ）＆ｉｎｆｉｎ；（／ｍｏ）?（／ｍｒｏｗ）?（ｍｏ）＆ｉｎｆｉｎ；（／ｍｏ）?（／ｍｕｎｄｅｒｏｖｅｒ）?（ｍｕｎｄｅｒｏｖｅｒ）?（ｍｏ）＆Ｉｎｔｅｇｒａｌ；（／ｍｏ）?（ｍｒｏｗ）?（ｍｏ）－（／ｍｏ）?（ｍｏ）＆ｉｎｆｉｎ；（／ｍｏ）?（／ｍｒｏｗ）?（ｍｏ）＆ｉｎｆｉｎ；（／ｍｏ）?（／ｍｕｎｄｅｒｏｖｅｒ）?（ｍｉ）Ｘ（／ｍｉ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｉ）ｔ（／ｍｉ）?（ｍｏ），（／ｍｏ）?（ｍｉ）ｘ（／ｍｉ）?（ｍｏ），（／ｍｏ）?（ｍｉ）ｙ（／ｍｉ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（ｍｓｕｐ）?（ｍｉ）ｅ（／ｍｉ）?（ｍｒｏｗ）?（ｍｏ）－（／ｍｏ）?（ｍｉ）ｉ（／ｍｉ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｉ）ｗｔ（／ｍｉ）?（ｍｏ）＋（／ｍｏ）?（ｍｓｕｂ）?（ｍｉ）ｋ（／ｍｉ）?（ｍｉ）ｘ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｉ）ｘ（／ｍｉ）?（ｍｏ）＋（／ｍｏ）?（ｍｓｕｂ）?（ｍｉ）ｋ（／ｍｉ）?（ｍｉ）ｙ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｉ）ｙ（／ｍｉ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（／ｍｒｏｗ）?（／ｍｓｕｐ）?（ｍｉ）ｄｔｄｘｄｙ（／ｍｉ）?（ｍｏ）－（／ｍｏ）?（ｍｏ）－（／ｍｏ）?（ｍｏ）－（／ｍｏ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｎ）１（／ｍｎ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（／ｍｒｏｗ）?（／ｍａｔｈ）式中Ｆ（ｗ，ｋｘ，ｋｙ）是正交子集记录的三维频率波数谱，Ｘ（ｔ，ｘ，ｙ）是正交子集地震记录，ｗ是时间方向的频率，ｋｘ是沿炮线方向的空间频率，ｋｙ是沿检波线方向的空间频率；７）根据上述４）和５）得到的数值，带入以下公式，得到ＦＫＫ域压制面波的滤波器响应Ｆ′（ｗ，ｋｘ，ｋｙ）；（ｍａｔｈ）??（ｍｒｏｗ）?（ｍｓｕｐ）?（ｍｉ）Ｆ（／ｍｉ）?（ｍｏ）＆ｐｒｉｍｅ；（／ｍｏ）?（／ｍｓｕｐ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（...
【技术特征摘要】
【专利技术属性】
技术研发人员：詹毅，赵波，
申请(专利权)人：中国石油天然气集团公司，中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：11
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。