0_正演部分_公式推导.docx

1.1对于无吸收边界的正演方程由二维声波方程：式中为在模型中(x,z)点t时刻的声波波场值，为在模型中(x,z)点的介质速度，为在模型中(x,z)点的震源在t时刻的值分别记、所对应的傅里叶变换为、对式(1-1)进行傅里叶变换可得对速度模型进行网格离散化(且使在x,z方向的网格间距相等，均为)，在角频率为时用代表，用代表，用代表利用有限差分，在忽略高阶无限小的情况下可得由于任一点的波场值可用此点和它上下左右四个点的加权平均来近似，也即对式(1-3),(1-4)式加权相加并利用式(1-5)可得1.2对于添加PML吸收边界的正演方程为了减少人为边界所造成的边界反射，我们可以添加PML层来进行吸收，在PML层中由于没有震源，声波波场方程变为：也可表示成如下形式若令则存在的拆分使得均为x,z,t的函数、满足下式对(1-7)、(1-8)式添加衰减系数可得对(1-9)、(1-10)式进行傅叶变换可得(记、、、的傅里叶变换分别为、、、)，则可得令、则上面四式可化为对(1-11)、(1-12)式分别进行有限差分，忽略高阶无限小，在角频率下用代表用代表，代表可得利用及(1-13)、(1-14)、(1-15)、(1-16)式可得同理，对原坐标进行45度旋转后，可得类似于(1-17)式的方程其中指旋转了45度后的x,z方向对式(1-17)、(1-18)式加权相加，并利用式(1-5)可得在(1-19)式的基础上添加一个震源项，可得在正常区域中由于无衰减也即、、、的值均为1，式(1-20)可化成与式(1-6)完全相同的形式若相邻点的波场值变化不大，在一定误差(也即零阶近似下)允许的情况下可得：若在速度模型中速度模型的速度值在小区域中变化不大，在一定误差(也即零阶近似下)允许的情况下可得：于是(1－20)式可化为：对于旋转后的x,z方向的衰减值可以直接用x,z方向的衰减值来代替式(1-21)的整体可以表示成线性方程：AP = S由(1-21)式可以求得线性方程中矩阵A的具体形式矩阵A对速度模型中的任意网格点的速度值求偏导可得若没有使用速度值的零阶近似则可得上式中表示元素所在网格点在总网格点中的索引值，上式中没有说明在其为与PML层交界时的具体形式由矩阵A的具体形式可知对于矩阵，其中除了上面所列的五个元素外，所有其余的的元素值均为0