钻孔化孔弹性应力分析benchmark

钻孔化孔弹性应力分析benchmark

本文介绍体积元（FVM）求解水力化三场耦合方程组的方法和结果简单分析。通过与解析解的结果对比验证程序的有效性(原始数据以及理论解由H.D Cheng提供）

1.问题描述

均质地层钻孔，孔内流体与地层压力相同，孔内流体与地层流体有浓度差，且假设钻孔内流体浓度与地层浓度平衡。初始地应力为0，即仅仅观察由于孔压产生的井眼周围应力改变（so called 诱导应力）

2.问题解析

该问题是一个几何轴对称平面应变边值问题，对于几何轴对称边值问题，可以很容易的在柱坐标系下将弹性力学中势函数改为需要的形式。该问题解析解形式简单，方便求解。

3.假设条件

 1.孔内流体与孔隙流体为稀溶液
 2.多孔介质饱和流体
 3.固体骨架满足线弹性小变形

假设条件分析：
稀溶液假定是为了线性化控制方程。饱和地层流体假定是为了避免大变形问题。因为干燥泥岩遇水膨胀分为两个阶段，晶间水化膨胀和渗透水化膨胀。前者对水化分散起主导作用，具有很高的膨胀压，是粘土矿物晶体片层强吸水导致的。渗透水化膨胀是由扩散双电层的同性电荷云相互排斥导致的，很明显这是一个可逆过程。晶间膨胀是个不可逆过程，对于孔隙材料的结构破坏是颠覆性的。真实岩层由于矿物分布非均质性，这种非均匀膨胀会产生裂纹，在连续介质框架内已然无法讨论（有人强行使用，但笔者认为不合理）。

4.控制方程

应力场：
$$G{\nabla }^2\vec{u}+\left(K_c+\frac{G}{3}\right)\nabla (\nabla \cdot u$$