專利名稱：一種大地電磁阻抗測量方法
技術(shù)領(lǐng)域：
本發(fā)明屬于地球物理勘探技術(shù)，具體涉及大地電磁勘測信號數(shù)據(jù)處理，是一種大地電磁阻抗測量方法。
背景技術(shù)：
地球物理勘探大地電磁勘測的方法，通常是通過大地電磁(MT)電場或磁場傳感器采集數(shù)據(jù)，輸出的電場或磁場信號通過濾波和放大電路進(jìn)入A/D轉(zhuǎn)換電路輸出數(shù)字信號，經(jīng)濾波、抽樣、增益改正等處理，輸出MT信號計(jì)算大地電磁阻抗，并將計(jì)算結(jié)果用圖形方式顯示出來，用來反映勘探目標(biāo)的情況。
要獲得可靠的MT測量數(shù)據(jù)，MT阻抗計(jì)算是關(guān)鍵。由于MT法采用天然電磁場信號作為場源，MT法所測量的兩個正交電場分量和磁場分量不可避免地受到各種干擾。在MT阻抗計(jì)算中，以磁場信號作為輸入信號，以電場信號作為輸出信號。MT阻抗計(jì)算相當(dāng)于求取兩個輸入信號激勵(正交磁場)兩個輸出信號響應(yīng)(正交電場)時的傳遞函數(shù)。最早的阻抗計(jì)算方法是最小二乘法，該方法假設(shè)MT信號的噪音服從高斯分布。但上述假設(shè)在干擾較大的情況很難滿足，飛點(diǎn)時時出現(xiàn)在MT測量結(jié)果中。針對上述情況，Sutarno和Vozoff于1991年提出了大地電磁阻抗的魯棒(Robust)方法(D.Sutarno and K.Vozoff，Phase-smoothed robust M-estimation of magnetotelluricimpedance functions，Geophysics，v56，n12，p.1999-2007，1991)。Egbert和Livelybrooks后來進(jìn)一步發(fā)展了魯棒方法(Gary D.Egbert and Dean W.Livelybrooks，Single station magnetotelluric impedance estimationCoherence weighting and the regression M-estimate，Geophysics，V61，n4，p964-970，1996.)。魯棒方法有2個缺點(diǎn)一、MT信號噪音服從高斯分布時，魯棒方法計(jì)算結(jié)果不是最優(yōu)的，是次優(yōu)；二、計(jì)算量大，計(jì)算時間比最小二乘法多了數(shù)倍。另一方面，當(dāng)干擾很大時，魯棒方法也無能為力。因?yàn)轸敯粲?jì)算法要求所有測量數(shù)據(jù)中至少有50％以上的點(diǎn)不是飛點(diǎn)，才能獲得比較可靠的測量結(jié)果。在上述情況下，對測量的MT信號進(jìn)行濾波，就成了干擾很大時很有前途的一種數(shù)據(jù)處理方法。然而，MT信號是隨機(jī)信號，在測量的整個頻段上都以隨機(jī)信號存在，用普通的濾波方法根本不可行。因?yàn)楦蓴_和MT信號一樣表現(xiàn)為隨機(jī)信號。為了解決上述問題，Trad和Travassos于2000年提出用小波閥值濾波法壓制MT磁信號中的干擾(Daniel O.Trad and JandryM.Travassos，Wavelet filtering of magnetotel luric data，Geophysics，V65，n2，P482-491，2000.)。具體方法是用小波變換對MT信號進(jìn)行多尺度分解后，對一個或多個尺度上的小波分解結(jié)果采用閥值濾波。上述方法也有三個缺點(diǎn)，一是MT測量中大干擾出現(xiàn)的概率很小，90％以上的測量數(shù)據(jù)用魯棒估計(jì)法或最小二乘法處理就可以滿足要求，在上述情況下采用小波壓噪，不僅浪費(fèi)計(jì)算時間，而且會濾掉有用信號。因?yàn)镸T信號分布在整個測量頻段上，無論采用何種方法濾波，都會對有用信號造成影響。二是確定閥值濾波的閥值沒有可計(jì)算的依據(jù)，只是分成多個層次后通過人的經(jīng)驗(yàn)來確定閥值大小，人為因素很大，在實(shí)踐中應(yīng)用存在困難。三是計(jì)算量大，是最小二乘法所需時間的數(shù)十倍。
另外，目前對MT數(shù)據(jù)質(zhì)量判斷依據(jù)是電場信號和磁場信號的相干度。
以磁場信號為輸入信號，電場信號為輸出信號，電場和磁場的相干度Гxy2(z)定義如下Γxy2(z)=|Φxy(z)|2Φxx(z)Φyy(z)]]>式中，Φxy(z)，Φxx(z)和Φyy(z)分別表示電場和磁場信號的互功率譜、磁場信號自功率譜和電場信號自功率譜。
通常，相干度接近1，表示數(shù)據(jù)質(zhì)量好，相干度接近0，表示數(shù)據(jù)質(zhì)量差。然而，當(dāng)電場信號和磁場信號受到同源近場干擾時，相干度基本等于1，只利用相干度來判斷數(shù)據(jù)質(zhì)量有時會得到錯誤的結(jié)果。例如，在變電站、城市和礦山附近，其干擾往往體現(xiàn)為強(qiáng)同源近場干擾，這時相干度接近1，利用相干度來選擇數(shù)據(jù)就會導(dǎo)致強(qiáng)干擾數(shù)據(jù)被保留，從而導(dǎo)致最差的MT阻抗測量結(jié)果。

發(fā)明內(nèi)容
本發(fā)明目的在于提供一種克服目前采用單一MT阻抗測量技術(shù)的結(jié)果有偏差缺點(diǎn)，能根據(jù)大地電磁信號干擾的類型和數(shù)據(jù)質(zhì)量，確定MT阻抗測量步驟，能提高M(jìn)T阻抗計(jì)算的速度的一種大地電磁阻抗測量方法。
本發(fā)明采用如下技術(shù)方案大地電磁阻抗測量方法，采用通常的方法采集電場或磁場數(shù)據(jù)，通過濾波和放大后轉(zhuǎn)換為數(shù)字信號，采用以下步驟1)對采集的信號采用最小二乘法計(jì)算大地電磁阻抗，并計(jì)算頻點(diǎn)電場信號和磁場信號相干度，當(dāng)75％以上測量頻點(diǎn)相干度為0.7-1且阻抗中相位在寬頻率范圍內(nèi)均不為175-180度，進(jìn)行步驟4)，否則進(jìn)行步驟2)；2)采用魯棒法計(jì)算大地電磁阻抗，計(jì)算電場信號和磁場信號相干度，如果75％以上測量頻點(diǎn)相干度為0.7-1并且阻抗中相位在寬頻率范圍內(nèi)均不為175-180度進(jìn)行步驟4)，否則進(jìn)行步驟3)；3)根據(jù)干擾類型在時間域或頻率域自動采用小波閥值壓噪法對大地電磁信號進(jìn)行濾波壓噪，再利用魯棒法計(jì)算大地電磁阻抗；4)計(jì)算大地電磁阻抗并將結(jié)果記錄或以圖形方式顯示。
所述的小波閥值壓噪法在有固定頻率干擾時，則在頻率域進(jìn)行小波閥值壓噪。
所述的小波閥值壓噪法在有瞬時干擾時，在時間域進(jìn)行壓噪。
所述的寬頻率范圍是通過連續(xù)有多個頻點(diǎn)阻抗相位大于175度。
本發(fā)明實(shí)現(xiàn)了大地電磁阻抗計(jì)算結(jié)果最優(yōu)，避免了次優(yōu)結(jié)果和錯誤結(jié)果出現(xiàn)。有利于提高大地電磁數(shù)據(jù)的測量精度和可靠性。在大地電磁噪音服從高斯分布，或干擾比較小時，計(jì)算速度比魯棒法計(jì)算快數(shù)倍。當(dāng)魯棒估計(jì)失敗時，自動采用小波閥值法壓制干擾，比單一的魯棒估計(jì)效果好。避免了人為因數(shù)對大地電磁阻抗計(jì)算的影響。通過是否連續(xù)有多個頻點(diǎn)阻抗相位大于175度來識別近場干擾，避免了單一采用相關(guān)度判斷MT數(shù)據(jù)質(zhì)量在近場存在時判斷錯誤這一缺點(diǎn)。


圖1為本發(fā)明流程圖。
圖2為本發(fā)明數(shù)據(jù)質(zhì)量分析與評價流程圖。
具體實(shí)施例方式
本發(fā)明在將最小二乘法、魯棒方法和小波壓噪有機(jī)結(jié)合起來，分步實(shí)現(xiàn)MT阻抗的計(jì)算，獲取最優(yōu)結(jié)果。
本發(fā)明采用通常的方法采集電場或磁場數(shù)據(jù)，通過濾波和放大后轉(zhuǎn)換為數(shù)字信號，步驟是1)對采集的信號采用最小二乘法計(jì)算大地電磁阻抗，并計(jì)算頻點(diǎn)電場信號和磁場信號相干度，當(dāng)75％以上測量頻點(diǎn)相干度為0.7-1且阻抗中相位在寬頻率范圍內(nèi)均不為175-180度，則輸出MT阻抗數(shù)據(jù)，完成MT抗計(jì)算過程，否則；2)采用魯棒法計(jì)算大地電磁阻抗，計(jì)算電場信號和磁場信號相干度，如果75％以上測量頻點(diǎn)相干度為0.7-1并且阻抗中相位在寬頻率范圍內(nèi)均不為175-180度則輸出MT阻抗數(shù)據(jù)，完成MT抗計(jì)算過程，否則；3)根據(jù)干擾類型在時間域或頻率域自動采用小波閥值壓噪法對大地電磁信號進(jìn)行濾波壓噪，再利用魯棒法計(jì)算大地電磁阻抗；計(jì)算大地電磁阻抗并將結(jié)果記錄或以圖形方式顯示，所述的小波閥值壓噪法在有固定頻率干擾時，則在頻率域進(jìn)行小波閥值壓噪。
所述的小波閥值壓噪法在有瞬時干擾時，在時間域進(jìn)行壓噪。
所述的寬頻率范圍是通過連續(xù)有多個頻點(diǎn)阻抗相位大于175度。
現(xiàn)結(jié)合說明書附圖1說明本發(fā)明實(shí)現(xiàn)MT阻抗計(jì)算基本過程。
<1>輸入MT時間域數(shù)據(jù)從MT時間域數(shù)據(jù)文件中讀出MT時間序列；<2>MT阻抗最小二乘法阻抗計(jì)算對輸入的MT數(shù)據(jù)進(jìn)行最小二乘法MT阻抗計(jì)算；<3>判斷MT阻抗數(shù)據(jù)質(zhì)量，如果數(shù)據(jù)質(zhì)量符合要求，則依次執(zhí)行步驟<10>和步驟<11>，否則執(zhí)行步驟<4>；<4>對MT數(shù)據(jù)進(jìn)行魯棒方法MT阻抗計(jì)算；<5>判斷MT阻抗數(shù)據(jù)質(zhì)量，如果數(shù)據(jù)質(zhì)量符合要求，則依次執(zhí)行步驟<10>和步驟<11>，否則執(zhí)行步驟<6>；<6>對MT數(shù)據(jù)進(jìn)行小波閥值法去噪處理；<7>對MT數(shù)據(jù)進(jìn)行魯棒法MT阻抗計(jì)算；<8>判斷MT阻抗數(shù)據(jù)質(zhì)量，如果數(shù)據(jù)質(zhì)量符合要求，則依次執(zhí)行步驟<10>和步驟<11>，否則依次執(zhí)行步驟<9>和步驟<11>；<9>確認(rèn)MT數(shù)據(jù)質(zhì)量不符合要求，并報告分析結(jié)果；<10>輸出MT阻抗數(shù)據(jù)；<11>完成MT阻抗計(jì)算過程。
本發(fā)明對MT數(shù)據(jù)質(zhì)量判斷除了依據(jù)電場信號和磁場信號的相干度以外，還通過直接分析計(jì)算的大地電磁阻抗來判斷數(shù)據(jù)質(zhì)量。
當(dāng)存在同源近場干擾時，大地電磁阻抗表現(xiàn)為阻抗相位在較寬的頻率范圍內(nèi)接近180度。利用上述特點(diǎn)，可以判斷是否存在同源近場干擾，進(jìn)而判斷MT數(shù)據(jù)質(zhì)量。
現(xiàn)結(jié)合說明書附圖2，說明判斷MT數(shù)據(jù)質(zhì)量的方法。
<1>輸入MT頻率域數(shù)據(jù)和MT阻抗數(shù)據(jù)；<2>利用MT頻率域數(shù)據(jù)計(jì)算電場信號和磁場信號的相關(guān)度；<3>判斷相關(guān)度小于0.7的頻點(diǎn)是否超過總頻點(diǎn)的15％，如果為真，則依次執(zhí)行步驟<6>和步驟<7>，否則執(zhí)行下一步；<4>判斷阻抗相位連續(xù)大于175度的頻點(diǎn)是否小于總頻點(diǎn)15％，如果條件不滿足，則依次執(zhí)行步驟<6>和步驟<7>，否則執(zhí)行下一步；<5>作出MT數(shù)據(jù)質(zhì)量達(dá)到要求判斷，然后執(zhí)行步驟<7>；<6>作出MT數(shù)據(jù)質(zhì)量不符合要求判斷，然后執(zhí)行步驟<7>；<7>輸出MT數(shù)據(jù)質(zhì)量判斷結(jié)果。
通過上述方法判斷MT數(shù)據(jù)質(zhì)量，避免了單獨(dú)使用相關(guān)度無法識別近場干擾這一缺點(diǎn)。
本發(fā)明可以在V5-2000大地電磁數(shù)據(jù)采集器上實(shí)現(xiàn)。
權(quán)利要求
1.一種大地電磁阻抗測量方法，采用通常的方法采集電場或磁場數(shù)據(jù)，通過濾波和放大后轉(zhuǎn)換為數(shù)字信號，其特征在于還采用以下步驟1)對采集的信號采用最小二乘法計(jì)算大地電磁阻抗，并計(jì)算頻點(diǎn)電場信號和磁場信號相干度，當(dāng)75％以上測量頻點(diǎn)相干度為0.7-1且阻抗中相位在寬頻率范圍內(nèi)均不為175-180度，進(jìn)行步驟4)，否則進(jìn)行步驟2)；2 )采用魯棒法計(jì)算大地電磁阻抗，計(jì)算電場信號和磁場信號相干度，如果75％以上測量頻點(diǎn)相干度為0.7-1并且阻抗中相位在寬頻率范圍內(nèi)均不為175-180度進(jìn)行步驟4)，否則進(jìn)行步驟3)；3)根據(jù)干擾類型在時間域或頻率域自動采用小波閥值壓噪法對大地電磁信號進(jìn)行濾波壓噪，再利用魯棒法計(jì)算大地電磁阻抗；4)計(jì)算大地電磁阻抗并將結(jié)果記錄或以圖形方式顯示。
2.根據(jù)權(quán)利要求1所述的一種大地電磁阻抗測量方法，其特征在于所述的小波閥值壓噪法在有固定頻率干擾時，則在頻率域進(jìn)行小波閥值壓噪；在有瞬時干擾時，則在時間域進(jìn)行壓噪。
3.根據(jù)權(quán)利要求1所述的一種大地電磁阻抗測量方法，其特征在于所述的小波閥值壓噪法在有瞬時干擾時，在時間域進(jìn)行壓噪。
4.根據(jù)權(quán)利要求1所述的一種大地電磁阻抗測量方法，其特征在于所述的寬頻率范圍是連續(xù)有多個頻點(diǎn)阻抗相位大于175度。
全文摘要
本發(fā)明涉及大地電磁勘測信號數(shù)據(jù)處理技術(shù)，是一種大地電磁阻抗測量方法。采用步驟是對采集的信號采用最小二乘法計(jì)算大地電磁阻抗，并計(jì)算頻點(diǎn)電場信號和磁場信號相干度，當(dāng)75％以上測量頻點(diǎn)相干度為0.7－1且阻抗中相位在寬頻率范圍內(nèi)均不為175－180度，計(jì)算大地電磁阻抗并將結(jié)果記錄或以圖形方式顯示。否則采用魯棒法計(jì)算大地電磁阻抗或根據(jù)干擾類型在時間域或頻率域自動采用小波閥值壓噪法對大地電磁信號進(jìn)行濾波壓噪，本發(fā)明有利于提高大地電磁數(shù)據(jù)的測量精度和可靠性，避免了人為因數(shù)對大地電磁阻抗計(jì)算的影響和單一采用相關(guān)度判斷MT數(shù)據(jù)質(zhì)量在近場存在時判斷錯誤。
文檔編號G01V3/00GK101071179SQ200610078739
公開日2007年11月14日 申請日期2006年5月11日 優(yōu)先權(quán)日2006年5月11日
發(fā)明者楊輪凱, 何展翔, 何蘭芳, 陳儒軍 申請人:中國石油集團(tuán)東方地球物理勘探有限責(zé)任公司