2. 油藏产能方程基础
要绘制IPR曲线，先要弄清曲线背后的物理方程。IPR曲线描述的是“井底流压变化时，油井能提供多少产量”，而这一关系的根源，是流体在多孔介质中的渗流规律。本节从达西定律出发，推导平面径向流的产能公式，引出采油指数的概念，并明确供给边界压力与井底流压之间的驱动关系。
一、达西渗流定律：产能方程的物理起点
达西（Darcy）通过砂柱实验发现，流体通过多孔介质的流量与两端压差成正比、与渗流路径长度成反比，比例系数由介质渗透率和流体黏度决定。其微分形式为：
v = -(K/μ) · dp/dr式中：v为渗流速度（m/s）；K为有效渗透率（m²，工程常用μm²或mD）；μ为流体黏度（Pa·s）；dp/dr为压力梯度（Pa/m）。负号表示流体沿压力降低方向流动，这是本节统一采用的符号约定——压力梯度取流动方向为正，速度方向与压力降落方向一致。
达西定律有三个关键前提：层流、单相、稳定渗流。油藏工程中的IPR分析通常建立在这三个假设之上，因此得到的是一条理想化的直线或曲线；当出现湍流、多相流或非达西效应时，实际曲线会偏离理论形态，这正是后续章节要修正的内容。
把微分形式的达西定律与连续性方程结合，并引入岩石和流体的压缩性，可得到渗流的控制方程。对于均质、等厚、圆形封闭或定压边界的油藏，在稳定渗流条件下，控制方程简化为径向的拉普拉斯方程，其解即为平面径向流产能公式。
二、平面径向流产能公式：从微分方程到产量表达式
考虑一个半径为re的圆形供给边界，边界压力为pe；井半径为rw，井底流压为pwf。流体从供给边界沿径向流入井底，渗流截面积为2πrh。将达西定律写成流量形式并对半径积分：
q = 2πKh(p_e - p_wf) / [μB·ln(r_e/r_w)]式中：q为地面产量（m³/s）；h为油层有效厚度（m）；B为原油体积系数（地下体积/地面体积）；re/rw为供给半径与井半径之比，工程上常取对数项ln(re/rw)≈7～8。若采用工程单位（q以m³/d、K以μm²、h以m、p以MPa、μ以mPa·s），公式可写为：
q = 0.0864 · 2πKh(p_e - p_wf) / [μB·ln(r_e/r_w)]这个公式揭示了三个重要事实。第一，产量与生产压差(pe−pwf)成正比，压差越大产量越高，这是IPR曲线横纵坐标关系的直接来源。第二，产量与Kh/μ成正比，即地层系数越大、流体越稀，同样的压差能采出更多油，这决定了IPR曲线的斜率。第三，产量与ln(re/rw)成反比，说明井径的影响是对数级的——把井径扩大一倍，产量增幅有限，这也解释了为什么靠扩眼提产往往不经济。
需要特别强调：公式中的pe是供给边界压力，代表油藏能量；pwf是井底流压，代表井筒对地层的“回压”。当pwf=pe时，压差为零，产量为零，对应IPR曲线与纵轴的交点；当pwf降低，产量沿曲线上升；理论上pwf=0时产量达到最大值，即最大潜在产量qmax。实际生产中pwf不可能降到零，因此qmax是外推得到的参考值。
三、采油指数：把产能“浓缩”成一个系数
把平面径向流公式改写为q = J·(pe−pwf)，其中：
J = 2πKh / [μB·ln(r_e/r_w)]J即采油指数（m³/(d·MPa)），表示单位生产压差下的日产油量。它把渗透率、厚度、黏度、体积系数、供给半径和井半径等参数“打包”成一个系数，是矿场上最常用的产能指标。J的物理意义很直观：J越大，IPR曲线越陡，同样的流压下降能换来更多产量；J越小，曲线越平缓，说明地层供液能力弱。
采油指数有两个实用价值。其一，快速估算产量：已知J和当前流压，可直接算出稳定产量，无需每次重算全套参数。其二，反演地层信息：通过测试获得q与pwf数据，回归出J，再结合已知的h、μ、B等参数，可反算有效渗透率K，用于判断储层改造效果或污染程度。例如某井测试得到J=12 m³/(d·MPa)，若理论计算值应为20 m³/(d·MPa)，则说明存在近井污染或完井不完善，需要酸化解堵。
四、供给边界压力与井底流压的关系
供给边界压力pe与井底流压pwf共同构成生产压差，二者关系可从三个层面理解。
能量层面：pe是油藏能量的“源头”，由边底水、溶解气膨胀或人工注水维持；pwf是井筒对地层的“出口回压”。pe−pwf就是驱动流体从边界流向井底的压差。pe下降意味着油藏能量衰竭，整条IPR曲线会向左下方移动；pwf下降则是主动放大压差、提高产量的操作手段。
约束层面：pwf不能无限降低。一方面受井筒液柱和泵吸入口压力限制，另一方面过低流压会引发脱气、出砂、水锥等问题。因此实际可用的生产压差存在下限约束，IPR曲线上真正可用的工作段往往只是曲线的一部分。
动态层面：随着开采进行，pe逐渐下降，而pwf由生产制度决定。若保持pwf不变，压差会随pe下降而减小，产量自然递减；若想稳产，就必须同步下调pwf。这一动态关系是后续“IPR与VLP联合分析”和“合理生产制度设计”的基础。
上图左侧为平面径向流的物理图像：压力从供给边界pe沿径向向井底pwf逐渐降落，等压线呈同心圆；右侧为对应的稳定渗流IPR直线，纵轴为井底流压、横轴为产量，直线斜率等于1/J。工作点落在直线上，其横坐标就是当前流压下的产量。需要说明的是，这条直线是单相稳定渗流的理想结果，实际油藏在低于饱和压力后会出现弯曲，那将是后续Vogel方程要处理的问题。
五、一个教学算例
假设某均质油藏：K=0.15 μm²，h=12 m，μ=2.0 mPa·s，B=1.15，re=300 m，rw=0.1 m，pe=28 MPa。先算对数项ln(300/0.1)=ln3000≈8.01。代入采油指数公式（工程单位）：
J = 2π × 0.15 × 12 / (2.0 × 1.15 × 8.01) ≈ 0.614 m³/(d·MPa)若当前井底流压pwf=22 MPa，则生产压差为6 MPa，稳定产量q=0.614×6≈3.7 m³/d。若把流压降到18 MPa，压差变为10 MPa，产量升至约6.1 m³/d。可见在直线段内，产量随流压下降呈线性增长。若继续降到10 MPa，理论产量约11.1 m³/d，但此时需检查是否低于饱和压力——若饱和压力为15 MPa，则10 MPa已进入两相流区，直线预测会偏高，必须改用Vogel关系修正。这一判断正是IPR绘制中“先判流态、再选方程”的关键步骤。
归纳本节要点：达西定律给出渗流速度与压力梯度的正比关系；平面径向流积分得到q与生产压差的线性产能公式；采油指数J把储层与流体参数浓缩为一个斜率系数；供给边界压力pe代表能量来源，井底流压pwf代表生产回压，二者之差是产量的直接驱动。掌握这四层关系，就具备了绘制单相稳定渗流IPR直线的完整基础。