于悦颖 李正楷 周昱辰 杨 驰

（中国南京 210014 江苏省地震局）

0 引言

作为一种全局优化算法，模拟退火算法因算法简单、可对目标函数全局寻优且不易陷入局部极小值，最初被应用于计算机辅助设计、图像处理、代码设计等领域，帮助解决组合优化问题（Van Laarhoven et al，1987），继而被广泛用于地球物理反演。1994 年Billings将模拟退火算法引入地震定位，发现可以得到更可靠和准确的全局最小值解。此后，利用模拟退火算法，Gao 等（2002）对1999 年9 月21 日中国台湾集集地震进行了重定位，窦立婷等（2018）、杨波等（2019）对模拟地震数据和少量大陆震例进行了地震定位，尹奇峰等（2019）对井中微地震事件进行了反演定位，均表明模拟退火算法能够对实际地震事件进行地震定位，具有地球物理实用价值。与Geiger（1912）提出的经典线性方法相比，模拟退火算法在解决地震定位问题时避免了线性方法处理过程中容易出现的高阶项省略、初始模型选定等问题，不会出现线性迭代中陷入局部极小值的情况（师学明等，2007；田玥等，2002）。此外，模拟退火算法不需要进行偏导数计算，可以使用更少的计算量定位到全局最小值附近，有效提高了计算效率和精度（Billings，1994）。

为验证模拟退火算法在江苏及邻区地震事件定位的适用性，文中利用P 波走时数据构建全局寻优模型，以2019 年4 月6 日江苏溧水ML3.3 地震为例，进行反演试算，确定合适的初始温度、扰动函数和降温策略，运用于2019 年江苏省测震台网中心记录的105 次ML≥1.8 地震事件重定位，并对定位误差进行分析。……