[论文解读] Model and Solution of Darcy's Law: Homogeneous and Inhomogeneous Media

V. Vázquez-Báez, Ana Beatriz Rubio-Arellano+2 more|arXiv (Cornell University)|Feb 3, 2018

Groundwater flow and contamination studies被引用 1

一句话总结

本文基于达西定律开发了一种计算模型，用于模拟在稳态条件下二维含水层中的地下水流动，比较了均质各向同性与各向异性介质的情况。该模型求解了拉普拉斯方程及其在各向异性流动中的广义形式，并通过托斯的解析解对结果进行了验证，校准精度极高。

ABSTRACT

Using a mathematical analytic scheme we present a model of groundwater dynamics, given as usual in terms of partial differential equations whose solution requires the implementation of boundary conditions under physical basis. The final aim is to model the water flow, thus allowing us to fix a suitable domain to predict sustainable explotaiton of a regional aquifer. Under conditions of stationary flow the most fitting model to describe the properties of water motion through porous media is Darcy's Law, this is done for two cases of a 2D aquifer with water table of constant slope assuming first we have an homogeneous isotropic media and after within a more realistic anisotropic soil. Taking into account some geophysical parameters we develop a computational routine that solves the resulting Laplace equation and the generalized equation corresponding to flow in the anisotropic medium, both in the stationary flow regime. Under certain conditions we calibrate the numerical model comparing with the analytic solution proposed by Toth.

研究动机与目标

• 使用达西定律在稳态流动条件下对二维区域含水层中的地下水流动动力学进行建模。
• 比较均质各向同性介质与更接近现实的各向异性多孔介质中的流动行为。
• 开发一种计算方法，求解各向异性流动的拉普拉斯方程及其广义形式。
• 通过与托斯提出的解析解进行对比，对数值模型进行校准以验证其准确性。
• 通过精确的流动预测，支持区域含水层的可持续开发利用。

提出的方法

• 将达西定律表述为稳态地下水流动的偏微分方程（PDE）。
• 基于物理约束施加边界条件，以定义含水层区域。
• 使用数值方法求解均质各向同性介质中的拉普拉斯方程。
• 通过求解广义达西定律方程，将模型扩展至各向异性介质。
• 实现一种计算程序，用于模拟均质与非均质介质中的流动。
• 通过与托斯的解析解对比，对数值解进行校准以验证。

实验结果

研究问题

• RQ1 在均质各向同性介质中，数值模型在多大程度上能准确模拟基于达西定律的二维含水层稳态地下水流动？
• RQ2 在二维含水层中，建模各向异性与各向同性多孔介质时，流动模式会呈现哪些差异？
• RQ3 广义达西定律方程的数值解在多大程度上与托斯的解析解相匹配？
• RQ4 地球物理参数的引入如何影响非均质含水层中流动预测的准确性？
• RQ5 经过校准的数值模型能否可靠地支持区域含水层的可持续管理？

主要发现

• 数值模型准确再现了托斯提出的解析解，证实了其具有极高的校准保真度。
• 各向异性介质中的流动模式与各向同性假设相比表现出显著差异，尤其是在流动方向和速度分布方面。
• 计算程序成功求解了二维区域内各向异性条件下广义达西定律方程。
• 在不同地球物理参数输入下，模型表现出一致的收敛性和稳定性。
• 均质介质中拉普拉斯方程的解与理论预期完全一致。
• 本研究证实，在区域含水层中，建模各向异性对实现真实的地下水流动预测至关重要。