👨🎓 作者及单位
-
作者:未在摘要中明确列出(从致谢推测为某研究生)
-
学位:理学硕士(Master of Science)
-
专业:石油工程
-
单位:美国德克萨斯A&M大学(Texas A&M University)
-
研究组:MCERI(Model Calibration Efficient Reservoir Imaging)研究组
-
完成年份:未明确标注(参考文献中有2018年文献)
📖 中文摘要(核心内容)
研究背景
近年来,美国中西部(如俄克拉荷马州、德克萨斯州)的地震活动显著增加,主要归因于废水注入等人类活动。诱发地震的机制涉及孔隙压力变化、应力变化、应变变化以及断层活化。传统的油藏模拟软件通常不支持地质力学效应(如孔隙压力变化导致的应力应变响应),而有限元分析软件(如Abaqus)具备地质力学能力,但可能仅限于单相流。
研究目标
本研究旨在耦合流体流动与地质力学效应,使用商业软件Abaqus(有限元法) 和CMG(有限差分法) 对注水井-生产井系统(目标储层下存在断层)进行模拟对比,研究:
-
孔隙压力、应力变化、应变变化
-
不平衡注采对基底应变变化和诱发地震矩的影响
-
两种数值方法(有限元 vs 有限差分)的对比
研究方法
-
CMG:采用顺序迭代耦合(双向耦合),通过“储层孔隙度”耦合参数在流体流动模块与地质力学模块之间交换压力、应变、渗透率等信息
-
Abaqus:采用全耦合有限元法,同时求解压力与位移未知量(C3D8P单元)
-
失效准则:摩尔-库仑准则判断断层滑动
-
地震矩计算:基于应变变化、单元体积和刚度张量,使用Kanamori关联式转换为矩震级
主要结论
-
耦合流体流动与地质力学是必要的,特别是对于复杂稳定性现象(如诱发地震)
-
不平衡注采活动(产率 > 注率)会导致基底应变变化增加和诱发地震活动增强
-
应变变化对生产模式更敏感(由于压实效应),相比注入模式影响更大
-
诱发地震与应变变化的相关性比与孔隙压力变化更强
-
Abaqus与CMG结果对比:
-
两者趋势一致,但存在定量差异(主要由于断层表征差异:Abaqus为节点输出,CMG为单元中心输出)
-
移除断层后,两者结果更接近,说明断层表征方式是差异的主要来源
-
🖥️ CMG软件应用情况总结
| 项目 | 说明 |
|---|---|
| 使用软件 | CMG(未明确具体模块,应为IMEX或GEM,支持地质力学耦合) |
| 耦合方法 | 顺序迭代耦合(双向耦合) — 先求解流体流动方程,再将压力传递给地质力学模块,计算位移、应力、应变,更新孔隙度和渗透率,迭代至收敛 |
| 耦合参数 | “储层孔隙度” ϕ∗=(1−ϵv)ϕϕ∗=(1−ϵv)ϕ,连接流体流动与地质力学 |
| 渗透率更新 | 基于体积应变的经验关系:ln(k/k0)=Cnϵvln(k/k0)=Cnϵv |
| 模型规模 | 二维简化模型(Azle案例简化版),网格 62×1×2262×1×22,单注单采 |
| 失效准则 | 摩尔-库仑准则(通过材料输入间接实现,非CMG内置功能) |
| 计算效率 | 未明确对比,但有限差分法(CMG)比有限元法(Abaqus)更适合大尺度问题 |
| 与Abaqus对比 | CMG有限差分为体积中心输出,Abaqus有限元为节点输出,需后处理平均后对比 |
🧾 主要结论
物理机制
-
诱发地震的机制:流体注入 → 孔隙压力增加 → 有效正应力降低 → 断层滑动;或流体采出 → 压实 → 应变变化 → 断层滑动
-
不平衡注采的影响:产率/注率比值越大,基底应变变化和诱发地震矩越大
-
应变变化 vs 孔隙压力变化:诱发地震与应变变化的相关性更强,而非仅与孔隙压力变化相关
数值方法对比(CMG vs Abaqus)
| 对比项 | CMG(有限差分) | Abaqus(有限元) |
|---|---|---|
| 耦合方式 | 顺序迭代(双向) | 全耦合 |
| 输出位置 | 单元中心 | 节点 |
| 断层表征 | 单元属性 | 节点/单元边界 |
| 早期孔隙压力响应 | 单调变化 | 可能出现非单调(Mandel-Cryer现象) |
| 与Abaqus一致性 | 趋势一致,定量有差异 | — |
敏感性分析
-
生产-注入比(1, 2, 5, 10):比值越大,应变变化和地震矩越大
-
断层渗透率(0.003 ~ 30 md):断层渗透率增加 → 注采区连通性增加 → 压力差减小 → 应变变化减小 → 趋于稳态
📊 关键数据摘要
| 参数 | 值 |
|---|---|
| 模型尺寸 | 10×10×8 km(62×1×22网格) |
| 目标层渗透率 | 30 md |
| 断层渗透率(基准) | 0.003 md |
| 注/采速率(基准) | 416.67 m³/day(注),833.34 m³/day(采) |
| 模拟时间 | 300天 ~ 5年 |
| 杨氏模量(目标层) | 6e4 kPa |
| 泊松比 | 0.2 |
| 内摩擦角 | 30° |
| 黏聚力 | 2.0e4 kPa |
🔮 未来工作建议
-
扩展到全尺度Azle案例:2个注入井 + 70个生产井,考虑更复杂的地质力学效应
-
开发自主代码:基于基本方程、离散化和耦合技术,实现完全灵活性
-
研究压实破坏:关注最大主应力增加导致的破坏
-
指导行业最佳实践:为废水注入提供地质力学指导,尤其是在有地震历史的地区
📚 关键词(中文)
诱发地震、流体流动-地质力学耦合、CMG、Abaqus、有限差分、有限元、顺序耦合、全耦合、摩尔-库仑、断层活化、应变变化、孔隙压力、不平衡注采、Azle
🙏 致谢
-
导师:Dr. Datta-Gupta, Dr. King, Dr. Banerjee
-
研究组成员:Jaeyoung Park, Rongqiang Chen, Xu Xue, Atsushi Lino
-
资助:LSAMP-BTD奖学金(Texas A&M University)
📄 论文来源
-
机构:Texas A&M University
-
研究组:MCERI(Model Calibration Efficient Reservoir Imaging)
⚠️ 特别说明
-
本论文是已整理的14篇文献中唯一一篇以地质力学-流体流动耦合为核心的研究
-
CMG在本研究中作为有限差分对比工具,而非主要模拟平台(与Abaqus对比)
-
论文详细推导了CMG和Abaqus的离散化方程,可作为耦合模拟的参考
