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

【文章來(lái)源】：工程地球物理學(xué)報(bào). 2020,17(06)

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

【部分圖文】：

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

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

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

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



本文編號(hào)：3284117