二维注热开采天然气水合物模型控制方程

📌 1. 模型假设

为保证模型可解与稳定，采用如下简化：

• 单相等效流体（忽略气/水速度差异）
• 各相局部热平衡（LTE）
• 多孔介质各向同性
• 水合物分解为动力学控制
• 忽略毛细压力与重力

🧩 2. 控制方程

2.1 质量守恒（Darcy 流）

连续性方程

```math id=”eq_mass” \frac{\partial (\phi \rho)}{\partial t} + \nabla \cdot (\rho \mathbf{u}) = m_g

---

### Darcy 定律
```math id="eq_darcy"
\mathbf{u} = -\frac{k}{\mu} \nabla p

渗透率演化

```math id=”eq_perm” k = k_0 \left( \frac{\phi}{\phi_0} \right)^3

---

## 2.2 能量守恒（多孔介质）

```math id="eq_energy"
(\rho C_p)_{eff} \frac{\partial T}{\partial t}
+ \rho C_p \mathbf{u} \cdot \nabla T
= \nabla \cdot (k_{eff} \nabla T) + Q_h

等效体积热容

```math id=”eq_ceq” (\rho C_p)_{eff} = (1-\phi)\rho_s C_s + \phi \rho_f C_f

---

### 等效导热系数

```math id="eq_keq"
k_{eff} = (1-\phi)\lambda_s + \phi \lambda_f

2.3 水合物分解动力学

```math id=”eq_kinetic” \frac{dS_h}{dt} = -k_{react} \cdot \max\left( \frac{p - p_{eq}}{p_{eq}}, 0 \right) S_h

---

## 2.4 分解热源项
```math id="eq_Qh"
Q_h = \rho_h \cdot k_{react} \cdot \max\left( \frac{p - p_{eq}}{p_{eq}}, 0 \right) S_h \cdot \Delta H

🔍 3. 方程推导要点

3.1 能量方程推导

从多相体系能量守恒： ```math id=”eq_energy_origin” \sum_i \frac{\partial (\phi \rho_i S_i C_i T)}{\partial t}

• \sum_i \nabla \cdot (\rho_i C_i \mathbf{u}i T) = \nabla \cdot (k{eff} \nabla T) + Q_h ``` —

在假设：

• 单相等效流速：(\mathbf{u}_i \approx \mathbf{u})
• 局部热平衡

可化简为：

```math id=”eq_energy_simplified” (\rho C_p)_{eff} \frac{\partial T}{\partial t}

• \rho C_p \mathbf{u} \cdot \nabla T = \nabla \cdot (k_{eff} \nabla T) + Q_h ``` —

3.2 热源项来源

水合物分解：

```math id=”eq_reaction” Hydrate \rightarrow Gas + Water

单位体积反应速率：

```math id="eq_rate"
r = \rho_h \cdot \frac{dS_h}{dt}

乘以反应焓：

```math id=”eq_Qh_derivation” Q_h = r \cdot \Delta H

---

# ⚙️ 4. COMSOL 对应关系

| 物理方程     | COMSOL接口                      |
| -------- | ----------------------------- |
| Darcy 定律 | Darcy's Law                   |
| 能量方程     | Heat Transfer in Porous Media |
| 分解动力学    | Domain ODE                    |
| 热源项      | Heat Source                   |

---

### 关键耦合

```comsol
Velocity field = [dl.u, dl.v]
Q = Qh

📊 5. 初始与边界条件

初始条件

```math id=”eq_init” T = T_0,\quad p = p_0,\quad S_h = S_{h0}

---

## 边界条件

### 注热边界

```math id="eq_bc_T"
T = T_{inj}

生产井

```math id=”eq_bc_p” p = p_{out}

---

### 其余边界

```math id="eq_bc_other"
\mathbf{n}\cdot \nabla T = 0,\quad \mathbf{n}\cdot \mathbf{u} = 0

🚀 6. 模型扩展方向

6.1 双相流模型

```math id=”eq_two_phase” \mathbf{u}_g \neq \mathbf{u}_w

---

## 6.2 孔隙率演化

```math id="eq_phi"
\phi = \phi_0 + \alpha (S_{h0} - S_h)

6.3 THM耦合

• 热（T）
• 渗流（p）
• 力学（σ）
• 分解（Sh）

✅ 7. 总结

该模型构建了一个注热驱动水合物分解的最小物理闭环：

text id="eq_flow" 注热 → 升温 → 分解 → 孔隙变化 → 渗流 → 继续传热