压缩质面法在海区地壳厚度估算中的应用与思考

压缩质面法的壳厚估算，有点意思，适合海区这类数据比较稀的区域。它的思路其实挺朴素的：用地壳的密度和地表重力异常反推厚度，数学模型也不复杂，写个fortran脚本就能搞定。你要是做地壳结构反演，这方法挺省事的，不用非得跑全波形反演那套大阵仗。

压缩质面法的计算，关键在密度分区设定。一般用上下地壳+地幔三层结构就够了，密度值定合理点，结果差不了太多。嗯，记得单位要统一，尤其重力单位经常坑新人。起来，建议先把布格重力异常数据拉平，再开始算。

要找灵感或者参考思路，可以看看这篇造山带壳厚估算的综述，里头讲得蛮清楚。你如果顺带要做点CFD或面元法模拟，下面这些也可以一起逛逛：

• CFD 面元法
• 响应面法优化
• 面元法 fortran

，这种方法图个快、省、稳。对精度要求不是变态的项目，压缩质面法真的还挺香的。如果你刚好在做海区构造解释，可以试试看。