孙 强，郑建常，王 鹏

（1.淄博市应急管理局，山东 淄博 255086；2.山东省地震局，济南 250014；3.东营市地震监测中心，山东 东营 257000）

0 引言

目前使用地震波资料求取地壳构造应力场的方法主要有两种：一是使用地震的震源机制解来反演应力场；二是使用P 波初动符号计算综合震源机制解的方法获得应力场方向[1]。使用P 波初动求取单个地震震源机制解需要较多的P 波初动信息，并且记录地震波的台站需要合理的覆盖，否则计算出的震源机制解会有较大的误差。再者有许多小微震的震源机制解由于数据原因并不能准确求取，这样通过震源机制解求解应力场的数据量就会大大减少，降低了解的可靠性。而通过P 波初动计算综合震源机制解方法则可以利用众多小微震的P 波初动信息，增大了输入的信息量，提高了结果的置信度。

利用大量小震P 波初动符号推断应力场的方法最早是由Aki[2]提出的，该方法充分运用大量不能确定震源机制解的小微震的P 波初动信息，因此获得较广泛的应用。李钦祖等[3]运用该法对单个台站的小震资料分别求取了红山台和沙城台的区域应力场方向，并给出了其置信范围；许忠淮等[4]使用该方法运用多个台站和多个地震资料求取了区域应力场的方向，并通过合成数据的方法进行验证，确定了该方法能得到稳定的应力场；薛宏运等[5]运用该方法求解了鄂尔多斯周围地区分区现代地壳应力场并进行了分析；汪素云等[6]运用该法求取了中国东部大陆分区的地震构造应力场并进行了分析；……