[發(fā)明專利]一種應(yīng)用于可控源電磁法相關(guān)辨識(shí)中的降噪方法、系統(tǒng)及可讀存儲(chǔ)介質(zhì)

在審

• 申請(qǐng)?zhí)枺?/td>: 202110038837.5
• 申請(qǐng)日: 2021-01-12
• 公開（公告）號(hào): CN112765546A
• 公開（公告）日: 2021-05-07
• 發(fā)明（設(shè)計(jì)）人: 湯井田;楊智;肖曉
• 申請(qǐng)（專利權(quán)）人: 中南大學(xué)
• 主分類號(hào): G06F17/15
• 分類號(hào): G06F17/15;G06F17/16
• 代理公司: 長(zhǎng)沙市融智專利事務(wù)所(普通合伙) 43114
• 代理人: 姚瑤

1.一種應(yīng)用于可控源電磁法相關(guān)辨識(shí)中的降噪方法，其特征在于：包括如下步驟：

設(shè)定發(fā)射信號(hào)的時(shí)長(zhǎng)T與工頻噪聲周期T_c的關(guān)系；

連續(xù)發(fā)出發(fā)射信號(hào)傳送至大地系統(tǒng)，并基于T與T_c的關(guān)系確定等效發(fā)射信號(hào)，其中，基于等效發(fā)射信號(hào)對(duì)應(yīng)的等效接收信號(hào)中不包含工頻噪聲；

若發(fā)射信號(hào)的時(shí)長(zhǎng)T為工頻噪聲周期T_c的整數(shù)倍，等效發(fā)射信號(hào)為連續(xù)發(fā)射的相鄰兩組發(fā)射信號(hào)相減得到；若發(fā)射信號(hào)的時(shí)長(zhǎng)T為工頻噪聲周期T_c的整數(shù)倍加二分之一個(gè)T_c，等效發(fā)射信號(hào)為連續(xù)發(fā)射的相鄰兩組發(fā)射信號(hào)相加得到。

2.根據(jù)權(quán)利要求1所述的降噪方法，其特征在于：所述發(fā)射信號(hào)的時(shí)長(zhǎng)T等于發(fā)射信號(hào)的有效激勵(lì)時(shí)間與零激勵(lì)時(shí)間之和。

3.根據(jù)權(quán)利要求2所述的降噪方法，其特征在于：若發(fā)射信號(hào)為m序列信號(hào)，發(fā)射信號(hào)的有效激勵(lì)時(shí)間以及零激勵(lì)時(shí)間對(duì)應(yīng)的補(bǔ)零的個(gè)數(shù)如下所示：

T_v＝(N_p*L_m*N_spb)/f_s

式中，N_pz為零激勵(lì)時(shí)間對(duì)應(yīng)的補(bǔ)零的個(gè)數(shù)，T_pz為零激勵(lì)時(shí)間，T_v為發(fā)射信號(hào)的有效激勵(lì)時(shí)間，N_p為m序列的周期數(shù)，L_m為一個(gè)m序列周期的碼元數(shù)，N_spb為每個(gè)碼元的采樣數(shù)，f_s為采樣頻率，K為正整數(shù)；

其中，case1表示發(fā)射信號(hào)的時(shí)長(zhǎng)T為工頻噪聲周期T_c的K倍，case2表示發(fā)射信號(hào)的時(shí)長(zhǎng)T為工頻噪聲周期T_c的K倍加二分之一個(gè)T_c。

4.根據(jù)權(quán)利要求1所述的降噪方法，其特征在于：還包括利用等效發(fā)射信號(hào)和等效接收信號(hào)辨識(shí)大地脈沖響應(yīng)，其中，等效發(fā)射信號(hào)需滿足：

case1：當(dāng)發(fā)射信號(hào)的時(shí)長(zhǎng)T為工頻噪聲周期T_c的K倍時(shí)，連續(xù)發(fā)射的相鄰兩組發(fā)射信號(hào)x₁(t)、x₂(t)分別表示為：x₁(t)＝x_m(t)、x₂(t)＝-x_m(t)，對(duì)應(yīng)得到的等效發(fā)射信號(hào)為m序列2x_m(t)；

case2：當(dāng)發(fā)射信號(hào)的時(shí)長(zhǎng)T為工頻噪聲周期T_c的K倍加二分之一個(gè)T_c時(shí)，連續(xù)發(fā)射的相鄰兩組發(fā)射信號(hào)x₁(t)、x₂(t)分別表示為：x₁(t)＝x_m(t)、x₂(t)＝x_m(t)，對(duì)應(yīng)得到的等效發(fā)射信號(hào)為m序列2x_m(t)，其中，x_m(t)為m序列。