公式推导_反演部分.docxVIP

由于所用模型的参数只考虑速度大小，可用列矩阵M来表示模型的速度参数信息(共有n个，并用表示矩阵的第i个元素)，由在实际中观测到的数据加上由其它未观测处的值在角频率为时可按模型参数的排列方式组成列矩阵，由于建立成数学模型后，可表示在此模型上的模拟值P与M间的关系成：建立目标函数： ，其中H指矩阵的共轭转置其中T指矩阵的转置，*指矩阵的复共轭故可得：其中J是雅可比矩阵，且有，指雅可比矩阵第i行第j列的元素,J可表示成如下形式：由正演过程满足：AP = S对上述方程对求导，由于震源S和M无关，可得可得故由(4)可得，的第i个行元素为由于指第j个元素为1其余元素皆为0的列矩阵，而[]指按内部元素的坐标排列所形成的矩阵故上式可得由(4)式代入(2)式可得由于原来求正演时的方法可得此时A是对称的故(6)式可化为对目标在当前模型附近按泰勒展开，并忽略二阶以上的高阶项可得上式中H是黑塞矩阵，且故可得对(10)式对求导可得其中E是单位矩阵，由(8)可得H是对称的故若在此处可取得极小值，则得上式为0，可得由于在(10)式中第一项远大于第二项，故可得考滤到黑塞矩阵通常是病态的或奇异的，所以需要对它进行一些正则化处理，最简单的正则化方式是在近似黑塞后面加上一个衰减项，也即由于近似黑塞矩阵是一个对角占优的矩阵，所以可以简化成只取对角元来代替，并添加一个因子(可看作是一个步长)，可得现说明求近似黑塞矩阵的主对角元素的方法由上面的式(5)公式推导_正演部分的式(22)可得的第i个对角元素上式是对于单炮点的情形，如果对于多炮点由于模型比较大，则上式中j的可取范围就会很大，此时可以进行一些合理限定，如限定为所有的接收点所在的网格点，则上式可变成 (13-3)利用式(13-2)，可使(13-3)化为(13-4)上式中中的指在模型中下标为i的速度模型的值，而整体则表示对列矩阵进行处理而使新得的矩阵中只有第i个量保持不变，而使其余的值全部置0；是在设定的初始模型上正演的结果，为了使所选用的震源更接近实际震源，通常对其进行处理，也即添加了一个震源修正因子求梯度的方法由式(2)得由式(7)得公式推导_正演部分的式(22)式(14-2)可化为考虑到多炮点可得求步长因子的方法由上面的分析可知应该大于零，可用一定的方法确定两个值、使得与三者满足上式的下标用来标示所用的利用抛物线插值法即可求得使最小的值