4. 测井资料在出砂预测中的应用
出砂预测的起点，是判断“哪一段储层可能出砂”。岩心分析最直接，但取心成本高、覆盖率低；而几乎每口井都有成套的测井曲线，它们连续记录了井筒剖面上岩石的声学、密度、电性和放射性响应。把这些响应翻译成岩石力学参数与强度指标，就能在无岩心的层段上给出连续的出砂倾向剖面。本节的核心任务，是建立“测井曲线—岩石力学参数—出砂判别指标”这条转换链，并说明如何用声波、密度、电阻率三类曲线圈定易出砂层段。
一、三类曲线各自提供什么信息
声波时差（DT 或 AC）是纵波穿过单位长度岩石所需时间，单位 μs/m 或 μs/ft。它主要受岩石骨架、孔隙度和孔隙流体影响：孔隙度越高、岩石越疏松，时差越大。由纵波时差可计算动态杨氏模量、体积模量和泊松比，是强度评价的主线曲线。若同时有横波时差（DTS），还能得到更可靠的动态泊松比；只有纵波时，常用经验关系估算横波。
密度曲线（DEN/RHOB）给出体积密度，单位 g/cm³。它一方面用于计算孔隙度（配合骨架密度和流体密度），另一方面直接参与上覆岩层压力（垂向应力）的积分计算。密度曲线的另一个重要作用是识别“低密度、高声波时差”的疏松砂岩组合——这正是出砂的高风险信号。
电阻率曲线（RT/LLD 等）本身不直接反映力学强度，但它能区分流体性质与泥质含量：油层电阻率高、水层低；同时结合自然伽马（GR）可估算泥质含量，而泥质含量又影响岩石强度与胶结程度。电阻率还常用于判断层内流体界面，帮助确定哪些层段在生产压差下更容易被水侵蚀、携砂产出。
二、从曲线到力学参数的计算链
实际工作中通常按以下顺序推进，每一步都对应可检查的输入与输出：
- 预处理：检查井径（CAL）是否规则，剔除扩径段；对声波、密度做环境校正；用 GR 或 SP 划分渗透层与泥岩隔层。
- 计算孔隙度与泥质含量：由密度曲线得 φD=(ρma−ρb)/(ρma−ρf)，由声波得 φS；用 GR 指数求泥质含量 Vsh。
- 计算动态弹性参数：由纵波时差与密度求动态杨氏模量 Ed、动态泊松比 νd，再按地区经验折减为静态参数用于强度评价。
- 估算强度：用泥质含量、孔隙度和声波时差建立抗压强度、内摩擦角、黏聚力的经验式，得到连续的强度剖面。
- 计算地应力与出砂判别指标：对密度曲线积分得上覆压力，结合孔隙压力与构造应力估算水平应力，再计算临界生产压差或出砂指数。
下面给出一个可参考的示例计算片段（参数为教学假设值，非实测）：
# 示例：由测井曲线计算动态弹性参数（教学假设数据）
# 输入：DT=320 us/m, DEN=2.10 g/cm3, 换算系数 K=1e6
# 纵波速度
Vp = 1e6 / DT # ≈ 3125 m/s
# 动态杨氏模量（GPa）
Ed = DEN * Vp**2 / 1e9 # ≈ 20.5 GPa
# 若 DTS=560 us/m，则 Vs=1786 m/s
# 动态泊松比
nu = (Vp**2 - 2*Vs**2) / (2*(Vp**2 - Vs**2)) # ≈ 0.26
得到 Ed≈20.5 GPa、ν≈0.26 后，若本地区岩心标定关系为 Es≈0.6Ed，则静态模量约 12.3 GPa，属于偏疏松砂岩区间，应列入重点观察层段。这类“算得出、对得上”的中间结果，比单纯看曲线形态更有判别力。
三、用曲线组合圈定出砂层段
单一曲线容易误判，工程上更可靠的做法是建立多曲线组合的判别规则。常见的高风险组合是：声波时差高（疏松）、密度低、泥质含量中等偏低、电阻率显示含油或含水；低风险组合则是声波时差低、密度高、泥质含量高或胶结致密。把这几条曲线按深度对齐后，可以在剖面上直接标出“甜点段”和“风险段”。
上图把四条信息放在同一深度轴上：声波时差向右偏大、密度向左偏小、电阻率向右偏高的层段，被标为高风险；致密泥岩和隔层则相反。实际解释时，还应叠加自然伽马判断泥质、叠加井径排除扩径干扰，并把判别结果与邻井出砂史、生产压差对照，形成本井的出砂风险剖面。
四、结果判读与使用边界
测井判别给出的是“相对倾向”，不是绝对结论。判读时建议关注三点：一是曲线组合是否自洽，若声波时差高但密度也高，可能是含重矿物或井眼影响，需复核；二是判别层段是否与射孔层、生产压差大的层段重合，重合度越高，实际出砂概率越大；三是用已有出砂井做回判，检验本地区判别门槛是否合适。只有经过本地化刻度与回判的判别规则，才能作为筛管选型和生产压差设计的输入依据。
把本节方法落到操作上，可以形成一张“单井出砂倾向判别表”：按层段记录深度区间、DT、DEN、RT、GR、计算孔隙度、泥质含量、动态杨氏模量、判别结论与置信度。这张表既是出砂预测的中间成果，也是后续防砂筛管选型时确定“哪一段需要挡砂、挡什么粒径”的直接依据。