步長(zhǎng)自適應(yīng)的有限差分復(fù)雜地表波場(chǎng)數(shù)值模擬
本文選題:自適應(yīng)網(wǎng)格 切入點(diǎn):波動(dòng)方程 出處:《地球物理學(xué)進(jìn)展》2017年03期 論文類(lèi)型:期刊論文
【摘要】:復(fù)雜地表?xiàng)l件下的有限差分地震波場(chǎng)的數(shù)值模擬,由于受到低速層和地表起伏的限制,模型速度分布范圍變大,一般使用精細(xì)的差分網(wǎng)格來(lái)抑制頻散,提高模擬分辨率,但精細(xì)網(wǎng)格會(huì)顯著增加計(jì)算成本.為了能有效地解決這一問(wèn)題,本文提出一種步長(zhǎng)自適應(yīng)有限差分波動(dòng)方程數(shù)值模擬方法.(1)新方法根據(jù)模型中的介質(zhì)速度分布,對(duì)不同的速度區(qū)域采用與該速度匹配的空間步長(zhǎng),實(shí)現(xiàn)對(duì)模型空間網(wǎng)格的步長(zhǎng)自適應(yīng)精細(xì)劃分.對(duì)于速度分布范圍大的復(fù)雜地表模型,新方法不僅能夠極大地減少模型的網(wǎng)格節(jié)點(diǎn)數(shù),同時(shí)又能提高波場(chǎng)的時(shí)間采樣步長(zhǎng),減少時(shí)間采樣數(shù),提高計(jì)算效率.(2)推導(dǎo)了不同步長(zhǎng)邊界網(wǎng)格節(jié)點(diǎn)Laplace算子的二階有限差分表達(dá)式,避免了在這些結(jié)點(diǎn)進(jìn)行插值計(jì)算產(chǎn)生的假擾動(dòng)和數(shù)值不穩(wěn)定問(wèn)題.(3)為了降低有限差分產(chǎn)生的數(shù)值頻散,本文在常規(guī)的差分方程中增加了一頻散校正項(xiàng),能有效地衰減了高波數(shù)成分,抑制了數(shù)值頻散.對(duì)復(fù)雜近地表的波場(chǎng)數(shù)值模擬結(jié)果表明,本文提出的步長(zhǎng)自適應(yīng)新方法能夠有效減少網(wǎng)格節(jié)點(diǎn)數(shù)和時(shí)間采樣數(shù),極大地提高計(jì)算效率,計(jì)算量比常規(guī)粗網(wǎng)格增加一些,但效果能夠達(dá)到了常規(guī)精細(xì)網(wǎng)格的模擬結(jié)果.
[Abstract]:The numerical simulation of finite difference seismic wave field under complex surface condition is limited by the low velocity layer and surface fluctuation, and the velocity distribution range of the model becomes larger. The fine difference grid is generally used to suppress dispersion and improve the simulation resolution. In order to solve this problem effectively, a step size adaptive finite difference wave equation numerical simulation method, I. e., a new method based on the velocity distribution of the medium in the model, is presented in this paper. For different velocity regions, the space step size matched with the velocity is adopted to realize the adaptive fine partition of the model space mesh. For the complex surface model with large velocity distribution, The new method can not only greatly reduce the number of mesh nodes in the model, but also increase the time sampling step size of the wave field and reduce the number of time samples. The second order finite-difference expression of the Laplace operator of grid nodes with different step sizes is derived. In order to reduce the numerical dispersion generated by finite difference, a dispersion correction term is added to the ordinary difference equation. The numerical simulation results show that the proposed step size adaptive method can effectively reduce the number of mesh nodes and the number of time samples, and the numerical simulation results show that the proposed method can effectively reduce the number of mesh nodes and the number of time samples, and the numerical simulation results show that the proposed step size adaptive method can effectively reduce the number of mesh nodes and the number of time samples. The computational efficiency is greatly improved, and the computational complexity is a little more than that of the conventional coarse mesh, but the simulation results of the conventional fine mesh can be achieved.
【作者單位】: 西安理工大學(xué)自動(dòng)化與信息工程學(xué)院;
【基金】:國(guó)家自然科學(xué)基金資助重大項(xiàng)目課題(41390454)資助
【分類(lèi)號(hào)】:P631.4
【相似文獻(xiàn)】
相關(guān)期刊論文 前10條
1 魏文斌;袁俊泉;劉偉明;劉清國(guó);;基于正交規(guī)范化的小波的有限差分表示(英文)[J];數(shù)學(xué)雜志;2008年03期
2 陳楚琳,陳華杰;連續(xù)梁的有限差分解[J];力學(xué)與實(shí)踐;1989年06期
3 周家紀(jì);劉竟波;王鐸;;有限差分的快速算法[J];物探化探計(jì)算技術(shù);1993年02期
4 張育英,張維全;二階偏微分方程的有限差分解法及其在氣動(dòng)力計(jì)算中的應(yīng)用[J];甘肅工業(yè)大學(xué)學(xué)報(bào);1992年02期
5 S.A.Magnier;;采用最小網(wǎng)格的有限差分法[J];石油物探譯叢;1995年03期
6 C·L·艾肯;D·A·哈斯廷斯;張桂青;;激發(fā)極化的數(shù)字模擬[J];物化探電子計(jì)算技術(shù);1981年01期
7 張德榮,袁陵濤,魯百年;一類(lèi)具有磁場(chǎng)效應(yīng)的Schr銉dinger方程組的有限差分解[J];高校應(yīng)用數(shù)學(xué)學(xué)報(bào)A輯(中文版);1987年02期
8 苑志江;金良安;田恒斗;盧yN斌;;海洋拖曳系統(tǒng)運(yùn)動(dòng)方程的變步長(zhǎng)有限差分?jǐn)?shù)值仿真[J];指揮控制與仿真;2013年01期
9 王候庭;姚國(guó)治;劉書(shū)成;;厚板焊接溫度場(chǎng)的隱顯交變有限差分計(jì)算方法[J];大連鐵道學(xué)院學(xué)報(bào);1984年01期
10 馬紅磊,王燕昌,李進(jìn);變步長(zhǎng)的有限差分線法[J];寧夏大學(xué)學(xué)報(bào)(自然科學(xué)版);2000年02期
相關(guān)會(huì)議論文 前10條
1 陳勇;韓波;陳小宏;;有限差分注入法在時(shí)間推移地震中的應(yīng)用[A];中國(guó)地球物理學(xué)會(huì)年刊2002——中國(guó)地球物理學(xué)會(huì)第十八屆年會(huì)論文集[C];2002年
2 席豐;;爆炸與火載荷作用下鋼梁響應(yīng)的有限差分計(jì)算模型[A];中國(guó)力學(xué)學(xué)會(huì)學(xué)術(shù)大會(huì)'2009論文摘要集[C];2009年
3 唐暉;李小軍;;顯式積分格式對(duì)波動(dòng)傳播規(guī)律的影響[A];中國(guó)地震學(xué)會(huì)第11次學(xué)術(shù)大會(huì)論文摘要集[C];2006年
4 林庚浩;馬天寶;寧建國(guó);;三維前處理有限差分網(wǎng)格剖分技術(shù)的研究[A];中國(guó)計(jì)算力學(xué)大會(huì)'2010(CCCM2010)暨第八屆南方計(jì)算力學(xué)學(xué)術(shù)會(huì)議(SCCM8)論文集[C];2010年
5 徐朝繁;張先康;段永紅;;三維射線有限差分混合地震模型及應(yīng)用研究[A];2001年中國(guó)地球物理學(xué)會(huì)年刊——中國(guó)地球物理學(xué)會(huì)第十七屆年會(huì)論文集[C];2001年
6 梁文全;楊長(zhǎng)春;王彥飛;劉紅偉;;用于聲波方程數(shù)值模擬的時(shí)間-空間域有限差分系數(shù)確定新方法[A];中國(guó)科學(xué)院地質(zhì)與地球物理研究所2013年度(第13屆)學(xué)術(shù)論文匯編——油氣資源研究室[C];2014年
7 孫林潔;印興耀;吳國(guó)忱;;基于PM_L邊界條件的高倍可變網(wǎng)格有限差分?jǐn)?shù)值模擬方法[A];中國(guó)地球物理2010——中國(guó)地球物理學(xué)會(huì)第二十六屆年會(huì)、中國(guó)地震學(xué)會(huì)第十三次學(xué)術(shù)大會(huì)論文集[C];2010年
8 田小波;吳慶舉;曾融生;;有限差分合成接收函數(shù)[A];中國(guó)地球物理學(xué)會(huì)年刊2002——中國(guó)地球物理學(xué)會(huì)第十八屆年會(huì)論文集[C];2002年
9 張金海;王衛(wèi)民;符力耘;姚振興;;傅立葉有限差分算子和廣義屏算子的對(duì)比研究[A];中國(guó)地球物理學(xué)會(huì)第22屆年會(huì)論文集[C];2006年
10 蘇麗娟;;雙邊空間分?jǐn)?shù)階對(duì)流—擴(kuò)散方程的一種有限差分解法[A];中國(guó)力學(xué)學(xué)會(huì)學(xué)術(shù)大會(huì)'2009論文摘要集[C];2009年
相關(guān)博士學(xué)位論文 前2條
1 呂桂霞;拋物方程有限差分并行算法理論[D];吉林大學(xué);2004年
2 胡平;探地雷達(dá)數(shù)值模擬技術(shù)的應(yīng)用研究[D];中國(guó)地質(zhì)大學(xué)(北京);2005年
相關(guān)碩士學(xué)位論文 前10條
1 劉埕伊;基于樹(shù)圖法和有限差分法的多叉樹(shù)美式期權(quán)定價(jià)模型研究[D];哈爾濱工業(yè)大學(xué);2016年
2 李芳;2.5-D聲波方程高階有限差分?jǐn)?shù)值分析[D];吉林大學(xué);2007年
3 王麗麗;結(jié)構(gòu)網(wǎng)格虛擬單元浸入邊界有限差分基本無(wú)振蕩格式[D];南京航空航天大學(xué);2012年
4 孫旭;對(duì)比源方法在時(shí)間推移地震中的應(yīng)用[D];哈爾濱工業(yè)大學(xué);2013年
5 王紅;變網(wǎng)格步長(zhǎng)聲波方程有限差分?jǐn)?shù)值模擬[D];中南大學(xué);2013年
6 楊瑩;二維地震波場(chǎng)有限差分法數(shù)值模擬研究[D];中國(guó)地質(zhì)大學(xué)(北京);2009年
7 赫英紅;帶有阻尼項(xiàng)的波動(dòng)方程有限差分解的長(zhǎng)時(shí)間行為[D];黑龍江大學(xué);2012年
8 田多勤;一維耗散的Schr(o|¨)dinger方程有限差分解的長(zhǎng)時(shí)間行為[D];黑龍江大學(xué);2012年
9 周忠雙;基于MPI的一類(lèi)橢圓型方程有限差分區(qū)域分解算法的并行實(shí)現(xiàn)[D];山東大學(xué);2013年
10 趙良;拋物方程的有限差分區(qū)域分解算法[D];山東大學(xué);2008年
,本文編號(hào):1559489
本文鏈接:http://sikaile.net/kejilunwen/kuangye/1559489.html