泰州网站设计,郑州网络运营平台有哪些,铜仁市住房和城乡建设部网站,如何做品牌营销在三维空间中#xff0c;弹性波方程#xff08;Navier 方程#xff09;是一个矢量方程。施图姆–刘维尔#xff08;S-L#xff09;理论的应用#xff0c;通常是通过分离变量法#xff0c;将复杂的矢量场分解为一系列互相正交的标量场或矢量基函数的过程。 在地震学中&am…在三维空间中弹性波方程Navier 方程是一个矢量方程。施图姆–刘维尔S-L理论的应用通常是通过分离变量法将复杂的矢量场分解为一系列互相正交的标量场或矢量基函数的过程。在地震学中这最典型的应用是处理球对称地球模型如 PREM 模型的简正振型Normal Modes或分层介质中的波传播。弹性波方程的算子化各向同性均匀介质中的三维弹性波方程为ρ∂2u∂t2(λμ)∇(∇⋅u)μ∇2u\rho \frac{\partial^2 \mathbf{u}}{\partial t^2} (\lambda \mu) \nabla (\nabla \cdot \mathbf{u}) \mu \nabla^2 \mathbf{u}ρ∂t2∂2u​(λμ)∇(∇⋅u)μ∇2u其中u(x,t)\mathbf{u}(\mathbf{x}, t)u(x,t)是位移矢量。为了应用 S-L 理论我们考虑定态解简谐波令u(x,t)U(x)e−iωt\mathbf{u}(\mathbf{x}, t) \mathbf{U}(\mathbf{x}) e^{-i\omega t}u(x,t)U(x)e−iωt得到空间部分的特征方程−ω2ρU(λμ)∇(∇⋅U)μ∇2U- \omega^2 \rho \mathbf{U} (\lambda \mu) \nabla (\nabla \cdot \mathbf{U}) \mu \nabla^2 \mathbf{U}−ω2ρU(λμ)∇(∇⋅U)μ∇2U这可以写成线性算子形式LUω2U\mathcal{L} \mathbf{U} \omega^2 \mathbf{U}LUω2U。这里的ω2\omega^2ω2就对应 S-L 问题中的本征值。分离变量与标量化对于三维问题我们通常利用**亥姆霍兹分解Helmholtz Decomposition**将位移矢量U\mathbf{U}U分解为标量位ϕ\phiϕ纵波 P和矢量位Ψ\boldsymbol{\Psi}Ψ横波 SU∇ϕ∇×Ψ,∇⋅Ψ0\mathbf{U} \nabla \phi \nabla \times \boldsymbol{\Psi}, \quad \nabla \cdot \boldsymbol{\Psi} 0U∇ϕ∇×Ψ,∇⋅Ψ0以球坐标系(r,θ,ϕ)(r, \theta, \phi)(r,θ,ϕ)为例对于球对称介质位移场可以分解为类扭转振荡Toroidal modes只有横波分量对应SHSHSH波。类球振荡Spheroidal modes包含纵波和垂直横波分量对应P−SVP-SVP−SV波。径向方程标准的 S-L 问题以类扭转振荡为例其位移场可以表示为UTW(r)[r^×∇Ylm(θ,ϕ)]\mathbf{U}_T W(r) [ \hat{\mathbf{r}} \times \nabla Y_{lm}(\theta, \phi) ]UT​W(r)[r^×∇Ylm​(θ,ϕ)]其中YlmY_{lm}Ylm​是球谐函数。将此形式代入波动方程角度部分被球谐函数抵消剩下的径向部分W(r)W(r)W(r)满足如下方程ddr[μr4dWdr][ω2ρr4−(l−1)(l2)μr2]W0\frac{d}{dr} \left[ \mu r^4 \frac{dW}{dr} \right] \left[ \omega^2 \rho r^4 - (l-1)(l2)\mu r^2 \right] W 0drd​[μr4drdW​][ω2ρr4−(l−1)(l2)μr2]W0如果我们对比标准的 S-L 形式ddr[p(r)dWdr][q(r)λρ(r)]W0\frac{d}{dr} [p(r) \frac{dW}{dr}] [q(r) \lambda \rho(r)] W 0drd​[p(r)drdW​][q(r)λρ(r)]W0可以清晰地看到本征值λω2\lambda \omega^2λω2频率的平方。系数函数p(r)μ(r)r4p(r) \mu(r) r^4p(r)μ(r)r4。权重函数ρSL(r)ρ(r)r4\rho_{SL}(r) \rho(r) r^4ρSL​(r)ρ(r)r4。势函数项q(r)−(l−1)(l2)μ(r)r2q(r) -(l-1)(l2)\mu(r) r^2q(r)−(l−1)(l2)μ(r)r2。边界条件与奇异性在地球物理应用中S-L 问题必须配合边界条件才能定解自由表面条件在地球表面rRrRrR应力为零即Trμ(dWdr−Wr)0T_r \mu (\frac{dW}{dr} - \frac{W}{r}) 0Tr​μ(drdW​−rW​)0。内部连续性在介质分界面处如核幔边界位移和应力必须连续。地心奇点在r0r0r0处p(0)0p(0)0p(0)0这正是你图片中提到的奇异边界点。根据 S-L 理论我们要求解在r0r0r0处保持有限。为什么这样做应用价值正交性关系S-L 理论保证了不同频率的振型Un\mathbf{U}_nUn​和Um\mathbf{U}_mUm​在全地球体积VVV上关于质量密度ρ\rhoρ正交∭Vρ(x)Un∗⋅UmdVδnm\iiint_V \rho(\mathbf{x}) \mathbf{U}_n^* \cdot \mathbf{U}_m dV \delta_{nm}∭V​ρ(x)Un∗​⋅Um​dVδnm​这意味着我们可以把任何地震位移记录分解为这些固有振型的叠加。合成地震图Synthetic Seismograms在 CPS 等软件中计算理论地震图的核心逻辑就是先通过求解 S-L 方程组得到本征值色散曲线和本征函数振型深度分布然后进行模态叠加Mode Summation。反演地球内部结构如果我们观测到了地球的固有频率ωobs\omega_{obs}ωobs​我们可以根据 S-L 方程建立ω\omegaω对参数μ(r),ρ(r)\mu(r), \rho(r)μ(r),ρ(r)的敏感核从而反推地球深部的物理性质。