發(fā)布時(shí)間：2018-07-31 09:04

【摘要】：近年來(lái),隨著瞬變電磁技術(shù)的不斷發(fā)展,瞬變電磁儀器的不斷進(jìn)步,我們可以采集到越來(lái)越精確可靠的瞬變電磁野外數(shù)據(jù)。然而我們隊(duì)瞬變電磁法的研究仍不深入,沒(méi)能很好地發(fā)揮作用。目前瞬變電磁反演,主要以一維為主,二維三維反演仍不成熟。為了更加深入地研究瞬變電磁法,我們對(duì)時(shí)間域有限差分三維瞬變電磁正演進(jìn)行深入研究,但是常規(guī)地CPU上的計(jì)算耗時(shí)太長(zhǎng),模型太大的話計(jì)算時(shí)間長(zhǎng)到無(wú)法忍受。隨著近些年GPU并行計(jì)算的深入發(fā)展,形成了包括CUDA、OpenGL、OpenAcc 等一大批 GPU 并行技術(shù)。但是由于 CUDA、OpenGL等傳統(tǒng)的并行技術(shù)要與硬件結(jié)構(gòu)高度結(jié)合,因而程序編寫難度較大,并且對(duì)于已經(jīng)寫好的程序移植性較差,所以技術(shù)普及應(yīng)用受到限制。而OpenAcc技術(shù)解決了傳統(tǒng)并行技術(shù)在這些方面難題,并取得了較大的進(jìn)步。本文利用新型的GPU并行編程技術(shù)OpenAcc對(duì)時(shí)域有限差分三維瞬變電磁正演算法進(jìn)行改進(jìn),極大縮短了三維正演的計(jì)算時(shí)間,增大了計(jì)算效率,并針對(duì)其進(jìn)行了如下研究:首先,應(yīng)用Yee (1966)提出的交錯(cuò)網(wǎng)格方案和Wang Hohmann (1993)提出的DuFort-Frankel方案對(duì)均勻網(wǎng)格剖分下的麥克斯韋方程組的FDTD離散形式進(jìn)行推倒,并采用軟源直接注入的方式施加激勵(lì)源,推導(dǎo)得到有源媒質(zhì)中的MaxWell有限差分方程組。對(duì)GPU并行技術(shù)的發(fā)展與優(yōu)勢(shì)進(jìn)行簡(jiǎn)要介紹,重點(diǎn)介紹了 OpenAcc并行技術(shù)的執(zhí)行模式與優(yōu)勢(shì),應(yīng)用實(shí)際實(shí)例證明OpenAcc并行技術(shù)相較于傳統(tǒng)CUDA技術(shù)的優(yōu)勢(shì)以及選取其作為本文研究工具的原因;其次,對(duì)FDTD三維正演算法進(jìn)行分析,找出程序最耗時(shí)的部分,應(yīng)用OpenAcc并行技術(shù)對(duì)其進(jìn)行并行化處理。并應(yīng)用經(jīng)典的均勻半空間模型、層狀模型以及三維異常體模型對(duì)并行程序的準(zhǔn)確性進(jìn)行驗(yàn)證;最后,對(duì)并行程序的加速性能進(jìn)行評(píng)價(jià),通過(guò)試驗(yàn),并行程序的最大加速比可以達(dá)到20x左右。通過(guò)大量的模型試驗(yàn)得到在本文所用的計(jì)算環(huán)境下并行程序的計(jì)算上限,分析網(wǎng)格剖分方式對(duì)并行程序加速性能的影響,尋找最適合的模型剖分方法。最后在更加實(shí)際的有覆蓋的地質(zhì)模型、復(fù)雜的地質(zhì)模型進(jìn)行計(jì)算,驗(yàn)證了并行程序在絕大多數(shù)地質(zhì)模型下準(zhǔn)確性極高,加速性性能也達(dá)到滿意效果。
[Abstract]:In recent years, with the continuous development of transient electromagnetic technology and the continuous progress of transient electromagnetic instruments, we can collect more and more accurate and reliable transient electromagnetic field data. However, the study of transient electromagnetic method in our team is still not thorough and has not worked well. At present, transient electromagnetic inversion is mainly one-dimensional, two-dimensional inversion is still immature. In order to study the transient electromagnetic method more deeply, we study the 3-D transient electromagnetic forward modeling in finite difference time domain, but the computation time on the conventional CPU is too long, and if the model is too large, the calculation time is too long to bear. With the development of GPU parallel computing in recent years, a large number of GPU parallel technologies have been formed, such as CUDAN OpenGL OpenAcc and so on. However, because the traditional parallel technology such as CUDA-OpenGL is highly integrated with the hardware structure, it is difficult to program writing, and the portability of the written program is poor, so the popularization and application of the technology is limited. OpenAcc technology solves these problems in traditional parallel technology and has made great progress. In this paper, a new GPU parallel programming technique, OpenAcc, is used to improve the 3-D transient electromagnetic forward algorithm of finite-difference time-domain, which greatly shortens the computing time and increases the computational efficiency. Using the staggered grid scheme proposed by Yee (1966) and the DuFort-Frankel scheme proposed by Wang Hohmann (1993), the FDTD discrete form of Maxwell equations under uniform mesh division is pushed down, and the excitation source is applied by direct injection of soft source. The MaxWell finite difference equations in active media are derived. The development and advantages of GPU parallel technology are briefly introduced, and the execution mode and advantages of OpenAcc parallel technology are emphatically introduced. A practical example is used to prove the advantage of OpenAcc parallel technology compared with traditional CUDA technology and the reason why it is chosen as the research tool in this paper. Secondly, the 3D forward algorithm of FDTD is analyzed to find out the most time-consuming part of the program. OpenAcc parallel technology is applied to parallelize it. The accuracy of the parallel program is verified by the classical uniform half-space model, layered model and 3D abnormal volume model. Finally, the acceleration performance of the parallel program is evaluated and tested. The maximum speedup of parallel programs can reach about 20x. Through a large number of model tests, the upper limit of parallel programs is obtained in the computing environment used in this paper. The effect of grid generation on the acceleration performance of parallel programs is analyzed, and the most suitable model generation method is found. Finally, in the more practical overlay geological model and the complex geological model, it is verified that the parallel program is highly accurate under most geological models, and the speedup performance is satisfactory.
【學(xué)位授予單位】：中國(guó)地質(zhì)大學(xué)(北京)
【學(xué)位級(jí)別】：碩士
【學(xué)位授予年份】：2017
【分類號(hào)】：P631.325

【參考文獻(xiàn)】

相關(guān)期刊論文 前7條

1 李建慧;朱自強(qiáng);魯光銀;曾思紅;;回線源瞬變電磁法的三維正演研究[J];地球物理學(xué)進(jìn)展;2013年02期

2 藍(lán)朝暉;胡希偉;江中和;;Interaction of Electromagnetic Waves with Two-Dimensional Metal Covered with Radar Absorbing Material and Plasma[J];Plasma Science and Technology;2008年06期

3 李進(jìn)璽;程引會(huì);周輝;吳偉;;用傳輸線和時(shí)域有限差分法計(jì)算電纜X射線響應(yīng)[J];強(qiáng)激光與粒子束;2007年12期

4 胡淑欣;朱守正;;含有左手材料的T結(jié)波導(dǎo)分析[J];材料導(dǎo)報(bào);2007年S2期

5 許愛(ài)國(guó);;用FDTD分析微帶結(jié)構(gòu)過(guò)程的優(yōu)化[J];微波學(xué)報(bào);2007年S1期

6 薛國(guó)強(qiáng);李貅;底青云;;瞬變電磁法理論與應(yīng)用研究進(jìn)展[J];地球物理學(xué)進(jìn)展;2007年04期

7 閆玉波,葛寧,鄭美艷,葛德彪,田春明;網(wǎng)絡(luò)并行FDTD方法分析電大目標(biāo)電磁散射[J];電子學(xué)報(bào);2003年06期

相關(guān)博士學(xué)位論文 前3條

1 毛玉蓉;時(shí)間域電磁響應(yīng)三維正演計(jì)算及GPU實(shí)現(xiàn)[D];長(zhǎng)江大學(xué);2014年

2 孫懷鳳;隧道含水構(gòu)造三維瞬變電磁場(chǎng)響應(yīng)特征及突水災(zāi)害源預(yù)報(bào)研究[D];山東大學(xué);2013年

3 關(guān)珊珊;基于GPU的三維有限差分直升機(jī)瞬變電磁響應(yīng)并行計(jì)算[D];吉林大學(xué);2012年

相關(guān)碩士學(xué)位論文 前10條

1 劉二萌;瞬變電磁時(shí)域有限差分正演模擬[D];長(zhǎng)安大學(xué);2015年

2 吳學(xué)銀;瞬變電磁法三維有限差分正演研究[D];中國(guó)地質(zhì)大學(xué)(北京);2015年

3 牛延強(qiáng);聲波正演模擬的OpenACC實(shí)現(xiàn)研究[D];燕山大學(xué);2015年

4 劉堯;基于有限差分的瞬變電磁三維數(shù)值模擬[D];中國(guó)地質(zhì)大學(xué)(北京);2015年

5 李雁冰;面向異構(gòu)多核處理器并行化編譯中的數(shù)據(jù)優(yōu)化[D];解放軍信息工程大學(xué);2014年

6 陳呈;面向MIC平臺(tái)的OpenACC實(shí)現(xiàn)與優(yōu)化關(guān)鍵技術(shù)研究[D];國(guó)防科學(xué)技術(shù)大學(xué);2013年

7 郭春輝;基于GPU的電力系統(tǒng)并行計(jì)算的研究[D];山東大學(xué);2013年

8 趙云威;矩形回線源瞬變電磁法三維有限差分正演模擬[D];中南大學(xué);2012年

9 賴俊良;基于CUDA的超聲彈性成像系統(tǒng)的研究與實(shí)現(xiàn)[D];電子科技大學(xué);2012年

10 吳長(zhǎng)莉;基于MPI和OpenMP的三維FDTD并行算法的研究[D];華中科技大學(xué);2009年

，

本文編號(hào)：2155095