發(fā)布時(shí)間：2018-01-29 00:55

  本文關(guān)鍵詞： 曲線(xiàn)坐標(biāo) 貼體網(wǎng)格 有限差分 各向異性介質(zhì) 固體-液體耦合介質(zhì)　出處：《中國(guó)科學(xué)技術(shù)大學(xué)》2017年博士論文　論文類(lèi)型：學(xué)位論文

【摘要】：在地震學(xué)和與它的相關(guān)領(lǐng)域中,有限差分方法是一種重要的模擬地震波在復(fù)雜介質(zhì)中傳播的數(shù)值計(jì)算方法;它具有簡(jiǎn)單直接計(jì)算效率高等優(yōu)點(diǎn),已經(jīng)被廣泛地應(yīng)用于地震勘探、強(qiáng)地面運(yùn)動(dòng)預(yù)測(cè)、全波形反演、地震各向異性介質(zhì)的研究、合成理論地震圖等方面。傳統(tǒng)的有限差分方法采用矩形網(wǎng)格離散計(jì)算區(qū)域,并且使用臺(tái)階近似的方法來(lái)描述復(fù)雜的地形;這種對(duì)地形的近似既降低了差分方法實(shí)施邊界條件的精度又會(huì)讓它在計(jì)算中產(chǎn)生虛假的散射波,從而不能正確地計(jì)算地形對(duì)地震波傳播的影響;雖然對(duì)它改進(jìn)的方法在不斷地被提出,但是這些方法的適用性還有待深入地研究。本文所研究的曲線(xiàn)網(wǎng)格有限差分方法采用貼體網(wǎng)格描述復(fù)雜地形;這種離散地形的方法不僅可以讓有限差分方法避免因?yàn)榕_(tái)階近似所導(dǎo)致的虛假散射波,而且還可以方便它準(zhǔn)確地實(shí)施自由表面邊界條件,從而正確地計(jì)算出復(fù)雜地形對(duì)地震波傳播的影響。地震波傳播過(guò)程的數(shù)值模擬是研究地震波在各向異性介質(zhì)中傳播規(guī)律的重要手段,對(duì)地震勘探和天然地震的研究有著重要的意義。傳統(tǒng)的有限差分方法采用交錯(cuò)網(wǎng)格中心差分格式;這種方法雖然可以高效率地計(jì)算出地震波在各向同性介質(zhì)中的傳播過(guò)程,但是當(dāng)任意各向異性介質(zhì)出現(xiàn)在計(jì)算區(qū)域時(shí),它需要對(duì)各向異性介質(zhì)中的地震波場(chǎng)分量進(jìn)行插值從而導(dǎo)致很大的計(jì)算誤差。本文所研究的曲線(xiàn)網(wǎng)格有限差分方法采用同位網(wǎng)格DRP/opt MacCormack差分格式;這種方法不僅可以避免傳統(tǒng)有限差分方法對(duì)各向異性中的地震波場(chǎng)分量的插值所導(dǎo)致的近似,而且還可以結(jié)合貼體網(wǎng)格技術(shù)和牽引力鏡像技術(shù)處理復(fù)雜地形,從而正確地計(jì)算出地震波在任意各向異性介質(zhì)中的傳播過(guò)程以及地形對(duì)各向異性介質(zhì)中地震波傳播的影響。準(zhǔn)確地模擬地震波通過(guò)固體-液體界面的物理過(guò)程對(duì)于研究海洋中的T震相和斯通利波、海底地形對(duì)地震波的散射作用以及解釋海底地震儀記錄的振動(dòng)數(shù)據(jù)等具有重要的意義。傳統(tǒng)的有限差分方法對(duì)復(fù)雜的幾何形狀的固體-液體界面的離散和固體-液體邊界條件的實(shí)施都做了相當(dāng)程度的近似,不能準(zhǔn)確地模擬起伏的固體-液體界面對(duì)地震波傳播的影響。本文所研究的曲線(xiàn)網(wǎng)格有限差分方法可以準(zhǔn)確地描述固體-液體界面的復(fù)雜幾何形狀,更加準(zhǔn)確直接地實(shí)施邊界條件,為研究復(fù)雜幾何形狀的固體-液體界面對(duì)地震波傳播的影響提供了一種更加可靠的計(jì)算工具。本文將集成了貼體網(wǎng)格、同位網(wǎng)格DRP/opt MacCormack差分格式和牽引力鏡像等技術(shù)的曲線(xiàn)網(wǎng)格有限差分方法分別推廣到了二維和三維情況下含有復(fù)雜地形的各向異性介質(zhì)和含有復(fù)雜固體-液體界面形狀的固體-液體耦合介質(zhì)的地震波傳播計(jì)算問(wèn)題中,并且將曲線(xiàn)網(wǎng)格有限差分方法計(jì)算的結(jié)果分別同廣義反射/透射系數(shù)方法和譜元方法計(jì)算的結(jié)果進(jìn)行了對(duì)比驗(yàn)證。驗(yàn)證結(jié)果表明,本文所研究的曲線(xiàn)網(wǎng)格有限差分方法可以正確并且有效地應(yīng)用于上述地震波傳播問(wèn)題。
[Abstract]:In the field of seismology and with it, the finite difference method is a calculation method for simulation of seismic wave propagation in complex media is important in numerical; it has simple calculation and high efficiency, has been widely used in seismic exploration and prediction of strong ground motion, full waveform inversion, seismic anisotropy the synthetic seismogram. The traditional finite difference method using rectangular grid discrete computational domain, and use the method of step approximation to describe the complex topography of the terrain; the approximate difference method not only reduces the accuracy of the boundary condition and the implementation will allow it to produce a false scattering wave in the calculation, thus could not correctly calculate topographic effects on seismic wave propagation; although the method is put forward to improve it constantly, but the applicability of these methods still need further study in this research. The curvilinear grid finite difference method is used to describe the complex terrain grid; this method of discrete terrain can not only make the finite difference method to approximate the steps to avoid because of the false scattering wave, but also it can easily and accurately implement the free surface boundary conditions, so as to correctly calculate complex topographic effects on seismic wave propagation the process of seismic wave propagation. Numerical simulation is an important method for studying seismic wave propagation in anisotropic media, it is very important to study the seismic exploration and natural earthquake. The traditional finite difference method on a staggered grid center difference scheme; although this method can efficiently calculate the propagation of seismic waves in in isotropic media, but when arbitrary anisotropic media appear in the computational domain, it takes on seismic wave field in anisotropic medium. The amount of interpolation leads to large calculation error. This paper studies the curvilinear grid finite difference method with collocated grid DRP/opt MacCormack difference scheme; this method can not only avoid the traditional finite difference method to approximate the result of anisotropic seismic wave field in component interpolation, but also can be combined with grid technology and traction mirror technology to deal with complex terrain, so as to correctly calculate the process of seismic wave propagation in arbitrary anisotropic medium and the topographic effects on seismic wave propagation in anisotropic medium. Accurately simulating seismic wave through the solid - liquid interface for the study of physical processes in the ocean T phase and Stone wave, has an important the significance of seabed terrain scattering effect on seismic wave and ocean bottom seismometers recorded vibration data. The traditional finite difference method for complex The geometry of the solid - liquid interface and discrete solid liquid boundary conditions are done with a considerable degree of approximation, cannot accurately simulate the effect of solid - liquid interface on the ups and downs of the seismic wave propagation. The curvilinear grid finite difference method can accurately describe the complex geometry of the solid - liquid interface. A more direct implementation of boundary conditions, and provides a more reliable computational tool for studying the complex geometry of the solid - liquid interface effect on seismic wave propagation. This paper will integrate body fitted grid parity grid DRP/opt MacCormack difference scheme and traction image technology of curvilinear grid finite difference method are extended to in the case of two-dimensional and three-dimensional complex terrain and complex anisotropic medium containing solid - liquid interface shape of solid liquid coupling medium. The calculation of the seismic wave propagation, and the curve of grid finite difference calculation results respectively with generalized reflection / transmission coefficient method and spectral element method calculation results have been verified. The verification results show that the curve grid finite difference method can correctly and effectively applied to the problem of seismic wave propagation.

【學(xué)位授予單位】：中國(guó)科學(xué)技術(shù)大學(xué)
【學(xué)位級(jí)別】：博士
【學(xué)位授予年份】：2017
【分類(lèi)號(hào)】：P631.4

【參考文獻(xiàn)】

相關(guān)期刊論文 前2條

1 ;Seismogram Synthesis in Multi-layered Half-space Part Ⅰ. Theoretical Formulations[J];Earthquake Research in China;1999年02期

2 謝小碧,姚振興;二維不均勻介質(zhì)中點(diǎn)源P-SV波響應(yīng)的有限差分近似算法[J];地球物理學(xué)報(bào);1988年05期

，

本文編號(hào)：1472043