王议迎，丁仁伟,2，3，李建平,4，赵俐红,2，3，赵 硕，张硕伟

(1. 山东科技大学 地球科学与工程学院，山东 青岛 266590; 2. 青岛海洋科学与技术国家实验室，山东 青岛 266237; 3. 海洋矿产资源评价与探测技术国家功能实验室，山东 青岛 266237;4. 山东省物化探勘查院，山东 济南 221116)

微地震监测已经成为当今油藏地球物理描述和二氧化碳地质储存检测等方面的关键技术手段[1]。在非常规油气开采过程中，为了提高油气产量,在储层附近进行高压注液，从而产生微裂缝、诱发微地震。实时绘制微震震源，能更好地控制注入过程从而提高储层开发的效率[2]。因此，微震的自动检测对于监测诱导水力压裂过程至关重要。

随着微震记录数据规模的增大，手动拾取微震到时变得更具挑战性和不可靠，且可能会产生若干假警报或出现到时误拾现象。因此，各种自动或半自动的方法被提出并广泛用于微地震震相拾取。现有的微震震相自动拾取方法主要有长短时窗平均能量比法(short-term average/long-term average, STA/LTA)[3]，该算法简单、计算效率高、适于实时处理。在过去的几十年里，许多研究人员对该方法进行了改进，通过构建不同的特征函数或构建能反映信号振幅强度与频率变化的权值函数等方式获取更准确的微震到时拾取结果[4-10]。自回归(auto regressive, AR)是另一种常用的到时自动拾取方法[11]，该方法基于统计特性的差异将微震记录分为真实微震信号到达前和到达后两个部分，用以表征微震记录。在此基础上，Akaike[12]提出了Akaike信息准则(Akaike information criterion, AIC)方法，将微震数据建模为AR过程，将AIC最小值设定为微震事件的起始时间。……