准噶尔盆地石炭系主要是一套中深层火山岩(图 1),地震响应特征以低频为主(2800ms~3600ms)。受强能量面波等影响,信噪比较低。其上覆侏罗系发育一套煤层,对下伏地层反射有很强的屏蔽作用,使石炭系有效波高频成分迅速衰减,有效波的频率偏低。由于石炭系构造较复杂,地层速度横向变化大,造成石炭系顶界反射波难以成像(图 2红色虚线处)。目前针对石炭系勘探目标,地震采集技术已经实现从拓展高频到同时拓展高、低频的转变,随着低频可控震源采集技术的应用,中深层地震数据的品质得到明显提升[1-7]。
在地震数据处理方面,传统的叠前面波压制处理会造成有效波低频信息的损失,本文研究选择保护低频的十字排列域噪声压制方法改善石炭系资料成像品质。同时自主研发了层控低频增强处理方法,大大提高目的层段地震资料的信噪比。
2 方法原理 2.1 十字排列域FKK噪声压制三维地震采集中,地震数据按文件号记录方式进行记录,即一个单炮数据对应一个文件号,一次炮点激发有多个检波点接收,以单条炮线与区块内所有相关检波点对应的地震道进行分选,形成的子集是十字交叉子集,对该子集进行三维傅里叶变换,利用傅里叶变换域中信号和规则干扰可以分离的特点,进行视速度滤波处理,达到叠前真三维压制规则干扰的效果[8-11]。图 3是在炮域(图 3a)和十字排列域(图 3b)中噪声压制效果对比,从单炮去噪效果可见,十字排列域能有效区分低频噪声与一次波低频端信息,在提高信噪比的同时,最大程度保留有效反射的低频成分。
首先通过快速傅里叶正变换将地震数据从时间域变换到频率域,然后设计一个频率域低频增强因子,将其与频率域地震数据相乘,再通过快速傅里叶反变换到时间域,即得到低频增强后的地震数据,具体流程如图 4所示。
设计的低频增强因子频率与强度关系如图 5所示。f1、f2和A是根据地震资料直接给定的参数,该低频增强因子的作用是:大于f2的强度为1,对应振幅保持不变,小于f1的强度为A,对应振幅增强A倍,f2和f1之间由增强1倍逐渐过渡到增强A倍。
对空间上剧烈变化的目的层段进行低频增强处理,需要运用时间随空间变化的插值处理技术,本文针对三维地震数据采用了二元插值方法[12-15]。
假设给定的矩形域上n×m个节点在x方向与y方向上的坐标分别为
x0<x1<…<xn-1与y0<y1<…<yn-1,对应的函数值为
$z_{k j}=z\left(x_{k}, y_{j}\right) k=0, 1, \cdots, n-1 ; j=0, 1, \cdots, m-1 $ | (1) |
在两个方向上,以(u, v)为中心插值点,左右各取4个坐标,即
$ \begin{array}{c} x_{p}<x_{p+1}<x_{p+2}<x_{p+3} <u<x_{p+4}<x_{p+5}< \\ x_{p+6}<x_{p+7} \end{array} $ | (2) |
$ y_{q} <y_{q+1} <y_{q+2}<y_{q+3}<v<y_{q+4}<y_{q+5}< \\ y_{q+6}<y_{q+7} $ | (3) |
再用二元插值公式
$ z(u, v)=\sum\limits_{i=p}^{p+7} \sum\limits_{j=q}^{q+7}\left(\prod\limits_{k=p \atop k \neq i}^{p+7} \frac{u-x_{k}}{x_{i}-x_{k}}\right)\left(\prod\limits_{l=q \atop l \neq j}^{q+7} \frac{v-y_{l}}{y_{j}-y_{l}}\right) z_{i j} $ | (4) |
计算插值点(u,v)处的函数近似值。
二元全区间插值技术能够对三维地震资料进行中深层空变低频增强,从而提高中深层地震数据信噪比及成像质量(图 6)。
上述低频补偿和低频保护去噪处理技术被应用到准噶尔盆地滴南地区资料处理中,处理成果显示,石炭系成像质量得到明显提升,石炭系内幕断裂体系展布特征清晰(图 7b)。基于新资料的解释在该区石炭系发现多个有利勘探目标,部署的多口探井均获得成功。
本文总结出一套针对准噶尔盆地石炭系火山岩低频保护的处理技术流程,将该流程应用于新疆准噶尔盆地多个目标区块,地震资料石炭系成像品质得到明显提升,取得如下结论。
(1) 通过以十字排列域滤波为主的噪声压制叠前配套技术的应用,能够在压制噪声的同时较好地保护有效低频信息,为石炭系成像提供基础。
(2) 二元全区间插值技术使低频能量补偿技术有效应用到地震层位起伏变化剧烈的地震数据,并且可以应用到不同地震目标的成像处理中。
(3) 层控频率域低频补偿技术在新疆准噶尔盆地多块目标区推广应用并取得效果,该技术是提高地震资料深层信噪比卓有成效的特色技术。
[1] |
陈新发, 匡立春, 查明. 火山岩油气成藏机理与勘探技术——以准噶尔盆地为例[M]. 北京: 科学出版社, 2014. CHEN Xinfa, KUANG Lichun, ZHA ming. Volcanic Oil and Gas Accumulation Mechanism and Explora-tion Technology in the Junggar Basin as An Example[M]. Beijing: Science Press, 2014. |
[2] |
管振清, 郭树祥, 单联瑜. 利用低频信息改善地震成像质量[J]. 断块油气田, 2009, 16(1): 39-41. GUAN Zhenqing, GUO Shuxiang, SHAN Lianyu. Improving the quality of seismic imaging using low-frequency information[J]. Fault-Block Oil & Gas Field, 2009, 16(1): 39-41. |
[3] |
韩站一, 尚新民, 孙成禹, 等. 火成岩下伏地层地震反射波场特征研究[J]. 石油物探, 2009, 48(3): 226-231. HAN Zhanyi, SHANG Xinmin, SUN Chengyu, et al. Study on characteristics of seismic reflection wave field under igneous rock[J]. Geophysical Prospecting for Petroleum, 2009, 48(3): 226-231. |
[4] |
王仁冲, 徐怀民, 邵雨, 等. 准噶尔盆地陆东地区石炭系火山岩储层特征[J]. 石油学报, 2008, 29(3): 350-355. WANG Renchong, XU Huaimin, SHAO Yu, et al. Reservoir characteristics of Carboniferous volcanic rocks in Ludong area of Junggar Basin[J]. Acta Pe-trolei Sinica, 2008, 29(3): 350-355. |
[5] |
赵宁, 石强. 裂缝孔隙型火山岩储层特征及物性主控因素——以准噶尔盆地陆东-五彩湾地区石炭系火山岩为例[J]. 天然气工业, 2012, 32(10): 14-23. ZHAO Ning, SHI Qiang. Characteristics of fractured pore-type volcanic reservoirs and the main controlling factors of physical properties: A case study of the Carboniferous volcanic rocks in Ludong-Wucaiwan area of Junggar Basin[J]. Natural Gas Industry, 2012, 32(10): 14-23. |
[6] |
罗静兰, 邵红梅, 张成立. 火山岩油气藏研究方法与勘探技术综述[J]. 石油学报, 2003, 24(1): 31-32. LUO Jinglan, SHAO Hongmei, ZHANG Chengli. Review of research methods and exploration technolo-gies for volcanic reservoirs[J]. Acta Petrolei Sinica, 2003, 24(1): 31-32. |
[7] |
王卫华. 速度横向变化对地震成像的影响[J]. 石油地球物理勘探, 2007, 42(6): 720-726. WANG Weihua. Impact of lateral velocity changes on seismic imaging[J]. Oil Geophysical Prospecting, 2007, 42(6): 720-726. |
[8] |
张军华. 地震资料去噪方法——原理、算法、编程及应用[M]. 山东东营: 中国石油大学出版社, 2011. ZHANG Junhua. Denoising Method of Seismic Data: Principle, Algorithm, Programming and Application[M]. Dongying, Shan-dong: China University of Petroleum Press, 2011. |
[9] |
苏贵仕, 丁成震. 叠前随机噪声衰减及其应用技巧[J]. 石油地球物理勘探, 2014, 49(增刊1): 87-92. SU Guishi, DING Chengzhen. Random noise attenua-tion on prestack seismic data and its application skills[J]. Oil Geophysical Prospecting, 2014, 49(S1): 87-92. |
[10] |
张白林, 潘树林, 尹成. 地震资料数字处理方法[M]. 北京: 石油工业出版社, 2011. ZHANG Bailin, PAN Shulin, YIN Cheng. Digital Pro-cessing Method of Seismic Data[M]. Beijing: Petroleum In-dustry Press, 2011. |
[11] |
徐海, 樊太亮, 赵夏, 等. 基于矢量中值滤波的ISS高密度地震数据去噪方法[J]. 石油地球物理勘探, 2013, 48(3): 340-344. XU Hai, FAN Tailiang, ZHAO Xia, et al. A deblending method based on vector median filter for ISS high density seismic data[J]. Oil Geophysical Pro-specting, 2013, 48(3): 340-344. |
[12] |
严又生. 一种非线性滤波去噪方法[J]. 石油地球物理勘探, 1998, 33(1): 104-108. YAN Yousheng. A non-linear filtering and denoising method[J]. Oil Geophysical Prospecting, 1998, 33(1): 104-108. |
[13] |
林娟, 罗勇, 张龙, 等. GeoEast处理系统在准南复杂高陡构造地区的应用[J]. 石油地球物理勘探, 2014, 49(增刊1): 60-65. LIN Juan, LUO Yong, ZHANG Long, et al. GeoEast seismic data processing in complex steep structure areas, Junggar Basin[J]. Oil Geophysical Prospecting, 2014, 49(S1): 60-65. |
[14] |
张军华, 张在金, 张彬彬, 等. 地震低频信号对关键处理环节的影响分析[J]. 石油地球物理勘探, 2016, 51(1): 54-62. ZHANG Junhua, ZHANG Zaijin, ZHANG Binbin, et al. Analysis of the influence of seismic low-frequency signals on key processing steps[J]. Oil Geophysical Prospecting, 2016, 51(1): 54-62. |
[15] |
王金龙, 胡治权. 三维锥形滤波方法研究及应用[J]. 石油地球物理勘探, 2012, 47(5): 705-711. WANG Jinlong, HU Zhiquan. Research and Applica-tion of 3D cone filtering method[J]. Oil Geophysical Pro-specting, 2012, 47(5): 705-711. |
[16] |
伍建, 王润秋, 魏加明, 等. 可控震源地震数据谐波滤方法[J]. 石油地球物理勘探, 2014, 49(1): 47-52. WU Jian, WANG Runqiu, WEI Jiaming, et al. Har-monic filtering method of controllable seismic data[J]. Oil Geophysical Prospecting, 2014, 49(1): 47-52. |
[17] |
宋玉琢, 王维红, 初海红. GeoEast保幅处理技术在大庆油田龙26区块的应用[J]. 石油地球物理勘探, 2014, 49(增刊1): 41-47. SONG Yuzhuo, WANG Weihong, CHU Haihong. Application of GeoEast amplitude-preserving processing technology in Long 26 block of Daqing Oilfield[J]. Oil Geophysical Prospecting, 2014, 49(S1): 41-47. |