第六章弹性波波动方程及其解.pptVIP

§6.6 有源地震波的传播 构成半空间的理想弹性体的密度、拉梅系数已知。当t=0时，弹性体处理静止状态；当t0时，弹性体受到均匀压力P(t)的作用，不计体力。求t0时，弹性体内的位移场及应力场。 分析：此问题同样是个一维问题 u1=u2=0 u3=u3(x3,t) x1 x3 当t0时，弹性体在x3=0的表面上受力为t1=t2=0,t3=P(t)。 当t=0时，弹性体处于平衡、静止状态 由已知的位移场，即可求出应变张量场和应力张量场。 由上式可知，对应特定的时刻t来说，对于x3vt的 区域，弹性体还没有开始振动；对于x3≤vt的区域， 弹性体质点位移由上式给出。 无界弹性体中球腔振源产生的波 无界的理想弹性体内有如图所示的空球腔，半径为a，受力前处于静止状态。今突然在球腔壁上突然施加均匀压力p0H(t)，其中p0为常数，H(t)为阶跃函数。求在该弹性体内产生的位移场。 在这种情况下使用直角坐标系是否合适？ 此时球坐标系的方程如何表达？ 当t0, r=a时 将t=0时，位移和速度为“0”的初始条件代入 6.7 各向异性介质中平面波的传播 各向异性介质中地震波传播的一般表达式 Christoffel方程 各向异性介质中地震波传播的速度 Christoffel方程普遍成立的前提条件是其系数行列式为零 VTI介质的群速度与相速度 X1 X3 o θ f vpha vgroup 6.8 波动方程的数值解 各向同性波动方程 当不计体力, 只讨论地震波沿ox3轴传播时,方程可简化为: 只讨论纵波时, 有: 差分格式理论 时间一阶导数 时间二阶导数 差分格式理论 整理得 写成离散格式为: 已知条件 二维波动方程的有限差分解 二维二分量各向异性介质有限差分解 根据各向异性介质的虎克定律和运动微分方程求解，此时解由5个方程组成的方程组。 高精度地震波波动方程有限差分解 当接收的是ox2轴分量时的情况又怎样？ 此时接收的波场中只有横波分量，没有纵波分量！ i必须等于2 三维单垂向分量波动方程(3D1VC) 三维单垂向分量是目前常规三维地震勘探数据采集的观测方式； 即在地表平面上采集地震数据的观测方式； 取地表平面为ox1x2x3平面，所接收的分量为u3，平行ox3轴垂直向下。 此时，有： X1 X2 X3 O 上式即为三维单垂向波动方程。分析上式可知，接收点的响应的地震波场中即有纵波成份，同时也包含横波成份。 一维三分量波动方程 一维三分量是目前VSP测井中常规数据采集方式； 即在一维线上激发与接收，所接收的分量u=(u1,u2,u3) X1 X2 X3 O 平面波仅沿ox3轴传播 平面波沿ox2轴传播 平面波沿ox1轴传播 沿过原点任意射线传播时 当平面波的传播方向与ox1x2x3坐标系的任意一个坐标轴都不重合时，此时的平面波方程比较复杂。 首先，以波的传播方向为ox3轴建立一个新坐标系； 其次，就新坐标系内将平面波方程表达出来； 最后，通过坐标旋转关系求出在旧坐标系内的表达式。 §6.5 波动方程的一般解 平面波假设 平面波只是一种理想化模型。波前面离开波源足够远时，可以把波前面近似地看作平面，叫做平面波。 平面波分类 根据波函数对时间的依赖关系：脉冲型、简谐型； 根据振幅随场点坐标的变化：均匀平面波、非均均平面波； 一维波动方程的解 由于平面被的波前面是一系列互相平行的平面，因此，同一波前面上各点振动情况完全相同，所以对于沿某一方向传播的平面波，我们可以选择一个坐标系，使得波就象在其中某一个坐标轴方向上传播一样。 行波法解 令波传播方向沿ox1轴，则ui=ui(x1,t)； 则此时的波动方程为： 由达朗贝尔方程，得一维波动方法为： 上式即为一维波动方程的通解，f1,f2为任意函数。 上式意义 f1(x-vt)表示由波源出发，以速度v，沿ox轴方向传播的行波；? f2(x+vt)表示由无穷远处出发，以速度v，沿ox轴负方向传播的行波；? 坐标轴上某点处的波是由上述两种行波干涉波叠加的结果； 上述行波法只考虑了波的传播，没有考虑波的初始扰动和边界情况等因素；因此，所得解没有太多实用价值。 地震波传播的边界条件和初始条件 边界条件 给定了弹性体在其边界面上所满足的条件，即称之为边界条件。 位移边界条件：在弹性体的所有边界面上都给定了位移Su。此类问题称为第一类边值问题。 应力边界条件：在弹性体的所有边界面上都给定了表面力St。此类问题称为第二类边值问题。 混合边界条件：弹性体的边界面可以分成两部分，一部分给定了位移Su，另一部分给定了表面力St。此类问题称为第三