孔隙度数值反演计算

计算物理 ›› 2010, Vol. 27 ›› Issue (3): 413-422.

孔隙度数值反演计算

刘福平^1,2, 曹跃祖¹, 王安玲¹, 杨长春²

1. 北京印刷学院, 北京 102600;
2. 中国科学院地质与地球物理研究所, 北京 100029

收稿日期:2009-03-23 修回日期:2009-07-07 出版日期:2010-05-25 发布日期:2010-05-25
作者简介:刘福平(1960-),男,山东莒县,教授,博士,主要从事油藏数值模拟、电磁测井理论及数值计算方面的研究和教学工作,北京印刷学院等离子体研究室102600.
基金资助:
北京市自然基金(1083011);北京市属市管高等学校人才强教计划资助项目

Numerical Inversion of Porosity

LIU Fuping^1,2, CAO Yuezu¹, WANG Anling¹, YANG Changchun²

1. Beijing Institute of Graphic Communication, Beijing 102600, China;
2. Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China

Received:2009-03-23 Revised:2009-07-07 Online:2010-05-25 Published:2010-05-25

摘要/Abstract

摘要： 由渗流微分方程定解问题和Peaceman方程给出了网格压力、井底压力对网格孔隙度的导数,利用三维渗流方程压强数值解计算井底压力对网格孔隙度的导数;采用共轭梯度法实现孔隙度均匀(或分块均匀)分布油藏模型的反演计算.算例表明,经过8~10次迭代后反演结果与真值的最大相对误差在0.03%以内,反演收敛于真值.

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

Abstract: 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.

Key words: porosity inversion, numerical calculation, optimization, sensitivity coefficients, convergent property

中图分类号:

P631

刘福平, 曹跃祖, 王安玲, 杨长春. 孔隙度数值反演计算[J]. 计算物理, 2010, 27(3): 413-422.

LIU Fuping, CAO Yuezu, WANG Anling, YANG Changchun. Numerical Inversion of Porosity[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2010, 27(3): 413-422.

[1]	陈上, 崔国民, 张春伟, 段欢欢. 一种热集成系统的动态多智能体微分进化算法[J]. 计算物理, 2016, 33(3): 349-357.
[2]	薛社生, 李守先, 刘阳. 一种气体低温吸附模型与计算[J]. 计算物理, 2014, 31(4): 412-416.
[3]	刘福平, 王安玲, 杨长春. 地层渗透参数反演计算[J]. 计算物理, 2013, 30(6): 879-885.
[4]	曾凡海, 李常品. 时间分数阶亚扩散方程的高阶差分方法[J]. 计算物理, 2013, 30(4): 491-500.
[5]	游荣义, 黄晓菁. 纳米结构金属表面的近场电势分布[J]. 计算物理, 2013, 30(4): 582-586.
[6]	刘全, 王瑞利, 林忠, 刘希强. 流体力学拉氏程序收敛性及数值计算不确定度初探[J]. 计算物理, 2013, 30(3): 346-352.
[7]	王海, 王海鹰, 周伟国, 周海珠. 供热管网中管段阻力系数的辨识方法[J]. 计算物理, 2013, 30(3): 422-432.
[8]	程钰锋, 聂万胜. 电弧放电等离子体诱导激波的计算[J]. 计算物理, 2012, 29(2): 213-220.
[9]	贺宏, 李会雄, 冯永昌. 超燃冲压发动机等截面隔离段中流动的LBM模拟[J]. 计算物理, 2012, 29(2): 234-238.
[10]	吴磊, 刘全金, 章礼华, 王其申. 外伸梁差分离散系统的模态反问题[J]. 计算物理, 2010, 27(3): 407-412.
[11]	王廷春, 张鲁明, 陈芳启. 求解耦合SchrÖdinger方程组的数值逼近算法[J]. 计算物理, 2008, 25(2): 179-185.
[12]	刘福春, 金明星, 高欣, 路秀秀, 丁大军. 静电六极装置电场分布的数值计算方法及在N₂O分子转动态选择、取向中的应用[J]. 计算物理, 2008, 25(2): 230-234.
[13]	张毅锋, 邓小刚, 毛枚良, 陈坚强. 一种可压缩流动的高阶加权紧致非线性格式(WCNS)的加速收敛方法[J]. 计算物理, 2007, 24(6): 698-704.
[14]	殷海荣, 宫玉彬, 魏彦玉, 黄民智, 路志刚, 王文祥. 虚边界元法与强流电子枪的电子轨迹模拟[J]. 计算物理, 2006, 23(6): 647-654.
[15]	王菊芬, 蒲家宁, 孟浩龙. 输油管道油流静电带电数值计算[J]. 计算物理, 2006, 23(5): 614-620.