|
|
Prestack seismic stochastic inversion method based on spatial co-simulation |
CAO Yamei1,2, ZHOU Hui1,2, YU Bo3, ZHANG Yuangao4, DAI Shili4 |
1. National Key Laboratory of Petroleum Resources and Engineering, China University of Petroleum(Beijing), Beijing 102249, China; 2. CNPC Key Laboratory of Geophysical Exploration, China University of Petroleum(Beijing), Beijing 102249, China; 3. School of Earth Sciences, Northeast Petroleum University, Daqing, Heilongjiang 163318, China; 4. Exploration Department of Daqing Oilfield Company Ltd, Daqing, Heilongjiang 163453, China |
|
|
Abstract Sequential stochastic simulation is generally used to characterize reservoir heterogeneity both in the iterative stochastic inversion and linear Bayesian stochastic inversion. Most sequential simulation methods rely on variograms or training images to describe the spatial correlation of model parameters. In addition, the simulation results are required to be calculated point by point, which makes parallel computing difficult and reduces computational efficiency. Therefore, a conditioned fast Fourier transform moving average (FFT-MA) is introduced into the linear inversion framework, and a prestack seismic stochastic inversion method based on spatial co-simulation is proposed. Firstly, the posterior probability distribution of elastic parameters is obtained by integrating seismic data and low-frequency well logging information under the Bayesian framework. Then, the probability field is generated according to the FFT-MA algorithm. Well logging data is taken as conditional data to conduct Bayesian posterior probability field co-simulation. High-resolution prestack stochastic inversion results of elastic parameters constrained by well logging and seismic data are thus obtained. No iteration and update of model parameters is required by the method, which greatly improves the computational efficiency of stochastic inversion. Finally, the validity of the proposed method is demonstrated by numerical model examples and practical data application cases. The numerical model examples show that the proposed method has significant advantages over conventional methods in terms of high-resolution reservoir prediction and computational efficiency. Small-scale reservoir characterizations can be explored stably and accurately. The practical data application cases show that the high-resolution inversion results obtained by the proposed method match well with well logging data. The practicality of stochastic inversion in characterizing quantitatively thin reservoirs is greatly improved.
|
Received: 28 March 2023
|
|
|
|
|
[1] |
撒利明,杨午阳,姚逢昌,等.地震反演技术回顾与展望[J].石油地球物理勘探,2015,50(1):184-202.SA Liming,YANG Wuyang,YAO Fengchang,et al.Past,present,and future of geophysical inversion[J].Oil Geophysical Prospecting,2015,50(1):184-202.
|
[2] |
DOWNTON J E.Seismic Parameter Estimation from AVO Inversion[D].Department of Geology and Geophysics,University of Calgary,2005.
|
[3] |
TIHONOV A N.Solution of incorrectly formulated problems and the regularization method[J].Soviet Math Dokl,1963,DOI:https://doi.org/10.11499/sicejl1962.36.468.
|
[4] |
裴松,印兴耀,李坤.全域正则化快速匹配追踪稀疏地震反演方法[J].石油地球物理勘探,2022,57(6):1400-1408,1426.PEI Song,YIN Xingyao,LI Kun.Sparse seismic inversion method based on full-domain regularized fast matching pursuit[J].Oil Geophysical Prospecting,2022,57(6):1400-1408,1426.
|
[5] |
ZHANG R,SEN M K,PHAN S,et al.Stochastic and deterministic seismic inversion methods for thin-bed resolution[J].Journal of Geophysics and Engineering,2012,9(5):611-618.
|
[6] |
张枫,贾学成,张晓敏,等.相控反演薄储层预测技术在鄂尔多斯盆地东缘的应用[J].石油地球物理勘探,2022,57(增刊1):215-222.ZHANG Feng,JIA Xuecheng,ZHANG Xiaomin,et al.Application of thin reservoir prediction technology based on facies-controlled inversion at eastern margin of Ordos Basin[J].Oil Geophysical Prospecting,2022,57(S1):215-222.
|
[7] |
郭同翠,姜明军,纪迎章,等.叠前地质统计学反演在页岩甜点和薄夹层预测中的应用——以西加拿大盆地W区块为例[J].石油地球物理勘探,2020,55(1):167-175.GUO Tongcui,JIANG Mingjun,JI Yingzhang,et al.The application of prestack geostatistical inversion in the prediction of shale sweet spots and thin interbeds:a case study of Block W in Western Canada Basin[J].Oil Geophysical Prospecting,2020,55(1):167-175.
|
[8] |
PEREIRA P,BORDIGNON F,AZEVEDO L,et al.Strategies for integrating uncertainty in iterative geostatistical seismic inversion[J].Geophysics,2019,84(2):R207-R219.
|
[9] |
HAAS A,DUBRULE O.Geostatistical inversion:a sequential method of stochastic reservoir modeling constrained by seismic data[J].First Break,1994,12(11):561-569.
|
[10] |
郭强,雒聪,刘红达,等.自适应优化参数模拟退火的叠前地震联合反演方法[J].石油地球物理勘探,2023,58(3):670-679.GUO Qiang,LUO Cong,LIU Hongda,et al.Prestack seismic hybrid inversion based on simulated annealing algorithm with adaptive optimization parameters[J].Oil Geophysical Prospecting,2023,58(3):670-679.
|
[11] |
PEREIRA Â,NUNES R,AZEVEDO L,et al.Geostatistical seismic inversion for frontier exploration[J].Interpretation,2017,5(4):T477-T485.
|
[12] |
LIU M,GRANA D.Stochastic nonlinear inversion of seismic data for the estimation of petroelastic properties using the ensemble smoother and data reparamete- rization[J].Geophysics,2018,83(3):M25-M39.
|
[13] |
YU B,ZHOU H,WANG L,et al.Prestack Bayesian statistical inversion constrained by reflection features[J].Geophysics,2020,85(4):R349-R363.
|
[14] |
HANSEN T M,JOURNEL A G,TARANTOLA A,et al.Linear inverse Gaussian theory and geostatistics[J].Geophysics,2006,71(6):R101-R111.
|
[15] |
YU B,ZHOU H,WANG L,et al.Prestack Bayesian linearized inversion with decorrelated prior information[J].Mathematical Geosciences,2021,53(3):437-464.
|
[16] |
DE FIGUEIREDO L P,GRANA D,ROISENBERG M,et al.Gaussian mixture Markov chain Monte Carlo method for linear seismic inversion[J].Geophysics,2019,84(3):R463-R476.
|
[17] |
李祺鑫,罗亚能,张生,等.高分辨率波阻抗贝叶斯序贯随机反演[J].石油地球物理勘探,2020,55(2):389-397.LI Qixin,LUO Yaneng,ZHANG Sheng,et al.High-resolution Bayesian sequential stochastic inversion[J].Oil Geophysical Prospecting,2020,55(2):389-397.
|
[18] |
樊鹏军,马良涛,王宗俊,等.地质统计学反演中变差函数地质含义及求取方法探讨[J].地球物理学进展,2017,32(6):2444-2450.FAN Pengjun,MA Liangtao,WANG Zongjun,et al.Variogram geological implication and its calculating method discussing for geostatistical inversion[J].Progress in Geophysics,2017,32(6):2444-2450.
|
[19] |
LIU X,LI J,CHEN X,et al.Stochastic inversion of facies and reservoir properties based on multi-point geostatistics[J].Journal of Geophysics and Engineering,2018,15(6):2455-2468.
|
[20] |
LE RAVALEC M,NOETINGER B,HU L Y.The FFT moving average (FFT-MA) generator:an efficient numerical method for generating and conditioning Gaussian simulations[J].Mathematical Geology,2000,32(6):701-723.
|
[21] |
印兴耀,刘婵娟,王保丽.基于混合遗传算法的叠前随机反演方法[J].中国石油大学学报(自然科学版),2017,41(4):65-70.YIN Xingyao,LIU Chanjuan,WANG Baoli.Pre-stack stochastic inversion based on hybrid genetic algorithm[J].Journal of China University of Petroleum (Edition of Natural Science),2017,41(4):65-70.
|
[22] |
杨修伟,朱培民,毛宁波,等.基于FFT-MA的随机介质建模方法[J].地球物理学报,2018,61(12):5007-5018.YANG Xiuwei,ZHU Peimin,MAO Ningbo,et al.Random medium modeling based on FFT-MA[J].Chinese Journal of Geophysics,2018,61(12):5007-5018.
|
[23] |
DE FIGUEIREDO L P,SCHMITZ T,LUNELLI R,et al.Direct multivariate simulation-a stepwise conditional transformation for multivariate geostatistical simulation[J].Computers & Geosciences,2021,DOI:10.1016/j.cageo.2020.104659.
|
[24] |
王保丽,印兴耀,丁龙翔,等.基于FFT-MA谱模拟的快速随机反演方法研究[J].地球物理学报,2015,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.
|
[25] |
AKI K,RICHARDS P G.Quantitative Seismology (2nd Edition)[M].University Science Books,2002.
|
[26] |
BORDIGNON F L,DE FIGUEIREDO L P,AZEVEDO L,et al.Hybrid global stochastic and Bayesian linearized acoustic seismic inversion methodo- logy[J].IEEE Transactions on Geoscience and Remote Sensing,2017,55(8):4457-4464.
|
[27] |
OLIVER D S.Moving averages for Gaussian simulation in two and three dimensions[J].Mathematical Geology,1995,27(8):939-960.
|
|
|
|