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