• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    淺部頻率域電磁勘探方法綜述

    2015-03-01 01:34:48湯井田任政勇周聰張林成原源肖曉
    地球物理學(xué)報(bào) 2015年8期
    關(guān)鍵詞:場源反演電磁

    湯井田, 任政勇, 周聰, 張林成, 原源, 肖曉

    1 中南大學(xué)有色金屬成礦預(yù)測與地質(zhì)環(huán)境監(jiān)測教育部重點(diǎn)實(shí)驗(yàn)室, 長沙 410083 2 中南大學(xué)地球科學(xué)與信息物理學(xué)院, 長沙 410083

    ?

    淺部頻率域電磁勘探方法綜述

    湯井田1,2, 任政勇1,2, 周聰1,2, 張林成1,2, 原源1,2, 肖曉1,2

    1 中南大學(xué)有色金屬成礦預(yù)測與地質(zhì)環(huán)境監(jiān)測教育部重點(diǎn)實(shí)驗(yàn)室, 長沙 410083 2 中南大學(xué)地球科學(xué)與信息物理學(xué)院, 長沙 410083

    適用于近地表(2000 m以內(nèi))勘探的頻率域電磁法主要有音頻大地電磁法(audio-frequency magnetotellurics, AMT), 無線電大地電磁法(radio-magnetotellurics, RMT), 可控源音頻大地電磁法(controlled source audio-frequency magnetotellurics, CSAMT),廣域電磁法(Wide Field Electromagnetic Method, WFEM).本文擬從最新的數(shù)據(jù)采集技術(shù)、數(shù)據(jù)處理技術(shù)、正反演算法、實(shí)例等四個方面,論述適用于淺部勘探的AMT, RMT, CSAMT和WFEM方法的國內(nèi)外最新進(jìn)展,總結(jié)目前AMT, RMT, CSAMT和WFEM方法遇到的困難,對潛在的發(fā)展方向提出建議.綜述表明: (1) 張量測量、多站陣列、多站疊加可提高AMT、RMT和 CSAMT數(shù)據(jù)的質(zhì)量.利用近區(qū)數(shù)據(jù)WFEM法可獲得良好的效果.國產(chǎn)與國外儀器在質(zhì)量方面的差距正在逐步縮小.(2)數(shù)學(xué)形態(tài)濾波技術(shù)、Hilbert-Huang變換等可有效分離出有用的數(shù)據(jù),局部畸變?nèi)匀皇秦酱鉀Q的難題,需要更為深入的研究.(3)矢量有限元與非結(jié)構(gòu)網(wǎng)格的出現(xiàn)大幅度提高了有限元處理復(fù)雜電磁問題模擬的精度與應(yīng)用范圍,成為目前電磁正演的首選工具.完全非線性反演算法仍然局限于1D、2D問題,共軛梯度法和高斯牛頓算法等為解決3D問題的發(fā)展趨勢.地質(zhì)約束的引入和多數(shù)據(jù)聯(lián)合反演可以減小反演的非唯一性.各向異性的反演為目前反演研究的熱點(diǎn)之一.(4)野外數(shù)據(jù)解釋的正確性嚴(yán)重依賴于對地下結(jié)構(gòu)先期的維性判別,在2D特性不明顯、3D特性明顯時,需要采用3D進(jìn)行反演解釋.

    近地表勘探; 頻率域電磁法; AMT; RMT; CSAMT; WFEM

    1 引言

    隨著人們對周圍生活環(huán)境認(rèn)識的進(jìn)一步加強(qiáng),以及對礦產(chǎn)資源需求的進(jìn)一步提高,要求更為可靠的探索地下淺部結(jié)構(gòu)的地球物理勘探方法.適用于淺部勘探的電磁勘探方法主要有音頻大地電磁法(AMT),射線大地電磁法(RMT),可控源音頻大地電磁法(CSAMT),廣域電磁法(WFEM).AMT和RMT的場源分別對應(yīng)于磁暴和雷電現(xiàn)象、潛艇信號.由于距離足夠遠(yuǎn),磁暴、雷電和潛艇天線產(chǎn)生的電磁波可近似為平面波.CSAMT的場源為人工電偶極子或者磁偶極子,數(shù)據(jù)采集一般在遠(yuǎn)區(qū),場源可近似為平面波.因此,AMT, RMT和CSAMT 可歸納為基于平面波場源的電磁勘探方法.WFEM的場源與CSAMT類似,但測量不局限于遠(yuǎn)區(qū).

    AMT, RMT,CSAMT和WFEM數(shù)據(jù)采集的頻率范圍略有不同.AMT的工作頻率一般為0.1 Hz~100 kHz; RMT的工作頻率一般為10 kHz~250 kHz; CSAMT和WFEM的工作頻率一般為0.1 Hz~10 kHz.在地表,我們一般采集隨時間變化的電磁場信號E(t)和H(t)(t為時間變量),通過信號處理獲得頻率域電磁場E(ω)和H(ω)(ω為頻率變量).電磁場信號一般為地下電導(dǎo)率(σ)的函數(shù),電場與磁場的比值Z(ω,σ)=E(ω)/H(ω)(通常稱Z為阻抗)建立了頻率與地下電性結(jié)構(gòu)的關(guān)系.利用地表采集的阻抗Z,反演可獲得地下電導(dǎo)率分布.基于廣域電阻率概念的廣域電磁法在地質(zhì)勘探中獲得了越來越多的成功應(yīng)用.

    本文擬從數(shù)據(jù)采集和處理技術(shù)、反演算法和實(shí)例等四個方面來論述淺部勘探AMT,RMT, CSAMT和WFEM方法的國內(nèi)外最新進(jìn)展,給出AMT,RMT, CSAMT和WFEM勘探方法目前遇到的困難,并指明潛在的發(fā)展方向.

    2 數(shù)據(jù)采集技術(shù)

    數(shù)據(jù)采集技術(shù)的發(fā)展主要體現(xiàn)在數(shù)據(jù)采集方案和硬件設(shè)備兩方面.

    2.1 采集方案

    AMT, RMT和CSAMT的一個共同點(diǎn)為在源的遠(yuǎn)區(qū)進(jìn)行數(shù)據(jù)采集.在遠(yuǎn)區(qū),入射源場可以近似為平面波.遠(yuǎn)區(qū)可采集的數(shù)據(jù)為水平電場、水平磁場和垂直磁場.由于空氣的電導(dǎo)率一般為0,地表的垂直電場分量為0,因此在頻率域電磁法一般不采集地表下的垂直電場.如果測站分布設(shè)計(jì)相對較密,測點(diǎn)距較小,高密度的數(shù)據(jù)足以保證較高分辨率模型的獲取,用于標(biāo)識地下橫向電性變化的垂直磁場的測量也可省略.

    由于天然場源信號弱,極化方向隨機(jī),極易受各類噪聲干擾,這些特征使得高質(zhì)量AMT數(shù)據(jù)的獲取極為困難.為提高數(shù)據(jù)質(zhì)量,合理的設(shè)計(jì)采集方案是首要任務(wù).在天然源電磁法觀測方案的發(fā)展歷程中,有三個關(guān)鍵節(jié)點(diǎn)值得注意.一是張量觀測方案的提出(Sims et al., 1971),通過定義水平電場與磁場分量的二階阻抗張量來估計(jì)地下電性分布.該方案使得天然場的觀測真正走向?qū)嵱?,其典型的采集裝置如圖1所示.二是遠(yuǎn)參考法(Gamble et al., 1979)的提出,通過在噪聲較低的區(qū)域布置參考測站,有效壓制了測區(qū)的噪聲干擾,改善了數(shù)據(jù)質(zhì)量.該方案自提出后即成為天然場源特別是低頻信號觀測的標(biāo)準(zhǔn)程序.三是電磁陣列剖面(Electromagnetic Array Profiling,EMAP)(Torres-Verdín and Bostick, 1992)觀測的提出,為壓制局部異常體對觀測數(shù)據(jù)的畸變影響提供了解決思路.該方案及其處理程序也成為2D電磁剖面處理中重要的步驟.學(xué)者們現(xiàn)已將其擴(kuò)展至三維空間觀測層面,如小面元(何繼善等, 2010c)、多站疊加(Jiang and Xu, 2013)等觀測方式.

    人工源電磁法(CSAMT/RMT)由于源的引入提高了數(shù)據(jù)信噪比,且由于場源極化方向固定,使得矢量和標(biāo)量觀測兩種方案均可使用.實(shí)際施工中,為提高生產(chǎn)效率,CSAET成為一種常用的觀測方案(用一個場源在一個方向上只做電場觀測的標(biāo)量方案).幾種典型的觀測方案如圖2所示.

    然而,CSAET等觀測方案顯然并不適用于地質(zhì)條件復(fù)雜的地區(qū).針對此,學(xué)者們提出了許多改進(jìn)方案,這些改進(jìn)主要體現(xiàn)在發(fā)送端.綜合來看,可分為以下幾個方面.一是改進(jìn)信號源的發(fā)送形式,如人工源極低頻電磁波技術(shù)(趙國澤等,2010; 卓賢軍等,2011; 卓賢軍和陸建勛,2010),在高電阻率地區(qū)建立一個固定的大功率發(fā)射臺,用人工方法發(fā)射0.1~300 Hz超低/極低頻電磁信號.一些學(xué)者討論了考慮電離層的影響時該方法的響應(yīng)(底青云等,2008; 付長民等,2012; 陳小斌和趙國澤,2009; 李帝銓等,2011; 徐志鋒和吳小平,2010).二是改進(jìn)發(fā)送信號的編碼形式,主要以何繼善(2010c)提出的偽隨機(jī)信號編碼為代表,羅維斌等(2012)研究了逆重復(fù)M序列偽隨機(jī)信號編碼方式.三是改進(jìn)發(fā)送源的布設(shè)方法,如“L”型發(fā)射源(王顯祥等,2014)、十字型源(孟慶奎,2013)以及多偶極發(fā)射裝置(雷達(dá)等,2014)等.

    圖1 張量觀測方案(Christopherson, 2012; Ren, 2012)Fig.1 Configuration for tensor measurements (Christopherson, 2012; Ren, 2012)

    圖2 CSAMT法的張量、矢量和標(biāo)量觀測方案(湯井田和何繼善,2005)(a)(b)用兩個場源,且每個場源做5次的全張量測量; (c)每個場源觀測3個分量的部分張量測量; (d)單一場源5次觀測的矢量測量;(e)單一場源2次觀測的標(biāo)量測量;(f)用一個場源在一個方向上只做電場觀測的標(biāo)量CSAET法,在少量點(diǎn)上測磁場用以把電場數(shù)據(jù)轉(zhuǎn)換為近似的電阻率.Fig.2 The configuration of CSAMT tensor, vector and scalar measurements (Tang and He, 2005)(a) and (b) Two sources, a set of five tensor data is measured for each source. (c) A set of three tensor data is measured for each source. (d) A set of five vector data is measured for each source. (e) A set of two scalar data is measured for each source. (f) Only the electric data is measured along the profile for one source. At several measurement sites, the magnetic field is also used to transfer the measured electric field to the apparent resistivity.

    值得一提的是,何繼善(2010a,2010b)提出了一種新的頻率域方法,即廣域電磁法(WFEM).廣域電磁法具有獨(dú)特的觀測方式,其僅僅測量單分量電場,避免了磁場干擾對測量數(shù)據(jù)的影響.另外,在發(fā)送端采用偽隨機(jī)信號編碼,并且可同時發(fā)送包含多頻率的信號,提高發(fā)送效率;在接收端,進(jìn)行單分量陣列觀測,同時接收多個頻率的地電響應(yīng)信號,提高精度和速度.

    2.2 儀器設(shè)備

    電磁場觀測設(shè)備分為發(fā)送端與接收端兩大部分,對于天然場源電磁法,發(fā)送端部分可省略.發(fā)送端主要包括發(fā)電機(jī)、逆變器及信號發(fā)送控制器等組成部分,接收端主要包括電極、磁場傳感器、GPS、電源、信號放大器及采集控制器等組成部分.隨著技術(shù)的進(jìn)步,儀器設(shè)備也在不斷更新?lián)Q代.綜合來看,未來電磁勘探的觀測技術(shù)及儀器設(shè)備的趨勢主要為以下幾個方面:

    (1)采集系統(tǒng)的集成化,實(shí)現(xiàn)一套系統(tǒng)的多種電磁數(shù)據(jù)采集.如加拿大鳳凰公司推出的V-8系統(tǒng),可以進(jìn)行AMT/CSAMT方法觀測;美國Zonge公司的GDP-32II系統(tǒng),可以進(jìn)行AMT/CSAMT/觀測;德國Metronix公司的GMS-07e系統(tǒng),可以進(jìn)行AMT/CSAMT/RMT觀測等.

    (2)全三維采集,數(shù)據(jù)記錄由單測線的剖面測量逐步轉(zhuǎn)變?yōu)槿臻g陣列三維采集;如Geometrics 公司的Geode EM 3D系統(tǒng)支持3D網(wǎng)絡(luò)式布設(shè),最多可測量240道數(shù)據(jù);Quantec Geoscience公司的Orion 3D系統(tǒng)支持開展全3D勘探等.

    (3)頻率范圍的拓展,更高精度的傳感器使得采集頻率得以向更低和更高的范圍拓展.德國Metronix公司的GMS-07e系統(tǒng)測量頻帶可達(dá)0.00001 Hz~250 kHz;美國Zonge公司的GDP-32II系統(tǒng)測量頻帶范圍可達(dá)0.000121 Hz~8 kHz.此外,更大的發(fā)送電流、更加輕便的裝備、更加持久的續(xù)航能力以及更加快捷的實(shí)時處理能力也是儀器發(fā)展一直努力的目標(biāo).

    與國外相對成熟的商業(yè)儀器相比,國內(nèi)電磁法儀器的發(fā)展相對滯后,主要停留在試驗(yàn)階段.我國多個科研單位和高校進(jìn)行了頻率域電磁法儀器的研制,如吉林大學(xué)(程德福等,2004)、中國地質(zhì)科學(xué)院物化探研究所(林品榮等,2007,2010)、中國科學(xué)院地質(zhì)與地球物理研究所(底青云等,2013; 真齊輝等,2013).

    基于全區(qū)視電阻率概念(湯井田和何繼善, 1994a), 中南大學(xué)成功研制出了JSGY-2廣域電磁儀(何繼善, 2010c),其工作頻率范圍為3/256 Hz~8192 Hz.該儀器可在遠(yuǎn)區(qū)至近區(qū)的廣大區(qū)域進(jìn)行觀測,拓展了觀測范圍.JSGY-2廣域電磁儀已經(jīng)廣泛應(yīng)用于油氣藏、頁巖氣、地?zé)?、煤田采空區(qū)、金屬礦、水資源等探測領(lǐng)域.因此,國產(chǎn)與國外儀器在質(zhì)量方面的差距正在逐步縮小,差距主要體現(xiàn)在制造標(biāo)準(zhǔn)、工藝水平和商業(yè)化水平等方面.

    3 數(shù)據(jù)處理技術(shù)

    3.1 天然場源AMT方法

    AMT方法的場源為天然雷電,記錄的數(shù)據(jù)為不同采樣率的時間域信號.為獲得可解釋的阻抗視電阻率及相位,需要進(jìn)行時間域數(shù)據(jù)處理、時頻轉(zhuǎn)換、阻抗估計(jì)及局部畸變校正等幾個步驟.提高數(shù)據(jù)質(zhì)量的努力也主要從這四個方面展開.

    3.1.1 時間域數(shù)據(jù)處理

    對于時間域數(shù)據(jù)處理,無論是數(shù)據(jù)刪選還是信噪分離,信噪識別是成敗關(guān)鍵.Weckmann 等(2005)提出可利用能量、阻抗在復(fù)平面上的分布、阻抗估計(jì)誤差、相干度和極化方向等多種參數(shù)進(jìn)行綜合判斷.湯井田等(2012b)指出,強(qiáng)干擾數(shù)據(jù)的時間域特征為幅值大,形態(tài)規(guī)則.Escalas等(2013)使用小波分析了天然場源和人工源不同的時頻域極化特征,并利用其差異進(jìn)行信噪識別.

    3.1.2 時頻轉(zhuǎn)換

    時頻轉(zhuǎn)換是將時間域信號轉(zhuǎn)換到頻率域的過程.快速傅里葉變換(Fast Fourier Transform,F(xiàn)FT)(Cooley and Tukey,1965; Cooley et al.,1969)是最為常用的方法,其優(yōu)點(diǎn)是計(jì)算穩(wěn)定,易于實(shí)現(xiàn),因此已成為天然場源電磁法數(shù)據(jù)處理的經(jīng)典程序之一.但它也存在明顯的不足,如時窗函數(shù)固定、對信號時間定位能力差、頻率混淆效應(yīng)與泄露效應(yīng)無法根除等,因此,一些不同的譜估計(jì)方法相繼被引入電磁數(shù)據(jù)的處理中,如高階統(tǒng)計(jì)(王書明和王家映,2004)、最大熵譜分析法(Tzanis and Beamish,1989)、S變換(陳海燕等,2012;景建恩等,2012)等.這些方法對于特定問題表現(xiàn)出了優(yōu)于FFT的性能,但尚未成為標(biāo)準(zhǔn)處理程序.小波變換是一種具有多分辨率、多尺度的數(shù)學(xué)方法,具有良好的時頻局域性和自適應(yīng)性,可對信號的不同頻率成分進(jìn)行分析,在信號處理領(lǐng)域已多有應(yīng)用,眾多學(xué)者研究了其對電磁數(shù)據(jù)的處理效果(Zhang and Paulson,1997; 何蘭芳等,1999; 徐義賢和王家映,2000;Trad and Travassos,2000; 劉宏等,2004;嚴(yán)家斌和劉貴忠,2007),證明了其有效性,但該方法依賴于小波基函數(shù)的選取,這也成為其難點(diǎn).天然電磁場信號是非線性、非平穩(wěn)信號(Chant and Hastie,1992),噪聲同樣常具非平穩(wěn)性(Banks,1998),因此,非線性、非平穩(wěn)信號分析工具是最佳選擇.湯井田等(2008)將Hilbert-Huang變換引入電磁數(shù)據(jù)處理中,基于經(jīng)驗(yàn)?zāi)B(tài)分解(Empirical Mode Decomposition, EMD)與希爾伯特譜分析,對大地電磁數(shù)據(jù)進(jìn)行處理,證明了該方法的有效性,其后該方法被進(jìn)一步完善和推廣(Cai et al.,2009; Chen et al.,2010; Cai,2013,2014; Neukirch and Garcia,2014),顯示出良好的應(yīng)用前景(見圖3).

    3.1.3 阻抗估計(jì)

    張量阻抗估計(jì)是獲取高質(zhì)量轉(zhuǎn)換函數(shù)、阻抗視電阻率及相位等解釋參數(shù)的重要步驟.針對不同的噪聲類型、含噪程度,主要的張量阻抗估計(jì)方法經(jīng)歷了從最小二乘估計(jì)到穩(wěn)健阻抗估計(jì),從引入單參考道到多道數(shù)據(jù)的利用等階段.

    最小二乘(Sims et al.,1971)是估計(jì)阻抗張量的經(jīng)典方法.在噪聲環(huán)境相對平靜或以隨機(jī)噪聲為主的地區(qū),該方法是有效的,但是當(dāng)存在有色噪聲或相關(guān)噪聲時,最小二乘估計(jì)結(jié)果往往會出現(xiàn)嚴(yán)重偏離.如何在含噪情況下獲取高質(zhì)量的頻域阻抗參數(shù)成為電磁研究者們關(guān)心的問題.

    穩(wěn)健阻抗估計(jì)是一種主流的提高阻抗估計(jì)質(zhì)量的方法. 各種穩(wěn)健阻抗估計(jì)都是通過自適應(yīng)地賦權(quán)來降低奇異值在阻抗估計(jì)中的地位.一般來說,穩(wěn)健估計(jì)的成功依賴于整個數(shù)據(jù)段中性態(tài)良好的數(shù)據(jù)占主要地位,當(dāng)噪聲貫穿大部分甚至整個觀測時段時,穩(wěn)健估計(jì)方法就失效了(Egbert and Booker,1986; Chave et al. ,1987; Chave and Thomson,1989; Egbert and Livelybrooks,1996; Chave and Thomson,2003; Smirnov,2003).Jones等(1989)、Sutarno(2008)及湯井田等(2013a)分別對穩(wěn)健阻抗估計(jì)方法進(jìn)行了綜述和對比.

    圖3 HHT矯正振動干擾的實(shí)例(湯井田等,2008)(a) Hx的EMD 分解殘差; (b) Hx原始信號; (c) 振動干擾矯正后的Hx信號. Fig.3 A field example of the vibration interference which is corrected by the HHT method (Tang et al., 2008)(a) EMD decomposition of the Hx component; (b) The original Hx signal; (c) The corrected Hx signal.

    遠(yuǎn)參考道估計(jì)是提高阻抗估計(jì)質(zhì)量的另一種可靠方法.自Gamble 等(1979)首次提出遠(yuǎn)參考方法以來,該方法已成為天然場源電磁法的標(biāo)準(zhǔn)程序,并取得了長足的發(fā)展.利用測點(diǎn)磁場與參考點(diǎn)磁場之間的相關(guān)性,可以實(shí)現(xiàn)含噪數(shù)據(jù)的剔選(Ritter et al.,1998; Varentsov,2006; Sokolova et al.,2005).Larsen等(1996)提出了利用參考道實(shí)現(xiàn)信噪分離的思路,當(dāng)遠(yuǎn)參道中不含噪聲時,可以估計(jì)出測站磁場中的天然場信號項(xiàng),進(jìn)而獲得平面波阻抗的估計(jì).Oettinger 等 (2001)進(jìn)一步證明,當(dāng)參考道中含噪時,引入第二個遠(yuǎn)參考道可以得到無偏的平面波阻抗估計(jì).

    多道數(shù)據(jù)的利用是學(xué)者們努力的方向之一. 利用測站間電磁場分量的相互關(guān)系對含噪時段數(shù)據(jù)進(jìn)行預(yù)測是一種常見的思路,如大地電流-大地電磁(Telluric-Magnetotelluric,T-MT)方案(García and Jones, 2005), “偽遠(yuǎn)參考”(Muoz and Ritter,2013)及組合站間轉(zhuǎn)換函數(shù)(Estimation of Local transfer-functions by Combining Interstation Transfer-functions,ELICIT)(Campanyà et al., 2014)等一批基于測站間“轉(zhuǎn)換函數(shù)”的采集處理方案.這些方案不僅改善了MT的數(shù)據(jù)質(zhì)量,而且拓展了數(shù)據(jù)采集的思路.然而,參考道的數(shù)據(jù)質(zhì)量決定了此類方法的適用范圍,另一些學(xué)者提出了算法方面的改進(jìn),如Kappler (2012)提出了IARWR(Intersite Activity Ratio Wiener-filter Replacement)時間域去噪方法.王輝等(2014) 提出利用多道同步依賴關(guān)系進(jìn)行阻抗估計(jì)的方法,Cui 等(2013)引入獨(dú)立分量分析(Independent Component Analysis,ICA)的方法,Varentsov 等(2003)提出利用多道遠(yuǎn)參考聯(lián)合使用的方法,另外,Egbert和Booker(1989)、Egbert(1997,2002)系統(tǒng)闡述了利用陣列數(shù)據(jù)進(jìn)行信噪分離的阻抗估計(jì)方法,在其基礎(chǔ)上,Smirnov 和Egbert (2012)發(fā)展了陣列數(shù)據(jù)的主成分分析穩(wěn)健估計(jì)方法(Principal Component Analysis,PCA).

    此外,利用多頻點(diǎn)數(shù)據(jù)進(jìn)行光滑處理也是提高數(shù)據(jù)質(zhì)量的有效手段之一.如利用視電阻率、相位數(shù)據(jù)的轉(zhuǎn)換關(guān)系(Weidelt,1972),可以實(shí)現(xiàn)阻抗相位參數(shù)對畸變視電阻率曲線的校正(楊生等,2001).Parker提出的Dplus及Rhoplus反演方法(Parker,1980; Parker et al.,1981; Parker et al.,1996)是這類處理的優(yōu)秀代表.該方法考慮了阻抗視電阻率、相位數(shù)據(jù)的關(guān)聯(lián)性以及數(shù)據(jù)在頻率域的連續(xù)性,其擬合數(shù)據(jù)來自于所謂D+或R+模型的響應(yīng),具有明確的物理意義,在頻域數(shù)據(jù)的一維分析與擬合中應(yīng)用廣泛(Beamish and Travassos,1992; Fischer and Weibel,1991; Spratt et al.,2005).

    3.1.4 局部畸變校正

    當(dāng)?shù)乇泶嬖谌S局部電性不均勻體時,會因積累電荷而使測量電場產(chǎn)生局部畸變,進(jìn)而影響計(jì)算所得的阻抗數(shù)據(jù).這種局部畸變的影響表現(xiàn)在兩個方面(Bahr,1988,1991; Berdichevsky and Dimitriev,1976):一是電流型畸變(galvanic distortion),以準(zhǔn)靜態(tài)形式影響觀測的電場,使阻抗視電阻率數(shù)據(jù)在所有頻點(diǎn)上產(chǎn)生所謂“靜態(tài)位移”,或稱為“靜態(tài)效應(yīng)”;二是感應(yīng)型畸變(inductive distortion),會同時影響阻抗視電阻率、相位數(shù)據(jù).對于實(shí)測數(shù)據(jù)而言,局部畸變往往是不可避免的,即使是信噪比很高的人工源電磁法數(shù)據(jù),因此,對于局部畸變的校正研究引起了學(xué)者們的極大興趣. 阻抗張量分解技術(shù)將觀測值視為區(qū)域阻抗與局部畸變兩部分,具有較為清晰的理論基礎(chǔ),主要的分解算法有Swift 方法(Swift,1967)、Bahr分解( Bahr,1988; Bahr,1991; Prácser and Szarka,1999)、GB 分解(Groom and Bailey,1989,1991; Jones and Groom,1993)、Chave-Smith電磁場全畸變分解法(Chave and Smith, 1994; Chave et al., 2004)、多測點(diǎn)-多頻點(diǎn)分解方法(陳小斌等,2014; McNeice and Jones,2001)以及相位張量分解(Caldwell et al.,2004; Weaver et al.,2006; Booker,2014)等. 針對靜態(tài)效應(yīng)問題,除了上述阻抗張量分解方法外,還發(fā)展了一系列校正方法, 這些方法可歸為三大類,即利用相位或磁場分量測量結(jié)果(閻述和陳明生,1996)、利用空間域(Torres-Verdín and Bostick,1992)或波數(shù)域(宋守根等,1995)數(shù)字濾波方法和聯(lián)合反演方法(Sternberg et al.,1988)等.

    值得指出,目前實(shí)際應(yīng)用較多的局部畸變校正方法主要為阻抗張量分解與空間濾波算法.

    3.2 人工源電磁法

    CSAMT、WFEM和RMT為人工源電磁法,其處理方式與AMT法有許多相似之處,但也存在明顯的差異,如發(fā)射信號可以固定在特定的頻率或組合頻率,提高該頻率的信號能量,在接收端可以獲得很高的信噪比,進(jìn)而得到高質(zhì)量的數(shù)據(jù),同時因場源極化方向固定,在不同的空間位置不同方向上的信號能量差異巨大,當(dāng)采集區(qū)域限制在特定范圍內(nèi)時,可以簡化觀測方案及波阻抗的計(jì)算方法.與天然場源方法相比較,人工源電磁法的數(shù)據(jù)處理較為簡單,但因?yàn)樵吹囊?,也產(chǎn)生了一系列與源有關(guān)的問題(何繼善,1991; 湯井田和何繼善,2005; 高文,1991).

    非平面波效應(yīng)是人工場源效應(yīng)之一,其主要特征是會造成所謂“過渡區(qū)”和“近區(qū)”卡尼亞視電阻率、相位數(shù)據(jù)畸變.對此,首先需要討論人工源條件下,觀測場區(qū)的劃分問題.何繼善(1991)、湯井田和何繼善(2005)對這一問題進(jìn)行了詳細(xì)的論述,指出場區(qū)的劃分與地下電阻率、觀測頻率以及發(fā)收距等相關(guān),可用綜合參數(shù)“電距離”來進(jìn)行衡量.岳瑞永和徐義賢(2004)、陳明生和閆述(2005)、林威(2009)、程輝和李帝銓(2014)分別給出了不同的定量場區(qū)劃分方法,為畸變數(shù)據(jù)的識別提供了依據(jù).

    在過渡區(qū)及近區(qū),卡尼亞視電阻率及相位的畸變常常會帶來錯誤的認(rèn)識,通過定義適當(dāng)?shù)碾娮杪蕝?shù)以合理利用過渡區(qū)數(shù)據(jù)是應(yīng)對措施之一,如Das(1995)、曹昌祺(1978)、黃皓平和樸化榮(1992)、殷長春和樸化榮(1991)、詹少全等(2011)以及馮兵等(2013)均做過不同的嘗試.實(shí)際工作中,合理設(shè)計(jì)采集方案,盡量保證數(shù)據(jù)處于遠(yuǎn)區(qū)是避免非平面波效應(yīng)影響的最佳策略.湯井田等(2013b)、王若等(2004)討論了不同條件下適宜的觀測方案. 陳衛(wèi)營和薛國強(qiáng)(2014)則給出了適用于全場域的有效趨膚深度快速估算公式.

    陰影和場源附加效應(yīng)是另一類由場源所引起的問題. 目前對該問題的研究多處于數(shù)值模擬階段(Boschetto and Hohmann,1991;閆述和陳明生,2004; 王若等,2009; 湯井田等,2012c).陳明生和閆述(2005)通過實(shí)例表明,陰影和場源附加效應(yīng)是由地層波所攜帶的源與觀測點(diǎn)之間、源下方真實(shí)地質(zhì)情況的信息決定,是可以加以利用的.有效的陰影和場源附加效應(yīng)的處理方案仍有待后續(xù)研究,是今后的一個研究方向.

    湯井田和何繼善(1993,1994a,1994b)從電磁場統(tǒng)一的角度,系統(tǒng)地討論了頻率電磁測深中幾種等效全區(qū)視電阻率的原理、定義和具體算法,并提出形式更為簡潔的Ey視電阻率(湯井田等,2011).基于等效全區(qū)視電阻率概念(湯井田和何繼善, 1994a)的WFEM法,可不局限于遠(yuǎn)區(qū),從而拓展了觀測范圍.基于單分量的測量, 廣域電磁法有效壓制電磁場干擾對測量數(shù)據(jù)的影響,從而提高抗干擾能力.

    4 正反演算法發(fā)展

    4.1 正演

    數(shù)值模擬方法具有普遍適用性,不僅僅適用于本論文研究的 AMT,RMT, CSAMT和WFEM問題,還適用于其他問題,如MT,CSEM, 直流電法等.為了更好地為讀者敘述數(shù)值模擬技術(shù)的發(fā)展趨勢,下面的闡述著眼于其在整個地電磁場模擬的發(fā)展.通常,數(shù)值模擬方法可分為四類:體積積分法(Raiche, 1974; Hohmann, 1975; Weidelt, 1975; Avdeev et al., 2002; Zhdanov et al., 2006),邊界積分法(Parry and Ward, 1971; Doherty, 1988; Liu and Becker, 1992; Xu et al., 1997; Liu and Lamontagne, 1998; Ren et al., 2012),有限差分法(Mackie et al., 1994; Aprea et al., 1997; Smith, 1996; Haber et al., 2000; Newman and Alumbaugh, 2002; Hou et al., 2006; Streich, 2009)和有限單元法(Badea et al., 2001; Mitsuhata and Uchida, 2004; Key and Weiss, 2006; Nam et al., 2007).對于有限差分、體積分方程法的優(yōu)缺點(diǎn),請查看已有的綜述文章Avdeev (2005), B?erner (2010) 和Everett(2011).需要特別提出的,Penz等(2013)提出了利用有限差分解決復(fù)雜地電模型的新思路,擴(kuò)展了傳統(tǒng)有限差分的應(yīng)用范圍. 本節(jié)只闡述適用于任意復(fù)雜模型的有限單元法.

    從Coggon(1971)首次將有限單元法應(yīng)用到電磁模擬以來,該領(lǐng)域有了突飛猛進(jìn)的發(fā)展,技術(shù)的進(jìn)步可大概歸納為微分控制方程、有限元空間、網(wǎng)格離散化技術(shù)、線性求解器等四個方面.電磁模擬可基于電場方程(Nam et al., 2007; Farquharson and Miensopust, 2011)、磁場方程(Siripunvaraporn et al., 2002;Franke et al., 2008)和不同種類的勢場方程. 對于低頻問題,由于空氣的電導(dǎo)率為0,磁場方程在空氣中未定義.因此,磁場方程一般很少用于求解低頻地球物理電磁場問題.低頻情況下,電場方程在空氣中變得非常的病態(tài),從而導(dǎo)致高條件數(shù)的離散系統(tǒng)矩陣.高條件數(shù)一方面降低了求解的精度,另外一方面增加了迭代求解器的迭代次數(shù),從而增加了求解時間.常用的勢場方程包含H-Ω式,A-Ω式 (Haber et al., 2000; Badea et al., 2001)和T-Ω式(Mitsuhata and Uchida, 2004).勢場方程導(dǎo)致的離散系統(tǒng)矩陣的條件數(shù)一般較低,從而降低了對求解器的要求, 在處理電導(dǎo)率不連續(xù)界面上的邊界條件時,勢場方程存在一定的復(fù)雜性.

    有限元空間常常包含節(jié)點(diǎn)型基函數(shù)和向量基函數(shù)(Nédélec, 1986; Farquharson and Miensopust, 2011).節(jié)點(diǎn)型基函數(shù)構(gòu)成連續(xù)的有限元空間,要求電場或者磁場在地下為連續(xù)函數(shù).在電性不連續(xù)界面,如電導(dǎo)率跳躍界面處,電場在跳躍界面的法向分量不連續(xù).因此節(jié)點(diǎn)型基函數(shù)構(gòu)造出的電磁場有限元解,強(qiáng)加了法向連續(xù)性條件,從而使得有限元解不準(zhǔn)確,或者偏離真實(shí)解.向量基函數(shù)可解決節(jié)點(diǎn)型基函數(shù)產(chǎn)生的上述問題. 向量基函數(shù)構(gòu)造的有限元解只要求電場的切向連續(xù)性,保留電場的法向不連續(xù)性;另外一方面,能夠保證單元內(nèi)場的散度為零.因此,向量基函數(shù)構(gòu)造的有限元解不存在偽解.由于其優(yōu)越的性能,向量基函數(shù)已逐步代替節(jié)點(diǎn)型基函數(shù),成為地球物理電磁場數(shù)值模擬的標(biāo)準(zhǔn)方式.

    網(wǎng)格離散化技術(shù)經(jīng)歷了從結(jié)構(gòu)化網(wǎng)格到非結(jié)構(gòu)化網(wǎng)格的過渡(Key and Weiss, 2006; Franke et al., 2008).早期的電磁場模擬一般采用簡單的結(jié)構(gòu)化網(wǎng)格,如四邊形、六面體.隨著研究的深入,人們探究復(fù)雜地電磁場模型欲望的加強(qiáng),結(jié)構(gòu)化網(wǎng)格的缺點(diǎn)越來越突出.結(jié)構(gòu)化網(wǎng)格很難真實(shí)地模擬復(fù)雜的地下結(jié)構(gòu),如任意復(fù)雜的地形情況. 隨著計(jì)算幾何技術(shù)的突飛猛進(jìn),產(chǎn)生了非??焖?、可靠地用三角形和四面體來任意精確逼近復(fù)雜地電模型的技術(shù). 有限元本身對單元類型沒有要求,自然地,基于三角形與四面體的非結(jié)構(gòu)化網(wǎng)格越來越多地應(yīng)用于地電磁場的有限元模擬. 非結(jié)構(gòu)化網(wǎng)格能夠非常容易地實(shí)現(xiàn)網(wǎng)格的自適應(yīng)加密,因此,非結(jié)構(gòu)化有限元模擬越來越重視網(wǎng)格自適應(yīng)加密技術(shù),生成高精度的數(shù)值結(jié)果.典型的應(yīng)用為: Key和Weiss (2006)應(yīng)用自適應(yīng)三角形網(wǎng)格研究了復(fù)雜的2D 大地電磁問題. Li和Key (2007)應(yīng)用自適應(yīng)三角形網(wǎng)格研究了復(fù)雜的2D 海洋CSEM問題(圖4). Ren 和Tang(2010)應(yīng)用自適應(yīng)四面體網(wǎng)格研究了3D DC問題. Schwarzbach等(2011)應(yīng)用自適應(yīng)四面體網(wǎng)格研究了3D CSEM問題.Ren等(2013) 應(yīng)用自適應(yīng)四面體網(wǎng)格研究了3D MT和RMT 問題(圖5).目前, 非結(jié)構(gòu)網(wǎng)格在有限元模擬中逐漸成為主流.

    線性求解器方面也展現(xiàn)出較大的進(jìn)步. 電磁場問題經(jīng)有限元離散化后生成大規(guī)模的稀疏矩陣.由于計(jì)算機(jī)內(nèi)存的限制以及直接法的落后,早期的研究者一般采用迭代法來完成線性方程組的求解.然而,求解低頻電磁場問題時,線性方程的條件數(shù)往往較大,導(dǎo)致過多的迭代求解步驟,從而增加了計(jì)算時間.因此,如何設(shè)計(jì)低成本、高效的迭代預(yù)處理器成為了迭代法的發(fā)展方向.值得提及的是,Grayver 和Burg (2014)提出了一種基于輔助空間的高效預(yù)處理,能夠有效地減小迭代次數(shù).目前,直接求解器方面獲得了非常大的突破,產(chǎn)生了有效的直接法求解器,如基于LU矩陣分解法和OpenMP的PARDISO求解器(Schenk和G?rtner, 2004),基于波前法和MPI的MUMPS求解器(Amestoy et al., 2002).Streich (2009)測試了MUMPS求解器的優(yōu)越性能,求解了超過百萬未知數(shù)的大型地電磁場問題.地電問題一般為多源問題,即含多個右端項(xiàng).直接法能夠快速地處理多個右端項(xiàng)問題,因此獲得了越來越多的重視,有逐步取代迭代求解器的發(fā)展趨勢.另外為了求解超大規(guī)模的問題,基于MPI并行的求解技術(shù)也獲得了大量的嘗試和應(yīng)用.Puzyrev 等(2013) 利用區(qū)域分解和MPI求解了大規(guī)模的CSEM 3D 問題(見圖6);Ren等(2014)成功應(yīng)用了一種新的FEDI-DP區(qū)域分解技術(shù)計(jì)算了大規(guī)模MT問題.

    4.2 反演

    反演算法是一個優(yōu)化過程,目前常用的反演算法可歸納為兩類:非線性反演算法和線性迭代算法.

    非線性反演算法包括人工神經(jīng)網(wǎng)格法、模擬退火法、蒙特卡羅法等.非線性反演算法最大的優(yōu)點(diǎn)是可避開局部極大、極小值,搜索出全局最優(yōu)模型.然而,非線性反演算法需要大量的樣本,根據(jù)每個樣本的測試結(jié)果,進(jìn)行有方向性的樣本篩選或樣本升級,搜索出符合野外實(shí)測數(shù)據(jù)的樣本. 在模型搜索過程中,需要大量的正演計(jì)算,由于每一次正演需要消耗大量的計(jì)算時間,非線性反演算法目前僅僅適合1D或者小規(guī)模2D反演問題.

    線性迭代算法主要包括高斯-牛頓算法及其變種、共軛梯度算法及其變種.高斯-牛頓算法及其變種是目前流行的反演方法,廣泛應(yīng)用于MT反演、CSEM反演、DC反演等.靈敏度矩陣是制約反演速度的關(guān)鍵因素,由于靈敏度矩陣是密實(shí)矩陣,對于3D大規(guī)模問題,計(jì)算和存儲靈敏度矩陣是困難的.共軛梯度算法及其變種通過減小目標(biāo)函數(shù)梯度來達(dá)到反演收斂的目的,由于不需要引入靈敏度矩陣,共軛梯度算法及其變種的計(jì)算消耗相對小.關(guān)于反演算法在地球物理電磁法的應(yīng)用,請查看已有的綜述文章Avdeev (2005), B?erner (2010) 和Everett(2011).

    總體來說,對于2D和3D問題,可行的反演算法為線性迭代算法.由于上述最優(yōu)化問題求解技術(shù)發(fā)展已經(jīng)相當(dāng)成熟,地球物理學(xué)者往往關(guān)注如何提高迭代中的正演計(jì)算效率,因此4.1節(jié)闡述的正演算法的發(fā)展,直接促進(jìn)了反演計(jì)算的效率和發(fā)展.

    圖4 2D 海洋 CSEM問題的非結(jié)構(gòu)化三角形網(wǎng)格的自適應(yīng)加密 (Li and Key, 2007)Fig.4 Adaptive refinements for 2D marine CSEM problems with triangles (Li and Key, 2007)

    圖5 RMT問題中三維梯形模型的非結(jié)構(gòu)化四面體網(wǎng)格的自適應(yīng)加密,為單元誤差(Ren et al., 2013)Fig.5 Adaptive refinements for 3D RMT problems with tetrahedral elements (Ren et al., 2013)

    圖6 基于區(qū)域分解的大型CSEM 3D 問題的MPI計(jì)算, 不同顏色代表不同的子區(qū)域,最上藍(lán)色為空氣區(qū)域,網(wǎng)格包含2.8 百萬個四面體,16CPU上計(jì)算時間為1.5 min (Puzyrev et al., 2013)Fig.6 Large-scale parallel computation with MPI. The grid contains 2.8 millions elements, with a run-time of 1.5 minutes using 16 CPUs (Puzyrev et al., 2013)

    圖7 CSEM各向異性3-D反演結(jié)果圖(Newman et al., 2010)Fig.7 Inversion results for CSEM anisotropic problems (Newman et al., 2010)

    圖8 2D RMT反演結(jié)果,頂圖忽略位移電流造成了錯誤的高阻區(qū)域.底圖未忽略位移電流,結(jié)果與測井結(jié)果吻合 (Kalscheuer et al., 2008)Fig.8 2D RMT inversion results. The top figure shows the case of ignoring the displacements which leads to an incorrect highly resistive area. The bottom figure shows that when considering the displacement currents, the inverted result agrees well with logging results (Kalscheuer et al., 2008)

    圖9 巴西米納斯吉拉斯州xxx測線二維反演結(jié)果圖(Terra and Menezes, 2012)Fig.9 2D inversion results from an AMT profile in the Minas Gerais State of Brazil (Terra and Menezes, 2012)

    反演多解性促進(jìn)了聯(lián)合反演的發(fā)展,因?yàn)槁?lián)合反演可減低結(jié)果的多解性,提高結(jié)果的精度和可靠性. Hu等(2009)發(fā)表了2.5D頻率域電磁反演和2D地震聯(lián)合反演;Candansayar 和Tezkan(2008)采用最小二乘法進(jìn)行了RMT和DC的2D聯(lián)合反演;Bala和Pieta(2010)采用蒙特卡羅法進(jìn)行了相似的1D DC/AMT聯(lián)合反演.另外,隨著計(jì)算資源的大幅度提高,各向異性的反演越來越成為一種研究熱點(diǎn),例如Newman等(2010)進(jìn)行了CSEM數(shù)據(jù)各向異性3D反演(見圖7).

    位移電流在RMT反演過程中的應(yīng)用越來越得到重視.在高阻環(huán)境和高頻率情況下,位移電流并不遠(yuǎn)遠(yuǎn)小于傳導(dǎo)電流.對于淺地表問題,特別是在地形起伏嚴(yán)重的區(qū)域,位移電流應(yīng)該考慮在反演過程中,忽略位移電流會造成反演結(jié)果錯誤.比如,Kalscheuer等(2008)詳細(xì)地分析了位移電流在RMT反演中的重要性,給出了忽略位移電流會造成反演結(jié)果的不合理性(見圖8).

    另外,超級計(jì)算機(jī)促使了反演的工業(yè)化實(shí)現(xiàn).高性能計(jì)算平臺的發(fā)展經(jīng)歷了單核到多核,再到近年來的GPU.Alumbaugh等(1996)、Alumbaugh和Newman(1997), Newman和Alumbaugh (1997, 1999)首先將高性能計(jì)算平臺應(yīng)用在大規(guī)模電磁計(jì)算中,完成了真實(shí)的3D地質(zhì)模型的電磁模擬.Commer和Newman(2008)、Newman等(2010)采用了32768個CPU進(jìn)行了3D CSEM數(shù)據(jù)成像的大規(guī)模并行計(jì)算,計(jì)算時間為幾天,奠定了大規(guī)模的3D電磁問題成像計(jì)算在工業(yè)中應(yīng)用的基礎(chǔ).Weiss和Schultz (2011)提供了使用單一GPU進(jìn)行大規(guī)模電磁數(shù)據(jù)模擬的可行性研究.

    5 實(shí)例應(yīng)用

    隨著數(shù)據(jù)采集技術(shù)、信號處理技術(shù)、正反演算法的快速發(fā)展,頻率域電磁法得到了人們越來越多的關(guān)注,其應(yīng)用范圍越來越寬廣.由于其信號頻率范圍的特性,上述三類頻率域電磁法主要應(yīng)用于礦產(chǎn)資源勘探、地?zé)豳Y源或地下水勘探、環(huán)境與工程地球物理等淺部地球物理問題.

    5.1 礦產(chǎn)資源

    隨著我國國民經(jīng)濟(jì)的快速發(fā)展,各種礦產(chǎn)資源的需求量越來越大.近年來,面對礦產(chǎn)資源的日益枯竭,開展中深部(大于500 m)找礦迫在眉睫(夏時斌等,2012).最新的成礦理論研究和深部定位預(yù)測驗(yàn)證結(jié)果均表明,地下500~1500 m深度見礦范例眾多,表明大陸深部蘊(yùn)藏著潛力巨大的礦產(chǎn)資源(葉益信等,2011),因此如何準(zhǔn)確、有效地開展中深部和各類隱伏區(qū)的探礦找礦,已成為現(xiàn)階段及未來的重大研究任務(wù).

    音頻大地電磁法和可控源音頻大地電磁法已逐步應(yīng)用在隱伏礦及深部礦的勘察中,并取得了不錯的效果(Sheard et al., 2005;Zhdanov, 2010;Vallée et al., 2011;Smith, 2014).在AMT和CSAMT兩種方法中,AMT為天然場源,抗干擾能力差,但其裝備簡便,易于操作,適合在地形復(fù)雜,干擾相對較小的地區(qū)開展勘察;而CSAMT采用人工源,信號強(qiáng),抗干擾能力大,但源的布設(shè)及采集較為復(fù)雜,適用于地形簡單但干擾相對較大的地區(qū).

    AMT和CSAMT受高阻屏蔽作用小,對低阻體反映靈敏.圖9給出了一個典型的AMT勘探低阻角礫云母橄欖巖的實(shí)例(La Terra and Menezes, 2012),角礫云母橄欖巖是原生金剛石的成礦標(biāo)志,而其電性特征表現(xiàn)為低阻.針對AMT死頻帶問題,La Terra等人采用發(fā)射天線,在距測點(diǎn)四倍趨膚深度外發(fā)射人工信號,以此提高頻帶信噪比及數(shù)據(jù)質(zhì)量.對AMT數(shù)據(jù)進(jìn)行了可靠的數(shù)據(jù)處理,并獲得了各個測線的二維反演結(jié)果.為進(jìn)一步驗(yàn)證二維反演結(jié)果的可靠性,利用WSINV3DMT程序(Siripunvaraporn and Sarakorn 2011),作者以二維反演結(jié)果為初始模型,進(jìn)行了3D反演(圖9b).三維反演結(jié)果與二維反演結(jié)果的一致性驗(yàn)證了反演結(jié)果的可靠.

    近年來,廣域電磁法以其探測深度大,分辨率高,抗干擾能力強(qiáng)、儀器輕便、觀測速度快,勘探范圍大等特點(diǎn),在深部資源勘探中發(fā)揮著越來越重要的作用. 胡艷芳和何繼善(2015)在安徽銅陵地區(qū)開展了強(qiáng)干擾礦區(qū)廣域電磁法的有效性試驗(yàn),為保證數(shù)據(jù)采集的質(zhì)量,進(jìn)行了強(qiáng)干擾區(qū)最佳發(fā)射電流的試驗(yàn)研究(電流從100~10 A),結(jié)果表明當(dāng)發(fā)射電流大于20 A時,數(shù)據(jù)信噪比才能得以保證. 圖10為銅陵強(qiáng)干擾礦區(qū)某測線2D反演結(jié)果,可以看出礦區(qū)地層的似層狀體現(xiàn)非常明顯,礦體橫向和深度范圍內(nèi)出現(xiàn)了低阻異常,并識別出了兩條斷層,與已知鉆孔信息一致.經(jīng)地質(zhì)、物探等綜合信息,推斷2號異常區(qū)還存在成礦的空間,該處低阻異常為礦體或礦化引起.

    5.2 地?zé)豳Y源或地下水勘探

    地?zé)豳Y源是集“熱、礦、水”三位一體的清潔可再生能源,目前已被廣泛應(yīng)用于取暖空調(diào)、生活熱水、醫(yī)療等許多領(lǐng)域,具有巨大的開發(fā)前景(周厚芳等,2004;黃力軍等,2004;廖忠禮等,2006).

    地?zé)豳Y源一般位于充滿流體的斷層或裂縫中,這些流體富含鹽分,具有良好的導(dǎo)電性,相對圍巖表現(xiàn)為低阻特性,這些特征使得電法及電磁法在地?zé)豳Y源的勘探中應(yīng)用效果顯著.一般地?zé)豳Y源主要分為巖漿型與斷陷盆地型兩種類型,巖漿型地?zé)豳Y源分布范圍小,溫度高,但是埋深淺,可采用常規(guī)直流電阻率法進(jìn)行勘探. 斷陷盆地型地?zé)豳Y源是目前利用最廣的地?zé)豳Y源,具有分布廣、中低溫、埋深大等特點(diǎn),可利用勘探深度較大的AMT和CSAMT以及MT進(jìn)行勘探. 對于中淺部地?zé)豳Y源,常選用AMT或CSAMT.

    Barcelona等(2013)給出了一個AMT勘探地?zé)豳Y源的實(shí)例,為了確定阿根廷西北部地區(qū)地下電阻率的分布,圈定出與溫泉相關(guān)的區(qū)域,綜合考慮測區(qū)的地形、水文、地質(zhì)及構(gòu)造等信息,選用AMT法成功進(jìn)行了地?zé)豳Y源的前期勘探.在該實(shí)例中,由于測區(qū)位于山區(qū),地形起伏大,區(qū)域維性相對比較復(fù)雜,所以根據(jù)區(qū)域與局部的地質(zhì)特征先確定出可能的地電響應(yīng),然后對AMT數(shù)據(jù)進(jìn)行了維性分析,分別做了一維和二維反演,綜合實(shí)際水文、地質(zhì)等先驗(yàn)信息,二維反演解釋結(jié)果更為合理可靠,成功獲得了地下一定深度內(nèi)的電阻率分布圖.

    Fu等(2012)給出了一個CSAMT勘探地下水資源的實(shí)例,在北京郊區(qū)強(qiáng)干擾條件下成功進(jìn)行了地下水的勘探工作.在該實(shí)例中,為保證強(qiáng)干擾條件下的數(shù)據(jù)質(zhì)量,通過加大發(fā)射功率、延長偶極子長度、增加測站采集時間等措施來提高數(shù)據(jù)信噪比;在數(shù)據(jù)處理時,利用HHT(湯井田等, 2009,蔡劍華等, 2011)進(jìn)行了信噪分離.圖11 給出了1D反演結(jié)果及地質(zhì)解釋示意圖,通過與已知的實(shí)際鉆孔資料對比,清晰地查明了高阻石灰石巖層的埋藏深度,并根據(jù)含水帶低阻特性,成功圈出了地下400~800 m處的含水破碎帶.

    總的來說,在地?zé)峄虻叵滤Y源的勘探中,一般常根據(jù)測區(qū)地形、干擾情況以及地?zé)崧癫厣疃鹊葘?shí)際因素選取合適的電磁方法,對采集的數(shù)據(jù)進(jìn)行可靠的數(shù)據(jù)處理,獲得有用的信號.然后,針對測區(qū)的地形,進(jìn)行地形校正,分析采集數(shù)據(jù)的維性,確定地下水資源的流動方向(走向),在二維假設(shè)成立的情況下,才能夠運(yùn)用二維反演算法.對于二維假設(shè)不滿足的情況(地下水流動一般呈不規(guī)則狀態(tài)),三維反演是必要的.

    5.3 環(huán)境與工程地球物理

    環(huán)境與工程地球物理,其任務(wù)是通過各種高精度的地球物理探測,了解近地表地下空間內(nèi)的物質(zhì)分布、運(yùn)動規(guī)律等,從而為各類工程和中淺層資源的設(shè)計(jì)、施工、檢修、維護(hù)、開發(fā)提供依據(jù).

    圖10 安徽銅陵強(qiáng)干擾礦區(qū)某線廣域電磁法二維反演結(jié)果P2d,二疊系上統(tǒng)大隆組;T,三疊系;C,石炭系;P2l,二疊系上統(tǒng)龍?zhí)督M;P1,二疊系下統(tǒng)地層.Fig.10 2D inversion result in Tongling P2d,Dalong formation of Upper Permian; T, Triassic; C, Carboniferous; P2l, Longtan formation of Upper Permian; P1, formation of Under Permian.

    圖11 北京郊區(qū)某測線CSAMT反演結(jié)果及地質(zhì)解釋示意圖(Fu et al., 2013)Fig.11 The CSAMT inversion results and its geological interpretation along the profile located on the outskirts of the Beijing city (Fu et al., 2013)

    RMT是近20年才逐步發(fā)展的一種淺地表地球物理勘探方法(Müller et al., 1994; Tezkan et al., 1996; Tezkan et al., 2000; Tezkan et al., 2005; Bastani, 2001; Bastani et al., 2009; Bastani et al., 2013; Kalscheuer et al., 2008; Ren et al., 2012).RMT 法測量頻段為10000 Hz~250 kHz,勘探深度為數(shù)米到數(shù)十米,目前多應(yīng)用于淺地表工程、資源勘探及環(huán)境監(jiān)測.在RMT應(yīng)用中,對于淺地表問題,特別是在起伏地形嚴(yán)重的區(qū)域,位移電流應(yīng)該考慮在反演過程中,不適當(dāng)?shù)睾雎晕灰齐娏?,會造成反演結(jié)果錯誤(Kalscheuer et al., 2008).

    近年來,隨著計(jì)算機(jī)的快速發(fā)展,RMT三維反演也已經(jīng)逐漸應(yīng)用于實(shí)際生產(chǎn)中,并在數(shù)據(jù)解釋方面發(fā)揮著重要作用.Newman等(2003)采用RMT法對德國一地下垃圾場是否發(fā)生滲漏污染等情況進(jìn)行了勘探,在該實(shí)例中,對RMT數(shù)據(jù)進(jìn)行了二維和三維反演結(jié)果的對比.結(jié)果表明,二維反演結(jié)果過高估計(jì)了垃圾埋藏深度及范圍,得出了錯誤結(jié)論,而三維反演結(jié)果能夠精確地圈定垃圾坑的邊界及深度,判斷垃圾坑是否發(fā)生滲漏污染(圖12).

    另一方面,為了壓制反演的多解性,進(jìn)一步提高RMT反演結(jié)果的可靠性, RMT與多種方法的聯(lián)合反演得到了發(fā)展和應(yīng)用.Candansayar和Tezkan(2008)從RMT和DC聯(lián)合反演的方法原理、合成理論數(shù)據(jù)以及野外實(shí)際數(shù)據(jù)的聯(lián)合反演等方面做了詳細(xì)的介紹,并成功地查明了德國凱爾彭實(shí)驗(yàn)基地?cái)鄬游恢眉s在地下16 m處.文中指出地表淺部的小尺度目標(biāo)體對RMT數(shù)據(jù)有較大的影響(易產(chǎn)生靜位移),但單獨(dú)的RMT反演結(jié)果不能分辨出淺地表的小尺度目標(biāo)體,并增加了反演的多解性.直流電阻率法對地表淺部目標(biāo)體有很好的反映,因此將RMT聯(lián)合DC進(jìn)行反演,不僅清晰查明了地表淺部小尺度目標(biāo)體以及深部的電性特征,同時極大地減小了反演結(jié)果的非唯一性,提高了反演結(jié)果模型分辨率,在此次實(shí)際勘查中取得了良好的應(yīng)用效果.

    圖12 測線反演結(jié)果對比圖(a)二維反演結(jié)果; (b)三維反演結(jié)果. 鉆孔資料顯示垃圾坑基底深度約為13 m (Newman et al., 2003)Fig.12 Comparison of the inversion results(a) 2D inversion result; (b) 3D inversion result. Landfill drilling data show the basement has a depth of 13 meters (Newman et al., 2003)

    6 結(jié)論

    基于目前最新的文獻(xiàn),本文從數(shù)據(jù)采集技術(shù)、信號處理技術(shù)、反演算法和實(shí)例分析等四個方面分別論述了頻率域電磁法的發(fā)展現(xiàn)狀,主要的結(jié)論可概述為:

    (1)數(shù)據(jù)采集與儀器.數(shù)據(jù)采集方案由單點(diǎn)張量采集向多點(diǎn)同步、多參考道、三維陣列的方向發(fā)展,主動源電磁法在信號發(fā)送方式、編碼方式、場源布設(shè)方法以及接收端的觀測方案等方面涌現(xiàn)出諸多新嘗試,提高了數(shù)據(jù)采集質(zhì)量和效率. 隨著儀器開發(fā)技術(shù)的日趨成熟,采集系統(tǒng)已逐步實(shí)現(xiàn)多方法集成、全三維采集,元器件精度的提高也使得數(shù)據(jù)采集頻率范圍、采集精度及效率等獲得了大幅擴(kuò)展和提高.尤其可喜的是,隨著人力與資金的大幅投入,國產(chǎn)儀器與國外儀器在質(zhì)量方面的差距正在逐步縮小.

    (2)數(shù)據(jù)處理技術(shù)方面.由于天然電磁場強(qiáng)度較低,AMT數(shù)據(jù)的信噪比通常不高.為了獲得或者分離出有用的信號,可采用一些新的信噪分離技術(shù),如數(shù)學(xué)形態(tài)濾波技術(shù)、Hilbert-Huang變換等.張量阻抗估計(jì)經(jīng)歷了從最小二乘估計(jì)到穩(wěn)健阻抗估計(jì),從引入單參考道到多道數(shù)據(jù)的利用等階段.局部畸變?nèi)匀皇秦酱鉀Q的難題,目前的主要校正方法為阻抗張量分解、空間濾波及聯(lián)合反演等.CSAMT和RMT方法由于人工源的引入,信噪比得到大幅提高,其數(shù)據(jù)處理技術(shù)的發(fā)展集中在非平面波效應(yīng)、陰影和場源附加效應(yīng)的識別與壓制等方面.值得提出的是,全區(qū)電阻率以及廣域電磁法的提出,有效地利用了非遠(yuǎn)區(qū)的數(shù)據(jù),為非平面波效應(yīng)問題的解決提供了途徑.

    (3)正反演算法方面.為了求解復(fù)雜AMT/RMT/CSAMT/WFEM問題,有限元法成為了不可或缺的正演計(jì)算技術(shù).為了進(jìn)一步提高有限單元法的計(jì)算精度和效率,基于邊的矢量有限元和可以離散任意復(fù)雜模型的非結(jié)構(gòu)化網(wǎng)格越來越成為主流.另外,由于直接求解器領(lǐng)域的在求解大型線性方程組的快速發(fā)展,對于多源地球物理問題,直接法有取代傳統(tǒng)的迭代法的趨勢.相對于正演算法的快速發(fā)展,由于其非線性,地球物理領(lǐng)域的反演算法本身的發(fā)展相對滯后,基本上為對經(jīng)典的非線性最優(yōu)化算法的應(yīng)用,如共軛梯度法、高斯牛頓算法等.由于其超大的計(jì)算量,完全非線性反演算法如遺傳算法等仍然處于1D、2D問題.與此不同,基于線性化的反演算法,如共軛梯度法、高斯牛頓算法等,3D模型的反演已經(jīng)變得相對實(shí)用.另外,為了提高反演算法的精度,聯(lián)合反演,如與地震數(shù)據(jù)等其他地球物理數(shù)據(jù)的聯(lián)合反演, 目前正獲得越來越多的應(yīng)用.對于一些簡單的問題,電阻率的各向異性的反演也獲得了一些重視.RMT問題的頻率較高,需要考慮位移電流,為了獲得正確的RMT反演結(jié)果,必須改正現(xiàn)有的代碼,使其包含位移電流項(xiàng).由于WFEM法出現(xiàn)較晚,需要開發(fā)更為可靠的正反演算法.

    (4)實(shí)際應(yīng)用.目前上述頻率域電磁法在不同的地球物理問題中得到了廣泛的應(yīng)用,如礦產(chǎn)資源勘探、地?zé)峄虻叵滤Y源勘探、環(huán)境與工程地球物理和油氣勘探等.由于地下結(jié)構(gòu)可能存在2D、3D特性,在應(yīng)用反演算法進(jìn)行數(shù)據(jù)解釋時,需要對地下結(jié)構(gòu)的維性進(jìn)行正確的判別.在2D特性不明顯,或者3D特性明顯時,需要采用3D反演算法.另外,即使對于2D特性明顯的區(qū)域,3D反演不僅可以驗(yàn)證結(jié)果,還可以輔助2D反演,獲得一些更為精細(xì)的局部異常結(jié)構(gòu).因此,三維采集與反演已逐漸成為主流.此外,為了獲得更為可信的反演結(jié)果,聯(lián)合反演的重要性日益明顯.

    Alumbaugh D L, Newman G A, Prevost L, et al. 1996. Three-dimensional, wideband electromagnetic modeling on massively parallel computers.RadioScience, 31(1): 1-23.

    Alumbaugh D L, Newman G A. 1997. Three-dimensional massively parallel electromagnetic inversion-part II. Analysis of a crosswell electromagnetic experiment.GeophysicalJournalInternational, 128(2): 355-363.

    Amestoy P R, Duff I S, L′Excellent J Y, et al. 2002. A fully asynchronous multi-frontal solver using distributed dynamic scheduling.SIAMJournalonMatrixAnalysisandApplications, 23(1): 15-41.Aprea C, Booker J R, Torquil Smith J. 1997. The forward problem of electromagnetic induction: accurate finite-difference approximations for two dimensional discrete boundaries with arbitrary geometry.GeophysicalJournalInternational, 129(1): 29-40.Avdeev D B, Kuvshinov A V, Pankratov O V, et al. 2002. Three dimensional induction logging problems, Part I: an integral equation solution and model comparisons.Geophysics, 67(2): 413-426.

    Avdeev D B. 2005. Three-dimensional electromagnetic modelling and inversion from theory to application.SurveysinGeophysics, 26(6): 767-799. Badea E A, Everett M E, Newman G A, et al. 2001. Finite-element analysis of controlled-source electromagnetic induction using Coulomb-gauged potentials.Geophysics, 66(3): 786-799.

    Bahr K. 1988. Interpretation of the magnetotelluric impedance tensor: regional induction and local telluric distortion.J.Geophys., 62: 119-127.

    Bahr K. 1991. Geological noise in magnetotelluric data: a classification of distortion types.PhysicsoftheEarthandPlanetaryInteriors, 66(1-2): 24-38.

    Bastani M. 2001. EnviroMT: A new controlled source/radio magnetotelluric system [Ph. D. thesis]. Uppsala: Acta Universitatis Upsaliensis. Bastani M, Malehmir A, Ismail N, et al. 2009. Delineating hydrothermal stockwork copper deposits using controlled-source and radio-magnetotelluric methods: A case study from northeast Iran.Geophysics, 74(5): B167-B181. Bastani M, Persson, L, Beiki M, et al. 2013. A radio magnetotelluric study to evaluate the extents of a limestonequarry in Estonia.GeophysicalProspecting, 61(3): 678-687.Beamish D, Travassos J M. 1992. The use of the D+solution in magnetotelluric interpretation.JournalofAppliedGeophysics, 29(1): 1-19. Berdichevsky M N, Dimitriev V T. 1976. Basic principles of interpretation of magnetotelluric sounding curves. ∥Ada S A ed. Geoelectric and Geothermal Studies, (East-Central Europe, Soviet-Asia), KAPG Geophysical Monograph, Akadémiai Kiadó, Budapest. 165-221. Bielecka M, Danek T, Wojdyla M, et al. 2009. Neural networks application to reduction of train caused distortions in magnetotelluric measurement data.SchedaeInformaticae, 17-18: 75-86.

    B?erner R U. 2010. Numerical modelling in geo-electromagnetics: advances and challenges.SurveysinGeophysics, 31(2): 225-245.

    Booker J R. 2014. The magnetotelluric phase tensor: a critical review.SurveysinGeophysics, 35(1): 7-40.

    Boschetto N B, Hohmann G W. 1991. Controlled-source audio frequency magnetotelluric responses of three-dimensional bodies.Geophysics, 56(2): 255-264.

    Cai J H, Tang J T, Hua X R, et al. 2009. An analysis method for magnetotelluric data based on the Hilbert-Huang transform.ExplorationGeophysics, 40(2): 197-205.

    Cai J H, Tang J T, Wang X C. 2011. Human noise elimination for magnetotelluric data based on empirical mode decomposition.JournalofCentralSouthUniversity(ScienceandTechnology) (in Chinese), 42(6): 1786-1790. Cai J H. 2013. Magnetotelluric response function estimation based on Hilbert-Huang transform.PureandAppliedGeophysics, 170(11): 1899-1911.

    Cai J H. 2014. A combinatorial filtering method for magnetotelluric time-series based on Hilbert-Huang transform.ExplorationGeophysics, 45(2): 63-73.

    Caldwell T G, Bibby H M, Brown C. 2004. The magnetotelluric phase tensor.GeophysicalJournalInternational, 158(2): 457-469.

    Campanyà J, Ledo J, Queralt P, et al. 2014. A new methodology to estimate Magnetotelluric (MT) Tensor relationships: Estimation of Local transfer-functions by Combining Interstation Transfer-functions (ELICIT).GeophysicalJournalInternational, 198(1): 484-494.Emin Candansayar M, Tezkan B. 2008. Two-dimensional joint inversion of radiomagnetotelluric and direct current resistivity data.GeophysicalProspecting, 56(5): 737-749.

    Cao C Q. 1978. The apparent resistivity for layered earth.ChineseJ.Geophys. (ActaGeophysicaSinica) (in Chinese), 21(3): 248-261.

    Chant I J, Hastie L M. 1992. Time-frequency analysis of magnetotelluric data.GeophysicalJournalInternational, 111(2): 399-413.

    Chave A D, Thomson D J, Ander M E. 1987. On the robust estimation of power spectra, coherences, and transfer functions.JournalofGeophysicalResearch:SolidEarth(1978—2012), 92(B1): 633-648. Chave A D, Thomson D J. 1989. Some Comments on magnetotelluric response function estimation.JournalofGeophysicalResearch:SolidEarth(1978—2012), 94(B10): 14215-14225.

    Chave A D, Smith J T. 1994. On electric and magnetic galvanic distortion tensor decompositions.JournalofGeophysicalResearch, 99(B3): 4669-4682.

    Chave A D, Thomson D J. 2003. A bounded influence regression estimator based on the statistics of the hat matrix.JournaloftheRoyalStatisticalSociety:SeriesC(AppliedStatistics), 52(3): 307-322.

    Chave A D, Luther D S, Meinen C S. 2004. Correction of motional electric field measurements for galvanic distortion.JournalofAtmosphericandOceanicTechnology, 21(2): 317-330.

    Chen H Y, Wei W B, Jing J E, et al. 2012. The generalized S-transform and its application in the processing of magnetotelluric sounding data.ProgressinGeophysics(in Chinese), 27(3): 988-996, doi: 10.6038/j.issn.1004-2903.

    Chen J, Jegen-Kulcsar M, Moorkamp M. 2010. The empirical mode decomposition and hilbert instantaneous spectrum on MT data Processing. ∥20th International Electromagnetic Induction Workshop, Giza, Egypt.

    Chen M S, Yan S. 2005. Analytical study on field zones, record rules, shadow and source overprint effects in CSAMT exploration.ChineseJ.Geophys. (in Chinese), 48(4): 951-958.

    Chen W Y, Xue G Q. 2014. Effective skin depth of whole EM field due to a grounded wire source.ChineseJ.Geophys. (in Chinese), 57(7): 2314-2320, doi: 10.6038/cjg20140725.

    Chen X B, Zhao G Z. 2009. Study on the transmitting mechanism of CSELF waves: response of the alternating current point source in the uniform space.ChineseJ.Geophys. (in Chinese), 52(8): 2158-2164, doi: 10.3969/j.issn.0001-5733.2009.08.026.

    Chen X B, Cai J T, Wang L F, et al. 2014. Refined techniques for magnetotelluric data processing and two-dimensional inversion (IV): Statistical image method based on multi-site, multi frequency tensor decomposition.ChineseJ.Geophys. (in Chinese), 57(6): 1946-1957, doi: 10.6038/cjg20140625.

    Cheng D F, Wang J, Li X P, et al. 2004. Development in the research of HSSMT instrument.ProgressinGeophysics(in Chinese), 19(4): 778-781, doi: 10.3969/j.issn.1004-2903.2004.04.012.

    Cheng H, Li D Q. 2014. A new electromagnetic field division method for horizontal electric dipole.JournalofCentralSouthUniversity(ScienceandTechnology) (in Chinese), 45(7): 2250-2258.

    Christopherson K R. 2012. Magnetotellurics (MT): Technique, Interpretation, and Application. http:∥www.searchanddiscovery.com/documents/geophysical/christopherson/index.htm.Coggon J H. 1971. Electromagnetic and electrical modeling by the finite element method.Geophysics, 36(1): 132-155.

    Commer M, Newman G A. 2008. New advances in three-dimensional controlled-source electromagnetic inversion.GeophysicalJournalInternational, 172(2): 513-535. Cooley J W, Tukey J W. 1965. An algorithm for the machine calculation of complex Fourier series.MathematicsofComputation, 19(90): 297-301.

    Cooley J W, Lewis P A, Welch P D. 1969. The fast Fourier transform and its applications.IEEETransactionsonEducation, 12(1): 27-34. Cui J L, Deng M, Jing J E, et al. 2013. Using independent component analysis to process magnetotelluric data.AppliedMechanicsandMaterials, 295-298: 2795-2798.

    Das U C. 1995. Apparent resistivity curves in Controlled-source electromagnetic sounding directly reflecting true resistivities in a layered earth.Geophysics, 60(1): 53-60.

    Di Q Y, Wang M Y, Wang R, et al. 2008. Study of the long dipole and large power electromagnetic field.ChineseJ.Geophys. (in Chinese), 51(6): 1917-1928.

    Di Q Y, Fang G Y, Zhang Y M. 2013. Research of the Surface Electromagnetic Prospecting (SEP) system.ChineseJ.Geophys. (in Chinese), 56(11): 3629-3639, doi: 10.6038/cjg20131104.

    Doherty J. 1988. EM modelling using surface integral equations.GeophysicalProspecting, 36(6): 644-668.

    Egbert G D, Booker J R. 1986. Robust estimation of geomagnetic transfer functions.GeophysicalJournalInternational, 87(1): 173-194.

    Egbert G D, Booker J R. 1989. Multivariate analysis of geomagnetic array data: 1. The response space.JournalofGeophysicalResearch:SolidEarth(1978—2012), 94(B10): 14227-14247. Egbert G D, Livelybrooks D W. 1996. Single Station magnetotelluric impedance estimation: coherence weighting and the regression M-estimate.Geophysics, 61(4): 964-970.

    Egbert G D. 1997. Robust Multiple-station magnetotelluric data processing.GeophysicalJournalInternational, 130(2): 475-496.

    Egbert G D. 2002. Processing and interpretation of electromagnetic induction array data.SurveysinGeophysics, 23(2-3): 207-249.

    Escalas M, Queralt P, Ledo J, et al. 2013. Polarization analysis of magnetotelluric time series using a wavelet-based scheme: a method for detection and characterization of cultural noise sources.PhysicsoftheEarthandPlanetaryInteriors, 218: 31-50.

    Everett M E. 2011. Theoretical developments in electromagnetic induction geophysics with selected applications in the near surface.SurveysinGeophysics, 33(1): 29-63.

    Farquharson C G, Miensopust M P. 2011. Three-dimensional finite-element modelling of magnetotelluric data with a divergence correction.JournalofAppliedGeophysics, 75(4): 699-710.

    Feng B, Wang J L, Zhou X W, et al. 2013. Application of full-region apparent resistivity of CSAMT Ex in exploration.CoalGeology&Exploration(in Chinese), 41(6): 78-82.

    Fischer G, Weibel P. 1991. A new look at an old problem: Magnetotelluric modelling of 1-D structures.GeophysicalJournalInternational, 106(1): 161-167.

    Fontes S L, Harinarayana T, Dawes G, et al. 1988. Processing of noisy magnetotelluric data using digital filters and additional data selection criteria.PhysicsoftheEarthandPlanetaryInteriors, 52(1): 30-40.

    Franke A, Borner R U, Spitzer K. 2008. Is there a most efficient formulation of the three-dimensional magnetotelluric boundary value problem? ∥19th Workshop on Electromagnetic Induction in the Earth, Beijing, China.

    Fu C M, Di Q Y, Xu C, et al. 2012. Electromagnetic fields for different type sources with effect of the ionosphere.ChineseJ.Geophys. (in Chinese), 55(12): 3958-3968, doi: 10.6038/j.issn.0001-5733.2012.12.008.

    Fu C M, Di Q Y, An Z G. 2013. Application of the CSAMT method to groundwater exploration in a metropolitan environment.Geophysics, 78(5): B201-B209.

    Gamble T, Goubau W M, Clarke J. 1979. Magnetotellurics with a remote magnetic reference.Geophysics, 44(1): 53-68.

    Gao W. 1991. Line source effects on magnetotelluric responses.ChineseJ.Geophys. (in Chinese), 34(2): 210-215.

    García X, Jones A G. 2005. A new methodology for the acquisition and processing of Audio-Magnetotelluric (AMT) data in the AMT dead band.Geophysics, 70(5): G119-G126.

    Grayver A V, Bürg M. 2014. Robust and scalable 3D geo-electromagnetic modelling approach using the finite element method.GeophysicalJournalInternational, 198(1): 110-125.

    Groom R W, Bailey R C. 1989. Decomposition of magnetotelluric impedance tensors in the presence of local Three-Dimensional galvanic distortion.JournalofGeophysicalResearch:SolidEarth(1978—2012), 94(B2): 1913-1925.

    Groom R W, Bailey R C. 1991. Analytic investigations of the effects of Near-Surface Three-dimensional galvanic scatters on MT tensor decompositions.Geophysics, 56(4): 496-518.

    Haber E, Ascher U M, Aruliah D A, et al. 2000. Fast simulation of 3D electromagnetic problems using potentials.JournalofComputationalPhysics, 163(1): 150-171.

    Hattingh M. 1989. The use of Data-adaptive filtering for noise removal on magnetotelluric data.PhysicsoftheEarthandPlanetaryInteriors, 53(3-4): 239-254.

    He J S. 1991. Controlled Source Audio-frequency Magnetotellurics (in Chinese). Changsha: Central South University Press.

    He J S. 2010a. Closed addition in a three-element set and 2nsequence pseudo-random signal coding.JournalofCentralSouthUniversity(ScienceandTechnology) (in Chinese), 41(2): 632-627.

    He J S. 2010b. Wide field electromagnetic sounding methods.JournalofCentralSouthUniversity(ScienceandTechnology) (in Chinese), 41(3): 1065-1072. He J S. 2010c. Wide Field Electromagnetic Method and Pseudo Random Signal Method (in Chinese). Beijing: Higher Education Press.

    He L F, Wang X B, Wang C X. 1999. Improving the s/n ratio of MT data by wavelet analysis.JournalofChengduUniversityofTechnology(in Chinese), 26(3): 299-302.

    Hohmann G W. 1975. Three-dimensional induced polarization and electromagnetic modeling.Geophysics, 40(2): 309-324.

    Hou J S, Mallan R K, Torres-Verdín C. 2006. Finite-difference simulation of borehole EM measurements in 3D anisotropic media using coupled scalar-vector potentials.Geophysics, 71(5): G225-G233.

    Hu W Y, Abubakar A, Habashy T M. 2009. Joint electromagnetic and seismic inversion using structural constraints.Geophysics, 74(6): R99-R109.

    Hu Y F, He J S. 2015. Comparison of wide field electromagnetic method with CSAMT method in strong interferential mining area [Ph. D. thesis] (in Chinese). Changsha: Central South of University.

    Huang H P, Piao H R. 1992. Full-wave apparent resistivity from vertical magnetic dipole frequency sounding on a layered earth.ChineseJ.Geophys. (in Chinese), 35(3): 389-395.

    Huang L J, Lu G F, Liu R D, et al. 2004. The application of electromagnetic sounding method in deep geothermal investigation.Geophysical&GeochemicalExploration(in Chinese), 28(6): 493-495. Jiang L, Xu Y X. 2013. Multi-station superposition for magnetotelluric signal.StudiaGeophysicaEtGeodaetica, 57(2): 276-291.

    Jing J E, Wei W B, Chen H Y, et al. 2012. Magnetotelluric sounding data processing based on generalized S transformation.ChineseJ.Geophys. (in Chinese), 55(12): 4015-4022, doi: 10.6038/j.issn.00015733. Jones A G, Chave A D, Egbert G D, et al. 1989. A comparison of techniques for magnetotelluric response function estimation.JournalofGeophysicalResearch:SolidEarth(1978—2012), 94(B10): 14201-14213.

    Jones A G, Groom R W. 1993. Strike-angle determination from the magnetotelluric impedance tensor in the presence of noise and local distortion: Rotate at your peril!GeophysicalJournalInternational, 113(2): 524-534.

    Kalscheuer T, Pedersen L B, Siripunvaraporn W. 2008. Radiomagnetotelluric two-dimensional forward and inverse modelling accounting for displacement currents.GeophysicalJournalInternational, 175(2): 486-514.

    Kappler K N. 2012. A data variance technique for automated despiking of magnetotelluric data with a remote reference.GeophysicalProspecting, 60(1): 179-191.

    Key K, Weiss C. 2006. Adaptive finite-element modeling using unstructured grids: the 2D magnetotelluric example.Geophysics, 71(6): G291-G299.La Terra E F, Menezes P T L. 2012. Audiomagnetotelluric 3D imaging of the Regis kimberlite pipe, Minas Gerais, Brazil.JournalofAppliedGeophysics, 77: 30-38.

    Larsen J C, Mackie R L, Manzella A, et al. 1996. Robust smooth magnetotelluric transfer functions.GeophysicalJournalInternational, 124(3): 801-819. Lei D, Zhang G H, Huang G Y, et al. 2014. The application of tensor CSAMT method.ChineseJournalofEngineeringGeophysics(in Chinese), 11(3): 286-294.

    Li Y G, Key K. 2007. 2D marine controlled-source electromagnetic modeling: Part 1—an adaptive finite-element algorithm.Geophysics, 72(2): WA51-WA62.

    Liao Z L, Zhang Y J, Chen W B, et al. 2006. Available persist exploitation and utilization of geothermal resources.ChinaMiningMagazine(in Chinese), 15(10): 8-11.

    Lin P R, Zheng C J, Shi F S, et al. 2006. The research of integrated electromagnetic method system.ActaGeologicaSinica(in Chinese), 80(10): 1539-1548.

    Lin P R, Guo P, Shi F S, et al. 2010. A study of the techniques for large-depth and multi-functional electromagnetic survey.ActaGeoscienticaSinica(in Chinese), 31(2): 149-154.

    Lin W. 2009. Transition region electromagnetic field characteristics of the CSAMT method.Geophysical&GeochemicalExploration(in Chinese), 33(2): 148-150. Liu E H, Lamontagne Y. 1998. Geophysical application of a new surface integral equation method for EM modeling.Geophysics, 63(2): 411-423.

    Liu G M, Becker A. 1992. Evaluation of terrain effects in AEM surveys using the boundary element method.Geophysics, 57(2): 272-278.

    Liu H, He L F, Wang X B, et al. 2004. Adaptability of wavelet analysis in MT noise-suppressed processing.OilGeophysicalProspecting(in Chinese), 39(3): 338-341.

    Luo W B, Li Q C, Tang J T. 2012. Coded source electromagnetic sounding method.ChineseJ.Geophys. (in Chinese), 55(1): 341-349, doi: 10.6038/j.issn.0001-5733.2012.01.035.

    Mackie R L, Smith J T, Madden T R. 1994. Three-dimensional electromagnetic modeling using finite difference equations: The magnetotelluric example.RadioScience, 29(4): 923-935.

    Manoj C, Nagarajan N. 2003. The application of artificial neural networks to magnetotelluric time-series analysis.GeophysicalJournalInternational, 153(2): 409-423.

    Mcneice G W, Jones A G. 2001. Multisite, multifrequency tensor decomposition of magnetotelluric data.Geophysics, 66(1): 158-173.

    Meng Q K. 2013. Preliminary study of tensor CSAMT data acquisition and processing technology [Master′s thesis] (in Chinese). Beijing: Chinese Academy of Geological Sciences.

    Mitsuhata Y, Uchida T. 2004. 3D magnetotelluric modeling using the T-Ω finite-element method.Geophysics, 69(1): 108-119.

    Nam M J, Kim H J, Song Y, et al. 2007. 3D magnetotelluric modelling including surface topography.GeophysicalProspecting, 55(2): 277-287. Nédélec J. 1986. A new family of mixed finite elements in R3.NumerischeMathematik, 50(1): 57-81.

    Neukirch M, Garcia X. 2014. Nonstationary magnetotelluric data processing with instantaneous parameter.JournalofGeophysicalResearch:SolidEarth, 119(3): 1634-1654.Newman G A, Alumbaugh D L. 1997. Three-dimensional massively parallel electromagnetic inversion—I. theory.GeophysicalJournalInternational, 128(2): 345-354.

    Newman G A, Alumbaugh D L. 1999. 3-D electromagnetic modeling and inversion on massively parallel computers. ∥ Oristaglio M N, Spies B R eds. Three-Dimensional Electromagnetics. Society of exploration geophysicists, Geophysical Developments No. 7, Tulsa OK, 299-321. Newman G A, Alumbaugh D. 2002. Three-dimensional induction logging problems, part 2: a finite-difference solution.Geophysics, 67(2): 484-491. Newman G A, Recher S, Tezkan B, et al. 2003. 3D inversion of a scalar radio magnetotelluric field data set.Geophysics, 68(3): 791-802.

    Newman G A, Commer M, Carazzone J J. 2010. Imaging CSEM data in the presence of electrical anisotropy.Geophysics, 75(2): F51-F61.

    Oettinger G, Haak V, Larsen J C. 2001. Noise reduction in magnetotelluric Time-series with a new signal-noise separation method and its application to a field experiment in the Saxonian granulite massif.GeophysicalJournalInternational, 146(3): 659-669.

    Parker R L. 1980. The inverse problem of electromagnetic induction: existence and construction of solutions based on incomplete data.JournalofGeophysicalResearch:SolidEarth(1978—2012), 85(B8): 4421-4428.

    Parker R L, Whaler K A. 1981. Numerical methods for establishing solutions to the inverse problem of electromagnetic induction.JournalofGeophysicalResearch:SolidEarth(1978—2012), 86(B10): 9574-9584.

    Parker R L, Booker J R. 1996. Optimal One-dimensional inversion and bounding of magnetotelluric apparent resistivity and phase measurements.PhysicsoftheEarthandPlanetaryInteriors, 98(3-4): 269-282.

    Parry J R, Ward S H. 1971. Electromagnetic scattering from cylinders of arbitrary cross-section in a conductive half-space.Geophysics, 36(1): 67-100.

    Penz S, Chauris H, Donno D, et al. 2013. Resistivity modelling with topography.GeophysicalJournalInternational, 194(3): 1486-1497.

    Prácser E, Szarka L. 1999. A correction to Bahr′s “Phase Deviation” method for tensor decomposition.Earth,PlanetsandSpace, 51(10): 1019-1022.

    Puzyrev V J, Koldan J, de la Puente J, et al. 2013. A parallel finite-element method for three-dimensional controlled-source electromagnetic forward modelling.GeophysicalJournalInternational, 193(2): 678-693.Raiche A P. 1974. An integral equation approach to three-dimensional modelling.GeophysicalJournaloftheRoyalAstronomicalSociety, 36(2): 363-376.

    Ren. 2012. New Developments in Numerical Modeling of Broadband Geo-electromagnetic Fields. Zurich: ETH Zurich Press.

    Ren Z Y, Tang J T. 2010. 3D direct current resistivity modeling with unstructured mesh by adaptive finite-element method.Geophysics, 75(1): H7-H17

    Ren Z Y, Kalscheuer T, Greenhalgh S, et al. 2012. Boundary element solutions for broadband 3D geo-electromagnetic problems accelerated by multi-level fast multipole method.GeophysicalJournalInternational, 192(2): 473-499.

    Ren Z Y, Kalscheuer T, Greenhalgh S, et al. 2013. A goal-oriented adaptive finite-element approach for plane wave 3D electromagnetic modelling.GeophysicalJournalInternational, doi: 10.1093/gji/ggt154. Ren Z Y, Kalscheuer T, Greenhalgh S, et al. 2014. A finite-element based domain-decomposition approach for plane wave 3D electromagnetic modeling.Geophysics, 79(6): E255-E268.

    Ritter O, Junge A, Dawes G. 1998. New equipment and processing for magnetotelluric remote reference observations.GeophysicalJournalInternational, 132(3): 535-548.

    Schenk O, G?rtner K. 2004. Solving unsymmetric sparse systems of linear equations with PARDISO.FutureGenerationComputerSystems, 20(3): 475-487.

    Schwarzbach C, B?rner R U, Spitzer K. 2011. Three-dimensional adaptive higher order finite element simulation for geo-electromagnetics—a marine CSEM example.GeophysicalJournalInternational, 187(1): 63-74. Sheard S N, Ritchie T J, Christopherson K R, et al. 2005. Mining, environmental, petroleum, and engineering industry applications of electromagnetic techniques in geophysics.SurveysinGeophysics, 26(5): 653-669.Sims W E, Bostick Jr F X, Smith H W. 1971. The estimation of magnetotelluric impedance tensor elements from measured data.Geophysics, 36(5): 938-942.

    Siripunvaraporn W, Egbert G, Lenbury Y. 2002. Numerical accuracy of magnetotelluric modeling: a comparison of finite difference approximations.Earth,PlanetsandSpace, 54(6): 721-725.

    Siripunvaraporn W, Sarakorn W. 2011. An efficient data space conjugate gradient Occam′s method for three-dimensional magnetotelluric inversion.GeophysicalJournalInternational, 186(2): 567-579.

    Smirnov M Y, Egbert G D. 2012. Robust principal component analysis of electromagnetic arrays with missing data.GeophysicalJournalInternational, 190(3): 1423-1438. Smirnov M Y. 2003. Magnetotelluric data processing with a robust statistical procedure having a high breakdown point.GeophysicalJournalInternational, 152(1): 1-7. Smith J T. 1996. Conservative modeling of 3-D electromagnetic fields, Part II: bi-conjugate gradient solution and an accelerator.Geophysics, 61(5): 1319-1324.Smith R. 2014. Electromagnetic induction methods in mining geophysics from 2008 to 2012.SurveysinGeophysics, 35(1): 123-156.

    Sokolova E Y, Varentsov M I, Group E P W. 2005. The RRMC Technique fights highly coherent EM noise. ∥ Dtsch. Geophys. Ges. Potsdam, Germany, 124-136.

    Song S G, Tang J T, He J S. 1995. Wavelets analysis and the recognition, separation and removal of the static shift in electromagnetic soundings.ChineseJ.Geophys. (in Chinese), 38(1): 120-128.

    Spratt J E, Jones A G, Nelson K D, et al. 2005. Crustal structure of the India-Asia collision zone, Southern Tibet, from INDEPTH MT investigations.PhysicsoftheEarthandPlanetaryInteriors, 150(1-3): 227-237. Sternberg B K, Washburne J C, Pellerin L. 1988. Correction for the static shift in magnetotellurics using transient electromagnetic soundings.Geophysics, 53(11): 1459-1468.

    Streich R. 2009. 3D finite-difference frequency-domain modeling of controlled-source electromagnetic data: direct solution and optimization for high accuracy.Geophysics, 74(5): F95-F105.

    Sutarno D. 2008. Development of robust magnetotelluric impedance estimation: a review.IndonesianJournalofPhysics, 16(3): 79-89.

    Swift Jr C M. 1967. A magnetotelluric investigation an electrical conductivity anomaly in the southwestern United States [Ph. D. thesis]. Massachusetts: Massachusetts Institute of Technology.Tang J T, He J S. 1993. Effective resistivity defined by the electromagnetic fields induced by horizontal electric dipole.JournalofCentralSouthInstituteofMiningandMetallurgy(in Chinese), 24(2): 137-142.

    Tang J T, He J S. 1994a. A new method to define the full-zone resistivity in horizontal electric dipole frequency soundings on a layered earth.ChineseJ.Geophys. (in Chinese), 37(4): 543-552.

    Tang J T, He J S. 1994b. Impedance real part equivalent resistivity in frequency electromagnetic sounding of horizontal electric double source on horizontal multilayer media.Geophysical&GeochemicalExploration(in Chinese), 18(2): 92-96.

    Tang J T, He J S. 2005. Methods and Applications of CSAMT (in Chinese). Changsha: Central South University Press.

    Tang J T, Hua X R, Cao Z M, et al. 2008. Hilbert-Huang transformation and noise suppression of magnetotelluric sounding data.ChineseJ.Geophys. (in Chinese), 51(2): 603-610.

    Tang J T, Cai J H, Ren Z Y, et al. 2009. Hilbert-Huang transform and time frequency analysis of magnetotelluric signal.JournalofCentralSouthUniversity(ScienceandTechnology) (in Chinese), 40(5): 1399-1405

    Tang J T, Zhou C, Zhang L C. 2011. A new apparent resistivity of CSAMT defined by electric fieldy-direction.JournalofJilinUniversity(EarthScienceEdition) (in Chinese), 41(2): 552-558.

    Tang J T, Li J, Xiao X, et al. 2012a. Mathematical morphology filtering and noise suppression of magnetotelluric sounding data.ChineseJ.Geophys. (in Chinese), 55(5): 1784-1793, doi: 10.6038/j.issn.0001-5733.2012.05.036.

    Tang J T, Xu Z M, Xiao X, et al. 2012b. Effect rules of strong noise on magnetotelluric (MT) sounding in the Luzong ore cluster area.ChineseJ.Geophys. (in Chinese), 55(12): 4147-4159, 10.6038/j.issn.0001-5733.2012.12.027.

    Tang J T, Ge W N. 2012. Shadow and source overprint effects in 3D CSAMT.ComputingTechniquesforGeophysical&GeochemicalExploration(in Chinese), 34(1): 19-26.

    Tang J T, Zhou C, Wang X Y, et al. 2013. Deep electrical structure and geological significance of Tongling ore district.Tectonophysics, 606: 78-96.

    Tang J T, Zhang C, Xiao X, et al. 2013a. Comparison of methods for magnetotelluric impedance estimation.TheChineseJournalofNonferrousMetals(in Chinese), 23(9): 2351-2358.

    Tang J T, Zhou C, Xiao X. 2013b. Selection of minimum transmit-receive distance of CSAMT on complicated media.TheChineseJournalofNonferrousMetals(in Chinese), 23(6): 1681-1693.

    Tezkan B, H?rdt A, Gobashy M. 2000. Two-dimensional radiomagnetotelluric investigation of industrial and domestic waste sites in Germany.JournalofAppliedGeophysics, 44(2-3): 237-256.

    Tezkan B, Georgescu P, Fauzi U. 2005. A radiomagnetotelluric survey on an oil-contaminated area near the Brazi Refinery, Romania.GeophysicalProspecting, 53(3): 311-323.

    Tezkan B, Goldman M, Greinwald S, et al. 1996. A joint application of radiomagnetotellurics and transient electromagnetics to the investigation of a waste deposit in Cologne (Germany).JournalofAppliedGeophysics, 34(3): 199-212. Torres-Verdín C, Bostick Jr F X. 1992. Principles of spatial surface electric field filtering in magnetotellurics: Electromagnetic Array Profiling (EMAP).Geophysics, 57(4): 603-622. Trad D O, Travassos J M. 2000. Wavelet filtering of magnetotelluric data.Geophysics, 65(2): 482-491. Turberg P, Müller I, Flury F. 1994. Hydrogeological investigation of porous environments by radio magnetotelluric-resistivity (RMT-R 12-240 kHz).JournalofAppliedGeophysics, 31(1-4): 133-143.

    Tzanis A, Beamish D. 1989. A high-resolution spectral study of audiomagnetotelluric data and noise interactions.GeophysicalJournalInternational, 97(3): 557-572.

    Vallée M A, Smith R S, Keating P. 2011. Metalliferous mining geophysics-State of the art after a decade in the new millennium.Geophysics, 76(4): W31-W50.

    Varentsov I M, Sokolova E Y, Martanus E R, et al. 2003. System of electromagnetic field transfer operators for the BEAR array of simultaneous soundings: methods and results.Izvestiya-PhysicsoftheSolidEarth, 39(2): 118-148.

    Varentsov I M. 2006. Arrays of simultaneous electromagnetic soundings: design, data processing and analysis. ∥ Electromagnetic Sounding of the Earth′s Interior, Vol. 40. Methods in Geochemistry and Geophysics. Amsterdam: Elsevier, 259-273. Wang H, Wei W B, Jin S, et al. 2014. Removal of magnetotelluric noise based on synchronous time series relationship.ChineseJ.Geophys. (in Chinese), 57(2): 531-545, doi: 10.6038/cjg20140218.

    Wang R, Wang M Y, Lu Y L. 2004. CSAMT observation system study in high mountain and steep gorge area.ProgressinGeophysics(in Chinese), 19(1): 125-130, doi: 10.3969/j.issn.1004-2903.2004.01.018.

    Wang R, Di Q Y, Wang M Y, et al. 2009. Research on the effect of 3D body between transmitter and receivers on CSAMT response using Integral Equation method.ChineseJ.Geophys. (in Chinese), 52(6): 1573-1582, doi: 10.3969/j.issn.0001-5733.2009.06.019.

    Wang S M, Wang J Y. 2004. Application of higher-order statistics in magnetotelluric data processing.ChineseJ.Geophys. (in Chinese), 47(5): 928-934.

    Wang X X, Di Q Y, Xu C. 2014. Characteristics of multiple sources and tensor measurement in CSAMT.ChineseJ.Geophys. (in Chinese), 57(2): 651-661, doi: 10.6038/cjg20140228.

    Weaver J T, Agarwal A K, Lilley F E M. 2006. The relationship between the magnetotelluric tensor invariants and the phase tensor of Caldwell, Bibby, and Brown.ExplorationGeophysics, 37(3): 261-267.Weckmann U, Magunia A, Ritter O. 2005. Effective noise separation for magnetotelluric single site data processing using a frequency domain selection scheme.GeophysicalJournalInternational, 161(3): 635-652.

    Weidelt P. 1972. The inverse problem of geomagnetic induction.J.Geophys., 38: 257-289.

    Weidelt P. 1975. Electromagnetic induction in three-dimensional structures.J.Geophys., 41: 85-109.

    Weiss C J, Schultz A. 2011. An evaluation of parallelization strategies for low-frequency electromagnetic induction simulators using staggered grid discretizations. ∥ American Geophysical Union Fall Meeting Conference Proceedings, Informatics Session, San Francisco.

    Xia S B, Min G, Qiu L, et al. 2012. A study of MT method on the prospecting for deep ore deposits.ChineseJournalofEngineeringGeophysics(in Chinese), 9(5): 531-536.

    Xu S Z, Ruan B Y, Zhou H, et al. 1997. Numerical modeling of 3-D terrain effect on MT field.ScienceinChinaSeriesD:EarthSciences, 40(3): 269-275.

    Xu Y X, Wang J Y. 2000. Power spectrum estimation for magnetotelluric signal based on continuous wavelet transform.ChineseJ.Geophys. (in Chinese), 43(5): 677-683.

    Xu Z F, Wu X P. 2010. Electromagnetic fields excited by the horizontal electrical dipole on the surface of the ionosphere-homogeneous earth model.ChineseJ.Geophys. (in Chinese), 53(10): 2497-2506, 10.3969/j.issn.0001-5733.2010.10.023.

    Yan J B, Liu G Z. 2007. Like-impulse electromagnetic noise processing based wavelet transform.CoalGeology&Exploration(in Chinese), 35(5): 61-65.

    Yan S, Chen M S. 1996. Static offset and correction in frequency domain electromagnetic sounding.OilGeophysicalProspecting(in Chinese), 31(2): 238-247. Yan S, Chen M S. 2004. Shadow and source overprint effects in CSAMT sounding.OilGeophysicalProspecting(in Chinese), 39(S1): 8-10.

    Yang S, Bao G S, Zhang S Y. 2001. The correction to aberrant appearance resistivity curve by using impedance phase data in magnetotelluric method.Geology&Prospecting(in Chinese), 37(6): 42-45.Ye Y X, Deng J Z, Li M, et al. 2011. Application status and vistas of electromagnetic methods to deep ore prospecting.ProgressinGeophysics(in Chinese), 26(1): 327-334, doi: 10.3969/j.issn.10042903.2011.01.039.

    Yin C C, Piao H R. 1991. A study of the definition of apparent resistivity in electromagnetic sounding.Geophysical&GeochemicalExploration(in Chinese), 15(4): 290-299.

    Yue R Y, Xu Y X. 2004. Study on electric azimuth anisotropy of controlled source in frequency domain and Near-field effects.OilGeophysicalProspecting(in Chinese), 39(3): 342-347.

    Zhan S Q, Qian M P, Feng J J. 2011. The CSAMT near-field correction method for electric field apprarent resistivity forward iterative fitting.Geophysical&GeochemicalExploration(in Chinese), 35(5): 663-665.

    Zhang Y, Paulson K V. 1997. Enhancement of signal-to-noise ratio in natural-source transient magnetotelluric data with wavelet transform.PureandAppliedGeophysics, 149(2): 405-419.

    Zhao G Z, Wang L F, Tang J, et al. 2010. New experiments of CSELF electromagnetic method for earthquake monitoring.ChineseJ.Geophys. (in Chinese), 53(3): 479-486, doi: 10.3969/j.issn.0001-5733.2010.03.002.

    Zhdanov M S, Lee S K, Yoshioka K. 2006. Integral equation method for 3D modeling of electromagnetic fields in complex structures with inhomogeneous background conductivity.Geophysics, 71(6): G333-G345.

    Zhdanov M S. 2010. Electromagnetic geophysics: Notes from the past and the road ahead.Geophysics, 75(5): 75A49-75A66.

    Zhen Q H, Di Q Y, Liu H B. 2013. Key technology study on CSAMT transmitter with excitation control.ChineseJ.Geophys. (in Chinese), 56(11): 3751-3760, doi: 10.6038/cjg20131116.

    Zhou H F, Liu C, Shi K F. 2003. A review of study on geothermal resources exploration.ProgressinGeophysics(in Chinese), 18(4): 656-661, doi: 10.3969/j.issn.1004-2903.2003.04.014.

    Zhuo X J, Lu J X. 2010. Application and prospect of WFEM to resource exploration and earthquake predication.ShipScienceandTechnology(in Chinese), 32(6): 3-7.

    Zhuo X J, Lu J X, Zhao G Z, et al. 2011. The extremely low frequency engineering project using WFEM for underground exploration.EngineeringScience(in Chinese), 13(9): 42-50.

    附中文參考文獻(xiàn)

    蔡劍華, 湯井田, 王先春. 2011. 基于經(jīng)驗(yàn)?zāi)B(tài)分解的大地電磁資料人文噪聲處理. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 42(6): 1786-1790.

    曹昌祺. 1978. 水平分層大地的交流視電阻率. 地球物理學(xué)報(bào), 21(3): 248-261.

    陳海燕, 魏文博, 景建恩等. 2012. 廣義 S 變換及其在大地電磁測深數(shù)據(jù)處理中的應(yīng)用. 地球物理學(xué)進(jìn)展, 27(3): 988-996, doi: 10.6038/j.issn.1004-2903.

    陳明生, 閆述. 2005. CSAMT勘探中場區(qū)、記錄規(guī)則、陰影及場源復(fù)印效應(yīng)的解析研究. 地球物理學(xué)報(bào), 48(4): 951-958.

    陳衛(wèi)營, 薛國強(qiáng). 2014. 接地導(dǎo)線源電磁場全域有效趨膚深度. 地球物理學(xué)報(bào), 57(7): 2314-2320, doi: 10.6038/cjg20140725.

    陳小斌, 趙國澤. 2009. 關(guān)于人工源極低頻電磁波發(fā)射源的討論—均勻空間交流點(diǎn)電流源的解. 地球物理學(xué)報(bào), 52(8): 2158-2164, doi: 10.3969/j.issn.0001-5733.2009.08.026.

    陳小斌, 蔡軍濤, 王立鳳等. 2014. 大地電磁資料精細(xì)處理和二維反演解釋技術(shù)研究(四)——阻抗張量分解的多測點(diǎn)-多頻點(diǎn)統(tǒng)計(jì)成像分析. 地球物理學(xué)報(bào), 57(6): 1946-1957, doi: 10.6038/cjg20140625.

    程德福, 王君, 李秀平等. 2004. 混場源電磁法儀器研制進(jìn)展. 地球物理學(xué)進(jìn)展, 19(4): 778-781, doi: 10.3969/j.issn.1004-2903.2004.04.012.

    程輝, 李帝銓. 2014. 一種新的水平電偶極子電磁波場區(qū)劃分方法. 中南大學(xué)學(xué)報(bào) (自然科學(xué)版), 45(7): 2250-2258.

    底青云, 方廣有, 張一鳴. 2013. 地面電磁探測系統(tǒng)(SEP)研究. 地球物理學(xué)報(bào), 56(11): 3629-3639, doi: 10.6038/cjg20131104.

    底青云, 王妙月, 王若等. 2008. 長偶極大功率可控源電磁波響應(yīng)特征研究. 地球物理學(xué)報(bào), 51(6): 1917-1928.

    馮兵, 王珺璐, 周祥文等. 2013. CSAMT探測中電場Ex全區(qū)視電阻率定義及應(yīng)用. 煤田地質(zhì)與勘探, 41(6): 78-82.

    付長民, 底青云, 許誠等. 2012. 電離層影響下不同類型源激發(fā)的電磁場特征. 地球物理學(xué)報(bào), 55(12): 3958-3968, doi: 10.6038/j.issn.0001-5733.2012.12.008.

    高文. 1991. 大地電磁感應(yīng)的場源效應(yīng). 地球物理學(xué)報(bào), 34(2): 210-215. 何繼善. 1991. 可控源音頻大地電磁法. 長沙: 中南工業(yè)大學(xué)出版社. 何繼善. 2010a. 三元素集合中的自封閉加法與 2n系列偽隨機(jī)信號編碼. 中南大學(xué)學(xué)報(bào) (自然科學(xué)版), 41(2): 632-627.

    何繼善. 2010b. 廣域電磁測深法研究. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 41(3): 1065-1072.

    何繼善. 2010c. 廣域電磁法和偽隨機(jī)信號電法. 北京: 高等教育出版社.

    何蘭芳, 王緒本, 王成祥. 1999. 應(yīng)用小波分析提高M(jìn)T資料信噪比. 成都理工學(xué)院學(xué)報(bào), 26(3): 299-302.

    胡艷芳, 何繼善. 2015. 強(qiáng)干擾礦區(qū)廣域電磁法與CSAMT法對比[碩士論文]. 長沙:中南大學(xué).

    黃皓平, 樸化榮. 1992. 水平多層大地上垂直磁偶極頻率測深的全波視電阻率. 地球物理學(xué)報(bào), 35(3): 389-395.

    黃力軍, 陸桂福, 劉瑞德等. 2004. 電磁測深方法在深部地?zé)豳Y源調(diào)查中的應(yīng)用. 物探與化探, 28(6): 493-495.

    景建恩, 魏文博, 陳海燕等. 2012. 基于廣義S變換的大地電磁測深數(shù)據(jù)處理. 地球物理學(xué)報(bào), 55(12): 4015-4022, doi: 10.6038/j.issn.00015733.

    雷達(dá), 張國鴻, 黃高元等. 2014. 張量可控源音頻大地電磁法的應(yīng)用實(shí)例. 工程地球物理學(xué)報(bào), 11(3): 286-294.

    李帝銓, 底青云, 王妙月. 2011. “地-電離層”模式有源電磁場一維正演. 地球物理學(xué)報(bào), 54(9): 2375-2388, doi: 10.3969/j.issn.0001-5733.2011.09.021.

    廖忠禮, 張予杰, 陳文彬等. 2006. 地?zé)豳Y源的特點(diǎn)與可持續(xù)開發(fā)利用. 中國礦業(yè), 15(10): 8-11.

    林品榮, 郭鵬, 石福升等. 2010. 大深度多功能電磁探測技術(shù)研究. 地球?qū)W報(bào), 31(2): 149-154.

    林品榮, 鄭采君, 石福升等. 2007. 電磁法綜合探測系統(tǒng)研究. 地質(zhì)學(xué)報(bào), 80(10): 1539-1548.

    林威. 2009. CSAMT法過渡區(qū)電磁場的特征. 物探與化探, 33(2): 148-150.

    劉宏, 何蘭芳, 王緒本等. 2004. 小波分析在 MT 去噪處理中的適定性. 石油地球物理勘探, 39(3): 338-341.

    羅維斌, 李慶春, 湯井田. 2012. 編碼電磁測深. 地球物理學(xué)報(bào), 55(1): 341-349, doi: 10.6038/j.issn.0001-5733.2012.01.035.

    孟慶奎. 2013. 張量CSAMT數(shù)據(jù)采集與處理技術(shù)初步研究[碩士論文]. 北京: 中國地質(zhì)科學(xué)院.

    宋守根, 湯井田, 何繼善. 1995. 小波分析與電磁測深中靜態(tài)效應(yīng)的識別, 分離及壓制. 地球物理學(xué)報(bào), 38(1): 120-128.

    湯井田, 何繼善. 1993. 水平電偶極子源的場分量等效電阻率. 中南礦冶學(xué)院學(xué)報(bào), 24(2): 137-142.

    湯井田, 何繼善. 1994a. 水平電偶源頻率測深中全區(qū)視電阻率定義的新方法. 地球物理學(xué)報(bào), 37(4): 543-552.

    湯井田, 何繼善. 1994b. 水平多層介質(zhì)上水平電偶源頻率電磁測深的阻抗實(shí)部等效電阻率. 物探與化探, 18(2): 92-96.

    湯井田, 何繼善. 2005. 可控源音頻大地電磁法及其應(yīng)用. 長沙: 中南大學(xué)出版社.

    湯井田, 化希瑞, 曹哲民等. 2008. Hilbert-Huang 變換與大地電磁噪聲壓制. 地球物理學(xué)報(bào), 51(2): 603-610.

    湯井田, 蔡劍華, 任政勇等. 2009. Hilbert-Huang 變換與大地電磁信號的時頻分析. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 40(5): 1399-1405.

    湯井田, 周聰, 張林成. 2011. CSAMT電場y方向視電阻率的定義及研究. 吉林大學(xué)學(xué)報(bào)(地球科學(xué)版), 41(2): 552-558.

    湯井田, 李晉, 肖曉等. 2012a. 數(shù)學(xué)形態(tài)濾波與大地電磁噪聲壓制. 地球物理學(xué)報(bào), 55(5): 1784-1793, doi: 10.6038/j.issn.0001-5733.2012.05.036.

    湯井田, 徐志敏, 肖曉等. 2012b. 廬樅礦集區(qū)大地電磁測深強(qiáng)噪聲的影響規(guī)律. 地球物理學(xué)報(bào), 55(12): 4147-4159, 10.6038/j.issn.0001-5733.2012.12.027.

    湯井田, 葛偉男. 2012. 三維 CSAMT 中的陰影和場源附加效應(yīng). 物探化探計(jì)算技術(shù), 34(1): 19-26.

    湯井田, 張弛, 肖曉等. 2013a. 大地電磁阻抗估計(jì)方法對比. 中國有色金屬學(xué)報(bào), 23(9): 2351-2358.

    湯井田, 周聰, 肖曉. 2013b. 復(fù)雜介質(zhì)條件下CSAMT最小發(fā)收距的選擇. 中國有色金屬學(xué)報(bào), 23(6): 1681-1693.

    王輝, 魏文博, 金勝等. 2014. 基于同步大地電磁時間序列依賴關(guān)系的噪聲處理. 地球物理學(xué)報(bào), 57(2): 531-545, doi: 10.6038/cjg20140218.

    王若, 王妙月, 盧元林. 2004. 高山峽谷區(qū)CSAMT觀測系統(tǒng)研究.地球物理學(xué)進(jìn)展, 19(1): 125-130, doi: 10.3969/j.issn.1004-2903.2004.01.018.

    王若, 底青云, 王妙月等. 2009. 用積分方程法研究源與勘探區(qū)之間的三維體對CSAMT觀測曲線的影響. 地球物理學(xué)報(bào), 52(6): 1573-1582, doi: 10.3969/j.issn.0001-5733.2009.06.019. 王書明, 王家映. 2004. 高階統(tǒng)計(jì)量在大地電磁測深數(shù)據(jù)處理中的應(yīng)用研究. 地球物理學(xué)報(bào), 47(5): 928-934.

    王顯祥, 底青云, 許誠. 2014. CSAMT的多偶極子源特征與張量測量. 地球物理學(xué)報(bào), 57(2): 651-661, doi: 10.6038/cjg20140228. 夏時斌, 閔剛, 邱林等. 2012. 大地電磁 (MT) 法在深部礦產(chǎn)勘探中的分辨率探討. 工程地球物理學(xué)報(bào), 9(5): 531-536.

    徐義賢, 王家映. 2000. 基于連續(xù)小波變換的大地電磁信號譜估計(jì)方法. 地球物理學(xué)報(bào), 43(5): 677-683.

    徐志鋒, 吳小平. 2010. 電離層-均勻地球模型中地表水平電偶極子激發(fā)的電磁場. 地球物理學(xué)報(bào), 53(10): 2497-2506, 10.3969/j.issn.0001-5733.2010.10.02.

    嚴(yán)家斌, 劉貴忠. 2007. 基于小波變換的脈沖類電磁噪聲處理. 煤田地質(zhì)與勘探, 35(5): 61-65.

    閻述, 陳明生. 1996. 頻率域電磁測深的靜態(tài)偏移及校正方法. 石油地球物理勘探, 31(2): 238-247.

    閆述, 陳明生. 2004. CSAMT測深中的陰影和場源復(fù)印效應(yīng)問題. 石油地球物理勘探, 39(增刊): 8-10.

    楊生, 鮑光淑, 張少云. 2001. MT 法中利用阻抗相位資料對畸變視電阻率曲線的校正. 地質(zhì)與勘探, 37(6): 42-45.

    葉益信, 鄧居智, 李曼等. 2011. 電磁法在深部找礦中的應(yīng)用現(xiàn)狀及展望. 地球物理學(xué)進(jìn)展, 26(1): 327-334, doi: 10.3969/j.issn.10042903.2011.01.039. 殷長春, 樸化榮. 1991. 電磁測深法視電阻率定義問題的研究. 物探與化探, 15(4): 290-299.

    岳瑞永, 徐義賢. 2004. 可控源頻率域電方位各向異性與近場效應(yīng)研究. 石油地球物理勘探, 39(3): 342-347.

    詹少全, 錢美平, 馮戔戔. 2011. CSAMT 全區(qū)視電阻率電場正演迭代擬合近場校正方法. 物探與化探, 35(5): 663-665.

    趙國澤, 王立鳳, 湯吉等. 2010. 地震監(jiān)測人工源極低頻電磁技術(shù)(CSELF)新試驗(yàn). 地球物理學(xué)報(bào), 53(3): 479-486, doi: 10.3969/j.issn.0001-5733.2010.03.002.

    真齊輝, 底青云, 劉漢北. 2013. 勵磁控制的CSAMT發(fā)送機(jī)若干技術(shù)研究. 地球物理學(xué)報(bào), 56(11): 3751-3760, doi: 10.6038/cjg20131116.

    周厚芳, 劉闖, 石昆法. 2004. 地?zé)豳Y源探測方法研究進(jìn)展. 地球物理學(xué)進(jìn)展, 18(4): 656-661, doi: 10.3969/j.issn.1004-2903.2003.04.014.

    卓賢軍, 陸建勛. 2010. “極低頻探地工程”在資源探測和地震預(yù)測中的應(yīng)用與展望. 艦船科學(xué)技術(shù), 32(6): 3-7.

    卓賢軍, 陸建勛, 趙國澤等. 2011. 極低頻探地(WEM)工程. 中國工程科學(xué), 13(9): 42-50.

    (本文編輯 何燕)

    Frequency-domain electromagnetic methods for exploration of the shallow subsurface: A review

    TANG Jing-Tian1,2, REN Zheng-Yong1,2, ZHOU Cong1,2, ZHANG Lin-Cheng1,2, YUAN Yuan1,2, XIAO Xiao1,2

    1KeyLaboratoryofMetallogenicPredictionofNon-FerrousMetalsandGeologicalEnvironmentMonitor,MinistryofEducation,CentralSouthUniversity,Changsha410083,China2InstituteofAppliedGeophysics,SchoolofGeosciencesandInfo-Physics,CentralSouthUniversity,Changsha410083,China

    Four electromagnetic (EM) methods are widely used in exploration of the shallow subsurface (above depth 2000 meters), which are the audio-frequency magnetotellurics (AMT), radio-magnetotellurics (RMT), controlled source audio-frequency magnetotellurics (CSAMT) and the wide-field electromagnetic methods (WFEM). The purpose of this review is to present the advances of these technologies from five aspects, i.e. data collection, data processing algorithms, forward modeling, inversion algorithms and case studies. We also point out their current difficulties and suggest the potential development trends.To address the issues above, we have made a detailed investigation to the recent literature and chosen the most significant papers published.(1) Data acquisition and instrument. Along with the growing power of EM data collection instruments, the tensor measurement has been commonly adopted. In order to obtain more reliable data, instrument array or multi-instruments have been widely used in the field. By large funds for developing instruments, the quality of domestic instruments is gradually approaching that of the foreign equipment. (2) Data processing algorithms. The time series of AMT data usually contains different types of noise, such as artificial square waves and industry noise around target zones. Therefore, the signal-to-noise ratio of AMT data is usually low. In order to obtain useful signals, we can utilize the robust digital morphological filtering technology. In transforming the time domain signal into the frequency domain, the Hilbert Huang transform is a suitable choice. Recently, the main reliable approach to estimate the tensor impedance is simultaneously using multi-channel multi-frequency data. As for the CSAMT and RMT methods, the state of art in the data acquisition is to identify near field and far field zones. A new idea is the wide-band electromagnetic acquisition technique, which presents a unified way to compute the apparent resistivity in both near field and far field zones. Compared to the AMT, RMT and CSAMT methods, the WFEM method can work better in the near field zone. (3) Forward modeling and inversion. In order to compute the EM responses for complicated AMT/RMT/CSAMT problems, the finite element method has become the indispensable tool. The edge-based finite-element methods have replaced the nodal-based finite-element methods. To have a capability of dealing with complicated topography, the unstructured grids instead of traditional structured grids have become more important. The high efficiency of direct solvers such as MUMPS and Pardiso have played an important role in solving multi-source EM problems. As for the inversion, to simulate 3D cases, two inversion algorithms, which are the nonlinear conjugate gradient method and Gauss-Newton algorithm, have been applied to successfully performing 3D inversion of AMT/RMT/CSAMT data. To further improve the inversion quality, the joint inversion scheme becomes more popular. (4) Applications. The above AMT/RMT/CSAMT methods are widely used in mineral resources exploration, geothermal resources or groundwater surveys, and environmental and engineering geophysics. In inverting the AMT/RMT/CSAMT data, we should analyze the possible dimensionality of the subsurface structure. When they show a 2D structure, the results by 2D inversion are reliable. However, for 3D cases, inversion results and geological interpretations must be based on 3D inversion algorithms.

    Shallow subsurface exploration; Frequency-domain electromagnetic methods; AMT; RMT; CSAMT; WFEM

    國家自然科學(xué)基金(41104071,41574120,41174105),國家高技術(shù)研究發(fā)展計(jì)劃(863)項(xiàng)目(2014AA06A602),國家重點(diǎn)基礎(chǔ)研究發(fā)展規(guī)劃(青年973)項(xiàng)目(2015CB060201)資助.

    湯井田,男,1965年生,教授,博士生導(dǎo)師,主要研究方向:電磁法正反演研究及相關(guān)數(shù)據(jù)處理.E-mail:jitang@mail.csu.edu.cn

    10.6038/cjg20150807.

    10.6038/cjg20150807

    P631

    2014-12-22,2015-06-25收修定稿

    湯井田, 任政勇,周聰?shù)? 2015. 淺部頻率域電磁勘探方法綜述.地球物理學(xué)報(bào),58(8):2681-2705,

    Tang J T, Ren Z Y, Zhou C,et al. 2015. Frequency-domain electromagnetic methods for exploration of the shallow subsurface: A review.ChineseJ.Geophys. (in Chinese),58(8):2681-2705,doi:10.6038/cjg20150807.

    猜你喜歡
    場源反演電磁
    例談求解疊加電場的電場強(qiáng)度的策略
    基于深度展開ISTA網(wǎng)絡(luò)的混合源定位方法
    信號處理(2022年10期)2022-11-16 00:50:56
    反演對稱變換在解決平面幾何問題中的應(yīng)用
    基于矩陣差分的遠(yuǎn)場和近場混合源定位方法
    三維多孔電磁復(fù)合支架構(gòu)建與理化表征
    基于低頻軟約束的疊前AVA稀疏層反演
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    掌握基礎(chǔ)知識 不懼電磁偏轉(zhuǎn)
    一種識別位場場源的混合小波方法
    疊前同步反演在港中油田的應(yīng)用
    国产精品国产三级国产av玫瑰| 久久韩国三级中文字幕| 久久久a久久爽久久v久久| 日本色播在线视频| 精品人妻熟女av久视频| 欧美三级亚洲精品| 一级毛片久久久久久久久女| 波多野结衣巨乳人妻| 国产精品人妻久久久影院| 亚洲av不卡在线观看| 久久久午夜欧美精品| 国产爱豆传媒在线观看| 人妻少妇偷人精品九色| 久久精品91蜜桃| 亚洲高清免费不卡视频| 精品少妇黑人巨大在线播放 | 亚洲精品日韩av片在线观看| 国产片特级美女逼逼视频| 国产三级中文精品| 日韩欧美精品免费久久| 免费无遮挡裸体视频| 嫩草影院新地址| 国产午夜精品论理片| 国产一区二区亚洲精品在线观看| 久久久a久久爽久久v久久| 亚洲综合色惰| 国产久久久一区二区三区| 亚洲精品日韩av片在线观看| 久久久国产成人精品二区| av中文乱码字幕在线| 十八禁国产超污无遮挡网站| 一区二区三区高清视频在线| 亚洲欧美清纯卡通| 在线a可以看的网站| 性插视频无遮挡在线免费观看| 日韩,欧美,国产一区二区三区 | 久99久视频精品免费| 我的老师免费观看完整版| 久久国内精品自在自线图片| 亚洲成人av在线免费| 久久天躁狠狠躁夜夜2o2o| av天堂中文字幕网| 男女视频在线观看网站免费| 国产精品一区二区三区四区免费观看 | 国内揄拍国产精品人妻在线| 欧美极品一区二区三区四区| 欧美bdsm另类| 国产精品久久电影中文字幕| 一a级毛片在线观看| 亚洲久久久久久中文字幕| 国内揄拍国产精品人妻在线| 老司机福利观看| 色噜噜av男人的天堂激情| 蜜桃久久精品国产亚洲av| 国产伦精品一区二区三区视频9| 色av中文字幕| 草草在线视频免费看| 老师上课跳d突然被开到最大视频| 日韩欧美一区二区三区在线观看| 午夜精品国产一区二区电影 | 波多野结衣高清无吗| 久久久久性生活片| 亚洲真实伦在线观看| 国产亚洲欧美98| 日本在线视频免费播放| 校园春色视频在线观看| 18禁裸乳无遮挡免费网站照片| 变态另类丝袜制服| 国产三级在线视频| h日本视频在线播放| 啦啦啦啦在线视频资源| 一进一出好大好爽视频| 成人一区二区视频在线观看| 99九九线精品视频在线观看视频| 特级一级黄色大片| 91久久精品国产一区二区成人| 免费av不卡在线播放| 在现免费观看毛片| 日本五十路高清| 亚洲欧美日韩高清在线视频| 九九爱精品视频在线观看| 真实男女啪啪啪动态图| 秋霞在线观看毛片| 99国产精品一区二区蜜桃av| 我的女老师完整版在线观看| 午夜福利在线在线| 亚洲专区国产一区二区| 久久热精品热| 久久久久久久久久久丰满| 亚洲无线在线观看| .国产精品久久| 中文字幕精品亚洲无线码一区| 97人妻精品一区二区三区麻豆| 级片在线观看| 丰满乱子伦码专区| 一个人看的www免费观看视频| 亚洲欧美日韩高清专用| 自拍偷自拍亚洲精品老妇| 国产久久久一区二区三区| 国产黄a三级三级三级人| 日本熟妇午夜| 在现免费观看毛片| 亚洲精品日韩av片在线观看| 日本一本二区三区精品| 噜噜噜噜噜久久久久久91| 欧美不卡视频在线免费观看| 日本欧美国产在线视频| 亚洲欧美日韩无卡精品| 成人美女网站在线观看视频| 男女之事视频高清在线观看| 全区人妻精品视频| 夜夜夜夜夜久久久久| av免费在线看不卡| 亚洲,欧美,日韩| 国产亚洲精品久久久com| 亚洲自拍偷在线| 国产在线精品亚洲第一网站| 欧美性感艳星| 日本-黄色视频高清免费观看| 成年免费大片在线观看| 亚洲av美国av| 国产成人freesex在线 | 美女被艹到高潮喷水动态| 亚洲av熟女| 亚洲精品亚洲一区二区| 欧美极品一区二区三区四区| 国国产精品蜜臀av免费| 国产精品1区2区在线观看.| 最近2019中文字幕mv第一页| 亚洲欧美成人综合另类久久久 | 91午夜精品亚洲一区二区三区| 蜜臀久久99精品久久宅男| 热99在线观看视频| 国产精品国产三级国产av玫瑰| 不卡一级毛片| 男人狂女人下面高潮的视频| 国产欧美日韩精品亚洲av| 黄色一级大片看看| 一级黄片播放器| 亚洲自偷自拍三级| 亚洲av二区三区四区| 国产精品一二三区在线看| 天堂动漫精品| 校园春色视频在线观看| 色综合站精品国产| 淫妇啪啪啪对白视频| 五月玫瑰六月丁香| 国产av麻豆久久久久久久| 日韩av在线大香蕉| 神马国产精品三级电影在线观看| 国产综合懂色| 一区二区三区四区激情视频 | 国产高清视频在线观看网站| 国产大屁股一区二区在线视频| 欧美性猛交╳xxx乱大交人| 亚洲电影在线观看av| 99热这里只有精品一区| 我要搜黄色片| 日韩av不卡免费在线播放| 国产精品一区二区免费欧美| 69av精品久久久久久| 亚洲精品在线观看二区| 搡老熟女国产l中国老女人| 悠悠久久av| 亚洲欧美日韩无卡精品| 婷婷六月久久综合丁香| 亚洲内射少妇av| 午夜福利在线观看吧| 少妇丰满av| 国产午夜精品论理片| 国语自产精品视频在线第100页| 成人二区视频| 狂野欧美激情性xxxx在线观看| 精品久久久精品久久久| 日本wwww免费看| 看免费成人av毛片| 男女啪啪激烈高潮av片| 国产色婷婷99| 春色校园在线视频观看| 最近2019中文字幕mv第一页| 亚洲国产欧美在线一区| 久久免费观看电影| 视频区图区小说| 日日撸夜夜添| 一区二区av电影网| 精品国产一区二区久久| 一级毛片 在线播放| 日韩中文字幕视频在线看片| 七月丁香在线播放| 日本-黄色视频高清免费观看| 欧美一级a爱片免费观看看| 中文乱码字字幕精品一区二区三区| 亚洲国产欧美在线一区| 久久99一区二区三区| 免费少妇av软件| 中文乱码字字幕精品一区二区三区| 又粗又硬又长又爽又黄的视频| 寂寞人妻少妇视频99o| 久久国产精品男人的天堂亚洲 | 国产永久视频网站| 一区在线观看完整版| 少妇被粗大猛烈的视频| 国产 一区精品| 久久国产乱子免费精品| 美女主播在线视频| 日韩精品有码人妻一区| 一级毛片久久久久久久久女| 老司机影院成人| 青青草视频在线视频观看| 69精品国产乱码久久久| 亚洲精品国产av蜜桃| 最黄视频免费看| 女的被弄到高潮叫床怎么办| 国产真实伦视频高清在线观看| 亚洲精品成人av观看孕妇| 日韩一本色道免费dvd| 最近的中文字幕免费完整| 这个男人来自地球电影免费观看 | 国产精品一区二区性色av| 18禁动态无遮挡网站| 六月丁香七月| 国产成人精品婷婷| 免费高清在线观看视频在线观看| 国产欧美亚洲国产| 老女人水多毛片| 欧美精品亚洲一区二区| 久久久久久久精品精品| 亚洲国产av新网站| 在线看a的网站| 91精品一卡2卡3卡4卡| 欧美 日韩 精品 国产| 国产免费一区二区三区四区乱码| 国产男人的电影天堂91| 午夜影院在线不卡| 热re99久久国产66热| 桃花免费在线播放| 日韩av不卡免费在线播放| 欧美日本中文国产一区发布| av在线app专区| 少妇 在线观看| 国精品久久久久久国模美| 亚洲精品色激情综合| 免费人妻精品一区二区三区视频| 国产女主播在线喷水免费视频网站| 国产精品三级大全| 国产一区二区在线观看日韩| 九九在线视频观看精品| 久久毛片免费看一区二区三区| 久热久热在线精品观看| 大香蕉97超碰在线| 青春草亚洲视频在线观看| 久久久久精品久久久久真实原创| 日本黄大片高清| 欧美激情国产日韩精品一区| 老女人水多毛片| 久久精品久久久久久久性| 国语对白做爰xxxⅹ性视频网站| av一本久久久久| 精品国产一区二区久久| 久久热精品热| 人人妻人人澡人人爽人人夜夜| 午夜免费观看性视频| 在线亚洲精品国产二区图片欧美 | √禁漫天堂资源中文www| 免费播放大片免费观看视频在线观看| 天天操日日干夜夜撸| 国产黄片美女视频| 人人妻人人看人人澡| 极品少妇高潮喷水抽搐| 国产69精品久久久久777片| 亚洲国产最新在线播放| 亚洲欧洲精品一区二区精品久久久 | 国产视频内射| 2021少妇久久久久久久久久久| 免费大片18禁| 久久久国产一区二区| 亚洲欧洲国产日韩| 亚洲精品久久久久久婷婷小说| 99精国产麻豆久久婷婷| 欧美人与善性xxx| 精品视频人人做人人爽| 成人二区视频| 乱系列少妇在线播放| 性色avwww在线观看| 日本av手机在线免费观看| 国产男女内射视频| 精品一区二区免费观看| 91aial.com中文字幕在线观看| 国产免费一级a男人的天堂| 久久久久久久久久人人人人人人| 高清视频免费观看一区二区| 黄色欧美视频在线观看| 日韩,欧美,国产一区二区三区| 免费播放大片免费观看视频在线观看| 一级片'在线观看视频| 久久97久久精品| 国产 精品1| 老熟女久久久| 国产伦理片在线播放av一区| 亚洲性久久影院| 免费看不卡的av| 国产欧美日韩精品一区二区| 一个人免费看片子| 十八禁高潮呻吟视频 | 欧美xxⅹ黑人| 国产有黄有色有爽视频| 老司机影院毛片| 日韩熟女老妇一区二区性免费视频| 色5月婷婷丁香| 赤兔流量卡办理| 日韩一区二区视频免费看| 又黄又爽又刺激的免费视频.| 99热这里只有精品一区| 亚洲无线观看免费| 欧美三级亚洲精品| 国产一区二区在线观看日韩| 亚洲欧洲日产国产| 亚洲婷婷狠狠爱综合网| 天堂俺去俺来也www色官网| 亚洲第一av免费看| 午夜福利在线观看免费完整高清在| 黑丝袜美女国产一区| 精品国产一区二区三区久久久樱花| 国产黄片美女视频| .国产精品久久| 成人亚洲欧美一区二区av| 国产精品一区二区在线观看99| 三级经典国产精品| 成人免费观看视频高清| 亚洲中文av在线| 三级经典国产精品| 亚洲成色77777| 国产av精品麻豆| 肉色欧美久久久久久久蜜桃| 一级毛片黄色毛片免费观看视频| 九九在线视频观看精品| 你懂的网址亚洲精品在线观看| 国产有黄有色有爽视频| av.在线天堂| 成年美女黄网站色视频大全免费 | 国产免费又黄又爽又色| 国产成人精品一,二区| 国产日韩欧美在线精品| 亚洲三级黄色毛片| 欧美97在线视频| 国产伦精品一区二区三区四那| 日韩电影二区| 97超碰精品成人国产| 国产伦在线观看视频一区| 久久久久久久久久成人| 成年人免费黄色播放视频 | 美女脱内裤让男人舔精品视频| 国产熟女欧美一区二区| 亚洲电影在线观看av| 国产免费一级a男人的天堂| 伊人亚洲综合成人网| 久久久午夜欧美精品| 在线观看免费日韩欧美大片 | 人人妻人人澡人人爽人人夜夜| 色网站视频免费| 亚洲欧美精品自产自拍| 在线观看av片永久免费下载| 亚洲av国产av综合av卡| 国产高清国产精品国产三级| 国产亚洲最大av| 人人妻人人添人人爽欧美一区卜| 日本爱情动作片www.在线观看| 美女福利国产在线| 欧美 亚洲 国产 日韩一| 日韩亚洲欧美综合| 最近的中文字幕免费完整| 亚洲av综合色区一区| 内地一区二区视频在线| 国产一区有黄有色的免费视频| 日韩一区二区视频免费看| 亚洲国产精品999| 日韩熟女老妇一区二区性免费视频| 插逼视频在线观看| 天天躁夜夜躁狠狠久久av| 国产精品偷伦视频观看了| 3wmmmm亚洲av在线观看| 最近最新中文字幕免费大全7| 国产亚洲午夜精品一区二区久久| 美女xxoo啪啪120秒动态图| 久久 成人 亚洲| 热99国产精品久久久久久7| 韩国av在线不卡| 精品久久久噜噜| 超碰97精品在线观看| 在线免费观看不下载黄p国产| 高清不卡的av网站| 国产精品熟女久久久久浪| 99热全是精品| 国产深夜福利视频在线观看| tube8黄色片| 国产精品不卡视频一区二区| 99国产精品免费福利视频| 免费av不卡在线播放| av网站免费在线观看视频| 少妇人妻一区二区三区视频| 99热网站在线观看| 噜噜噜噜噜久久久久久91| 美女福利国产在线| 国产精品久久久久久精品电影小说| 麻豆精品久久久久久蜜桃| 久久ye,这里只有精品| 久久6这里有精品| 一本一本综合久久| 国产无遮挡羞羞视频在线观看| 亚洲国产精品国产精品| 国产伦精品一区二区三区四那| 欧美日本中文国产一区发布| 永久网站在线| 亚洲天堂av无毛| 欧美精品国产亚洲| 啦啦啦中文免费视频观看日本| 婷婷色综合www| 免费看日本二区| 青春草视频在线免费观看| 久久精品国产亚洲av涩爱| 永久网站在线| av福利片在线| 高清av免费在线| 日日摸夜夜添夜夜爱| 精品一区二区三区视频在线| 欧美 日韩 精品 国产| 男女无遮挡免费网站观看| 久久久久国产网址| 国产极品天堂在线| 成人无遮挡网站| 亚洲国产日韩一区二区| 亚洲av成人精品一二三区| 日韩免费高清中文字幕av| www.av在线官网国产| 自拍欧美九色日韩亚洲蝌蚪91 | 看十八女毛片水多多多| 最近最新中文字幕免费大全7| 亚洲精品日本国产第一区| 伦精品一区二区三区| 国产免费一级a男人的天堂| av不卡在线播放| 欧美 日韩 精品 国产| 三级经典国产精品| 黑人高潮一二区| 丝袜脚勾引网站| 在线播放无遮挡| 亚洲熟女精品中文字幕| 六月丁香七月| 亚洲av成人精品一区久久| 看十八女毛片水多多多| 黄片无遮挡物在线观看| 2018国产大陆天天弄谢| 日韩免费高清中文字幕av| 亚洲av不卡在线观看| 日日摸夜夜添夜夜添av毛片| 国产成人精品无人区| 99久久中文字幕三级久久日本| 中文字幕av电影在线播放| 国产日韩欧美视频二区| 黄色欧美视频在线观看| 91久久精品电影网| 国产成人免费观看mmmm| 日日爽夜夜爽网站| 亚洲av综合色区一区| 免费看不卡的av| 国内少妇人妻偷人精品xxx网站| 一本色道久久久久久精品综合| 男男h啪啪无遮挡| 国产精品蜜桃在线观看| 人人妻人人爽人人添夜夜欢视频 | 在线观看免费日韩欧美大片 | 亚州av有码| 自线自在国产av| 五月伊人婷婷丁香| 免费看av在线观看网站| 久久99一区二区三区| 国产男人的电影天堂91| 王馨瑶露胸无遮挡在线观看| 中文字幕免费在线视频6| 丁香六月天网| 成年人免费黄色播放视频 | 菩萨蛮人人尽说江南好唐韦庄| 91精品国产国语对白视频| 精品99又大又爽又粗少妇毛片| 亚洲久久久国产精品| 亚洲欧美精品专区久久| 一级毛片久久久久久久久女| 天堂俺去俺来也www色官网| a级片在线免费高清观看视频| 久久久亚洲精品成人影院| 一级a做视频免费观看| 乱系列少妇在线播放| 亚洲丝袜综合中文字幕| 22中文网久久字幕| 一级a做视频免费观看| 熟女人妻精品中文字幕| 国产欧美日韩精品一区二区| 欧美精品人与动牲交sv欧美| 男女边吃奶边做爰视频| 国产精品一区www在线观看| 秋霞伦理黄片| 亚洲精品视频女| 偷拍熟女少妇极品色| 丝袜在线中文字幕| 亚洲中文av在线| 精品一区二区免费观看| 人妻人人澡人人爽人人| 精品久久久噜噜| 美女大奶头黄色视频| 国产精品一区二区在线不卡| 亚洲欧美中文字幕日韩二区| a 毛片基地| 午夜视频国产福利| 免费人妻精品一区二区三区视频| 视频中文字幕在线观看| 中文字幕免费在线视频6| 成人亚洲欧美一区二区av| 97超视频在线观看视频| 男人爽女人下面视频在线观看| 日韩欧美一区视频在线观看 | 国产成人a∨麻豆精品| 在线观看人妻少妇| 交换朋友夫妻互换小说| 我要看日韩黄色一级片| 亚洲精品456在线播放app| 国产白丝娇喘喷水9色精品| 美女国产视频在线观看| 中文乱码字字幕精品一区二区三区| 午夜视频国产福利| 亚洲精品日韩av片在线观看| 亚洲高清免费不卡视频| 国产色婷婷99| a级毛片在线看网站| 最近中文字幕高清免费大全6| 老熟女久久久| 中文天堂在线官网| 美女xxoo啪啪120秒动态图| 亚洲人成网站在线播| 国产精品久久久久久久电影| 有码 亚洲区| 91精品国产九色| 你懂的网址亚洲精品在线观看| 欧美另类一区| 99热这里只有精品一区| 日本黄大片高清| 国产综合精华液| 日本91视频免费播放| 91午夜精品亚洲一区二区三区| 91在线精品国自产拍蜜月| 观看免费一级毛片| 日日啪夜夜撸| 嫩草影院入口| 国产在线男女| 亚洲美女视频黄频| 国产淫片久久久久久久久| 亚洲在久久综合| 精品午夜福利在线看| 婷婷色av中文字幕| 日本-黄色视频高清免费观看| 偷拍熟女少妇极品色| 精品亚洲成a人片在线观看| 久久久亚洲精品成人影院| 久久韩国三级中文字幕| 国产欧美日韩精品一区二区| 3wmmmm亚洲av在线观看| 中文乱码字字幕精品一区二区三区| 十八禁网站网址无遮挡 | 国产亚洲一区二区精品| 99视频精品全部免费 在线| 亚洲精品国产成人久久av| 精品亚洲乱码少妇综合久久| 亚洲精品视频女| 丰满少妇做爰视频| 日日摸夜夜添夜夜添av毛片| 国产成人免费观看mmmm| 久久精品久久久久久久性| 18+在线观看网站| 国产成人免费观看mmmm| 国产69精品久久久久777片| 日本-黄色视频高清免费观看| 视频中文字幕在线观看| 夜夜爽夜夜爽视频| 最新中文字幕久久久久| 久久久久国产精品人妻一区二区| 亚洲经典国产精华液单| 成人18禁高潮啪啪吃奶动态图 | 国产男人的电影天堂91| 国产91av在线免费观看| 亚洲精品亚洲一区二区| videossex国产| 一个人看视频在线观看www免费| 2018国产大陆天天弄谢| 亚洲av在线观看美女高潮| 中文欧美无线码| 熟女人妻精品中文字幕| 国产av精品麻豆| 国产av国产精品国产| 国产精品人妻久久久久久| 黑人高潮一二区| 高清毛片免费看| 国产一区二区三区av在线| 午夜激情福利司机影院| 亚洲美女黄色视频免费看| 成年av动漫网址| 国产欧美另类精品又又久久亚洲欧美| 久久久国产欧美日韩av| 日韩中文字幕视频在线看片| 日本欧美视频一区| 一边亲一边摸免费视频| 日韩视频在线欧美| 大片电影免费在线观看免费| 一级毛片aaaaaa免费看小| 黑人猛操日本美女一级片| 国产精品不卡视频一区二区| 黄色怎么调成土黄色|