1.一种基于初至时差和波形叠加的微地震震源定位方法，其特征在于，包括以下步骤：S1：输入速度模型；

S2：拾取并输入实际初至时间，以及读取地震数据；

S3：依据步骤S1中速度模型，计算可行解区域内所有网格点到每一个检波器的理论初至时间表；

S4：基于步骤S2中输入的实际初至时间和步骤S3中计算的理论初至时间表，依据公式(1)和(2)构造走时残差函数Tr；

上式中：M、N分别是观测到的横波和纵波初至的个数； 和 分别是第j个检波器记录的横波和纵波初至时间， 和 分别是与之对应的第j个检波器的横波及纵波理论计算时间； 和 分别是第i个检波器记录的横波和纵波初至时间， 和 分别是与之对应的第i个检波器的横波及纵波理论计算时间；Tshift是观测及理论初至之间的一个恒定的漂移量，以解决实际资料监测中微地震震源发震时间未知的难题；

S5：将步骤S2中读入的地震数据、走时目标函数对应微地震事件中第m个检波器记录的可识别纵波初至时间 和步骤S3中理论初至时间表结合，通过极性修正，然后依据公式(3)、(4)和(5)构造波形叠加函数Ews；

式中ΔP、ΔS分别是纵横波滑动时窗样点数；dt是采样间隔；n是检波器个数；

分别是假设震源点到第i个检波器的纵横波理论时间 与第m个检波器的纵波理论时间 之差， 分别是第i个检波器纵横波振幅的符号系数，ui是地震数据第i道振幅值；

S6：输入权重系数β，依据步骤S4中走时残差函数Tr和步骤S5中波形叠加函数Ews，按照公式(6)构造改进目标函数；

Q＝Tr-βEws                  (6)

S7：通过网格搜索法寻找改进目标函数的最小值，输出此时对应的最优解，即震源位置。

2.根据权利要求1所述的一种基于初至时差和波形叠加的微地震震源定位方法，其特征在于：在所述步骤S1中，速度模型可通过声波测井资料和地质分层资料建立初始速度模型，然后利用射孔资料进行速度模型校正。

3.根据权利要求1或2所述的一种基于初至时差和波形叠加的微地震震源定位方法，其特征在于：在所述步骤S3中，可行解区域内所有网格点是指以射孔点为中心，给定一个求解区域，然后对该区域进行网格化，每个网格点均为可能的震源点。

4.根据权利要求3所述的一种基于初至时差和波形叠加的微地震震源定位方法，其特征在于：在所述步骤S5中，极性修正可通过符号系数 实现，即判断检波器记录中时窗内振幅符号，需要正负号保持一致，此处保证波形的同向叠加。

5.根据权利要求4所述的一种基于初至时差和波形叠加的微地震震源定位方法，其特征在于：在所述步骤S6中，权重系数β>0，大小与微地震资料的信噪比情况成反比。

6.根据权利要求5所述的一种基于初至时差和波形叠加的微地震震源定位方法，其特征在于：在所述步骤S7中，在网格搜索定位中，假设网格点是震源位置时，Tr为最小，Ews为最大，即-βEws为最小，故改进目标函数Q为最小值。