4. 节理产状数据预处理

野外测得的节理产状数据在进入赤平投影和统计分析之前,几乎不可能直接使用。一个典型的路勘或边坡编录数据集里,往往混杂着记录笔误、量测口径不一致、同一测点重复录入、倾向与走向混填、倾角超出合理范围等问题。若不做预处理就直接投图,极点会散乱、优势方位会被伪数据带偏,后续的统计结论也就失去了可靠性。本节围绕四件事展开:数据清洗、异常值识别与剔除、产状格式统一、缺失数据处理,并给出可直接套用的处理流程与判断依据。

一、先统一参照约定,再谈数据对错

产状预处理的第一条纪律是明确坐标与角度约定,并在整份数据集中保持一致。本课程统一采用:倾向为从正北起顺时针量测的方位角,取值 0°~360°;倾角为节理面与水平面的夹角,取值 0°~90°,向下倾斜为正。走向则取节理面与水平面交线的方位角,取值 0°~360°,且走向与倾向相差 90°。需要特别注意的是,走向本身具有二义性——一条走向 30° 的线,也可以记作 210°,二者是同一条直线。因此凡涉及走向的记录,必须同时保留“倾向落在走向哪一侧”的信息,否则无法唯一还原倾向。这一约定是后面所有换算和异常判断的基础,切换参照系时必须显式说明。

二、数据清洗:字段级与记录级两道关

数据清洗解决的是“格式与完整性”层面的问题,不涉及地质合理性判断。建议按字段逐项检查:测点编号是否唯一、坐标是否落在工作区范围内、岩性字段是否为空、倾向与倾角是否为数值型、是否存在文本混入(如“N30E”“约45°”“310/”等)。常见脏数据包括:全角字符与半角字符混用、度分秒与十进制度混填、Excel 自动把“310/45”识别成日期、倾向写成“310°”带单位符号导致无法参与计算。

记录级清洗主要处理重复与冲突。同一测点、同一组节理被重复录入两次,会在统计中形成虚假的密度峰值,应按“测点编号+节理组号”去重,保留首次记录或经复核后的记录。若同一测点出现两条倾向相差很小但倾角明显不同的记录,不能简单删除,应回到原始编录表核对,判断是两条真实节理还是录入错误。

建议在清洗阶段保留原始值列与清洗后值列,并增加“清洗动作”备注列(如“去重”“单位换算”“字符修正”),使每一步修改可追溯,便于复核与复盘。

三、异常值识别:范围判断与地质合理性判断并用

异常值分两类。第一类是数值越界,属于硬性错误:倾向不在 0°~360°、倾角不在 0°~90°、倾角为负、倾向或倾角为空但被填成 0。这类记录可直接标记为无效。第二类是数值在合法范围内但地质上不合理,需要结合现场认识判断,例如:倾角小于 5° 的近水平节理在陡倾岩层区大量出现;同一测点相邻节理倾向突变超过 120° 且无构造解释;某测点记录数远超同批次其他测点(可能存在整列复制)。

对成批数据,可用统计方法辅助筛查。以倾向为例,可计算圆周均值与圆周标准差,将偏离均值超过 2~3 倍圆周标准差的记录列为可疑;对倾角可用箱线图规则,将超出 Q1−1.5IQR 与 Q3+1.5IQR 的记录标出。但统计筛查只能“提示可疑”,不能自动“判定错误”——一个真实的优势组完全可能落在统计尾部。因此推荐的处理顺序是:先做硬性范围校验并直接剔除越界记录,再对统计可疑记录逐条人工复核,复核结论写入备注列。剔除比例应记录在案,若单批次剔除超过 5%~10%,说明原始编录质量存在问题,应优先复查编录流程而非继续删数据。

节理产状数据预处理:清洗 → 异常识别 → 格式统一 → 缺失处理 原始编录数据 测点号 / 倾向 / 倾角 走向 / 岩性 / 坐标 ① 数据清洗 字符/单位修正 · 去重 字段完整性检查 ② 异常值识别 范围校验 + 统计筛查 + 地质合理性复核 ③ 格式统一 走向/倾角 → 倾向/倾角 越界记录:直接剔除 倾向∉[0,360] · 倾角∉[0,90] · 空值填0 统计可疑记录:逐条人工复核后决定 走向/倾角 → 倾向/倾角 换算 已知:走向 S(0°~360°)、倾角 D、倾向侧别 若倾向在走向顺时针一侧: 倾向 = (S + 90°) mod 360° 若倾向在走向逆时针一侧: 倾向 = (S + 270°) mod 360° ④ 缺失数据处理策略 · 倾向缺失、倾角完整 → 无法还原,标记无效,不参与投图 · 倾角缺失、倾向完整 → 不插补,单独列出待补测 · 走向完整、倾向缺失 → 按侧别信息换算,可保留 · 岩性/坐标缺失 → 不影响产状统计,标注后保留 · 禁止用组内均值插补产状,避免制造虚假优势方位 输出:统一为“倾向/倾角”的干净数据集 附清洗日志、剔除清单、待补测清单

四、产状格式统一:走向/倾角与倾向/倾角的互转

不同单位、不同年代的编录表常混用两种记法:倾向/倾角(如 310°/45°)与走向/倾角(如 40°/45°)。进入赤平投影前必须统一为倾向/倾角,因为投影网上的极点位置直接由倾向和倾角确定。换算的关键是走向的二义性:走向 S 对应两个可能的倾向,即 S+90° 与 S+270°(均对 360° 取模),必须依据“倾向落在走向的哪一侧”来确定。

举例说明:某节理记录为走向 40°、倾角 45°,编录备注“倾向南东”。南东方向约为 135°,而 40°+90°=130°,40°+270°=310°,显然 130° 更接近南东,故倾向取 130°,最终产状为 130°∠45°。若备注为“倾向北西”,则取 310°∠45°。若原始记录既无侧别信息、又只有走向,则该记录无法唯一还原,应归入待补测清单,而不能随意取一个方向。

反向换算同样要掌握:已知倾向 α、倾角 β,则走向可取 α−90° 或 α+90°(对 360° 取模),两者等价。实际工作中建议统一保留倾向/倾角作为主字段,走向仅作为辅助显示字段由程序自动生成,避免人工填写造成不一致。换算时还应注意角度归一化,例如 (40°+270°) mod 360° = 310°,而 (350°+90°) mod 360° = 80°,取模运算不可遗漏。

五、缺失数据处理:不插补优于乱插补

产状数据缺失的处理原则与一般统计数据不同。倾向和倾角共同决定一个三维方向,用组内均值插补会人为制造出大量集中在优势方位的“伪节理”,直接夸大优势组的密度,这是赤平投影分析中最隐蔽的错误来源之一。因此推荐策略是:倾向缺失而倾角完整、且无走向信息可换算的,标记为无效记录,不参与投图;倾角缺失而倾向完整的,不插补,单独列入待补测清单;走向完整、倾向缺失但有侧别信息的,按上节公式换算后保留;岩性、坐标等辅助字段缺失不影响产状统计,标注后保留即可。

缺失率本身也是质量指标。若某批次倾角缺失率超过 10%,应检查是否是编录表格设计问题(例如倾角列被合并单元格或格式限制),从源头修正比事后补数据更有效。

六、一份可直接套用的预处理记录表

为使处理过程可复核,建议每条记录至少保留以下字段。下表为示例结构,数值为教学假设,非实测数据。

测点号 原始记录 清洗动作 异常判断 统一后倾向/倾角 处理结论
J-012 310/45 无 范围合法 310°∠45° 保留
J-015 走向40°,倾角45°,倾向南东 格式换算 范围合法 130°∠45° 保留
J-018 375/50 无 倾向越界 — 剔除,回查原始表
J-021 120/— 无 倾角缺失 — 列入待补测
J-023 310/45(重复录入) 去重 与J-012重复 — 删除重复项

完成上述四步后,数据集应满足:全部产状以倾向/倾角表示、倾向在 0°~360°、倾角在 0°~90°、无重复测点、每条被剔除或待补测的记录都有明确理由和备注。此时的数据才具备投图条件,后续的极点图、等密度图和优势组划分才有可信的基础。需要强调的是,预处理不是一次性动作,而应形成日志:记录每批次的数据量、剔除数、剔除比例、换算条数和缺失条数,这些数字本身就是评价野外编录质量的重要依据。