[发明专利]一种大地电磁时间域标定方法及系统

权利要求书

1.一种大地电磁时间域标定方法，其特征在于，包括：

获取大地电磁测深仪中时间序列的系统响应、各频率的幅值校正量以及各频率的相位校正量；

根据所述幅值校正量和所述相位校正量确定标定的传输函数；

根据所述标定的传输函数采用傅里叶反变换方法确定标定的单位脉冲响应函数；

采用希尔伯特变换法确定所述系统响应的复数序列；

将所述单位脉冲响应函数与所述复数序列进行卷积操作，得到标定后的时间序列；

其中，

所述根据所述标定的传输函数采用傅里叶反变换方法确定标定的单位脉冲响应函数，具体包括：

根据所述标定的传输函数确定双边谱传输函数；所述双边谱传输函数的公式如下：

式中，H'(ω)为双边谱传输函数，H(ω)为标定的传输函数，ω₀为H(ω)频率范围上确界，ω为离散的角频率；

对所述双边谱传输函数采用三次样条插值法，得到插值后的传输函数；

对所述插值后的传输函数进行离散傅里叶反变换，得到单位脉冲响应函数；所述单位脉冲响应函数的公式如下：

式中，h(n)为单位脉冲响应函数，H(ω)为插值后的传输函数，N为H(ω)序列的长度，n为时间点自变量；

所述采用希尔伯特变换法确定所述系统响应的复数序列，具体包括：

对所述系统响应进行希尔伯特变换，得到系统响应对应的虚部序列；所述虚部序列的公式如下：

式中，y'(n)为系统响应对应的虚部序列，M为系统响应的序列长度，y(m)为第m个时间点的电磁场大小，n与m为不同的时间点，m为时间点卷积变量，n＝1,2,3,...,M，m＝1,2,3,...,M；

根据所述系统响应和所述系统响应对应的虚部序列确定所述系统响应的复数序列；所述复数序列的公式如下：

y(n)＝y(n)+iy'(n)

式中，y(n)为自变量为n的复数序列，y(n)为系统响应；

所述标定的传输函数公式如下：

H(ω)＝A(ω)e^iφ(ω) 0≤ω≤ω₀

式中，H(ω)为标定的传输函数，A(ω)为幅值校正量，φ(ω)为相位校正量；

所述标定后的时间序列公式如下：

式中，y₁(n)为标定后的时间序列，y(m)为卷积变量为m的复数序列，h(n-m)为单位脉冲响应函数。

2.根据权利要求1所述的大地电磁时间域标定方法，其特征在于，在所述将所述单位脉冲响应函数与所述复数序列进行卷积操作，得到标定后的时间序列，之后还包括：

将所述标定后的时间序列中前M项作为大地电磁时间域标定结果；所述M的个数与所述系统响应的序列长度相等。