為了計算一維地電模型的大地電磁響應,采用高精度的Chebyshev譜方法進行了數(shù)值模擬。首先,從磁場滿足的微分方程出發(fā),利用Chebyshev求導矩陣導出了一維大地電磁正演計算的線性方程組;其次,利用Matlab工具編寫了一維大地電磁響應的Chebyshev譜方法計算程序,同時通過計算均勻半空間模型的數(shù)值解,并與理論解析結(jié)果對比,驗證了正演算法的正確性;最后,通過對二層D型地電模型的大地電磁響應模擬,說明了Chebyshev譜算法的有效性,并總結(jié)了近地表單元剖分間距對大地電磁響應的影響。 

【文章來源】：工程地球物理學報. 2020,17(06)

【文章頁數(shù)】：7 頁

【部分圖文】：

一維網(wǎng)格的Chebyshev結(jié)點示意圖

由式(3)～式(6),可以直觀地給出Chebyshev求導矩陣DN的結(jié)構(gòu),如圖2所示。這樣,便得到了用于求一階導數(shù)的Chebyshev求導矩陣DN。DN是用于求一階導數(shù)的Chebyshev求導矩陣,要得到用于求二階導數(shù)的Chebyshev求導矩陣,可以直接對DN進行平方處理,即

加入邊界條件后的正演方程組(16)含有N+1個方程和N+1個未知數(shù)。這有別于有限差分法或有限單元法正演所形成的系數(shù)矩陣[19-25],具有非稀疏形式,如圖3所示。求解線性方程組(16)便能得到地下介質(zhì)剖分節(jié)點處的電場值,再近似計算大地電磁輔助場值,即可得到一維地電模型的視電阻率和相位。3.2 大地電磁響應計算

【參考文獻】：
期刊論文
[1]基于譜元法的二維大地電磁數(shù)值模擬[J]. 方小姣,李長偉,劉劍.  礦產(chǎn)與地質(zhì). 2019(05)
[2]利用非均勻網(wǎng)格有限差分法模擬一維大地電磁響應[J]. 童孝忠,吳思洋,程東俊.  工程地球物理學報. 2018(02)
[3]頻率—空間域有限差分正演迭代法實現(xiàn)[J]. 徐夢婕.  工程地球物理學報. 2017(06)
[4]一種新的三維大地電磁積分方程正演方法[J]. 任政勇,陳超健,湯井田,周峰,陳煌,邱樂穩(wěn),胡雙貴.  地球物理學報. 2017(11)
[5]一種新型有限差分網(wǎng)格剖分方法在大地電磁一維正演中的應用[J]. 張輝,唐新功.  物探與化探. 2015(03)
[6]帶地形的二維各向異性大地電磁測深數(shù)值模擬[J]. 田紅軍,佟鐵鋼.  工程地球物理學報. 2015(02)
[7]探地雷達極化探測時域有限差分法模擬效果分析[J]. 李雙喜,曾昭發(fā),李靜,宋世明,張景,張寶松.  工程地球物理學報. 2014(04)
[8]二維彈性波方程頻率空間域有限差分正演研究[J]. 潘威,成景旺,閆立爽.  工程地球物理學報. 2013(06)
[9]有限單元法與有限差分法在MT一維正演模擬中的對比分析[J]. 王濤,柳建新,童孝忠,曹創(chuàng)華,譚神湘.  物探化探計算技術(shù). 2013(05)
[10]VTI介質(zhì)交錯網(wǎng)格FCT有限差分數(shù)值模擬[J]. 張省,何兵壽,王玉鳳.  工程地球物理學報. 2012(05)



本文編號：3284117