[论文解读] Upscaling of unsaturated flow in fractured porous media

Florian List, Kundan Kumar+2 more|arXiv (Cornell University)|Jul 16, 2018

Advanced Mathematical Modeling in Engineering参考文献 48被引用 5

一句话总结

本文通过取薄裂缝宽度 ε → 0 的极限，严格推导了裂缝性孔隙介质中非饱和流的有效模型，表明所得到的尺度化模型取决于孔隙率和渗透率比值随 ε 的缩放方式。该分析根据指数 κ 和 λ 的不同，导出了不同的有效方程，例如当 λ < 1 且 κ > -1 时，裂缝界面处为一维 Richards 方程。数值验证显示，当 ε ≤ 10⁻³ 时，误差低于 10⁻⁵。

ABSTRACT

In this work, we consider a mathematical model for flow in a unsaturated porous medium containing a fracture. In all subdomains (the fracture and the adjacent matrix blocks) the flow is governed by Richards' equation. The submodels are coupled by physical transmission conditions expressing the continuity of the normal fluxes and of the pressures. We start by analyzing the case of a fracture having a fixed width-length ratio, called $\varepsilon &gt; 0$. Then we take the limit $\varepsilon o 0$ and give a rigorous proof for the convergence towards effective models. This is done in different regimes, depending on how the ratio of porosities and permeabilities in the fracture, respectively matrix scale with respect to $\varepsilon$, and leads to a variety of effective models. Numerical simulations confirm the theoretical upscaling results.

研究动机与目标

• 通过从薄裂缝几何结构严格推导，数学上证明在非饱和流模拟中使用低维裂缝模型的合理性。
• 根据孔隙率和渗透率比值随裂缝厚度 ε 的缩放方式，识别不同有效模型出现的条件。
• 在 ε → 0 极限下，建立 ε-问题（具有正宽度的裂缝）收敛到有效模型的证明。
• 为地下水流模拟中常用的离散裂缝建模方法提供物理解释一致的理论基础。
• 通过数值模拟验证理论结果，展示在 ε 较小时模型具有极高精度。

提出的方法

• 在基质块和宽度为 ε 的薄裂缝中均采用 Richards 方程建模，并在界面处保证压力和法向通量连续。
• 引入缩放假设：孔隙率比值 ∝ ε^κ，渗透率比值 ∝ ε^λ，其中 κ > -1 且 λ < 1。
• 应用严格的渐近分析与均质化技术，在 ε → 0 极限下推导有效模型。
• 根据 κ 和 λ 的取值，推导出多种有效模型，当 λ < 1 且 κ > -1 时，裂缝界面处为一维 Richards 方程。
• 采用弱形式和能量估计，证明 ε-问题的弱解存在唯一性，并证明其收敛到有效模型。
• 通过使用 van Genuchten–Mualem 参数的数值模拟，比较 ε-解与有效模型的差异，计算各子域上的 L² 误差。

实验结果

研究问题

• RQ1 当裂缝具有小但正的宽度时，在 ε → 0 限制下，根据孔隙率和渗透率比值随 ε 的缩放方式，会涌现出何种有效模型？
• RQ2 孔隙率和渗透率缩放中指数 κ 和 λ 所定义的不同参数区域，如何影响尺度化方程的结构？
• RQ3 对于由 Richards 方程控制的非退化非饱和流，能否严格证明 ε-问题收敛到有效模型？
• RQ4 在压力水头和饱和度方面，有效模型对原始 ε-问题的近似精度如何，特别是在 ε 较小时？
• RQ5 有效模型的数值解是否能捕捉到完整二维裂缝问题的关键特征，特别是流动动力学和饱和度分布？

主要发现

• 有效模型对缩放指数 κ 和 λ 依赖显著：当 λ < 1 且 κ > -1 时，裂缝表现为由简化 Richards 方程控制的一维界面。
• 当 ε ≤ 10⁻³ 时，裂缝和基质块中的 L² 误差均小于 10⁻⁵，与非线性求解器的容差相当，表明数值精度极高。
• 即使在 ε = 10⁻² 时，L² 误差仍低于 10⁻⁴，证实了尺度化模型在 ε → 0 时的鲁棒收敛性。
• ε-解中裂缝宽度方向上的平均压力水头与有效模型解高度吻合，验证了平均假设的合理性。
• 尽管界面处压力连续，由于水力传导率较低，裂缝仍保持低于基质块的饱和度，与模型预测一致。
• 该理论框架支持在实际模拟中使用离散裂缝模型，为其物理一致性提供了严格的数学基础。