瞬變電磁法資料解釋精度的改進(jìn)_第1頁(yè)
瞬變電磁法資料解釋精度的改進(jìn)_第2頁(yè)
瞬變電磁法資料解釋精度的改進(jìn)_第3頁(yè)
全文預(yù)覽已結(jié)束

下載本文檔

版權(quán)說(shuō)明:本文檔由用戶(hù)提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡(jiǎn)介

瞬變電磁法資料解釋精度的改進(jìn)

1瞬變電磁電性界面的識(shí)別瞬變電壓法(tem)便于操作,采集快,測(cè)量深度變化范圍大,生產(chǎn)成本低,對(duì)地形條件要求小等特點(diǎn)。在某些特殊情況下,如弱地震反射區(qū)、地形復(fù)雜地區(qū)等,可彌補(bǔ)地震勘探的不足。提高瞬變電磁法的解釋精度,擴(kuò)大其應(yīng)用范圍,一直是地球物理工作者不斷追求的目標(biāo),也是地質(zhì)、工程等給該方法提出的急待解決的課題。傳統(tǒng)的瞬變電磁法資料解釋手段單一,主要是通過(guò)視電阻率值這一參數(shù)來(lái)對(duì)地下電性分布情況作出判斷,這往往使人為因素增大。所以,有必要尋求一種更直觀地反映地下介質(zhì)分層情況的多元參數(shù)。本文就是在分析瞬變電磁在水平層狀介質(zhì)下電磁響應(yīng)的特性的基礎(chǔ)上,提出了以反射系數(shù)為特征量,進(jìn)行瞬變電磁電性界面的識(shí)別。即從野外實(shí)測(cè)的電壓衰減數(shù)據(jù)進(jìn)行一系列計(jì)算,求出反射系數(shù)序列,并以此為參數(shù)繪制斷面圖,從而判定地下電性界面。從對(duì)兩個(gè)三層模型和實(shí)測(cè)資料的處理結(jié)果看,用文中提出的方法繪制的剖面清晰直觀地顯示出了地下電性層分界面。2時(shí)間域與反射系數(shù)的關(guān)系瞬變電磁場(chǎng)的求解是通過(guò)在頻率域中將麥克斯韋方程組轉(zhuǎn)化為亥姆霍茲方程求解,然后利用域的轉(zhuǎn)換技術(shù)將頻率域的結(jié)果轉(zhuǎn)換到時(shí)間域來(lái)。對(duì)于水平層狀大地,在大回線發(fā)射諧變電流激發(fā)下,回線中心處頻率域電磁響應(yīng)為:式中a為回線半徑;I為發(fā)射電流;J1(λa)為一階貝塞爾函數(shù);μ0為磁導(dǎo)率;ω為圓頻率;dB(ω)dtdB(ω)dt為磁場(chǎng)微分量;Z(1)為水平層狀介質(zhì)地面波阻抗,由下列公式求得:式中σj為第j層的電導(dǎo)率。根據(jù)頻譜分析理論,在水平層狀介質(zhì)表面,回線源形成的瞬變電磁響應(yīng)在時(shí)間域的表達(dá)式:另外,在頻率域內(nèi),地面波阻抗與反射系數(shù)序列之間的關(guān)系可表達(dá)如下:Z(0?ω)=Z01[1+2∑m=1∞qme?2?iωt0√?m]Ζ(0?ω)=Ζ01[1+2∑m=1∞qme-2-iωt0?m]。(3)式中qm=ρm+1√?ρm√ρm+1√+ρm√qm=ρm+1-ρmρm+1+ρm,為第m個(gè)界面的反射系數(shù);ρm為第m層的電阻率值。在時(shí)間域,地面波阻抗與反射系數(shù)之間的關(guān)系,可由(3)式兩邊同時(shí)進(jìn)行域的變換后得到。式(2)、(3)是本文提出的計(jì)算反射系數(shù),進(jìn)行地下電性界面成像的兩個(gè)重要數(shù)學(xué)模型。對(duì)(3)式兩邊同時(shí)進(jìn)行付氏變換,可以得到在時(shí)間域內(nèi)的以反射系數(shù)序列為未知的線性方程組:R(t)=∑m=1∞qmPR(t)=∑m=1∞qmΡ,(4)式中R(t)=?12πi∫∞?∞R(ω)/ωe?iωtdω;式中R(t)=-12πi∫-∞∞R(ω)/ωe-iωtdω;P=?12πi∫∞?∞e?2?iωt0√ωe?iωtdω;Ρ=-12πi∫-∞∞e-2-iωt0ωe-iωtdω;R(ω)=Z(0?ω)Z01R(ω)=Ζ(0?ω)Ζ01為頻率域歸一化波阻抗。對(duì)某一時(shí)刻tj而言,R(tj)=∑m=1∞qmPR(tj)=∑m=1∞qmΡ。(5)用有限級(jí)數(shù)代替無(wú)窮級(jí)數(shù),上式可表述為:R(tj)≈∑m=1N1qmPR(tj)≈∑m=1Ν1qmΡ。(6)采用線性規(guī)劃法求解方程組(6),可得出反射系數(shù)。3反射系數(shù)剖面圖上述方法的實(shí)施步驟為:由實(shí)測(cè)或理論模型正演數(shù)據(jù)求出磁感應(yīng)強(qiáng)度微分參數(shù),并由此參數(shù)求出歸一化后的波阻抗,最終由(6)式求出反射系數(shù)。對(duì)于以上各步驟,已在計(jì)算機(jī)上編程實(shí)現(xiàn),并對(duì)理論模型和實(shí)測(cè)資料進(jìn)行了處理。設(shè)計(jì)了兩個(gè)理論模型進(jìn)行了數(shù)值模擬計(jì)算。一個(gè)是H型斷面,一個(gè)是A型斷面。圖1是兩個(gè)地?cái)嗝嫖镄詤?shù)示意圖;圖2是兩個(gè)模型的反射系數(shù)剖面圖。以計(jì)算出來(lái)的反射系數(shù)為參數(shù),繪制成經(jīng)過(guò)時(shí)深轉(zhuǎn)換后的剖面圖,如圖2所示。所用的時(shí)深轉(zhuǎn)換公式為:Vi=2μ0σi???√Vi=2μ0σi和hj=tiVi,由圖2a可看出在100m和200m處分別出現(xiàn)了一個(gè)電性分界面,由圖2b可知,在100m和500m處出現(xiàn)電性層分界。通過(guò)此圖可直觀地劃分地下介質(zhì)電性分層情況。圖3為某一測(cè)區(qū)某線的視電阻率等值線斷面圖,由這一剖面上可知:視電阻率值反映地下電性界面不明顯也不直觀,在此圖上能看出地下電性層分層情況,但只能憑感覺(jué)和經(jīng)驗(yàn)進(jìn)行分析判斷。圖4為繪制出的反射系數(shù)剖面圖,在此圖上明顯看出有3個(gè)界面,結(jié)合圖5(綜合解釋地質(zhì)剖面圖)可知:第一層深度1~50m,第二層深度50~80m,第三層深度為80~160m或更深,分別解釋為河漫灘沖積層,粘土和亞粘土層,粘土與礫石互層。4反射系數(shù)繪制的剖面識(shí)別效果文中提出用反射系數(shù)作為特征量來(lái)對(duì)地下電性界面進(jìn)行識(shí)別,給出了數(shù)學(xué)模型、實(shí)現(xiàn)步驟。經(jīng)過(guò)理論模型

溫馨提示

  • 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶(hù)所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶(hù)上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶(hù)上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶(hù)因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。

最新文檔

評(píng)論

0/150

提交評(píng)論