李 铂，崔 鑫，叶庆东，于 澄

（1.中国科学技术大学，合肥 230026；2.山东省地震局，济南 250014；3.中国地震局地球物理研究所，北京 100080）

0 引言

地震波速度结构的反演是地震学的经典问题，精准的速度模型是许多研究的基础性资料，如地震定位、大尺度构造应力场研究、地球动力学等。速度结构研究主要可分为一维速度结构和三维速度结构研究，一维速度结构经典的研究方法主要有Herglotz－Wiechert方 法、Gutenberg 反 演 方 法、τ（p）法［1］等，主要依据大量地震资料的走时曲线计算速度结构，目前应用较少。地震波速度结构研究有2个主要的难点，一是地震波速度结构和地震定位的耦合问题，两者会互相影响其误差和准确性；二是地震数据的结构问题，即地震台站深度普遍位于近地表，难于深入地下，从而难以在深度上进行有效约束。为了解决地震波速度结构和地震定位的耦合性问题，在Geiger定位方法的基础上，1976年，Crosson［2］首次提出了震源参数和速度结构联合反演方法，将速度结构作为未知参数与震源参数同时反演，由此弥补选取速度模型引起的误差，得到了广泛的应用。Aki［3－4］等人将该方法推广到三维速度结构，Pavlis［5］，Spencer［6］，刘福田［7－8］等人相继提出了参数分离法，大大提高了运算效率。Kissling［9］等于1994年提出可以基于走时残差均方根最小作为目标函数进行一维速度结构和震源位置联合反演，定位精度更高。

基于现在全球数字地震观测技术的发展和观测台网的大量建设，积累了大量高质量的地震波走时数据和数字波形记录，使得三维速度结构的反演方法取得了长足的进步。……