[發(fā)明專利]一種全自動的航空電磁數(shù)據(jù)調(diào)平方法及裝置

在審

• 申請?zhí)枺?/td>: 202111371807.2
• 申請日: 2021-11-18
• 公開（公告）號: CN114690253A
• 公開（公告）日: 2022-07-01
• 發(fā)明（設(shè)計(jì)）人: 張瓊;劉云清;顏飛;呂超
• 申請（專利權(quán)）人: 長春理工大學(xué)
• 主分類號: G01V3/38
• 分類號: G01V3/38
• 代理公司: 北京睿博行遠(yuǎn)知識產(chǎn)權(quán)代理有限公司 11297
• 代理人: 董自亮

1.一種全自動的航空電磁數(shù)據(jù)調(diào)平方法，其特征在于：包括以下步驟，

1).對待測區(qū)域進(jìn)行航測并記錄航測數(shù)據(jù)D，航測所得的測區(qū)數(shù)據(jù)中包含多條測線數(shù)據(jù)；

2).錄入航測數(shù)據(jù)，將航測數(shù)據(jù)中異常區(qū)域數(shù)據(jù)剔除，剔除后得到測區(qū)數(shù)據(jù)記為

3).所述的D_{ext_gru}經(jīng)過變分模態(tài)分解后，變分模態(tài)分解后將分解得到2個(gè)模態(tài)分量u₁(t)和u₂(t)；

4).求解調(diào)平誤差D_{ext_error}，模態(tài)分量u₁(t)代入下式

得到調(diào)平誤差D_{ext_error}

5).求解調(diào)平結(jié)果

D_result＝D-D_error，得到各測區(qū)數(shù)據(jù)調(diào)平后調(diào)平航測數(shù)據(jù)結(jié)果；

輸出測區(qū)數(shù)據(jù)調(diào)平結(jié)果D_result

6).依據(jù)步驟5)中輸出的調(diào)平航測數(shù)據(jù)結(jié)果D_result進(jìn)行成圖，得到調(diào)平后航空電磁數(shù)據(jù)圖。

2.根據(jù)權(quán)利要求1所述的一種全自動的航空電磁數(shù)據(jù)調(diào)平方法，其特征在于：步驟2)中，所述的將航測數(shù)據(jù)中異常區(qū)域數(shù)據(jù)剔除，還包括，

在測線數(shù)據(jù)的兩端數(shù)據(jù)分別進(jìn)行鏡像周期延拓，并利用格拉布斯準(zhǔn)則進(jìn)行異常區(qū)域數(shù)據(jù)的剔除處理；

其中數(shù)據(jù)延拓過程中，設(shè)定各測線兩端的延拓測點(diǎn)數(shù)為：

此時(shí)各測線包含測點(diǎn)為：

各測線數(shù)據(jù)經(jīng)鏡像周期延拓后，測區(qū)數(shù)據(jù)表示為：

3.根據(jù)權(quán)利要求2所述的一種全自動的航空電磁數(shù)據(jù)調(diào)平方法，其特征在于：數(shù)據(jù)延拓過程中計(jì)算測點(diǎn)數(shù)據(jù)與該測線數(shù)據(jù)均值的殘差與格拉布斯準(zhǔn)則，當(dāng)測點(diǎn)數(shù)據(jù)與該測線數(shù)據(jù)均值的殘差滿足時(shí)記為異常測點(diǎn)數(shù)據(jù)，并剔除且由該測線的中值數(shù)據(jù)代替：

為該測線數(shù)據(jù)的均值，g(num，a)為臨界值系數(shù)，取決于測量次數(shù)num和顯著性水平a，為測線數(shù)據(jù)的標(biāo)準(zhǔn)差；

測線數(shù)據(jù)剔除所述的異常測點(diǎn)數(shù)據(jù)后，表示為，

4.根據(jù)權(quán)利要求3所述的一種全自動的航空電磁數(shù)據(jù)調(diào)平方法，其特征在于：

并依次對測區(qū)各測線數(shù)據(jù)進(jìn)行異常區(qū)域數(shù)據(jù)的剔除處理，剔除處理后的測區(qū)數(shù)據(jù)記為

5.根據(jù)權(quán)利要求1所述的一種全自動的航空電磁數(shù)據(jù)調(diào)平方法，其特征在于：步驟3)中所述的D_{ext_gru}經(jīng)過變分模態(tài)分解，還包括

①變分模態(tài)分量u(t)、變分模態(tài)分量的中心頻率ω和拉格朗日乘子λ，代入變分模態(tài)分量計(jì)算式循環(huán)更新迭代，在第m次時(shí)候得到：

其中，δ(t)為狄拉克函數(shù)，*表示卷積；

②變換到譜域并積分可得到第m次迭代的變分模態(tài)分量：

③變分模態(tài)分量的中心頻率計(jì)算公式如下：

④變換到譜域，從而得到第m次迭代的變分模態(tài)分量的中心頻率

⑤拉格朗日乘子更新迭代公式如下：

⑥設(shè)置迭代中止條件，符合中止條件時(shí)迭代終止。

對所得到的對變分模態(tài)分量和進(jìn)行傅里葉反變換，得到時(shí)域的分解結(jié)果u₁(t)和u₂(t)，至此完成測區(qū)數(shù)據(jù)D_{ext_gru}的變分模態(tài)分解。

6.根據(jù)權(quán)利要求5所述的一種全自動的航空電磁數(shù)據(jù)調(diào)平方法，其特征在于：所述的變換到譜域具體的為，

將步驟①中得到的式變換為譜域式：

sgn為符號函數(shù)。

7.根據(jù)權(quán)利要求6所述的一種全自動的航空電磁數(shù)據(jù)調(diào)平方法，其特征在于：所述的積分為構(gòu)建埃爾米特矩陣，并將所述的譜域式的兩項(xiàng)結(jié)果寫為非負(fù)頻率上的半空間積分，同時(shí)令ω＝ω+ω₁，得到：