[發明專利]一種基于數據空間快速計算的聯合反演方法

有效

• 申請號: 201810323943.6
• 申請日: 2018-04-12
• 公開（公告）號: CN108229082B
• 公開（公告）日: 2019-12-06
• 發明（設計）人: 張镕哲;李桐林;鄧馨卉;莊嚴;侯宇健;石會彥;楊賀天;李少朋;白云天
• 申請（專利權）人: 吉林大學
• 主分類號: G06F17/50
• 分類號: G06F17/50
• 代理公司: 21229 沈陽維特專利商標事務所(普通合伙)
• 代理人: 楊群<國際申請>=<國際公布>=<進入國

本發明公開了一種基于數據空間快速計算的聯合反演方法，包括如下步驟：S1：建立重磁電震聯合反演目標函數；S2：建立電阻率、速度、密度和磁化率的初始模型；S3：正演計算和雅克比矩陣的求取；S4：計算交叉梯度偏導數和拉格朗日算子；S5：計算獲得模型空間迭代模型；S6：計算獲得數據空間迭代模型；S7：對獲得的數據空間迭代模型進行正演計算，并與觀測數據進行擬合差求取，如果擬合差增大、達到最大迭代次數或者擬合差達到閾值，則停止循環反演，輸出該數據空間聯合反演模型，否則返回S3。該基于數據空間快速計算的聯合反演方法計算速度快、占用內存小。

技術領域

本發明涉及多物性聯合反演的處理技術領域，特別提供了一種基于數據空間快速計算的聯合反演方法。

背景技術

綜合地球物理方法已經廣泛應用于深部地質探測和礦產、水、石油、地熱等資源勘查。而聯合反演是其綜合地球物理方法數據處理的主流方法，國外科學家已經提出了二種或多種物性參數之間結構耦合的聯合反演方法，國內也實現了二種或多種物性參數之間的結構耦合聯合反演方法。

目前，對開展大數據量的聯合反演數據處理已經成為重要的發展方向。現有的聯合反演多采用在模型空間下進行反演計算，但是當反演物性參數種類從二個增加到多個或者反演區域網格剖分數據量過大時，對計算機的計算速度和內存占用都將成為考驗，將會嚴重制約了綜合地球物理方法的發展。

發明內容

鑒于此，本發明的目的在于提供一種基于數據空間快速計算的聯合反演方法，以解決現有的聯合反演方法計算速度慢、占用內存大等問題。

本發明提供的技術方案是：一種基于數據空間快速計算的聯合反演方法，包括如下步驟：

S1：建立重磁電震聯合反演目標函數，所述聯合反演目標函數如式(1)：

約束條件：τ(m)＝0

其中：

C_d＝diag[C_d1,C_d2,C_d3,C_d4]，C_m＝diag[C_m1,C_m2,C_m3,C_m4]，

α＝[α₁,α₂,α₃,α₄]，