弹性介质中非线性地震波动力学：Klein–Gordon模型的解析研究

《Results in Engineering》：Nonlinear Seismic Wave Dynamics in Elastic Media: Analytical Investigation of the Klein–Gordon Model

　　

地震波传播在石油与天然气勘探、地震工程及环境监测等诸多地球物理应用中起着至关重要的作用。地震波的数学分析可提供关于地下结构的基本信息，并对理解弹性介质中非线性波动力学有重要贡献。在非线性演化模型中，Klein–Gordon（KG）方程因具备刻画各类数学物理分支中色散与非线性波现象的显著能力而受到广泛关注。本研究通过若干互补的解析途径，研究所关联的与地震波传播相关的非线性Klein–Gordon模型。具体而言，研究人员系统地采用F-展开法（F-expansion method）、完全判别系统法（complete discriminant system method）及能量平衡法（energy balance method）来推导并分析控制模型的精确波解。所应用的方法在适当的参数区间内生成了多种具有物理意义的波结构，包括亮孤子（bright solitons）、暗型孤子（dark-type solitons）、W型剖面及周期行波解。所得结果揭示了色散弹性介质中非线性地震波的丰富动力学行为。此外，文中给出图形可视化以说明所推导解的传播特性与结构性质。本研究的结果丰富了非线性Klein–Gordon型地震模型的解析理论，并可为非均匀弹性环境中非线性波传播现象的进一步研究提供有用框架。

研究背景方面，非线性偏微分方程（nonlinear partial differential equations, NPDEs）在科学工程复杂物理现象数学建模中具有基础地位，相较于线性模型，NPDEs能描述波相互作用、能量局域化、振幅依赖传播及相干结构涌现，已广泛应用于流体动力学、等离子体物理、非线性光学、凝聚态物理及地球物理科学等领域。在众多NPDE中，Klein–Gordon（KG）方程最初由Oskar Klein与Walter Gordon于1920年代末作为薛定谔方程的相对论性推广引入相对论量子力学，随后在场论、非线性光学、等离子体物理及固态物理中找到重要应用。KG方程在同一数学框架内纳入色散与非线性效应，可描述孤立波、周期波及扭结型解等多种非线性波结构，近年来广义KG方程在弹性与地球物理介质中非线性波传播建模受到相当关注。地震波传播在地球物理研究中地位显著，由地震、火山活动、构造运动及人工震源产生，传播过程提供地球内部结构、组成与地质建造力学特征信息。地震波一般分为体波与面波，体波含纵波（P-waves）与横波（S-waves），面波含Love波与Rayleigh波。弹性介质中地震波数学描述通常由连续介质力学导出的弹性动力学方程给出，关联位移场、应力张量与应变张量，在均匀各向同性假设下可化简为经典波方程。线性波方程虽能成功描述诸多地震现象，但难以捕捉真实地质环境中大振幅扰动、非线性材料响应及不同波模间复杂相互作用等特征，因此NPDEs为地震波演化与传播动力学提供更合适框架。目前针对非线性地震波模型的解析研究相对有限，多数现有研究聚焦于孤立精确解的构造，而未在非线性地震波传播背景下将解析求解程序与定性波分析及基于能量的表征相结合。为此，研究人员对耦合非线性Klein–Gordon系统进行解析与定性综合研究，系统采用F-展开法（F-expansion method, FEM）、完全判别系统法（complete discriminant system method, CDSM）及能量平衡法，推导精确解析解并通过能量平衡途径考察周期地震波的振荡传播属性，系统分析控制模型参数对波幅、传播速度、局域化及波剖面的影响，以丰富非线性Klein–Gordon型地震模型的解析理论并提供非线性地震波动力学的物理洞察。本文发表于《Results in Engineering》。

为开展研究，研究人员采用的关键技术方法主要包括：首先将三维弹性动力学公式约化为（1+1）维纵向地震P波传播波方程，引入含非线性应力-应变相互作用、振幅依赖响应及色散机制的非线性项F(V)，并在特定参数选取下得到耦合非线性Klein–Gordon系统；随后引入行波变换V(x,t)=v(ξ)exp(i(?α₂t+α₁x+Ω))，P(x,t)=p(ξ)，Q(x,t)=q(ξ)，ξ=x??t，将原系统化为关于v(ξ)的非线性常微分方程；继而分别应用F-展开法（FEM），将解表示为辅助函数Φ(ξ)的有限多项式并利用其满足的一阶非线性常微分方程在不同参数下对应Jacobi椭圆函数解，通过齐次平衡确定展开阶数并解代数方程组得系数；应用完全判别系统法（CDSM），将约化常微分方程化为(v′(ξ))²=v³+q₂v²+q₁v+q₀形式，借助三次多项式Ξ(v)的判别式Δ与D对根的配置分类，得到有理、三角、双曲及椭圆型等不同类别精确解；此外采用能量平衡法分析所获非线性波解的定性行为与传播特征。

研究结果部分保留原文小标题并依次说明如下：

2. Theoretical Framework and Model Description：研究人员从均匀各向同性弹性介质中运动控制方程出发，依据牛顿第二定律给出x方向位移分量动量方程，通过Hooke定律建立应力张量σ_ij=λε_kkδ_ij+2με_ij与应变张量ε_ij=(1/2)(?_iV_j+?_jV_i)的本构关系，结合动量平衡导出三维形式并化简为二维耦合系统，再在均匀介质假定下令Lamé系数λ、μ及密度?为常数，进一步约化为（1+1）维纵向地震P波模型?V_x,tt+F(V)=(2μ+λ)V_x,xx，其中非线性项F(V)可表为多种形式如αV、βV³、αV+βV|V²|等。略去非线性项得经典P波波速c=√((2μ+λ)/?)，对Poisson固体取λ=μ=3×10¹⁰ Pa、?=3×10³ kg/m³得c≈4.47 km/s，1 Hz下波长λ_c=4.47 km。恢复非线性情形F(V)≠0，将方程写为V_tt?c²V_xx+F(V)=0，并取F(V)=a(V,P,Q)+b|V|^n?1V得具n-散焦非线性的二维Klein–Gordon型方程V_tt?c²ΔV+a(V,P,Q)=?b|V|^n?1V，在特定参数a(V,P,Q)=V(P?Q)、b=?2/(1?ρ²)、n=3下化为耦合系统V_tt?c²V_xx+V(P?Q)=2/(1?ρ²)|V|²V，P_x+1/(ρ+1)(|V|²)_t+(ρ+1)P_t=0，Q_x?1/(ρ?1)(|V|²)_t+(ρ?1)Q_t=0（ρ≠±1），以此作为后续解析研究之控制模型。

3. Framework of the Applied Methods：研究人员概述所采用解析方法框架，对一般NPDEs系统作行波变换ξ=x??t，将原方程化为常微分方程组，分离实部与虚部得关于v(ξ)的二阶非线性常微分方程(?²?c²)v″+(c²α₁?α₂²+p?q)v?2/(1?ρ²)v³=0，并由虚部得α₂=α₁c²/?，同时由耦合关系得p=??/((ρ+1)(?(ρ+1)?1))v²，q=?/((1?ρ)(?(1?ρ)?1))v²，代入后得最终约化方程(?²?c²)v″+α₁²c²(1?c²/?²)v+(??/((1?ρ)(?(1?ρ)?1))+?/((ρ+1)(?(ρ+1)?1))?2/(1?ρ²))v³=0。随后分别阐述F-展开法步骤：设v(ξ)=∑_l=0^kA_lΦ^l(ξ)，Φ满足( dΦ(ξ)/dξ )²=PΦ⁴+QΦ²+R，依参数(P,Q,R)取不同得Jacobi椭圆函数解（如SN(ξ,m)、CN(ξ,m)、DN(ξ,m)等），通过齐次平衡定k，代入并令同幂系数为零解代数系统得系数A_l；进而简述完全判别系统法：将约化ODE化为(v′(ξ))²=v³+q₂v²+q₁v+q₀，定义Ξ(v)=v³+q₂v²+q₁v+q₀，引入判别式Δ=?27(2q₂³/27+q₀?q₁q₂/3)²?4(q₁?q₂²/3)³与D=q₁?q₂²/3，据Δ、D符号对多项式根的配置分类以得不同精确解。

4. Derivation of Soliton Solutions：研究人员对方程(3.8)作齐次平衡得k=1，故F-展开取v(ξ)=A₀+A₁Φ(ξ)。代入并配平系数解得簇解c=√(?Q) ? α₁，A₀=0，A₁=?( ? α₁ )√( ?( (ρ²?1)( ρ²?²??²+2 ??1 )( α₁²+Q )P )/( ρ²?²?2 ρ ?²??²+ρ ?+2 ??1 ) )。由此得Jacobi椭圆函数解V_1,1(x,t)=?( ? SN(x?? t,m) α₁ )√( ( (ρ⁴?²?2 ρ²?²+2 ρ²??ρ²+?²?2 ?+1 )m²(m²?α₁²+1) )/( ρ²?²?2 ρ ?²+ρ ???²+2 ??1 ) )×exp(i(?(α₁c²/?)t+α₁x+Ω))，及相应P_1,1(x,t)、Q_1,1(x,t)表达式，表明在F-展开法下可获得含SN(ξ,m)的精确波解，当m→1或m→0时可退化为双曲或三角函数解，包括亮孤子与周期解。

5. Energy Balance Method（原文述及将能量平衡法用于进一步分析所获解的性质）：研究人员通过能量平衡法考察所得非线性波解的振荡传播属性，分析周期地震波在参数变化下的能量分布与局域化特征，以补充解析解构造之外的定性认识。

6. Graphical Representations（原文述及给出结果波结构的图形可视化）：研究人员通过图形展示所得波结构（亮孤子、暗型孤子、W型剖面、周期行波）的传播特性与结构性质，以直观呈现不同参数区间下波幅、宽度及形态变化。

讨论与结论部分总结：研究人员通过对耦合非线性Klein–Gordon系统的解析与定性综合分析，系统应用F-展开法、完全判别系统法与能量平衡法，推导出多种精确波解，包括Jacobi椭圆函数解及其退化形式的亮孤子、暗型孤子、W型剖面和周期行波解，揭示色散弹性介质中非线性地震波的丰富动力学行为；参数分析表明波速、波幅、局域化特性与Lamé系数λ、μ、密度?、波数α₁、行波速度?及椭圆模数m等密切相关；能量平衡分析提供周期波振荡传播的定性认识；研究结果丰富了非线性Klein–Gordon型地震模型的解析理论，可为非均匀弹性环境中非线性波传播现象的进一步研究提供有用框架，对地球物理中地震波建模与解释具有理论与应用意义。