用途:Bruges 面向地震与岩石物理分析,涵盖界面反射率计算、弹性属性与流体替换、解析子波以及时间域合成地震道。其定位是反射率与岩石物理计算;波传播建模应使用波动方程求解器,LAS 测井输入应使用测井读取工具。
主要能力: 反射率:bruges.reflection 提供 zoeppritz、shuey、akirichards,入射角以度为单位。shuey 支持 return_gradient=True 返回截距与梯度对;该参数属于 shuey 而非 akirichards。Shuey 截距是精确法向入射阻抗差的近似,可用 (rho2vp2-rho1vp1)/(rho2vp2+rho1vp1) 校验。临界角及以外,Zoeppritz 结果可能为复数,需保留相位或明确所绘振幅。AVO 分类应结合极性约定、角度覆盖与岩性判断,截距/梯度符号对并非唯一流体诊断。 子波与合成道:bruges.filters 的 ricker 等返回具名元组(先振幅后时间),注意字段避免混淆。合成地震道需要规则采样的双程时间;须先将深度测井转换并重采样,并保留速度/基准面选择。生成井震标定前,应在单个已知界面上检查采样间隔、时间对齐、峰值极性与振幅。 岩石物理:模量与流体替换使用 moduli.mu、moduli.pr、moduli.vp、moduli.vs 和 avseth_fluidsub;bruges.rockphysics 中没有顶层 gassmann 或 castagna。Gardner 位于 bruges.petrophysics,默认返回 kg/m³。Gassmann 假设孔隙系统连通、低频压力平衡且固体骨架与剪切模量不变,应检查密度变化与剪切模量不变性(不变量不是 Vs 不变)。
输入约定:速度用 m/s、密度用 kg/m³、模量用 Pa,报告 GPa 时再除以 1e9。提供有限层属性;仅正密度与正剪切模量不足,各向同性固体还需 vp2 > (4/3)*vs2 以保证体积模量为正。缺失测井段应保持标记,不可将 NaN 当作零值岩石属性,也不可跨缺失段做卷积。
使用方式与环境:示例经 Bruges 0.5.4、NumPy 1.26.4、SciPy 1.17.1 验证(2026-09-14)。Bruges 0.5.4 导入 pkg_resources,测试环境使用 setuptools 80.9.0;建议在隔离环境安装科学依赖。scripts/avo_analysis.py 为经验证的 AVO 辅助脚本,其 CLI 速度用 m/s、密度用 g/cm³,与上述 SI 示例不同;精确反射依赖 Bruges,缺失导入会显式失败,临界角以外的绘图保留实部/虚部分量,分类仅为筛查启发式,不能仅凭梯度符号推断 IIp。
限制:上述示例未覆盖并认证 Bruges 全部操作;其他功能、收费与平台信息暂未提供。