《地球物理遗传反演方法》科普文学杂谈.pdfVIP

地球物理反演理论和方法 的发展与评价 引 言 有关地球深部 的信 息 ，绝大多数来源于对地表地球物理观测资料 的解释 。因此 ，任何一 种地球物理观测 ，最终都要求解反演 问题 ，即根据各种位场 、地震波地球物理观测资料 ，来 推测地球 内部 的结构形态和物质成分 ，定量计算各种地球物理参数 （朱介寿 ， 解反演 问题是一 门广泛应用于各个科学技术领域 的基本研 究方法 ，凡涉及到观测与演绎 的科学技术 ，都会碰到这类 问题 。例如 ，天体物理学 中根据光谱 的测定研 究宇宙 的物质组成 及星体 的演化 ；天体力学 中根据人造卫星瞬时位置及速度 的观测数据计算其轨道要素 ；气象 学 中根据气象观测数据研 究大气层运动规律 ；大地测量学 中根据大地测量数据推算近代应力 场 ；地球物理学 中根据地表 的地球物理观测数据推算地下物质 的结构 、形态及成分等等 ，这 些都是反演 问题 的例子 。可见 ，相 当多的科学技术 问题 ，都 以不 同的方式 、在不 同程度上涉 及反 演 问题 。 解反演 问题 的古典方法是最小二乘法及统计学 中的回归、参数估计等 。近 年来 ，由于 在计算 中广泛应用 了信息论 、线性及非线性规划 、广义逆理论及最优化方法等一些数学工具 ， 使得地球物理反演 问题在理论和方法上都有 了重大进展 。 任何反演 问题 的计算都是建立在正演 问题 的基础之上 。对于某个地球物理 问题 ，只有得 到了它的数学解答 （解析 的或数值 的），即能计算其正演 问题时，才有可能解反演 问题 。对于 那些成 因或机制 尚不清楚 的地球物理 问题 ，是谈不上解反演 问题 的。所谓正演 问题 ，是 由理 论模型参数计算观测数据 的理论值 。而反演 问题 ，是 由实 际观测数据推算所有可能 的模型参 数 。 反演 问题是地球物理学研 究 中一个重要而又 困难 的问题 ，最主要 的困难是反演结果存在 着多解性 （不唯一性 ）。这种不唯一性不是 由于反演方法或技巧上 的缺陷引起 的，而是地球物 理 问题本身存在 的固有 困难 。困难之一是观测资料不完备 。由于地球物理观测 只能在地球表 面或其上空进行 ，无法直接获得反 映地球 内部物质 的充足 的信 息 。这就使得反演结果有 多种 可能性 ，而每一种可能性都能满足不充分 的观测数据 。地球物理观测信 息的不足是无法用数 学技巧来弥补 的。困难之二是任何地球物理观测都存在干扰和误差 。这种受畸变 的观测数据 使得反演计算不稳定 ，即观测数据 中少量 的错误可 以导致反演结果 的很大变动 。因此 ，在解 地球物理反演 问题 时 ，要根据 已有资料 ，尽可能多地 引入先验信 息 、附加约束条件 ，以减少 解 的不唯一性 。 在不 同的科学技术领域 中，对反演 问题有不 同的提法 。如 ，根据数学 的观 点 ，反演是在 给定 的子空 间中寻找与 函数空 间中的某 点最靠近 的点 ；根据统计学 的观 点 ，反演是 回归和参 数估计 问题 ；根据信息论 的观点，反演又可看成滤波或过程鉴别 问题 。尽管术语和提法不 同， 但它们都表示 同一 内容 ，即从观测数据计算地球 内部物质 的各种参数 。 下面对地球物理反演理论和方法的发展做一简要 回顾和归纳 ，并对现有 的 （传统的）一 般的反演方法做一定性的评价，其中不涉及为解决专 门问题而设计的特殊反演方法。 地球物理反演理论和方法的发展 在地球物理学发展 的初期 ，由于观测 资料 比较粗糙 ，一般 只对观测数据做 定性 的或粗 略 的定量解释 。从 世纪初到 年代 ，随着地球物理方法 的发展 ，相继提 出了各种解反演 问 题 的定量计算方法 。但这些方法大多是孤立地为解决某一特定 问题而提 出来 的。例如 ，重 、磁 位场解释 中常用到 的特征 点法 、积分法就是为解释某一孤立位场异常而提 出的一些简单反演 计算方法 。又如 ，地震学 中根据地球表面体波走 时观测来反演速度 随深度变化 的 公式 ，也是在假设速度是深度 的单调函数时的一种积分运算 。这些方法 的特 点是 ，仅 能提供有关平均化 的参量 ，如异常物体 的质量或埋藏深度等 。 近 年来 ，随着科学技术 的迅速发展 ，地球物理资料无论在空间分布范围上 （从大陆、 大气层扩展到海洋和太空 ）还是数据 的可靠性及丰 富程度上 ，都较过去有很大改进 。数学方