油藏数值模拟是用数值方法模拟和再现油藏内流体渗流现象和渗流规律的一门学科,目前它已成为研究与制定开发方案的重要工具,同时也为提高采收率等石油工程问题提供一种可靠的科学依据。油藏数值模拟融合了油藏工程、渗流力学、计算数学、计算机等多学科知识,其目的在于求解油藏内压力和各种流体的饱和度等多个物理量在油藏中的分布。在实际油藏开发过程中,大多需要揭示油水在油藏中的运移规律,而描述这种运移规律通常运用一些物理定律建立一组数学模型,其中应用最为广泛的是油水两相数学模型。油水两相数学模型实质上是由一组非线性抛物型偏微分方程、辅助方程、边界条件和初始条件等构成的定解问题。油水两相数学模型具有较强的非线性,这一性质极大地增加了求解难度,从稳定性角度考虑,在众多解决方法中,全隐式方法的稳定性最高,因此本文主要研究油水两相全隐式数学模型的求解方法。针对该问题,本文选用了块中心网格和角点网格两种方法对油藏进行网格剖分,采用全隐式方法对流动方程进行数值离散,由此得到非线性代数方程组,并利用Newton-Raphson方法将非线性代数方程线性化后迭代求解。在求解过程中,本文着重解决了网格传导系数计算问题、油藏相关未知量初始化问题、井筒中流体平均密度计算等问题。最后针对5个油藏生产案例,先采用本文算法进行数值计算,然后将其数值结果与Eclipse软件的计算结果进行对比,从而验证了本文算法的准确性与可行性。