1.示例模型示例模型的合成炮集記錄見圖4-20。模型參數(shù)見表4-1。圖4-20 示例模型的合成炮集記錄:60道,最小偏移距1m,道間距1m;采樣率1000Hz,記錄時間1.024s表4-1 二層理論模型參數(shù)2.程序部分代碼%%%%%%相關(guān)參數(shù)%%%%%%%%%%seis %地震數(shù)據(jù)...[繼續(xù)閱讀]
海量資源,盡在掌握
1.示例模型示例模型的合成炮集記錄見圖4-20。模型參數(shù)見表4-1。圖4-20 示例模型的合成炮集記錄:60道,最小偏移距1m,道間距1m;采樣率1000Hz,記錄時間1.024s表4-1 二層理論模型參數(shù)2.程序部分代碼%%%%%%相關(guān)參數(shù)%%%%%%%%%%seis %地震數(shù)據(jù)...[繼續(xù)閱讀]
將式(5-3)帶入柱坐標系統(tǒng),對于頻率為ω,相速度為Vr的軸對稱柱面瑞雷波,考慮其軸對稱性,由各層介質(zhì)的Navier方程引入Helmholtz勢分解得出以下的位移與應力關(guān)系式:式(5-4)中,ur,uθ,uz,τrz,τθz,τzz分別為位移應力矢量在柱坐標系下的各個分...[繼續(xù)閱讀]
考慮頻率為ω、相速度為Vr的地震波在圖5-1所示的n層水平層狀介質(zhì)中傳播,在除自由界面外的每個界面的位移分量和應力分量連續(xù),而自由界面處的兩個應力分量為零。這樣共有4n+2個邊界條件組成一個方程組,當且僅當方程組的系數(shù)行...[繼續(xù)閱讀]
圖5-2 二維水平層狀介質(zhì)模型考慮如圖5-2所示的層狀模型,其位移和應力[Pj(z) Sj(z) Tj(z)]對于瑞雷波(P-SV波)滿足下列關(guān)系:=(5-18)式中:k為波數(shù);ω為頻率;λ、μ為拉梅常數(shù);ρ為密度。其中:ζj=4μj(λj+μj)/(λj+2μj)ξj=1/(λj+2μj)令fi=T,則...[繼續(xù)閱讀]
與瑞雷波頻散方程推導方式類似,我們可以通過連續(xù)兩層上界面之間的關(guān)系,應力位移連續(xù)邊界條件、自由邊界條件及無窮遠輻射條件得到勒夫波的頻散函數(shù):其中,當1≤m≤n-2時,式中:Qm=krβmhm,rβm=-i。當m=n-1時且半空間的上界面為...[繼續(xù)閱讀]
將面波的頻散曲線方程看作為一個函數(shù):F(fj,VRj,VS,VP,d,h)=0或F(fj,VLj,VS,d,h)=0, j=1,2,3,…,m式中:fj為頻率點;VRj和VLj為該頻點的瑞雷波相速度;VS,VP,d,h分別為S波速度、P波速度、密度和層厚;m為總頻率點個數(shù)。把整個頻散曲線的求取過程看作...[繼續(xù)閱讀]
用一個簡單的兩層模型(表5-2)計算瑞雷波和勒夫波的頻散曲線,分別采用快速矢量傳遞算法和Knopoff方法,計算該二層模型的瑞雷波和勒夫波頻散曲線,并對其頻散曲線特征進行分析。其計算結(jié)果如圖5-4。在上述代碼中,取vp=[400,800],vs=[2...[繼續(xù)閱讀]
用一個兩層模型來驗證將復數(shù)根的實部當作瑞雷波相速度的正確性。該模型的表層厚度為5m,VS=420m/s,VP=900m/s,ρ=2000kg/m3,半空間VP=270m/s,VS=700m/s,ρ=2000kg/m3。圖5-7(a)為通過數(shù)值模擬得到的48道瑞雷波的炮集記錄。該炮集的最小偏移距...[繼續(xù)閱讀]
因為層數(shù)的多變性,多層模型情況有些復雜。采用較少層數(shù)的模型來代替多層模型來簡化問題,表5-3就用來驗證該想法的正確性。用三層模型來代替多層模型,將多層模型的頻散曲線作為實測數(shù)據(jù),并用Xia等(1999)討論的反演方法將其反演...[繼續(xù)閱讀]
當表層的橫波速度高于某些下伏地層時,由目前算法得到的瑞雷波相速度在其高頻部分會趨近除半空間以外的最低橫波層速度,而不是與表層橫波速度相關(guān)的某個值?;诓▌臃匠痰臄?shù)值模擬結(jié)果表明,瑞雷波的相速度在其波長遠小于...[繼續(xù)閱讀]