周嘉欣， 杨锴， 邵炜栋

1 中煤科工集团常州研究院有限公司， 江苏常州 213015 2 天地(常州)自动化股份有限公司， 江苏常州 213015 3 同济大学海洋地质国家重点实验室， 上海 200092 4 阿里巴巴集团， 杭州 311121

0 引言

各向异性在地下介质中广泛存在，在地震波传播中各向异性会同时影响信号的相位和振幅.对于大偏移距数据，不考虑各向异性会导致速度估计发生很大误差.这促使近年来学术界和工业界在研究基于地震波的反演、成像以及储层建模技术时，必须将算法从各向同性拓展到各向异性(Tsvankin et al.，2010).因此面向各向异性参数估计的层析成像技术成为现代地震反演与成像技术中的一个重要组成部分.然而，由于各向异性层析成像方法引入了更多未知参量，使得参数估计的不确定性和复杂度都大大增加，因此长期以来一直是一个非常有挑战的研究课题.

基于走时的各向异性层析成像已经有大量研究工作(Chapman and Pratt，1992；Zhou et al., 2008; Wang and Tsvankin, 2013；卢明辉等, 2005； 刘玉柱等，2014；黄光南等，2015).20世纪末以来，一种全新的层析成像算法——立体层析，开始为学术界和工业界所知.立体层析将局部相干同相轴在炮集与检波点道集内的射线参数(或P参数)水平分量，以及炮点坐标和检波点坐标都纳入到反演的数据空间之中，重新定义了层析反演的数据分量和模型分量，使得立体层析成为射线类层析反演方法中唯一一种可以同时反演速度结构、反射点位置与反射层局部形态的方法(Billette and Lambare, 1998).国内外学者在各向同性立体层……