2. 井数据加载与质量检查
井震标定的第一步不是做合成记录,而是把井筒数据“干净地”装进解释系统。合成记录的质量上限,几乎在数据加载阶段就已经被决定了:声波时差(DT/AC)决定反射系数序列的时间刻度,密度(RHOB)决定波阻抗的绝对值,伽马(GR)与自然电位(SP)则用于岩性识别和层位对比的辅助判断。如果曲线存在单位错位、深度偏移或异常尖峰,后续的反射系数计算和子波匹配都会建立在错误的基础上,标定结果看似“对上了”,实际上是把错误拟合成了假象。
曲线加载:先明确单位、采样与深度基准
加载前应逐条确认三类元数据。第一是单位:声波时差在国内常见 μs/m 与 μs/ft 两种,密度常见 g/cm³ 与 kg/m³,伽马常见 API 与 cps。单位不统一会导致数值量级相差数倍,例如 2.0 g/cm³ 若被当作 kg/m³ 读入,波阻抗计算将完全失真。第二是采样间隔:常规测井曲线多为 0.125 m 或 0.1 m 采样,加载时应保留原始采样率,不要在加载阶段做重采样,以免后续深度校正时丢失细节。第三是深度基准:明确曲线深度是井深(MD)还是垂深(TVD),以及是否经过补心海拔校正,整口井的所有曲线必须使用同一深度基准。
以某示例井为例,加载后的曲线头信息可整理如下(示例数据,仅用于说明字段组织方式):
| 曲线 | 单位 | 采样间隔 | 深度范围 | 用途 |
|---|---|---|---|---|
| DT(声波时差) | μs/m | 0.125 m | 1200–3100 m | 时深转换、反射系数 |
| RHOB(密度) | g/cm³ | 0.125 m | 1200–3100 m | 波阻抗计算 |
| GR(伽马) | API | 0.125 m | 1200–3100 m | 岩性识别、层位辅助 |
| SP(自然电位) | mV | 0.125 m | 1200–3100 m | 渗透层判断 |
曲线一致性检查:让不同曲线“对得上同一段地层”
一致性检查的核心是判断各曲线在同一深度上是否描述同一地质事实。常用做法有三:其一,将 GR 与 SP 叠合显示,检查渗透性砂岩段是否同时表现为 GR 低值、SP 负异常,若两者错位则说明某条曲线存在深度偏移;其二,将 DT 与 RHOB 交会或并排显示,在致密灰岩段两者应同时表现为低声波时差、高密度,若出现“声波正常而密度异常偏低”,往往是密度曲线井径影响未校正;其三,检查曲线之间的相关系数,在岩性稳定的井段内,DT 与 RHOB 通常呈负相关,若相关系数接近零或为正,应怀疑其中一条曲线单位或深度有误。
此外还要检查曲线覆盖范围是否一致。若 DT 只测量到 2900 m,而 RHOB 到 3100 m,则 2900 m 以下无法计算波阻抗,需要在加载后明确标注“有效计算井段”,避免在合成记录中误用外推段。
深度校正:以标志层为锚点做整体与分段平移
测井曲线深度与地震深度(或钻井分层深度)不一致是常见问题,来源包括电缆张力、仪器零深误差和井斜校正差异。校正方法是以区域稳定标志层为锚点,比较测井曲线响应特征与已知分层深度,求出深度偏移量。若整口井偏移量基本一致,做整体平移;若上下段偏移量不同,则分段平移并在段间做线性过渡,避免出现深度台阶。
例如某示例井在标志层 A(设计分层深度 1850.0 m)处,GR 曲线实际响应最低点位于 1852.5 m,则整体偏移量为 +2.5 m,应将全部曲线深度减去 2.5 m 后重新采样。校正后必须复查:标志层 B、C 的响应位置是否也同步对齐,若仅 A 对齐而 B 偏差增大,说明偏移量并非常量,应改为分段校正。
异常值处理:识别、剔除与合理插值
异常值主要分三类:尖峰(单点或数点突变)、平直段(仪器故障导致的常数段)和井径垮塌段引起的密度/声波失真。识别时可用滑动窗口统计,例如以 5 m 窗口计算中值与标准差,偏离中值超过 3 倍标准差的点标记为可疑。处理原则是“先标记、再判断、后插值”:尖峰直接剔除后用相邻有效值线性插值;平直段若长度小于 1 m 可插值,超过 3 m 应整段标记为无效,不得强行填补;垮塌段密度曲线应结合井径曲线(CAL)判断,CAL 明显扩径处 RHOB 偏低,应剔除并插值,同时记录无效井段。
下图概括了本节从加载到可用曲线的处理链路与判断关系:
处理记录与可复现性
完成上述四步后,应形成一份可复核的处理记录,至少包含:每条曲线的单位与采样、整体/分段深度偏移量及锚点标志层、异常值剔除的深度区间与插值方法、标记为无效的井段。这份记录既是质量检查的凭证,也是后续合成记录出现异常时回溯问题的依据。只有曲线集、深度基准和处理记录三者一致,井震标定才具备可靠的数据基础。