本文關(guān)鍵詞：求解波動(dòng)方程的低數(shù)值頻散方法及數(shù)值模擬，由筆耕文化傳播整理發(fā)布。

【摘要】：數(shù)值頻散和計(jì)算效率問題是求解波動(dòng)方程的關(guān)鍵和難點(diǎn)。針對(duì)這一問題,本論文應(yīng)用帕德逼近和近似解析離散化的思想,首先發(fā)展了求解波動(dòng)方程的帕德逼近方法(簡稱PAM)。該方法采用帕德逼近方法進(jìn)行時(shí)間離散,時(shí)間微分算子是一個(gè)有理分式,導(dǎo)致離散格式是一種隱式格式。為避免求解系數(shù)為塊三對(duì)角矩陣的大型線性代數(shù)方程組,我們采用差分算子顯式化的方法,將隱式格式變成顯式格式。在空間離散方面,采用近似解析離散化方法,使用函數(shù)及其空間梯度的線性組合共同逼近空間高階導(dǎo)數(shù),使得離散格式包含更多的波場(chǎng)信息。這些信息更加有利于提高地震波反演和地震偏移的精度和成像質(zhì)量。同時(shí),這種近似解析離散化算子半徑短,有更好的緊致性,遵循地震介質(zhì)的物理特性。將帕德逼近方法和近似解析離散化方法結(jié)合,可以達(dá)到時(shí)間四階精度,空間八階精度。另外,對(duì)數(shù)值離散格式中將會(huì)出現(xiàn)4階和5階微分算子,本文采用了算子分裂的方式,可以有效降低算法中所含微分算子的階數(shù),從而減少計(jì)算量。本文首先給出帕德逼近方法數(shù)值求解地震波方程的離散格式,對(duì)一維和二維離散格式穩(wěn)進(jìn)行了穩(wěn)定性、誤差和數(shù)值頻散等分析,同時(shí)與其它具有相同精度的數(shù)值方法,如8階Lax Wendoroff Correction(簡稱LWC)方法、8階Staggered Grid(簡稱SG)方法進(jìn)行計(jì)算效率和數(shù)值模擬效果的比較,說明帕德逼近方法的低數(shù)值頻散和高計(jì)算效率的特性。然后論文將帕德逼近方法應(yīng)用于求解三維地震波方程,在空間離散采取了4階精度和8階精度兩種不同的離散格式。在數(shù)值離散格式的理論分析和數(shù)值模擬方面與具有同階精度的LWC方法和SG方法進(jìn)行了相應(yīng)的比較研究,進(jìn)一步證明了帕德逼近方法模擬三維地震波傳播是可行的,且具有低數(shù)值頻散性和高計(jì)算效率等優(yōu)勢(shì)。最后論文還將龍格庫塔方法求解三維地震波方程推廣到空間八階精度,相比原來的四階方法,在提高數(shù)值精度的同時(shí),也降低了數(shù)值頻散。另一方面論文將帕德逼近方法推廣到數(shù)值求解雙相Biot方程。給出了雙相介質(zhì)模型中帕德逼近方法的數(shù)值格式,并進(jìn)行了含流體多孔隙介質(zhì)中地震波傳播模擬。從數(shù)值模擬雙相介質(zhì)模型的波場(chǎng)快照中可見清晰的快P波,慢P波和S波。同時(shí)與LWC方法等傳統(tǒng)方法進(jìn)行了比較,PAM方法壓制數(shù)值頻散效果更好。
【關(guān)鍵詞】：正演方法 帕德逼近 數(shù)值頻散 PML 數(shù)值模擬
【學(xué)位授予單位】：清華大學(xué)
【學(xué)位級(jí)別】：博士
【學(xué)位授予年份】：2015
【分類號(hào)】：O241.82
【目錄】：

• 摘要3-4
• Abstract4-9
• 第1章 引言9-21
• 1.1 問題背景9-10
• 1.2 基本彈性波方程10-13
• 1.2.1 彈性介質(zhì)中的平衡方程式10-11
• 1.2.2 均勻各向同性介質(zhì)中的彈性波方程11-12
• 1.2.3 均勻橫向各向同性介質(zhì)中的彈性波方程12-13
• 1.3 常用的數(shù)值模擬方法13-18
• 1.3.1 Lax-Wendroff修正格式14-16
• 1.3.2 位移-應(yīng)力交錯(cuò)網(wǎng)格格式16-17
• 1.3.3 近似解析離散化方法17-18
• 1.4 本文的主要內(nèi)容和安排18-21
• 第2章 帕德逼近方法和近似解析離散化算子21-32
• 2.1 帕德逼近方法21-23
• 2.2 近似解析離散化算子23-31
• 2.3 本章小結(jié)31-32
• 第3章 求解波動(dòng)方程的低數(shù)值頻散帕德逼近方法32-79
• 3.1 直接方法32-37
• 3.2 顯式PAM方法37-39
• 3.3 PAM穩(wěn)定性分析39-43
• 3.3.1 一維穩(wěn)定性條件39-41
• 3.3.2 二維穩(wěn)定性條件41-43
• 3.4 誤差分析43-48
• 3.4.1 理論誤差分析43
• 3.4.2 數(shù)值誤差43-48
• 3.5 數(shù)值頻散分析48-58
• 3.5.1 一維PAM方法的數(shù)值頻散48-51
• 3.5.2 二維PAM方法的數(shù)值頻散51-58
• 3.6 計(jì)算效率分析與數(shù)值結(jié)果比較58-71
• 3.6.1 計(jì)算效率分析58-63
• 3.6.2 數(shù)值算例63-71
• 3.7 匹配PAM的PML吸收邊界條件71-78
• 3.7.1 分裂波場(chǎng)PML的構(gòu)造74-76
• 3.7.2 分裂波場(chǎng)PML與PAM結(jié)合76-78
• 3.8 本章小結(jié)78-79
• 第4章 三維帕德逼近方法及地震波模擬79-107
• 4.1 三維帕德逼近空間四階方法81-98
• 4.1.1 四階三維帕德逼近方法的穩(wěn)定性條件及數(shù)值頻散82-87
• 4.1.2 誤差分析87-89
• 4.1.3 計(jì)算效率比較89-92
• 4.1.4 數(shù)值模擬算例92-98
• 4.2 8階空間精度的三維帕德逼近方法98-106
• 4.2.1 穩(wěn)定性條件及頻散分析98-101
• 4.2.2 計(jì)算效率比較101-104
• 4.2.3 數(shù)值模擬比較104-106
• 4.3 本章小結(jié)106-107
• 第5章 求解波動(dòng)方程的三維龍格庫塔方法107-118
• 5.1 八階龍格庫塔方法107-109
• 5.2 穩(wěn)定性條件109
• 5.3 計(jì)算效率109-112
• 5.4 數(shù)值模擬112-117
• 5.5 本章小結(jié)117-118
• 第6章 求解雙相Biot方程的帕德逼近方法118-131
• 6.1 雙相介質(zhì)中的彈性波方程118-119
• 6.2 雙相Biot方程的PAM方法119-122
• 6.3 數(shù)值模擬122-130
• 6.3.1 與其他數(shù)值方法對(duì)比122-128
• 6.3.2 雙層模型128-130
• 6.4 本章小結(jié)130-131
• 第7章 總結(jié)與展望131-133
• 7.1 總結(jié)131-132
• 7.2 后續(xù)研究工作與展望132-133
• 致謝133-134
• 參考文獻(xiàn)134-142
• 附錄PAM方法的增長矩陣142-149
• 個(gè)人簡歷、在學(xué)期間發(fā)表的學(xué)術(shù)論文與研究成果149

【參考文獻(xiàn)】

中國期刊全文數(shù)據(jù)庫 前10條

1 鐘擁;;地震屬性在古城地區(qū)碳酸鹽巖儲(chǔ)層預(yù)測(cè)中的應(yīng)用[J];石油工業(yè)計(jì)算機(jī)應(yīng)用;2014年01期

2 汪文帥;李小凡;魯明文;張美根;;基于多辛結(jié)構(gòu)譜元法的保結(jié)構(gòu)地震波場(chǎng)模擬[J];地球物理學(xué)報(bào);2012年10期

3 王美霞;楊頂輝;宋國杰;;二維SH波方程的半解析解及其數(shù)值模擬[J];地球物理學(xué)報(bào);2012年03期

4 宋國杰;楊頂輝;童平;練永生;;求解3D彈性波方程的并行WNAD方法及其TI介質(zhì)中的波場(chǎng)模擬[J];地球物理學(xué)報(bào);2012年02期

5 陳可洋;;完全匹配層吸收邊界條件研究[J];石油物探;2010年05期

6 肖開宇;胡祥云;;正演模擬技術(shù)在地震解釋中的應(yīng)用[J];工程地球物理學(xué)報(bào);2009年04期

7 王磊;楊頂輝;鄧小英;;非均勻介質(zhì)中地震波應(yīng)力場(chǎng)的WNAD方法及其數(shù)值模擬[J];地球物理學(xué)報(bào);2009年06期

8 陳可洋;;聲波完全匹配層吸收邊界條件的改進(jìn)算法[J];石油物探;2009年01期

9 LIU Lucy;;Seismic tomography of the Moon[J];Chinese Science Bulletin;2008年24期

10 馮英杰;楊長春;吳萍;;地震波有限差分模擬綜述[J];地球物理學(xué)進(jìn)展;2007年02期

中國博士學(xué)位論文全文數(shù)據(jù)庫 前6條

1 李靜爽;基于近似解析離散化算子的逆時(shí)偏移方法及其應(yīng)用[D];清華大學(xué);2014年

2 楊磊;含多相流彈性及粘彈性介質(zhì)中波頻散與衰減分析及其應(yīng)用[D];清華大學(xué);2014年

3 馬嘯;波動(dòng)方程保辛近似解析離散化算法研究[D];清華大學(xué);2013年

4 童平;地震層析成像方法及其應(yīng)用研究[D];清華大學(xué);2012年

5 宋國杰;并行WNAD算法及其波場(chǎng)模擬[D];清華大學(xué);2011年

6 陳山;求解波動(dòng)方程的龍格—庫塔型方法及其地震波傳播模擬[D];清華大學(xué);2010年

中國碩士學(xué)位論文全文數(shù)據(jù)庫 前5條

1 章曉;求解波動(dòng)方程的近似解析指數(shù)時(shí)間差分方法及其數(shù)值模擬[D];清華大學(xué);2013年

2 王美霞;雙相及黏彈性介質(zhì)中波傳播方程的保辛算法及其波場(chǎng)模擬[D];清華大學(xué);2012年

3 王念;隱式Runge-Kutta方法的顯式化算法及波場(chǎng)模擬[D];清華大學(xué);2009年

4 王磊;求解波動(dòng)方程的多步法及其數(shù)值模擬[D];清華大學(xué);2009年

5 宋國杰;三維彈性波方程的改進(jìn)近似解析離散化方法及波場(chǎng)模擬[D];清華大學(xué);2008年

  本文關(guān)鍵詞：求解波動(dòng)方程的低數(shù)值頻散方法及數(shù)值模擬，，由筆耕文化傳播整理發(fā)布。

本文編號(hào)：257003