基于sn算法的复合油藏模型求解-西南石油大学学报.pdf

西南石油大学学报（自然科学版） 2016 年6 月第38 卷第3 期 Journal of Southwest Petroleum University（Science TechnologyEdition） Vol.38 No. 3 Jun. 2016 DOI ：10.11885/j.issn.1674 5086.209 文章编号：1674 5086（2016）03 0114 07 中图分类号：TE331.1 文献标志码：A 基于SN 算法的复合油藏模型求解 * 李铁军 ，仝立媛，郭大立，崔林鹏，杨武剑 西南石油大学理学院，四川成都 610500 摘 要：针对井筒储集和表皮效应的复合油藏模型的解析解形式复杂、常规数值求解过程中易出现曲线震荡问题，开 展了对模型数值求解精度的研究。对拉氏空间的井底无因次压力的解运用Stehfest 数值反演，得到对应的实空间解， 分析解的形式可知，虚宗量Bessel 函数的求解是整个复合油藏模型数值计算的关键，给出了虚宗量Bessel 函数的分段 数值（SN）算法，其计算精度可与Matlab 计算结果的精度比拟，在此基础上得到的无穷大、定压、封闭3 种外边界条件 下的的曲线平滑完整，并且程序运行速度远远满足试井分析快速响应的实际需要，可以很好地适用于复合油藏的试井 分析。 关键词：曲线震荡；Stehfest 数值反演；分段数值算法；虚宗量Bessel 函数；平滑完整 TheSolutionofCompositeReservoirModelBasedonBesselFunctionof ImaginaryArgumentSNAlgorithm * LI Tiejun ,TONGLiyuan,GUODali,CUILinpeng,YANGWujian Facultyof Science，SouthwestPetroleumUniversity，Chengdu，Sichuan610500，China Abstract ：Analytical solution of composite reservoir model which considers wellbore storage and skin effect is complex in its form and curve shocks occur in the conventional numerical solving process. Therefore，it is necessary to make researches on accuracy of numerical solution. Firstly，the real space solution is obtained by using the stehfest numerical inversion to the analyticalsolutionofLaplacespacedimensionlesspressure. Afterthat，byanalyzingtheformofthesolution ，wedeterminethe solving of the Bessel function of imaginary argument is the key to numerical solution of the model. Then we present sectional numerical（SN）algorithmofBesse