2. 双重介质渗流的基本假设与连续性方程
裂缝性油藏与常规砂岩油藏最本质的差别，在于储集空间与渗流通道被“分工”到了两套介质上：基质岩块孔隙负责储油，但渗透率极低；裂缝网络渗透率很高，却只占很小体积。要把这种“两套介质、两种角色”的物理图像转化为可求解的数学模型，就必须先建立一组明确的基本假设，再在代表性单元体上写出质量守恒方程，最后用窜流项把两个系统耦合起来。本节按“假设—REV—守恒方程—窜流项”的顺序展开，并给出可直接用于编程的方程形式。
一、双重介质的基本假设
Warren–Root 模型是双重介质渗流最经典的理论框架，其基本假设可归纳为以下几条，每一条都直接决定后续方程的形态：
- 两套连续介质并存：基质系统和裂缝系统各自充满整个油藏空间，二者在空间上处处重叠，而不是“基质块之间夹着裂缝”的离散拼图。这是把裂缝“连续化”的关键一步。
- 基质岩块为均质、各向同性：每个岩块内部孔隙度 φm、渗透率 Km 视为常数，岩块形状用规则几何体（正交六面体或球体）近似。
- 裂缝系统为规则正交网络：由三组互相垂直、间距相等的裂缝面构成，裂缝开度 b、间距 a 为统计平均值，裂缝渗透率 Kf 远大于 Km。
- 基质向裂缝的流动为拟稳态窜流：岩块内部压力梯度被忽略，窜流量只与基质—裂缝之间的压差成正比，比例系数即形状因子 σ。
- 流体微可压缩、渗流服从达西定律：忽略毛管力和重力影响，等温渗流。
- 裂缝系统直接与井筒连通：油井只从裂缝产油，基质不直接向井供液，而是先窜流到裂缝再流入井底。
二、代表性单元体（REV）与双重孔隙度
连续性假设要求存在一个尺度：在该尺度上取体积平均，孔隙度、渗透率等参数趋于稳定，这个体积就是代表性单元体（REV）。对双重介质，REV 的选取有一个额外约束——它必须同时包含足够多的基质岩块和足够多的裂缝，即 REV 尺度要远大于裂缝间距 a，又要远小于整个油藏尺度。若取体积 V，则：
裂缝孔隙度 φf = Vf/V，基质孔隙度 φm = Vm/V，总孔隙度 φ = φm + φf。由于裂缝体积占比通常只有 0.1%～1%，故 φ ≈ φm，但裂缝渗透率可高出基质 2～4 个数量级，这正是“储渗分离”的定量体现。
对正交裂缝网络，若三组裂缝间距均为 a、开度均为 b，则裂缝孔隙度近似为 φf ≈ 3b/a；裂缝渗透率按平板流立方定律估算为 Kf ≈ b³/(12a)。这两个式子把几何参数（a、b）与物性参数（φf、Kf）直接挂钩，是后续取值的依据。
三、质量守恒方程的推导
在 REV 上对裂缝系统写质量守恒：单位时间流入与流出 REV 的流体质量差，加上基质窜流进入的质量，等于 REV 内流体质量的变化率。设流体密度 ρ、黏度 μ，裂缝压力 pf，则：
裂缝系统：
∂(φf·ρ/μ)/∂t = ∇·[ (Kf·ρ/μ)·∇pf ] + q + q_well
基质系统：
∂(φm·ρ/μ)/∂t = ∇·[ (Km·ρ/μ)·∇pm ] − q其中 q 为单位体积基质向裂缝的窜流量（质量/体积/时间），q_well 为井筒采出项（仅裂缝系统有）。引入微可压缩假设 ρ = ρ₀ec(p−p₀)，并定义综合压缩系数，可将上式线性化为常见的双重介质扩散方程形式：
裂缝：φf·ct·∂pf/∂t = ∇·(Kf/μ·∇pf) + α·(pm − pf)
基质：φm·ct·∂pm/∂t = ∇·(Km/μ·∇pm) − α·(pm − pf)式中 α 为窜流系数，α = σ·Km/μ，σ 为形状因子。注意符号约定：本文取“基质压力高于裂缝压力时 q 为正、由基质流向裂缝”，因此基质方程中窜流项为负、裂缝方程中为正，全文与图示保持一致。
四、窜流项的引入与形状因子
窜流项是双重介质模型区别于单重介质模型的灵魂。Warren–Root 给出的拟稳态窜流为：
q = σ·(Km/μ)·(pm − pf)，其中形状因子 σ 对正交裂缝网络取 σ = 4n(n+2)/a²，n 为裂缝组数。当 n = 3 时，σ = 60/a²；对球状岩块则常用 σ = 15/r²（r 为岩块半径）。
以一个教学示例说明量级：取裂缝间距 a = 1 m、基质渗透率 Km = 0.1 mD、原油黏度 μ = 2 mPa·s，则 σ = 60/1² = 60 m⁻²，窜流系数 α = σKm/μ ≈ 60 × 1×10⁻¹⁶ / (2×10⁻³) ≈ 3×10⁻¹² m/(Pa·s) 量级。若基质压力比裂缝高 1 MPa，则单位体积窜流强度约为 α·Δp ≈ 3×10⁻⁶ s⁻¹ 量级的压力响应速率——这正是“基质供液慢、裂缝压力恢复快”的数学根源。上述数值为示例假设，实际取值须以岩心分析与裂缝描述为准。