储层重力密度反演后验约束正则化方法
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
本文针对蒸汽辅助重力泄油(SAGD)生产中开发监测问题,发展了综合应用地震及重力数据反演储层密度的联合反演算法.通过测井数据建立纵波阻抗与密度的直接关系,并推导出这种关系下重力与纵波阻抗数据联合反演的计算方法,从而计算出蒸汽腔体密度分布规律.文中应用密度反演后验约束正则化方法,采用Tikhonov正则化模型,通过波阻抗数据作为约束进行联合反演,在算法上提高了稳定性,同时得到较高的反演精度.文中对SAGD生产中的理论模型进行了方法试算,并分析了算法的误差,最终应用于SAGD生产的实际数据中,通过最终反演结果分析,该方法取得了很好的应用效果.
We present the result of a seismic and gravity joint reservoir density inversion algorithm applied to monitoring the process of steam assisted gravity drainage(SAGD).The method can be used to analyze the density distribution of steam chambers in SAGD in the case of combination of the single well with horizontal wells.Gravimetric density and seismic impedance are combined as a posteriori constraint to formulate the Tikhonov regularization model.Minimization of the Tikhonov regularization model can yield a reasonable inversion result of the density distribution of steam chambers in SAGD.First,a regression model of the seismic impedance and the density through well logging data is built,which can be used to calculate the density constraint data through the seismic P wave data and gives the density constraint.Then,an optimization model based on L2 norm is established.In solving the minimization model,Euler equations along with proper choice of the regularization parameter are simultaneously solved.Utilizing the micro-gravimetric Bouguer gravity anomaly data,the joint inversion scheme based on seismic impedance data and gravimetric data is performed,which yields the density distribution of the reservoir.By analyzing properties of thereservoir density,the density distribution of steam chambers in SAGD can be drawn.Faced with the problem of development and monitoring during SAGD production,this paper investigates the joint inversion method based on the seismic and gravimetric data.The results obtained are as follows:(1)applying the P wave data as a constraint to invert the gravimetric density is proposed;(2)aposteriori regularization method is utilized to calculate the density,in particular,aposteriori choice of the regularization parameter is performed;(3)aquadratic regression model is built for the seismic impedance and the density,theoretical and field data applications reveal that the joint inversion algorithm is stable and can yield high precision of inversion results;(4)it is shown from the field data applications using our joint inversion method that the inverted low density distribution area coincides well with the known well gas injection position,and the connectivity of the low density area is similar to the forward model,hence proves the feasibility of the joint inversion algorithm in practice.We also remark that the joint inversion method relies on the proper initial model and high SNR data;this may be a requirement for practical usage of the joint inversion method.For the development and monitoring problem of SAGD production,we develop ajoint inversion method based on the seismic data and gravimetric data.Using aposteriori choice of the regularization parameter and a quadratic regression model for the seismic impedance and the density,an inversion procedure is performed by solving an Euler equation.The inversion results can well explain and draw the density distribution of steam chambers in SAGD.
引文
Basuyau C,Tiberi C.2011.Imaging lithospheric interfaces and 3Dstructures using receiver functions,gravity,and tomography ina common inversion scheme.Computers&Geosciences,37(9):1381-1390.
    Bulter R M.2004.The behavior of non-condensable gas in SAGD:A rationalization.Journal of Canadian Petroleum Technology,43(1):28-34.
    Coutant O,Bernard M L,Beauducel F,et al.2012.Joint inversionof P-wave velocity and density,application to La Soufrière ofGuadeloupe hydrothermal system.Geophysical Journal International,191(2):723-742.
    Dobróka M,Gyulaiá,Ormos T,et al.1991.Joint inversion ofseismic and geoelectric data recorded in an underground coalmine.Geophys.Prosp.,39(5):643-665.
    Gallardo L A,Meju M A.2004.Joint two-dimensional DC resistivityand seismic travel time inversion with cross-gradients constraints.J.Geophys.Res.,109(B3),doi:10.1029/2003JB002716.
    Hayashi K,Matsuoka T,Hatakeyama H.2005.Joint analysis of asurface-wave method and micro-gravity survey.Journal ofEnvironmental and Engineering Geophysics,10(2):175-184.
    Lines L R,Schultz A K,Treitel S.1988.Cooperative inversion ofgeophysical data.Geophysics,53(1):8-20.
    Liu S Q,Wang X C,Gao Y R,et al.2007.SAGD process with thecombination of vertical and horizontal wells in super-heavy oilreservoir.Petroleum Exploration and Development(in Chinese),34(2):234-238.
    Maceira M,Ammon C J.2009.Joint inversion of surface wavevelocity and gravity observations and its application to centralAsian basins shear velocity structure.J.Geophys.Res.,114(B2),doi:10.1029/2007JB005157.
    Meng W,Jia D,Xie J N,et al.2006.Optimization of geologicaldesign for SAGD process in super heavy oil reservoirs using acombination of vertical and horizontal wells.Journal of DaqingPetroleum Institute(in Chinese),30(2):44-47.
    Moorkamp M,Heincke B,Jegen M,et al.2011.A framework for3-D joint inversion of MT,gravity and seismic refraction data.Geophys.J.Int.,184(1):477-493.
    Onizawa S,Mikada H,Watanabe H,et al.2002.A method forsimultaneous velocity and density inversion and its application toexploration of subsurface structure beneath Izu-Oshimavolcano,Japan.Earth,Planets and Space,54(8):803-817.
    Parsons T,Blakely R,Brocher T.2001.A simple algorithm forsequentially incorporating gravity observations in seismic traveltime tomography.Int.Geol.Rev.,43(12):1073-1086.
    Savino J M,Rodi W L,Masso J F.1980.Simultaneous inversion ofmultiple geophysical data sets for earth structure.45th AnnualInternational Meeting,SEG.
    Sun Y H,Schuster G T.1992.Hierarchic optimizations forsmoothing and cooperative inversion.Annual Meeting Abstracts,Society of Exploration Geophysicists,745-748.
    Tikhotsky S,Achauer U.2008.Inversion of controlled-sourceseismic tomography and gravity data with the self-adaptivewavelet parametrization of velocities and interfaces.Geophys.J.Int.,172(2):619-630.
    Tondi R,De Franco R,Barzaghi R.2000.Sequential integratedinversion of refraction and wide-angle reflection travel times andgravity data for two-dimensional velocity structures.Geophys.J.Int.,141(3):679-698.
    Tondi R,De Franco R,Biella G.2003.Deep and shallow solidearth structures reconstructed with sequential integratedinversion(SII)of seismic and gravity data.Space ScienceReviews,108(1-2):115-130.
    Vozoff K,Jupp D L B.1977.Effective search for a buried layer:Anapproach to experimental design in geophysics.Society ofExploration Geophysicists,8(1):6-15.
    Wang Y F.2007.Computational Methods for Inverse Problems andTheir Applications(in Chinese).Beijing:Higher EducationPress.
    Wang Y F,Stepanova I E,Strakhov V N,et al.2011.InverseProblems in Geophysics and Solution Methods(in Chinese).Beijing:Higher Education Press.
    Wang Y F,Xiao T Y.2001.Fast realization algorithms fordetermining regularization parameters in linear inverseproblems.Inverse Problems,17(2):281-291.
    Yang H.1998.Basement density inversion using gravimetric andseismic data and the integrative interpretation.OGP(inChinese),33(4):496-502.
    Yang H,Dai S K,Song H B,et al.2002.Overview of jointinversion of integrated geophysics.Progress in Geophysics(inChinese),17(2):262-271.
    Yang L Q,Chen Y M,Wang H Y,et al.2007.Physical andnumerical simulation of steam assisted gravity drainage withvertical and horizontal well combination in extra heavy oilreservoir.Journal of China University of Petroleum(in Chinese),31(4):64-69.
    Zeyen H,Pous J.1993.3-D joint inversion of magnetic andgravimetric data with a priori information.Geophysical JournalInternational,112(2):244-256.
    Zhang G B,Shen N H,Wang X C,et al.1993.The programsystem on generalized linear cooperative inversion of thepotential anomalies.Journal of Changchun University ofEarth Sciences(in Chinese),23(2):197-204.
    刘尚奇,王晓春,高永荣等.2007.超稠油油藏直井与水平井组合SAGD技术研究.石油勘探与开发,34(2):234-238.
    孟巍,贾东,谢锦男等.2006.超稠油油藏中直井与水平井组合SAGD技术优化地质设计.大庆石油学院学报,30(2):44-47.
    王彦飞.2007.反演问题的计算方法及其应用.北京:高等教育出版社.
    王彦飞,斯捷潘诺娃I E,提塔连科V N等.2011.地球物理数值反演问题.北京:高等教育出版社.
    杨辉.1998.重力、地震联合反演基岩密度及综合解释.石油地球物理勘探,33(4):496-502.
    杨辉,戴世坤,宋海滨等.2002.综合地球物理联合反演综述.地球物理学进展,17(2):262-271.
    杨立强,陈月明,王宏远等.2007.超稠油直井-水平井组合蒸汽辅助重力泄油物理和数值模拟.中国石油大学学报(自然科学版),31(4):64-69.
    张贵宾,申宁华,王喜臣等.1993.位场广义线性综合反演系统的建立.长春地质学院学报,23(2):197-204.

版权所有:© 2023 中国地质图书馆 中国地质调查局地学文献中心