基于模拟退火粒子群优化算法的裂缝型储层各向异性参数地震反演

doi:10.13278/j.cnki.jjuese.20180255

吉林大学学报(地球科学版) ›› 2019, Vol. 49 ›› Issue (5): 1466-1476.doi: 10.13278/j.cnki.jjuese.20180255

基于模拟退火粒子群优化算法的裂缝型储层各向异性参数地震反演

罗腾¹, 冯晅¹, 郭智奇¹, 刘财¹, 刘喜武^2,3,4

1. 吉林大学地球探测科学与技术学院, 长春 130026;
2. 页岩油气富集机理与有效开发国家重点实验室, 北京 100083;
3. 中国石化页岩油气勘探开发重点实验室, 北京 100083;
4. 中国石化石油勘探开发研究院, 北京 100083

收稿日期:2018-10-08 出版日期:2019-09-26 发布日期:2019-10-10
通讯作者: 郭智奇(1980-),男,教授,博士生导师,主要从事岩石物理、地震各向异性、油气储层预测等研究,E-mail:guozhiqi@jlu.edu.cn E-mail:guozhiqi@jlu.edu.cn
作者简介:罗腾(1989-),男,博士研究生,主要从事地震各向异性正反演研究,E-mail:autumn2009lt@126.com
基金资助:
国家"十三五"重大专项（2017ZX05049-002）；国家自然科学基金石油化工联合基金项目（U1663207）；国家自然科学基金项目（41430322）

Seismic Inversion of Anisotropy Parameters of Fractured Reservoirs by Simulated Annealing and Particle Swarm Optimization

Luo Teng¹, Feng Xuan¹, Guo Zhiqi¹, Liu Cai¹, Liu Xiwu^2,3,4

1. College of GeoExploration Science and Technology, Jilin University, Changchun 130026, China;
2. National Key Laboratory of Corporation of Shale Oil/Gas Enrichment Mechanism and Effective Development, Beijing 100083, China;
3. SinoPEC Key Laboratory of Shale Oil/Gas Exploration and Production Technology, Beijing 100083, China;
4. SinoPEC Petroleum Exploration and Production Research Institute, Beijing 100083, China

Received:2018-10-08 Online:2019-09-26 Published:2019-10-10
Supported by:
Supported by National Key Projects During the 13th Five-Year Plan of China (2017ZX05049-002), NSFC and SinoPEC Joint Key Project (U1663207) and National Natural Science Foundation of China (41430322)

摘要/Abstract

摘要： 发育垂直定向排列裂缝的地下岩石可等效为具有水平对称轴的横向各向同性（horizontal transverse isotropic，HTI）介质。针对HTI介质模型，本文研究了裂缝型储层的各向异性参数地震振幅随方位角变化（amplitude variations with azimuth，AVAZ）的反演方法。首先，在地震AVAZ反演流程中，提出采用模拟退火粒子群优化算法实现裂缝型储层各向异性参数反演。之后，通过理论模型测试，验证了基于模拟退火粒子群优化算法的地震AVAZ反演的有效性。最后，将反演方法应用于四川盆地龙马溪组页岩气储层实际方位地震数据；在反演之前先利用傅里叶级数方法估计裂缝方位并对实际数据进行方位校正，以提供更准确的输入数据；通过计算得到的P波、S波各向异性参数可用于评价裂缝发育程度。反演结果表明，研究区域构造顶部裂缝较发育，与地质基本理论一致，验证了反演方法的合理性。

关键词: HTI, 裂缝, 储层, 各向异性, 地震AVAZ, 反演, 四川盆地, 页岩气

Abstract: Underground rocks with vertical aligned fractures can be equivalent to horizontal transversely isotropic (HTI) media. For HTI medium model, the authors studied the inversion method to estimate anisotropic parameter using the seismic amplitude varying with azimuth (AVAZ) of fractured reservoirs. Firstly, in the seismic AVAZ inversion process, the simulated annealing particle swarm optimization algorithm was proposed to invert the anisotropic parameters of fractured reservoirs; then, the validity of seismic AVAZ inversion based on simulated annealing particle swarm optimization algorithm was verified by theoretical model test; finally, the inversion method was applied to the field azimuth seismic data of Longmaxi shale gas reservoir in Sichuan basin. Before inversion, the Fourier series method was used to estimate the fracture orientation and correct the field data to provide more accurate input data. The anisotropic parameters of P and S waves obtained by calculation can be used to evaluate the degree of fracture development. The inversion results show that the fractures on the top of the study area are well developed, which is consistent with the basic geological theory. It verifies the rationality of the inversion method and proves that the inversion results can provide a basis for reservoir fracture identification.

Key words: HTI, fracture, reservoirs, anisotropy, seismic AVAZ, inversion, Sishuan basin, shale gas

中图分类号:

P631.4

罗腾, 冯晅, 郭智奇, 刘财, 刘喜武. 基于模拟退火粒子群优化算法的裂缝型储层各向异性参数地震反演[J]. 吉林大学学报(地球科学版), 2019, 49(5): 1466-1476.

Luo Teng, Feng Xuan, Guo Zhiqi, Liu Cai, Liu Xiwu. Seismic Inversion of Anisotropy Parameters of Fractured Reservoirs by Simulated Annealing and Particle Swarm Optimization[J]. Journal of Jilin University(Earth Science Edition), 2019, 49(5): 1466-1476.

参考文献

[1] Crampin S. Effective Anisotropic Elastic Constants for Wave Propagation Through Cracked Solids[J]. Geophysical Journal of the Royal Astronomical Society, 1983, 76(1):135-145.
[2] Liu E, Martinez A. Seismic Fracture Characterization[M]. Houten:EAGE Publications, 2012.
[3] 张世俊,杨慧珠,董渊,等. 遗传算法反演HTI介质各向异性参数[J].石油地球物理勘探,2002, 37(1):24-28. Zhang Shijun, Yang Huizhu, Dong Yuan, et al. Inversion of Anisotropic Parameters in HTI Medium by Genetic Algorithm[J]. Oil Geophysical Prospecting, 2002, 37(1):24-28.
[4] 杜启振,董渊,杨慧珠. HTI介质横波正常时差速度反演[J]. 中国石油大学学报(自然科学版),2002, 26(2):26-30. Du Qizhen, Dong Yuan, Yang Huizhu. Inversion of Normal Time Difference Velocity of Shear Wave in HTI Medium[J]. Journal of the University of Petroleum (Nature Science Edition), 2002, 26(2):26-30.
[5] 孙炜,何治亮,李玉凤,等. 基于HTI介质各向异性正演的裂缝预测属性优选[J]. 石油物探,2014, 53(2):223-231. Sun Wei, He Zhiliang, Li Yufeng, et al.Seismic Attributes Optimization for Fractures Prediction Based on the Anisotropy Forwarding of the HTI Media[J]. Geophysical Prospecting for Petroleum, 2014, 53(2):223-231.
[6] 窦喜英,韩立国,刘春成,等. 裂缝型HTI介质中的弹性阻抗[J]. 吉林大学学报(地球科学版),2012,42(4):1192-1198. Dou Xiying, Han Liguo, Liu Chuncheng, et al. Elastic Impedance in Fractured HTI Media[J]. Journal of Jilin University (Earth Science Edition), 2012, 42(4):1192-1198.
[7] 范晓敏,李舟波. 裂缝性碳酸盐岩储层声波时差曲线的波动和增幅分析[J]. 吉林大学学报(地球科学版),2007, 37(1):168-173. Fan Xiaomin, Li Zhoubo. Analysis of Fluctuation and Amplitude Increase in Sonic Logs Caused by Fractures in Carbonate Rocks[J]. Journal of Jilin University (Earth Science Edition), 2007, 37(1):168-173.
[8] 杜启振,孔丽云,韩世春. 裂缝诱导各向异性双孔隙介质波场传播特征[J]. 地球物理学报,2009, 52(4):1049-1058. Du Qizhen, Kong Liyun, Han Shichun. Wavefield Propagation Characteristics in the Fracture-Induced Anisotropic Double-Porosity Medium[J]. Chinese Journal of Geophysics, 2009, 52(4):1049-1058.
[9] 李春鹏,印兴耀,刘志国,等. 裂缝型储层预测的各向异性梯度反演方法研究[J]. 石油物探,2017, 56(6):835-840. Li Chunpeng, Yin Xingyao, Liu Zhiguo, et al. An Anisotropic Gradient Inversion for Fractured Reservoir Prediction[J]. Geophysical Prospecting for Petroleum, 2017, 56(6):835-840.
[10] Thomsen L. Weak Elastic Anisotropy[J]. Geophysics, 1986, 51(10):1954-1966.
[11] Ruger A. P-Wave Reflection Coefficients for Transversely Isotropic Models with Vertical and Horizontal Axis of Symmetry[J]. Geophysics, 1997, 62:713-722.
[12] Rüger A. Variation of P-Wave Reflectivity with Offset and Azimuth in Anisotropic Media[J]. Geophysics, 1998, 63(63):935-947.
[13] Mallick S, Meister L J, Chambers R E, et al. Determination of the Principal Directions of Azimuthal Anisotropy from P-Wave Seismic Data[J]. Geophysics, 1998, 63(2):692-706.
[14] 朱培民,王家映,於文辉,等. 用纵波AVO数据反演储层裂隙密度参数[J]. 石油物探,2001, 40(2):1-12. Zhu Peimin, Wang Jiaying, Yu Wenhui, et al. Inverting Reservoir Crack Density Using P-Wave AVO Data[J]. Geophysical Prospecting for Petroleum, 2001, 40(2):1-12.
[15] Bachrach R, Sengupta M, Salama A, et al. Reconstruction of the Layer Anisotropic Elastic Parameters and High-Resolution Fracture Characterization from P-Wave Data:A Case Study Using Seismic Inversion and Bayesian Rock Physics Parameter Estimation[J]. Geophysical Prospecting, 2009, 57(2):253-262.
[16] 张广智,陈怀震,印兴耀,等. 基于各向异性AVO的裂缝弹性参数叠前反演方法[J]. 吉林大学学报(地球科学版),2012, 42(3):845-851. Zhang Guangzhi, Chen Huaizhen, Yin Xingyao, et al. Method of Fracture Elastic Parameter Inversion Based on Anisotropic AVO[J]. Journal of Jilin University (Earth Science Edition), 2012, 42(3):845-851.
[17] 刘财,刘宇巍,冯晅,等. 基于方位相交的纵波AVA数据运用SVD反演HTI介质裂缝密度[J]. 吉林大学学报(地球科学版),2013, 43(5):1655-1662. Liu Cai, Liu Yuwei, Feng Xuan, et al. Invert Crack Density of HTI Media by Using SVD Based on PP-Wave AVA Data from Crossing Seismic Survey Lines[J]. Journal of Jilin University (Earth Science Edition), 2013, 43(5):1655-1662.
[18] 李博南,刘财,郭智奇. 基于等效介质模型和频变AVO反演的裂缝储层参数估算方法[J]. 吉林大学学报(地球科学版),2017, 47(1):234-244. Li Bonan, Liu Cai, Guo Zhiqi. Estimation of Fractured Reservoir Parameters Based on Equivalent Media Model and Frequncy-Dependent AVO Inversion[J]. Journal of Jilin University (Earth Science Edition), 2017, 47(1):234-244.
[19] Liu B, Yin X, Yang F. Anisotropic Parameters Inversion of Fractured Reservoir Based on Orthogonal Azimuth Reflection Coefficient Difference[C]//Beijing 2014 International Geophysical Conference & Exposition. Beijing:Society of Exploration Geophysicists, 2014:584-587.
[20] Bakulin A, Tsvankin I, Grechka V. Estimation of Fracture Parameters from Reflection Seismic Data:Part I:HTI Model Due to A Single Fracture Set[J]. Geophysics, 2000, 65(6):1788-1802.
[21] Kennedy J, Eberhart R. Particle Swarm Optimization[C]//International Conference on Neural Networks. Perth:IEEE, 2002:1942-1948.
[22] Kirkpatrick S, Gelatt C D, Vecchi M P. Optimization by Simulated Annealing[J]. Science, 1983(220):671-680.
[23] 高鹰,谢胜利. 基于模拟退火的粒子群优化算法[J]. 计算机工程与应用,2004, 40(1):47-50. Gao Ying, Xie Shengli. Particle Swarm Optimization Algorithms Based on Simulated Annealing[J]. Computer Engineering and Applications, 2004, 40(1):47-50.
[24] Rüger A. Reflection Coefficients and Azimuthal AVO Analysis in Anisotropic Media[M]. Tulsa:Society of Exploration Geophysicists, 2002.
[25] Sayers C M, Dean S. Azimuth-Dependent AVO in Reservoirs Containing Non-Orthogonal Fracture Sets[J]. Geophysical Prospecting, 2010, 49(1):100-106.
[26] Ikelle L T. Amplitude Variations with Azimuths (AVAZ) Inversion Based on Linearized Inversion of Common Azimuthal Sections[M]. Tulsa:Society of Exploration Geophysicists, 1996.
[27] Downton J, Russell H. Azimuthal Fourier Coefficients:A Simple Method to Estimate Fracture Parameters[C]//Seg Technical Program Expanded. Tulsa:Society of Exploration Geophysicists,2011:269-273.
[28] 郭彤楼,张汉荣. 四川盆地焦石坝页岩气田形成与富集高产模式[J]. 石油勘探与开发,2014, 41(1):28-36. Guo Tonglou, Zhang Hanrong. Formation and Enrichment Mode of Jiaoshiba Shale Gas Field, Sichuan Basin[J]. Petroleum Exploration and Development, 2014, 41(1):28-36.

Metrics

Viewed

Full text

Abstract

Cited

Shared

Discussed

基于模拟退火粒子群优化算法的裂缝型储层各向异性参数地震反演

Seismic Inversion of Anisotropy Parameters of Fractured Reservoirs by Simulated Annealing and Particle Swarm Optimization

PDF (PC)

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

Metrics

本文评价

推荐阅读 10

[1]	唐华风, 张亚东, 邹明倬, 陆国超, 户景松, 张津铭, 卞海越. 渤海湾盆地莱州湾凹陷K油田中生界火山岩储层定量评价[J]. 吉林大学学报(地球科学版), 2026, 56(3): 739-753.
[2]	李冠, 何文祥, 文志刚, 胡勇, 高小洋, 肖娜. 鄂尔多斯盆地陇东地区长71-2亚段不同源储组合的页岩油富集机理及富集模式[J]. 吉林大学学报(地球科学版), 2026, 56(3): 768-787.
[3]	陈泽熙, 徐志明, 李斌, 钟笠, 彭军, 张昆, 魏祥峰, 郝景宇. 基质型页岩甜点特征及评价方法：以四川盆地涪陵地区中侏罗统凉高山组陆相页岩为例[J]. 吉林大学学报(地球科学版), 2026, 56(3): 788-803.
[4]	王蓓东, 杨胜来 , 赵帅 , 邓惠 , 张钰祥 , 白浩言, 李沁怡, 姜艺. 高磨地区震旦系灯影组四段台内储层储集空间类型及其孔喉结构特征[J]. 吉林大学学报(地球科学版), 2026, 56(3): 804-817.
[5]	夏英杰, 陈华斌, 陈健, 姚明宇, 杨海. 不同施工参数下深层页岩水力压裂数值模拟[J]. 吉林大学学报(地球科学版), 2026, 56(3): 924-936.
[6]	张恩威, 孟庆涛, 唐佰强, 胡菲, 党微. 基于机器学习的TOC测井预测方法：以松辽盆地南部青山口组一段为例[J]. 吉林大学学报(地球科学版), 2026, 56(3): 1062-1075.
[7]	丁梦颜, 冯晅, 刘财. 裂缝对岩石非线性弹性特征的影响[J]. 吉林大学学报(地球科学版), 2026, 56(3): 1051-1061.
[8]	赵天硕, 宋超, , 刘财, , , , 徐雨歆. 物理信息神经网络地震走时层析成像程函方程因式分解方法 [J]. 吉林大学学报(地球科学版), 2026, 56(2): 694-702.
[9]	刘涛, 刘宗堡, 张可佳, 张岩, 张瑞雪, 刘晓文, 徐翠云. 基于深度学习的致密砂岩储层薄片图像生成与识别方法 [J]. 吉林大学学报(地球科学版), 2026, 56(2): 724-738.
[10]	周三栋, 张伟鑫, 成巧耘, 魏佳成, 周国栋, . 深部煤储层中煤层水可动性评价方法及动力学机制研究进展[J]. 吉林大学学报(地球科学版), 2026, 56(2): 397-412.
[11]	王健, 张媛瑗, 吴楠, 徐清海, 崔子岳, 刘显凤, 付清萌. 鄂尔多斯盆地志靖—安塞地区长7段夹层型页岩油储层特征及分类评价[J]. 吉林大学学报(地球科学版), 2026, 56(2): 453-468.
[12]	王振军. 油气储层裂缝预测进展[J]. 吉林大学学报(地球科学版), 2026, 56(2): 469-482.
[13]	李新琦, 霍飞宇, 赵志平, 崔普媛, 李祺. 渤中凹陷中生界火山岩规模性裂缝储层特征及成因#br#[J]. 吉林大学学报(地球科学版), 2026, 56(2): 483-496.
[14]	黄玉龙, 张皓, 郭强, 陈昌. 低阻玄武岩的矿物学特征及其油气地质意义——以辽河坳陷古近系为例 [J]. 吉林大学学报(地球科学版), 2026, 56(2): 511-521.
[15]	吴小奇, 刘全有, 王萍, 黎华继, 朱东亚, 李朋朋. 沉积盆地天然气藏中氦气贫化机理及资源潜力：以四川盆地新场气田为例[J]. 吉林大学学报(地球科学版), 2025, 55(6): 1837-1850.