[发明专利]一种基于克里金插值的多台站远洋磁力勘探方法

权利要求书

1.一种基于克里金插值的多台站远洋磁力勘探方法，其特征在于，包括以下步骤：

S1、计算测区中心点经纬度坐标(x₀,y₀)，获取磁力测量在测区作业的起止时间；

S2、设置搜索半径，获取半径范围内的地磁台站对应时间范围的日变曲线，设获取台站数量为n；

S3、确定距离最近台站并将其设为主站，各台站分别与主站进行互相关，计算相关性指标，删除相关性低的记录，计算个台站与主站记录的相位差，修正日变曲线；

S4、以测区中心为原点，转换台站位置为极坐标，计算内插指标，判断台站分布是否符合内插条件，若是，则进行步骤S5，若否，则扩大搜索半径后返回所述步骤S2；

S5、将各台站日变曲线记录进行转置，按时间分组，将同一时刻各台站日变曲线作为输入，基于克里金发进行插值，获取拟合的测区日变曲线。

2.根据权利要求1所述的一种基于克里金插值的多台站远洋磁力勘探方法，其特征在于，所述步骤S3具体包括：

S31、以距离测区最近台站为主站，设其日变记录为f₁，其他台站日变记录为f_i(i＝2,3,…,n)，所有台站记录采样间隔为Δt，将各台站与主站的日变记录依次做互相关：

其中，k为日变曲线采样点数，m＝-k,-k+1,…,k-1,k；

S32、进行归一化处理，获取相关性指标：

其中：

设置阈值，提出相关性低的台站记录，设剩余台站数量为n′；

S33、以主台站为基准，各台站与其互相关最大值max(R_1l),l＝2,3,…,n′对应的mΔt为相位差，对各台站日变曲线相位调整获得修正后的日变曲线f′_l(t)＝f_l(t+mΔt)，基于最小二乘法进行拟合，修正幅值差，

f″_l(t)＝af′_l(t)+b。

3.根据权利要求1所述的一种基于克里金插值的多台站远洋磁力勘探方法，其特征在于，所述步骤S4具体包括：

S41、设各地磁台站经纬度坐标为(x_s,y_s)，s＝1,2,…,n′，以测区中心点(x₀,y₀)为原点，将各台站坐标转换为极坐标(r_s,θ_s)，其中：

S42、对θ_s由小到大排序，获取有序数组计算相邻两站的夹角：

S43、获取内插指数：

其中，a_max＝max(a_w)，若采用的地磁台站分布不符合内插条件，即γ＝0，则扩大搜索半径r，重复步骤S2-S4，直至满足条件，即γ＝1。

4.根据权利要求1所述的一种基于克里金插值的多台站远洋磁力勘探方法，其特征在于，所述步骤S5具体包括：

将同一时刻各台站日变曲线作为输入，基于克里金法进行插值，获得拟合的测区日变曲线：

f₀(t)＝K({f_s(t)}),s＝1,2,…,n

其中，K为克里金算子，f_s(t)为输入的台站日变曲线。