M-PML技術是由Meza-Fajardo和Papageorgiou(2008)提出,用來解決傳統(tǒng)PML的不穩(wěn)定問題。M-PML的基本理念是,波在正交方向上根據(jù)多個阻尼系數(shù)同時衰減,每個阻尼系數(shù)之間是成比例的。例如圖3-6中的二維PML模型所示,x方向上的阻尼可以被定義為:圖...[繼續(xù)閱讀]

    傳統(tǒng)單軸PML技術對于高泊松比(>0.38)淺地表地質(zhì)模型,是不穩(wěn)定的。泊松比越高,傳統(tǒng)PML算法發(fā)散越快。自由表面的存在造成了數(shù)值模擬時的不穩(wěn)定。自由表面復雜的波現(xiàn)象,是影響PML內(nèi)數(shù)值誤差快速累積的重要因素。對泊松比在0.1...[繼續(xù)閱讀]

    在一個二維各向同性介質(zhì)中,結(jié)合胡克定律,SH波的彈性動力學方程(Virieux,1984)為:(3-23)式中:vy(x,z,t)為時間域的速度場;ρ(x,z)為密度;μ(x,z)為剪切模量;σxy(x,z,t)和σzy(x,z,t)為剪應力;S(x,z,t)為震源。利用交錯網(wǎng)格高階有限差分對上述方程組進...[繼續(xù)閱讀]

    可通過對比解析解和數(shù)值解來驗證交錯網(wǎng)格有限差分算法模擬勒夫波的正確性。四分之一平面問題是兩個邊緣角度為90°的無限楔子的一個實際案例。Wait(1959)指出可以使用圖像理論找到這個問題的解。圖3-14　自由界面和剛性界面情...[繼續(xù)閱讀]

    通過兩層模型(參數(shù)見表3-4),對比其所得頻散曲線和理論相速度之間的關系,來驗證交錯網(wǎng)格FD代碼的正確性。表3-4　兩層模型參數(shù)空間網(wǎng)格由600×600節(jié)點,以及沿著模型左邊、右邊和底部過渡區(qū)的50個點組成,正如有限差分(每個波長至少...[繼續(xù)閱讀]

    f-k變換實質(zhì)上是一種二維傅里葉變換(Yilmaz,1987)。在時間域上,對一道地震信號做傅里葉變換,可以得到在時間上不同頻度(稱為頻率)的波動組分的振幅和相位信息。同樣地,在空間上,我們也可以對多道地震信號做類似傅里葉變換的數(shù)值...[繼續(xù)閱讀]

    τ-p變換法是由McMechan和Yedlin于1981年提出的。τ-p變換目前已廣泛地應用到地震資料處理的各個領域。它來源于圖像攝影理論中的拉東變換(RadonTransform),是離散化的拉東變換。τ-p變換是一種線性變換,它將時間-空間域的地震數(shù)據(jù)按不同...[繼續(xù)閱讀]

    相移法是由Park等在1998年提出的一種提取頻散曲線的方法。設時間-空間域的一個炮集記錄為U(x,t),對各道信號沿時間方向作一維傅里葉變換,其頻譜U(x,ω)可寫成如下兩項乘積的形式,即:U(x,ω)=p(x,ω)A(x,ω)(4-6)式中:p(x,ω)為相位譜;A(x,ω)為...[繼續(xù)閱讀]

    傾斜疊加法是由夏江海等在2007年提出的。該方法首先利用一個頻率掃描函數(shù)與炮集記錄(t-x域)進行卷積運算將頻率分解,把時間變換為頻率;接著利用類似于τ-p變換中傾斜疊加的方式得到每個速度的疊加能量值。由這兩步即得到了f...[繼續(xù)閱讀]

    高分辨率線性拉東變換法提取瑞雷波頻散曲線是由Luo等在2008年提出的。對二維波場d(x,t)的每一道做時域傅里葉變換得到d(x,f)后(圖4-18),每一個頻率的線性拉東變換可以由以下兩式計算:d(x,f)=m(p,f)ei2πfpx(4-15)和m(p,f)=d(x,f)e-i2πfpx(4-16)圖...[繼續(xù)閱讀]