[发明专利]一种基于格林函数的地震监测方法

权利要求书

1.一种基于格林函数的地震监测方法，其特征在于，包括以下步骤：

步骤一：地震监测区域划分为经度、纬度和深度的三维网格，网格所在位置为可能的震源位置；

步骤二：准备格林函数文件；依据一维速度模型、观测台站的位置坐标和步骤一的网格点的经纬度、深度坐标，准备好全部格林函数文件g_i，i∈[0,2]， 计算出格林函数的归一化系数；归一化系数是2g₀g₁、2g₁g₂和2g₀g₂；

步骤三：格林函数与观测波形进行互相关；假定震源在某一个三维网格位置上，将实际的连续观测记录从ENZ坐标旋转到RTZ坐标；将格林函数与实际的连续观测做互相关，并得到半归一化的互相关格林函数：u是某一分量的实际的观测波形；每一个分量都需要各自进行步骤三；

步骤四：对步骤三得到的半归一化的互相关格林函数进行加权叠加，并完成归一化；假定一个震源机制δ,λ，是断层面的走向，δ是断层面的倾角，λ是断层面的滑动角，对格林函数进行加权叠加，即可得平均互相关波形；首先，单分量的互相关波形可以由下面的公式(1)给出：

公式(1)中的已经由步骤二完成；公式(1)中的ai是权重系数，

如果u是Z或R分量，a₁＝sinφcos2δsinλ-cosφcosδcosλ，如果u是T分量，a₀＝0，a₁＝-(sinφ cosδ cosλ+cosφcos2δ sinλ)，

公式(1)中的v是理论波形，|v|由公式(2)得到：

在单分量的互相关波形得到后，叠加后再平均即可得平均互相关波形；将平均互相关波形中振幅达到了既定阈值的采样点的相关信息保存，平均互相关波形不保存；

步骤五：步骤四完成后在完成了一个震源位置上一个震源机制的疑似地震的探测；不断重复步骤四，完成同一个震源位置上不同震源机制的疑似地震的探测；

步骤六：重复步骤三至五，完成所有震源位置上所有震源机制的疑似地震的探测并保留相关信息；

步骤七：依据所有的疑似地震的记录，挑选出确定的探测到的地震的信息。

2.根据权利要求1所述的一种基于格林函数的地震监测方法，其特征在于：步骤中三维网格覆盖需要进行地震监测的区域；三维网格的网格点的分布应该均匀覆盖；如果有需要特别监测的区域需要加密网格点；网格点的间距选择经纬度方向上2km，深度方向上5km。

3.根据权利要求1所述的一种基于格林函数的地震监测方法，其特征在于：步骤七中的挑选方法是：首先选择出互相关值最大的结果作为确定的探测到的地震，然后找到互相关值第二大的结果，如果该结果的发震时刻和前面已经确定的地震的发震时刻不小于固定时窗，则也确定为探测到的地震，否则舍弃；不断重复这个过程，完成对疑似地震的挑选。