發(fā)布時間：2018-08-13 20:05

【摘要】：有限差分方法(Finite Difference Method,FDM)作為一種快速而精確的數(shù)值方法,在波動方程正演數(shù)值模擬領域得到了最為廣泛的應用,然而,當在不規(guī)則或者地表起伏不平的區(qū)域中模擬時,基于規(guī)則網(wǎng)格的有限差分法對波動方程離散時會產(chǎn)生階梯狀近似,影響模擬的精度,給有限差分方法的求解帶來困難。通過求解橢圓型偏微分方程生成的貼體網(wǎng)格為這種復雜區(qū)域中的有限差分波場模擬問題提供了一種有效的手段。借助貼體網(wǎng)格以及鏈式法則,將空間導數(shù)的計算從不規(guī)則的物理區(qū)域轉換到規(guī)則的計算域,在計算域中繼續(xù)使用有限差分方法。相比其他的不規(guī)則網(wǎng)格剖分方法,如插值法、映射法以及非結構網(wǎng)格,貼體網(wǎng)格在普遍性、精度和穩(wěn)定性方面都有優(yōu)勢。本文在復雜區(qū)域綜合應用了貼體網(wǎng)格、聲波方程正演模擬、逆時偏移(Reverse Time Migration,RTM)以及完全匹配層(Perfectly Matched Layer,PML)等技術。首先通過求解橢圓偏微分方程法生成貼體網(wǎng)格,根據(jù)鏈式法則代入順序不同,將二階位移形式的二維聲波方程改寫到曲線網(wǎng)格中,得到了兩種形式的貼體坐標二維聲波方程,第二種形式的波動方程相比第一種形式而言,形式對稱、更緊湊。為了壓制截斷邊界造成的人工邊界反射,推導了兩種形式的聲波方程對應的PML方程。在復雜起伏地表區(qū)域中,分部求和(Summation-by-Parts,SBP)有限差分方法的使用能夠確保曲線網(wǎng)格中非均勻介質數(shù)值模擬的穩(wěn)定性。SBP有限差分法特別適合求解含變系數(shù)的導數(shù)項。因此,對兩種形式的波動方程及PML針對性地采用了不同的差分方法:第一種形式中,導數(shù)項的系數(shù)可預先求得,采用中心差分方法進行數(shù)值離散;第二種形式中含變系數(shù),因此采用顯式的二階精度SBP有限差分方法進行離散,并采用Fourier譜分析方法討論了這種離散格式的穩(wěn)定性。與中心有限差分方法對比,分部求和有限差分方法的穩(wěn)定性更高。四階精度的有限差分方法在降低存儲需求和提高效率方面有很大的優(yōu)勢。在二階精度中心差分方法和SBP有限差分方法的基礎上,將兩種形式的貼體網(wǎng)格聲波方程和PML中空間導數(shù)項的離散均擴展到四階精度。采用Fourier譜分析方法研究了該四階SBP差分離散格式的穩(wěn)定性,得到了離散方程的穩(wěn)定性條件,并驗證了四階精度SBP有限差分方法比中心有限差分方法的穩(wěn)定性更高。同時,在同樣尺寸的模型中比較了二階和四階精度SBP有限差分方法的計算時間及存儲需求,驗證了四階方法的高精度、低頻散特性,顯示了四階方法在提高運算效率、降低存儲方面的能力。本論文將貼體網(wǎng)格的方法分別應用于地面地震和VSP(Vertical Seismic Profiling,垂直地震剖面)資料的逆時偏移研究,得到了曲線坐標系中各自對應的成像數(shù)據(jù)體,為解決實際地震采集中觀測系統(tǒng)不規(guī)則、地表起伏等情況下的數(shù)據(jù)處理難題提供一種有效的技術手段。
[Abstract]:As a fast and accurate numerical method, finite difference method (Finite Difference) has been widely used in the field of forward numerical simulation of wave equations. The finite difference method based on regular grid will produce ladder approximation when the wave equation is discretized, which will affect the accuracy of the simulation and make it difficult to solve the finite difference method. The body-fitted meshes generated by solving elliptic partial differential equations provide an effective means for the finite difference wave field simulation in this complex domain. With the help of body-fitted meshes and chain rules, the computation of spatial derivatives is transformed from irregular physical regions to regular computational domains, and the finite difference method is used in the computational domain. Compared with other irregular mesh generation methods, such as interpolation, mapping and unstructured meshes, body-fitted meshes have advantages in universality, accuracy and stability. In this paper, the techniques of body-fitted mesh, forward modeling of acoustic wave equation, inverse time migration (Reverse Time migration and perfectly matched layer (Perfectly Matched are applied in complex regions. Firstly, the body-fitted grid is generated by solving the elliptic partial differential equation. According to the different order of the chain rule, the two-dimensional acoustic wave equation in the form of second-order displacement is rewritten into the curvilinear grid, and two kinds of two-dimensional acoustic wave equations with body-fitted coordinates are obtained. The second form of wave equation is more symmetrical and compact than the first one. In order to suppress the artificial boundary reflection caused by truncated boundary, the PML equation corresponding to two kinds of acoustic wave equations is derived. In complex undulating surface regions, the use of the Summation-by-Partss-SBP finite difference method can ensure the stability of numerical simulation of non-uniform media in curved grids. The SBP finite difference method is particularly suitable for solving derivative terms with variable coefficients. In the first form, the coefficient of the derivative term can be obtained in advance, and the central difference method is used to discretize the two forms of wave equation and PML. Therefore, the explicit second-order precision SBP finite difference method is used for discretization, and the stability of the scheme is discussed by using the Fourier spectrum analysis method. Compared with the central finite difference method, the partial summation finite difference method is more stable. The fourth order finite-difference method has great advantages in reducing storage requirements and improving efficiency. Based on the second-order precision central difference method and the SBP finite difference method, the dispersion of the two kinds of body-fitted mesh acoustic equations and the spatial derivative term in PML are extended to the fourth order accuracy. The stability of the fourth order SBP difference discrete scheme is studied by using the Fourier spectrum analysis method. The stability conditions of the discrete equation are obtained, and the stability of the fourth order precision SBP finite difference method is proved to be higher than that of the central finite difference method. At the same time, the computation time and storage requirement of the second-order and fourth-order precision SBP finite difference method are compared in the model of the same size. The high precision and low frequency dispersion characteristics of the fourth-order method are verified, and it is shown that the fourth-order method is improving the computational efficiency. Reduce storage capabilities. In this paper, the method of body-fitted grid is applied to the inverse time migration of surface seismic and VSP (Vertical Seismic Profiling, vertical seismic profiles, and the corresponding imaging data volume in the curvilinear coordinate system is obtained. In order to solve the problem of data processing in the case of irregular observation system and surface undulation in actual seismic acquisition, this paper provides an effective technical means.
【學位授予單位】：中國石油大學(北京)
【學位級別】：博士
【學位授予年份】：2016
【分類號】：P631.4

【相似文獻】

相關期刊論文 前10條

1 薛桂霞,鄧世坤,劉秀娟;逆時偏移在探地雷達信號處理中的應用[J];煤田地質與勘探;2004年01期

2 楊勤勇;段心標;;逆時偏移技術發(fā)展現(xiàn)狀與趨勢[J];石油物探;2010年01期

3 張春燕;李振春;孫小東;;逆時偏移方法技術進展綜述[J];勘探地球物理進展;2010年05期

4 滕厚華;許衛(wèi)華;于曉莉;尚新民;賀紅蕾;;無反射遞推算法疊后逆時偏移的研究與應用[J];油氣地球物理;2011年04期

5 許璐;孟小紅;劉國峰;;逆時偏移去噪方法研究進展[J];地球物理學進展;2012年04期

6 王童奎;李瑩;郭愛華;石文武;楊曉利;孫明;潘歡;;逆時偏移技術在南堡1號構造中的應用研究[J];地球物理學進展;2012年06期

7 康智清;劉恩;司杰戈;陳康;;逆時偏移剖面假象產(chǎn)生機制及其消除方法[J];石油天然氣學報;2013年01期

8 胡明順;潘冬明;董守華;金紅娣;趙新;;基于散射理論的逆時偏移隨機邊界構建策略優(yōu)化分析[J];地球物理學進展;2013年04期

9 張叔倫;陳靜;張關泉;;逆時偏移的方向導數(shù)算法[J];石油地球物理勘探;1991年06期

10 Ruey-Chyuan Shih;A.R.Levander;陳炳文;;層剝離逆時偏移[J];國外油氣勘探;1995年04期

相關會議論文 前10條

1 杜啟振;秦童;;逆時偏移最佳匹配層邊界實現(xiàn)[A];中國地球物理·2009[C];2009年

2 陳康;;逆時偏移剖面的波數(shù)域分析及高階拉普拉斯濾波算法的實現(xiàn)[A];“地球·資源”全國博士生學術論壇會議論文摘要集[C];2011年

3 龍云;韓立國;韓利;鄧武斌;;單程波逆時偏移[A];中國地球物理學會第二十七屆年會論文集[C];2011年

4 李慶洋;李振春;黃建平;;基于雙變算法的精細逆時偏移成像[A];中國地球物理2013——第二十二專題論文集[C];2013年

5 曹曉莉;黃建平;李振春;;最小二乘逆時偏移高精度成像方法研究[A];中國地球物理2013——第二十二專題論文集[C];2013年

6 許璐;孟小紅;劉國峰;;成像后壓制逆時偏移噪聲[A];中國地球物理學會第二十七屆年會論文集[C];2011年

7 薛霆哠;張書杰;;逆時偏移中坡印廷向量濾波的改進方法[A];中國地球物理2013——第二十二專題論文集[C];2013年

8 楊仁虎;毛偉健;凌云;;逆時偏移儲層成像波形校正研究[A];中國地球物理2013——第二十專題論文集[C];2013年

9 曾同生;李振春;;基于疊前逆時偏移的速度分析方法研究[A];中國地球物理學會第二十七屆年會論文集[C];2011年

10 胡昊;劉伊克;王一博;;地震疊前逆時偏移不同邊界處理方式研究[A];中國地球物理學會第二十七屆年會論文集[C];2011年

相關重要報紙文章 前2條

1 王躍;物探院逆時偏移成像技術躋身國際前列[N];中國石化報;2013年

2 本報記者 劉洪宇;石油HPC走向“平衡計算”[N];中國計算機報;2011年

相關博士學位論文 前10條

1 王穎;貼體網(wǎng)格剖分的地震資料逆時偏移方法研究[D];中國石油大學(北京);2016年

2 李靜爽;基于近似解析離散化算子的逆時偏移方法及其應用[D];清華大學;2014年

3 程飛;多井井間三維波場數(shù)值模擬與逆時偏移成像研究[D];中國地質大學;2015年

4 馮蘭蘭;保辛幾何結構不變的逆時偏移方法研究[D];清華大學;2015年

5 楊佳佳;多分量地震波逆時偏移的關鍵技術研究[D];中國海洋大學;2015年

6 高新成;基于云計算的逆時偏移數(shù)據(jù)處理方法研究[D];東北石油大學;2016年

7 Naveed Akram;基于彈性波模擬的逆時偏移以及地震數(shù)據(jù)解釋[D];中國科學技術大學;2017年

8 王娟;TTI介質準縱波逆時偏移方法研究[D];中國石油大學(華東);2015年

9 趙旭;時域聲波方程逆時偏移成像及其角道集計算[D];吉林大學;2015年

10 蔡志成;基于Q補償?shù)酿椊橘|逆時偏移成像研究[D];中國地質大學;2017年

相關碩士學位論文 前10條

1 陳洪杰;基于聲波方程的數(shù)值模擬與逆時偏移方法研究[D];吉林大學;2010年

2 魏巍;基于疊前逆時偏移方法的圓管損傷識別[D];哈爾濱工業(yè)大學;2015年

3 張敬秋;交錯網(wǎng)格有限差分法地震資料逆時偏移[D];成都理工大學;2015年

4 劉詩竹;VSP數(shù)據(jù)成像方法研究[D];東北石油大學;2015年

5 柯璇;復雜構造疊前逆時偏移成像方法研究[D];東北石油大學;2015年

6 侯維麗;逆時偏移成像[D];東北石油大學;2015年

7 郭雪豹;VSP地震資料逆時偏移方法研究[D];東北石油大學;2015年

8 張千祥;純P波各向異性介質波場數(shù)值模擬及逆時偏移研究[D];吉林大學;2016年

9 于江龍;聲波逆時偏移及噪聲壓制[D];吉林大學;2016年

10 方修政;各向異性介質逆時偏移方法研究[D];東北石油大學;2016年

，

本文編號：2182038