郝识杰，黄周传，王良书，徐鸣洁，米 宁，于大勇

南京大学 地球科学与工程学院，南京 210023

1 引言

青藏高原是印度板块和欧亚板块碰撞形成的年轻的高原（图1；Yin and Harrison, 2000）。对于青藏高原的隆升机制存在热烈的讨论，前人提出许多端元模型，如地壳流模型（Clark and Royden,2000; Royden et al., 2008），刚性块体挤出模型（Tapponnier et al., 2001），岩石圈连续变形模型（Molnar and Tapponnier, 1975）等。获取精细的岩石圈地震波速度结构，对于讨论青藏高原的构造演化过程，具有重要的意义，广泛布设的密集地震观测台阵为获取更高精度的地球内部结构提供了可能。

“中国地震科学台阵南北地震带北段”观测项目（ChinArray 二期），在青藏高原东北缘及其相邻区域布设了675台流动地震仪。面波速度成像是利用密集台阵数据，对岩石圈结构进行约束的有效方法。传统的面波成像研究多基于射线假设，即假设地震台接收到的面波信号是沿震源与台站之间的大圆弧路径传播（Barmin et al., 2001）。在实际情况下，由于岩石圈波速结构存在横向变化，面波的传播方向会偏离大圆弧路径，导致基于射线假设的方法出现误差。Foster等（2014）利用类似聚束分析的方法分析了USArray台阵记录的地震事件的入射角度，发现入射角度的偏差可以达到15°，由此造成的相速度误差可以达到4%。Lin等（2009）提出了基于程函方程的成像方法，该方法首先计算面波的相位走时场，然后基于程函方程计算相速度的大小和方向。在该方法中，面波速度的方向从地震记录计算得来，理论上可以减少因面波传播偏离大圆弧路径导致的误差。……