[发明专利]一种高精度预测ASF的窄带离散分布抛物方程方法

摘要

本发明公开了一种高精度预测ASF的窄带离散分布抛物方程方法，具体为：首先对罗兰‑C电流时域信号进行采样，对采样后的信号作离散傅里叶变换，将其分解为多个频率电流分量；然后采用平地面公式计算每个频率电流分量在近区所辐射的磁场，并通过近远区边界上的磁场结果作为非均匀网格剖分的离散分布抛物方程方法的初始场计算远区磁场，从而得到每个频率电流分量在地表所产生的磁场；最后，采用基于滑动窗思想的傅里叶逆变换，从得到的多个频率电流分量在地表辐射的频域磁场恢复出时域磁场信号。本发明方法克服现有理论难于预测实际长距离罗兰‑C信号ASF分布的不足，与现有频域方法相比，预测精度明显提高，具有实用性强的特点。

主权项

一种高精度预测ASF的窄带离散分布抛物方程方法，其特征在于，具体按照以下步骤实施：步骤1：对实测罗兰‑C时域电流时域信号i(t)进行采样，采样后得到信号i(n)，其中，n＝0,1,2,…,Nt，将采样后的信号作离散傅里叶变换，分解为多个频率电流分量I(m)，其中，m＝0,1,2,…,Nf；步骤2：当m＝0时，利用平地面公式计算频率f(m)对应的电流分量I(m)在近区即区域I(ρ≤ρ0)产生的磁场分布步骤3：以步骤2计算的近区最大处即近远区边界ρ＝ρ0处的磁场作为远区即区域II(ρ≥ρ0)的激励源，用非均匀网格剖分的离散分布抛物方程方法用于求解远区的磁场分布步骤4：通过步骤2和步骤3可得到频率为f(m)的电流分量I(m)在地面接收点辐射的磁场强度；更新步骤2和步骤3中的m，对其赋值为m＝m+1，当f(m)≤fmax时，重复步骤2和步骤3即可求得罗兰‑C带宽中每个频率处电流分量I(m)在地面接收点辐射的磁场强度步骤5：基于滑动窗的思想，对步骤4中计算的磁场分布，采用离散傅里叶逆变换求解地表任意距离处接收到的时域罗兰‑C信号和滑动时间twindow(ρ)，并提取ASF(ρ)分布。