有效应力计算
煤矿开挖区的三维渗流仿真 煤矿开挖区模型。 计算了渗流速度场以及结构的应力场
煤矿开挖区的三维渗流仿真是个挺有意思的活。咱们今天聊点硬核的，直接上代码。先说说模型构建，这玩意儿得考虑煤层裂隙网络。用Python的pyvista库处理地质数据挺顺手：
import pyvista as pv mesh = pv.read('coal_seam.vtk') mesh.plot(show_edges=True, cpos='xz')这段代码读入VTK格式的煤层模型，用三维可视化确认结构完整性。注意cpos参数控制视角方向，这里选xz平面方便观察层理走向。遇到过实际案例里钻孔数据生成的网格存在非流形边，得用mesh.clean()做预处理。
渗流场的计算核心是达西定律。用FEniCS做有限元求解时，边界条件设置是关键。比如底部承压水层的处理：
from dolfin import * mesh = Mesh('coal.xml') V = FunctionSpace(mesh, 'P', 1) bc = DirichletBC(V, Constant(2e6), "bottom")这里2MPa是实测水压值。注意单元类型选P1线性单元，计算效率与精度的折中。求解器用Krylov迭代配合AMG预条件器，比直接求解快三倍不止。
应力场耦合这块更有意思。渗流-应力耦合方程得考虑有效应力原理，代码里这么实现：
sigma_eff = sigma_total - alpha*p*Identity(3) # 渗透率动态更新 k = k0 * exp(beta*(sigma_eff.tr()/3 - sigma0))其中alpha是Biot系数，beta来自实验室拟合参数。有个坑是当单元受压过大时渗透率可能突变，得加个clip限制变化幅度。
可视化结果时用ParaView的StreamTracer滤镜追踪渗流路径，对比应力云图发现：最大主应力方向与水流优势通道呈45°夹角，这解释了现场观测到的偏流现象。有个小技巧是把计算结果存成XDMF格式，时间序列数据能动态播放。
最后给个实用建议：在Ubuntu系统下用mpirun跑并行计算，16核工作站能把30万单元的模型计算时间从8小时压到40分钟。内存不够的话试试HDF5分块存储，配合Out-of-Core算法，亲测有效。