地震勘探1-理论.doc

第一章 地震勘探的基本理论 第一节 地震波在弹性介质中的传播规律 §1.1弹性介质中的基本波 （一）弹性波控制方程 从固体弹性理论可知，在均匀、各向同性的理想弹性介质中，三维波动方程可以用矢量表示为： （1-1-1） 式中向量u表示介质质点受外力F作用后的位移，称位移向量；向量F作用的外力，称为向量；常数、是介质的弹性常数，称拉梅(Lame)常数；常量是介质的密度；标量称体变系数，亦可表示为。 标符为拉普拉斯(Laplace)算子 如果对作用外力分别取散度(div)和旋度(rot)，则式（1-1-1）可分别写成 (1-1-2) (1-1-3) 式中：ω=rotu。 这说明如果对这种介质分别作用胀缩外力divF和旋转外力rotF的话，则在介质中分别存在二种扰动，胀缩力作用下产生由体变系数决定的介质体积相对胀缩的扰动，这就是纵波；在旋转外力作用下，则产生由向量ω决定的角度转动的扰动，这就是横波。这二种独立的扰动，分别以速度Vp和Vs传播。 和 (1-1-4) 如同重力场可用重力位，电场可用电位来描述一样，地震波场亦可用质点位移的位移位来描述。 根据亥姆霍兹（Helmholtz）涡流理论，任何一个矢量场，如果在定义域内有散度和旋度，则该矢量场可以用一个标量位的梯度场和一个矢量位的旋度场之和表示。则上式中的位移矢量u力矢量F可分别用位函数表示为 u=up+us=gradφ+rotψ F=Fp+Fs=gradΦ+rotΨ （1-1-5） 式中φ代表位移场的标量位；ψ代表位移场面的矢量位；Φ代表队标量力位；Ψ代表矢量力位。 将式（1-1-5）代入式（1-1-2）和（1-1-3）中，得到用位移位表示的波动方程： （1-1-6） （1-1-7） 分别称为纵波波动方程和横波波动方程。 在各种形式的弹性波中，平面波是一种特别简单的波。这种波以波面的形式在介质中传播，即在垂直于波传播方向的任一平面上，各点的振动是同相的。平面简谐波是波函数为简谐形式的平面波。因此，在研究波的传播问题时经常使用简谐波假定。 平面波波动方程的达朗贝尔（D`Alembert）解为： f（x,y,z,t）=f1（lx+my+nz-ct）+f2（lx+my+nz+ct） (1-1-8) 其中为波的传播方向，c为波速，f1表示沿正方向传播的平面简谐波，f2表示沿负方向传播的平面简谐波。 （二）弹性波的能量 地震波的传播实质是能量的传播。根据一般波动理论可知，波在介质中传播时的能量等于动能Eτ和位能Ep之和。假设波通过的介质体积为W，介质的密度为，对于简谐振动来说，波的振荡E可用下式表示： E=Eτ+Ep?A2f2W （1-1-9） 式中A表示波动的振幅，f表示波的频率。 上式说明波的能量和振幅平方、频率平方及介质的密度成正比。于是包含在介质中单位体积内的能量（称能量密度）亦应正比于振幅平方，即： ?=E/W?A2f2 （1-1-10） 定义单位时间通过单位面积的能量为波的能通量密度或波的强度I，因为实际地震勘探是在波前面的单位面积上观测波的能量信息的，如果时间dt内通过面积ds的能量为??V?dt?ds，则波的强度I为 （1-1-11） 式中V为速度。因此波的强度I正比于波的振幅平方。 （三）球面波 地震勘探中一般使用炸药振源，介质中发生的弹性振动从中心向四周传播。在均匀各项同性介质中这种波动过程具有中心对称性质，波面为球面，称为球面波，在球坐标系下，波动方程可以描述为： （1-1-12） 其达朗贝尔解为： （1-1-13） 其中第一项表示由中心向四周扩展的波，而第二项表示由无限远向中心汇集的波。当r、t固定，则复合变量（t—r/c）为一常数，波函数为定值。这样，在瞬间以为半径的球面上波场值相同，该球面为等相位面。波函数前的系数1/r表示波远离震源向外传播，其振幅不断衰减，且与到震源的距离成反比。事实上，随增大，波前面的面积越来越大，与r2成正比。从中心点震源所产生的波所据有的能量是一定的。在波向外传播的过程中，通过单位面积的能量将与r2成正比。又由于能量与振