[發明專利]一種火山巖型鈾礦床磁性界面反演方法

有效

• 申請號: 201510729695.1
• 申請日: 2015-10-30
• 公開（公告）號: CN106650192B
• 公開（公告）日: 2019-02-26
• 發明（設計）人: 喻翔;陳聰
• 申請（專利權）人: 核工業北京地質研究院
• 主分類號: G06F17/10
• 分類號: G06F17/10
• 代理公司: 核工業專利中心 11007
• 代理人: 高尚梅

1.一種火山巖型鈾礦床磁性界面反演方法，其特征在于：具體包括以下步驟：

步驟一、在勘查區開展高精度磁法測量，獲得磁測數據；

步驟二、在勘查區開展RTK測量，獲得勘查區的測點精確經緯度坐標和高程數據；

步驟三、收集勘查區地層、巖體磁化率數據并把勘查區地層、巖體磁化率數據轉換成CGS國際單位制；

步驟四、對步驟一獲得的磁測數據進行日變校正、正常場梯度改正、高度改正及基點改正，得到測點異常△T值；

步驟五、對步驟二得到的測點高程數據進行網格化處理，對步驟四所得的磁異常△T值進行數據平滑處理，得到光滑的磁異常曲線并對所有光滑處理后的磁異常△T值進行網格化處理；

步驟六、對步驟五所得到的磁異常△T網格數據進行化極處理，化極處理所需的磁傾角、磁偏角由步驟三所得出；

步驟七、開始創建反演初始模型，并把模型分為兩層，根據步驟三中不同巖體磁化率數據對兩層模型進行磁化率賦值，針對反演初始模型的兩層模型的交界面創建反演目標界面常量網格；

步驟八、對目標界面進行反演，在反演過程中加入余弦濾波器，余弦濾波器選擇低通濾波，且設置上限截止波長，下限截止波長，反演出磁性界面網格；

步驟九、修改余弦濾波器參數上限截止波長與下限截止波長；

步驟十、重復步驟八、步驟九，迭代次數為M次，直至擬合誤差達到5nT，獲得最終磁性界面網格。

2.根據權利要求1所述的一種火山巖型鈾礦床磁性界面反演方法，其特征在于：所述的步驟一中，磁測總精度應高于2nT，勘查區選在地形平坦的區域。

3.根據權利要求2所述的一種火山巖型鈾礦床磁性界面反演方法，其特征在于：所述的步驟五、步驟六中，所述的網格化處理中的差值方法為最小曲率法。

4.根據權利要求3所述的一種火山巖型鈾礦床磁性界面反演方法，其特征在于：所述的步驟七中，反演目標界面常量網格中，網格常量為-1000m、所需的反演觀測面選擇步驟五得出的高程網格。

5.根據權利要求4所述的一種火山巖型鈾礦床磁性界面反演方法，其特征在于：所述的步驟八中，對目標界面進行反演采用帕克反演方法采用下式進行計算：

式中為磁性界面起伏的傅里葉變換，ΔT為磁異常，J是上下界面磁化率差值，S為徑向波數，n為迭代次數。

6.根據權利要求5所述的一種火山巖型鈾礦床磁性界面反演方法，其特征在于：所述的步驟九中，對步驟八得到的磁性界面網格利用帕克正演公式進行正演計算，所述的帕克正演方法采用下式進行計算：

式中為磁性界面起伏的傅里葉變換，為反演出磁性界面起伏深度，J是上下界面磁化率差值，S為徑向波數，n為迭代次數。

7.根據權利要求6所述的一種火山巖型鈾礦床磁性界面反演方法，其特征在于：所述的步驟十中，M≥10。