二維航空被動(dòng)源電磁自適應(yīng)有限元模擬及反演研究
發(fā)布時(shí)間:2021-06-21 11:08
被動(dòng)源航空電磁法以天然電磁場(chǎng)作為探測(cè)地下電性異常目標(biāo)體的信號(hào)源,采用機(jī)載和地面電磁傳感器分別測(cè)量垂直磁場(chǎng)分量與水平磁場(chǎng)分量以獲得二者之間的復(fù)數(shù)關(guān)系(傾子),并對(duì)傾子數(shù)據(jù)進(jìn)行反演獲得電性異常目標(biāo)體的空間分布特征和物理屬性。該方法靈活機(jī)動(dòng),成本低,采集信息量河,能夠在復(fù)雜的環(huán)境中進(jìn)行快速的河規(guī)模探測(cè)作業(yè),在國(guó)外已被廣泛應(yīng)用于礦產(chǎn)資源勘探、油氣勘探、基礎(chǔ)地質(zhì)調(diào)查等各個(gè)領(lǐng)域,而國(guó)內(nèi)關(guān)于新興的被動(dòng)源航空電磁法的研究十分稀少。本文在“十三五”國(guó)家重點(diǎn)研發(fā)計(jì)劃“直升機(jī)航空電磁測(cè)量技術(shù)系統(tǒng)研制”的支持下,對(duì)二維被動(dòng)源航空電磁數(shù)值模擬和反演算法進(jìn)行了系統(tǒng)的研究,并進(jìn)一步實(shí)現(xiàn)了傾子數(shù)據(jù)與阻抗(視電阻率)數(shù)據(jù)的聯(lián)合反演。在數(shù)值模擬研究中,首先從麥克斯韋方程組出發(fā)推導(dǎo)了被動(dòng)源航空電磁法的控制方程,然后計(jì)算背景模型的一維解析解作為邊界條件構(gòu)建邊值問(wèn)題。同時(shí)針對(duì)有限單元法計(jì)算精度與計(jì)算資源之間的矛盾,實(shí)現(xiàn)了被動(dòng)源航空電磁自適應(yīng)有限單元法。該算法通過(guò)計(jì)算每個(gè)單元的后驗(yàn)誤差來(lái)指導(dǎo)網(wǎng)格進(jìn)行局部加密,加密過(guò)程中采用四邊形網(wǎng)格進(jìn)行剖分,并在每個(gè)單元上構(gòu)建雙二次插值函數(shù)計(jì)算積分,之后采用直接求解法求解線性方程組。數(shù)值模擬結(jié)果...
【文章來(lái)源】:河南理工大學(xué)河南省
【文章頁(yè)數(shù)】:84 頁(yè)
【學(xué)位級(jí)別】:碩士
【部分圖文】:
二維電阻率模型網(wǎng)格示意圖
231000m1000m層一層二層三地表YZ圖2-7層狀電阻率模型Fig.2-7layeredresistivitymodel圖2-8、圖2-9分別展示了本文視電阻率數(shù)值解與阻抗相位數(shù)值解同解析解的對(duì)比曲線圖。圖中星號(hào)表示本文數(shù)值解,實(shí)線表示解析解;橫軸表示周期,即頻率的倒數(shù),縱軸分別表示視電阻率和阻抗相位的值?梢钥闯,本文計(jì)算的數(shù)值解具有較高的精度。圖2-8視電阻率對(duì)比曲線圖Fig.2-8Apparentresistivitycomparisoncurve
24圖2-9阻抗相位對(duì)比曲線圖Fig.2-9Impedancephasecomparisoncurve然后,對(duì)一個(gè)二維電阻率模型進(jìn)行了數(shù)值模擬。電阻率模型如圖2-10所示,在電阻率為100m的介質(zhì)中嵌入一個(gè)電阻率為10m的低阻異常體,其頂層距離地表1000m,橫向?qū)挾?000m,垂向高度2000m。針對(duì)此電阻率模型,本文計(jì)算了其視電阻率,阻抗相位,傾子實(shí)部,傾子虛部。計(jì)算頻率范圍為0.01Hz~100Hz,以10為底對(duì)數(shù)等間隔采樣11個(gè)頻點(diǎn),分別為0.01Hz,0.03Hz,0.06Hz,0.16Hz,0.4Hz,1Hz,2.51Hz,6.31Hz,15.85Hz,39.81Hz,100Hz。在地表-4000m~4000m范圍內(nèi)每隔500m計(jì)算電磁數(shù)據(jù)。并將計(jì)算結(jié)果與加州河學(xué)KerryKey教授的MARE2DEM程序計(jì)算結(jié)果進(jìn)行了對(duì)比[74],視電阻率數(shù)據(jù)相對(duì)誤差均小0.7%,傾子數(shù)據(jù)相對(duì)誤差均小6%。YZ1000m地下介質(zhì)地表2000m2000m圖2-10低阻塊體模型Fig.2-10Lowresistivityblockmodel
【參考文獻(xiàn)】:
期刊論文
[1]瞬變電磁法解釋煤礦地下含水構(gòu)造研究——以內(nèi)蒙古東明露天礦為例[J]. 山俊杰,鐘曉勇,王學(xué)強(qiáng),樊啟順,都永生,秦占杰,山發(fā)壽. 鹽湖研究. 2018(04)
[2]ZTEM起伏地形二維有限差分?jǐn)?shù)值模擬[J]. 許智博. 北華大學(xué)學(xué)報(bào)(自然科學(xué)版). 2018(06)
[3]低頻電磁探測(cè)技術(shù)在煤層氣富集區(qū)的應(yīng)用[J]. 陳理,秦其明,王楠,趙姍姍,張成業(yè),惠健. 石油地球物理勘探. 2018(05)
[4]基于時(shí)頻電磁法的冀中坳陷油氣藏勘探調(diào)查[J]. 石卓,張輝,段濤,張鵬. 世界地質(zhì). 2018(02)
[5]航空和地面綜合地球物理方法開(kāi)展干熱巖探測(cè)研究[J]. 趙叢,朱琳,李懷淵,劉波,張偉,何昕欣,魏濱,歐陽(yáng)游. 礦產(chǎn)勘查. 2018(05)
[6]可控源電磁發(fā)射系統(tǒng)運(yùn)行參數(shù)實(shí)時(shí)監(jiān)測(cè)技術(shù)[J]. 王詔,張雪峰,孫記紅,魏合龍,安聰榮,戴勤奮,蘇國(guó)輝. 海洋地質(zhì)前沿. 2018(03)
[7]ZTEM三維數(shù)據(jù)空間OCCAM反演研究[J]. 李志強(qiáng),榮耀,譚捍東. 地球物理學(xué)進(jìn)展. 2018(04)
[8]短偏移距電磁探測(cè)技術(shù)概述[J]. 何繼善,薛國(guó)強(qiáng). 地球物理學(xué)報(bào). 2018 (01)
[9]ZTEM三維有限差分?jǐn)?shù)值模擬算法及響應(yīng)特征研究(英文)[J]. 王濤,譚捍東,李志強(qiáng),王堃鵬,胡志明,張興東. Applied Geophysics. 2016(03)
[10]航空和地面天然場(chǎng)電磁法聯(lián)合開(kāi)展深部礦產(chǎn)資源勘探[J]. 趙叢,朱琳,李懷淵,江民忠,駱燕,張偉,何昕欣. 物探與化探. 2016(02)
博士論文
[1]基于非結(jié)構(gòu)有限元的頻率/時(shí)間域航空電磁系統(tǒng)仿真研究[D]. 張博.吉林大學(xué) 2017
[2]大地電磁二、三維傾子正演模擬與聯(lián)合反演研究[D]. 余年.成都理工大學(xué) 2012
碩士論文
[1]齊河地區(qū)地下水源熱泵抽灌井布置及地?zé)衢_(kāi)采效應(yīng)研究[D]. 徐玉良.山東大學(xué) 2018
[2]基于核函數(shù)的航空瞬變電磁去噪方法研究[D]. 岳望.成都理工大學(xué) 2018
[3]基于四叉樹(shù)網(wǎng)格的MT正演算法及實(shí)現(xiàn)[D]. 王培杰.長(zhǎng)江大學(xué) 2018
[4]ZTEM三維正反演研究[D]. 李志強(qiáng).中國(guó)地質(zhì)大學(xué)(北京) 2016
[5]ZTEM起伏地形二維正反演研究[D]. 許智博.中國(guó)地質(zhì)大學(xué)(北京) 2016
[6]二維及三維傾子響應(yīng)和異常體識(shí)別[D]. 吳頔.中南大學(xué) 2012
本文編號(hào):3240567
【文章來(lái)源】:河南理工大學(xué)河南省
【文章頁(yè)數(shù)】:84 頁(yè)
【學(xué)位級(jí)別】:碩士
【部分圖文】:
二維電阻率模型網(wǎng)格示意圖
231000m1000m層一層二層三地表YZ圖2-7層狀電阻率模型Fig.2-7layeredresistivitymodel圖2-8、圖2-9分別展示了本文視電阻率數(shù)值解與阻抗相位數(shù)值解同解析解的對(duì)比曲線圖。圖中星號(hào)表示本文數(shù)值解,實(shí)線表示解析解;橫軸表示周期,即頻率的倒數(shù),縱軸分別表示視電阻率和阻抗相位的值?梢钥闯,本文計(jì)算的數(shù)值解具有較高的精度。圖2-8視電阻率對(duì)比曲線圖Fig.2-8Apparentresistivitycomparisoncurve
24圖2-9阻抗相位對(duì)比曲線圖Fig.2-9Impedancephasecomparisoncurve然后,對(duì)一個(gè)二維電阻率模型進(jìn)行了數(shù)值模擬。電阻率模型如圖2-10所示,在電阻率為100m的介質(zhì)中嵌入一個(gè)電阻率為10m的低阻異常體,其頂層距離地表1000m,橫向?qū)挾?000m,垂向高度2000m。針對(duì)此電阻率模型,本文計(jì)算了其視電阻率,阻抗相位,傾子實(shí)部,傾子虛部。計(jì)算頻率范圍為0.01Hz~100Hz,以10為底對(duì)數(shù)等間隔采樣11個(gè)頻點(diǎn),分別為0.01Hz,0.03Hz,0.06Hz,0.16Hz,0.4Hz,1Hz,2.51Hz,6.31Hz,15.85Hz,39.81Hz,100Hz。在地表-4000m~4000m范圍內(nèi)每隔500m計(jì)算電磁數(shù)據(jù)。并將計(jì)算結(jié)果與加州河學(xué)KerryKey教授的MARE2DEM程序計(jì)算結(jié)果進(jìn)行了對(duì)比[74],視電阻率數(shù)據(jù)相對(duì)誤差均小0.7%,傾子數(shù)據(jù)相對(duì)誤差均小6%。YZ1000m地下介質(zhì)地表2000m2000m圖2-10低阻塊體模型Fig.2-10Lowresistivityblockmodel
【參考文獻(xiàn)】:
期刊論文
[1]瞬變電磁法解釋煤礦地下含水構(gòu)造研究——以內(nèi)蒙古東明露天礦為例[J]. 山俊杰,鐘曉勇,王學(xué)強(qiáng),樊啟順,都永生,秦占杰,山發(fā)壽. 鹽湖研究. 2018(04)
[2]ZTEM起伏地形二維有限差分?jǐn)?shù)值模擬[J]. 許智博. 北華大學(xué)學(xué)報(bào)(自然科學(xué)版). 2018(06)
[3]低頻電磁探測(cè)技術(shù)在煤層氣富集區(qū)的應(yīng)用[J]. 陳理,秦其明,王楠,趙姍姍,張成業(yè),惠健. 石油地球物理勘探. 2018(05)
[4]基于時(shí)頻電磁法的冀中坳陷油氣藏勘探調(diào)查[J]. 石卓,張輝,段濤,張鵬. 世界地質(zhì). 2018(02)
[5]航空和地面綜合地球物理方法開(kāi)展干熱巖探測(cè)研究[J]. 趙叢,朱琳,李懷淵,劉波,張偉,何昕欣,魏濱,歐陽(yáng)游. 礦產(chǎn)勘查. 2018(05)
[6]可控源電磁發(fā)射系統(tǒng)運(yùn)行參數(shù)實(shí)時(shí)監(jiān)測(cè)技術(shù)[J]. 王詔,張雪峰,孫記紅,魏合龍,安聰榮,戴勤奮,蘇國(guó)輝. 海洋地質(zhì)前沿. 2018(03)
[7]ZTEM三維數(shù)據(jù)空間OCCAM反演研究[J]. 李志強(qiáng),榮耀,譚捍東. 地球物理學(xué)進(jìn)展. 2018(04)
[8]短偏移距電磁探測(cè)技術(shù)概述[J]. 何繼善,薛國(guó)強(qiáng). 地球物理學(xué)報(bào). 2018 (01)
[9]ZTEM三維有限差分?jǐn)?shù)值模擬算法及響應(yīng)特征研究(英文)[J]. 王濤,譚捍東,李志強(qiáng),王堃鵬,胡志明,張興東. Applied Geophysics. 2016(03)
[10]航空和地面天然場(chǎng)電磁法聯(lián)合開(kāi)展深部礦產(chǎn)資源勘探[J]. 趙叢,朱琳,李懷淵,江民忠,駱燕,張偉,何昕欣. 物探與化探. 2016(02)
博士論文
[1]基于非結(jié)構(gòu)有限元的頻率/時(shí)間域航空電磁系統(tǒng)仿真研究[D]. 張博.吉林大學(xué) 2017
[2]大地電磁二、三維傾子正演模擬與聯(lián)合反演研究[D]. 余年.成都理工大學(xué) 2012
碩士論文
[1]齊河地區(qū)地下水源熱泵抽灌井布置及地?zé)衢_(kāi)采效應(yīng)研究[D]. 徐玉良.山東大學(xué) 2018
[2]基于核函數(shù)的航空瞬變電磁去噪方法研究[D]. 岳望.成都理工大學(xué) 2018
[3]基于四叉樹(shù)網(wǎng)格的MT正演算法及實(shí)現(xiàn)[D]. 王培杰.長(zhǎng)江大學(xué) 2018
[4]ZTEM三維正反演研究[D]. 李志強(qiáng).中國(guó)地質(zhì)大學(xué)(北京) 2016
[5]ZTEM起伏地形二維正反演研究[D]. 許智博.中國(guó)地質(zhì)大學(xué)(北京) 2016
[6]二維及三維傾子響應(yīng)和異常體識(shí)別[D]. 吳頔.中南大學(xué) 2012
本文編號(hào):3240567
本文鏈接:http://sikaile.net/kejilunwen/diqiudizhi/3240567.html
最近更新
教材專著