發(fā)布時間：2021-04-09 11:39

　　利用衛(wèi)星測高數(shù)據(jù)解算重力場是獲取大范圍海域重力異常的主要途徑.然而,由于涉及到全球范圍的多重積分,傳統(tǒng)串行計算方法效率低下.針對上述問題,本文研究提出一種高效的海域重力異常計算方法.該快速方法利用逆Stokes公式進(jìn)行重力場解算,同時針對其中不同點(diǎn)位多重循環(huán)計算的效率低下問題提出了一種基于多核CPU技術(shù)的并行化改進(jìn)方法,該方法通過將程序合理地劃分為CPU不同線程,在運(yùn)算過程中合理調(diào)配計算機(jī)資源以完成海域重力異常的高效并行計算.文章利用EGM2008模型提供的大地水準(zhǔn)面數(shù)據(jù)進(jìn)行了快速算法測試,并將解算出的重力異常與船測重力異常進(jìn)行了比對,表明該快速算法能夠準(zhǔn)確高效地求取出研究區(qū)重力異常,計算效率相比較于傳統(tǒng)串行算法提升30%以上.本文提出的并行化方法對解算大范圍海域重力異常和精化地球重力場模型具有積極意義. 

【文章來源】：地球物理學(xué)進(jìn)展. 2020,35(06)北大核心CSCD

【文章頁數(shù)】：6 頁

【部分圖文】：

并行計算模式偽代碼

運(yùn)算過程中發(fā)現(xiàn),在計算點(diǎn)及其附近區(qū)域到計算點(diǎn)的理論距離值極小或等于零,使得逆Stokes公式中的積分核函數(shù)奇異(王瑞和李厚樸, 2010).考慮到該公式是對全球范圍進(jìn)行積分,因此該區(qū)域范圍對計算點(diǎn)的重力異常值的貢獻(xiàn)可忽略不計,即在計算點(diǎn)處結(jié)果直接取零值.同時由于積分區(qū)域較廣,數(shù)據(jù)較多,積分時需要對每個點(diǎn)進(jìn)行一次對全球的積分,這就導(dǎo)致積分運(yùn)算過程緩慢,耗費(fèi)大量時間.針對運(yùn)算效率低下問題,本文研究提出一種高效的基于多核CPU技術(shù)的并行化改進(jìn)方法,通過合理分配任務(wù),減少核間通信達(dá)到降低串行計算部分的比例和減少交互開銷的目的,充分發(fā)揮多核性能,提高運(yùn)算速度.圖3 3.75′×3.75′研究區(qū)海域大地水準(zhǔn)面

3.75′×3.75′研究區(qū)海域大地水準(zhǔn)面

【參考文獻(xiàn)】：
期刊論文
[1]MATLAB并行計算解決方案[J]. 姚尚鋒,劉長江,唐正華,戴迪.  計算機(jī)時代. 2016(09)
[2]南海重力異常特征及其顯著的構(gòu)造意義[J]. 李淑玲,孟小紅,郭良輝,姚長利.  現(xiàn)代地質(zhì). 2012(06)
[3]多核CPU和GPU加速分子動力學(xué)模擬[J]. 林江宏,林錦賢,呂暾.  計算機(jī)應(yīng)用. 2011(03)
[4]基于逆Stokes公式的測高重力反演中央?yún)^(qū)效應(yīng)計算[J]. 王瑞,李厚樸.  武漢大學(xué)學(xué)報(信息科學(xué)版). 2010(04)
[5]衛(wèi)星測高在確定地球重力場中的應(yīng)用[J]. 翟國君.  海洋測繪. 2002(02)
[6]中國近海衛(wèi)星測高數(shù)據(jù)的重力異常反演——I.Stokes公式逆運(yùn)算+2DFFT法[J]. 王海瑛,陸洋,王廣運(yùn).  高技術(shù)通訊. 2000(10)
[7]南海及鄰區(qū)重力場特征與地殼構(gòu)造區(qū)劃[J]. 張訓(xùn)華.  海洋地質(zhì)與第四紀(jì)地質(zhì). 1998(03)
[8]用T/P測高數(shù)據(jù)反演海域重力異常[J]. 歐陽永忠,翟國君.  解放軍測繪學(xué)院學(xué)報. 1998(01)
[9]衛(wèi)星測高在物理大地測量應(yīng)用中的若干問題[J]. 李建成,寧津生,晁定波.  武漢測繪科技大學(xué)學(xué)報. 1996(01)

碩士論文
[1]衛(wèi)星測高反演海洋重力應(yīng)用研究[D]. 馬宏達(dá).中國地質(zhì)大學(xué)(北京) 2017
[2]多核CPU上數(shù)據(jù)庫散列連接算法的研究與實(shí)現(xiàn)[D]. 王梓.西安電子科技大學(xué) 2014



本文編號：3127555