考虑孔隙流体的质量守恒,流体连续性方程描述了流体在多孔介质中的流动情况。以下是详细推导过程:
流体连续性方程
在多孔介质中,流体的质量守恒方程可以表示为:
∂
(
ϕ
ρ
f
)
∂
t
+
∇
⋅
(
ρ
f
v
f
)
=
0
\frac{\partial (\phi \rho_f)}{\partial t} + \nabla \cdot (\rho_f \mathbf{v_f}) = 0
∂t∂(ϕρf)+∇⋅(ρfvf)=0
其中:
- ϕ \phi ϕ 是孔隙度。
- ρ f \rho_f ρf是流体密度。
- v f \mathbf{v_f} vf 是流体的平均速度。
为了简化计算,假设流体密度 (\rho_f) 是常数,即流体是不可压缩的。则方程简化为:
ϕ
∂
ρ
f
∂
t
+
ρ
f
∂
ϕ
∂
t
+
ρ
f
∇
⋅
v
f
=
0
\phi \frac{\partial \rho_f}{\partial t} + \rho_f \frac{\partial \phi}{\partial t} + \rho_f \nabla \cdot \mathbf{v_f} = 0
ϕ∂t∂ρf+ρf∂t∂ϕ+ρf∇⋅vf=0
ρ
f
∂
ϕ
∂
t
+
ρ
f
∇
⋅
v
f
=
0
\rho_f \frac{\partial \phi}{\partial t} + \rho_f \nabla \cdot \mathbf{v_f} = 0
ρf∂t∂ϕ+ρf∇⋅vf=0
将流体密度 (\rho_f) 提取出来,得:
∂
ϕ
∂
t
+
∇
⋅
v
f
=
0
\frac{\partial \phi}{\partial t} + \nabla \cdot \mathbf{v_f} = 0
∂t∂ϕ+∇⋅vf=0
达西定律
达西定律描述了流体在多孔介质中的流动:
v
f
=
−
k
μ
∇
p
\mathbf{v_f} = -\frac{k}{\mu} \nabla p
vf=−μk∇p
其中:
- k k k是渗透率。
- μ \mu μ 是流体的黏度。
- p p p 是孔隙压力。
将达西定律代入流体连续性方程中:
∂
ϕ
∂
t
−
∇
⋅
(
k
μ
∇
p
)
=
0
\frac{\partial \phi}{\partial t} - \nabla \cdot \left( \frac{k}{\mu} \nabla p \right) = 0
∂t∂ϕ−∇⋅(μk∇p)=0