[發(fā)明專利]全張量重力梯度儀數(shù)據(jù)質(zhì)量評(píng)價(jià)和濾波處理方法

在審

• 申請(qǐng)?zhí)枺?/td>: 202010257866.6
• 申請(qǐng)日: 2020-04-03
• 公開(kāi)（公告）號(hào): CN111427096A
• 公開(kāi)（公告）日: 2020-07-17
• 發(fā)明（設(shè)計(jì)）人: 袁園
• 申請(qǐng)（專利權(quán)）人: 自然資源部第二海洋研究所
• 主分類(lèi)號(hào): G01V7/00
• 分類(lèi)號(hào): G01V7/00;G01V13/00
• 代理公司: 北京方韜法業(yè)專利代理事務(wù)所(普通合伙) 11303
• 代理人: 黨小林

1.一種全張量重力梯度儀數(shù)據(jù)質(zhì)量評(píng)價(jià)和濾波處理方法，其特征在于，包括：

通過(guò)奇偶網(wǎng)格分析法分別對(duì)全張量重力梯度儀測(cè)量的各個(gè)梯度分量進(jìn)行噪聲估計(jì)；

根據(jù)由噪聲估計(jì)結(jié)果確定的各重力梯度分量權(quán)值，以及拉普拉斯方程的傅里葉級(jí)數(shù)解，構(gòu)建加權(quán)聯(lián)合反演濾波方程；

利用最優(yōu)線性反演方法求解加權(quán)聯(lián)合反演濾波方程中的傅里葉級(jí)數(shù)系數(shù)；

利用計(jì)算得到的傅里葉級(jí)數(shù)系數(shù)，正演計(jì)算得到濾波結(jié)果。

2.根據(jù)權(quán)利要求1所述的全張量重力梯度儀數(shù)據(jù)質(zhì)量評(píng)價(jià)和濾波處理方法，其特征在于，通過(guò)奇偶網(wǎng)格分析法分別對(duì)全張量重力梯度儀測(cè)量的各個(gè)梯度分量進(jìn)行噪聲估計(jì)，分析其噪聲水平，包括：

奇偶網(wǎng)格分析法利用奇數(shù)測(cè)線和偶數(shù)測(cè)線分別對(duì)數(shù)據(jù)進(jìn)行網(wǎng)格，利用兩個(gè)獨(dú)立網(wǎng)格數(shù)據(jù)的和與差來(lái)實(shí)現(xiàn)對(duì)噪聲的估計(jì)。

3.根據(jù)權(quán)利要求2所述的全張量重力梯度儀數(shù)據(jù)質(zhì)量評(píng)價(jià)和濾波處理方法，其特征在于，利用兩個(gè)獨(dú)立網(wǎng)格數(shù)據(jù)的和與差來(lái)實(shí)現(xiàn)對(duì)噪聲的估計(jì)，包括：

根據(jù)如下公式對(duì)噪聲的均方根誤差進(jìn)行估計(jì)：

n_d(i,j)＝[s_o(i,j)+n_o(i,j)]-[s_e(i,j)+n_e(i,j)]＝n_o(i,j)-n_e(i,j)

其中，n_d(i,j)為奇偶網(wǎng)格差，s_o(i,j)為奇數(shù)網(wǎng)格信號(hào)成分，s_e(i,j)為偶數(shù)網(wǎng)格信號(hào)成分，n_o(i,j)為奇數(shù)網(wǎng)格噪聲成分，n_e(i,j)為偶數(shù)網(wǎng)格噪聲成分。

4.根據(jù)權(quán)利要求1所述的全張量重力梯度儀數(shù)據(jù)質(zhì)量評(píng)價(jià)和濾波處理方法，其特征在于，加權(quán)聯(lián)合反演濾波方程具有以下形式：

Ac＝d

其中，A為方程系數(shù)，c為要擬合的系數(shù),d為測(cè)量的重力梯度分量值。

5.根據(jù)權(quán)利要求4所述的全張量重力梯度儀數(shù)據(jù)質(zhì)量評(píng)價(jià)和濾波處理方法，其特征在于，利用最優(yōu)線性反演方法求解加權(quán)聯(lián)合反演濾波方程中的傅里葉級(jí)數(shù)系數(shù)，還包括：

對(duì)目標(biāo)函數(shù)添加正則化項(xiàng)。

6.根據(jù)權(quán)利要求5所述的全張量重力梯度儀數(shù)據(jù)質(zhì)量評(píng)價(jià)和濾波處理方法，其特征在于，正則化項(xiàng)具有如下形式：

其中，μ為正則化參數(shù)。

7.根據(jù)權(quán)利要求4所述的全張量重力梯度儀數(shù)據(jù)質(zhì)量評(píng)價(jià)和濾波處理方法，其特征在于，利用最優(yōu)線性反演方法求解加權(quán)聯(lián)合反演濾波方程中的傅里葉級(jí)數(shù)系統(tǒng)，還包括：

根據(jù)噪聲估計(jì)結(jié)果，定義不同重力梯度分量的方向加權(quán)矩陣。

8.根據(jù)權(quán)利要求7所述的全張量重力梯度儀數(shù)據(jù)質(zhì)量評(píng)價(jià)和濾波處理方法，其特征在于，方向加權(quán)矩陣具有如下形式：

W＝diag(1/σ₁,…,1/σ₆)

其中，σ_i為第i個(gè)數(shù)據(jù)的標(biāo)準(zhǔn)偏差。