地球物理反演的过去,现在和未来.docVIP

地球物理反演的过去，现在和未来—— 新世纪的回顾 Sven Treitel Larry Lines (美国TriDekon公司) （加拿大卡尔加里大学地质地球物理系） 引 言 自地球物理学这一专业创建之时，地球物理学家就一直在为寻求反演问题的解决方案而努力工作。解释人员通常是根据观测到的数据例如地震记录或位场记录来推断地下的构造特征。在处理这些观测数据的过程中形成的一个初始模型会驻留在解释员的大脑中，然后，解释员会利用这个记忆中的模型根据观测到的数据来重建地下构造特征。按照现代的说法，根据观测数据对地下特征所做出的推断与所谓的“反演问题”的解决方案应该是一致的。对此相比，“正演”所包括是对在一个假定的地下构造上记录的数据进行测定并且假设这一测定过程是符合物理学法则的。直到六十年代初期，地球物理学家的大脑里几乎想的都是在做地球物理反演。从那时起，我们知道了如何在不断发展的理论研究成果的指导下，将地球物理反演做的更加定量化和多样化，并且借助于计算机的能力去缩短从理论到实践的进程。但是，我们应该强调这样一个显而易见的事实，并不是理论和计算机算法能够轻易地成为最终的裁决者来确定反演结果是否有意义或有效，最终的裁决者是地球物理解释人员。也许我们的子孙后代在下一个世纪会写一篇文章宣称计算机已解决了反演问题，而无需人来充当裁决者。但到目前为止，被称之为“无约束的地球物理反演”仍然存在但只是一个梦想。 现在让我们来遵循上述对反演问题所做的非常概括的定义来继续进行讨论。我们在处理中心所应用的各种类似的算法都可以被认为是在各个过程中来反演地球物理数据。例如，地震偏移是试图根据地震记录来重构实际的地下构造几何形态（Garder，1985）。通过预测反褶积来衰减多次波（Peacock和Treitel，1969）或者通过对地层脉冲响应中的一次反射或多次反射进行模拟（Lines和Treitel，1984）都可以对地层的反射信号进行反演。AVO技术（Castagna和Backus(1993)所做的反演是根据地表振幅的观测值来反演岩石特性。反演可以处理不同类型的地球物理数据。人们可以同时或连续将不同的地球物理数据例如地震、位场、井眼数据与同一个地层模型进行拟合（Lines等1988）。其它各种反演的实例还可以列举很多。对于每一种反演，我们假定物理法则都是成立的。举例来说，对于地震反演而言，这个法则就是波动方程或者是诸多近似方程中的其中之一。这样，基于这个物理法则的算法就能使我们将观测到的数据反演为地下岩层特征，而我们也把这一结果放在首要的位置上。 我们所说的反演是一个定义非常概括的反演，这个定义没有沿用一般文献所讲，肯定与勘探地球物理学文献中的定义不尽相同。在概述了反演方法的理论背景之后，我们来讨论一些在勘探地球物理学中最常用的一些技术，这些技术一般都被视为是“反演方法”。然而，我们必须重申，这些方法的判定和划分常常是任意的，我们日常所做的许多事都归落在地球物理反演这个保护伞下。 理 论 背 景 反演可以被定义为获得地下模型的一个过程，而这个模型足以描述所观测到的数据体。在地球物理数据中，我们所观测到的东西可以被称之为是地下构造的物理特征，包括由于地震或者电磁源的激发而形成的构造的反射（或散射）波场以及构造的重力异常和磁场异常等等。现代反演理论的基础可以认为是Backus和Gilbert（1967，1968，1970）的早期研究工作中创立的。 反演过程与正演模拟密切相关。正演模拟利用了数学关系例如波动方程来人工合成一个地球模型对激发（如地震能量脉冲）的响应。这样的方法是按照一组参数来定义的，例如层速度和层密度。当然，至关重要的是如何选择一个能够足以描述观测物体的正演模拟方法。在地震技术中，正演模拟是用一个产生合成地震记录的算法来实现的，例如，地震射线追踪器、有限差分或有限元波动方程求解器。在重力数据中，正演模拟过程包含了一个编码，该编码根据假定的地下密度分布值来计算重力场。除了选择合适的数学模型之外，重要的是要知道应该选用多少种模型参数并且哪些参数是最有意义的。选择正确的模型取决于目前要解决何种勘探问题，例如水平层模型可能适合于堪萨斯中部的地质情况，但并不适用于怀俄明或阿尔伯达山前的逆冲断裂带。 反演或者说“反演模拟”是试图根据一组已知的地球物理观测数据来重建地下构造特征，因此，反演的要求是在允许的观测误差范围之内模型响应应该与观测结果相吻合。为此，选择一个“好的”模型是至关重要的。即使假设我们选择的模型已经足够好，也仍然会存在很多问题。实际上，Jackson（1972）已经巧妙地将反演描述为“是对不精确的、不充足的、不一致的数据所做的解释”。 为了进一步说明这些问题，我们最好利用一些符