地震数据提高分辨率处理监控评价技术

1、地震数据提高分辨率处理监控评价技术刁瑞期刊名称】岩性油气藏年(卷),期】2020(032)001【总页数】8页(P94-101) 【关键词】 地震数据; 高分辨率; 效果评价; 软件模块; 效果保真【作 者】 刁瑞【作者单位】 中国石化胜利油田分公司物探研究院 山东东营 257022【正文语种】 中 文【中图分类】 P631.40 引言薄互层岩性油气藏逐渐成为勘探开发的重点领域，薄层砂体、河道砂体等储层对分 辨率的要求越来越高 1。常规地震资料难以满足薄层砂体储层刻画精度的要求， 需要开展针对性的地震资料提高分辨率处理2。提高地震分辨率是地震资料处 理的重要环节之一，方法主要包括：分频带反褶积

2、、预测反褶积、零相位反褶积、反Q滤波、谱模拟反褶积、井约束提频和吸收补偿等方法3-8。各种提高地震分辨率处理的方法在实际地震资料处理中发挥了重要作用，并取得了良好的应用效 果。由于薄层砂体储层对提高地震分辨率处理方法提出了更高的要求，针对不同品 质的地震数据，应该采用何种提高分辨处理的方法和参数，以及地震数据的分辨率 应该提高到何种程度等问题，均制约着提高分辨率处理的效果。刘浩杰等9从 客观量化角度对地震资料分辨率与频谱特征参数的关系进行了研究。李曙光等 10从剖面和频谱特征对 3 种提高分辨率处理的方法进行了探讨。钱荣钧11 分析了分辨率的分类和定义，并给出了计算空间分辨率的方法。万欢等12

3、通 过频谱、信噪比和相似系数等参数，对比分析了 4 种提高分辨率处理技术，在叠 前地震资料处理中取得了较好的应用效果。在常规地震资料提高分辨率处理的过程中，由于缺乏系统的提高地震分辨率处理效 果的监控评价分析方法和软件模块，需要对提高分辨率处理的方法和参数等进行大 量的重复试验，并需要凭借处理人员的经验来判定方法及参数选取的正确性和可靠 性，人为因素影响较大。针对薄层砂体精细描述的地质需求，开展提高分辨率处理 效果的监控评价方法研究，并研发软件系统，优选提高分辨率处理方法、流程和参 数，全面可靠地监控提高分辨率处理的流程，实现提高分辨率处理效果的监控评价 分析，并在准噶尔盆地中部庄 3 井地区

4、进行推广应用，对提高分辨率处理进行全 流程的效果监控评价分析，以期提高分辨率处理后识别薄层砂体的能力，为后续的 油藏描述和井位部署奠定基础。方法原理如何实现提高地震分辨率处理方法、流程及参数的优化选择，一直是提高地震分辨 率处理所面临的难题。在常规的频谱和剖面对比分析的基础上，提出了多种属性计 算、组合优选的效果监控评价方法。首先，选用不同的提高地震分辨率处理的方法、流程及参数进行数据处理；然后， 通过数学公式计算目标层位地震数据的属性参数，包括：振幅、频率、相位、信噪 比、时频特征、相干特征、振幅比特征、波形特征、断层特征等；再次，根据计算 得到的属性数值进行综合分析，对提高分辨率处理的结果

5、进行效果量化评价；最后， 从不同测试结果中优选出最理想的效果，进而确定最佳的提高分辨率处理的方法、 流程及参数。通过技术(参数)差异性绘图，结合属性提取、合成记录标定等多因 素综合分析，可以快速高效地优选出最佳提高分辨率处理的结果，缩短了地震资料 提高分辨率处理的周期。频谱类方法 地震波在地层中传播，随着传播距离的增加，频率成分会发生变化，高频成分相对 于低频成分更容易衰减，包含不同频率成分的频谱则能够反映地震数据的分辨率 13。通过傅里叶变换算法可以求得地震数据的振幅谱和分贝谱，进而计算主 频、优势频宽、有效频宽、倍频程、低截频值、高截频值、不同频段能量值等属性， 并分析地震数据不同提高分辨

6、率处理的频谱属性，从而实现最优参数的优选。分贝 谱和振幅谱的计算公式分别为式中：d B(f)为分贝谱；|H(f)|为归一化后的振幅谱；H(3)为振幅谱，3=2 n f； h(t)为地震记录；f为频率，Hz；t为时间，ms。对地震数据不同提高分辨率处理结果进行傅里叶变换，在频谱中计算主频、频宽、 各频段能量值等属性，通过频谱类分析方法能够直观地定量评价提高分辨率处理的 效果。优势频宽以-18 dB为标准，有效频宽以-24 dB为标准。主频越高，有效频 带越宽，倍频程越大，则地震分辨率越高，识别薄层砂体的能力越强。随着“两宽高”地震技术的发展，越来越重视低频段信息的恢复和保护，特别是310 Hz

7、的低频段信息。低截频值越小，则表示低频段信息越丰富，倍频程也越大，反演结 果越可靠，因而地震分辨率越高。子波类方法 地震子波的提取方法主要包括：自相关方法、同态谱方法、谱模拟方法、高阶统计 量方法、高阶谱方法和非线性理论方法等。通过地震记录的高阶累积量可以近似得 到地震子波的高阶累积量，根据地震记录的功率谱可计算得到子波振幅谱，利用地震记录的双谱重构地震子波的相位。地震子波频谱的计算公式14为地震子波的相位谱与双谱之间的关系为式中：3为角频率，rad/s ;碰（3）为地震子波相位谱；屮3b（w1,w2）为根据地震 信号三阶累积量加时窗计算得到的双谱。在地震数据目标层位划取时窗提取地震子波，一方

8、面，可以计算地震子波的主频、 频宽、相位等，量化评价方法与频谱类方法的评价标准相一致，主频越高，有效频 带越宽，倍频程越大，则地震分辨率越高；另一方面，对提高分辨率处理前、后的 子波进行互相关和子波一致性分析，互相关数值越大，则越保幅、保真。图1 为 提高分辨率处理前、后的地震子波，提频前子波主频为32 Hz、频宽为48 Hz，提 频后子波主频为45 Hz、频宽为74 Hz，频带拓宽26 Hz。另外，在压缩地震子波、 提高分辨率处理的同时，地震资料子波保持了较高的相似性，提高分辨率处理具有 较高的保真性。图1提高分辨率处理前（蓝）、后（红）的地震子波Fig.1 Wavelet of seism

9、ic data before（blue）and after（red）improving high resolution processing信噪比类方法 信噪比是衡量地震资料品质的重要参数，以提高分辨率并保持一定的信噪比为标准， 通过分析提高分辨率处理前、后数据的信噪比，可以定量地评价处理效果，从而根 据信噪比分析结果优选最佳参数。信噪比分析方法包括：能量叠加法、频谱估计法、 功率谱估算法、相关法和特征值法等15。通过信噪比分析可以避免过多高频噪音的不利影响，从而确保扩展频带宽度的可信 度。由于不同层位或同一层位的地震波能量差异巨大，通常情况下在某一时窗内进 行信噪比分析。能量叠加法的信噪比计

10、算公式为式中：M为分析时窗的采样点数，个；N为地震道数；dij为地震记录振幅；Es 为有效信号能量；E为地震记录能量。频谱估算法的信噪比计算公式为式中：fL，fH，fc分别为信号有效频带的起始频率、截止频率、振幅谱截止频率， Hz ； X(f)为地震记录振幅谱；Xs(f)为地震信号振幅谱。通过信噪比分析对提高分辨率处理的结果进行效果评价，在多种方法达到相同的提 高分辨率处理效果时，通过信噪比的高低进行评价，优选信噪比最高的提频效果， 从而指导方法、流程或参数的优选。时频分析方法 时频分析方法解决了一维傅里叶变换中时间域和频率域分离的不足，可以同时分析 不同时刻包含的各频率成分以及不同频率成分随

11、传播时间的变化特征16-17。修正S变换的表达式为式中：GSTN(T,f)为二维时频谱；f为频率,Hz；t为时间，ms;入p均为修正S 变换调节参数，可以获得不同分辨率的时频分析结果。对比提高分辨率处理前、后的时频谱特征，一方面避免了傅里叶变换时窗大小的限 制，能够更加精细地分析频率变化特征；另一方面可以分析不同频率成分随时间的 变化情况，即不同频段能量的时间变化曲线，其计算公式为式中：t为时间,ms；f1和f2分别为频段的起始频率和截止频率，Hz ； GSTN(t,f)为地震记录的二维时频谱。进行提高分辨率处理后，高频段能量曲线的数值明显增大，表明时频谱中优势频带 和有效频带宽度均增加，而高

12、频段能量曲线的变化趋势保持不变，表明地震记录能 量得到了很好地恢复和补偿，保持了原始数据的频谱特征和时频特性。图2 为提高分辨率处理前、后二维时频谱对比。从图2 可以看出，提高分辨率处 理前、后频带分别为10 45 Hz和10 70 Hz,地震记录深层能量也得到了很好 的补偿，时频谱的能量匹配较好，保持了原始数据的频谱特征和时频特性。图2提高分辨率处理前（a）、后（b）时频谱对比Fig.2 Comparison of timefrequency spectrum of seismic data before（a）and after（b）high resolution processing实际资

13、料应用效果常用的编程语言主要有C/C+，Fortran，Matlab等17-18 ,此次采用C+ + 和Qt编程语言进行提高地震分辨率处理的效果评价软件研发，包括：频率分析、 波形分析、子波一致性分析、能量分析、信噪比分析、振幅谱分析、时频谱分析、 子波分析、沿层属性分析等模块功能。对原始地震记录进行提高分辨率处理，得到不同方法和参数的一系列结果，可以对 单个数据进行频谱分析、时频分析和子波分析等，也可以对一系列提高分辨率处理 的结果进行编号，计算属性特征的变化趋势，从而根据柱状图进行定性和定量评价 对单个数据进行评价分析，在地震剖面上划取时窗可以快速实现频谱分析、时频分 析和子波分析等功能。

14、图3 为地震剖面和二维时频谱分析结果，通过二维时频谱 能够直观地分析各频率成分随时间的变化情况。将一系列提高分辨率处理的数据进行编号，分析不同方法或参数的属性变化特征， 通过计算的属性值进行定量评价，从而优选出最佳提频方法或参数。以频谱类量化 评价方法为例，如图4 所示，频谱类属性包括：主频、频宽、低截频值、高截频 值、低频段能量和高频段能量、倍频程等。对不同提高分辨率处理的方法和参数的 测试结果进行编号，图4 中共计6 个处理参数。从主频属性分析结果显示，不同 处理参数的主频基本相同，编号46的主频略高，约为20 Hz ；从频宽属性、高 截频值属性、高频段能量属性的分析结果可看出，编号为 4

15、 和 5 的处理参数最佳， 频宽约为44 Hz，高截频值为50 Hz ；低截频值属性显示编号14处理效果较佳， 约为6 Hz ，但低频段能量属性显示编号15 的效果整体比较理想。 综合考虑主频、频宽、低截频值、高截频值、低频段能量和高频段能量，最佳提频 参数是编号4，主频为20 Hz，频宽为44 Hz，低截频值为6 Hz，高截频值为50 Hz ，倍频程为3。通过提高分辨率处理的效果监控评价方法，可以对提高地震分 辨率处理数据进行定量评价，从而优选最佳方法、流程或参数。图 3 地震剖面（a ）及时频谱（b）Fig.3 Seismic section（ a）and timefrequency sp

16、ectrum（b）图 4 频谱类方法监控评价结果 Fig.4 Monitoring and evaluation results of spectrum准噶尔盆地中部庄3井地区主体构造位于盆1井西凹陷，其东南是莫索湾凸起， 构造的主体形态为西北向东南渐深的单斜，主要目的层为三叠系、侏罗系、白垩系 等，目的层埋深为3 000 -6 000 m。地震资料存在纵横向分辨率低、反射轴杂乱、 层位追踪困难、空白反射较多、断点与断面不清楚等问题，造成低序级断层识别难 度大，储层精细描述落实存在一定风险，迫切需要高分辨率处理地震资料的支撑。 将提高地震分辨率处理效果的监控评价方法在庄3 井地区进行了应用，实

17、现了提 高分辨率处理全流程的监控，确保优选最佳方法和参数，从而获取最优的提高分辨 率处理数据。采用逐级提频的方式提高地震资料分辨率处理，以信噪比、子波一致性和沿层属性 为监控手段，最大程度地提升主频、拓展频宽。图5 为提高分辨率处理前、后的 剖面效果对比，通过分步逐级提频处理，避免信噪比的大幅降低。图6 为原始数 据、叠前提频和叠后提频的频宽属性效果对比，原始数据主频约为35 Hz，叠前提 频后主频约为40 Hz，叠后提频后主频约为50 Hz，主频值不断提升。从图5处 理前、后的地震剖面可以看出，提高分辨率处理成果数据解决了原始资料分辨率低 反射轴杂乱、层位难以追踪的问题，提高分辨率处理资料的

18、断裂特征更加明确，三 工河组反射特征更加清晰。从图6 中可以看出，提高分辨率分辨率处理后，频带 拓宽15 25 Hz,主频提升15 Hz以上。图 5 地震剖面效果对比 Fig.5 Comparison of seismic section图 6 频宽属性切片对比 Fig.6 Comparison of frequency band width slices 从合成记录标定结果（图7）可见，通过Z4井、Z110井、Z102井、Z109井的 合成记录标定，提高分辨率处理结果与井资料更加吻合，验证了提高分辨率处理结 果的可靠性，有助于下一步的井位部署。图 7 合成记录标定结果 Fig.7 Calib

19、ration results of synthetic records结论（1）针对提高地震分辨率处理效果监控的迫切需求，创新研究了提高地震分辨率 处理效果的监控评价技术，从9个方面的属性特征对不同处理方法、流程或参数 进行了量化评价，实现了提高分辨率处理效果评价由定性向量化的转变。（2）利用效果监控评价技术对提高分辨率处理进行全流程的监控，在每一步提高 分辨率处理中均进行了效果监控评价，最终建立了针对不同地质目标的最佳处理流 程，能够大幅提高薄层砂体的识别能力。（3）随着高密度地震技术的不断发展，地震数据量呈指数型增加，提高分辨率处 理效果监控评价技术既能够大幅提高处理效率，缩短处理周期，还

20、可以应用于海量 数据的质量监控和其他处理环节的效果评价中，具有广阔的应用前景。参考文献：相关文献】李庆忠岩性油气藏地震勘探若干问题的讨论(口)：关于垂向分辨率的认识岩性油气藏， 2008，20(3): 1-5.LI Q 乙Discussion on seismic exploration of lithologic reservoirs (H): Seismic vertical resolution.Lithologic Reservoirs，2008，20(3)：1-5. 2 姚健，郑敬贵，杜晓峰，等.辽东湾探区古近系浊积体针对性地震处理与识别技术.岩性油气藏， 2012，24(1)：69

21、-73.YAO J ， ZHENG J G，DU X F，et al.Seismic processing and identification technology for Paleogene turbidite in Liaodong Bay.Lithologic Reservoirs， 2012，24(1)：69-73.3 徐兴荣，苏勤，王劲松，等加权MPFI方法及其在三维连片处理中的应用岩性油气藏， 2019，31(1)：122-129.WANG X R，SU Q，WANG J S，et al.Weighted MPFI method and its application in 3D

22、 joint processing.Lithologic Reservoirs，2019，31(1)：122-129.4 刁瑞分频带预测反褶积方法研究断块油气田，2015 ,22 ( 1 ): 53-57.DIAO R.Method of separate frequency predictive deconvolution.Fault-Block Oil&Gas Field，2015，22(1)： 53-57. 5 王万里，李国发，桂金咏.混合相位子波有色反褶积.岩性油气藏， 2013 ， 25 ( 3 ) : 82- 86.WANG W L，LI G F，GUI J Y.Colored d

23、econvolution of mixed phase wavelet.Lithologic Reservoirs，2013，25(3):82-86.王孝，刘文卿，曾华会，等.复杂区分层约束近地表建模方法及应用.岩性油气藏，2018，30:68-73.WANG X，LIU W Q，ZENG H H，et al.Stratified constrained near-surface model building method and its application in complex surface area.Lithologic Reservoirs ， 2018，30(4):68-73.杨

24、应，杨巍，朱仕军基于EEMD的高分辨率层序t也层划分方法岩性油气藏，2018 , 30:59-67.YANG Y，YANG W，ZHU S J.Method for high-resolution sequence stratigraphy division based on Ensemble Empirical Mode Decomposition.Lithologic Reservoirs，2018，30(5):59-67.刁瑞，尚新民，芮拥军，等.时频域谱模拟反褶积方法研究.岩性油气藏，2013，25(1): 116-121.DIAO R，SHANG X M，RUI Y J，et al.

25、Study on spectral modeling deconvolution in time-frequency domain.Lithologic Reservoirs，2013，25(1):116-121. 9 刘浩杰，王延光.也震资料分辨率若干问题的探讨.油气也球物理， 2009 ， 7 ( 4 ) : 5-10.LIUH J，WANG Y G.Discussion of some question on resolution for seismic data.Petroleum Geophysics，2009，7(4):5-10. 10 李曙光，唐建明，徐天吉，等.几种提高也震资料

26、分辨率的方法及效果分析.勘探也球物理进 展，2010 ,33 ( 5 ): 323-327.LI S G , TANG J M , XU T J , et al.Methods for improving seismic data resolution.Progress in Exploration Geophysics，2010，33(5):323-327.钱荣钧.也震波分辨率的分类研究及偏移对分辨率的影响.石油也球物理勘探,2010,45:306-313.QIAN R J.Classification studies on seismic resolution and migrations

27、 affect to resolution.Oil Geophysical Prospecting，2010，45（2）：306-313.万欢，樊小意，刘涛，等.叠前地震资料提高分辨率处理方法及应用.地球物理学进展，2012， 27（1）：304-311.WAN H，FAN X Y，LIU T，et al.Methods and applications for improving pre-stack seismic data resolution.Progress in Geophysics ，2012，27（1）： 304-311.王永生，胡杰，张静，等.宽频勘探技术在柴北缘深层侏罗系勘探

28、中的应用.岩性油气藏， 2017，29（6）：101-107.WANG Y S，HU J，ZHANG J，et al.Wide-frequency prospecting technology and its application on deep-seated Jurassic exploration in northern margin of Qaidam Basin.Lithologic Reservoirs，2017，29（6）：101-107.崔庆辉，芮拥军，尚新民，等.混合相位地震子波提取及应用.石油物探，2011，50（5）： 481-486.CUI Q H，RUI Y J，SHANG X M，et al.Mixed-phase seismic wavelet extraction and its application.Geophysical Prospecting for Petroleum，2011，50（5）：481-486.张军华，臧胜涛，周振晓，等.地震资料信噪比量化计算及方法比较.石油地球物理勘探， 2009，44（4）：481-486.ZHANG J H，ZANG S T，ZHOU Z X，et al.Quantitative compu