第 12 卷 第 9 期 2012 年 3 月 1671—1815( 2012) 09-2128-04 科 學(xué) 技 術(shù) 與 工 程 Science Technology and Engineering Vol. 12 No. 9 Mar． 2012  2012 Sci. Tech. Engrg. 均勻半空間中心回線瞬變電磁法正演研究 晏沖為 李文堯 ( 昆明理工大學(xué) 國(guó)土資源與工程學(xué)院，昆明 650093) 摘 要 通過對(duì)均勻半空間中心回線裝置瞬變電磁法表達(dá)式分析，可以為瞬變電磁法的反演提供理論正演數(shù)據(jù)、可以為瞬變電磁法的反演提供一定的方法，同時(shí)還可以為野外生產(chǎn)設(shè)計(jì)與室內(nèi)數(shù)據(jù)篩選與剔除提供一定的參考。通過一個(gè)反演實(shí)例說(shuō)明在晚期微小的二次電位誤差將對(duì)視電阻率的計(jì)算產(chǎn)生很大的影響。編制了 matlab 程序計(jì)算各種工作裝置和視電阻率理論正演二次電位衰減到噪聲影響程度高的準(zhǔn)確時(shí)間。對(duì)工程設(shè)計(jì)、評(píng)價(jià)視電阻計(jì)算的準(zhǔn)確度及剔除受噪聲影響嚴(yán)重的數(shù)據(jù)有一定的參考性。 關(guān)鍵詞 中心回線瞬變電磁法 二次電位 衰減時(shí)間中圖法分類號(hào) P631. 34; 文獻(xiàn)標(biāo)志碼 A 2012 年 1 月 4 日收到 地質(zhì)調(diào)查項(xiàng)目( 1212011121273) 資助第一作者簡(jiǎn)介: 晏沖為( 1984—) ，男，云南富源人，碩士研究生，研究方向: 瞬變電磁法理論與儀器。E-mail:89534861@ qq. com。 瞬變電磁法( TEM) ［1，2］中用接收線圈接收的感應(yīng)電壓隨時(shí)間呈早期信號(hào)幅值高且衰變速度快，晚期信號(hào)很弱衰減速度慢的特點(diǎn)，均勻半空間中心回線瞬變電磁法的感應(yīng)電壓表達(dá)式可以為理論研究提供無(wú)干擾的正演數(shù)據(jù)，同時(shí)還能為反應(yīng)的研究提供一定的方法，均勻半空間的瞬變電磁場(chǎng)與地下電阻率之間存在著復(fù)雜的隱函數(shù)關(guān)系，無(wú)法獲得全期視電阻率的解析表達(dá)式，目前通常使用各種近似定義和精確定義的方法來(lái)計(jì)算早期與晚期視電阻率或全期視電阻率［3］，全期視電阻率的計(jì)算很多學(xué)者提出了各種成熟的計(jì)算方法［3—7］，其中陳清禮［7］基于在 0. 08 ms ＜ t ＜1 s 的范圍內(nèi)，感應(yīng)電動(dòng)勢(shì)隨電阻率的增大而單調(diào)下降的特性提出了全期視電阻率的二分搜索算法，搜索出理論感應(yīng)電動(dòng)勢(shì)與實(shí)際觀測(cè)的電動(dòng)勢(shì)相符的視電阻率即是基于正演研究提供反演的一種方法的例子，現(xiàn)通過一個(gè)實(shí)例來(lái)討論晚期干擾噪聲對(duì)視電阻率計(jì)算的影響，并通過 mat-lab 編寫了一個(gè)計(jì)算各種工作裝置和視電阻率理論正演二次電位衰減到自主設(shè)定的噪聲影響程度高的準(zhǔn)確時(shí)間的程序，對(duì)于前期工程設(shè)計(jì)、后期數(shù)據(jù) 處理解釋時(shí)剔除不合理晚期數(shù)據(jù)，對(duì)提高解釋準(zhǔn)確度有一定的參考意義。 1 均勻半空間下感應(yīng)電動(dòng)勢(shì)的解析表達(dá) 中心回線或重疊回線方式是時(shí)間域瞬變電磁法大多采用的工作方法，1983 年美國(guó)地球物理學(xué)家P. Raab 和 F. Frisehknecht［8］推導(dǎo)出了中心回線裝置的感應(yīng)電壓表達(dá)式為 V( t) = qIπ3/2ρ L3 ［3Φ( z) － ( 3z + 2z3) Φ · ( z) ］u( t) ( 1) 其中 q 為接收線圈的有效面積( m2) ，L 為發(fā)送回線邊長(zhǎng)( m) ，I 為發(fā)送電流強(qiáng)度( A) ，它為滿足 I( t) = 0， t ＜ 0 I0， t ≥{ 0 ( 2) 的階躍電流。Φ( z) 為概率積分， Φ( z) = ∫ z 0 2槡 π exp( － t2) dt ( 3) Φ · ( z) = 2 槡 π exp( － z2) ( 4) z 為瞬變場(chǎng)參數(shù)，τ 為擴(kuò)散參數(shù)( m) z = 2 槡 πL/τ ( 5) τ = 2π 2ρ

總結(jié)

以上是生活随笔為你收集整理的均匀半空间视电阻率MATLAB代码,均匀半空间中心回线瞬变电磁法正演研究的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。