隨著地表速度復雜區(qū)域勘探力度的加大,地震勘探作業(yè)區(qū)域的速度結構越來越復雜。特別地,在我國西部探區(qū)和山前帶區(qū)域,由于構造運動和風化作用的影響,地表構造特別復雜,低速異常普遍存在,由地表起伏引起的地震資料信噪比低、靜校正問題復雜以及觀測系統(tǒng)布置困難等問題正影響著復雜區(qū)域地震工作的開展。為了研究這一系列實際問題,必須對復雜地表下地震波的傳播規(guī)律有清醒的認識。在地震波傳播的研究中,地球物理界普遍認可地震波傳播數(shù)值模擬,也就是通常所說的正演模擬。正演模擬的算法很多且各有優(yōu)勢,本次選擇有限差分法和譜元法進行研究。有限差分法使用差分格式替代空間偏導項,該方法編程難度小、計算效率高以及對復雜地質(zhì)模型適用的優(yōu)點使其被廣泛采用。而譜元法作為一種高階有限元法,兼具有限元法剖分模型的貼體性和譜方法的計算效率,逐漸被更多人關注。本文先總結了復雜地表條件下地震波數(shù)值模擬方法,再從有限差分法和譜元法入手完成以下研究工作:1)通過應力鏡像法實施自由邊界條件后,使用交錯網(wǎng)格差分完成了起伏地表下彈性波方程有限差分法數(shù)值模擬;2)在彈性波方程弱形式的基礎上,使用Galerkin譜展開法、Lagrange插值基函數(shù)以及Gauss-Lobatto-Legendre積分方法,完成了起伏地表下彈性波方程譜元法數(shù)值模擬;3)對比有限差分法和譜元法的精度、計算效率和計算機內(nèi)存占用情況,并使用MPI+CUDA的異構并行技術完成譜元法的加速。4)使用譜元法對地表起伏以及地表速度復雜模型進行彈性波傳播數(shù)值模擬,并對波場特征做了詳細分析。通過研究發(fā)現(xiàn):1)譜元法的貼體網(wǎng)格劃分能克服矩形網(wǎng)格剖分帶來的角點數(shù)值散射,提高了數(shù)值模擬精度;2)有限差分法的計算時間和計算機內(nèi)存占用均小于譜元法,但是并不能改善其由于矩形網(wǎng)格劃分所帶來的數(shù)值散射,使用MPI+CUDA的異構并行技術能極大的提高計算效率,使較大模型的正演變得現(xiàn)實;3)在起伏地表模型下彈性波的傳播模式發(fā)生了很大的變化,散射波和多次波同時影響著記錄,這也是起伏地表下地震資料信噪比低的原因。
【學位單位】：西南石油大學
【學位級別】：碩士
【學位年份】：2018
【中圖分類】：P631.4
【部分圖文】：

圖２－３?ＰＭＬ吸收邊界模型試算快照ｔ＝６６０ｍｓ：??

圖２－５水平自由地表（地表過ｃｒＹＹ、ｃｒｚｚ方向）示意圖??２０??

圖２－６起伏自由地表（地表過（Ｊ？．、方向）示意圖，黑色實線為為地表，實心方形代表切??應力，實心圓形代表正應力，垂直速度平行于切應力，水平速度垂直于正應力??
【參考文獻】

相關期刊論文 前10條

1 劉有山;滕吉文;徐濤;劉少林;司薌;馬學英;;三角網(wǎng)格譜元法地震波場數(shù)值模擬[J];地球物理學進展;2014年04期

2 劉有山;滕吉文;劉少林;徐濤;;稀疏存儲的顯式有限元三角網(wǎng)格地震波數(shù)值模擬及其PML吸收邊界條件[J];地球物理學報;2013年09期

3 黃超;董良國;;可變網(wǎng)格與局部時間步長的交錯網(wǎng)格高階差分彈性波模擬[J];地球物理學報;2009年11期

4 黃超;董良國;;可變網(wǎng)格與局部時間步長的高階差分地震波數(shù)值模擬[J];地球物理學報;2009年01期

5 楊旭明;裴正林;;復雜近地表彈性波波場特征研究[J];石油地球物理勘探;2007年06期

6 董良國;郭曉玲;吳曉豐;馬在田;;起伏地表彈性波傳播有限差分法數(shù)值模擬[J];天然氣工業(yè);2007年10期

7 王祥春;劉學偉;;起伏地表二維聲波方程地震波場模擬與分析[J];石油地球物理勘探;2007年03期

8 王秀明;G.SERIANI;林偉軍;;利用譜元法計算彈性波場的若干理論問題[J];中國科學(G輯:物理學 力學 天文學);2007年01期

9 董良國;復雜地表條件下地震波傳播數(shù)值模擬[J];勘探地球物理進展;2005年03期

10 裴正林;任意起伏地表彈性波方程交錯網(wǎng)格高階有限差分法數(shù)值模擬[J];石油地球物理勘探;2004年06期

本文編號：2891029