關鍵詞:大地電磁 三次插值 自適應有限元法
摘要:大地電磁(MT)數值模擬中通常使用有限單元法,通過伽遼金(Galerkin)法將微分方程轉化為與其等價的泛函形式,對泛函求取極值并在單元上定義插值基函數,得到節點上電磁場值的線性方程組,最終形成大型復對稱稀疏矩陣。要達到較高的有限元計算精度,一般采用密集的網格或高次插值的方法,這樣做大大的減慢了正演的速度。結合兩者的優點利用三次插值和h-型自適應相結合的有限元法來實現MT的正演算法。首先從一個粗網格出發并利用三次插值,通過后驗誤差估計方法局部加密網格,在計算量較小的情況獲得較高的計算精度。這種方法可以針對目標區域和介質分界面發生突變處進行網格加密,不需要全局加密網格。最后通過對國際標準模型COMMEMI-2D1的模擬,分別比較二次插值與三次插值的自適應網格數量和數值模擬結果,證明了三次插值自適應有限元算法的可行性。
物探化探計算技術雜志要求:
{1}若需列出項目資助相關內容,均置于首頁底部,并于論文題目后標出星號。
{2}來稿不退,請自留底稿。如三個月未見通知,作者可另行處置。
{3}論文中插人圖表時,說明性文字要先出現,圖表附在說明性文字之后。圖的題目在圖之下,表的題目在表格之上。
{4}注釋:一律采用腳注,每頁編號,自為起止。
{5}文內體例要求,一般以“一、”、“(一)”、“1.”、“(1)”作為文章層次,通過簡短的小標題方式提煉主要觀點。
注:因版權方要求,不能公開全文,如需全文,請咨詢雜志社