这个推导过程是变分法（Calculus of Variations）在连续介质力学中的经典应用。作为地球物理专业的学生，理解这个推导能帮你从能量的角度直观把握波的产生。
我们从作用量泛函 III 的变分为零（δI=0\delta I = 0δI=0）开始，逐步展开推导。
- 建立变分表达式
首先，对作用量积分 III 施加变分算子 δ\deltaδ。由于积分号和变分算子可以互换，我们直接对被积函数（拉格朗日密度）取变分：
δI=∫t0t1∫0lδ(12ρut2−12μux2)dxdt=0\delta I = \int_{t_0}^{t_1} \int_0^l \delta \left( \frac{1}{2} \rho u_t^2 - \frac{1}{2} \mu u_x^2 \right) dx dt = 0δI=∫t0t1∫0lδ(21ρut2−21μux2)dxdt=0
根据全微分规则，对 utu_tut（对时间的导数）和 uxu_xux（对空间的导数）分别求偏导：
δI=∫t0t1∫0l(ρutδut−μuxδux)dxdt=0\delta I = \int_{t_0}^{t_1} \int_0^l \left( \rho u_t \delta u_t - \mu u_x \delta u_x \right) dx dt = 0δI=∫t0t1∫0l(ρutδut−μuxδux)dxdt=0
注意这里的关系：δut=∂∂t(δu)\delta u_t = \frac{\partial}{\partial t}(\delta u)δut=∂t∂(δu)，δux=∂∂x(δu)\delta u_x = \frac{\partial}{\partial x}(\delta u)δux=∂x∂(δu)。
- 利用分部积分法（Integration by Parts）
为了从积分中提取出关于 δu\delta uδu 的项，我们需要处理掉 δut\delta u_tδut 和 δux\delta u_xδux 中的导数符号。
第一项：处理时间导数项
对第一项关于 ttt 进行分部积分：
∫t0t1ρut∂∂t(δu)dt=ρutδut0t1⏟Boundary Term 1−∫t0t1∂∂t(ρut)δudt\int_{t_0}^{t_1} \rho u_t \frac{\partial}{\partial t}(\delta u) dt = \underbrace{\left \\rho u_t \\delta u \\right{t_0}^{t_1}}{\text{Boundary Term 1}} - \int_{t_0}^{t_1} \frac{\partial}{\partial t}(\rho u_t) \delta u dt∫t0t1ρut∂t∂(δu)dt=Boundary Term 1 ρutδut0t1−∫t0t1∂t∂(ρut)δudt
第二项：处理空间导数项
对第二项关于 xxx 进行分部积分：
∫0l−μux∂∂x(δu)dx=−μuxδu0l⏟Boundary Term 2+∫0l∂∂x(μux)δudx\int_0^l -\mu u_x \frac{\partial}{\partial x}(\delta u) dx = \underbrace{\left -\\mu u_x \\delta u \\right0^l}{\text{Boundary Term 2}} + \int_0^l \frac{\partial}{\partial x}(\mu u_x) \delta u dx∫0l−μux∂x∂(δu)dx=Boundary Term 2 −μuxδu0l+∫0l∂x∂(μux)δudx
- 处理边界条件
在变分原理中，我们通常假设路径在起始和结束点，以及空间的边界上是确定的：
时间边界：在 t0t_0t0 和 t1t_1t1 时刻，系统的状态 uuu 是给定的，因此变分 δu(x,t0)=δu(x,t1)=0\delta u(x, t_0) = \delta u(x, t_1) = 0δu(x,t0)=δu(x,t1)=0。这意味着 Boundary Term 1 消失。
空间边界：对于两端固定的弦（如图片提到的 u(0,t)=u(l,t)=0u(0,t)=u(l,t)=0u(0,t)=u(l,t)=0），其边界变分 δu(0,t)=δu(l,t)=0\delta u(0, t) = \delta u(l, t) = 0δu(0,t)=δu(l,t)=0。这意味着 Boundary Term 2 消失。
- 提取波动方程
将剩下的积分项重新合并，我们得到：
δI=∫t0t1∫0l(−ρutt+μuxx)δu dxdt=0\delta I = \int_{t_0}^{t_1} \int_0^l \left( -\rho u_{tt} + \mu u_{xx} \right) \delta u \, dx dt = 0δI=∫t0t1∫0l(−ρutt+μuxx)δudxdt=0
其中 utt=∂2u∂t2u_{tt} = \frac{\partial^2 u}{\partial t^2}utt=∂t2∂2u，uxx=∂2u∂x2u_{xx} = \frac{\partial^2 u}{\partial x^2}uxx=∂x2∂2u。
根据变分法基本引理（Fundamental Lemma of Calculus of Variations），由于 δu\delta uδu 是区间内任意的微小扰动，要使整个积分为零，括号内的项必须恒为零：
−ρutt+μuxx=0-\rho u_{tt} + \mu u_{xx} = 0−ρutt+μuxx=0
整理后即得标准的一维波动方程：
ρutt−μuxx=0\rho u_{tt} - \mu u_{xx} = 0ρutt−μuxx=0
物理总结
ρutt\rho u_{tt}ρutt：代表单位长度的惯性力（质量 ×\times× 加速度）。
μuxx\mu u_{xx}μuxx：代表单位长度的弹性恢复力（张力与曲率的乘积）。
这个推导展示了物理世界的一个深刻真理：波的传播本质上是介质在动能和势能之间不断交换并试图寻找作用量极小值过程的动力学表现。
既然你在研究泸州地区的诱发地震，需要我进一步解释这个推导如何扩展到三维各向异性介质（即地震勘探中更常用的位移方程）吗？