
地震波在地球内部的传播方式是描绘其隐藏结构的关键。几十年来,人们做了一个简化的假设:地球的地下是各向同性的,即波在所有方向上的传播速度都相同。然而,地质学的现实——层状的沉积物和应力诱发的裂缝——却表现为各向异性,其中波速等属性具有固有的方向性。这种方向性会扭曲我们对地下的认知,导致不准确的绘图和被误解的数据。主要的挑战在于描述这一现象的复杂性,它最多需要21个独立的弹性常数,而这些常数几乎不能提供任何物理直观性。
本文通过探讨Leon Thomsen的开创性工作来应对这一挑战。我们将深入研究一个已成为现代地球物理学基石的简化而强大的框架。在“原理与机制”一章中,您将了解到三个直观的汤姆森参数——epsilon (ε)、delta (δ) 和 gamma (γ)——并发现它们如何支配着各向异性波传播中奇特而迷人的物理现象。随后,“应用与跨学科联系”一章将展示这些参数不仅是理论构想,还被积极用于校正地震图像、推断岩石属性,并指导揭示地球奥秘的探索。首先,我们必须理解这些参数的基本原理以及它们所开启的波物理学新视角。
想象一下声波穿过一块完全均匀的钢块。无论它向哪个方向传播——向上、向下或侧向——速度都相同。从波的角度来看,宇宙在每个方向上都是一样的。这个优美而简单的性质被称为各向同性(isotropy)。很长一段时间里,这是物理学家和地球物理学家在模拟波在地球中传播方式时的默认假设。这是一个优雅的模型,并且在许多情况下都很有用。但像许多简单的模型一样,它也是极不完整的。
地球的地下并非一块均匀的钢块。它是由数百万年来沉积的岩层、构造应力造成的裂缝以及本身具有方向性排列的晶体交织而成的复杂织物。想一想一块木头。它很容易沿纹理劈开,但很难横着纹理劈开。它的强度是方向依赖的。这就是各向异性(anisotropy)的本质。穿行于地壳的地震波发现自己身处一个更像木头而非钢材的世界。一个平行于沉积岩层传播的波可能比试图横穿它们的波快得多。这种方向依赖性不仅改变了速度,还改变了波能量传播的路径,从而导致了各向同性模型永远无法预测的现象。要理解地球,我们必须首先学会各向异性的语言。
描述一个完全各向异性的材料可能是一场噩梦。在最一般的情况下,需要21个独立的弹性常数才能完全表征其行为。即使对于一个常见且更简单的情况,即横向各向同性(TI),其中存在一个特殊的对称轴(例如平坦沉积层中的垂直方向),我们仍然需要五个常数()来描述其物理特性。虽然这些数字包含了所有信息,但它们几乎不能提供任何物理直观性。如果比大10%,这意味着什么?
正是在这里,物理学家Leon Thomsen提供了一个豁然开朗的时刻。1986年,他提出了一种新语言来描述“弱”各向异性——在地质学中常见的一种情况,即方向差异存在但并不极端。他没有使用一堆晦涩的刚度常数,而是给了我们三个直观的无量纲数: (epsilon)、 (delta) 和 (gamma)。这些汤姆森参数衡量了与简单各向同性背景的偏差。如果它们都为零,我们就回到了完美的各向同性世界。如果它们是小的非零数,它们就能准确地告诉我们这个世界是如何各向异性的。
让我们来认识一下这个非凡的三参数组合:
Epsilon ():这是最直观的参数。它描述了纵波(P波)的各向异性。具体来说,是水平方向和垂直方向上P波刚度(近似于速度)的分数差异。其定义为 。一个正的意味着P波平行于材料“纹理”或层理的传播速度比横穿它们的速度要快。这是人们想了解各向异性岩石时,首先也是最显而易见的事情。
Gamma ():这个参数讲述了一个类似的故事,但针对的是两种横波中的一种,即SH波(质点水平运动)。其定义为 ,它量化了水平和垂直方向上SH波速度的分数差异。
Delta ():这是三个参数中最微妙,在许多方面也是最重要的一个。与和不同,并非由恰好沿着水平和垂直轴的速度定义。相反,它控制着P波在中间角度,特别是近垂直方向的行为。它的定义 ,涉及“离轴”刚度。这使得成为控制近垂直传播波的P波波前形状的关键参数。由于大部分地震勘探都涉及向下发送波并监听近垂直的回波,在获得正确的地震图像方面具有超乎寻常的重要性。它告诉我们波前是比我们预期的更接近椭圆形还是更接近箱形。
这些参数在抽象的底层物理学(刚度常数)和我们实际可以测量的事物(如不同方向的波速)之间架起了一座桥梁。它们为我们提供了一个简单而强大的公式,用于计算P波相速度作为与垂直对称轴夹角的函数:
这里,是垂直P波速度。这个优雅的表达式是开启各向异性波传播奇特而美丽世界的钥匙。
故事在这里发生了有趣的转折。在各向同性介质中,波的能量传播方向与波前推进方向相同。如果你向平静的池塘中投下一颗石子,圆形的涟漪会向外扩展,涟漪上每一部分的能量都沿径向向外传播,垂直于涟漪本身。能量流的方向(群速度)和波前法线方向(相速度)是完全一致的。
在各向异性介质中,情况不再如此。能量流的方向可以偏离波前法线方向。想象一队士兵正穿过一片泥泞的田地,但一侧的泥比另一侧更厚。为了保持队形(即“波前”)笔直,泥地较浅一侧的士兵必须放慢脚步。这队士兵看起来会发生转动,整个队伍前进的方向将与每个士兵面向的方向不同。
群速度角(能量路径或射线的方向)和相速度角(波前法线方向)之间的这种差异不仅仅是一个定性的概念;它是由汤姆森参数所描述的各向异性的一个直接且可计算的后果。对于弱各向异性,这个差异 可以用一个简单的关系式来近似:
这告诉我们,射线路径偏离波前法线的程度取决于和之间的相互作用。在一阶近似下,能量包的速度与相速度相同,但其方向不同。这对于旅行时层析成像等技术至关重要,在该技术中,我们测量能量在两点之间传播所需的时间。能量传播的路径由群速度而非相速度决定。
一个强有力的可视化方法是通过慢度面——即所有方向上相速度倒数()的图。在各向同性介质中,这个面是一个完美的球面。在各向异性介质中,它是一个扭曲的形状。这个面的神奇之处在于,群速度矢量——能量流的方向——总是垂直于慢度面。如果这个面是一个球面,任何一点的法线都沿着半径方向,因此群速度和相速度方向一致。如果这个面不是球面,法向矢量通常不会是径向的,这两个方向就会分离。
如果各向异性变得更强会发生什么?由和相互作用控制的P波慢度面,可能会以剧烈的方式发生扭曲。如果条件合适,这个面可能会出现“凹痕”或凹曲率区域。
这就是物理学变得真正狂野的地方。慢度面上的一个凹痕意味着存在多个不同的相速度方向(面上的点)共享同一个法线方向。在物理上,这意味着能量可以从一个震源沿着几条不同的路径传播到同一个接收器,并在不同的时间到达。这种现象被称为旅行时三叉波。想象一下一次地震或爆炸,但地震仪记录到P波分三次到达。这给地震成像带来了巨大的困惑,就像试图用一个能产生三重图像的镜头拍照一样。
波前上这些三叉波区域的边界是能量高度集中的点,称为尖角。它们是波前上尖锐的角状特征。
有没有办法知道这种美丽但混乱的行为是否会发生?同样,汤姆森参数提供了答案。关键是非椭圆度参数 (eta),定义为:
这是一个非凡的物理现象:两个参数 和 之间的一个简单关系,决定了波场是简单且唯一定义的,还是一个由多路径到达组成的复杂网络。
汤姆森参数远不止是一种方便的简写。它们代表了视角的深刻转变。它们在岩石弹性刚度张量的微观物理学与穿行其中波的宏观、可测量行为之间架起了一座桥梁。它们将各向异性波传播令人困惑的复杂性提炼为几个直观的原则,使我们能够预测、建模和理解一个并不像保龄球那么简单的世界。它们揭示了波速、能量路径和波前形状之间隐藏的关系,表明它们是同一潜在物理现实的不同侧面。这,本质上就是物理学的美丽与力量:找到支配复杂世界的简单、统一的规律。
在我们之前的讨论中,我们熟悉了汤姆森参数 、 和 。我们视其为一种绝妙直观的语言,用以描述岩石的内在构造如何能具有“纹理”,导致地震波在不同方向以不同速度传播。这似乎是一个微妙的、近乎学术性的观点。但正如我们即将看到的,这种微妙之处是通向广阔应用世界的钥匙,将这些简单的参数从纯粹的描述转变为强大的发现工具。有了这种语言,我们就能以一种前所未有的清晰度来审视地壳。
想象一下,你正试图绘制一张城市地图,但你的尺子会根据指向南北还是东西而伸缩。你的地图很快就会变得一团糟!这正是地球物理学家在绘制地球地下结构时面临的问题。我们测量地震波从震源出发,经岩层反射后返回接收器所需的时间。如果我们假设波速在所有方向上都相同(各向同性),而实际上并非如此,那么我们的深度计算就会出错。
这并非一个假设性的担忧。许多构成石油、天然气和水储层的岩石,如页岩,是由细小、扁平的颗粒在漫长的地质时期中沉降成水平岩层而形成的。这种层状结构造成了显著的垂向横向各向同性(VTI),其中波在水平方向的传播速度比垂直方向快。如果我们使用垂直传播速度来计算远离震源的岩层深度——此时波已经传播了相当长的水平距离——我们将会把该岩层放置得比实际更深。汤姆森参数使我们能够精确地量化这种效应并校正我们的图件。
当我们使用像地震偏移这样的复杂成像技术时,忽略各向异性的后果变得更加明显。偏移是一个计算过程,它将记录到的杂乱地震信号聚焦回其原始位置,从而创建出地下的剖面图像。当我们假设地球是各向同性的来进行偏移,但实际上它是VTI介质时,得到的图像就不够清晰。在所谓的角度域共成像道集(ADCIGs)上——这是一种诊断图,其中反射体的深度是根据波的角度绘制的——一个单一的平坦岩层会呈现为弯曲的。它会从中心向外“皱眉”(向下弯曲)或“微笑”(向上弯曲)。
这种弯曲不仅仅是一个错误,它是一个线索!这条曲线的形状是各向异性的直接标志。一个经典的结果表明,成像深度 作为角度 的函数,与真实深度 之间存在如下形式的表达式:
这个优美的公式告诉我们,小角度处的初始曲率由 控制,而在较大角度处偏离简单抛物线形状的程度则由组合 控制。因此,通过测量我们偏移数据中的“皱眉”形态,我们可以直接估算岩石的汤姆森参数,并用它们来正确地重新偏移数据,将曲线拉平成它应有的直线,从而使地下结构清晰聚焦。
地震波的旅程不仅仅记录在其传播时间中。它所携带的能量——以及在边界上反射的能量大小——本身就讲述了一个故事。反射振幅(或“亮度”)随炮检距(偏移距)变化的方式是一种强大的技术,称为振幅随偏移距变化(AVO)分析。
各向异性为这个故事增添了引人入胜的新篇章。当一个波撞击到两种不同VTI岩石之间的边界时,反射的强度不仅取决于它们各向同性属性(如密度和垂直速度)的差异,还取决于它们汤姆森参数的差异。在一个线性化近似中,反射系数的各向异性部分 可以用跨越边界的 和 的变化量 和 来表示:
这个方程 揭示了地震振幅中编码了关于岩石构造的信息。通过分析反射亮度随角度的变化,我们可以开始区分不同类型的岩石界面。例如,我们或许能够区分含气多孔砂岩和非多孔砂岩,因为气的存在可以改变岩石的各向异性响应。汤姆森参数给了我们解读这种微妙而关键信息的密码环。
我们已经看到,汤姆森参数在我们的地震数据上留下了它们的指纹。这自然引出了“反演问题”:我们能通过观测数据推断出这些参数吗?这正是地球物理学的真正艺术和科学所在。这并不像听起来那么简单,因为不同参数的影响可能会相互纠缠。
想象一下试图从旅行时中同时确定和。事实证明,你区分它们的能力关键取决于你的实验设计。如果你只使用小偏移距的数据,此时波近乎垂直传播,旅行时几乎完全只对敏感。你没有办法看到。如果你只使用非常大偏移距的数据,和的影响会以一种难以解耦的方式混合在一起。为了唯一地确定两者,你需要一个从近到远的宽范围偏移距。这揭示了一个深刻的原则:我们能了解世界的多少,与我们选择观察世界的方式密不可分。
在现代方法如全波形反演(FWI)中,这一挑战变得更大。FWI试图通过拟合记录到的地震波形的每一个波动来创建地球的高分辨率模型。这是一个极其复杂的优化问题。试图一次性求解背景速度()、和是注定要失败的。取而代之的是采用一种分层策略,这种策略巧妙地由波散射物理学所指导。我们首先求解控制模型大尺度、低波数特征的参数。这通常是速度,它对最稳健记录到的小角度散射数据最为敏感。一旦我们有了一个好的背景模型,我们接着求解对下一角度范围敏感的参数,即。最后,我们尝试解析,其特征在最宽、最难记录的角度处最强。这种从粗到精建立模型的策略是解决复杂反演问题的通用原则,其应用远不止地球物理学。
此外,对于给定的实验,有些东西可能根本无法知晓。描述水平偏振横波(SH波)各向异性的参数,在VTI介质中对P波或垂直偏振SV波的旅行时或振幅基本没有影响。这意味着任何数量的P波数据都无法约束;它对我们的实验是“不可见的”。要测量它,我们必须设计一个不同的实验,一个使用S波的实验。
一个好的测量不仅仅是一个数字,它是一个带有不确定性的数字。在我们运行了复杂的反演算法之后,我们对得到的和值的置信度有多高?这就是不确定性量化的领域。
一个思考这个问题的有力方式是通过贝叶斯推断的视角。我们从对参数的一个先验信念开始,这个信念可能来自地质知识或附近井的测量数据。这个先验信念具有很大的不确定性。然后,我们将这个先验信念与我们的地震数据信息结合起来。数据的作用是减少我们的不确定性,从而得到一个在真实值附近更尖锐峰值的后验信念。
这个框架使我们能够量化不同数据类型如何贡献于我们的知识。例如,仅使用有限角度范围的P波数据可能会使我们对和之间的关系存在高度不确定性和强烈的“权衡”(后验相关性)。但是,仅仅增加一个水平传播P波的测量,或一个单一的SH波测量,就可以显著减小不确定性并打破这种相关性,从而使我们对最终模型有更高的置信度。这正式化了一种直觉:更多、更多样化的数据会带来更准确的知识。
到目前为止,我们的故事都关乎波的速度。但是当波在岩石中传播时,它们也会损失能量,振幅会衰减——这个过程称为衰减。正如速度可以是各向异性的,衰减同样也可以。一个沿着岩石“纹理”传播的波可能比横穿它传播的波以不同的速率损失能量。
汤姆森框架可以被优雅地扩展来处理这种现象。通过允许参数和为复数,它们的虚部可以用来描述与方向相关的能量损失。这将地球物理模型与粘弹性和材料科学的深层物理学联系起来。品质因子,一个衡量材料传输波能量效率的指标,现在变成了,一个方向的函数。其角度变化由和的虚部控制。这是一个绝佳的例子,说明一个稳健的物理参数化如何能扩展以包含更复杂、更真实的物理世界。
因此,汤姆森参数远不止是简单的符号便利。它们是一种物理语言,使我们能够理解和校正地震图像中的失真,从地震振幅中读取关于岩石属性的新信息,智能地设计实验和反演策略,正式量化我们知识的局限,并将波速与更微妙的能量损失物理学联系起来。它们是现代定量地球物理学的基石,帮助我们照亮脚下行星黑暗而复杂的结构。