地球物理场论实习报告.docx

声波方程数值模拟实验报告姓名：专业年级：学号：授课教师：小模型大模型程序分析拓展总分中国海洋大学实验环境：编程软件：MATLAB R2010a C-free 5.0软件环境：Microsoft Windows 7 Ultimate，32-bit硬件环境：CPU：Intel(R) Core(TM) i3 M350 @2.27GHz 内存：2.00G 硬盘：500G实验一实验题目：声波方程数值模拟———小模型波前快照实验目的：认识地震子波，以雷克子波为例，对地震波波形有直观认识。认识波的传播和反射等。利用数值算法模拟波前快照。实验要求：1.应用声波方程作为正演模拟的波动方程；2.将所提供震源函数离散后绘图； 3.给定小模型的二维速度-深度模型，绘出图形来4.在小模型中，整个区域的速度值可设为常数，即只有一种介质，将震源点放在模型中间，分别记录两个时刻的波前快照（即该时刻区域内所有网格点的波场值）。第一时刻为地震波还未传播到边界上的某时刻，第二时刻为地震波已经传播到边界上的某时刻，体会其人工边界反射；实验原理：地震波场模拟即地震正演，是指已知模型结构，通过物理或数值计算的方法模拟该地质结构下的地震波的传播，最终合成地震记录，也可以认为其是野外数据采集过程的室内再现。物理模拟花费昂贵，人们一般采用比较经济的数值模拟技术。地震波场数值模拟是在给定数学模型（如弹性波方程，声波方程等）、震源和地下几何界面、物性参数（岩层密度、速度等）情况下，研究弹性波或声波的传播规律。在本次试验中，选用一个小模型，即200*200的范围，视为完全弹性模型，为一层介质，在介质中取速度为定值，本实验中取为2500，将震源放在模型中央，分别记录两个时刻的波前快照，（即区域内所有网格点的波场值）。第一时刻为地震波还未传播到边界上的某时刻，第二时刻为地震波已经传播到边界上的某时刻，体会其人工边界反射。地震波场数值模拟是在给定数学模型（如弹性波方程，声波方程等）、震源和地下几何界面、物性参数（岩层密度、速度等）情况下，研究弹性波或声波的传播规律。地震波波动方程数值模拟方法主要包括克希霍夫积分法、傅里叶变换法、有限元法和有限差分法等。相对于上述几种方法，有限差分法是一种更为快速有效的方法。虽然其精度比不上有限元法，但因其具有计算速度快，占用内存较小的优点，在地震学界受到广泛的重视与应用。实验内容：波动方程对于二维速度-深度模型，地下介质中地震波的传播规律可以近似地用声波方程描述：（1）其中v(x,z)是介质在点(x , z)处的纵波速度，本实验中v为常数2500，为描述速度位或者压力的波场，为震源函数，本实验选用雷克子波（零相位子波）上式中，t为时间，fm为中心频率，一般取为20-40HZ，本实验取20，为控制频带宽度的参数，一般取3-5,这里取4，在实际计算过程中，需把此震源函数离散，参与波场计算。由图像可知子波能量主要集中在前部，大致上t35的时候S(t)都为零值，故在绘制S(t)图像时k取100，来比较清晰得呈现雷克子波。本实验不涉及吸收边界条件的使用。这样声波方程数值模拟所需的（1）震源函数（2）地层速度（波速）（3）边界条件均定义完全。为求（1）式的数值解，必须将此式离散化（包括时间离散和空间离散，这里dt取0.002，dh取4，即用有限差分来逼近导数，用差商代替微商。为此，先把空间模型网格化，得到横向上x=i*dh,纵向上z=j*dh,t=n*dt，k时刻(i,j)点的波场值为。实验中波场本应定义三维数组，但考虑到内存问题，将其定义为二维数组，采用循环交换赋值实现三维数组的功能，这里利用泰勒展式，展到二阶略去高阶小量，整理得(i,j)点k时刻的二阶时间微商，对于空间微分，采用四阶精度差分格式分别在(i,j)点k时刻展开到四阶小量，消除四阶小量并解出二阶微分，由时间和空间的二阶微商整理得到波场值。}（2）一般来讲，差分时差分格式阶数越高，得到的波场值精度越高，但稳定性会受影响，差分格式的稳定性不仅与差分格式本身有关，也随着差分精度的不同而不同，而且与网格步长之比的大小有关，值得指出的是差分格式的稳定性与微分方程无关，本实验所采用的（2）式为条件稳定，稳定条件：，（2）式由差分方程近似替代微分方程所得，所以如果空间和时间采样间隔不当，就会产生频散现象，导致波形畸变，甚至派生出多个同相轴，为减少频散，Δh需满足Dablain的经验公式：为Nyquist频率，一般取为主频的两倍，G为每个波长所占的网格点数，时间、空间为两阶差分的情况G取8，而时间、空间为四阶差分的情况G取4。则由上面两式得v∈[1000，3061.8]二、震源函数离散后的图像及速度-深度模型由由上面定义的数值，用Matlab将震源函数成图，下面是成图的源程序：clf;f=20