液层覆盖指数梯度横观各向同性介质中三维波传播的解析解与基于完全匹配层（PML）的方法

《Soil Dynamics and Earthquake Engineering》：Analytical and PML-based solutions for 3D wave propagation in a liquid layer over an exponentially graded transversely isotropic medium

　　

本文研究了由液层覆盖非均匀横观各向同性基底所构成系统中的波传播问题。针对无界基底的辐射条件（radiation condition）满足，采用了两种方法：（i）严格的解析公式；（ii）近似吸收边界条件（absorbing boundary condition）。通过截断该非均匀半空间，实现了完全匹配层（Perfectly Matched Layer, PML）作为吸收边界条件。解析方法与PML方法得到的波场均表示为半无限线积分。文中采用围道积分定理（contour integration theorem），为处理积分路径上的大量奇异性提供了稳健手段。所提出解答的精度通过若干特定退化情形下与已有结果的对比得到验证。最后，给出了若干数值结果，用以比较两种方法所得解，并描述PML参数——如近场厚度与PML厚度——在吸收入射外传波中的敏感性。参数分析表明，流体动力压力传递（hydrodynamic pressure transfer）与流体扰动（fluid perturbation）对非均匀性、激励频率以及液层厚度具有敏感性。本文建立的公式为流体-基底相互作用问题中高级区域型数值技术的实现与验证提供了可靠的基准解，例如有限元法（Finite Element Method, FEM）与无网格法（Meshless methods）；其应用场景包括海洋地震学（marine seismology）、海啸灾害预测以及海洋工程。

该文发表于《Soil Dynamics and Earthquake Engineering》，围绕“液层覆盖指数梯度横观各向同性半空间中的三维波传播”这一典型流体-固体耦合动力学问题，建立了解析解与完全匹配层（Perfectly Matched Layer, PML）截断解两套并行框架，并通过数值对比评估PML吸收无界域辐射波的有效性。研究背景在于，地震波、动力荷载波以及海洋环境中的波动传播常发生于液体覆盖地基或岩体的体系中，如水库、河流、海洋及近海工程场景。此类问题不仅涉及液层与基底的耦合作用，还受到基底材料各向异性与深度非均匀性的显著影响。传统均匀、各向同性假设虽便于求解，但对于天然土体与岩体随深度变化的力学性质而言往往不够精确，因此有必要在地壳与场地模型中引入横观各向同性（transverse isotropy）与指数梯度非均匀性（exponential grading），以提升波传播分析的可靠性。

现有研究已经分别讨论了液层上覆弹性半空间、横观各向同性基底中的表面波传播、以及非均匀各向异性介质中的散射与界面波问题；同时，解析格林函数（Green’s functions）在场地-结构相互作用与边界元法（Boundary Element Method, BEM）中具有核心作用。然而，对于兼具液体耦合、横观各向同性、指数型非均匀分布与三维无界传播特征的体系，严格解析解与适用于数值实现的PML框架仍较为缺乏。另一方面，有限元法（Finite Element Method, FEM）等区域型数值方法分析半无限域时，必须对计算域进行截断；若远场边界处理不当，易产生非物理反射，影响计算精度。正因如此，研究人员开展本项研究，目的在于同时构建严格解析基准与可实施的PML截断模型，以服务于复杂流体-基底相互作用问题的分析、解释和数值方法验证。

研究人员考虑了一个有限厚度液层覆盖于指数梯度横观各向同性半空间之上的系统。液层在水平方向无限延伸，基底在深度方向为半无限体。研究中分别建立了解析求解路径与PML求解路径：在解析方法中，依据无界域中的物理可接受性与辐射条件，剔除不符合衰减要求的项，从而在激励平面上、下区域构造满足边界与连续条件的通解；在PML方法中，则借助复坐标伸缩（complex-valued coordinate stretching transformation）和阻尼函数，将半无限非均匀基底截断为有限计算域，并通过人工吸收层实现外传波衰减。两种方法所得波场最终均写成半无限线积分形式。由于积分核中存在贝塞尔函数（Bessel functions）引起的振荡与慢衰减，且积分路径上分布着多种奇异性，研究人员进一步采用围道积分定理对奇异项进行稳健处理，并使用自适应求积进行数值积分。

用于开展研究的主要关键技术方法包括：其一，基于液层-横观各向同性非均匀半空间耦合体系建立三维波动控制方程与边界条件；其二，在解析框架中通过波场分解、辐射条件筛选与半无限线积分表达构建闭式表示；其三，采用完全匹配层（PML）复坐标伸缩方法对半无限域进行截断建模；其四，利用围道积分定理处理积分路径奇异性，并采用Wolfram数值积分程序进行自适应计算；其五，通过与特定退化情形的已有结果对比完成模型验证，并进一步执行参数敏感性分析与PML参数优化评估。文中未涉及样本队列来源问题。

关于研究结果，论文主体可按照文中主要部分进行概括如下。

Inhomogeneous substrate
研究首先明确了问题的物理构型，即有限液层覆盖指数梯度横观各向同性基底。液层厚度以h表示，体系在水平方向无限延展，基底沿深度方向为半无限。通过这一建模，研究人员将天然地质介质中材料参数随深度变化的客观事实纳入波动分析框架，为后续讨论液层-基底耦合波场奠定基础。该部分的核心结论是：在横观各向同性与指数非均匀性共同作用下，体系的传播特征不能由简单均匀半空间模型替代，必须显式考虑深度梯度与各向异性效应。

Analytical approach
在解析方法部分，研究人员依据物理可接受性去除介质III中随深度增长而发散的指数项，只保留满足辐射条件的衰减波成分，并据此重写介质II与III中的通解。通过结合控制方程及各界面、激励面上的条件，建立待定系数方程组，进而获得液层及基底中的三维波场表达。该部分表明：对于此类无界流体-固体耦合问题，可通过严格解析程序构建满足辐射条件的半解析闭式解，为后续数值验证与物理解释提供高可信基准。

PML approach
在PML方法部分，研究人员将问题域在一定深度位置截断，并在截断边界附近引入有限厚度的PML区域。该层通过复坐标伸缩与阻尼函数使位移场和应力场在PML中逐渐衰减，以满足等效辐射条件并抑制人工反射。文中强调，PML在这里的主要作用并非替代解析法、也非单纯降低相对于解析积分的计算代价，而是提供一个适用于区域型数值离散方法的伴随框架，并以解析解为基准评估其有效性。该部分结论是：对该非均匀各向异性液层-基底体系，PML能够构造出可操作的有限域近似模型，但其效果依赖于PML厚度及近场范围等参数的合理设置。

Numerical scheme
由于两种方法最终均得到半无限线积分表达，而此类积分难以直接闭式计算，研究人员采用数值积分策略进行求解。积分困难主要来自贝塞尔函数引起的高振荡与缓慢衰减，因此使用自适应求积方法处理。结合围道积分定理，数值方案能够更稳定地跨越积分路径中的奇异性。该部分说明：在严格解析公式已建立的基础上，合理的数值积分技术是将理论表达转化为可计算结果的关键环节。

Validations
在验证部分，研究人员通过若干退化情形与既有结果进行比较，以检验解析解的正确性。比较采用无量纲频率ω₀=aωρ_s/c₄₄，其中a为无量纲长度参数，可代表激励深度s或液层厚度h。验证结果表明，所建立解析公式在特殊情形下与文献结果一致，从而证明本文推导与实现具有可靠性。该部分给出的直接结论是：无论在建模假设退化为既有典型模型时，还是在数值实现层面，本文方法均保持良好精度。

Discussions
在讨论部分，研究人员首先基于解析结果考察不同模型参数对波场的影响，再通过与解析结果对照识别PML参数的优选范围。文中指出，流体动力压力传递与流体扰动对非均匀性、激励频率以及液层厚度表现出明显敏感性。这表明，基底材料性质随深度变化、外部激励特征及液层几何尺度共同控制体系响应，不同参数组合将改变液层与基底之间的动力耦合强度。此外，通过与解析解比较，研究人员进一步评估了近场厚度与PML厚度对外传波吸收性能的影响，从而揭示PML在非均匀横观各向异性介质中发挥吸收作用的机制。该部分的重要意义在于：一方面，它为理解液层-基底系统中的波动物理提供了参数化依据；另一方面，它为今后在有限元法与无网格法中选取PML计算参数提供了基准参考。

从全文的讨论来看，研究的核心贡献在于将流体耦合、横观各向同性、指数型非均匀分布以及PML截断技术统一于同一三维波传播框架中。以往研究虽分别涉及液层上覆基底、各向异性半空间、非均匀地基及PML远场处理，但较少同时处理这些因素。本文因此在理论建模和数值基准层面都具有较强价值。特别是，解析解并不仅用于给出某一特定工况下的结果，更重要的是作为严格参考标准，用于验证区域型数值方法在复杂无界流体-固体耦合问题中的精度与稳定性。对于海洋地震学、海啸灾害预测、近海工程以及一般流体-土相互作用分析，本文公式更适合作为基础性基准工具与模型开发平台，而不是替代大尺度通用数值模拟。

论文结论部分可译为如下内容：
本研究针对有限液层覆盖非均匀各向异性基底中的波传播问题，建立了解析公式与基于完全匹配层（PML）的方法。借助围道积分定理有效处理了解中的奇异性，并通过与特殊情形的比较验证了其精度。参数分析表明，流体动力压力传递与流体扰动对非均匀性、激励频率和液层厚度具有敏感性。

总体而言，该研究在《Soil Dynamics and Earthquake Engineering》发表，形成了一套兼具理论严谨性与数值实用性的三维波传播分析框架。其重要意义主要体现在三方面：其一，为液层覆盖非均匀横观各向同性半空间中的波传播提供了可靠解析基准；其二，为PML在复杂地球物理无界域中的应用效果提供了可量化评估；其三，为有限元法、无网格法及其他区域型先进数值技术处理流体-基底相互作用问题提供了验证基础与建模参照。