C4 · 出版卷 14
从岩石物性到地球物理响应
Contrast、geometry、depth、noise、survey configuration
学习目标
本章建立全书共用的物理语言。学习者应能用“物性差异假设”替代“岩石名称预期”,区分固有物性与有效调查响应,预测几何和距离对振幅与波长的影响,识别测量支持域,并在解释异常前完成简单正演检验。
需要纠正的核心直觉是:仪器并不识别岩性。仪器测量物理场、电位、计数率、电压、加速度、到时或其派生量。只有相关物性相对围岩存在差异、占据足够体积、位于敏感度足迹内并产生大于采集与环境变化的信号时,材料才可被探测。
物性差异,而不是材料标签
密度、磁化率、剩磁磁化强度、电导率、极化率、介电响应、放射性元素丰度、地震波速和衰减是物性或有效参数,但都不与某一地质名称一一对应。同名岩石可因矿物比例、孔隙、流体、温度、组构、蚀变和风化而具有很宽的物性范围;不同岩石也可能具有相同的有效物性。
必须相对于围岩或参考域定义差异。若地质体和围岩密度为 \rho_b 与 \rho_h,重力相关差异为 \Delta\rho=\rho_b-\rho_h。若二者电导率相同,即使地质接触很重要,电法也没有可探测差异。具有相同端元数值的渐变边界与突变边界会产生不同响应。
物性值必须附带测量条件。各向异性材料具有方向性;电响应可能依赖频率,激电依赖时间窗,地震速度依赖传播方向和压力,磁响应依赖外场和剩磁,放射性响应依赖地表状态与水分。实验室小样与原位体积的支持域和条件并不必然相同。
几何、深度与距离
响应取决于物性差异的位置与形状。体积、厚度、倾角、走向长度、顶面深度和边界粗糙度可与物性幅值相互补偿。一个小而强的地质体可能类似于大而弱的地质体。深部紧凑源一般比同类浅源产生更宽、幅度更低的地表响应,但这只是一般敏感关系,不是唯一的深度法则。
位场具有宽广空间影响,且通常随距离迅速衰减。波动与扩散方法还取决于发射—接收几何、传播路径和带宽。井中工具沿钻孔轨迹对局部体积取样。放射性测量主要来自浅表层,并被空气、水和覆盖层衰减。因此,“深度”在不同方法中可能指模型顶深、敏感度中心、勘探深度、特征扩散尺度、穿透程度、首波路径或工具探测半径。
应始终在明确坐标系中绘制几何,包括测量位置、传感器高度或深度、发射与接收方向、物性体与围岩、地形或钻孔轨迹,以及空气层或覆盖层。只有平面异常不能确定倾角。
正演响应与敏感度
正演模型根据物性分布和采集几何计算观测,简记为 \mathbf d^{\mathrm{pred}}=F(\mathbf m,\mathbf g)。算子可以是解析公式、数值积分、有限体积或有限元解、射线计算或卷积。其目的不只是生成逼真的图像,而是检验所提模型能否在明确物理条件下产生观测到的符号、振幅、波长、方向与带宽。
敏感度表示模型参数变化会使观测改变多少。小扰动下,\Delta\mathbf d\approx J\Delta\mathbf m,其中 J 是敏感度矩阵。敏感度高并不保证唯一性:J 的两列可能相似,使参数相互补偿。敏感度低则说明参数受约束很弱,即使模型图像很平滑也不例外。
详细反演前先做尺度试验。每次只改变一个因素:物性差异减半、深度加倍、宽度变窄、走向旋转、提高传感器高度、增大覆盖层电导率或移除某一分量。比较数据空间中的差异,而不仅是模型空间的图像。若差异小于现实不确定度,调查不能区分这些备选模型。
噪声、干扰与测量支持域
噪声包括随机仪器变化、计时误差、导航误差、平台运动、环境变化和未充分建模的背景。干扰是预期地质假设之外的结构化信号,例如基础设施、地形、近地表非均质、其他地质体或处理伪影。把所有未解释变化都称为“噪声”会掩盖可诊断过程。
测量支持域是对某个数据点有实质贡献的空间、时间与方向范围。测点坐标并不是地下的点样。支持域依赖物理核函数、传感器尺寸、积分时间、平台运动、发射波形与处理。即使地图坐标相同,若传感器高度、频率、装置长度或时间门不同,两个数值的支持域也会不同。
应按决策尺度估计总不确定度。固定点重复只测出其中一部分。复测、往返线、控制线、基站传感器、校准检查与独立导航可暴露其他分量。目标差异应与总不确定度以及目标波长处的结构化干扰比较。
合成算例
设一个合成紧凑地质体具有物性差异 p 和特征体积 V。仅为尺度检验,假定远场响应振幅与 pV/r^3 成正比,其中 r 是传感器至地质体中心的距离。模型 A 取 p=4 个任意单位、V=1{,}000\ \mathrm{m^3}、r=50\ \mathrm m,相对振幅为
$A_A=\frac{4(1{,}000)}{50^3}=0.032.$
模型 B 的体积为两倍、物性差异为一半且距离相同,响应仍为 0.032。距离增至 100 m 时,两模型响应均为 0.004,只剩原来的八分之一。这说明物性差异—体积等效与快速距离衰减,但不意味着所有地球物理方法都采用该核函数。
若相关波长处的合成标准不确定度为 0.0015,50 m 和 100 m 处信号—不确定度比分别约为 21.3 与 2.7。如果两个模型的差异只有 0.001,尽管单个异常可见,数据仍不能区分模型。设计问题不只是“能否探测异常”,而是“能否分开影响决策的备选解释”。
物性—响应审计流程
- 明确地质决策与竞争假设。
- 为每个假设赋予带单位和条件的物性分布。
- 定义围岩物性并计算差异,而不是只用岩石标签。
- 在明确坐标系中画出地质体、地形、传感器、发射源和接收器。
- 选择物理正演关系并记录其假设。
- 预测符号、振幅、波长、方向与带宽。
- 扰动物性差异、体积、深度、形状和采集几何。
- 将模型差异与总不确定度及干扰比较。
- 指出等效模型以及最能区分它们的观测。
练习与复习
- 列出三种使同名岩石具有不同电性的机制。
- 画出可能产生相似宽缓异常的浅部弱体与深部强体。
- 计算合成算例在 75 m 处的响应,并与给定不确定度比较。
- 解释为何单点仪器重复不能量化全部调查不确定度。
- 对倾斜物性边界,列出可约束倾角而非仅探测边界的观测。
复习问题:实际测量的是何种物性?相对于哪个围岩?哪个体积对数据有贡献?哪些模型参数相互补偿?关键模型差异是否大于不确定度?哪项独立测量可打破等效性?
来源
- 常见地球物理方法采用的物性关系,支撑方法与物性的区分。
- 地表地球物理方法的应用与局限,提供多方法官方概览。
- 磁化率定义与测量条件,说明物性测量方法和支持域必须记录。