[论文解读] Use of a speed equation for numerical simulation of hydraulic fractures

A. M. Linkov|arXiv (Cornell University)|Aug 31, 2011

Hydraulic Fracturing and Reservoir Analysis参考文献 2被引用 14

一句话总结

本文提出一种速度方程，通过将边界条件从裂纹尖端移开，改善了模拟水力裂缝时的数值稳定性，避免了因可忽略的滞后而导致的病态问题。e-正则化方法利用该速度方程，并将边界条件设置在前缘稍后方，实现了计算成本可控的精确且稳定的模拟。

ABSTRACT

The paper treats the propagation of a hydraulically driven crack. We explicitly write the local speed equation, which facilitates using the theory of propagating interfaces. It is shown that when neglecting the lag between the liquid front and the crack tip, the lubrication PDE yields that a solution satisfies the speed equation identically. This implies that for zero or small lag, the boundary value problem appears ill-posed when solved numerically. We suggest e - regularization, which consists in employing the speed equation together with a prescribed BC on the front to obtain a new BC formulated at a small distance behind the front rather than on the front itself. It is shown that - regularization provides accurate and stable results with reasonable time expense. It is also shown that the speed equation gives a key to proper choice of unknown functions when solving a hydraulic fracture problem numerically.

研究动机与目标

• 解决在水力裂缝模拟中，当流体前缘与裂纹尖端之间的滞后量可忽略或为零时，所引发的数值病态问题。
• 开发一种在该类条件下稳定且精确的数值模拟方法。
• 通过利用速度方程，为数值解中未知函数的选择提供系统性方法。

提出的方法

• 推导裂纹尖端的局部速度方程，该方程基于流体压力和断裂力学，控制裂纹扩展速度。
• 通过将边界条件从裂纹尖端移至其后方一小段距离，引入e-正则化方法，并利用速度方程定义新的边界条件。
• 在新边界条件下求解润滑方程，确保解恒等满足速度方程。
• 将速度方程作为关键约束，用于指导数值格式中未知函数的选择。
• 采用有限差分或有限元方法求解改进后的两点边值问题，提升数值稳定性。
• 通过数值实验验证方法，结果表明在不同条件下均具有收敛性和稳定性。

实验结果

研究问题

• RQ1 如何解决水力裂缝模拟中因滞后量趋近于零而引发的数值不稳定性？
• RQ2 速度方程在确保裂纹扩展模型的一致性和稳定性方面发挥什么作用？
• RQ3 将边界条件设定在裂纹尖端之后而非尖端本身，是否能带来更准确和稳定的解？
• RQ4 速度方程如何辅助在数值格式中合理选择未知函数？
• RQ5 在模拟中使用速度方程时，计算成本与精度之间的权衡如何？

主要发现

• 当滞后量为零或极小时，润滑偏微分方程（PDE）本身恒等满足速度方程，表明标准公式在数值上趋于病态。
• e-正则化通过将边界条件移至前缘之后，成功实现了模拟的稳定，避免了在裂纹尖端直接出现的数值问题。
• 该方法在计算成本合理的情况下产生精确结果，且在不同参数范围内表现出强鲁棒性。
• 速度方程为数值格式中选择合适未知函数提供了理论基础，提升了解的一致性。
• 数值实验表明，所提方法即使在滞后量极小或为零时，仍能保持稳定与收敛。
• 该方法适用于广泛的水力裂缝问题，显著提升了数值模拟的可靠性。