声波方程正演模拟1复习过程

1声波方程数值模拟

――地球物理场论

基础Ⅰ期末(qīmò)作业（1）任课教师：宋鹏第一页，共45页。2一、地震勘探基本原理

二、波动方程类型(lèixíng)及其局限性

三、数值算法类型(lèixíng)及其局限性

四、声波方程的有限差分法数值模拟内容提纲第二页，共45页。3一、地震(dìzhèn)勘探基本原理＊▽▽▽▽▽xt第三页，共45页。4 同相轴为双曲线，即反射(fǎnshè)波的时距曲线为双曲线，反射(fǎnshè)波一个同相轴可带来一个地层的信息。实际地下介质非常复杂，所得到的炮集记录也包含更多的地下信息。实际的炮集记录见图1-1和1-2。第四页，共45页。5图1-1陆上(lùshànɡ)某区实际地震记录第五页，共45页。6图1-2海上某区实际地震(dìzhèn)记录第六页，共45页。7广义的地震反演即是从地震炮集记录出发，经过复杂的去噪、速度分析以及偏移成像处理等手段(shǒuduàn)得到反映地下的地质结构的地震剖面。实际的地震剖面见图1-3和1-4。第七页，共45页。8图1-3 陆上(lùshànɡ)某区地震剖面第八页，共45页。9图1-4海上某区地震(dìzhèn)剖面第九页，共45页。10地震波场模拟即地震正演，是指已知模型结构，通过物理或数值计算的方法模拟该地质结构下的地震波的传播，最终合成地震记录，也可以认为其是野外数据采集过程的室内再现。物理模拟花费昂贵，人们(rénmen)一般采用比较经济的数值模拟技术。地震波场数值模拟是在给定数学模型（如弹性波方程，声波方程等）、震源和地下几何界面、物性参数（岩层密度、速度等）情况下，研究弹性波或声波的传播规律。第十页，共45页。11 一、地震勘探基本原理

二、波动(bōdòng)方程类型及其局限性

三、数值算法类型及其局限性

四、声波方程的有限差分法数值模拟内容提纲第十一页，共45页。12二、波动(bōdòng)方程类型及其局限性1、声波(shēnɡbō)方程：一阶压力(yālì)-速度方程组：二阶标量声波方程：第十二页，共45页。13二、波动(bōdòng)方程类型及其局限性能够描述且只能(zhīnénɡ)描述纵波的传播规律，包括直达波、反射波、透射波、折射波等，但不能描述转换波传播规律。需要的已知条件包括：1）震源函数2）地层速度(sùdù)/密度3）边界条件第十三页，共45页。142、弹性(tánxìng)波方程：第十四页，共45页。15能够描述纵、横波的传播(chuánbō)规律，包括直达波、反射波、透射波、折射波以及转换波等。需要的已知条件包括：1）震源函数(hánshù)2）地层速度或根据方程的类型需要提供的地层的其它弹性参数3）边界条件第十五页，共45页。163、粘声波/弹性波方程 前面(qiánmian)讨论的是理想弹性介质，波在其中传播时，没有能量的损耗，介质中应力和应变关系严格遵循胡克定律（这种理想介质称虎克固体），但波在实际介质中传播时，是有能量损耗的，这就是所谓的弹性波吸收。波在传播过程中，实际介质的不同部位之间会出现某种摩擦力，称为内摩擦力或粘滞力。这种力导致机械能向其他形式能量转换，最终转化为热能消耗掉。第十六页，共45页。17 在地震勘探中，地震波传播的实际介质是十分复杂的。在一定条件下，即震源作用时间短，作用力微小，地球介质可以看作完全弹性模型，但随着地震勘探技术的发展，勘探精度要求提高，面临复杂地质目标时，要求地震勘探采用更加符合实际的介质模型进行研究。粘弹性介质模型更符合实际。 但是到目前为止，在地震资料反演处理中应用(yìngyòng)最多的还是声波方程，弹性波以及粘弹性波方程的应用(yìngyòng)还只是停留在模拟层次上。第十七页，共45页。18 一、地震(dìzhèn)勘探基本原理

二、波动方程类型及其局限性

三、数值算法类型及其局限性

四、声波方程的有限差分法数值模拟内容提纲第十八页，共45页。19三、数值算法(suànfǎ)类型及其优缺点地震波波动方程数值模拟方法主要包括克希霍夫积分法、傅里叶变换法、有限元法和有限差分法等。克希霍夫积分法引入射线追踪过程(guòchéng)，本质上是波动方程积分解的一个数值计算，在某种程度上相当于绕射叠加。该方法计算速度较快，但由于射线追踪中存在着诸如焦散、多重路径等问题，故其一般只能适合于较简单的模型，难以模拟复杂地层的波场信息。第十九页，共45页。20傅里叶变换法是利用空间的全部信息对波场函数进行三角函数插值，能更加精确地模拟地震波的传播规律，同时(tóngshí)，利用快速傅里叶变换(FFT)进行计算，还可以提高运算效率，其主要优点是精度高，占用内存小，但缺点是计算速度较慢，对模型的适用性差，尤其是不适应于速度横向变化剧烈的模型。第二十页，共45页。21波动方程有限元法的做法是：将变分法用于单元分析，得到单元矩阵，然后将单元矩阵总体求和得到总体矩阵，最后求解总体矩阵得到波动方程的数值解；其主要优点是理论上可适宜于任意地质体形态的模型，保证复杂地层形态模拟的逼真性，达到很高的计算精度，但有限元法的主要问题(wèntí)是占用内存和运算量均较大，不适用于大规模模拟，因此该方法在地震波勘探中尚未得到广泛地应用。第二十一页，共45页。22相对于上述几种方法(fāngfǎ)，有限差分法是一种更为快速有效的方法(fāngfǎ)。虽然其精度比不上有限元法，但因其具有计算速度快，占用内存较小的优点，在地震学界受到广泛的重视与应用。第二十二页，共45页。23 一、地震勘探(kāntàn)基本原理

二、波动方程类型及其局限性

三、数值算法类型及其局限性

四、声波方程的有限差分法数值模拟内容提纲第二十三页，共45页。24四、声波方程(fāngchéng)的有限差分法数值模拟对于二维速度-深度模型，地下介质中地震波的传播规律可以近似地用声波方程(fāngchéng)描述：是介质在点（x,z）处的纵波(zònɡbō)速度，为描述速度位或者压力的波场，为震源函数。(4-1)第二十四页，共45页。25空间(kōngjiān)模型网格化（如图4-1所示）：图4-1差分(chàfēn)网格划分示意图第二十五页，共45页。26网格间隔(jiàngé)长度，时间(shíjiān)采样步长表示(biǎoshì)(i,j)点k时刻的波场值第二十六页，共45页。27 时间二阶、空间二阶差分格式(géshi)推导如下：将在(i,j)点k时刻(shíkè)用Taylor展式展开：(4-2)(4-3)第二十七页，共45页。28将上两式相加，略去(lüèqù)高阶小量，整理得(i,j)点k时刻的二阶时间微商为：(4-4)第二十八页，共45页。29同理可得(i,j)点k时刻的二阶空间(kōngjiān)微商分别为：(4-5)(4-6)第二十九页，共45页。30这就实现了用网个点波场值的差商代替(dàitì)了偏微分方程的微商，将上三个式子代入(4-1)式中得：(4-7)式中为介质速度的空间(kōngjiān)离散值是空间(kōngjiān)离散步长为时间采样率第三十页，共45页。31同理可得空间(kōngjiān)四阶精度的差分格式为：第三十一页，共45页。32

同理可得时间二阶、空间四阶精度的声波(shēnɡbō)方程差分格式为：(4-8)第三十二页，共45页。33为震源(zhènyuán)函数，一般使用一个理论(lǐlùn)的雷克型子波代替，即：为时间(shíjiān)为中心频率，一般取为20-40HZ为控制频带宽度的参数，一般取2-5确定震源位置第三十三页，共45页。344.1稳定性条件(tiáojiàn)对于特定的偏微分方程只有特定的几种有限(yǒuxiàn)差分格式是无条件或有条件稳定的，（4-7）、（4-8）式即是已被证明的有条件稳定格式，其稳定性条件分别为：(4-9)(4-10)第三十四页，共45页。354.2频散关系式同时，在差分计算过程中，如果空间和时间采样间隔不当，就会导致波形畸变，甚至派生出多个同相轴，这种现象称为频散现象。偏微分方程本身没有频散，网格(wǎnɡɡé)频散是由于差分方程近似替代微分方程引起的。当波场按照波动方程所表示的微分方程传播时，波场的传播速度就是波动方程中的速度，但当波场按照波动方程离散化后的差分方程传播时，波场的传播速度就不再是波动方程中的速度了，而是与波的频率和波数有关的函数，具有不同频率和波数的波有不同的传播速度，因而在传播过程中会出现频散，发生畸变，且随走时的增加而增加。第三十五页，共45页。36Dablain给出了一个能有效减少(jiǎnshǎo)频散的经验公式为：为Nyquist频率(pínlǜ)，一般取为主频的两倍G为每个波长所占的网格点数，时间、空间为两阶(liǎnɡjiē)差分的情况G取8，而时间、空间为四阶差分的情况G取4。(4-9)第三十六页，共45页。37当所给震源函数、空间(kōngjiān)网格间隔、时间采样间隔以及地下介质的波速满足稳定性条件及频散关系式时，就可以应用(4-7)式递推求得所给的速度—深度模型内任意时刻、任何采样点的波场值，并可最终得到一个人工合成的地震记录。第三十七页，共45页。38*4.3 边界条件在地震波场正演模拟中，必须引入人工边界来界定计算区域。人工边界若不做特殊处理，就会随着(suízhe)波场的递推计算在边界上产生虚假反射波从而扰乱波场，人工边界的处理是地震波场正演数值模拟的一个重要课题。本次作业不涉及边界条件的使用，可通过增大模型来避免边界反射干扰，有精力的同学可通过查阅资料获得(huòdé)解决边界问题的方法。第三十八页，共45页。39本次大作业的具体要求为：1、应用声波方程作为正演模拟的波动方程；2、将所提供震源函数离散后绘图；3、给定两个二维速度-深度模型（一个小模型；一个大模型），绘出图形来；4、对于小模型，整个区域的速度值可设为常数，即只有一种介质，将震源点放在模型中间(zhōngjiān)，分别记录两个时刻的波前快照（即该时刻区域内所有网格点的波场值）。第一时刻为地震波还未传播到边界上的某时刻，第二时刻为地震波已经传播到边界上的某时刻，体会其人工边界反射；第三十九页，共45页。405、对于大模型，定义为水平层状速度模型（至少两层）；做两个实验，一是将震源点放在区域表层任一点，记录下某些时刻的波前快照，体会地震波在两种介质的分界面上传播规律；二是合成一个(yīɡè)地震记录，即记录下与震源同一深度点的各点所有时刻的波场值，并指出记录上的同向轴分别对应哪些波？6、要附有实验报告，格式不固定，