The methodology of a post-stack stochastic seismic inversion with the co-constraint of deterministic inversion
ZHANG Fengqi1,2,3, LIU Junzhou1,2,3, LIU Lanfeng1,2,3, SHI Lei1,2,3, HAN Lei1,2,3
1. State Key Laboratory of Shale Oil and Gas Enrichment Mechanisms and Effective Development, Beijing 100083, China; 2. Sinopec Key Laboratory of Seismic Elastic Wave Technology, Beijing 100083, China; 3. Sinopec Petroleum Exploration and Production Research Institute, Beijing 102206, China
Abstract:Based on the previous work, a post-stack stochastic seismic inversion with the co-constraint of deterministic inversion was proposed. This new algorithm centered on sequential Gibbs sampling and an extended Metropolises-Hastings (M-H) algorithm. Owing to the introduction of sequence stratigraphy grids, geostatistics, structures, and sedimentary modes were integrated into the stochastic seismic inversion adaptively, thereby saving the effort of calculating the local dip of the strata or implementing complicated coordinate transformations. The high-frequency components (HFCs) in the stochastic seismic inversion result still featured large uncertainty, which resulted in big differences in reservoir characteristics among different implementations of stochastic seismic inversion. Given this problem, this paper combined sequential Gibbs sampling with collocated cokriging and introduced the co-constraint of deterministic inversion to restrict the candidate solution space and ultimately to reduce the uncertainty of the HFCs in the stochastic seismic inversion result. The following conclusions were obtained: ① Compared with deterministic inversion, stochastic seismic inversion can produce high-resolution results. The vertical resolution of the results is affected by the vertical variation whereas the lateral continuity of the results is affected by the lateral variation. ② Compared with uniform seismic grids, sequence strati-graphy grids are more suitable for stochastic inversion with a lateral variation constraint because of their integration with structures and sedimentary modes. With the help of a resampling matrix, sample simulation is carried out on sequence stratigraphy grids and convolutional forward modeling is carried out on uniform seismic grids. The whole inversion process not only meets the constraints of the structures and sedimentary modes but also conforms to geophysical principles. ③ Due to the introduction of the co-constraint of deterministic inversion, the candidate solution space is further restricted and the correlation between the results of stochastic inversion and deterministic inversion is enhanced. Furthermore, the uncertainty of the HFCs in the stochastic inversion results is reduced. Trial calculations with real data reveal that this new algorithm is verified by a comparison among four implementations of stochastic inversion.
张丰麒, 刘俊州, 刘兰锋, 时磊, 韩磊. 确定性反演协同约束的叠后随机地震反演方法[J]. 石油地球物理勘探, 2021, 56(5): 1137-1149.
ZHANG Fengqi, LIU Junzhou, LIU Lanfeng, SHI Lei, HAN Lei. The methodology of a post-stack stochastic seismic inversion with the co-constraint of deterministic inversion. Oil Geophysical Prospecting, 2021, 56(5): 1137-1149.
Joumel A, Alabert F.New method for reservoir mapping[J]. Journal of Petroleum Technology, 1990, 42(2):212-218.
[2]
Hass A, Dubrule O.Geostatistical inversion:a sequential method for stochastic reservoir modeling constrained by seismic data[J]. First Break, 1994, 12(11):561-569.
[3]
Dubrule O.Geostatistical reservoir characterization constrained by 3D seismic data[J]. Petoleum Gesocience, 1998, 4(2):121-128.
[4]
Debeye H, Sabbah E.Stochastic inversion[C]. SEG Technical Program Expanded Abstracts, 1996, 15:1212-1215.
[5]
Sams M S, Atkins D, Siad N, et al.Stochastic inversion for high resolution reservoir characterization in the Central Sumatra Basin[C]. SPE Asia Pacific Improved Oil Recovery Conference, 1999, SPE-57260-MS.
[6]
Mosegaard K, Tarantola A.Monte Carlo sampling of solutions to inverse problems[J]. Journal of Geophy-sical Research, 1995, 100(B7):12431-12447.
[7]
Hansen T M, Cordua K S, Mosegaard K.Inverse problems with non-trivial priors:efficient solution through sequential Gibbs sampling[J]. Computational Geosciences, 2012, 16(3):593-611.
张广智, 王丹阳, 印兴耀, 等.基于MCMC的叠前地震反演方法研究[J]. 地球物理学报, 2011, 54(11):2926-2932.ZHANG Guangzhi, WANG Danyang, YIN Xingyao, et al. Study on prestack seismic inversion using Markov Chain Monte Carlo[J]. Chinese Journal of Geophy-sics, 2011, 54(11):2926-2932.
[10]
王保丽, 印兴耀, 丁龙翔, 等.基于FFT-MA谱模拟的快速随机反演方法研究[J]. 地球物理学报, 58(2):664-673.WANG Baoli, YIN Xingyao, DING Longxiang, et al.Study of fast stochastic inversion based on FFT-MA spectrum simulation[J]. Chinese Journal of Geophy-sics, 2015, 58(2):664-673.
[11]
刘兴业, 李景叶, 陈小宏, 等.联合多点地质统计学与序贯高斯模拟的随机反演方法[J]. 地球物理学报, 2018, 61(7):2998-3007.LIU Xingye, LI Jingye, CHEN Xiaohong, et al.A stochastic inversion method integrating multi-point geostatistics and sequential Gaussian simulation[J]. Chinese Journal of Geophysics, 2018, 61(7):2998-3007.
[12]
赵晨, 张广智, 张佳佳, 等.基于Metropolis优化的叠前全局迭代地质统计学反演方法[J]. 地球物理学报, 2020, 63(8):3116-3130.ZHAO Chen, ZHANG Guangzhi, ZHANG Jiajia, et al. Prestack global iteration geostatistical inversion method based on metropolis sampling algorithm[J]. Chinese Journal of Geophysics, 2020, 63(8):3116-3130.
[13]
张丰麒, 金之钧, 盛秀杰, 等.贝叶斯三参数低频软约束同步反演[J]. 石油地球物理勘探, 2016, 51(5):965-975.ZHANG Fengqi, JIN Zhijun, SHENG Xiujie, et al. Bayesian prestack three-term inversion with soft low-frequency constraint[J]. Oil Geophysical Prospecting, 2016, 51(5):965-975.
[14]
姚凌青, 潘懋, 成秋明, 等.三维Kriging方法中的变异函数套合[J]. 地球科学——中国地质大学学报, 2009, 34(2):294-298.YAO Lingqing, PAN Mao, CHENG Qiuming, et al.Nested overlap of variograms in 3D Kriging[J]. Earth Science:Journal of China University of Geosciences, 2009, 34(2):294-298.
[15]
Buland A, Omre H.Bayesian linearized AVO inversion[J]. Geophysics, 2003, 68(1):185-198.
[16]
Xu W, Tran T, Srivastava M, et al.Integrating seismic data in reservoir modelling:The collocated cokriging alternative[C]. The 67th Annual Technical Conference and Exhibition of the Society of Petroleum Engineers, 1992, SPE-24742-MS.
[17]
Soares A.Direct sequential simulation and cosimulation[J]. Mathematical Geology, 2001, 33(8):911-926.
[18]
周爽爽, 印兴耀, 裴松, 等.地震波形约束的蒙特卡洛-马尔科夫链随机反演方法[J]. 石油地球物理勘探, 2021, 56(3):543-554, 592.ZHOU Shuangshuang, YIN Xingyao, PEI Song, et al. Monte Carlo-Markov Chain stochastic inversion constrained by seismic waveform[J]. Oil Geophysical Prospecting, 2021, 56(3):543-554, 592.