如何计算上覆应力？ 从密度曲线到 Sᵥ 剖面
学完这节课，你将带走 3 个能力
1把密度测井正确积分为 Sᵥ
2处理海水柱、浅层缺失和单位换算
3用分块近似与图形质量控制复核结果
上覆应力看似只是“岩石重量”，却是孔隙压力、水平应力与井壁稳定计算的地基。最常见的错误，是把底部一个密度乘以整段深度，或者在海上井中忘记海水柱和深度基准。
这节课从 ρgh 的单位换算开始，分别建立连续积分与五块近似，并用 Excel 曲线把错误显性化。
01 / KNOWLEDGE
上覆应力是什么：不是当前点的密度乘深度
垂向上覆应力 Sv 是目标深度以上全部物质重量在单位面积上的贡献。密度随深度变化时，应沿垂向积分：
Sv(z) = ∫0z ρb(z′)g dz′
因此，直接用目标点密度 ρ(z)×z 会把整个上覆柱都当成与底部一样致密，通常造成系统误差。正确做法是逐采样点积分，或把密度曲线分块后求和。
Sv 是总应力。进入有效应力或压实分析时，还要减去孔隙压力：σ′v=Sv−αPp。不要把上覆应力与有效上覆应力混用。
02 / KNOWLEDGE
第一道换算：密度曲线如何变成孔隙度与静水孔压
若把岩石简化为基质与孔隙流体两端元，体积密度为：
ρb = ρma(1−φ) + ρfφ
本例采用 ρma=2.7 g/cm³、ρf=1.0 g/cm³，因此：
φ = (2.7 − ρlog) / 1.7
在 5,750 ft 的数据检查中，Barnett 页岩孔隙度约 0.14383，墨西哥湾样本约 0.29765。这只是简化密度孔隙度；若矿物密度、气体或井径效应显著，需要岩性校正。
03 / KNOWLEDGE
把 SI 单位压缩成一个实用系数：0.4335
1 g/cm³=1000 kg/m³，1 ft=0.3048 m，1 Pa=0.000145038 psi。将这些单位代入 ρgh，可得到方便的油田单位积分式：
Sv(psi) = 0.4335 · ∫ρb(g/cm³) dz(ft)
Pp,h(psi) = 0.4335 · ρf(g/cm³) · TVD(ft)
0.4335 是每 1 g/cm³、每 1 ft 对应的压力增量（取标准重力时约值）。如果改用海水密度 1.03 g/cm³，静水梯度约为 0.447 psi/ft。
深度必须用真垂深 TVD，而不是斜深 MD；密度单位必须确认是 g/cm³、kg/m³ 还是 lb/gal。一个单位列错，整条应力曲线都会按固定倍数偏移。
04 / KNOWLEDGE
陆上与海上井：海水柱、未测浅层和测井段要分开积分
陆上井可直接对地表以下的体积密度积分。海上井则至少包含海水柱和海床以下沉积物；若浅部没有密度测井，还要为缺失段建立压实趋势或分块密度模型。
Sv,offshore = 0.4335[ρwzw + Σρ̄iΔzi]
海上深度基准尤其容易混乱：海平面、转盘面、泥线和测井零点必须统一。空气间隙没有有效上覆质量，海水段用水密度，海床以下才用沉积物体积密度。
05 / KNOWLEDGE
五块近似怎么做：用面积守恒替代“看起来差不多”
当浅层密度缺失，或需要快速手算时，可把密度—深度曲线分成若干块。每块选平均密度 ρ̄i 与厚度 Δzi，累加矩形面积：
Sv ≈ 0.4335 · Σρ̄iΔzi
好的分块应让矩形面积尽量接近原曲线面积。密度梯度突然变化、地层界面或海床位置应设为块边界。结果要与逐点积分比较；偏差过大时增加分块数或改用梯形法。
06 / KNOWLEDGE
Excel 工作流：中间量必须可见，图形必须能暴露错误
建议把 Depth、Δz、Density、Porosity、Pore Pressure、ΔSv、Cumulative Sv 和 Sv/z 分列，不要把所有换算塞进一个长公式。这样一旦出现异常跳点，可以迅速定位是密度、深度间隔还是单位系数。
07 / KNOWLEDGE
结果怎么报告：数值之外还要给出缺失浅层与密度模型
上覆应力往往是孔压、应力多边形和井壁稳定模型的基础输入，因此不确定性会向下游传播。报告中至少应说明密度测井起始深度、海水深度、浅层补偿方法、积分方法、重力常数与单位。
一条可信的 Sv 曲线，不是“公式算得出来”，而是每一段上覆质量都能解释、每一次单位换算都能复核。
一页结论
垂向上覆应力必须对目标深度以上的体积密度积分：Sv=∫ρgdz。油田单位可写成 0.4335∫ρ(g/cm³)dz(ft)；海上井要分开处理海水柱、缺失浅层与测井段，并用连续积分、分块近似和深度曲线相互校验。
建议收藏：把方法留在手边
拿到任何 Sᵥ 曲线，先问三个问题：密度从哪里开始、海水和浅层怎么补、积分用的是 TVD 还是 MD。
把复杂问题，拆成可以复核的工程判断。
公式、单位、假设和结果均已在正文中对应呈现。
