姜瑞忠，倪庆东，张春光，张凤博，郜益华，蔺靖淇

（1.中国石油大学（华东）石油工程学院，山东青岛 266580；2.中海石油（中国）有限公司上海分公司，上海 200335；3.中海油研究总院有限责任公司，北京 100028）

稠油油藏储量在未动用油藏中占比较大，具有较高的开发潜能［1-2］。基于稠油油藏黏度高、渗流阻力大的特点［3］，考虑变启动压力梯度的数值模拟研究用于稠油油藏开发［4］。张跃雷通过测量不同黏度稠油的流变性，对启动压力梯度与塑性黏度之间的关系做出了进一步的确定［5］。王泊等分别用数学模型以及实验的方法对非线性渗流启动压力梯度进行研究［6-7］。张芹月提出启动压力梯度与流度之间存在幂律关系［8］。张代燕等通过对实验数据进行回归，得到了有关启动压力梯度的经验公式［9］，并且提出稠油流动的非线性程度可用启动压力梯度来表征［10］。刘英宪等基于压敏效应得出启动压力梯度与生产井产量之间的关系［11-13］，同时严利咏等也提出启动压力梯度和压力敏感系数对地层压力的影响较为显著［14-15］。通过对特低渗透油藏的研究，李承龙等考虑压敏效应与启动压力梯度会对生产指标有显著影响［16-19］。目前稠油油藏非线性渗流数值模拟方法在模型求解过程中，未同时考虑储层的非均质性以及压敏效应和启动压力梯度变化［20-22］，一些简单的模型难以解决现场复杂区块的生产状况。

为此，笔者基于储层的非均质性以及稠油油藏在多孔介质中表现出非达西渗流的特点，应用启动压力梯度表征的油相视黏度与流动修正系数，描述稠油油藏黏度异常的现象，将压力梯度处理成与流度相关的变量并考虑其方向性，同时在考虑压敏效应的基础上建立全隐式数值模型。……