学完这节课,你将带走 3 个能力
1把 Δtp、Δts 转成 Vp、Vs 与动态 E
2按地层调用声波和密度经验式
3用双路径差值识别数据或模型问题
从测井估算 UCS,看似只是往 Excel 里填公式,真正决定结果可信度的却是三件事:单位是否一致、地层公式是否用对、动态参数是否经过静态校准。
这节课沿两条独立路径计算 UCS,并把二者的差异变成质量控制工具,而不是简单 平均掉。
01 / KNOWLEDGE
无侧限抗压强度 UCS 描述岩石在无侧向约束下达到峰值破坏时的轴向应力。声波和密度测井并不直接测 UCS:声波给出波传播速度,密度给出体积密度;我们先计算动态弹性模量或孔隙度,再用按地层标定的经验关系映射到 UCS。
动态弹性模量来自高频、小应变波传播,通常高于岩心静态模量。未经岩心校准,动态 E 与 UCS 的经验关系只能称为“估算”,不能替代实验强度。
02 / KNOWLEDGE
慢度 Δt 是单位距离的传播时间,速度是它的倒数。若输入单位为 μs/ft,转换到 m/s 时:
Vp = 304800 / Δtp ; Vs = 304800 / Δts
Δt:μs/ft;V:m/s。304800=0.3048×106。
有了 Vp、Vs 和体积密度 ρ,动态泊松比、剪切模量与杨氏模量为:
ν = (Vp2−2Vs2) / [2(Vp2−Vs2)]
E = ρVs2(3Vp2−4Vs2) / (Vp2−Vs2)
密度若以 g/cm³ 输入,先乘 1000 变为 kg/m³;公式输出 Pa,再除以 106 变为 MPa。
03 / KNOWLEDGE
用两端元混合关系估算孔隙度:
φ = (ρmatrix−ρlog) / (ρmatrix−ρfluid)
本例取 ρmatrix=2.7 g/cm³、ρfluid=1.0 g/cm³,因此 φ=(2.7−ρlog)/1.7。若算出 φ<0 或 φ>1,应先查岩性、气效应、井径扩大和密度校正,不能直接送入幂函数。
04 / KNOWLEDGE
下表整理了本例使用的地层边界与相关式。E 和 UCS 均以 MPa 计,φ 为小数:
这些系数是特定岩性与样本集的回归结果,不是普适物理定律。换盆地、换矿物组成、换成熟度后,指数和截距都可能变化;至少应检查原始样本范围与 R²。
Excel 中最稳妥的做法,是先按深度生成“Formation”列,再用分段 IF 或查找表调用对应公式。边界点必须提前规定归属,避免同一深度被重复或遗漏。
05 / KNOWLEDGE
将两套 UCS 统一为 MPa 后,取绝对差;若题目要求 psi,再乘 145.038。三个检查点的结果如下:
差值大并不自动说明某一条曲线“错了”。声波对裂缝、气体和各向异性敏感,密度对井径扩大、重晶石和矿物组成敏感;两者分离的位置,反而提示需要查看井径、伽马、含气性和岩性。
06 / KNOWLEDGE
连续 UCS 可进入井壁稳定、出砂、压裂与压实模型;但真正可决策的曲线,必须知道每一步的单位、适用域和校准来源。
一页结论
声波路径先由 Δt 得到 Vp、Vs 与动态 E,密度路径先由 ρb 得到 φ;两者再进入分地层 UCS 相关式。公式必须在正确单位和适用深度内使用,并用岩心静态强度完成最终校准。
建议收藏:把方法留在手边
以后再看到一条“漂亮的 UCS 曲线”,先问四件事:慢度单位、密度单位、地层边界和岩心校准。
把复杂问题,拆成可以复核的工程判断。
公式、单位、假设和结果均已在正文中对应呈现。