2_反演部分_公式推导.docx

2.1反演公式推导由于所用模型的参数只考虑模型的速度值，所以可用由模型的速度参数所组成的列矩阵M来表示模型(矩阵M共有n个元素，并用表示矩阵的第i个元素)，由在实际中观测到的数据加上在未设置接收装置的网格点上的数据在角频率为时的频谱值可按模型参数的排列方式组成列矩阵，建立数学模型后，可表示在此模型上的模拟值P与M间的关系为：反演的过程就是使速度模型逐步逼近实际模型、理论数据逐步逼近观测数据的过程，我们可以建立目标函数来衡量它的逼近程度建立目标函数： ，其中H指矩阵的共轭转置其中T指矩阵的转置，*指矩阵的复共轭故可得：其中J是雅可比矩阵，且有，指雅可比矩阵第i行第j列的元素,J可表示成如下形式：由正演过程满足：AP = S对上述方程对求导，由于震源S和M无关，可得可得故式(2-3)可表示为由(2-4)可得，的第i个行元素为其中是单位矩阵，由于可表示为指第j个元素为1其余元素皆为0的列矩阵，而[]指按内部元素的坐标排列所形成的矩阵故上式可得由(2-4)式代入(2-2)式可得由于原来求正演时的方法可得此时A是对称的故(2-6)式可化为(可能不是对称)对目标在当前模型附近按泰勒展开，并忽略二阶以上的高阶项可得上式中H是黑塞矩阵，且故可得对(2-7)式对求导可得由(2-8)式可知H是一个对称的矩阵，故可得若在处可取得极小值，则得上式为0，可得由于在(2-9)式中第一项远大于第二项，故可忽略掉第二项，从而得到考滤到黑塞矩阵通常是病态的或奇异的，故需要对它进行一些正则化处理，最简单的正则化方式是在近似黑塞后面加上一个衰减项，用此方法可得由于近似黑塞矩阵是一个对角占优的矩阵，所以可以简化成只取对角元来代替，由于又进行了黑塞矩阵近似故应添加一个因子(可看作是一个步长)，也即