A posteriori regularization method for elastic impe-dance inversion
Yang Xiao1,2,3, Zhu Houqin4, Wang Yanfei1,2,3
1. Key Laboratory of Petroleum Resource, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China; 2. University of Chinese Academy of Sciences, Beijing 100049, China; 3. Institutions of Earth Science, Chinese Academy of Sciences, Beijing 100029, China; 4. Research Institute of Petroleum Exploration and Development, PetroChina, Beijing 100083, China
Abstract:The rock brittleness index is one of the key factors for shale gas reservoir evaluation.In comparison with the post stack acoustic impedance inversion and AVO inversion,the elastic impedance inversion can obtain more abundant,stable and reliable inversion results of elastic parameters,which is beneficial to the selection of shale sensitivity brittleness index.However,the elastic impedance inversion is an ill posed problem in the sense of Hadamard,that is,the existence,uniqueness and stability of the solution cannot be met at the same time.It is necessary to use the regularization me-thod and with appropriate optimization techniques to improve the stability and accuracy of solutions.In order to overcome the ill-posedness of the elastic impedance inversion,a Tikhonov regularization inversion optimization model based on L2 norm constraint is established,which is supplemented by a posteriori optimal regular parameter selection method,and a dual filter factor regularization algorithm for minimizing problems is proposed.Nume-rical experiments on theoretical model and real data show that the dual filter regularization algorithm is feasible and promising.
张大伟.《页岩气发展规划(2011-2015年)》解读.天然气工业,2012,32(4):6-8.Zhang Dawei.Important measures for promoting the shale gas exploration and development in China:An explanation of The Shale Gas Development Programming from 2011 to 2015.Natural Gas Industry,2012,32(4):6-8.
[2]
王健,石万忠,舒志国等.富有机质页岩TOC含量的地球物理定量化预测.石油地球物理勘探,2016,51(3):596-604.Wang Jian,Shi Wanzhong,Shu Zhiguo et al.TOC content quantitative prediction in organic-rich shale.OGP,2016,51(3):596-604.
[3]
苏建龙,屈大鹏,陈超等.叠前地震反演方法对比分析——焦石坝页岩气藏勘探实例.石油地球物理勘探,2016,51(3):581-588.Su Jianlong,Qu Dapeng,Chen Chao et al.Application of pre-stack inversion technique for shale gas in Jiaoshiba gas field.OGP,2016,51(3):581-588.
[4]
Matt B,Bill G.Special techniques tap shale gas.Exploration and Production in Hart Energy,2007,80(3):89-93.
[5]
Altindag R.Correlation of specific energy with rock brittleness concepts on rock cutting.Journal of the South African Institute of Mining and Metallurgy,2003,103(3):163-171.
[6]
Yagiz S.Assessment of brittleness using rock strength and density with punch penetration test.Tunnelling and Underground Space Technology,2009,24(1):66-74.
[7]
Goodway B,Perez M,Varsek J et al.Seismic petro-physics and isotropic-anisotropic AVO methods for unconventional gas exploration.The Leading Edge,2010,29(12):1500-1508.
[8]
Guo Z,Li X,Liu C et al.A shale rock physics model for analysis of brittleness index,mineralogy and porosity in the Barnet Shale.Journal of Geophysics and Engineering,2013,10(2):25006-25015.
[9]
Shuey R T.A simplification of the Zoeppritz equa-tions.Geophysics,1985,50(4):609-614.
[10]
Cambois G.Can P-wave AVO be quatitative? The Leading Edge,2000,19(11):1246-1251.
[11]
Connolly P.Elastic impedance.The Leading Edge,1999,18(4):438-452.
[12]
Shao M L,George A M.Elastic impedance inversion of multichannel seismic data from unconsolidated se-diments containing gas hydrate and free gas.Geophy-sics,2004,69(1):164-179.
[13]
David N W,Patrick A C,Roger L R et al.Extended elastic impedance for fluid and lithology prediction.Geo-physics,2002,67(1):63-67.
[14]
刘百红,李学云,李建华等.弹性阻抗反演在储层含气性预测中的应用.石油地球物理勘探,2012,47(增刊1):72-77.Liu Baihong,Li Xueyun,Li Jianhua et al.Gas-bearing reservoir prediction by elastic impedance inversion.OGP,2012,47(S1):72-77.
[15]
王彦飞,斯捷潘诺娃I E,提塔连科V N等.地球物理数值反演问题.北京:高等教育出版社,2011.
[16]
Backus G E,Gilbert J F.Numerical application of a formulism for geophysical inverse problem.Geophy-sics Journal International,1967,13(1-3):247-276.
Morozov V A.On regularization of ill-posed problems and selection of regularization parameter.Springer,New York,1984.
[20]
Tihhonov A N,Goncharsky A V.Ill-Posed Problems in the Natural Sciences.MIR,Moscow,1987.
[21]
Engl H W,Hanke M,Neubauer.Regularization of Inverse Problems.Kluwer,Dordrecht,1966.
[22]
王彦飞.反演问题的计算方法及其应用.北京:高等教育出版社,2007.
[23]
Tikhonov A N,Arsenin V Y.Solutions of Ill-Posed Problems.Division of John Wiley and Sons,Halsted Press,Washington D C,1977.
[24]
Hanke M,Hansen P C.Regularization methods for large-scale problems.Surveys on Mathematics for Industry,1993,3(4):253-315.
[25]
Golub G H,Heath M,Wahba G.Generalized cross-validation as a method for choosing a good ridge parameter.Technometrics,1979,21(2):215-223.
[26]
李振春,李闯,黄建平等.基于先验模型约束的最小二乘逆时偏移方法.石油地球物理勘探,2016,51(4):738-744.Li Zhenchun,Li Chuang,Huang Jianping et al.Regularized least-squares reverse time migration with prior model.OGP,2016,51(4):738-744.
[27]
Hansen P C.Analysis of discrete ill-posed problems by means of the L-curve.SIAM Review,1992,34(4):561-580.
[28]
Hansen P C,O' Lear D P.The use of the L-curve in the regularization of discrete ill-posed problems.SIAM Journal of Scientific Computations,1993,14(6):1487-1503.
[29]
Castagna J P,Smith S W.Comparison of AVO indicators:a modeling study.Geophysics,1994,59(12):1849-1855.