石油地球物理勘探  2022, Vol. 57 Issue (s2): 29-34  DOI: 10.13810/j.cnki.issn.1000-7210.2022.S2.005
0
文章快速检索     高级检索

引用本文 

张春燕, 李亚旭, 王晓琳, 侯山红. GeoEast融合数据子波一致性处理技术及其应用研究. 石油地球物理勘探, 2022, 57(s2): 29-34. DOI: 10.13810/j.cnki.issn.1000-7210.2022.S2.005.
ZHANG Chunyan, LI Yaxu, WANG Xiaolin, HOU Shanhong. Research on wavelet consistency processing methods of GeoEast and its application in emerging data. Oil Geophysical Prospecting, 2022, 57(s2): 29-34. DOI: 10.13810/j.cnki.issn.1000-7210.2022.S2.005.

本项研究受国家自然科学基金项目“面向海洋深水资料的全波场最小二乘偏移方法研究”(41874164)资助

作者简介

张春燕  高级地球物理师,1986年生。2008年毕业于中国石油大学(华东),获勘查技术与工程专业学士学位;2011年毕业于中国石油大学(华东),获地球探测与信息技术专业硕士学位。目前就职于中国石油杭州地质研究院计算机应用研究所,主要从事地震方法研究及资料处理

张春燕, 浙江省杭州市西湖区西溪路920号中国石油杭州地质研究院, 310023。Email: zhangchy-jl@petrochina.com.cn

文章历史

本文于2022年6月15日收到,最终修改稿于同年9月14日收到
GeoEast融合数据子波一致性处理技术及其应用研究
张春燕 , 李亚旭 , 王晓琳 , 侯山红     
①. 中国石油杭州地质研究院,浙江杭州 310023;
②. 中油测井地质研究院,陕西西安 710000;
③. 中国石油吉林油田公司地球物理勘探研究院,吉林松原 138000
摘要:地震资料处理经常需要进行不同采集参数数据间的融合处理,其中子波一致性处理对改善成像质量、提高空间分辨率至关重要。利用静校正技术消除数据时差差异后,利用GeoEast系统的特色处理模块,可以建立针对不同采集条件下数据的子波一致性处理流程。首先,在保真去噪的基础上,利用振幅处理模块完成能量一致性处理,特别是利用分频振幅补偿技术改善近地表造成的高频能量衰减问题;然后,通过地表一致性反褶积处理消除近地表子波鸣震反射;最后,基于重叠区数据的特征差异设计整形因子,进一步消除激发子波差异。实际资料处理结果证明,经过融合一致性处理后,资料整体信噪比和分辨率都有所提高,不同震源子波类型趋于一致,空间子波的频率和相位特征均一致,利用标准反射层的沿层属性验证了子波的频率、相位的一致性。
关键词子波一致性    融合处理    子波差异    相位特征    
Research on wavelet consistency processing methods of GeoEast and its application in emerging data
ZHANG Chunyan , LI Yaxu , WANG Xiaolin , HOU Shanhong     
①. PetroChina Hangzhou Institute of Geology, Hangzhou, Zhejiang 310023, China;
②. Research Institute of Geology of CNPC Logging, Xi'an, Shaanxi 710000, China;
③. Geophysical Prospecting Research Institute, Jilin Oilfield Company, PetroChina, Songyuan, Jilin 138000, China
Abstract: In seismic data processing, it often needs to process emerging data of different acquisition parameters, among which wavelet consistency processing plays a vital role in improving imaging quality and spatial resolution. After using static correction technology to eliminate the time difference of data, we employ the characteristic processing modules of the GeoEast system to establish the technical flow of wavelet consistency processing for data acquired under different conditions. On the basis of fidelity denoising, firstly, the amplitude processing module is used to complete the energy consistency processing, especially the use of the frequency-division amplitude compensation technology for improvement in the high-frequency energy absorption problem caused by near surface. Then, through surface-consistent deconvolution, the near-surface reverberation and reflection of wavelets are eliminated. Finally, a shaping factor is designed on the basis of the characteristic difference in the overlapping data of different areas to further eliminate the difference in excitation wavelets. The actual data processing results show that after the emerging consistency processing, the overall signal-to-noise ratio and resolution of the data are improved. The source wavelets of different types tend to be consistent, and the frequency and phase characteristics of spatial wavelets are consistent. The frequency and phase consistency of wavelets is verified by the along-layer attribute of the standard reflection layer.
Keywords: wavelet consistency    emerging processing    wavelet difference    phase feature    
0 引言

目前在地震资料采集中,不同采集条件下的数据融合处理越来越常见。不同的激发、接收条件使得数据的子波特征(能量、频率、相位等)空间差异巨大,这些差异会给后续偏移成像带来困扰,造成画弧、虚假构造等现象,降低资料的空间分辨率,甚至可能导致错误的解释结果[1]。因此,解决不同采集条件下数据间激发子波的能量、频率、波形及相位等差异问题,使激发子波具有较好的空间一致性,是消除采集脚印、提高地震资料信噪比和分辨率的关键[2]

1 工区概况

本文针对准噶尔盆地西北缘的两块三维工区数据进行连片融合处理,这两片工区数据分别采集于2015年和2016年。工区炮点位置及激发井深见图 1,可以看到炮点位置分布较规范,且两个工区炮点位置有重叠,这对于子波一致性处理是有利的。

图 1 炮点井深分布 蓝色井深为0,对应可控震源激发;其他对应井炮激发;①~⑤为图 4中选取的分析炮位置。

连片工区位于环玛湖构造带,地表条件复杂。西部为滩涂区,中部为盐碱地,东部为公益林区。从表层结构看,低降速带厚度为3~24m,厚度变化较快,因此对地震数据的吸收作用在空间上同样存在较大差异。复杂的近地表条件造成地表局部变观较多,炮集数据间的噪声类型、能量、主频等差异较大。

2 GeoEast子波一致性处理技术

融合处理过程中,消除不同震源类型子波的时间差异是必不可少的基础工作,其处理效果决定了最终成像质量。本文利用GeoEast特色模块,建立了一套融合数据子波一致性处理技术流程(图 2)。静校正作为必须的常规处理步骤,主要作用是消除不同区块数据间的时差差异,避免将其视为数据间的子波相位差异。必须指出的是,基准面静校正必须是融合计算,剩余静校正可分区块进行,也可以融合处理,依实际情况而定。

图 2 GeoEast系统子波一致性处理技术流程

不同激发震源、检波器具有不同的激发接收响应,导致地震子波存在相位、频率和能量方面的差异。近地表的吸收因地表条件的不同而出现较大差异,进而造成地震子波传播过程中的振幅(能量)、频率和相位等再一次被改造,特征差异进一步加大(图 3),特别是高频能量衰减及不同地表条件带来的近地表鸣震是子波一致性处理的重中之重[3]图 3中以单炮统计自相关的共炮点叠加近似为地震子波,可看出可控震源(左)和井炮数据(右)能量级别不同,即便在井炮数据区域(右侧),单炮间统计子波差异也是非常明显的。

图 3 处理前质控线炮集统计子波特征差异分析图

图 1中井炮数据区选择5炮(图中四角星位置),分别位于不同地表条件,对单炮数据进行频谱统计分析(图 4)。可以看出,单炮数据间能量差异明显,能量差为6~16倍。从图中还可以看出,资料的优势频带为4~40Hz,炮间频率差异非常明显,单炮的优势主频、频宽差异较大。

图 4 选取的图 1中不同位置单炮统计频谱分析

首先,利用时频域分频补偿的处理思想,补偿不同地表条件对高频能量的吸收,恢复激发子波的高频成分,消除不同数据间激发子波的频率差异,并结合地表一致性振幅补偿处理,消除激发能量的空间差异,使不同采集条件下的数据激发能量趋于一致[4-5];其次,利用地表一致性反褶积消除近地表子波鸣震响应,使不同类型激发子波的旁瓣能量得到衰减或消除,并对子波的主瓣进行有效压缩,提高资料的分辨率;最后,应用整形处理消除残余的子波差异。

3 应用效果

本次融合处理的数据的采集条件差异并不大,且施工年代接近,检波器类型也无变化。激发类型包括井炮和可控震源,因此首先要利用小相位化消除震源类型差异带来的相位差异。虽然通常认为可控震源激发的地震子波是零相位的,实际上并非如此。受近地表吸收作用的影响,大部分可控震源激发的地震子波都是混合相位的,因此最小相位化不能完全消除两种震源子波的相位差异。凌云等[6]认为在地表强吸收情况下,必须利用可控震源反褶积才能将可控震源激发子波转化为最小相位子波。只有在吸收作用不严重的情况下才可将可控震源激发的子波视作零相位子波。分析图 5所示本区单炮形态,可知近地表吸收作用适中,可将可控震源子波近似为零相位子波。

图 5 工区可控震源单炮形态示意图

张伟等[7]通过模型数据正演研究,深入分析了能量一致性对偏移成像的影响,发现能量差异不仅会造成剖面画弧,还会使剖面的分辨率发生横向变化。连片处理中消除区块内和区块间的能量、频率差异对最终成果的品质起到了决定性作用[8-13]。利用时频域衰减补偿和地表一致性振幅补偿可消除近地表低降速带及地震波传播过程中带来的能量衰减,恢复地震波场的高频能量。图 6为原始数据和经振幅一致性处理的激发能量空间分布图,可以看出,不同震源激发带来的振幅差异得到了较好消除。

图 6 原始数据(上)和处理后数据(下) 单炮记录统计能量空间分布图

消除能量差异后,还需要利用地表一致性反褶积消除近地表带来的子波鸣震差异。从图 7所示质控线子波一致性处理后炮集统计子波分析可见,应用子波一致性处理后,炮间激发子波的旁瓣特征趋于一致,子波一致性进一步得到增强。

图 7 质控线子波一致性处理后炮集统计子波分析 图中两条曲线分别为炮点高程和井深,可控震源井深为0,可用于区分震源类型

将两种震源数据直接拼接在一起,得到图 8所示剖面。由图可见,拼接处两种震源数据依然存在明显差异,即不同类型炮集数据存在残余相位波形差异[14-18]。将两种震源数据进行融合,得到图 9所示剖面,可见剖面上没有明显的数据差异,这是因为高覆盖次数的可控震源数据占主体,叠加时几乎完全压制了井炮数据。但是,剖面上仍能看到视主频存在空间变化。应用GeoEast软件的DataClip模块可以快速寻找到两种数据的重叠区域。在本次处理中,寻找到两种震源数据的重叠区后,利用重叠区叠加剖面,以井炮激发数据为期望输出提取整形因子,对可控震源数据进行整形处理,消除残余的子波差异。整形后的拼接剖面见图 10,可见可控震源和井炮数据的统计子波达到一致,资料信噪比和分辨率得到明显提高。

图 8 质控线不同震源类型数据拼接剖面 箭头左侧为可控震源数据,右侧为井炮数据;箭头所指为数据交接重叠区域。图 9图 10同。

图 9 质控线融合数据叠加剖面

图 10 子波整形后质控线不同震源类型数据拼接剖面

基于子波一致性处理后的数据得到偏移剖面,再对成果资料作进一步质控分析。基于井信息特征和区域统层认识,拾取区域内的标志性反射层,作为地震属性的质控层位。因为储层特征一般变化较大,所以储层不能选作质控层位。沿标志性反射层提取属性切片,得到工区沿层瞬时频率和瞬时相位图(图 11)。由图可见,从平面属性响应上看不到采集痕迹,说明本次子波一致性处理中井炮和可控震源激发子波的一致性较好,可以认为区域子波一致,数据具备最优的空间分辨率。

图 11 工区标志性反射层位属性切片 左:沿层瞬时频率;右:沿层瞬时相位
4 结论

在不同区块地震数据的融合处理中,最大限度改善空间子波一致性对提高地震成果空间分辨率至关重要。融合处理中不同震源类型子波差异较大,应在消除时差的基础上针对波形、振幅、相位、频率等参数进行一致性处理,消除近地表激发、接收、低降速带吸收衰减等因素对子波能量、频率、相位的影响,是子波一致性处理的技术要点。

首先,应用分频振幅补偿技术消除激发震源类型及激发、接收条件不同所带来的空间子波的振幅差异,通过处理质控分析可以看出该技术方法切实有效。

然后,针对近地表造成的子波鸣震响应,利用地表一致性反褶积,以一致性作为评价标准,最大限度地衰减鸣震虚反射,并进行合理压缩,提高资料分辨率。处理过程中,不能完全依赖资料处理评价手段,要利用标准反射层的相关属性特征验证子波的空间变化,以监控采集脚印是否被消除,空间子波一致性是否满足要求。

本文中涉及可控震源和井炮的融合处理,通常做法是将可控震源激发子波作为最小相位进行处理,将最小相位化作为主要技术手段。随着可控震源的推广和应用,针对混合相位的子波处理技术将会得到发展,这对提高处理成果质量起到关键作用。

参考文献
[1]
羊屋三维处理、解释一体化方法研究组. 三维地震数据的分析和监测方法研究[J]. 石油地球物理勘探, 2002, 37(5): 433-440.
Yangwu Research Group. Yangwu 3-D processing and interpretation integration[J]. Oil Geophysical Prospecting, 2002, 37(5): 433-440. DOI:10.3321/j.issn:1000-7210.2002.05.001
[2]
俞寿朋. 高分辨率地震勘探[M]. 北京: 石油工业出版社, 1993.
[3]
凌云. 大地吸收衰减分析[J]. 石油地球物理勘探, 2001, 36(1): 1-8.
LING Yun. Analysis of attenuation by earth absorption[J]. Oil Geophysical Prospecting, 2001, 36(1): 1-8. DOI:10.3321/j.issn:1000-7210.2001.01.001
[4]
凌云, 高军, 吴琳. 时频空间域球面发散与吸收补偿[J]. 石油地球物理勘探, 2005, 40(2): 176-182, 189.
LING Yun, GAO Jun, WU Lin. Compensation for spherical dispersion and absorption in time-frequency-space domain[J]. Oil Geophysical Prospecting, 2005, 40(2): 176-182, 189. DOI:10.3321/j.issn:1000-7210.2005.02.018
[5]
凌云研究组. 叠前相对保持振幅、频率、相位和波形的地震数据处理与评价研究[J]. 石油地球物理勘探, 2004, 39(5): 543-552.
LING Yun Yanjiuzu. Study of seismic data processing and appreciation based on prestack relative preservation of amplitude, frequency, phase and waveform[J]. Oil Geophysical Prospecting, 2004, 39(5): 543-552. DOI:10.3321/j.issn:1000-7210.2004.05.010
[6]
凌云, 周熙襄. 可控震源自适应地表一致性反褶积[J]. 石油地球物理勘探, 1994, 29(3): 306-315, 356.
LING Yun, ZHOU Xixiang. Surface-consistent adaptive deconvolution for vibroseis data[J]. Oil Geophy-sical Prospecting, 1994, 29(3): 306-315, 356.
[7]
张伟, 于建华, 李云明, 等. 能量一致性对叠前时间偏移成像影响分析[J]. 地球物理学进展, 2012, 27(4): 1665-1670.
ZHANG Wei, YU Jianhua, LI Yunming, et al. The analysis of energy consistency for pre-stack time migration imaging[J]. Progress in Geophysics, 2012, 27(4): 1665-1670.
[8]
赵桂玲, 陈海峰, 高现俊, 等. 叠前偏移能量一致性处理方法及应用[J]. 地球物理学进展, 2016, 31(3): 1295-1299.
ZHAO Guiling, CHEN Haifeng, GAO Xianjun, et al. Prestack migration energy consistency processing method and application[J]. Progress in Geophysics, 2016, 31(3): 1295-1299.
[9]
裴江云, 张丽艳, 王丽娜, 等. 松辽盆地深层地震资料叠前时间偏移连片处理技术研究[J]. 地球物理学报, 2011, 54(2): 294-303.
PEI Jiangyun, ZHANG Liyan, WANG Lina, et al. Multi-survey joint pre-stack time migration proce-ssing study of deep beds in Songliao basin[J]. Chinese Journal of Geophysics, 2011, 54(2): 294-303. DOI:10.3969/j.issn.0001-5733.2011.02.004
[10]
左海, 魏庚雨, 何小松, 等. 三维地震资料叠前连片处理技术[J]. 石油地球物理勘探, 2008, 43(增刊1): 29-35.
ZUO Hai, WEI Gengyu, HE Xiaosong, et al. Block-jointed together processing technique of 3-D prestack seismic data[J]. Oil Geophysical Prospecting, 2008, 43(S1): 29-35.
[11]
曾华会, 王小卫, 苏勤, 等. 混合震源高精度匹配处理技术在中国西部复杂障碍区的应用[J]. 石油地球物理勘探, 2021, 56(3): 476-484.
ZENG Huahui, WANG Xiaowei, SU Qin, et al. Research and application of mixed sources high-precision matching processing technology in the complex obstacle area in western China[J]. Oil Geophysical Prospecting, 2021, 56(3): 476-484.
[12]
蒋立, 陈勇, 肖艳玲, 等. 地表过渡带近地表Q补偿与地表一致性反褶积处理效果对比研究[J]. 石油物探, 2018, 57(6): 870-877.
JIANG Li, CHEN Yong, XIAO Yanling, et al. A comparison of near-surface Q compensation and surface-consistent deconvolution in the near-surface transition zone[J]. Geophysical Prospecting for Petroleum, 2018, 57(6): 870-877.
[13]
干大勇, 高宪伟, 朱丹, 等. 提高浅层地震成像品质的处理技术——以川中地区沙溪庙组为例[J]. 石油地球物理勘探, 2020, 55(增刊1): 56-63.
GAN Dayong, GAO Xianwei, ZHU Dan, et al. Processing workflow for improving shallow seismic image quality on GeoEast platform: a case study on Shaxi-miao Formation in central Sichuan Basin[J]. Oil Geophysical Prospecting, 2020, 55(S1): 56-63.
[14]
和海雷, 史庆阳, 伊鸿斌, 等. 黄土山地可控震源与井炮联合激发配套技术研究[C]. 中国石油学会2021年物探技术研讨会论文集, 2021, 188-191.
[15]
吴德明, 戴海涛, 李国强, 等. 黄土塬井震混采资料子波一致性处理[C]. 中国石油学会2021年物探技术研讨会论文集, 2021, 568-571.
[16]
陈杨, 任勇, 丁在宇, 等. 复杂过渡带地震资料成像[C]. SPG/SEG南京2020年国际地球物理会议论文集(中文), 2020, 506-509.
[17]
刘传奇, 周建科, 李宾, 等. 利用Ricker子波相移特性估计地震资料剩余相位[J]. 石油地球物理勘探, 2021, 56(2): 249-256.
LIU Chuanqi, ZHOU Jianke, LI Bin, et al. Estimation of residual phase of seismic data by using the phase shift characteristics of Ricker wavelet[J]. Oil Geophysical Prospecting, 2021, 56(2): 249-256.
[18]
李晓峰, 王晓涛, 林娟, 等. 一致性处理技术在准噶尔盆地三维连片处理中的应用[J]. 石油地球物理勘探, 2018, 53(增刊2): 33-41.
LI Xiaofeng, WANG Xiaotao, LIN Juan, et al. Consistent processing for 3D seismic data merging in Junggar Basin[J]. Oil Geophysical Prospecting, 2018, 53(S2): 33-41.