4. 双重介质渗流数学模型建立

双重介质模型的建立，本质上是把“裂缝快速导流、基质缓慢供液”这一物理图像翻译成一组可求解的偏微分方程。它需要四类要素：每个系统的运动方程（流速与压力梯度的关系）、状态方程（流体与岩石物性随压力如何变化）、由质量守恒导出的渗流微分方程，以及使方程组有唯一解的定解条件。本节按这四类要素依次展开，并给出一个可直接对照的方程结构图。

一、运动方程：两个系统各自的达西定律

设基质系统压力为 p_m，裂缝系统压力为 p_f，取 z 轴向上为正、重力加速度为 g。两个系统各自满足达西定律：

基质：  v_m = -(K_m/μ) · (∇p_m - ρg∇z)
裂缝：  v_f = -(K_f/μ) · (∇p_f - ρg∇z)

其中 K_m 为基质渗透率（量级常在 0.01～1 mD），K_f 为裂缝系统等效渗透率（可达几十至上千 mD），μ 为流体黏度。两者相差 2～4 个数量级，正是“裂缝导流、基质供液”的数学体现。注意两个系统共用同一个流体黏度，但渗透率张量可以不同；若裂缝有优势方位，K_f 应写成张量形式，简化为各向同性时退化为标量。

二、状态方程：流体压缩性与岩石压缩性

对微可压缩流体，密度随压力线性化：ρ = ρ₀[1 + C_L(p - p₀)]，C_L 为流体压缩系数。岩石方面，双重介质有两个不同的压缩性：基质孔隙压缩系数 C_m 与裂缝压缩系数 C_f。裂缝性储层中 C_f 往往比 C_m 大一个量级以上，因为裂缝的闭合/张开比基质孔隙更敏感。综合储集系数常写为：

基质系统：  C_tm = C_L + C_m
裂缝系统：  C_tf = C_L + C_f

这两个系数分别进入两个系统的累积项，决定了压力波在各自系统中的传播速度与储量释放能力。

三、窜流项：两个系统之间的耦合

双重介质模型的关键耦合是基质向裂缝的窜流。最常用的是拟稳态窜流（Warren-Root 形式）：

q_ex = α · (K_m/μ) · (p_m - p_f)

α 为形状因子，量纲为 1/L²，对规则正交裂缝网络，α = 4n(n+2)/L²（n 为裂缝法向维数，L 为基质块特征尺寸）。当 p_m > p_f 时窜流为正，即基质向裂缝供液；反之裂缝向基质倒灌。这一项是双重介质区别于单重介质的核心，也是后面所有算例中“压力导数出现下凹”的物理来源。

四、完整渗流微分方程组

将运动方程代入质量守恒方程，并利用状态方程线性化，得到两个耦合的渗流微分方程（忽略重力项、以体积系数形式表达）：

基质：  ∇·[(K_m/μ)∇p_m] - α(K_m/μ)(p_m - p_f) = (φ_m C_tm) ∂p_m/∂t
裂缝：  ∇·[(K_f/μ)∇p_f] + α(K_m/μ)(p_m - p_f) = (φ_f C_tf) ∂p_f/∂t

式中 φ_m、φ_f 分别为基质与裂缝的孔隙度。两个方程通过窜流项 α(K_m/μ)(p_m−p_f) 双向耦合：基质方程中它是汇项，裂缝方程中它是源项，符号相反、量值相等，保证总质量守恒。若考虑井筒生产，还需在裂缝方程右端加入源汇项 q_w（生产井为负、注入井为正）。

为便于无量纲分析，常引入窜流系数 λ = α·(K_m/K_f)·r_w² 和储容比 ω = φ_fC_tf/(φ_fC_tf+φ_mC_tm)。λ 控制窜流发生的早晚，ω 控制裂缝储量的占比，这两个参数是双重介质试井解释的核心，也是后续算例中要反复调整的对象。

五、定解条件：初始条件与边界条件

上述两个方程是二阶偏微分方程，必须有初始条件和边界条件才能定解。初始条件通常取原始地层压力：p_m(x,y,z,0) = p_f(x,y,z,0) = p_i。若储层经过长期开采，也可用当前实测压力分布作为初始场。

边界条件分内边界与外边界。内边界描述井筒：定产时给定井底流量 q_w（裂缝系统承担主要产能），定压时给定井底流压 p_wf。外边界有三种常见类型：无穷大边界（p → p_i 当 r → ∞）、封闭边界（∂p/∂n = 0）、定压边界（p = p_i）。实际建模时，若边界距离大于 3～5 倍探测半径，可按无穷大处理；若存在断层遮挡，则用封闭边界。

六、一个可对照的教学示例

假设某裂缝性碳酸盐岩油藏：p_i = 30 MPa，K_m = 0.5 mD，K_f = 200 mD，φ_m = 0.18，φ_f = 0.004，C_L = 8×10⁻⁴ MPa⁻¹，C_m = 5×10⁻⁴ MPa⁻¹，C_f = 2×10⁻³ MPa⁻¹，基质块尺寸 L = 2 m，n = 3。则 α = 4×3×5/4 = 15 m⁻²，λ = 15×(0.5/200)×0.1² ≈ 3.75×10⁻⁴，ω = 0.004×2.8×10⁻³/(0.004×2.8×10⁻³ + 0.18×1.3×10⁻³) ≈ 0.045。这组参数意味着裂缝储量仅占约 4.5%，窜流发生较晚，压力导数曲线会出现明显的“下凹—平台”双重介质特征。上述数值为教学假设，用于说明参数如何进入方程，不代表任何实际油田的实测结果。