本文選題：節(jié)點(diǎn)密度 + 重力正演　； 參考：《中國地質(zhì)大學(xué)(北京)》2016年博士論文

【摘要】：正演是反演的基礎(chǔ),重力正演計(jì)算效率和質(zhì)量直接決定了反演方法的實(shí)用效果�；谥绷⒘骟w模型剖分的重力異常正演中存在嚴(yán)重的冗余計(jì)算,降低了模型正演效率,為此本研究引入網(wǎng)格節(jié)點(diǎn)密度的概念,對(duì)原有重力正演公式進(jìn)行簡(jiǎn)化,得到基于網(wǎng)格節(jié)點(diǎn)密度的重力正演公式。不同計(jì)算環(huán)境下的計(jì)算速度對(duì)比實(shí)驗(yàn)結(jié)果表明,基于網(wǎng)格節(jié)點(diǎn)密度的重力正演方法有效消除了冗余計(jì)算,提高了正演計(jì)算效率,加速比可以達(dá)到4到8倍。通過對(duì)邊界節(jié)點(diǎn)的特殊處理,基于網(wǎng)格節(jié)點(diǎn)密度的重力正演方法能夠正確計(jì)算無限延伸,或具有真實(shí)密度模型的重力異常。此外,網(wǎng)格節(jié)點(diǎn)密度定義為塊體密度的差分形式,具有良好的稀疏性,這在一定程度上有利于提高計(jì)算效率,節(jié)省物理內(nèi)存消耗。利用上述正演計(jì)算方案,本文正演計(jì)算了Marmousi2密度模型、BP2004密度模型和SEG/EAGE三維鹽丘密度模型的重力及重力梯度異常。在密度模型結(jié)構(gòu)簡(jiǎn)單情況下,網(wǎng)格節(jié)點(diǎn)密度模型較塊體密度模型具有更好的稀疏性。為此,本文在基于網(wǎng)格節(jié)點(diǎn)密度模型正演研究基礎(chǔ)上,進(jìn)一步開展了基于節(jié)點(diǎn)密度模型的稀疏反演方法研究。在反演過程中,采用柯西范數(shù)模型約束來得到稀疏節(jié)點(diǎn)密度模型解,通過引入深度加權(quán)約束改善模型深度分辨率。模型試驗(yàn)表明,該反演方法能夠真實(shí)有效恢復(fù)密度異常體,異常體邊界清晰,位置準(zhǔn)確,有效改善了最小模型解、光滑模型解、聚焦約束模型解等反演結(jié)果分辨率低、異常體邊界模糊的問題。此外,為改善重力物性反演深度分辨率,有效利用地震資料提供的地層深度信息,本文探索研究了地震層位約束的網(wǎng)格節(jié)點(diǎn)密度反演方法。通過在網(wǎng)格節(jié)點(diǎn)密度模型中引入地震層位加權(quán),有效提高了反演結(jié)果的深度分辨率。在反演密度模型解釋方面,提出通過構(gòu)建真實(shí)密度模型,真實(shí)反映研究區(qū)地下密度結(jié)構(gòu)信息,避免因異常分離導(dǎo)致的剩余密度異常解釋上的誤差。利用基于網(wǎng)格節(jié)點(diǎn)密度的重力正反演方法,結(jié)合曾母研究區(qū)區(qū)域地質(zhì)、巖石物性等資料,開展研究區(qū)三維真實(shí)密度模型構(gòu)建工作。利用已有地震層位信息構(gòu)建背景密度模型,利用基于網(wǎng)格節(jié)點(diǎn)密度模型的重力正演剝離方法,獲取研究區(qū)剩余重力異常,之后通過基于節(jié)點(diǎn)密度的稀疏反演方法獲得研究區(qū)剩余密度異常,從而成功構(gòu)建研究區(qū)真密度結(jié)構(gòu)。
[Abstract]:Forward modeling is the basis of inversion. The efficiency and quality of gravity forward modeling directly determine the practical effect of inversion method. The gravity anomaly forward modeling based on vertical hexahedron model has serious redundancy calculation, which reduces the forward modeling efficiency. In this paper, the concept of grid node density is introduced to simplify the original gravity forward modeling formula. The forward gravity formula based on grid node density is obtained. The experimental results show that the gravity forward modeling method based on grid node density can effectively eliminate redundant computation and improve the efficiency of forward calculation. The speedup ratio can reach 4 to 8 times. Through special processing of boundary nodes, gravity forward modeling based on grid node density can correctly calculate infinite extension or gravity anomaly with real density model. In addition, the grid node density is defined as the difference form of block density, which has good sparsity, which is helpful to improve the computational efficiency and save physical memory consumption to some extent. In this paper, the gravity and gravity gradient anomalies of Marmousi2 density model, BP2004 density model and SEG / EAGE three-dimensional salt dome density model are calculated by using the forward modeling scheme mentioned above. When the density model is simple, the mesh node density model is more sparse than the block density model. Therefore, based on the forward modeling of grid node density model, the sparse inversion method based on node density model is further studied in this paper. In the inversion process, Cauchy norm model constraints are used to obtain the sparse node density model solution, and the depth resolution of the model is improved by introducing depth-weighted constraints. The model tests show that the inversion method can restore the density anomaly body effectively, the boundary of the anomaly body is clear, the location is accurate, and the resolution of the inversion results such as minimum model solution, smooth model solution and focus constrained model solution are improved effectively, and the resolution of inversion results is low. The problem of fuzzy boundary of abnormal body. In addition, in order to improve the depth resolution of the inversion of gravity physical properties and make effective use of the depth information provided by seismic data, the inversion method of grid node density with seismic horizon constraints is studied in this paper. The depth resolution of inversion results is improved effectively by introducing seismic horizon weighting into the grid node density model. In the interpretation of inversion density model, it is proposed to construct a real density model to truly reflect the underground density structure information in the study area, so as to avoid the error in the interpretation of residual density anomaly caused by abnormal separation. Using the method of gravity forward and inverse modeling based on grid node density and combining the data of regional geology and petrophysical properties of Zengmu study area, the 3D real density model of the study area is constructed. The background density model is constructed by using the existing seismic horizon information and the gravity forward stripping method based on the grid node density model is used to obtain the residual gravity anomalies in the study area. Then the residual density anomaly is obtained by sparse inversion method based on node density, and the true density structure of the studied area is constructed successfully.
【學(xué)位授予單位】：中國地質(zhì)大學(xué)(北京)
【學(xué)位級(jí)別】：博士
【學(xué)位授予年份】：2016
【分類號(hào)】：P631.1

【參考文獻(xiàn)】

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

1 楊建輝;王亮;張家德;楊勝發(fā);;黔東南重磁異常與金異常特征及找礦意義[J];物探與化探;2013年05期

2 彭淼;譚捍東;姜枚;錢輝;譚嘉言;;基于交叉梯度耦合的大地電磁與地震走時(shí)資料三維聯(lián)合反演[J];地球物理學(xué)報(bào);2013年08期

3 劉銀萍;王祝文;杜曉娟;劉菁華;許家姝;;基于Extrapolation Tikhonov正則化算法的重力數(shù)據(jù)三維約束反演[J];地球物理學(xué)報(bào);2013年05期

4 陳少華;朱自強(qiáng);魯光銀;曹書錦;;重力梯度張量的預(yù)條件共軛梯度法反演[J];中南大學(xué)學(xué)報(bào)(自然科學(xué)版);2013年02期

5 熊莉娟;李三忠;索艷慧;劉鑫;余珊;程世秀;薛友辰;安慧婷;戴黎明;馬云;王霄飛;;南海南部新生代控盆斷裂特征及盆地群成因[J];海洋地質(zhì)與第四紀(jì)地質(zhì);2012年06期

6 陳召曦;孟小紅;郭良輝;劉國峰;;基于GPU并行的重力、重力梯度三維正演快速計(jì)算及反演策略[J];地球物理學(xué)報(bào);2012年12期

7 郭良輝;孟小紅;石磊;陳召曦;;優(yōu)化濾波方法及其在中國大陸布格重力異常數(shù)據(jù)處理中的應(yīng)用[J];地球物理學(xué)報(bào);2012年12期

8 祁光;呂慶田;嚴(yán)加永;吳明安;劉彥;;先驗(yàn)地質(zhì)信息約束下的三維重磁反演建模研究——以安徽泥河鐵礦為例[J];地球物理學(xué)報(bào);2012年12期

9 丘學(xué)林;趙明輝;徐輝龍;李家彪;阮愛國;郝天珧;游慶瑜;;南海深地震探測(cè)的重要科學(xué)進(jìn)程:回顧和展望[J];熱帶海洋學(xué)報(bào);2012年03期

10 陳召曦;孟小紅;郭良輝;;重磁數(shù)據(jù)三維物性反演方法進(jìn)展[J];地球物理學(xué)進(jìn)展;2012年02期

，

本文編號(hào)：2057606