📝 摘要
本研究通过油藏表征、油藏模拟和数据分析，研究提高成熟油田原油产量的多种途径，重点关注西德克萨斯二叠纪盆地。研究为优化水驱和CO₂提高采收率、最大化增产措施效益以及确定残余油带最佳开发策略提供了指导。针对电容电阻模型在大规模成熟油田应用中遇到的求解时间过长、收敛困难等问题，研究提出了多项数值改进：推导了CRMP目标函数的解析梯度向量和Hessian矩阵，采用全局优化算法验证了目标函数的凸性，并引入参数缩放策略。这些改进使求解时间减少90%以上（以Mallet Unit为例，从数千小时降至约3小时）。在增产候选井筛选方面，通过对Slaughter油田四个租区144口井的生产数据分析和数值模拟验证，发现增产前产油量与增产增量产油量之间存在显著正相关（相关系数0.4-0.9，p值<0.05），提出了以历史产油量作为增产候选井快速筛选标准的策略。在残余油带方面，基于Hubbert水动力学理论建立了Permian盆地San Andres油藏尺度的ROZ形成与开发模型，模拟了Altered Hydrodynamic Flow Fields过程（10万年时间尺度），匹配了现场观测的OWC倾角（30英尺/英里）和饱和度剖面，评估了六种CO₂注入策略，发现MPZ与ROZ同步开发且配合充足CO₂供应时效果最优（CO₂封存量98.9 Bscf，累计产油11.9 MMstb，净CO₂利用率8.3 Mscf/stb）。研究还讨论了ROZ中CO₂封存的水处理需求和泄漏风险降低等工程问题。
🔍 关键词
成熟油田；电容电阻模型（CRM）；增产候选井筛选；残余油带（ROZ）；CO₂提高采收率；水动力学；历史拟合
🖥️ CMG软件应用情况总结
| 项目 | 内容 |
|---|---|
| 所用软件 | CMG IMEX（黑油模拟器）；CMG WinProp（相态计算） |
| 模型类型 | 黑油模型（IMEX）；组分PVT模型（WinProp） |
| 主要应用 | 1) 验证BHP缺失对CRMP连接性影响（Sacroc油田模型，8P6I）；2) 增产候选井筛选策略验证（反倾构造模型，86口生产井+44口注入井）；3) ROZ形成与CO₂-EOR模拟（2D San Andres模型，100×1×60网格） |
| 模型参数 | 平均孔隙度0.12-0.13，平均渗透率15-17 mD，Kv/Kh=0.1，油藏温度105-120°F，初始压力2000-3000 psi |
| 主要输出 | 累计产油量、含水率、CO₂封存量、净CO₂利用率、饱和度剖面、OWC倾角 |
论文中明确指出：CMG IMEX用于黑油模拟，CMG WinProp用于生成黑油PVT表（第3.5节、第4.3节）。致谢中感谢CMG提供学术许可。
📊 研究内容与方法
| 研究方向 | 核心方法 | 主要贡献/结论 |
|---|---|---|
| 电容电阻模型（CRM） | 推导CRMP解析梯度与Hessian；全局优化验证凸性；参数缩放；步进式拟合 | 求解时间减少90%以上；BHP缺失时连接性误差<20%；可应用于>500井的大规模油田 |
| 增产候选井筛选 | 相关性分析（文献数据+144口现场井+数值模拟）；两年前后增量产油量定义 | 增产前产油量与增量产油量显著正相关（r=0.4-0.9）；高产量井应是优先候选 |
| 残余油带CO₂-EOR | AHFF过程模拟（10⁵年）；6种开发策略对比（CO₂注入量140/340 Bscf，同步/顺序开发） | OWC倾角匹配Hubbert公式（30英尺/英里）；同步开发+充足CO₂最优（累计油11.9 MMstb）；ROZ具有额外CO₂封存潜力 |
CRM的关键改进：
-
推导了解析梯度向量（∂F/∂Q、∂F/∂τ、∂F/∂f）和Hessian矩阵
-
计算时间从数千小时（Mallet Unit）降至约3小时
-
建立连接性与时间常数的物理约束（τ>0，f≥0，∑f≤1，半径影响范围约束）
增产候选井筛选的核心数据：
-
Slaughter油田四个租区：Mallet Unit（23口井，r=0.42，p=0.04）、East Mallet Unit（10口井，r=0.62，p=0.05）、W.A. Coons（12口井，r=0.88，p=0.0002）、F.L. Woodley（10口井，r=0.55，p=0.10）
-
模拟验证中“高预产”方案（优先增产高产量井）捕捉了62%的总增量潜力，接近“完美信息”方案的65%
ROZ模型关键参数与结果：
-
模型尺寸：10000×600×600 ft，网格100×1×60（Δx=100 ft，Δz=10 ft）
-
平均水平渗透率15 mD，Kv/Kh=0.01-1.0（基准0.5）
-
势面坡度7.5英尺/英里，放大因子3.8，OWC倾角30英尺/英里
-
MPZ与ROZ同步开发+340 Bscf CO₂（Case 6）：累计油11.9 MMstb，CO₂封存98.9 Bscf，净利用率8.3 Mscf/stb
ROZ开发策略对比（六种方案）：
| 方案 | CO₂注入量 | 开发策略 | 累计油 | CO₂封存 | 净利用率 |
|---|---|---|---|---|---|
| Case 1 | 60 Bscf | 仅MPZ | 5.8 MMstb | 40.3 Bscf | 6.8 Mscf/stb |
| Case 2 | 140 Bscf | 仅MPZ | 8.0 MMstb | 62.5 Bscf | 7.8 Mscf/stb |
| Case 3 | 140 Bscf | 顺序MPZ→ROZ | 8.4 MMstb | 55.9 Bscf | 6.7 Mscf/stb |
| Case 4 | 140 Bscf | 同步MPZ-ROZ | 8.2 MMstb | 56.8 Bscf | 7.0 Mscf/stb |
| Case 5 | 340 Bscf | 顺序MPZ→ROZ | 10.1 MMstb | 88.5 Bscf | 8.8 Mscf/stb |
| Case 6 | 340 Bscf | 同步MPZ-ROZ | 11.9 MMstb | 98.9 Bscf | 8.3 Mscf/stb |
ROZ中CO₂封存的关键考虑：
-
封存能力评估：MPZ约0.4 Bscf/ft，ROZ约0.15 Bscf/ft（减少60%，因高渗透通道导致垂向波及效率低）
-
盐水处理需求：MPZ-ROZ开发日均产水1.35 Mstb，比仅MPZ开发（0.45 Mstb）增加200%
-
泄漏风险降低：大多数井仅完井至MPZ（OWC以上），未穿透ROZ，降低了通过废弃井泄漏的风险
关于Hubbert倾斜公式的讨论（附录A）：
-
推导了势函数Φ = gz + p/ρ，建立了水动力学条件下油水界面的倾斜公式tanθ_c = [ρ_w/(ρ_w-ρ_o)]·(dh_w/dx)
-
通过敏感性分析验证了储层非均质性、毛管压力和垂向/水平渗透率比对OWC倾角的影响较小，而势梯度是主导因素
-
模拟结果与Hubbert公式预测值吻合良好，验证了公式在复杂储层条件下的适用性
✅ 主要结论
-
CRM的数值改进（解析梯度与Hessian、参数缩放）使大规模油田（>500井）应用成为可能，求解时间减少90%以上。
-
BHP数据缺失时，CRMP识别的连接性（尤其大连接性）误差<20%，仍可用于工程决策。
-
产油历史是增产潜力的有效预测指标，“高产量井优先”策略的筛选效率接近完美信息场景。
-
水动力学（AHFF）可解释Permian盆地San Andres储层中ROZ的形成，模型匹配现场OWC倾角和饱和度剖面。
-
MPZ与ROZ同步开发且配合充足CO₂供应时，EOR和封存效果最优；CO₂投资有限时应优先顺序开发。
-
降解作用显著降低ROZ的CO₂封存能力（从98.9 Bscf降至75.9 Bscf），成岩作用略有改善（升至103.2 Bscf）。
-
ROZ开发需充分考虑盐水处理需求（比仅MPZ开发增加200%）和泄漏风险降低等工程因素。
🏛️ 作者及单位信息
| 作者 | 单位 |
|---|---|
| Ali Jamali | 德克萨斯理工大学（Texas Tech University），Bob L. Herd石油工程系 |
导师：Dr. Amin Ettehadtavakkol（委员会主席）
学位：哲学博士（石油工程）