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

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

    在一個二維各向同性介質中,結合胡克定律,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)為震源。利用交錯網格高階有限差分對上述方程組進...[繼續閱讀]

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

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

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

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

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

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

    高分辨率線性拉東變換法提取瑞雷波頻散曲線是由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)圖...[繼續閱讀]