孔隙度数值反演计算  

Numerical Inversion of Porosity

作　　者：刘福平[1,2] 曹跃祖[1] 王安玲[1] 杨长春[2] 

机构地区：[1]北京印刷学院,北京102600 [2]中国科学院地质与地球物理研究所,北京100029

出　　处：《计算物理》2010年第3期413-422,共10页Chinese Journal of Computational Physics

基　　金：北京市自然基金(1083011);北京市属市管高等学校人才强教计划资助项目

摘　　要：由渗流微分方程定解问题和Peaceman方程给出了网格压力、井底压力对网格孔隙度的导数,利用三维渗流方程压强数值解计算井底压力对网格孔隙度的导数;采用共轭梯度法实现孔隙度均匀(或分块均匀)分布油藏模型的反演计算.算例表明,经过8～10次迭代后反演结果与真值的最大相对误差在0.03%以内,反演收敛于真值.Derivatives of well-bore and grid-block pressure to grid porosity are derived from initial boundary value problems of fluid equations in porous media and inner boundary condition.The derivatives are computed with numerical solution of 3D fluid equation in porous media.A inverse iterative equation is constructed by conjugate method.Numerical inversion of porosity distribution are realized for uniform or block uniform distribution reservoirs.It shows that the relative error of inverse porosity is less than three ten thousandths for a homogenous porosity reservoir with about 8-10 iterations.Relative error of inverse porosity is less than one thousandth for a block homogenous porosity reservoir.Computed porosity automatically approachs to the real value in a few iterations.

关 键 词：孔隙度反演 数值计算 最优化 敏感系数 收敛性 

分 类 号：P631[天文地球—地质矿产勘探]