本文選題：有限差分　切入點：高階交錯網格　出處：《中國科學技術大學》2017年碩士論文　論文類型：學位論文

【摘要】：有限差分算法是地震波場數值模擬中的重要工具,其具有實現簡單、計算效率高、對內存需求小等特點,在天然地震與勘探地球物理領域中得到了廣泛的應用。隨著地球科學的進一步發(fā)展,研究模型變得更加復雜與龐大,因此需要更高效的正演算法。使用高階交錯網格格式可以在使用較粗網格的情況下,有效的壓制色散誤差,進而降低計算內存需求節(jié)省計算資源。但是,較大的網格步長會造成有限差分網格對復雜內部界面的表征誤差增大。實際的地震模型中存在著大量的彈性系數間斷界面,因為網格離散介質的局限性,很難在算法中準確地滿足內部界面邊界條件。尤其在高階格式的情況下,由于使用了同低階格式相比更大的空間步長,因此這種問題會更加嚴重。本文使用高階交錯網格有限差分算法分別討論了處理介質間斷面的不同方法,并分析了界面誤差的數值特征,為高精度的地震波數值模擬提供了計算參數的選擇依據。本文第一部分,大致介紹了有限差分正演算法的基本情況,并對高階交錯網格與有限差分界面誤差研究現狀做了簡短的總結。第二部分,主要描述波動方程有限差分正演的基本實現,包括網格布局、高階格式、震源與邊界條件處理。第三部分引入內部界面邊界條件,顯式地分析了有限差分的界面數值誤差,同時說明了目前常用的四種介質參數取值方法:直接取值法、格點平均法、體積分平均法、正交各向異性介質法。在第四部分,通過數值實驗,分別展示了不同介質處理方法的震相對比,與界面誤差受不同因素的變化特征。最后,通過總結分析上述的誤差特征,給出高階交錯網格有限差分在實際正演計算中的網格參數選擇建議。
[Abstract]:The finite difference algorithm is an important tool in the numerical simulation of seismic wave field. It has the advantages of simple implementation, high computational efficiency and low memory requirement. It has been widely used in the field of natural earthquake and exploration geophysics. With the further development of geoscience, the research model becomes more complex and huge. Therefore, a more efficient forward algorithm is needed. Using higher-order interleaved grid schemes can effectively suppress dispersion errors in the case of coarse grids, thereby reducing the computational memory requirements and saving computing resources. The larger step size of the mesh will cause the finite difference mesh to increase the representation error of the complex internal interface. There are a large number of elastic coefficient discontinuous interfaces in the actual seismic model because of the limitation of the grid discrete medium. It is difficult to satisfy the boundary condition of internal interface accurately in the algorithm, especially in the case of high order scheme, because the space step size is larger than that of lower order scheme. Therefore, this problem will be more serious. In this paper, the different methods of dealing with dielectric discontinuity are discussed by using the high order staggered grid finite difference algorithm, and the numerical characteristics of the interface errors are analyzed. This paper provides a basis for the selection of calculation parameters for high precision seismic wave numerical simulation. In the first part of this paper, the basic situation of the finite difference forward algorithm is introduced. In the second part, the basic realization of finite difference forward modeling of wave equation is described, including grid layout, high order scheme, and so on. In the third part, the interface numerical error of finite difference is analyzed explicitly by introducing the boundary condition of internal interface. At the same time, four commonly used methods for determining the parameters of medium are described: direct method, lattice average method, and so on. The volume average method and the orthotropic medium method. In part 4th, the seismic characteristics of different medium processing methods are shown by numerical experiments, and the interface errors are changed by different factors. By summing up and analyzing the above error characteristics, the paper gives some suggestions on the selection of grid parameters in the actual forward calculation of high order staggered grid finite difference.
【學位授予單位】：中國科學技術大學
【學位級別】：碩士
【學位授予年份】：2017
【分類號】：P631.4;P315.8

【相似文獻】

相關期刊論文 前10條

1 袁修貴,黃誠;二階導數的小波多分辨表示與有限差分表示比較[J];數學理論與應用;2004年03期

2 魏文斌;袁俊泉;劉偉明;劉清國;;基于正交規(guī)范化的小波的有限差分表示(英文)[J];數學雜志;2008年03期

3 陳楚琳,陳華杰;連續(xù)梁的有限差分解[J];力學與實踐;1989年06期

4 周家紀;劉竟波;王鐸;;有限差分的快速算法[J];物探化探計算技術;1993年02期

5 張育英,張維全;二階偏微分方程的有限差分解法及其在氣動力計算中的應用[J];甘肅工業(yè)大學學報;1992年02期

6 S.A.Magnier;;采用最小網格的有限差分法[J];石油物探譯叢;1995年03期

7 C·L·艾肯;D·A·哈斯廷斯;張桂青;;激發(fā)極化的數字模擬[J];物化探電子計算技術;1981年01期

8 張德榮,袁陵濤,魯百年;一類具有磁場效應的Schr銉dinger方程組的有限差分解[J];高校應用數學學報A輯(中文版);1987年02期

9 苑志江;金良安;田恒斗;盧yN斌;;海洋拖曳系統(tǒng)運動方程的變步長有限差分數值仿真[J];指揮控制與仿真;2013年01期

10 王候庭;姚國治;劉書成;;厚板焊接溫度場的隱顯交變有限差分計算方法[J];大連鐵道學院學報;1984年01期

相關會議論文 前10條

1 陳勇;韓波;陳小宏;;有限差分注入法在時間推移地震中的應用[A];中國地球物理學會年刊2002——中國地球物理學會第十八屆年會論文集[C];2002年

2 席豐;;爆炸與火載荷作用下鋼梁響應的有限差分計算模型[A];中國力學學會學術大會'2009論文摘要集[C];2009年

3 唐暉;李小軍;;顯式積分格式對波動傳播規(guī)律的影響[A];中國地震學會第11次學術大會論文摘要集[C];2006年

4 林庚浩;馬天寶;寧建國;;三維前處理有限差分網格剖分技術的研究[A];中國計算力學大會'2010（CCCM2010）暨第八屆南方計算力學學術會議（SCCM8）論文集[C];2010年

5 徐朝繁;張先康;段永紅;;三維射線有限差分混合地震模型及應用研究[A];2001年中國地球物理學會年刊——中國地球物理學會第十七屆年會論文集[C];2001年

6 梁文全;楊長春;王彥飛;劉紅偉;;用于聲波方程數值模擬的時間-空間域有限差分系數確定新方法[A];中國科學院地質與地球物理研究所2013年度（第13屆）學術論文匯編——油氣資源研究室[C];2014年

7 孫林潔;印興耀;吳國忱;;基于PM_L邊界條件的高倍可變網格有限差分數值模擬方法[A];中國地球物理2010——中國地球物理學會第二十六屆年會、中國地震學會第十三次學術大會論文集[C];2010年

8 田小波;吳慶舉;曾融生;;有限差分合成接收函數[A];中國地球物理學會年刊2002——中國地球物理學會第十八屆年會論文集[C];2002年

9 張金海;王衛(wèi)民;符力耘;姚振興;;傅立葉有限差分算子和廣義屏算子的對比研究[A];中國地球物理學會第22屆年會論文集[C];2006年

10 蘇麗娟;;雙邊空間分數階對流—擴散方程的一種有限差分解法[A];中國力學學會學術大會'2009論文摘要集[C];2009年

相關博士學位論文 前2條

1 呂桂霞;拋物方程有限差分并行算法理論[D];吉林大學;2004年

2 胡平;探地雷達數值模擬技術的應用研究[D];中國地質大學（北京）;2005年

相關碩士學位論文 前10條

1 劉埕伊;基于樹圖法和有限差分法的多叉樹美式期權定價模型研究[D];哈爾濱工業(yè)大學;2016年

2 朱強;內部復雜速度界面處理對有限差分地震波模擬準確性影響研究[D];中國科學技術大學;2017年

3 李芳;2.5-D聲波方程高階有限差分數值分析[D];吉林大學;2007年

4 王麗麗;結構網格虛擬單元浸入邊界有限差分基本無振蕩格式[D];南京航空航天大學;2012年

5 孫旭;對比源方法在時間推移地震中的應用[D];哈爾濱工業(yè)大學;2013年

6 王紅;變網格步長聲波方程有限差分數值模擬[D];中南大學;2013年

7 楊瑩;二維地震波場有限差分法數值模擬研究[D];中國地質大學（北京）;2009年

8 赫英紅;帶有阻尼項的波動方程有限差分解的長時間行為[D];黑龍江大學;2012年

9 田多勤;一維耗散的Schr(o|¨)dinger方程有限差分解的長時間行為[D];黑龍江大學;2012年

10 周忠雙;基于MPI的一類橢圓型方程有限差分區(qū)域分解算法的并行實現[D];山東大學;2013年

，

本文編號：1587184