首页/文章/ 详情

如何建立一维机械地球模型?从测井曲线到可校准的地应力剖面

53分钟前浏览0

高级地质力学 · PGE 383 / L10

如何建立一维机械地球模型?从测井曲线到可校准的地应力剖面

机械地球模型不是把几条测井曲线叠在一起,而是把密度、声波、孔压、弹性、强度和应力放进同一深度与不确定性框架。

这一讲,你会真正弄清楚

01 由Vp、Vs和密度计算动态弹性参数
02 理解动态—静态转换为何必须校准
03 构建并质控一维MEM的完整流程

一维MEM沿井深描述岩石力学性质、孔隙压力和地应力随真垂深的变化。它服务于井壁稳定、完井和压裂等问题,但每一个输出都继承输入数据、经验转换和边界假设的误差。

图1:弹性参数可来自实验室静态或动态测试,也可由现场声波资料估算;来源不同不应混为同一尺度。  

图01 弹性参数可来自实验室静态或动态测试,也可由现场声波资料估算;来源不同不应混为同一尺度。

01 / 学习笔记

声波测井如何变成动态模量?

各向同性均匀介质中,P波与S波速度分别与纵波模量M和剪切模量G相关:M=ρVₚ²,G=ρVₛ²。由这两个量可换算动态E、ν和K。测井常给慢度Δt而非速度,必须先统一单位并取倒数。

图2:P波和S波对弹性参数的敏感性不同;缺少可靠S波时,不能无说明地声称得到了完整E与ν。  

图02 P波和S波对弹性参数的敏感性不同;缺少可靠S波时,不能无说明地声称得到了完整E与ν。

各向同性动态参数

Gdyn=ρVs²
Kdyn=ρ(Vp²−4Vs²/3)
νdyn=(Vp²−2Vs²)/[2(Vp²−Vs²)]
Edyn=2Gdyn(1+νdyn)

ρ用kg/m³、速度用m/s时模量单位为Pa。公式假设局部各向同性、线弹性和正确的波型识别。

图3:由ρ、Vp、Vs可以建立动态模量曲线;每一步都需保留原始曲线质量标志。  

图03 由ρ、Vp、Vs可以建立动态模量曲线;每一步都需保留原始曲线质量标志。

教学例题:ρ=2400 kg/m³、Vₚ=4000 m/s、Vₛ=2300 m/s。计算得到G_dyn=12.696 GPa、ν_dyn≈0.253、K_dyn≈21.472 GPa、E_dyn≈31.82 GPa。可用E=2G(1+ν)与K=E/[3(1−2ν)]交叉核验。

02 / 学习笔记

为什么动态模量通常不能直接当静态模量?

图4:剪切模量随应变幅值和加载条件变化;声波测量处在极小应变、高频一端。  

图04 剪切模量随应变幅值和加载条件变化;声波测量处在极小应变、高频一端。

声波传播的应变很小、频率高,裂纹和孔隙未必有足够时间发生与准静态加载相同的滑移或流体交换。实验室静态加载的应变幅值更大、时间尺度更长。因此动态模量常高于工程静态模量,但比例不是普适常数。

图5:不同测量方法覆盖不同应变幅值与频率;把一条经验转换用于所有岩性会掩盖真实差异。  

图05 不同测量方法覆盖不同应变幅值与频率;把一条经验转换用于所有岩性会掩盖真实差异。

经验转换只能写成待校准模型

Estatic=F(Edynamic,岩性,围压,孔隙流体,频率,应变幅值)

课程练习可指定例如E_static=0.65E_dynamic用于计算,但这一系数属于特定教学或标定条件,不能直接推广为行业常数。

图6:静—动态相关需要实验室或现场配对数据拟合,并检查残差、截距和适用范围。  

图06 静—动态相关需要实验室或现场配对数据拟合,并检查残差、截距和适用范围。

03 / 学习笔记

一维MEM需要哪些数据层?

模型层      
常见输入      
主要输出或作用      
几何与深度      
井轨迹、MD—TVD、地面/海面基准      
统一垂深和地层位置      
覆盖与孔压      
密度、浅层模型、泥浆与压力测试      
Sᵥ、Pₚ、有效应力      
弹性与强度      
Vp、Vs、岩心试验、矿物与经验式      
静态E、ν、UCS等      
应力方向与大小      
成像测井、LOT/DFIT、井壁破坏      
SHmax方位、Shmin及校准      
应用边界      
井眼轨迹、泥浆、压裂层段      
井壁或裂缝预测输入      
图7:沿深度填充弹性参数和应力曲线时,应保留层界、空缺段与数据分辨率。  

图07 沿深度填充弹性参数和应力曲线时,应保留层界、空缺段与数据分辨率。

04 / 学习笔记

从原始曲线到模型曲线的九步流程

1 · 深度匹配

完成环境校正、MD—TVD转换、不同测井趟次对齐和重复段处理。

2 · 岩性与质量控制

利用井径、自然伽马、密度校正等识别扩径、坏孔和异常声波。

3 · 总垂直应力

用体积密度沿TVD积分,补齐浅部和海水段并记录假设。

4 · 孔隙压力

结合压力测试、泥浆记录、速度或压实趋势建立模型,区分直接测量与间接预测。

5 · 动态弹性

由ρ、Vp、Vs计算,并筛除Vp≤Vs或导致非物理ν的记录。

6 · 静态转换与强度

按岩性分区用配对试验校准,不把相关式外推到未覆盖区间。

7 · 水平应力基线

选择各向同性或VTI本构及边界,计算无构造与构造应变情景。

8 · 独立观测校准

用泄漏/微压裂、井径、成像破坏和区域方向等分别约束。

9 · 不确定性传播

为关键输入设置范围或分布,报告应力窗口而非只有一条确定曲线。

图8:静态与动态模量的差异会直接传播到构造应变应力项,相关系数应按层段校准。  

图08 静态与动态模量的差异会直接传播到构造应变应力项,相关系数应按层段校准。

05 / 学习笔记

机械分层不是把曲线“画得更平滑”

机械层应反映参数和边界的真实变化尺度。过度平滑会抹掉能阻挡裂缝高度的小层,逐点保留全部噪声又会制造虚假的应力尖峰。应综合测井分辨率、岩性界面、岩心和具体工程尺度分层。

图9:地震或多井资料可帮助把井旁一维参数扩展到空间,但空间插值不能替代横向地质解释。  

图09 地震或多井资料可帮助把井旁一维参数扩展到空间,但空间插值不能替代横向地质解释。

图10:地层力学差异会影响裂缝几何;颜色剖面只有在标尺、参数和边界明确时才可用于比较。  

图10 地层力学差异会影响裂缝几何;颜色剖面只有在标尺、参数和边界明确时才可用于比较。

一维模型的边界

断层邻域、盐体边缘、复杂构造和强生产扰动需要二维或三维模型。一维MEM仍可提供初始条件和参数剖面,但不应承担它无法表达的横向应力重分布。

06 / 学习笔记

交付一份可审计MEM,应带什么?

  • 每条输出曲线的公式、单位、数据源和版本;每个经验系数的标定样本与适用区间。
  • 深度基准、坐标、总/有效应力、压缩符号和Biot系数约定。
  • 坏孔、缺失段、插值段和人工编辑标记,不能只交一张最终图。
  • 校准点的误差、模型残差及未用于拟合的验证证据。
  • 用于井壁或压裂应用时,另附强度准则、井底压力、方向与作业场景。

练习 / 学习笔记

把漂亮曲线变成可信模型

练习1|只有Vp和密度,能唯一求E和ν吗?

参考答案:不能。各向同性下仍缺少独立剪切信息,通常需要Vs或另一项弹性约束;经验估算必须标注。

练习2|Estatic=0.65Edynamic能否作为所有井的固定规则?

参考答案:不能。它可作为特定练习或经本区数据验证的经验式,需按岩性、状态和尺度校准。

练习3|MEM和测井综合图有什么根本区别?

参考答案:MEM包含明确物理模型、边界、参数转换、校准和不确定性;曲线并排本身不构成力学模型。

总结 / 学习笔记

把本讲压缩成几条判断原则

  • 动态模量由ρ、Vp、Vs在明确假设下换算。
  • 静—动态转换是需校准的经验关系,不是固定物理常数。
  • MEM必须统一深度、坐标、单位及总—有效应力定义。
  • 机械分层要匹配数据分辨率和工程问题尺度。
  • 一维结果必须由独立应力证据校准并报告不确定性。

课程资料:D. Nicolas Espinoza,UT Austin,PGE 383 Advanced Geomechanics(Fall 2020),L10。插图取自所提供的本讲材料,中文教学结构、补充解释和明确标注的教学例题为整理内容;本系列不是官方中文版。

来源:在北漂的石油人
UM裂纹材料控制试验
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2026-10-24
最近编辑:53分钟前
现代石油人
博士 签名征集中
获赞 44粉丝 175文章 1639课程 1
点赞
收藏
作者推荐
未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈