基于测井数据无监督异常检测的奥陶系风氧化带识别

Identification of weathered zones at the top of the ordovician limestone using unsupervised anomaly detection of well-logging data

  • 摘要:
    背景 准确识别奥陶系顶部致密风氧化带对华北型煤田的水害防治工程至关重要。然而,针对这一特定类型的风氧化带,现有识别方法常面临标定数据有限、测井响应多解性强及传统经验判识主观性大等难题。
    方法 基于奥陶系顶部风氧化带的地层赋存规律与水文地质演化特征,构建由顶界面地质先验划定、风氧化带水文地质特征测井响应判识以及底界面验证分析组成的风氧化带3阶步进判识体系,明确其在自然伽马(GR)、深侧向电阻率(LLD)、自然电位(SP)以及体密度(DEN)等测井指标的组合响应特征及地层约束。在此基础上,将风氧化带识别转化为深部原岩背景约束下的异常检测问题,提出融合多尺度卷积、注意力机制与Transformer的自编码器模型MS-TransAE,采用“深部原岩训练、浅部异常识别”策略,驱动模型学习未风化原岩的测井特征模式,并结合极值理论自适应确定异常阈值,进而实现奥陶系致密风氧化带的识别。
    结果和结论 对河北东庞矿21口钻孔开展方法应用,识别结果表明矿区奥陶系顶部潜在风氧化带厚度为0~26.35 m,平均厚度17.69 m,模型识别结果的绝对误差均小于1.5 m。各模块消融实验表明,多尺度卷积、多尺度注意力、测井曲线形态特征及岩性分类头对模型的重构性能、阈值稳定性和正则化学习方面均具有重要作用。对比实验结果表明,MS-TransAE模型的识别结果较其他常规无监督异常检测模型在平均EMA上降低了77.3%~83.0%,在平均IoU上提升了20.7%~46.9%。基于所识别的风氧化带层位,对东庞煤矿区域治理注浆层位进行了率定,注浆层位高程−819.55~−289.92 m,换算至奥陶系顶界面以深0~69.04 m,平均深度26.98 m。该结果可为华北型煤田同类条件下奥陶系顶部致密风氧化带的识别及底板水害治理提供参考。

     

    Abstract:
    Background Accurate identification of the dense weathered zone at the top of the Ordovician strata is essential for mine-water hazard prevention and control in North China coalfields. However, for this specific type of weathered zone, the applicability of existing identification methods remains limited, mainly due to insufficient calibration data, ambiguous logging responses, and the subjectivity of conventional empirical interpretation.
    Methods Based on the stratigraphic occurrence patterns and hydrogeological evolution characteristics of the weathered zone at the Ordovician top, a three-stage stepwise identification framework was established. The framework consists of geological-prior delineation of the top boundary, logging-response interpretation of hydrogeological characteristics, and verification of the bottom boundary. The combined responses of natural gamma ray (GR), deep lateral resistivity (LLD), spontaneous potential (SP), and bulk density (DEN), together with stratigraphic constraints, were clarified. On this basis, weathered-zone identification was formulated as an anomaly detection problem constrained by the unweathered bedrock background. A multi-scale convolutional, attention-enhanced, Transformer-integrated autoencoder model, termed MS-TransAE, was proposed. By adopting a “unweathered bedrock for training and shallow strata for anomaly identification” strategy, the model learns the logging patterns of unweathered bedrock and adaptively determines the anomaly threshold using extreme value theory.
    Results and Conclusions  The proposed method was applied to 21 boreholes in the Dongpang Mine, Hebei Province. The identified thickness of the potential weathered zone at the Ordovician top ranges from 0 to 26.35 m, with an average thickness of 17.69 m. The absolute errors of the identified bottom boundaries are all less than 1.5 m. Ablation experiments show that multi-scale convolution, multi-scale attention, logging-curve morphological features, and the lithology classification head play important roles in improving reconstruction performance, threshold stability, and regularized learning. Comparative experiments further show that, compared with conventional unsupervised anomaly detection models, MS-TransAE reduces the average EMA by 77.3%–83.0% and improves the average IoU by 20.7%–46.9%. Based on the identified weathered-zone horizons, the grouting horizons for regional grouting treatment in the Dongpang Coal Mine were calibrated. The calibrated grouting elevations range from −819.55 to −289.92 m, corresponding to depths of 0–69.04 m below the Ordovician top, with an average depth of 26.98 m. These results provide a reference for identifying dense weathered zones at the top of the Ordovician strata and for floor water-hazard control under similar geological conditions in North China coalfields.

     

/

返回文章
返回