平面S波在飽和多孔熱彈性介質(zhì)邊界的反射問(wèn)題研究

柳鴻博, 周鳳璽,2, 郝磊超

(1. 蘭州理工大學(xué) 土木工程學(xué)院, 甘肅 蘭州 730050; 2. 西部土木工程防災(zāi)減災(zāi)教育部工程研究中心, 甘肅 蘭州 730050)

0 引言

隨著石油工程、化工、路面工程和核廢料管理等領(lǐng)域的快速發(fā)展,飽和多孔熱彈性介質(zhì)中的波的傳播和反射問(wèn)題已經(jīng)引起了人們的廣泛關(guān)注。Biot[1-2]首先建立了在等溫和高低頻條件下流體飽和多孔彈性介質(zhì)中波傳播的理論框架。隨后國(guó)內(nèi)外眾多學(xué)者[3-9]借助經(jīng)典的Biot波動(dòng)理論從不同角度對(duì)不耦合熱的兩相飽和多孔彈性介質(zhì)中波的傳播特性進(jìn)行了大量的研究工作。與等溫條件下的多孔介質(zhì)相比較,飽和多孔熱彈性介質(zhì)由于孔隙中的流-固兩相的相互作用以及熱和力之間的相互耦合,故而其波動(dòng)響應(yīng)問(wèn)題變得更為復(fù)雜。Pecker等[10]研究了溫度對(duì)液飽和多孔彈性介質(zhì)中體波的傳播特性的影響,并給出了相應(yīng)的方程,但沒(méi)有考慮面波的傳播行為。Singh[11]基于Youssef提出的熱彈性理論,研究了飽和多孔熱彈性介質(zhì)中面波的傳播特性,但該模型比較復(fù)雜,其中一些物理參數(shù)的意義不夠明確。劉干斌等[12,13]基于Biot波動(dòng)理論及熱彈性動(dòng)力理論,建立了飽和多孔彈性介質(zhì)熱-流-固完全耦合的動(dòng)力響應(yīng)模型及控制方程,對(duì)飽和多孔彈性介質(zhì)的熱動(dòng)力響應(yīng)進(jìn)行了研究。

雖然目前國(guó)內(nèi)外學(xué)者[14-22]針對(duì)飽和多孔彈性介質(zhì)中波的傳播和反射問(wèn)題做了一些研究,但對(duì)彈性波在飽和多孔熱彈性介質(zhì)中的反射問(wèn)題涉及較少。因此,對(duì)于彈性波在飽和多孔熱彈性介質(zhì)表面的反射問(wèn)題的研究顯得十分重要。本文在熱-流-固耦合動(dòng)力響應(yīng)模型的基礎(chǔ)上,建立了飽和多孔熱彈性介質(zhì)的熱-流-固耦合波動(dòng)方程,通過(guò)引入勢(shì)函數(shù)并結(jié)合位移矢量和邊界條件,推導(dǎo)得到了向飽和多孔熱彈性介質(zhì)邊界入射平面S波后各反射波振幅反射率的理論表達(dá)式,最后通過(guò)參數(shù)分析討論了土體熱膨脹系數(shù)、入射頻率和入射角對(duì)各反射波的影響。

1 飽和多孔熱彈性介質(zhì)波動(dòng)方程

飽和多孔介質(zhì)中熱彈性本構(gòu)關(guān)系可表示為[12]:

σij=λeδij+2μεij-αpwδij-Kacθδij

(1)

式中:δij為Kronecker符號(hào);σij為總應(yīng)力;εij為應(yīng)變;e為體應(yīng)變;λ和μ為多孔彈性介質(zhì)的Lamb常數(shù);pw為孔隙水壓力;α=1-K/Ks(α≤ 1)為飽和多孔介質(zhì)的壓縮系數(shù),Ks為固體顆粒的體積模量(GPa);K=λ+2μ/3為土體在排水條件下的體積模量(GPa);ac為土體的熱膨脹系數(shù)(℃-1);θ=T-T0為溫度變化量,T和T0分別為絕對(duì)溫度和初始溫度(℃)。

忽略體積力后,飽和多孔介質(zhì)的運(yùn)動(dòng)方程可以表示為:

(2)

式中:ui和wi分別表示i方向的位移和流體的相對(duì)位移;ρ=(1-n)ρs+nρw為飽和多孔介質(zhì)的密度,ρs和ρw分別為固體顆粒和流體的密度(kg/m3),n為孔隙率。

應(yīng)變關(guān)系為:

(3)

考慮溫度的影響下,飽和多孔介質(zhì)的滲流方程為[12]:

(4)

式中:M為Biot模量,1/M=n/Kw+(α-n)/Ks,Kw為流體的體積模量(GPa);ξ為流體相對(duì)應(yīng)變;au=naw+(1-n)as-(1-α)ac,as和aw分別為固體顆粒和孔隙水的熱膨脹系數(shù)(℃-1)。

非等溫條件下飽和多孔介質(zhì)的液相連續(xù)方程為[12]:

(5)

式中:kl為飽和多孔介質(zhì)的滲透系數(shù)(m/s);DT為熱力滲透作用系數(shù)(m2/s·℃),反映飽和多孔介質(zhì)中溫度梯度對(duì)滲流的影響;g為重力加速度(m/s2)。

根據(jù)廣義的熱彈性理論,飽和多孔介質(zhì)中的熱傳導(dǎo)方程可以表示為[23]:

(6)

式中:k=(1-n)ks+nkw為飽和多孔熱彈性介質(zhì)的熱傳導(dǎo)系數(shù),ks和kw分別為固體顆粒和孔隙中流體的熱傳導(dǎo)系數(shù)(J/s·m·℃);C=(1-n)ρscs+nρwcw為土體比熱(J/m3·℃),cs和cw分別為固體顆粒和孔隙流體的比熱(J/kg·℃);t0為松弛時(shí)間(s)。式(6)為溫度控制方程,與運(yùn)動(dòng)控制方程式(2),(4)和式(5)相互耦合。由于經(jīng)典的溫度控制方程采用Fourier熱傳導(dǎo)定律,因此溫度控制方程為擴(kuò)散方程,表示熱以無(wú)限大速度傳播,與物理事實(shí)不相符,而式(6)考慮了熱波(T波)的有限傳播速度,描述了一個(gè)有限速度的彈性波與有限速度的熱波的耦合系統(tǒng)。

結(jié)合式(1)～(6),經(jīng)過(guò)推導(dǎo)可以得到飽和多孔熱彈性介質(zhì)的波動(dòng)方程為:

(7a)

(7b)

(7c)

式中:

由平面波條件e=divu和ξ=-divw可得到位移矢量形式的飽和多孔熱彈性介質(zhì)的波動(dòng)方程為:

μ▽2u+a11▽(▽·u)-a12▽(▽·w)+

(8a)

a21▽(▽·u)-a22▽(▽·w)+a23▽?duì)?

(8b)

為便于分析,引進(jìn)兩相介質(zhì)位移矢量的Helmholtz分解形式,即:

u=▽?duì)譻+▽×Hs

(9a)

w=▽?duì)譿+▽×Hw

(9b)

式中:ψs和ψw分別為固體骨架、孔隙水的標(biāo)量勢(shì)函數(shù);Hs和Hw分別為固體骨架、孔隙水矢量勢(shì)函數(shù)。其中:

ψs=Bsexp[i(kpnr-ωt)]

(10a)

ψw=Bwexp[i(kpnr-ωt)]

(10b)

Hs=Csexp[i(ksnr-ωt)]

(10c)

Hw=Cwexp[i(ksnr-ωt)]

(10d)

設(shè)T波的勢(shì)函數(shù)方程為:

θ=BTexp[i(kpnr-ωt)]

(11)

將式(9a)、(9b)和(11)代入式(8a)～(8c),并對(duì)方程兩端進(jìn)行散度和旋度運(yùn)算,將波動(dòng)方程解耦得到以下特征方程:

(12a)

(12b)

其中,式(12a)和(12b)分別表示飽和多孔熱彈性介質(zhì)中的壓縮波和剪切波的特征方程,兩個(gè)特征方程中的元素分別表示為:

式中:kp和ks分別表示飽和多孔熱彈性介質(zhì)中P波和S波的復(fù)波數(shù)。

對(duì)于式(12a)可以解出6個(gè)不同的復(fù)波數(shù)kp=Re(kp)+i Im(kp),式(12b)可以解出兩個(gè)不同的復(fù)波數(shù)ks=Re(ks)+iIm(ks),其中Re和Im分別為實(shí)部和虛部,Re反映常規(guī)波數(shù),Im反映波的衰減特性。由于振幅沿著波傳播的方向衰減,則Im(kp)>0,Im(ks)>0,故而kp只有3個(gè)有意義的復(fù)根,其中兩個(gè)為壓縮P1波和P2波的復(fù)波數(shù)kp1和kp2,第三個(gè)為T波(熱波)的復(fù)波數(shù)kp3,而ks只有1個(gè)有意義的復(fù)根,即為剪切S波的復(fù)波數(shù)。

2 S波的反射問(wèn)題

2.1 熱彈性波的勢(shì)方程

設(shè)飽和多孔熱彈性地基中頻率為ω的平面S波以任意角度αrs入射至邊界處,在邊界處將產(chǎn)生四種反射波(反射P1波、反射P2波、反射T波和反射S波),如圖1所示。

圖1 平面S波的反射示意圖Fig.1 Reflection diagram of plane S-wave

對(duì)于飽和多孔熱彈性介質(zhì),即z≥0的部分,則有入射S波的勢(shì)函數(shù):

(13a)

反射S波的勢(shì)函數(shù):

(13b)

反射P波的勢(shì)函數(shù):

反射T波的勢(shì)函數(shù):

式中:krs和vrs分別表示入射S波的波數(shù)和波速;kfs和vfs分別表示反射S波的波數(shù)和波速;kfp1、kfp2和kfp3分別表示反射P1波、P2波和T波的波數(shù);vfp1、vfp2和vfp3分別表示反射P1波、P2波和T波的波速。

將式(10)和(11)代入式(12),可以推導(dǎo)出飽和多孔熱彈性介質(zhì)中的P波、T波和S波在不同介質(zhì)中的勢(shì)函數(shù)幅值比例關(guān)系如下:

(14a)

(14b)

式中:γ=1,2,3。

2.2 邊界條件

考慮地表為自由透水、絕熱邊界,則在地表z=0處有以下邊界條件:

(15)

其中

αM▽2ψw-(αauM+Kac)θ

(16a)

(16b)

pw=-αM▽2ψs+M▽2ψw+auMθ

(16c)

將式(13a)～(13d)代入式(16a)～(16c)并結(jié)合邊界條件式(15),略去公共項(xiàng)后可得如下矩陣形式的各幅值之間的關(guān)系式:

(17)

式中:上標(biāo)T表示轉(zhuǎn)置,[D]為4×4階矩陣,[E]為4×1階矩陣,其中的元素分別為:

(αauM+Kac)δT1;

(αauM+Kac)δT2;

(αauM+Kac)δT3;

(18)

式中:nfp1、nfp2、nfT和nfs分別表示反射P1波、反射P2波、反射T波和反射S波的振幅反射率。

3 數(shù)值計(jì)算與討論

表1 介質(zhì)的材料參數(shù)

圖2分別表示在不同入射頻率條件下地基邊界處各反射波的振幅反射率隨平面S波入射角度的變化曲線。圖2表明,入射角的變化對(duì)各波的振幅反射率均有較明顯的影響。當(dāng)平面S波的入射角為0°,即平面S波垂直入射至飽和多孔熱彈性介質(zhì)的邊界時(shí),只存在反射S波,且其振幅反射率為1。隨著入射角的增大,反射P1波、反射P2波和反射T波的振幅反射率均隨之增大,當(dāng)入射角達(dá)到某一臨界角時(shí)其振幅反射率開(kāi)始減小。而反射S波的振幅反射率隨著入射角的增大而減小,且已有研究表明當(dāng)入射角達(dá)到某一值時(shí),其振幅反射率將增大。此外,由圖2可以看出平面S波的入射頻率對(duì)各反射波的振幅反射率亦有較大影響,主要表現(xiàn)為四種反射波的振幅反射率均隨著入射頻率的增大而增大。同時(shí),可以看出在同一入射角和入射頻率下,四種反射波的振幅反射率按由大到小的次序排列依次為反射S波、反射P1波、反射P2波和反射T波。

圖2 不同頻率下各波振幅反射率隨入射角的變化Fig.2 Variation of the amplitude reflection ratio of different reflected wave with incident angle under different frequencies

圖3分別給出了在不同熱膨脹系數(shù)的飽和多孔熱彈性介質(zhì)的邊界上四種反射波隨平面S波入射角度大小的變化曲線。圖3中四種反射波的振幅反射率隨平面S波入射角的變化曲線與圖2計(jì)算結(jié)果保持一致。由圖3可知,飽和多孔介質(zhì)的熱膨脹系數(shù)對(duì)四種反射波的振幅反射率均有不同程度的影響,其中熱膨脹系數(shù)對(duì)反射P1波的影響最為顯著?？梢钥闯?當(dāng)熱膨脹系數(shù)au從2.0×10-4℃-1增大到4.0×10-4℃-1時(shí),四種反射波的振幅反射率均有不同程度的增大,且反射P1波和反射S波的振幅反射率的增大幅度與反射P2波和反射T波的振幅反射率的增大幅度相差多個(gè)數(shù)量級(jí)。這是由于在飽和多孔介質(zhì)中,熱膨脹系數(shù)對(duì)于P1波有著很大影響,而其對(duì)P2波、T波和S波的影響基本可以忽略不計(jì)[23]。

圖3 不同熱膨脹系數(shù)下各波振幅反射率隨入射角的變化Fig.3 Variation of the amplitude reflection ratio of different reflected wave with incident angle under different thermal expansion coefficients

4 結(jié)語(yǔ)

本文基于多孔介質(zhì)理論和廣義熱彈性模型,研究平面S波由飽和多孔熱彈性介質(zhì)向其邊界入射并發(fā)生反射情形下,入射波頻率、土體熱膨脹系數(shù)和入射角對(duì)四種彈性反射波的振幅反射率的影響。通過(guò)數(shù)值算例并結(jié)合參數(shù)分析,得到以下結(jié)論:

(1) 平面S波的入射角對(duì)四種反射波的振幅反射率均有較大影響。當(dāng)平面S波垂直入射時(shí),只存在反射S波。

(2) 平面S波入射頻率的變化對(duì)四種反射波的振幅反射率亦有顯著影響,隨著入射波頻率的增大,四種反射波的振幅反射率亦將增大。

(3) 四種反射波的振幅反射率均隨著飽和土體的熱膨脹系數(shù)的增大而呈不同程度的增大。