[论文解读] On the Homogenization of Geological Fissured Systems With Curved non-periodic Cracks

Fernando A. Morales|arXiv (Cornell University)|Dec 14, 2013

Advanced Mathematical Modeling in Engineering参考文献 3被引用 4

一句话总结

本文提出了一种新颖的渐近均质化框架，用于具有弯曲、非周期性裂隙的地质裂隙介质，采用混合变分公式。通过将渗流阻力系数按比例缩放并将裂隙坍缩为二维流形，推导出具有界面切向流的耦合达西系统，实现了无需对裂隙网络施加几何限制假设的稳健数值处理。

ABSTRACT

We analyze the steady fluid flow in a porous medium containing a network of thin fissures i.e. width $\mathcal{O}(ε)$, where all the cracks are generated by the rigid translation of a continuous piecewise $C^{1}$ functions in a fixed direction. The phenomenon is modeled in mixed variational formulation, using the stationary Darcy's law and setting coefficients of low resistance $\mathcal{O}(ε)$ on the network. The singularities are removed performing asymptotic analysis as $ε ightarrow 0$ which yields an analogous system hosting only tangential flow in the fissures. Finally the fissures are collapsed into two dimensional manifolds.

研究动机与目标

• 开发用于模拟具有复杂、非周期性、弯曲裂隙的地质介质中稳态流体流动的数学框架。
• 克服先前均质化方法的局限性，这些方法要求裂隙几何形状具有周期性或规则性。
• 通过渐近分析在裂隙宽度 ε → 0 时消除奇异性，实现高效数值模拟。
• 推导出一个适定的尺度化模型，其中裂隙被坍缩为具有切向流的二维流形。
• 建立一种自然处理基质与裂隙之间流体交换的公式，无需在函数空间上强制施加耦合约束。

提出的方法

• 基于达西定律的混合变分公式，其中裂隙的流动阻力系数按 O(ε) 缩放。
• 在 ε → 0 时应用渐近分析，以消除由细裂隙引起的几何和材料奇异性。
• 采用几何变换，其中裂隙通过在固定方向上刚性平移一个 C¹ 连续曲线生成。
• 将三维裂隙网络坍缩为二维界面，得到基质和界面之间耦合的达西方程系统。
• 引入一种受 Showalter (2000) 启发的新记法，以高效编码复杂的界面几何形状和流动耦合。
• 通过先验估计推导极限问题，并证明解强收敛于尺度化系统。

实验结果

研究问题

• RQ1 如何将地质裂隙介质的均质化方法扩展到非周期性、弯曲裂隙网络，且无需施加严格的几何假设？
• RQ2 何种渐近缩放和变分公式可消除细裂隙系统中的奇异性，同时保持物理解释一致性？
• RQ3 如何在尺度化系统中一致地建模岩石基质与裂隙之间的流体交换条件？
• RQ4 当裂隙被坍缩为具有切向流的二维流形时，极限问题的结构是什么？
• RQ5 所提出的方法能否推广到高维情形，并扩展至时变或非对称裂隙几何形状？

主要发现

• 渐近分析通过将流动阻力缩放为 O(ε) 有效消除了原始问题中的奇异性，得到适定的极限系统。
• 极限问题包含一个耦合的达西系统，其中裂隙被表示为具有切向流的二维流形，流体交换由涉及压力和法向速度跃迁的界面条件控制。
• 该公式允许任意、非周期性、弯曲的裂隙，这些裂隙由 C¹ 曲线的刚性平移生成，显著拓宽了传统周期性或对称几何形状的适用范围。
• 该方法在 ε → 0 时实现了解的强收敛，确保尺度化模型能准确捕捉物理行为。
• 混合变分公式通过避免函数空间上的约束，实现了高效的数值实现，并能高效编码复杂的界面几何形状。
• 该方法可推广至 N 维设置，并可利用解析半群理论扩展至时变问题，尽管当前公式在高梯度裂隙情况下仍存在挑战。