基于SN算法的复合油藏模型求解

针对井筒储集和表皮效应的复合油藏模型的解析解形式复杂、常规数值求解过程中易出现曲线震荡问题,开展了对模型数值求解精度的研究。对拉氏空间的井底无因次压力的解运用Stehfest数值反演,得到对应的实空间解,分析解的形式可知,虚宗量Bessel函数的求解是整个复合油藏模型数值计算的关键,给出了虚宗量Bessel函数的分段数值(SN)算法,其计算精度可与Matlab计算结果的精度比拟,在此基础上得到的无穷大、定压、封闭3种外边界条件下的的曲线平滑完整,并且程序运行速度远远满足试井分析快速响应的实际需要,可以很好地适用于复合油藏的试井分析。