層狀巖體濾波特性研究

第一作者王觀石男，博士，副教授，1977年生

層狀巖體濾波特性研究

王觀石1，龍平1，胡世麗1,2

(1.江西理工大學(xué)建筑與測繪工程學(xué)院，江西贛州341000；2.中國地質(zhì)大學(xué)工程技術(shù)學(xué)院，北京100083)

摘要：用振幅、相位透射、反射系數(shù)描述應(yīng)力波在結(jié)構(gòu)面的波形變化，充分考慮應(yīng)力波在層間、層內(nèi)多重透射、反射條件下建立層狀巖體中透射、反射系數(shù)計算模型；分別計算的層狀巖體透射系數(shù)及波形與已有研究結(jié)果吻合，表明該模型計算結(jié)果合理。討論結(jié)構(gòu)面剛度、巖層厚度及結(jié)構(gòu)面數(shù)量對層狀巖體濾波特性影響。結(jié)果表明，結(jié)構(gòu)面剛度及巖層厚度是影響彈性層狀巖體濾波特性的主要因素，層狀巖體具有多尺度梳狀濾波器特性，隨巖層厚度與波長之比γ增加，振幅透、反射系數(shù)成兩種尺度周期性變化，由層內(nèi)多重透、反射形成大尺度周期γ=0.5，并存在一個帶通、一個帶阻，在帶通內(nèi)又存在多個小尺度周期，由層間多重透、反射形成；層狀巖體具有多帶通濾波特性，結(jié)構(gòu)面剛度增加，巖體帶通上、下限頻率同時增加，且上限頻率增加幅度較下限大，隨頻率增加，帶通寬度近似成負(fù)指數(shù)減小；所有巖層厚度相同時帶通中心頻率約為巖層一階固有頻率的整數(shù)倍，各巖層厚度不相等時層狀巖體僅在各巖層共同諧振頻率處表現(xiàn)為帶通。

關(guān)鍵詞：層狀巖體；多重透反射；梳狀濾波器；低通濾波

基金項目：國家自然科學(xué)基金(41462009,51104069)；2014年江西理工大學(xué)科研基金(NSFJ2014-G06)

收稿日期：2014-11-27修改稿收到日期：2015-04-27

中圖分類號:TU443

文獻(xiàn)標(biāo)志碼：A

DOI:10.13465/j.cnki.jvs.2015.20.023

Abstract:The waveform change of stress wave propagating across the interfaces can be described by the transmission and reflection coefficients of its amplitude and phase. The calculation model for the transmission and reflection coefficients of stress wave propagating through stratified rock mass was established on the basis of considering fully the conditions of propagating through multiple interlayers as well as the intrastratal transmissions and reflections. The transmission coefficient and transmission waveform through stratified rock mass were calculated using the established model, and the results are of good agreement with the previous research results. The analysis results show that the model is reasonable. The effects of interface stiffness, thickness of stratum and interface number on the filter property of stratified rock mass were discussed. The study results indicate that interface stiffness and thickness of stratum are the main influential factors on filter property of elastic stratified rock mass. And the stratified rock mass has the multi-scale comb filtering properties. The amplitude transmission and reflection coefficients have two-scale periodic change with the increase of the ratio of stratum thickness to wavelength. The large-scale period, for intrastratal multiple transmissions and reflections of stress wave across stratified rock mass, is 0.5. There exist a pass band and a stop band in a large-scale period. There are as well multiple small-scale periods in a pass band caused by multiple transmissions and reflections between interlayers. The stratified rock mass possesses the multiple band-pass filtering property. The upper limit frequency and lower limit frequency of pass band increase simultaneously with the increased of interface stiffness. The increased extent of upper limit frequency is more than that of lower limit frequency. The width of pass band decreases approximately in a negative exponential form with the increase of frequency. The center frequency of pass band is about an integral multiple of the first order inherent frequency of stratum when each stratal thickness is equal. But, the stratified rock mass has the band-pass property only at its resonance frequency when each stratum has unequal thickness.

Filter property of stratified rock mass

WANGGuan-shi1,LONGPing1,HUShi-li1,2(1. Jiangxi University of Science and Technology, School of Architectural and Surveying Engineering, Ganzhou 341000，China；2. China University of Geosciences, School of Engineering and Technology, Beijing 100083，China)

Key words:stratified rock mass; multiple transmissions and reflections; comb filter; low-pass filtering

層狀巖體在各類工程建設(shè)中經(jīng)常遇到，研究應(yīng)力波在其中的傳播過程及規(guī)律具有重要意義[1-2]。單一結(jié)構(gòu)面透射、反射系數(shù)作為研究應(yīng)力波在層狀巖體中傳播規(guī)律基礎(chǔ)，已建立巖層緊密連接模型、巖層彈簧連接模型、軟弱夾層模型、結(jié)構(gòu)面非線性變形模型及結(jié)構(gòu)面黏彈性變形模型等[3-4]。與單一結(jié)構(gòu)面相比，應(yīng)力波在層狀巖體傳播過程復(fù)雜得多，主要原因為應(yīng)力波會發(fā)生層內(nèi)及層間多重透射、反射，產(chǎn)生大量相干波形，導(dǎo)致波場計算非常復(fù)雜。

針對精確疊加多重透射、反射波計算的困難，已提出等效介質(zhì)法。研究結(jié)果表明，巖層厚度與波長之比遠(yuǎn)小于1時，層狀介質(zhì)波動現(xiàn)象方可用等效介質(zhì)理論描述[5]。而該方法忽略應(yīng)力波在結(jié)構(gòu)面間的多重透射、反射作用，不能完全體現(xiàn)節(jié)理巖體的頻率相關(guān)性，因而在某些頻率范圍內(nèi)對介質(zhì)的動力特性反映不準(zhǔn)確[6-7]，故嘗試建立充分考慮多重透射、反射條件的系數(shù)計算公式。Pyrak-Nolte等[8]利用位移不連續(xù)模型研究應(yīng)力波在多個平行且剛度相等結(jié)構(gòu)面間的傳播過程，認(rèn)為結(jié)構(gòu)面間距與入射波波長相比較大時，可不考慮應(yīng)力波多重透、反射；結(jié)構(gòu)面間距與波長相比較小時多重反、透射影響較顯著，并在不考慮多重透、反射條件下給出層狀巖體總透射系數(shù)計算公式。Cai等[9]在考慮多重反射前提下將位移不連續(xù)模型與特征值方法相結(jié)合，研究應(yīng)力波垂直入射多個平行結(jié)構(gòu)面時的衰減問題。結(jié)果表明，結(jié)構(gòu)面間距與波長之比ζ對層狀巖體總透射系數(shù)具有重要影響，隨ζ增加總透射系數(shù)先增加后減小，最后趨于穩(wěn)定值，且與基于一次透射波計算的透射系數(shù)相等，存在兩個臨界值ζcri1與ζcri2。文獻(xiàn)[10]研究表明，兩個臨界值與K/Zω有關(guān)，K/Zω增大臨界值減小。LI等[11]通過引入“虛振源”概念，討論應(yīng)力波在層狀巖體中的傳播特性，將透射波分為直達(dá)波與多重透射波兩部分，建立等間距平行節(jié)理巖體的透射系數(shù)計算模型。

應(yīng)力波穿過一組平行節(jié)理巖體時波形會發(fā)生明顯變化[12]，數(shù)值計算、理論及試驗研究結(jié)果表明，應(yīng)力波穿過層狀巖體或隨機(jī)介質(zhì)模型，形成特征頻率及特征波形, 不同厚度巖層會對某一頻率的波產(chǎn)生明顯放大效應(yīng)[13]，波場橫向中心頻率、縱向中心頻率及能量相對值均隨介質(zhì)自相關(guān)長度增加而降低，隨機(jī)介質(zhì)模型特征頻率隨自相關(guān)長度倒數(shù)成正比關(guān)系[14]。特征頻率及波形本質(zhì)上一致，前、后者分別在頻、時域上反映波場信息。解釋此現(xiàn)象，不僅需充分考慮應(yīng)力波在層狀巖體中層內(nèi)、層間多重透射、反射，亦需建立可考慮任意巖層厚度、結(jié)構(gòu)面剛度的數(shù)學(xué)模型。

應(yīng)力波穿過層狀巖體總透射波為所有透射波的疊加，每一列透射波不僅受巖層幾何參數(shù)(巖層厚度)、物理力學(xué)參數(shù)(密度、巖層彈模和結(jié)構(gòu)面剛度)及入射波頻率、相位影響，且受多重反射波影響，給出嚴(yán)謹(jǐn)?shù)耐干?、反射系?shù)計算公式有一定難度；但可進(jìn)行不同程度簡化，如假設(shè)巖層為等厚度或巖層物理力學(xué)參數(shù)相等。

本文用振幅透射、反射系數(shù)與相位透射、反射系數(shù)描述應(yīng)力波在結(jié)構(gòu)面的波形變化，將應(yīng)力波穿過含m條平行結(jié)構(gòu)面過程分解為應(yīng)力波下行入射與上行入射(m-1)條結(jié)構(gòu)面及單一結(jié)構(gòu)面兩種，充分考慮應(yīng)力波多重透射、反射，建立水平層狀巖體中透射、反射系數(shù)計算模型，分析結(jié)構(gòu)面剛度、巖層厚度及結(jié)構(gòu)面數(shù)量對層狀巖體透射系數(shù)影響規(guī)律。

1層狀巖體濾波的計算模型

1.1P波穿過線性變形結(jié)構(gòu)面振幅及相位變化

考慮水平層狀巖體，巖層間為薄結(jié)構(gòu)面，建立坐標(biāo)系，z軸向下為正向，坐標(biāo)原點(diǎn)在巖層1、2間結(jié)構(gòu)面1上，隨z軸坐標(biāo)值增加巖層及結(jié)構(gòu)面編號從1依次增加，第1個巖層與最后一個巖層均為半無限體。當(dāng)一列圓頻率為ω的P波垂直入射至層狀巖體，在結(jié)構(gòu)面處產(chǎn)生反射P波及透射P波時，傳播方向與z軸正向一致的稱為下行波，反之為上行波。P波入射至結(jié)構(gòu)面后，設(shè)結(jié)構(gòu)面發(fā)生線性變形，結(jié)構(gòu)面兩側(cè)應(yīng)力連續(xù)，位移不連續(xù)量等于應(yīng)力與節(jié)理剛度之比，不考慮多重透射、反射時P波在結(jié)構(gòu)面i的總透射、反射系數(shù)[15]分別為

(1)

(2)

式中：下標(biāo)i為結(jié)構(gòu)面及巖層編號；T_i為P波在結(jié)構(gòu)面i的總透射系數(shù)；Ri為P波在結(jié)構(gòu)面i的總反射系數(shù)；j為虛數(shù)單位；ω為入射波圓頻率；Ki為結(jié)構(gòu)面i的法向剛度；Zi為第i巖層P波阻抗。

由式(1)知，簡諧P波通過線性變形結(jié)構(gòu)面的總透射系數(shù)為復(fù)數(shù)形式，表明結(jié)構(gòu)面同時影響透射波振幅及相位。為便于分析結(jié)構(gòu)面對透射波振幅譜、相位譜影響，將總透射系數(shù)模定義為振幅透射系數(shù)，將總透射系數(shù)幅角主值定義為相位透射系數(shù)，反射波定義相同。在結(jié)構(gòu)面i處振幅透射、相位透射系數(shù)分別為

(3)

(4)

(5)

(6)

P波在結(jié)構(gòu)面i處的振幅及相位反射系數(shù)為

(7)

(8)

(9)

(10)

由式(3)～式(10)對比看出，對單一線性變形結(jié)構(gòu)面，入射波從結(jié)構(gòu)面任一側(cè)入射時振幅反射系數(shù)與相位透射系數(shù)不發(fā)生變化，但振幅透射系數(shù)與相位反射系數(shù)發(fā)生變化。

1.2層狀巖體透射及反射系數(shù)

在討論及計算公式基礎(chǔ)上，充分考慮多重透射、反射條件下P波穿過層狀巖體的透、反射系數(shù)計算。本文據(jù)各列透、反射波振幅及相位變化，設(shè)已知應(yīng)力波垂直穿過(m-1)條平行結(jié)構(gòu)面的透射、反射系數(shù)，將應(yīng)力波垂直穿過m條平行結(jié)構(gòu)面情況轉(zhuǎn)化為多次穿過(m-1)條結(jié)構(gòu)面與單一結(jié)構(gòu)面組合情況，從而建立應(yīng)力波垂直穿過m條及(m-1)條平行結(jié)構(gòu)面的透射、反射系數(shù)計算模型。

圖1　P波在層狀巖體中傳播過程 Fig.1 P-wave propagation process across stratified rock mass

圖2　下行入射P波時透、反射簡圖 Fig.2 Transmission and reflection diagram of downward incident P-wave

圖3　上行入射P波時透、反射簡圖 Fig.3 Transmission and reflection diagram of upgoing incident P-wave

(11)

設(shè)第i巖層厚度、波速為hi、ci，第i巖層走時i=hi/ci(不考慮第1層巖層走時，1=0)，前(m-1)個巖層走時及τi，不考慮第1及最后1層巖層走時，m≤2時T=0。由分析知，p波在層狀巖體傳播過程中振幅變化由振幅透射系數(shù)計算，相位變化由相位透射系數(shù)計算，并考慮走時引起的相位。設(shè)入射p波振幅為u0、圓頻率為ω、初相位為φ0的正弦波，寫出p2、p6、p10的計算公式為

(12)

(13)

(14)

(15)

(16)

(17)

(18)

(19)

(20)

(21)

(22)

(23)

給定各巖層阻抗、厚度、波速及結(jié)構(gòu)面剛度，由式(16)、(18)、(22)、(24)可計算應(yīng)力波穿過m條結(jié)構(gòu)面的透射、反射系數(shù)，即此4式構(gòu)成層狀巖體中透射、反射系數(shù)計算模型。

2層狀巖體濾波模型驗證

各巖層厚度相等時，定義γ為巖層厚度h與波長λ之比，即γ=h/λ。據(jù)遞推公式，采用歸納法可論證獲得振幅透射系數(shù)關(guān)于γ的周期為0.5(論證過程在此不詳述)。結(jié)構(gòu)面數(shù)量為2時法向剛度均為5.4 GPa/m，巖層密度均為2 700 kg/m3，各巖層P波波速均為4 000 m/s，入射波頻率為50 Hz。改變巖層厚度，用本文遞推公式計算獲得振幅透射系數(shù)與γ的關(guān)系，并與文獻(xiàn)[16]結(jié)果比較，見圖4。由圖4看出，兩種計算結(jié)果吻合較好，表明本文應(yīng)力波穿過層狀巖體計算模型合理。

圖4　γ對振幅透、反射系數(shù)影響 Fig.4 Effect of γ on amplitude transmission and reflection coefficients

一個半正弦波垂直入射至結(jié)構(gòu)面數(shù)量為5的層狀巖體時，入射波頻率為100 Hz，振幅為1.0 m/s，結(jié)構(gòu)面剛度均為3.5 GPa/m，巖層密度為2 650 kg/m3，巖層波速為5 746.2 m/s，巖層厚度為5.75 m。據(jù)本文層狀巖體振幅、相位透射系數(shù)計算獲得半正弦波穿過層狀巖體波形，并與文獻(xiàn)[17]結(jié)果比較，見圖5。由圖5看出，采用兩種方法計算所得波形非常一致，表明本文的遞推公式合理。

圖5　半正弦波穿過層狀巖體波形變化 Fig.5 Waveform change of half sine wave propagation across the stratified rock mass

3層狀巖體濾波特性

3.1層狀巖體多帶通濾波

設(shè)基本參數(shù)為：巖層密度ρ=2 500 kg/m3，巖層P波阻抗Z=8.0 MPa·s/m。當(dāng)入射波頻率為50 Hz、結(jié)構(gòu)面數(shù)量為10時，結(jié)構(gòu)面剛度均為4 GPa/m。改變巖層厚度，計算所得振幅透射系數(shù)與γ關(guān)系見圖6，相位透射系數(shù)見圖7。由兩圖看出，結(jié)構(gòu)面間距與各巖層物理力學(xué)參數(shù)相同時，層狀巖體具有多尺度周期帶通濾波特性，振幅透射、反射系數(shù)隨巖層厚度與波長之比成周期性變化；周期為0.5，在每周期內(nèi)存在一個帶通、一個帶阻，在帶通內(nèi)又存在多個小透射系數(shù)峰并在帶通內(nèi)等間距分布，個數(shù)與巖層數(shù)相等，反射系數(shù)類似。本文計算模型中，設(shè)第1巖層與最后1個巖層均為半無限體，不考慮第1巖層上邊界及最后1個巖層下邊界反射，因此當(dāng)結(jié)構(gòu)面數(shù)量為10時，參與多重透射、反射的巖層數(shù)量為9個，與圖6、圖7 中小透射系數(shù)峰個數(shù)相等。

分析表明，對任意巖層數(shù)量，振幅透射系數(shù)關(guān)于γ的周期為0.5，若考慮層內(nèi)多重透射、反射過程，巖層厚度為半波長及整數(shù)倍時，振幅透射系數(shù)最大，即層內(nèi)透射系數(shù)成周期性變化；若考慮層間多重透射、反射，即各巖層總厚度為半波長及整數(shù)倍時，振幅透射系數(shù)最大，亦成周期性變化。對含10個結(jié)構(gòu)面的層狀巖體，若同時考慮層內(nèi)與層間的多重透射、反射，每一巖層厚度增加半波長的1/9時，各巖層總厚度增加半波長，因?qū)娱g多重透射、反射使振幅透射系數(shù)出現(xiàn)最大值；所有巖層厚度均增加半波長時，層內(nèi)與層間的多重透射、反射使振幅透射系數(shù)出現(xiàn)最大值。此時層狀巖體的振幅透射系數(shù)才出現(xiàn)一個周期。相位透射系數(shù)分析過程類似?？傊?，振幅透射系數(shù)關(guān)于γ的大周期由層內(nèi)多重透射、反射決定，小透射系數(shù)峰形成的小周期由層間多重透射、反射決定。

巖層厚度為2 m時(其它參數(shù)同前)，改變?nèi)肷洳l率，計算獲得振幅透射系數(shù)與頻率的關(guān)系見圖8。由圖8看出，在頻域內(nèi)層狀巖體的帶通區(qū)間分兩部分，即低通區(qū)間與高頻帶通區(qū)間。低通區(qū)間透射系數(shù)大并成非對稱波浪式分布，共9個谷，第2個及其后的帶通組成高頻帶通區(qū)間，該區(qū)間內(nèi)透射系數(shù)成對稱波浪式分布，共9個峰；層狀巖體表現(xiàn)梳狀濾波器，隨頻率增加帶通寬度減小，前3個帶通寬度分別為300 Hz、88 Hz、48 Hz。

圖6 γ對振幅透射、反射系數(shù)影響Fig.6Effectofγonamplitudetransmissionandreflectioncoefficients圖7 γ對相位透射、反射系數(shù)影響Fig.7Effectofγonphasetransmissionandreflectioncoefficients圖8 振幅透射系數(shù)與頻率關(guān)系Fig.8Relationshipbetweenamplitudetransmissioncoefficientandfrequency

3.2結(jié)構(gòu)面剛度對層狀巖體濾波影響

結(jié)構(gòu)面數(shù)量為15、巖層厚度均為2.0 m、巖層密度阻抗為基本參數(shù)、所有結(jié)構(gòu)面剛度相等時，改變結(jié)構(gòu)面剛度，定義中心頻率為帶通上、下限頻率的平均值，計算獲得帶通寬度與結(jié)構(gòu)面剛度及中心頻率關(guān)系見圖9，用相位解纏繞運(yùn)算所得相位透射系數(shù)與結(jié)構(gòu)面剛度、頻率關(guān)系見圖10。由圖9看出，頻率增加帶通寬度近似成負(fù)指數(shù)減小，結(jié)構(gòu)面剛度增加帶通下限頻率略有增加，上限頻率增加明顯，使帶通寬度增加，中心頻率也有所增加；頻率越小帶通寬度及中心頻率增加越明顯，如結(jié)構(gòu)面剛度從2 GPa/m增加到15 GPa/m，第2個帶通下限頻率從800.6 Hz增加到804.1 Hz，而上限頻率從847.2 Hz增加到1 058.5 Hz，下限頻率增加3.5 Hz，上限頻率增加211.3 Hz，帶寬增加207.8 Hz，但第7個帶通下限頻率從4 800.4 Hz增加4 800.8 Hz，而上限頻率從4 808 Hz增加到4 861.7 Hz，帶寬增加53.3 Hz，較第2帶通帶寬增加幅度小得多。由圖10看出，結(jié)構(gòu)面對應(yīng)力波傳播具有延遲作用，其剛度越小相位透射系數(shù)越大，即延遲越大，且隨頻率增加相位延遲量增加。研究結(jié)果表明[18]，隨剛度增加諧振頻率向高頻方向漂移，本文結(jié)果與其一致。

圖9 巖體帶通寬度與結(jié)構(gòu)面剛度及中心頻率關(guān)系Fig.9Relationshipbetweenband-passwidthofrockmassandinterfacestiffnessandcenterfrequency圖10 相位透射系數(shù)與結(jié)構(gòu)面剛度及頻率關(guān)系Fig.10Relationshipbetweenphasetransmissioncoefficientsandinterfacestiffnessandfrequency圖11 結(jié)構(gòu)面剛度漸變時振幅透射系數(shù)與頻率關(guān)系Fig.11Relationshipbetweenamplitudetransmissioncoefficientsandfrequencyforgradualchangesininterfacesstiffness

以結(jié)構(gòu)面剛度遞減為例，研究結(jié)構(gòu)面剛度不相等時層狀巖體中透射系數(shù)變化規(guī)律。結(jié)構(gòu)面數(shù)量為5時其剛度分別為5 GPa/m、4 GPa/m、3 GPa/m、2 GPa/m、1 GPa/m(其它參數(shù)同前)，計算獲得漸變剛度下振幅透射系數(shù)與頻率關(guān)系，并與結(jié)構(gòu)面剛度均為3 GPa/m時的振幅透射系數(shù)比較，見圖11。由圖11看出，結(jié)構(gòu)面剛度逐漸減小時層狀巖體仍表現(xiàn)為多帶通濾波特性，在第2個及以后帶通寬度略有增加，但振幅透射系數(shù)明顯減小，且?guī)▋?nèi)各透射系數(shù)峰的間距不相等。如第3個帶通，結(jié)構(gòu)面剛度均為3 GPa/m時各透射系數(shù)峰值頻率分別為2 402.4 Hz、2 408.8 Hz、2 416.3 Hz、2 422.6 Hz，各透射系數(shù)峰關(guān)于中心頻率成對稱分布，而結(jié)構(gòu)面剛度遞減時帶通內(nèi)各透射系數(shù)峰值頻率分別為2 401.9 Hz、2 407.2 Hz、2 414.8 Hz、2 426.4 Hz，成非等間距分布。

3.3巖層厚度對層狀巖體濾波的影響

結(jié)構(gòu)面數(shù)量為15，剛度為4 GPa/m，巖層密度、阻抗為基本參數(shù)，在巖層厚度相同條件下改變巖層厚度，計算獲得帶通寬度與巖層厚度及中心頻率關(guān)系，見圖12，采用相位解纏繞運(yùn)算獲得相位透射系數(shù)與巖層厚度及頻率關(guān)系，見圖13。由圖12看出，巖層厚度不同帶通中心頻率也不同。當(dāng)巖層厚度為2 m時按參數(shù)計算得巖層一階固有頻率為800 Hz，前7個帶通中心頻率分別為134.8 Hz、845.2 Hz、1 624.6 Hz、2 416.6 Hz、3 212.6 Hz、4 010.1 Hz及4 808.4 Hz，第2及以后帶通約為巖層一階固有頻率的整數(shù)倍，頻率越大帶通中心頻率與該頻率整數(shù)倍誤差越??；由于1 600 Hz為4種巖層厚度的公共諧振頻率，故4種等巖層厚度層狀巖體在1 600 Hz及整數(shù)倍頻率處具有共同帶通；巖層厚度越小諧振頻率越高，帶通個數(shù)越少，破碎巖體濾波能力越強(qiáng)。由圖13看出，巖層厚度越大巖體對應(yīng)力波的相位延遲越大。

對結(jié)構(gòu)面厚度逐漸遞減的巖體，設(shè)應(yīng)力波穿過5層層狀巖體，5巖層厚度分別為5 m、4 m、3 m、2 m、1 m(其它參數(shù)同前)，計算得厚度漸變條件下振幅透射系數(shù)與頻率關(guān)系，并與巖層厚度均為3 m時振幅透射系數(shù)比較，見圖14。由圖14看出，與等厚度層狀巖體相比，巖層厚度逐漸減小時層狀巖體仍表現(xiàn)為多帶通濾波特性，但帶通頻率范圍發(fā)生變化，厚度漸變巖體在539.7～613.8 Hz范圍內(nèi)透射系數(shù)為0，只在所有巖層的公共諧振頻率1600 Hz及整數(shù)倍頻率處振幅透射系數(shù)不為0，因此巖層厚度分布對巖體濾波特性具有重要影響；巖層厚度漸變時帶通數(shù)減少，帶通寬度減小，帶通內(nèi)透射系數(shù)減小，各透射系數(shù)峰值頻率成非等間距分布。

圖12 帶通寬度與巖層厚度及中心頻率關(guān)系Fig.12Relationshipbetweenband-passwidthandstratumthicknessandcenterfrequency圖13 相位透射系數(shù)與巖層厚度及頻率關(guān)系Fig.13Relationshipbetweenphasetransmissioncoefficientsandstratumthicknessandfrequency圖14 結(jié)構(gòu)面厚度漸變時振幅透射系數(shù)與頻率關(guān)系Fig.14Relationshipbetweenamplitudetransmissioncoefficientandfrequencyforgradualchangesininterfacesthickness

3.4結(jié)構(gòu)面數(shù)量對層狀巖體濾波的影響

結(jié)構(gòu)面數(shù)量不同，層狀巖體帶通濾波特性會發(fā)生變化。結(jié)構(gòu)面剛度為4 GPa/m、巖層厚度2 m時改變結(jié)構(gòu)面數(shù)量，計算得巖體帶通寬度與結(jié)構(gòu)面數(shù)量及中心頻率關(guān)系，見圖15。圖15表明，結(jié)構(gòu)面數(shù)量增加帶通寬度增加，但增加幅度較小。隨頻率增加帶通寬度增加幅度減小。如結(jié)構(gòu)面由5增加到15時，第2帶通下限頻率從808 Hz減小到801 Hz，而上限頻率從883增加到889 Hz，帶寬增加13 Hz；第8帶通下限頻率從5 601 Hz減小到5 600 Hz，上限頻率從5 607 Hz增加到5614 Hz，帶寬增加2 Hz，較第2帶通，該帶寬增加幅度小得多。

圖15　巖體帶通寬度與結(jié)構(gòu)面數(shù)量及中心頻率關(guān)系 Fig.15 Relationship between band-pass width and interface number and center frequency

4結(jié)論

用振幅、相位的透射、反射系數(shù)描述應(yīng)力波在結(jié)構(gòu)面的波形變化，建立應(yīng)力波垂直穿過m條平行結(jié)構(gòu)面透射、反射系數(shù)計算模型，計算層狀巖體透射系數(shù)及波形，討論結(jié)構(gòu)面剛度、巖層厚度及結(jié)構(gòu)面數(shù)量對層狀巖體濾波特性影響，結(jié)論如下：

(1)層狀巖體具有多尺度梳狀濾波器特性，振幅透、反射系數(shù)隨巖層厚度與波長之比成兩種尺度周期性變化。由層內(nèi)多重透射、反射形成的大尺度周期為0.5，其內(nèi)存在帶通、帶阻各一。帶通內(nèi)存在的多個透射系數(shù)峰由層間多重透射、反射形成，對等厚度、等剛度層狀巖體，透射系數(shù)峰在帶通內(nèi)成等間距分布，個數(shù)與巖層數(shù)相等。

(2)頻域上層狀巖體具有多帶通濾波特性，結(jié)構(gòu)面剛度增加，帶通上、下限頻率同時增加，且上限頻率增幅較下限大，致帶通寬度及中心頻率同時增加，頻率越高增幅越小；頻率增加帶通寬度近似成負(fù)指數(shù)減??；巖體結(jié)構(gòu)面剛度相等時帶通內(nèi)各透射系數(shù)峰近似等間距分布，不相等時該系數(shù)峰成非等間距分布，且數(shù)值減小。

(3)所有巖層厚度相等時帶通中心頻率約為巖層一階固有頻率的整數(shù)倍，頻率越高兩者誤差越小；巖層厚度越小其諧振頻率越高，巖體帶通數(shù)越少，巖體濾波能力越強(qiáng)；各巖層厚度不相等時，層狀巖體僅在各巖層的共同諧振頻率處表現(xiàn)為帶通，即巖層厚度分布對巖體濾波具有重要影響。

(4)結(jié)構(gòu)面數(shù)量增加，帶通上限頻率小幅增加，下限頻率小幅減小，各帶通寬度小幅度增加，即結(jié)構(gòu)面數(shù)量對彈性層狀巖體影響較小。

參考文獻(xiàn)

[1]黃書嶺,徐勁松,丁秀麗,等. 考慮結(jié)構(gòu)面特性的層狀巖體復(fù)合材料模型與應(yīng)用研究[J]. 巖石力學(xué)與工程學(xué)報，2010,29(4):743-756.

HUANG Shu-ling，XU Jin-song，DING Xiu-li，et al. Study of layered rock mass composite model based on characteristics of structural plane and its application[J]. Chinese Journal of Rock Mechanics and Engineering，2010,29(4):743-756.

[2]程志寶,石志飛,向宏軍. 層狀周期結(jié)構(gòu)動力衰減域特性研究[J].振動與沖擊,2013,32(9):178-182.

CHENG Zhi-bao,SHI Zhi-fei,XIANG Hong-jun. Vibration attenuation zones of a laminated periodic structure[J]. Journal of Vibration and Shock,2013,32(9):178-182.

[3]Pyrak-Nolte L J，Myer L R，Cook N G W. Transmission of seismic waves across single natural fractures[J]. Journal of Geophysical Research: Solid Earth,1990,95(B6):8617-8638.

[4]Carcione J M. Reflection and transmission of qP-qS plane waves at a plane boundary between viscoelastic transversely isotropic media[J]. Geophysical Journal International,1997,129(3):669-680.

[5]Marion D,Mukerji T,Mavko G. Scale effects on velocity dispersion: from ray to effective medium theories in stratified media[J]. Geophysics,1994,59(4):1613-1619.

[6]房學(xué)謙,胡超,黃文虎. 非均勻復(fù)合材料板中剪切波傳播的研究[J].工程力學(xué),2008,25(2):66-70.

FANG Xue-qian，HU Chao，HUANG Wen-hu. Propagation of shear waves in a heterogeneous material plate[J]. Engineering Mechanics,2008,25(2):66-70.

[7]Frazer L N. SH propagation in rocks with planar fractures[J]. Geophysical Journal International，1995，122(1)：33-62.

[8]Pyrak-Nolte L J，Myer L R，Cook N G W. Anisotropy in Seismic velocities and amplitudes from multiple parallel fractures[J]. Journal of Geophysical Reseach,1990,95(11)：345-358.

[9]Cai J G，Zhao J. Effects of multiple parallel fractures on apparent attenuation of stress waves in rock masses[J]. International Journal of Rock Mechanics and Mining Sciences,2000,37(1):661-682.

[10]胡世麗,羅嗣海,王觀石. 結(jié)構(gòu)面間距對透射波的影響[J]. 礦冶工程,2011,31 (4):5-9.

HU Shi-li，LUO Si-hai，WANG Guan-shi. Influence of interface interval on transmission wave[J]. Mining and Metallurgical Engineering,2011,31 (4):5-9.

[11]Li J, Ma G, Zhao J. An equivalent viscoelastic model for rock mass with parallel joints[J]. Journal of Geophysical Research,2010, 115(B3):1-10.

[12]俞縉,宋博學(xué),錢七虎. 含單組節(jié)理的雙重非線性彈性巖石介質(zhì)中的P波傳播規(guī)律[J]. 巖石力學(xué)與工程學(xué)報，2012,31(12):2400-2451.

YU Jin,SONG Bo-xue，QIAN Qi-hu. Study of propagation of P-waves in dual nonlinear elastic rock medium with one set of joints[J]. Chinese Journal of Rock Mechanics and Engineering,2012,31(12):2400-2451.

[13]言志信,張森,張學(xué)東,等. 順層巖質(zhì)邊坡地震動力響應(yīng)及地震動參數(shù)影響研究[J]. 巖石力學(xué)與工程學(xué)報,2011, 30(增2):3522-3528.

YAN Zhi-xin，ZHANG Sen，ZHANG Xue-dong，et al. Study of dynamic response of bedding rock slope under earthquake and influence of ground motion parameters[J]. Chinese Journal of Rock Mechanics and Engineering，2011,30(S2):3522-3528.

[14]奚先,姚姚. 二維彈性隨機(jī)介質(zhì)中的波場特征[J]. 地球物理學(xué)進(jìn)展,2005,20(1):147-154.

XI Xian,YAO Yao. Wave field characters of 2-D elastic random medium[J]. Progress in Geophysics，2005，20(1)：147-154.

[15]Pyrak-Nolte L J. Seismic response of fractures and the interrelations among fractures[J]. International Journal of Rock Mechanics and Mining Sciences,1996,33(8):787-802.

[16]Perino A. Wave propagation through discontinuous media in rock engineering[D]. Politecnico Di Torino,2011.

[17]Li J C，Li H B，Ma G W，et al. A time-domain recursive method to analyse transient wave propagation across rock joints[J]. Geophysical Journal International，2012，188(2)：631-644.

[18]吳斌,邱兆國,何存富. 基于弱粘接結(jié)構(gòu)中諧振頻率的研究[J]. 應(yīng)用力學(xué)學(xué)報,2011,28(4):412-415.

WU Bin, QIU Zhao-guo, HE Cun-fu. Research on the resonance frequency in weak bonded structures[J]. Chinese Journal of Applied Mechanics,2011,28(4):412- 415.