[發明專利]一種基于多重分形克里金法的局部地磁圖構建方法

有效

• 申請號: 201210237318.2
• 申請日: 2012-07-10
• 公開（公告）號: CN102841385A
• 公開（公告）日: 2012-12-26
• 發明（設計）人: 趙玉新;常帥;劉廠;張振興;沈志峰
• 申請（專利權）人: 哈爾濱工程大學
• 主分類號: G01V3/40
• 分類號: G01V3/40
• 代理公司: 北京永創新實專利事務所 11121
• 代理人: 趙文利

1.一種基于多重分形克里金法的局部地磁圖構建方法，其特征在于，包括以下幾個步驟：

步驟一：去除實測數據中的主磁場成分；

判斷實測地磁數據中是否包含地球主磁場成分，如果包含，利用地磁場模型計算各測點主磁場強度，并在實測數據中將其去除；

步驟二：選定估值區域；

根據實測地磁數據的分辨率以及地磁異常場的空間特性，以待插值點為中心，選定估值區域；

步驟三：選定高斯模型描述地磁異常場的變異特性，擬合高斯模型中的未知參數，并利用此模型計算估值區域內各實測點之間以及各實測點和待插值點之間的半方差；

在估值區域內，利用實測地磁數據，計算不同的分離距離對應的變異函數值，從而得到分離距離與變異函數值序列對，根據此序列對，利用擬合法推算出高斯模型中的未知參數，并利用此模型計算估值區域內各實測點之間以及實測點和待插值點之間的半方差；

步驟四：計算估值區域內各測點的權重系數；

利用步驟三中得出的半方差值，通過解空間相關系數矩陣方程求出各實測點的權重系數；

步驟五：利用多重分形理論推算待插值點鄰域內的測度表達式；

在估值區域內，利用測度與尺度的在雙對數坐標系中的線性關系，用最小二乘擬合法得出插值點鄰域內的分形維數，進而推算出以待插值點為中心的小正方形的測度與整個估值區域測度的關系式；

步驟六：根據步驟四中利用克里金法求出的各實測點的權重系數，推算出待插值點鄰域內的多重分形克里金插值方程；

利用步驟四中求出的估值區域內各實測點的權重系數，求出整個估值區域測度表達式，并根據步驟五中得到的關系式，推算出待插值點在其小鄰域內的多重分形克里金法插值方程；

步驟七：采用交叉驗證法驗證插值的精度并構建地磁圖；

設X_i為測區內一實測點，將此實測點暫時去除，利用步驟二至步驟六所述方法對此點進行插值，得到這一點的地磁異常估計值?，按步驟二至步驟六得到測區內所有實測點的地磁異常估計值，并采用均方根預測誤差、平均估計誤差百分比、相對均方差三個指標對插值結果進行評估；

設定均方根預測誤差、平均估計誤差百分比、相對均方差三個指標的數值，如果插值結果符合指標要求，則根據插值后的地磁數據構建地磁基準圖；否則，返回步驟二，重新進行插值。