一种多参数全波形反演策略.docxVIP

一种多参数全波形反演策略

1参数化方式的影响

fwi是近30年开发的高分辨率地震数据处理方法。利用地震观测数据集的波形信息提取地下媒体的弹性参数（如速度、密度、弹性模量等）。全波形反演不仅可以精确地偏移成像（刘国峰等，2012），而且可以应用于存储层的定量描述领域（单瑞等，2011；李国平等，2011）。这在油气勘探和开发中发挥了重要作用。

由于全波形反演利用了地震波场的全部信息，理论上来讲，应该尽可能多的考虑影响波传播的因素，特别是将大孔径、宽方位数据用于全波形反演时，AVO/AVA效应和方位各向异性效应等需要特别考虑，这就要求将声波近似单参数全波形反演扩展到多参数全波形反演，包括纵横波速度、密度、吸收衰减效应、各向异性参数等．

然而，各参数之间的相互耦合，增加了多参数全波形反演的非线性程度．分析不同参数之间的相互关系，研究不同参数化方式对反演的影响程度，是建立多参数反演策略的有效途径．Tarantola(1986）和Mora(1987）指出，在P波震源反射地震波形反演中，弹性模量是一种比较差的参数化方式；对近偏移距数据（短波长成分）来说，阻抗是相对较好的参数化方式，而对远偏移距数据（长波长成分）来说，速度是比较好的参数化方式，但是，不管利用哪种数据信息，密度都不能被很好地估计．Kohn等（2012）指出，在弹性波多参数波形反演中，不管利用哪种参数化方式，弹性模量、纵横波速度、阻抗都可以被较好地反演出来，而密度反演结果较差；但总的来说，利用速度、阻抗参数化方式得到的反演结果要优于模量参数化方式得到的反演结果．Forgues和Lambaré（1997）通过辐射模式分析指出，密度很难利用全波形反演得到．因此，大部分学者在反演过程中，将密度固定成某一常数或利用速度通过经验关系换算得到密度．

密度是岩石物性参数的重要组成部分，如果能通过反演同时得到比较可靠的纵横波速度、密度资料，这对储层评价、岩性解释和油藏描述等都具有重要作用．弹性波多参数全波形反演的研究实例证明，虽然可以得到相对可靠的纵横波速度，但密度并不能被很好地估计（Choietal.,2008）．因此，需要对密度反演进行深入研究，提出有利于密度反演的多参数全波形反演策略．Jeong等（2012）提出了一种频率域密度反演策略，首先将密度固定为任意常数，通过单参数反演得到弹性模量，进而将换算得到的纵横波速度作为初始模型，再同时反演速度和密度，有效提高了密度的反演精度．Prieux等（2013）先利用大偏移距数据单参数反演速度参数，再利用全偏移距数据同时反演速度、密度参数，有效提高了密度反演的稳定性，但密度反演结果仍然不够理想．他们同时指出，利用速度-密度参数化方式反演得到的速度、密度和阻抗结果要好于利用速度-阻抗参数化方式得到的反演结果．

本文利用速度-密度参数化方式，通过辐射模式及目标函数敏感度分析，研究了速度、密度在反演过程中的相互影响，提出了一种有利于速度、密度同时反演的策略，并通过球状异常模型和Marmousi-2模型的数值实验验证了该策略的有效性．

2非负局部约束的非负局部约束

全波形反演通过拟合模拟数据umod和观测数据uobs，使波场残差δu达到最小，进而获取地下介质的物性参数m．目标函数可以表达为

可以利用拟牛顿法更新模型参数mn，使目标函数E达到最小：

其中，为目标Δ函数E对模型参数m的梯度，α为更新步长，为近似Hessian．常用的拟牛顿方法为有限内存BFGS(L-BFGS）法，它利用一定数量的梯度残差向量和模型残差向量来迭代构造正定矩阵Bk，并将它作为公式（2）中Hessian逆矩阵Hk-1的近似．矩阵更新公式为（NocedalandWright,2006):

选定某个初始矩阵Bk0，对（3）式做m次递推可得：

因而，只需记录m（一般取3～20，本文取m=5）个向量对{si,yi}，i=k-m，…k-1，就可以按公式（4）构造出近似矩阵，具体算法可参考NocedalandWright(2006).

频率域变密度声波方程可以表达为

其中，κ为模量，ρ为密度，ω为角频率，P(x，ω）为频率域压力波场，S（ω）是炮点xs处的震源函数．目标函数对模量和密度的梯度可以表达为（Tarantola,1984;PrattandWorthington,1990):

其中，Pf为正传波场，Pb为反传残差波场．梯度可以看成是正传波场和反传残差波场零延迟的互相关，而各参数之间辐射模式的不同导致了梯度在具体表达形式上的差异．

由于模量-密度并不是一种好的参数化方式（Tarantola,1986;Mora,1987;Kohnetal.,2012），本文选取速度-密度为模型参数．目标函数对其它参数的梯度可以利用链式法则（Mora,1987）求得：

利用模量、速