面波多模式頻散曲線正演PPT培訓(xùn)課件

面波多模式頻散曲線正演一、問(wèn)題的提出二、瑞雷波的基本原理和頻散方程三、固體層狀介質(zhì)中瑞雷波頻散曲線 5/20/20235/20/20235/20/20235/20/2023WaveTypePercentoftotalenergyRayleighShearCompression

672675/20/202360年代初，Hoykallen－半波長(zhǎng)解釋法，將穩(wěn)態(tài)瑞雷波法首先應(yīng)用于地基勘察80年代初，日本VIC株式會(huì)社－GR-810佐藤式全自動(dòng)地下勘探系統(tǒng)1989年，楊成林自行研制了一套穩(wěn)態(tài)瑞雷波勘探系統(tǒng)，用于第四系地層分層和地基處理效果評(píng)價(jià)瑞雷波勘探的歷史回顧

20世紀(jì)50年代初，開始利用人工地震波中的瑞雷波進(jìn)行工程勘察穩(wěn)態(tài)瑞雷波法5/20/20231986年,Nazarian利用表面波譜分析方法（SASW）對(duì)高速公路路面及路基進(jìn)行了探測(cè),為瞬態(tài)瑞雷波法在工程中的廣泛應(yīng)用奠定了基礎(chǔ)1996年，劉云楨等自行研制了SWS瞬態(tài)面波多道數(shù)據(jù)采集處理系統(tǒng),將其應(yīng)用于機(jī)場(chǎng)工程勘察、淺層煤田勘探、地下煤巷探測(cè)等方面的工作中1996年,李哲生應(yīng)用瞬態(tài)多道瑞雷波法對(duì)某建筑物場(chǎng)地地層分層進(jìn)行勘察1998年,肖柏勛等研制的LXII巖土工程質(zhì)量檢測(cè)分析儀中包含的瞬態(tài)多道瑞雷波勘探的子系統(tǒng)瞬態(tài)瑞雷波法5/20/2023數(shù)據(jù)的采集和處理正演問(wèn)題(快速、穩(wěn)定的算法)反演問(wèn)題(線性及非線性算法)瑞雷面波干擾波(在反射地震中)勘探波(頻散特性)5/20/2023多模問(wèn)題數(shù)值計(jì)算存在多種模式什么特征怎樣解釋怎樣與實(shí)測(cè)中的一條頻散曲線相對(duì)應(yīng)瑞雷波反演局部線性化方法：阻尼最小二乘方法。對(duì)初始模型要求較高非線性反演方法：遺傳算法、神經(jīng)網(wǎng)絡(luò)算法、蒙特卡羅法等。計(jì)算量大反演速度慢。5/20/2023一、問(wèn)題的提出二、瑞雷波的基本原理和頻散方程三、固體層狀介質(zhì)中瑞雷波頻散曲線 5/20/20232.1面波的形成面波的形成可以用一個(gè)簡(jiǎn)單的例子來(lái)說(shuō)明。下圖是地表下某一厚度為H的覆蓋層的模型，如果一體波的簡(jiǎn)諧波在層內(nèi)的傳播滿足全反射的條件，射線路徑為ABCDEF，則根據(jù)波前面與射線垂直的性質(zhì)，虛線CF即可代表波由C、D、E、F經(jīng)兩次反射后到達(dá)F的地震波波前，也可代表由A傳播到C的后續(xù)振動(dòng)的波陣面。如兩者波程差CDEF正好等于波長(zhǎng)的整數(shù)倍，則它們們完全同相，其合成或疊加屬于相長(zhǎng)干涉，并形成一種沿界面行進(jìn)的次生波，這種波的能量更強(qiáng)且主要能量集中在地表附近。

DHAEBCF5/20/2023彈性波動(dòng)方程：縱橫波波動(dòng)方程：瑞雷波存在條件：

α為縱波速度β為橫波速度,c為瑞雷波速度5/20/20232.2瑞雷波的傳播及其特點(diǎn)（1）沿地表水平方向傳播，振幅沿垂向方向按指數(shù)規(guī)律衰減的一種非均勻平面波。（2）瑞雷波沿水平方向傳播的速度小于橫波速度,（3）其能量分布僅限于自由表面垂向方向?yàn)閮杀度鹄撞úㄩL(zhǎng)范圍的薄層介質(zhì)內(nèi)。（4）其垂向位移分量較之水平位移分量超前相位（5）瑞雷波質(zhì)點(diǎn)位移的軌跡是一個(gè)逆進(jìn)的橢圓，如下面兩圖所示：5/20/2023（6）若半無(wú)限均勻彈性介質(zhì)含有一個(gè)疏松薄覆蓋層，則在自由表面存在具有頻散的瑞雷波，傳播速度是頻率的函數(shù)（7）在半無(wú)限均勻彈性介質(zhì)中的瑞雷波速度與頻率無(wú)關(guān)，即無(wú)頻散。（8）由實(shí)際觀測(cè)和數(shù)學(xué)模型研究知道，瑞雷波能量強(qiáng)，頻率低。（9）在兩層或者多層分層介質(zhì)中，瑞雷波的相速度都隨頻率變化，這種情況稱為面波的頻散。5/20/2023頻散：波初始擾動(dòng)的每一個(gè)簡(jiǎn)諧成分都以不同的速度前進(jìn)，從而初始擾動(dòng)的形狀將在前進(jìn)中發(fā)生變化，擾動(dòng)將分散成一系列波，這即為所謂的頻散現(xiàn)象。

正頻散：逆頻散：群速度：相速度：2.3瑞雷波的頻散和相速度a圖：粘彈性固體中頻散導(dǎo)致的波形變化b圖：能量轉(zhuǎn)移到后面的周期導(dǎo)致波形變化5/20/20232.3瑞雷波頻散曲線正演方法1、Thomson－Haskell法1953年，Haskell在Thomson的基礎(chǔ)上，通過(guò)相鄰兩界面的傳遞矩陣公式以及自由表面邊界條件和無(wú)窮遠(yuǎn)處的輻射條件，導(dǎo)出了層狀介質(zhì)中平面瑞雷波的頻散方程，人們把這種方法稱為Thomson－Haskell法。Knopoff在Haskell的矩陣方法基礎(chǔ)上，提出了一種新的求取頻散曲線的方法。與Haskell不同的是Knopoff得出了4n＋2階的頻散方程行列式，此行列式的求解使用了Knopoff分解法（一種循環(huán)的Laplace行分解法）。1967年，Randall對(duì)Knopoff分解進(jìn)行了計(jì)算，成功地避免了Thomson－Haskell方法中存在的數(shù)值不穩(wěn)定性問(wèn)題，但是此方法過(guò)于復(fù)雜，不利于廣泛應(yīng)用。

對(duì)于層狀介質(zhì)中面波頻散曲線的正演問(wèn)題，多年來(lái)大量學(xué)者進(jìn)行了研究，并提出了各種方法，取得了好的效果，也為其它復(fù)雜介質(zhì)中面波的正演問(wèn)題奠定了基礎(chǔ)。其主要的方法:5/20/20232、Schwab－Knopoff法Schwab和Knopoff對(duì)Knopoff分解在數(shù)值上進(jìn)行了進(jìn)一步的研究和改進(jìn)，形成了完整的Schwab－Knopoff方法1970年，Schwab通過(guò)變換使得頻散函數(shù)以特別簡(jiǎn)單的代數(shù)形式表達(dá)出來(lái)，并編制了計(jì)算量小的程序，這種方法稱為快速Schwab－Knopoff方法3、δ矩陣法1963年，Pestel和Leckie介紹了δ矩陣法1965年，Thrower和Dunkin首次將δ矩陣方法應(yīng)用于瑞雷波的頻散方程推導(dǎo)，避免了數(shù)值精度丟失的問(wèn)題。1970年，Watson應(yīng)用簡(jiǎn)化的δ矩陣求解瑞雷波頻散問(wèn)題，減小了計(jì)算量。1996年，Buchen和Ben－Hadon提出快速δ矩陣方法1998年，Ivasson用更簡(jiǎn)便的方法同樣得到了快速δ矩陣算法。5/20/20234、Abo－Zena法1979年，Abo-Zena為解決頻散方程的數(shù)值不穩(wěn)定性問(wèn)題，通過(guò)一系列的4×4階反對(duì)稱矩陣的循環(huán)計(jì)算得到了瑞雷波的頻散方程1979年，Menke解釋了Abo－Zena提出的反對(duì)稱矩陣的物理含義1968年，Ben－Menahem等提出B、P、C坐標(biāo)系，大大簡(jiǎn)化了層狀介質(zhì)中三維聲波場(chǎng)問(wèn)題的求解過(guò)程1982年，李幼銘等在B、P、C坐標(biāo)下采用的δ矩陣運(yùn)算用五個(gè)形式簡(jiǎn)單的矩陣連乘形式重新組織了Abo－Zana算法，使運(yùn)算步驟大量簡(jiǎn)化1996年，張碧星簡(jiǎn)化了傳遞矩陣的乘積運(yùn)算2001年，凡友華等從柱面波出發(fā)，討論了頻散函數(shù)的求取問(wèn)題。2002年，ShuangXZhang對(duì)各向異性介質(zhì)中瑞雷波的頻散特征進(jìn)行了系統(tǒng)研究。5、RT矩陣法1974年及1979年，Kennett等人提出并發(fā)展了RT矩陣祛，此方法使用了反射和透射子矩陣技術(shù)，巧妙地避開了數(shù)值不穩(wěn)定性問(wèn)題。但這種方法要求所有的計(jì)算均使用復(fù)數(shù)運(yùn)算形式，故計(jì)算量較大1993年，陳曉非在Luco和Apsel的廣義反射－透射系數(shù)方法基礎(chǔ)上提出了一種系統(tǒng)、有效的算法來(lái)求解固體層狀介質(zhì)中的面波振型5/20/2023…….xz自由表面空氣或者液體無(wú)限厚底層

u，w沿x,z方向的位移，波數(shù)切向應(yīng)力法向應(yīng)力

密度；

厚度；縱波波速橫波波速

Lame常數(shù)；對(duì)于第m層介質(zhì)：Knopoff快速計(jì)算法5/20/2023當(dāng)m+1為偶數(shù)時(shí)有：

其中:當(dāng)m+1為奇數(shù)時(shí)有:其中:

（a）（b）5/20/2023

無(wú)液體表層存在液體表層頻散函數(shù)是由以下實(shí)數(shù)初值遞推的：

將（c）代入（a）和（b）中反復(fù)遞推，直到頻散函數(shù)計(jì)算到（n-1）界面為止。（c）5/20/2023最后的頻散函數(shù)為5/20/2023頻散曲線算法流程框圖輸出頻散曲線給對(duì)應(yīng)模式下賦初值和迭代的精度E給定模型的參數(shù)以及需要計(jì)算的頻率輸出相應(yīng)頻率下的瑞雷波速度計(jì)算頻散函數(shù)Fi5/20/20232.5多模式面波頻散曲線的計(jì)算1、選取一定的相速度初值，給定一段相應(yīng)的頻率步長(zhǎng)和范圍，運(yùn)用線性插值法來(lái)計(jì)算得到與頻率相對(duì)應(yīng)的瑞雷波相速度，并且通過(guò)在計(jì)算的過(guò)程中不斷的修正得到的瑞雷波速度來(lái)提高計(jì)算瑞雷波頻散函數(shù)的精度和計(jì)算頻率。2、先利用上述方法計(jì)算基階模式的瑞雷波速度，得到基階模式的頻散曲線3、將基階模式的瑞雷波速度加上相應(yīng)的步長(zhǎng)，重復(fù)上面的步驟，得到第二階模式的瑞雷波速度4、以及類推，可以得到第3階、第4階、······第n階模式的頻散曲線，直到第n階模式的頻散曲線的截止頻率大于相應(yīng)的終止頻率為止5、一般情況下，階數(shù)越高，得到面波在高頻頻率段下的位移越小，相對(duì)應(yīng)的能量也就越弱。所以根據(jù)一般的實(shí)際需要和意義，我們一般計(jì)算到3－5階瑞雷波頻散曲線模式來(lái)分析其相應(yīng)的瑞雷波頻散曲線的特征5/20/2023編程計(jì)算需要注意的關(guān)鍵問(wèn)題：1、如何選擇初始速度。在基階模式頻散曲線的求取中，我們一般使用的初始速度為最小橫波速度的某個(gè)百分比（比如80％）；而在高階模式頻散曲線的求取中，我們使用前一階模式的瑞雷波速度加上相應(yīng)的步長(zhǎng)。2、頻散方程計(jì)算過(guò)程中的復(fù)數(shù)問(wèn)題。上述頻散方程在求取過(guò)程中全部使用了實(shí)數(shù)運(yùn)算，這樣大大減少了運(yùn)算量，也便于在計(jì)算機(jī)上實(shí)現(xiàn)編程，但是事實(shí)證明：在使用頻散方程求取瑞雷波速度的過(guò)程中，由于迭代速度過(guò)大（大于某一層介質(zhì)的橫波速度）導(dǎo)致運(yùn)算過(guò)程中仍然會(huì)出現(xiàn)虛數(shù)運(yùn)算的情況，這時(shí)我們對(duì)相應(yīng)的頻散方程中的各個(gè)元素，采用分別討論編程計(jì)算的方法避免虛數(shù)的出現(xiàn)。3、對(duì)高階模計(jì)算過(guò)程中速度步長(zhǎng)的問(wèn)題。速度步長(zhǎng)過(guò)大，可能導(dǎo)致漏根以至于漏掉了某一高階模式；而速度步長(zhǎng)過(guò)小，迭代次數(shù)過(guò)多又會(huì)導(dǎo)致計(jì)算量加大。5/20/2023一、問(wèn)題的提出二、瑞雷波的基本原理和頻散方程三、固體層狀介質(zhì)中瑞雷波頻散曲線5/20/2023橫波速度遞增的固體介質(zhì)模型，它模擬了工程勘探中最為常見(jiàn)的速度遞增型的地質(zhì)情況，如存在下伏高速基巖的黃土地區(qū)以及表面為土壤或者砂，下面為第四系砂、泥沉積，再往下為壓實(shí)底層的灘涂地區(qū)地質(zhì)結(jié)構(gòu)等。3.1.1橫波速度隨深度遞增的兩層模型

先討論上層介質(zhì)厚度不同時(shí)兩層模型的瑞雷波頻散曲線模型1模型號(hào)厚度h(m)1400025002500h122695000300028002743縱波速度橫波速度瑞雷波速度密度β(m/s)h250030002000h15/20/2023(a)(b)(c)(d)基階模式的零頻極限為高頻極限為

高階模式在截止頻率處的相速度為高階模式的高頻極限為

隨著瑞雷波模式的增加，各高階模式的截止頻率也隨之變大隨著上層介質(zhì)厚度的增加，較低頻率范圍內(nèi)的瑞雷波頻散曲線的模式階數(shù)越來(lái)越多5/20/2023橫波速度遞減的固體介質(zhì)模型，它在實(shí)際生產(chǎn)中大多都為道路型結(jié)構(gòu)的地球物理模型（表層為水泥混凝土或者瀝青之類的高速介質(zhì)）3.1.2橫波速度隨深度增加而遞減的兩層模型模型5模型組號(hào)厚度h(m)5500030002800h127434000250025002269縱波速度橫波速度瑞雷波速度密度其頻散曲線只有一個(gè)模式，且頻散曲線均表現(xiàn)為逆頻散，并且都存在截止頻率零頻極限為

截止頻率處的瑞雷波相速度都為當(dāng)上層介質(zhì)的厚度變小時(shí)，瑞雷波頻散曲線模式的截?cái)囝l率隨著變大，但是在相同頻率下的瑞雷波相速度變小C<β2的解釋：頻散方程的實(shí)數(shù)解對(duì)應(yīng)著瑞雷波，復(fù)數(shù)解對(duì)應(yīng)著泄漏模式波，對(duì)于長(zhǎng)距離的傳播，泄漏模式波的影響很小，我們?cè)诖藢⑵浜雎驭?m/s)h250030002000h15/20/2023模型6模型組號(hào)厚度h(m)6600031002800h128845000300025002743縱波速度橫波速度瑞雷波速度密度該模型的頻散曲線同樣表現(xiàn)為逆頻散，并且只有單一模式，但是其頻散曲線不存在截止頻率，頻散曲線的相速度的高頻極限為在較大的頻率范圍內(nèi)（小于1400Hz）當(dāng)上層介質(zhì)的厚度變小時(shí)，相同頻率下的瑞雷波相速度也變小。β(m/s)h300031002900h15/20/2023模型7模型組號(hào)厚度h(m)7500030002800h127435990296045002763縱波速度橫波速度瑞雷波速度密度曲線首先隨頻率的增大而遞增，然后在約200Hz附近后又隨頻率的增大而遞減；即在約200Hz以前為逆頻散，而200Hz以后則為正頻散，且頻散曲線在200Hz附近出現(xiàn)一個(gè)極大值點(diǎn)零頻極限為高頻極限為當(dāng)上層介質(zhì)的厚度變小時(shí)，整個(gè)頻散曲線的形態(tài)向高頻方向移動(dòng)β(m/s)h296030002900h15/20/2023模型組號(hào)厚度h(m)7’500030002800h127435990296025002763縱波速度橫波速度瑞雷波速度密度模型7密度變化此時(shí)的瑞雷波頻散曲線呈現(xiàn)出來(lái)的頻散曲線形態(tài)相反，原來(lái)表現(xiàn)為逆頻散的，現(xiàn)在表現(xiàn)為正頻散，而原來(lái)表現(xiàn)為逆頻散的，現(xiàn)在則表現(xiàn)為正頻散現(xiàn)象，并且此時(shí)瑞雷波的頻散曲線呈現(xiàn)一個(gè)極小值。當(dāng)表層介質(zhì)的厚度變小時(shí)，頻散曲線的極大值明顯向高頻方向移動(dòng)。頻散曲線在零頻極限與高頻極限的相速度趨勢(shì)情況與原來(lái)是相同的。5/20/2023模型組號(hào)厚度h(m)7－A20001400ρ1

h1274324001350ρ22763縱波速度橫波速度瑞雷波速度密度(a)(b)(c)(d)1、說(shuō)明密度在一定頻率和速度的情況下，對(duì)頻散曲線的形態(tài)影響也是比較大的。2、在一般的情況下，大多數(shù)學(xué)者都以(a)、(b)、(d)三個(gè)密度參數(shù)下的頻散曲線來(lái)研究其形態(tài)的變化3、從一定程度上說(shuō)明了瑞雷波頻散曲線的復(fù)雜以及利用瑞雷波進(jìn)行勘察研究的困難程度5/20/2023模型8：中間層厚度不同時(shí)三層模型3.2.1橫波速度隨深度遞增的三層模型上節(jié)我們已經(jīng)討論過(guò)兩層介質(zhì)情況下表層厚度不同時(shí)的頻散曲線形態(tài)的變化，因此，在三層厚度不同時(shí)橫波速度遞增的水平層狀介質(zhì)中，我們主要研究中間層厚度變化時(shí)所引起的頻散曲線形態(tài)特征的變化。模型組號(hào)厚度h(m)840002500250062269500030002800h227436000350030003213縱波速度橫波速度瑞雷波速度密度β(m/s)h2500300020003500h1h25/20/2023高頻極限的速度則隨著階數(shù)的變大而變大，并且互不相同。可以發(fā)現(xiàn)，隨著中間層介質(zhì)厚度的遞增，頻散曲線具有水平段的趨勢(shì)越來(lái)越明顯。相鄰高階模式在某一頻段處的間隔越來(lái)越小，即相同頻率段范圍內(nèi)，瑞雷波模式階數(shù)越來(lái)越多,通過(guò)我們加密相應(yīng)的頻率段,發(fā)現(xiàn)它們并不相交(a)(b)(c)(d)基階模式的零頻極限為

高頻極限為各高階模式均存在截止頻率，高階模式在截止頻率處的速度為最底層介質(zhì)的橫波速度

（d）圖第4、5階模式在某段頻率處的頻散曲線5/20/20233.2.2橫波速度隨深度而遞減的三層模型模型9模擬黃土地區(qū)混凝土地面下的等效模型模型組號(hào)厚度h(m)9600031003000h13213500030002800h227434000250025002269縱波速度橫波速度瑞雷波速度密度頻散曲線均只有一個(gè)瑞雷波模式，瑞雷波模式存在截?cái)囝l率，它是因?yàn)槿鹄撞ㄋ俣冗_(dá)到最底層介質(zhì)橫波速度的緣故，當(dāng)中間層和上層介質(zhì)變薄時(shí)，相同頻率下的相速度都變小零頻極限為β(m/s)h2500300020003100h1h25/20/2023模型10模型組號(hào)厚度h(m)10600031002200h12884500030002800h227435990296045402763縱波速度橫波速度瑞雷波速度密度頻散曲線的瑞雷波模式也只有一個(gè)，此時(shí)的瑞雷波頻散曲線不存在截止頻率，并且存在正頻散和逆頻散兩種現(xiàn)象。在較低頻率段（約小于10Hz）和較高的頻率段（約大于300Hz）時(shí)，瑞雷波頻散曲線表現(xiàn)為正頻散現(xiàn)象，而在大約10Hz到300Hz的頻率范圍內(nèi)，瑞雷波的頻散模式為逆頻散現(xiàn)象。而且可以看出，在300Hz左右存在一個(gè)極大值，而在10Hz左右存在一個(gè)極小值高頻極限為當(dāng)中間層變薄時(shí)，整個(gè)頻散曲線的形態(tài)變化不大，但是當(dāng)上層和中間層厚度都為3m時(shí)，瑞雷波頻散曲線的極大值變大，極小值變小，并且逆頻散的范圍變大零頻極限為h2960300029003100β(m/s)h1h25/20/20233.2.3含有低速夾層時(shí)的三層固體模型(a)存在低速夾層時(shí)的固體介質(zhì)模型，它模擬了地層中出現(xiàn)地裂縫、空洞、暗穴、塌陷以及破碎帶，還有諸如油氣勘探中的儲(chǔ)層、高強(qiáng)度材料中的裂隙、氣泡等情況。模型組號(hào)厚度h(m)11400025002500h12269235011201060h210485000300028002743縱波速度橫波速度瑞雷波速度密度存在低速夾層的瑞雷波頻散曲線的復(fù)雜性更為明顯，基階模式表現(xiàn)為正頻散和逆頻散現(xiàn)象，分別在大約50Hz和300Hz表現(xiàn)為局部極小值和局部極大值?；A模式相速度的零頻極限為

高階模式都存在截止頻率，且在截止頻率處的瑞雷波相速度

h2500300010001120β(m/s)h1h2高頻極限趨近中間層介質(zhì)的瑞雷波速度5/20/2023(b)(c)在約1075Hz附近，模型11中間層厚度為2m的瑞雷波第4階和第5階模式的頻散曲線并不相交上層厚度和中間層介質(zhì)的厚度變化時(shí)其頻散曲線的形態(tài)總體變化不大，但在(b)圖中的高階模式的頻散曲線，在很大頻段出現(xiàn)瑞雷波相速度相等的水平段區(qū)域，并在水平段的速度為第一層介質(zhì)的瑞雷波速度5/20/2023模型8：速度遞增型模型模型11：存在低速夾層模型在橫波速度遞增的模型中，基階模式的瑞雷波能量占主要地位，左圖中高階模式水平段的出現(xiàn)并不影響實(shí)際勘探中得到的頻散曲線，但在模型11的頻散曲線中，張碧星等）計(jì)算出水平段在地表的能量在其頻率范圍內(nèi)是最強(qiáng)的，這也說(shuō)明當(dāng)模型存在低速夾層時(shí)，頻散曲線在水平段瑞雷波的能量就主要集中在第一層介質(zhì)內(nèi)，符合高頻瑞雷波探測(cè)深度受限的物理機(jī)制。同時(shí)，配合基階零頻與高頻極限，再借助“水平段”特征，容易提取底層瑞雷波相速度和橫波速度及其關(guān)系，進(jìn)一步求出中間夾層的橫波速度，有助于對(duì)目的夾層的其它彈性參數(shù)進(jìn)行計(jì)算。

5/20/2023

代表了壓實(shí)路面（地面）下伏軟弱層（或空腔）的工程情況模型12模型組號(hào)厚度h(m)1260003100220062884235011201060h210485000300028002743縱波速度橫波速度瑞雷波速度密度當(dāng)中間層的厚度變小時(shí)的頻散曲線則更為復(fù)雜，如右圖所示，除了在相同頻段內(nèi)瑞雷波頻散曲線的導(dǎo)波模式減少以外，瑞雷波的各階模式都在水平方向上被“拉長(zhǎng)”，然后不同模式之間瑞雷波頻散曲線的相交趨勢(shì)更加猛烈，第1階與第2階，第2階與第3階，第3階與第4階等之間都有近似相交的趨勢(shì)β(m/s)3000310010001120hh1h25/20/2023模型12中第1階模式和第2階模式處的頻散曲線，可以發(fā)現(xiàn)它們并不相交

模型組號(hào)厚度h(m)135000300028006274323501120106031048400025002500226914500030002800h12743235011201060310485000300028002743縱波速度橫波速度瑞雷波速度密度h模型13β(m/s)300010001120h模型14h1h2β(m/s)2500300010001120h1h25/20/2023

模型13模型14與前面兩種存在低速夾層時(shí)的瑞雷波頻散曲線類似，說(shuō)明低速夾層對(duì)面波頻散的影響占主導(dǎo)地位。模型14實(shí)際上是構(gòu)筑物澆灌時(shí)內(nèi)部存在質(zhì)量問(wèn)題的情況，如橋面、橋墩、樁基等，其頻散曲線的理論研究可以為超聲波檢測(cè)提供技術(shù)依據(jù)。

5/20/2023

當(dāng)?shù)貙又谐霈F(xiàn)低速軟夾層的時(shí)候，頻散曲線不是從逐漸過(guò)渡到也就是說(shuō)，當(dāng)頻率較高時(shí)，實(shí)地勘探中接受道的肯定不止是只有基階模式，而是存在高階模式，也就是說(shuō)，接收到的信號(hào)肯定是多階瑞雷波模式的疊加。另外，此時(shí)模型相速度的頻散曲線都存在近似水平段，此時(shí)的瑞雷波相速度為

，而近似水平段又不是同一個(gè)模式，當(dāng)頻率增大時(shí)，接收到的導(dǎo)波必定以某一模式為主，在另一頻率段又以另一模式為主，這可能由于各導(dǎo)波間的能量存在跳躍，從而導(dǎo)致了“之”字形頻散現(xiàn)象。5/20/20233.2.4含有高速夾層時(shí)的三層固體模型

模型15存在高速夾層時(shí)的固體介質(zhì)模型，它模擬了疏松表層黃土區(qū)、泥灘區(qū)以及松散路面公路地區(qū)，探測(cè)下伏高速地質(zhì)不均勻體時(shí)的情況。

模型組號(hào)厚度h(m)15500030002500h12743600031002800h228844000250022002269縱波速度橫波速度瑞雷波速度密度零頻極限

截止頻率處的瑞雷波速度為最底層介質(zhì)的橫波速度

當(dāng)中間層厚度和上層介質(zhì)的厚度減小時(shí)，頻散曲線的截止頻率向高頻方向移動(dòng)。并且在相同頻率下的瑞雷波相速度變小。

β(m/s)3000310010002500hh1h25/20/2023模型16模型組號(hào)厚度h(m)16500030002500h12743600031002800h228845990296024502763縱波速度橫波速度瑞雷波速度密度瑞雷波的頻散曲線存在著正頻散和逆頻散現(xiàn)象，且都在100Hz左右出現(xiàn)局部極小值，并隨著頻率的增大出現(xiàn)一個(gè)極大值，零頻極限為高頻極限為當(dāng)中間層的厚度h2減小時(shí)，相速度的曲線幅值降低，曲線向高頻方向移動(dòng)，當(dāng)h1變小，相速度的曲線幅值增大，且其峰值位置向高頻方向移動(dòng)β(m/s)3000310010002960hh1h25/20/2023

模型17模型組號(hào)厚度h(m)17400025002200h12269600031002800h228845000300025002743縱波速度橫波速度瑞雷波速度密度這個(gè)模型是唯一存在多個(gè)瑞雷波模式的中間層為高速夾層的三層模型，各瑞雷波模式之間互不相交，并且都表現(xiàn)為正頻