地震波能量补偿的并行反Q滤波方法研究

doi:10.11885/j.issn.1674-5086.2021.02.03.03

西南石油大学学报(自然科学版) ›› 2023, Vol. 45 ›› Issue (1): 24-32.DOI: 10.11885/j.issn.1674-5086.2021.02.03.03

地震波能量补偿的并行反Q滤波方法研究

张全^1,2, 王一品¹, 张伟¹, 彭博¹, 胥林³

1. 西南石油大学计算机科学学院, 四川成都 610500;
2. 电子科技大学信息与通信工程学院, 四川成都 611731;
3. 西南石油大学信息学院, 四川南充 637001

收稿日期:2021-02-03 发布日期:2023-02-24
通讯作者: 张全,E-mail:zhangquan@swpu.edu.cn
作者简介:张全,1985年生,男,汉族,甘肃武威人,讲师,博士,主要从事地震信号处理、高性能计算及深度学习等方面的研究。E-mail:zhangquan@swpu.edu.cn
王一品,1996年生,男,汉族,河南漯河人,硕士,主要从事图像处理及计算机图形学等方面的研究。E-mail:wangyipin@sobey.com
张伟,1994年生,男,汉族,甘肃武威人,硕士,主要从事地震信号处理、并行计算及目标检测等方面的研究。E-mail:zhangwei_swpu2019@163.com
彭博,1980年生,男,汉族,四川南充人,教授,博士,主要从事医学图像与信号处理、医学超声成像及基于GPU并行计算等方面的研究。E-mail:bopeng@swpu.edu.cn
胥林,1977年生,男,汉族,四川资中人,教授,硕士,主要从事油田及企事业信息化建设、人工智能及数据挖掘与分析等方面的研究。E-mail:xulin@swpu.edu.cn
基金资助:
油气藏地质及开发工程国家重点实验室开放基金（PLN2022-51，PLN2021-21，PLN2021-25）

Research on Parallel Inverse Q Filtering Methods for Seismic Wave Energy Compensation

ZHANG Quan^1,2, WANG Yipin¹, ZHANG Wei¹, PENG Bo¹, XU Lin³

1. School of Computer Science, Southwest Petroleum University, Chengdu, Sichuan 610500, China;
2. School of Information and Communication Engineering, University of Electronic Science and Technology of China, Chengdu, Sichuan 611731, China;
3. School of Information, Southwest Petroleum University, Nanchong, Sichuan 637001, China

Received:2021-02-03 Published:2023-02-24

摘要/Abstract

摘要： 在石油勘探地震资料处理中，反Q滤波方法能有效地对地震波进行振幅补偿和相位校正，为地震反演和储层预测提供更准确的信息。对于大规模的地震道集数据处理，反Q滤波方法在CPU计算平台上执行时间较长，影响了地震解释的效率。分析发现，反Q滤波方法大量时间消耗在振幅相位补偿与短时傅里叶变换。在GPU平台上，首先，对振幅相位补偿部分进行并行化；其次，对批量短时傅里叶变换用CUFFT库进行加速；最后，对批量短时傅里叶变换进一步优化并将其应用于反Q滤波方法。实验结果表明，相比CPU计算环境，基于CUFFT库的反Q滤波并行算法效率提升了3.9倍，优化后的批量短时傅里叶变换进一步将效率提升了12%。

关键词: 反Q滤波, 振幅补偿, 傅里叶变换, 并行计算, 计算统一设备体系结构

Abstract: In seismic data processing of petroleum exploration, the inverse Q filtering method can effectively perform amplitude compensation and phase correction on seismic waves to provide more accurate information for seismic inversion and reservoir prediction. In large-scale seismic data processing, the inverse Q filtering method takes longer operation time under the CPU computing platform, which affects the efficiency of seismic interpretation. After analysis, it is found that the inverse Q filtering method consumes a lot of time in the short-time Fourier transform and calculates the amplitude and dispersion compensation terms. On the GPU platform, we first parallelizes the amplitude and dispersion compensation calculations, and accelerates the batch short-time Fourier transform with the CUFFT library, and then further optimizes the batch short-time Fourier transform and applies it to the inverse Q filtering method. The results show that compared with the CPU computing environment, the efficiency of the inverse Q filtering parallel algorithm based on the CUFFT library is improved by 3.9 times, and the optimized batch short-time Fourier transform further improves the efficiency of the parallel inverse Q filtering method by 12%.

Key words: Inverse Q filtering, amplitude compensation, Fourier transform, parallel computing, CUDA

中图分类号:

TE132
P631

张全, 王一品, 张伟, 彭博, 胥林. 地震波能量补偿的并行反Q滤波方法研究[J]. 西南石油大学学报(自然科学版), 2023, 45(1): 24-32.

ZHANG Quan, WANG Yipin, ZHANG Wei, PENG Bo, XU Lin. Research on Parallel Inverse Q Filtering Methods for Seismic Wave Energy Compensation[J]. Journal of Southwest Petroleum University(Science & Technology Edition), 2023, 45(1): 24-32.

0
/ / 推荐

导出引用管理器 EndNote|Ris|BibTeX

链接本文: http://journal15.magtechjournal.com/Jwk_xnzk/CN/10.11885/j.issn.1674-5086.2021.02.03.03

http://journal15.magtechjournal.com/Jwk_xnzk/CN/Y2023/V45/I1/24

参考文献

[1] FUTTERMAN W I. Dispersive body waves[J]. Journal of Geophysical Research, 1962, 67(13):5279-5291. doi:10.1029/JZ067i013p05279
[2] HALE D. An inverse Q-filter[R]. San Francisco:Stanford University, 1981.
[3] HARGREAVES N D, CALVERT A J. Inverse Q-filtering by fourier transform[J]. Geophysics, 1991, 56(4):519-527. doi:10.1190/1.1443067
[4] WANG Yanghua. A stable and efficient approach of inverse Q filtering[J]. Geophysics, 2002, 67(2):657-663. doi:10.1190/1.1468627
[5] WANG Yanghua. Quantifying the effectiveness of stabilized inverse Q filtering[J]. Geophysics, 2003, 68(1):337-345. doi:10.1190/1.1543219
[6] WANG Yanghua. Inverse Q-filter for seismic resolution enhancement[J]. Geophysics, 2006, 71(3):V51-V60. doi:10.1190/1.2192912
[7] WANG Yanghua. Seismic inverse Q filtering[M]. Oxford:Blackwell Publishing, 2008.
[8] WANG Yanghua. Stable Q analysis on vertical seismic profiling data[J]. Geophysics, 2014, 79(4):D217-D225. doi:10.1190/geo2013-0273.1
[9] YAN Hongyong, LIU Yang. Estimation of Q and inverse Q filtering for prestack reflected PP and converted PS-waves[J]. Applied Geophysics, 2009, 6(1):59-69. doi:10.1007/s11770-009-0009-y
[10] 严红勇,刘洋,于茜. 利用VSP提高地面地震资料的分辨率[J]. 石油天然气学报,2010,32(3):65-68. doi:10.3969/j.issn.1000-9752.2010.03.014 YAN Hongyong, LIU Yang, YU Qian. Improving the resolution of poststack surface seismic data by VSP[J]. Journal of Oil and Gas Technology, 2010, 32(3):65-68. doi:10.3969/j.issn.1000-9752.2010.03.014
[11] 严红勇,刘洋,赵前华,等. 一种提高VSP分辨率的反Q滤波方法[J]. 石油地球物理勘探,2011,46(6):873-880. YAN Hongyong, LIU Yang, ZHAO Qianhua, et al. A method for improving VSP resolution by inverse Q filtering[J]. Oil Geophysical Prospecting, 2011, 46(6):873-880.
[12] 刘财,冯晅,张瑾. 稳定的迭代法反Q滤波[J].石油地球物理勘探,2013,48(6):890-895. LIU Cai, FENG Xuan, ZHANG Jin. A stable inverse Q filtering using the iterative filtering method[J]. Oil Geophysical Prospecting. 2013, 48(6):890-895.
[13] 张固澜,贺振华,王熙明,等. 地震波频散效应与反Q滤波相位补偿[J]. 地球物理学报,2014(5):1655-1663. doi:10.6038/cjg20140528 ZHANG Gulan, HE Zhenhua, WANG Ximing, et al. Seismic wave dispersion effects and inverse Q-filter phase compensation[J]. Oil Geophysical Prospecting, 2014(5):1655-1663. doi:10.6038/cjg20140528
[14] 张固澜,林进,王熙明,等. 一种自适应增益限的反Q滤波[J]. 地球物理学报,2015,58(7):2525-2535. doi:10.6038/cjg20150727 ZHANG Gulan, LIN Jin, WANG Ximing, et al. A self-adaptive approach for inverse Q-filtering[J]. Chinese Journal of Geophysics, 2015, 58(7):2525-2535. doi:10.6038/cjg20150727
[15] 张固澜. 基于反Q滤波和IGST技术的地震高分辨处理研究[D]. 成都:成都理工大学,2016. ZHANG Gulan. High-resolution seismic data processing methods based on inverse Q-filter and IGST[D]. Chengdu:Chengdu University of Technology, 2016.
[16] 施羽,黄饶,周怀来,等. 一种变增益限反Q滤波法[C].北京:2017中国地球科学联合学术年会,2017. SHI Yu, HUANG Rao, ZHOU Huailai, et al. A variable gain-limited inverse-Q filtering method[C]. Beijing:2017 China Geosciences Joint Academic Annual Conference, 2017.
[17] 牛聪,施羽,周怀来,等. 基于等效加权模型的单步延拓反Q滤波方法研究及应用[J]. 地球物理学进展,2019,34(1):174-178. doi:10.6038/pg2019BB0583 NIU Cong, SHI Yu, ZHOU Huailai, et al. Research and application of one-step extended inverse Q filtering method based on equivalent weighted model[J]. Progress in Geophysics. 2019, 34(1):174-178. doi:10.6038/pg2019BB0583
[18] 张舒,褚艳利. GPU高性能运算之CUDA[M]. 北京:中国水利水电出版社,2009. ZHANG Shu, CHU Yanli. CUDA for GPU high-performance computing[M]. Beijing:China Water Resources and Hydropower Press, 2009.
[19] 刘钦,佟小龙. GPU/CPU协同并行计算(CPPC)在地震勘探资料处理中的应用[R]. 北京:北京吉星吉达公司,2008. LIU Qin, TONG Xiaolong. Application of GPU/CPU cooperative parallel computing (CPPC) in seismic exploration data processing[R]. Beijing:Beijing Jixing Jeddah Company, 2008.
[20] CEVAHIR A, NUKADA A, MATSUOKA S. Fast conjugate gradients with multiple GPUs[C]. Los Angeles:Computational Science-ICCS 2009, 9th International Conference, 2009. doi:10.1007/978-3-642-01970-8_90
[21] 张清,谢海波,赵开勇,等. 一种使用GPU加速地震叠前时间偏移的方法[J]. 信息技术与网络安全,2011,30(10):87-90. doi:10.3969/j.issn.1674-7720.2011.10.027 ZHANG Qing, XIE Haibo, ZHAO Kaiyong, et al. A method of accelerating pre-stack time migration using GPU[J]. Information Technology and Network Security, 2011, 30(10):87-90. doi:10.3969/j.issn.1674-7720.2011.10.027
[22] 袁崇鑫,邓飞. 基于GPU的波动方程正演模拟的实现[J]. 电脑知识与技术,2014,10(18):4333-4337. YUAN Chongxin, DENG Fei. GPU-based wave equation forward modeling to achieve[J]. Computer Knowledge and Technology, 2014, 10(18):4333-4337.
[23] 吴吉忠,左虎. 叠前衰减补偿时间偏移及GPU实现[J]. 石油地球物理勘探,2019,54(1):84-92. doi:10.13810/j.cnki.issn.1000-7210.2019.01.010 WU Jizhong, ZUO Hu. Attenuation compensation in prestack time migration and its GPU implementation[J]. Oil Geophysical Prospecting, 2019, 54(1):84-92. doi:10.13810/j.cnki.issn.1000-7210.2019.01.010
[24] 张全,林柏栎,彭博,等. 基于CUDA的地震相干体并行算法[J]. 地质与勘探,2020,56(1):147-153. doi:10.12134/j.dzykt.2020.01.014 ZHANG Quan, LIN Baiyue, PENG Bo, et al. Seismic coherence parallel algorithm based on CUDA[J]. Geology and Exploration, 2020, 56(1):147-153. doi:10.12134/j.dzykt.2020.01.014
[25] 张全,鲍华,饶长辉,等. GPU平台二维快速傅里叶变换算法实现及应用[J]. 光电工程,2016,43(2):69-75. doi:10.3969/j.issn.1003-501X.2016.02.012 ZHANG Quan, BAO Hua, RAO Changhui, et al. Realization and application of two-dimensional fast fourier transform algorithm based on GPU[J]. Opto-Electronic Engineering, 2016, 43(2):69-75. doi:10.3969/j.issn.1003-501X.2016.02.012
[26] VOLKOV V, KAZIAN B. Fitting FFT onto the G80 architecture[M]. Berkeley:University of California, 2008.
[27] 赵丽丽,张盛兵,张萌,等. 基于CUDA的高速FFT计算[J]. 计算机应用研究,2011,28(4):1556-1559. doi:10.3969/j.issn.1001-3695.2011.04.100 ZHAO Lili, ZHANG Shengbing, ZHANG Meng, et al. High performance FFT computation based on CUDA[J]. Application Research of Computers, 2011, 28(4):1556-1559. doi:10.3969/j.issn.1001-3695.2011.04.100

地震波能量补偿的并行反Q滤波方法研究

Research on Parallel Inverse Q Filtering Methods for Seismic Wave Energy Compensation

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 2

编辑推荐

Metrics

本文评价

[1]	伍轶鸣李勇李保柱姚军 . 双重介质油藏数值模拟并行算法研究[J]. 西南石油大学学报(自然科学版), 2009, 31(6): 80-84.
[2]	易宗富雷晓. GFT的性质与快速算法[J]. 西南石油大学学报(自然科学版), 1993, 15(1): 1-9.