基于小波變換時能密度法的隧道空洞充填物識別

張 勝，何文超，黎永索，胡 達，蔡 鑫

(1.湖南城市學院 土木工程學院，湖南 益陽 413000; 2.休斯敦大學 機械工程系，德克薩斯州 休斯敦 77204; 3.長沙理工大學 土木工程學院，湖南 長沙 410114)

近年來，隨著國家基礎設施建設的大力發展和高新技術的不斷應用，全國范圍內交通工程建設呈現一派繁榮景象。隧道工程因其具有改善線形、克服高程或地形障礙、保護生態環境等優勢，已成為交通工程建設中最重要的控制性工程之一。截至2018年底，中國隧道建設包括鐵路隧道、公路隧道以及城市軌道交通隧道共計36 103 km，正在施工的各類隧道約20 000 km，計劃修建的隧道約20 000 km[1]。目前，我國在隧道工程規模、工程質量與修建技術等方面已成為世界上發展最為迅速的國家。盡管隧道工程建設發展迅速，但在隧道修建過程中仍存在許多難題，其中，行業內公認的隧道建設難點之一是隧道掘進工作面前方的不良地質條件。由于修建隧道所處位置的地質情況較為復雜且地形起伏較大，在修建前期，限于現有的地質勘察技術手段，難以對工程區域內復雜地形地貌進行準確且全面地掌握，且隧道沿線出現的一些不良地質情況(如溶洞、斷層與破碎帶等)對后期的修建、維護與運營帶來了巨大的安全隱患[2-4]。

隨著國內外學者在地質探測領域的技術突破，使得無損探測技術發展迅速且廣泛應用于實際工程。隧道施工超前地質預報方法有地質分析法、超前鉆探法、地震反射波法、地質雷達法、瞬變電磁法、陸地聲納法、地震波層析成像法等，每種方法都有各自的優、缺點及適用范圍[5-8]。地質雷達因其具有分辨率高、圖像直觀以及對施工干擾少等優點，被廣泛應用于隧道短距離超前地質預報中(預測掘進工作面前方30 m左右)[9-10]。地質雷達信號特征分析是研究地質雷達電磁波在介質中傳播規律的前提與基礎，也是對地質雷達信號圖譜特征進行定量解釋的最有效手段。地質雷達信號分析方法普遍采用基于傅里葉理論的儀器自帶分析軟件，而傅里葉理論沒有時頻局部化分析能力，無法對隧道掘進工作面前方的地質情況進行精確定位。近年來，隨著科學技術的不斷進步，小波變換、S變換、復信號分析以及希爾伯特-黃變換等技術在地質雷達信號閾值去噪與特征提取等方面具備了良好的信號解譯能力[11-15]。由于小波變換具有許多優良的特性，并且理論比較成熟，因而在工程領域應用最為廣泛。然而，采用經典小波對地質雷達信號進行小波時頻局部化分析時，選取何種尺度的小波基通常依賴于技術人員的經驗，雖然由經驗選取的小波基能夠對地質雷達信號進行一般化處理，但是小波基天然無法準確選取的事實，阻礙了小波時頻局部化分析技術在地質雷達信號處理中的應用。因此，如何利用最新的信號處理方法，構造與地質雷達信號特征相似度高的最優小波基，將是小波理論運用于地質雷達信號處理中亟需突破的難題。

筆者以此為切入點，在小波變換與奇異性檢測原理的基礎上，構造與地質雷達發射子波波形契合度高的小波基，然后將該小波基添加到MATLAB小波分析工具箱中，提出一種新的基于雷達小波基的小波變換時能密度法，將其應用于隧道空洞充填物地質雷達正演模擬與室內模型實驗采集信號的定量識別與分析，以期為隧道短距離超前地質預報定量識別不良地質體提供技術支持。

1 識別分析方法

1.1 波形分析法

地質雷達超前地質預報的工作原理是利用地質雷達發射天線發射高頻電磁波以脈沖的形式傳播至隧道掘進工作面前方的巖土介質中，遇到溶洞、斷層與破碎帶等界面時會發生反射，反射波由地質雷達接收天線所接收，最后將在地質雷達主機顯示器上獲得時-距剖面圖像。通過對地質雷達時-距剖面進行處理、分析與解釋，達到對隧道掘進工作面前方的地質情況進行短距離超前地質預報的目的[8-10]。

地質雷達回波反射信號可以表示為

x(t)=A(t)cos[2πf0t+θ(t)]

(1)

式中，x(t)為地質雷達回波反射信號;A(t)為振幅函數;t為記錄時間;f0為中心頻率;θ(t)為相位函數。

地質雷達發射的脈沖電磁波遇到巖土層中存在2種不同介質的分界面時，由于相對介電常數的差異，將在界面上產生反射現象，其規律符合反射定律，產生的電磁波能量大小取決于反射系數Γ。

(2)

式中，ε1，ε2分別為分界面上、下介質的相對介電常數。由式(2)可知，隧道圍巖與充填物的相對介電常數的差異，直接影響了地質雷達圖像的振幅大小、頻率高低與相位起伏，這為地質雷達信號特征解譯提供了良好的理論基礎。

1.2 小波變換模極大值法

小波變換模極大值法最突出的特征是選取信號的奇異點。信號奇異點與突變部分往往帶有明顯的個性化特征，它是一個信號區別于其他信號的重要特征之一。利用小波變換方法分析局部奇異性特征時，小波變換分解系數取決于信號在某一點相鄰區域內的特征信號以及小波變換所選用的尺度。將任意信號f(t)的連續小波變換表示為Wf(a,b)[16]，若對b0的某一鄰域內的任意一點b，有

|Wf(a0,b)|≤|Wf(a0,b0)|

(3)

則稱(a0,b0)為f(t)在尺度a0上的小波變換模極大值點，|Wf(a0,b0)|為f(t)在尺度a0上的小波變換模極大值。

如果信號在某點或某區間內可微的次數越高，那么該信號越平滑，奇異性越弱。在數字信號處理中，常用李氏指數來表征信號的局部奇異性特征。信號f(t)在區間(t1,t2)上具有一致李氏指數α的充要條件為:存在常數k>0，使得對所有的t∈(t1,t2)，信號f(t)的小波變換Wf(a,b)滿足如下不等式:

|Wf(a,b)|≤kaα

(4)

由式(3)～(4)可以得出，信號產生突變的特征點就是小波變換的模極大值點。由于不同尺度下的小波變換模極大值，具有突出不同時頻局部特征信號的能力。因此，采用小波變換模極大值法對信號奇異性特征進行分析時，最重要的前提是選取一個最佳的尺度，在該尺度下才能最大限度顯示信號的奇異性與突變成分。

1.3 小波基的構造與添加

小波基函數的選取在小波變換分析中具有不唯一性。采用不同小波基分析同一信號，所得結果也會各有差異。構造適合被分析信號特征選用的最優小波基，是目前小波分析中亟需解決的難題。從理論上講，構造一個新的小波基的前提是只要滿足小波基的允許條件。因此，可根據地質雷達發射子波f(t)=t2e-atsinω0t的特點，其中，a為衰減系數；ω0為中心角頻率，由發射子波來進行自適應波形匹配，構造與地質雷達發射子波波形相似度高的自適應雷達小波基，將該小波基添加到小波分析工具箱中，以期對地質雷達信號進行小波變換時作為小波基進行選用。

1.4 小波變換時能密度法

設信號f(t)連續小波變換系數的幅度平方在時間尺度平面上的加權積分與信號的能量成比例關系，則滿足下列能量守恒關系[17-18]，即

(5)

式中，Cψ為小波基容許條件;R為實數。

由式(5)可知，小波變換幅度平方的積分同被分析信號的能量成正比。在分析與處理非平穩信號時，由于海森堡(Heisenberg)測不準原理的限制，在時-頻相空間中某一點的瞬時能量密度無法得知。因此，在某一指定時刻瞬時頻率的說法是不存在的。如果把|Wf(a,b)|2/Cψa2認為是時間-尺度平面上的能量密度函數，那么可將|Wf(a,b)|2ΔaΔb/Cψa2看作是以尺度a和時間b為中心、尺度間隔為Δa、時間間隔為Δb的能量。根據能量的概念，式(5)可以改寫成

(6)

式中，

(7)

小波變換中，尺度因子a在某種意義上對應于頻率。因此，式(7)表示信號所有頻帶的能量隨時間b的分布情況，稱為時能密度函數。

2 空洞充填物的正演模擬

2.1 時域有限差分法

地質雷達正演模擬是地質雷達圖像解譯的基礎，解譯人員要想對地質雷達信號進行清楚掌握與定量解釋，就必須事先了解電磁波在隧道掘進工作面前方的傳播規律與圖譜特征。時域有限差分法以其具有存儲空間小與計算效率高等優勢，成為了地質雷達正演模擬的最主要方法[19-21]。

在無源場區域內，Maxwell方程的2個旋度可以表示為

(8)

式中，H為磁場強度，A/m;E為電場強度，V/m;ε為媒質的介電常數;σ為電導率，S/m;t為時間s;μ為相對磁導率，H/m;σm為等效磁導率，w/m。

時域有限差分法利用二階精度的中心差分形式，將Maxwell方程中的2個旋度由微分轉成差分，電場與磁場在時間順序上交替抽樣，彼此相差半個時間步長。因此，二維電磁波的時域有限差分方程可以表示為

(9)

(10)

(11)

式中，c為電磁波在真空中的傳播速度，c=3×108m/s。

2.2 空洞充填物的正演模擬

假設隧道掘進工作面前方為半無限空間連續的均勻介質，電磁波反射與折射均在二維平面內進行。為了研究隧道空洞內充填不同介質的地質雷達正演模擬圖像的反射特征，設計了如圖1所示的隧道空洞充填物的地電模型。

圖1 隧道空洞充填物的地電模型

隧道掘進工作面前方空洞充填物的地電模型參數設置為:① 區域范圍:12 m×20 m，左上角為坐標原點，橫坐標為隧道掘進工作面水平距離，縱坐標為探測深度;② 目標體:長方形空洞，尺寸為2 m×1 m，空洞左側與模擬區域邊緣相距5 m，埋深為2 m;③ 圍巖灰巖的相對介電常數為6，電導率為0.001 S/m，導磁率為1;④ 空洞充填物分別為空氣、冰、干黏土與濕黏土。表1給出了隧道掘進工作面前方空洞內不同充填物的幾何與物理參數。

表1 隧道空洞充填物的幾何與物理參數

Table 1 Geometric and physical parameters of tunnel cavity fillings

充填物介電常數電導率/(S·m-1)導磁率電磁波速度/(m·ns-1)空氣1010.3000冰31.4×10-410.1732干黏土81×10-3～0.1110.1061濕黏土121×10-3～0.1110.0866

2.3 地質雷達正演模擬特征分析

根據隧道超前地質預報實際檢測所使用的參數，模擬時采用的中心頻率為100 MHz，邊界吸收條件為完全匹配層，激勵源采用Ricker子波，網格的空間步長為0.05 m，采樣步長為0.05 m，采樣道數為227，總采樣時間為200 ns。

由于隧道開挖過程中，特別是灰巖地區，隧道掘進工作面前方經常會遇到空洞，且空洞內含有不同的充填物。為了更好反映地質雷達正演模擬的實際情況，本次模擬將在空洞中分別填充空氣(ε=1)、冰(ε=3且ε<6)、干黏土(ε=8且ε>6)與濕黏土(ε=12且ε>6)。采用100 MHz頻率天線對如圖1所示的地電模型空洞內充填不同介質進行正演模擬，得到正演模擬雷達響應圖像，并將其減去完整灰巖的地質雷達正演模擬圖像，結果如圖2所示。處理過程中對每種工況進行標記，如“y1”對應充填物為空氣，“y2”對應充填物為冰，“y3”對應充填物為干黏土，“y4”對應充填物為濕黏土。從圖2可以看出:① 地質雷達電磁波遇到充填物空氣、冰、干黏土與濕黏土等不同介質時，界面處均存在反射，空洞充填物與灰巖的相對介電常數差異較小時，電磁波回波信號的反射不明顯;② 電磁波由空氣進入灰巖，經灰巖再進入充填物，最后進入灰巖，各界面反射波依次被接收天線所接收。當空洞充填物的相對介電常數小于或大于灰巖的相對介電常數時，電磁波會在空洞充填物的第一或第二界面的界面處發生相位反相。而從圖2所示的地質雷達時間剖面上無法識別電磁波遇到界面產生反射波的瞬時位置(究竟處于波峰還是波谷)。因此，需采用其他方法作進一步分析。

圖2 隧道空洞充填物的地質雷達信號時間剖面

3 空洞充填物正演模擬的定量識別

3.1 波形分析法的定量識別

由圖3可知，根據空洞充填物與灰巖相對介電常數的關系，通過計算反射系數來確定反射界面位于反射波波峰還是波谷，從而得到電磁波在不同空洞充填物上的第一、二反射界面，分別如圖3所示的a1，a2，b1，b2，c1，c2和d1，d2。依次計算不同空洞充填物單道信號兩特征點的時間差，計算結果分別為11.21，11.21，19.23和23.59 ns。也就是說，空洞充填物的厚度分別為1.68，0.97，1.02和1.02 m。采用波形分析法識別不同空洞充填物的相對誤差分別為68%，-3%，2%和2%，說明后三者的識別精度高。但須事先知道空洞充填物與灰巖相對介電常數的關系，即電磁波遇到異種界面時瞬時反射的位置究竟是反射波波峰還是波谷的位置。否則，難以準確確定反射界面在單道信號上的兩特征點。

3.2 小波變換模極大值法的定量識別

小波變換模極大值法是一種經典小波去噪方法，信號模極大值點的位置對應于信號奇異點。從理論上講，小波變換尺度選取得越小，模極大值點與信號奇異點的對應位置越準確，但小尺度下的小波變換系數易受噪聲的干擾，產生偽奇異點。尺度較大時，可使噪聲相對平滑，模極大值點位置易于尋找，但會使模極大值點位置產生一定的偏移。因此，在采用小波變換模極大值法識別信號的奇異點時，應選取適宜的尺度以避免小波變換所產生的交迭干擾。

2）本實驗中的混合菌群在原油降解的前期對中鏈、長鏈烴降解效果較好；而在降解的后期對短鏈烴的降解效果較強.

大量試驗與研究結果表明，離散小波變換選取合適小波基時，Daubechies系列小波基因其具有良好的緊支撐性、光滑性與近似對稱性等優點被廣泛應用于地質雷達信號的處理與分析中[22-25]。經過比較Daubechies系列小波變換模極大值法提取信號奇異點的效果后，選取Db4小波作為小波變換模極大值法提取地質雷達信號奇異點的小波基。通過分別對圖3所示的地質雷達單道信號進行Db4小波變換模極大值法分析，得到如圖4所示結果。

圖3 隧道空洞充填物的地質雷達單道信號

圖4 地質雷達單道信號Db4小波變換模值(a=10)

由圖4可知，空洞不同充填物的第一、二反射界面在地質雷達正演模擬單道信號Db4小波變換模極大值曲線上的兩特征點分別為a1，a2，b1，b2，c1，c2和d1，d2。兩特征點的時間差分別為8.26,13.45，21.81和26.41 ns，即空洞充填物的厚度分別為1.24,1.16,1.16和1.14 m，則識別空洞充填物的相對誤差分別為24%，16%，16%和14%，相對誤差均大于10%。由此說明，小波分析工具箱中已有的小波基分析地質雷達信號誤差較大，亟需構造新的小波基以提高地質雷達探測深度的準確率。

3.3 雷達小波變換模極大值法的定量識別

選擇合適尺度的小波基是對地質雷達信號進行小波分析與處理的重要環節，其原因是不同小波基分析同一信號會產生不同的結果。在小波分析工具箱中擁有眾多小波基可供挑選，但不能保證找到與地質雷達特征信號相吻合的小波基。因此，需構造與地質雷達特征信號相似度高的自適應雷達小波基，將其添加到小波分析工具箱中，以供進行小波變換時選用。

根據前文1.3節所列方法構造雷達小波基，對空洞不同充填物的地質雷達正演模擬單道信號進行雷達小波變換模極大值法分析，結果如圖5所示。

圖5 地質雷達單道信號雷達小波變換模值(a=10)

由圖5可知，空洞不同充填物的上、下反射界面在地質雷達正演模擬單道信號雷達小波變換模極大值曲線上的兩特征點分別為a1，a2，b1，b2，c1，c2和d1，d2。根據兩特征點的時間差8.61，11.09,19.22和23.58 ns，可以得到空洞不同充填物的厚度分別為1.29，0.96，1.02和1.02 m，即識別結果的相對誤差分別為29%，-4%，2%和2%。

3.4 小波變換時能密度法的定量識別

采用已構造的雷達小波基作為小波變換用的小波基，在MATLAB語言平臺工作環境下編寫小波變換時能密度分析程序，對空洞不同充填物的地質雷達正演模擬單道信號進行小波變換時能密度法分析，結果如圖6所示。

由圖6可知，空洞不同充填物的上、下反射界面在地質雷達正演模擬單道信號小波變換時能密度曲線上的兩特征點分別為a1，a2，b1，b2，c1，c2和d1，d2。通過計算兩特征點的時間間隔，結果依次為8.49，11.09，19.22和23.58 ns，即空洞充填物的厚度分別為1.27，0.96，1.02和1.02 m，由此可以說明識別結果的相對誤差分別為29%，-4%，2%和2%。

為了比較波形分析法、Db4小波變換模極大值法、雷達小波變換模極大值法和小波變換時能密度法的識別效果，將不同方法得到的測量值與相對誤差匯總于表2。

由表2可知，當空洞充填物為空氣時，不同方法得到的識別結果相對誤差都比較大，依據電磁波在空氣中的傳播速度0.3 m/ns，可以得到電磁波在空氣中的波長為3 m。由于充填物的厚度僅為1 m，可能發生波形混疊，導致無法識別空洞充填物的第二界面;空洞充填物為冰、干黏土與濕黏土時，電磁波在相應介質中的波長分別為1.73，1.06與0.87 m。波形分析法、雷達小波變換模極大值法與小波變換時能密度法的識別效果均較好，但波形分析法需預先了解圍巖前方的介質屬性，通過計算電磁波在界面上反射系數的正負，才能確定空洞充填物第一、二反射界面在電磁波反射信號上的具體位置;通過對圖5與圖6的比較不難看出，小波變換時能密度法的分辨率要比雷達小波變換模極大值法的高，表明小波變換時能密度法具有更好突出地質雷達信號奇異點的位置，從而驗證了將地質雷達發射子波作為地質雷達信號小波分析用的小波基是切實可行的，解決了適合地質雷達信號特征的小波基構造、添加與實現等問題，為分析地質雷達信號小波變換時選擇小波基提出了一種新的方法。

圖6 地質雷達單道信號小波變換時能密度曲線

表2 不同方法識別空洞充填物的測量值與相對誤差

4 實驗研究

4.1 方案設計

在眾多研究方法中，室內模型實驗是最直觀與最可靠的研究手段之一。隧道空洞充填物模型實驗中，模型尺寸與材料是影響室內模型實驗結果好壞的重要因素。根據已有相關文獻與資料，隧道模型尺寸相似比應盡量控制在15～30。采用地質雷達100 M天線對隧道掘進工作面前方地質情況進行短距離超前地質預報預測的探測深度一般為30 m，本次實驗模型箱相似比(深度方向)取為20，即模型箱的高度為1.5 m。地質雷達電磁波在空氣中的傳播速度為0.3 m/ns，即波長為3 m。根據模型尺寸相似比可將空洞的尺寸設置為10～20 cm，本次實驗空箱的尺寸取為10 cm，除去壁厚與加工誤差，空箱的實際尺寸為9.25 cm。地質雷達用于識別隧道掘進工作面前方空洞位置的精確程度主要取決于電磁波的傳播速度，即圍巖的介電常數。因此，應優先選擇與灰巖相對介電常數相同或相近的介質作為模型的實驗材料。干砂的相對介電常數為4～6，與灰巖的介電常數較為接近，具有重復使用、試制材料垃圾少以及成本低等優點，本次模型背景材料選擇干砂。為便于砂箱物理模型成型，砂箱四周采用混凝土圍制而成，長、寬和高分別設置為4，3和1.5 m，以盡可能減少混凝土邊界對地質雷達采集信號的干擾，具體如圖7所示。砂箱中心埋設空箱模擬空洞，空洞充填物分別為空氣、冰、干黏土與濕黏土的實物圖如圖8所示。

4.2 地質雷達信號的采集

由于空箱的實際尺寸為9.25 cm，在空洞充填物地質雷達檢測的物理模型實驗中，采用意大利RIS地質雷達，配以1 600 MHz天線。RIS地質雷達具有輕質便攜、天線屏蔽效果好與抗干擾能力強等優點。主機技術參數如下:掃描速度4 761掃/s，脈沖重復頻率400 kHz，采樣點數128～8 192，A/D轉換16 bit，疊加數1～32 768，分辨率5 psec，動態范圍>160 dB，信噪比>160 dB。按照奈奎斯特采樣定理，采樣頻率至少應為反射波最高頻率的2倍，但在頻率比僅僅為2時，雷達信號失真很明顯。地質雷達在實際應用過程中，為了達到良好的探測效果，其采樣頻率應為天線主頻的10倍以上，本次使用1 600 MHz天線所設采樣頻率為32 GHz。根據天線的中心頻率，將實驗記錄時間設為16 ns和疊加次數為512。

圖7 隧道空洞充填物的模型箱

圖8 隧道空洞不同充填物的實物

依據空洞充填物埋設在模型中的具體位置已知，地質雷達檢測時沿一條直線勻速前進。地質雷達數據采集過程中對每種工況進行標記，如“y5”對應充填物為空氣，“y6”為冰，“y7”為干黏土，“y8”為濕黏土。采用地質雷達所攜帶的分析軟件對采集的原始信號進行靜校正、去直流點漂移、增益以及帶通濾波等處理，得到空洞不同充填物的實測地質雷達信號時間剖面如圖9所示。

圖9 隧道空洞充填物的實測地質雷達信號時間剖面

由圖9可知，地質雷達電磁波遇到空洞、冰、干黏土與濕黏土等介質時，界面處均存在反射，但不能顯示具體的量值，尤其是空洞充填物的第一界面位于電磁波波形的波峰還是波谷無法得知。因此，需采用其他方法作進一步分析。

4.3 空洞充填物的定量識別

從如圖9所示的隧道空洞充填物實測地質雷達信號時間剖面中抽取最中間的1條信號，獲得1條反映空洞充填物典型特征的單道信號。采用前文已構造的雷達小波作為小波變換用的小波基，在MATLAB語言環境下運行小波變換時能密度法分析程序，對空洞不同充填物單道信號進行小波變換時能密度法分析，結果如圖10所示。

從圖10可以看出，空洞不同充填物的上、下反射界面在實測單道信號小波變換時能密度曲線上的兩特征點分別為e1，e2，f1，f2，g1，g2與h1,h2。通過計算兩特征點的時間間隔(分別0.908 4，1.065 0，1.911 0與2.412 1 ns)，可以得到空洞充填物的厚度分別為13.63，9.36，9.87與8.92 cm，則識別結果的相對誤差為47.35%,1.19%，3.35%與-3.57%。盡管采用地質雷達1 600 MHz高頻天線(探測深度淺、精度高)，但仍無法識別10 cm左右空洞，究其原因是電磁波在空氣中的波長為18.75 cm，大于空洞的實際尺寸。采用小波變換時能密度法識別空洞充填物為冰、干黏土與濕黏土的精度高，相對誤差均小于5%。由此表明，小波變換時能密度法能成功識別空洞充填物第一、二界面地質雷達反射信號奇異點的位置，為隧道施工現場超前地質預報圖像定量識別提供技術支持。

圖10 實測單道信號小波變換時能密度曲線

5 結 論

(1)波形分析法雖能有效識別空洞充填物的尺寸大小，但需通過計算反射系數的正或負值等先驗知識來確定空洞充填物的界面特征點處于波峰還是波谷，由于空洞充填物的介電常數存在大于或小于圍巖介電常數的情況，導致判斷過程較繁瑣，因此不適用于隧道超前地質預報精細化定量檢測。

(2)盡管小波變換模極大值法容易得到信號突變的模極大值點，但選取不同小波基，可能會產生截然不同的結果。通過地質雷達正演模擬分析結果得知，小波變換模極大值法識別結果的誤差約15%，誤差相對較大，其原因是小波基的時頻特征與被分析信號時頻特征的相似性不高。

(3)在小波變換原理的基礎上，從與地質雷達反射子波波形相匹配出發，成功構造了地質雷達小波基。雷達小波變換模極大值法與小波變換時能密度法的識別效果均較好，但小波變換時能密度法的分辨率更高，壓制隨機干擾的能力更強，更重要的是不需優選最優尺度。當空洞充填物的尺寸大于電磁波波長時，小波變換時能密度法識別結果的相對誤差均小于5%。