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