基于計算機圖像處理技術的斷級配瀝青混合料三維空隙率計算

王 雙，Mischelle Esguerra

(1.哈爾濱劍橋學院 電氣與信息工程學院 哈爾濱市 150000； 2.菲律賓萊西姆大學 信息技術學院，馬尼拉，4200)

0 引言

斷級配瀝青混合料主要用于瀝青路面磨耗層，可分為多空隙開級配排水式瀝青混合料(OGFC)與密實式斷級配瀝青混合料(SMA、Ruflex、CMHB、SAC等)。斷級配瀝青混合料在礦料組成中增加碎石比例以提高其高溫穩定性與抗滑性，并摻入較多礦粉以減少空隙率，降低其透水性[1]。其中SAC級配瀝青混合料在上世紀80年代由沙慶林院士提出，通過多條試驗路鋪筑逐步形成了整套礦料級配設計和檢驗方法[2]。

劉國耘[3]基于粒子干涉理論，提出了間斷級配瀝青混合料的體積設計方法。劉朝暉、彭波等[4-5]通過室內試驗表明SAC系列瀝青混合料具有良好的路用性能。

瀝青混合料的空隙率與空隙分布特征直接影響其各方面性能，SAC瀝青混合料技術指標要求空隙率為3%～5%，CMHB技術指標要求空隙率小于3%，可見空隙率對間斷級配瀝青混合料的性能影響顯著，因此斷級配瀝青混合料空隙率的計算至關重要。

傳統的瀝青混合料空隙率表干法受試驗環境和操作者影響較大，隨著計算機圖像技術在瀝青混合料研究中廣泛應用，基于CT的圖像技術可以識別瀝青混合料內部細觀結構并對其體積參數進行分析。Masad、Wang等[6-7]最先采用CT技術研究SGC旋轉壓實和Westrack試驗路芯樣瀝青混合料的空隙二維分布特征。裴建中等[8]對OGFC-13的豎向空隙分布特征進行分析，結果表明CT得出的空隙二維特征圖像可用來定性描述其分布特征。Wang等[9]對AC-16、SMA-16、OGFC-16內部空隙二維分布特征進行分析，表明SGC成型試件空隙呈“兩頭大、中間小”，試件徑向呈“外大內小”的特征。

雖然CT技術在瀝青混合料內部空隙研究領域已得到較好的運用，但上述研究在圖像處理中仍局限于斷層圖像的二維分析，未能反映空隙的三維特征及連通狀態，因此需要對瀝青混合料的CT圖像進行三維重構，并利用三維圖像計算瀝青混合料的空隙率。

為此，采用馬歇爾標準擊實與靜壓兩種成型方法制作SAC-16圓柱體試件，應用CT設備對試件進行掃描并獲取其斷層圖像，采用Mimics圖像處理系統進行空隙分離及三維重構，研究空隙的三維分布及連通狀態，并對空隙體積及空隙率進行了計算。

1 CT試驗

1.1 試驗材料

采用遼河90#基質瀝青，技術指標見表1。石灰巖集料和礦粉產自哈爾濱阿城區，SAC-16瀝青混合料級配組成如圖1所示。

表1 瀝青技術指標

圖1 SAC-16礦料級配曲線

采用馬歇爾試驗確定SAC-16瀝青混合料最佳油石比為3.9%，空隙率為4.2%，毛體積密度為2.460g/cm3。按《公路工程瀝青及瀝青混合料試驗規程》(JTG E20-2011)規定的馬歇爾標準擊實和靜壓兩種方法成型試件，馬歇爾擊實試件(SAC-16M)尺寸為Φ101.6×63.5mm，靜壓試件(SAC-16P)尺寸為Φ100×100mm。

1.2 CT掃描方案

研究采用Siemens螺旋CT設備對SAC-16M和SAC-16P試件進行掃描，設備掃描參數見表2。沿試件高度掃描寬度為1mm，SAC-16M和SAC-16P試件橫斷面CT掃描圖像數量分別為60張和120張。沿試件直徑方向環形掃描，每20°掃描一層，厚度設置為1mm。SAC-16M和SAC-16P試件均得到18張縱斷面掃描圖像，掃描位置如圖2所示。

表2 X-CT設備掃描參數

圖2 CT掃描示意圖

CT掃描后獲取的SAC-16M連續斷層圖像如圖3所示。

圖3 SAC-16M試件斷層掃描圖片

2 空隙三維重構

2.1 CT閾值分析

X射線穿過非均勻物質后，其強度將會按指數形式衰減，X射線通過物體后強度Iout與入射強度Iin關系式為：

3.利用課文內容，完善學生心理素質。語文課同其他學科相比，有一個獨特的特點，即思想性和文學性完美的統一。它不像政治課以直接的思想訓導引起學生的厭倦，而是把思想教育不知不覺寓于藝術美感之中。它們都極容易感染青年學生，都是使青年學生心理健康向上的催化劑，她們一旦被學生的心理所接受，不但能使學生形成崇高而偉大的理想，正確的人生觀，而且會使學生的心理素質達到高層次的境界。

(1)

式(1)中，l為X射線在均勻物質中的傳播距離；μ為物質對X射線的衰減系數。

由此可知，某物質的密度可通過X射線的衰減系數來體現，相比于物質密度，研究人員更為關注非均勻物質間的密度差異，因此引入CT值的概念，表述為：

(2)

由式(2)可知某物質CT值與其密度成正比關系。

為了將瀝青混合料中的空隙與集料、瀝青等材質進行分離，首先應準確獲取空隙的CT值范圍。本研究應用Profile line采集空隙徑向“邊緣-中心-邊緣”范圍內的CT值變化。圖4為SAC-16P試件第50斷層圖像中6個空隙的CT值隨徑向長度“邊緣-中心-邊緣”的變化規律。

圖4 空隙的CT值隨徑向長度變化

從圖4可以看出，空隙邊緣處CT值遠大于中心，空隙兩側邊緣至內部存在比較明顯的過渡區，長度為1～2mm，CT值由邊緣處正值在1～2mm范圍內迅速下降至-1000HU左右。每斷層圖像至少采集6個空隙的CT值，統計發現空隙CT值介于-1024～578HU之間，進而選取578HU作為空隙的CT閾值上限進行后續處理。

理論上空氣的衰減系數近似為0，代入式(2)可得空氣CT值應為-1000HU，而本研究中空隙內部CT值在-1000HU左右，邊緣處均高于-1000HU，主要是由于瀝青混合料中空隙、集料、瀝青三相材料間混合狀態的復雜性，空隙與高密度材料的迅速過渡必然使得空隙的CT值也存在一定的過渡區域。

2.2 空隙分離

采用578HU作為空隙的CT閾值上限，對SAC-16M、SAC-16P斷層圖像中空隙進行分離，圖5(a)、圖5(b)分別為SAC-16M試件第38斷層和SAC-16P試件第46斷層空隙分離后的圖像。

圖5 斷層空隙分離圖像

2.3 三維圖像重構

采用像素閾值法對空隙表面輪廓進行軸向三維重構，由于試驗設備和掃描方式的原因，本研究中試件頂部和底部的部分斷層圖像未能清晰地顯示空隙與集料、瀝青的像素區間，分離時會將部分集料誤認為空隙，導致三維重構與體積組成分析出現誤差。為此，三維重構前去除了上述部分斷層圖像，研究對象高度也隨之降低，SAC-16M試件高度降至54mm，SAC-16P試件高度降至94mm，圖6為調整后的空隙三維重構圖像。

3 空隙率計算

采用Mimics圖像處理系統中的體積計算功能對SAC-16M和SAC-16P試件的三維空隙率進行計算，表3為由三維重構圖像計算的空隙體積參數。

表3 三維空隙體積計算

從表3可以看出，由于成型時擊實功的影響，靜壓法成型試件的空隙率大于馬歇爾擊實試件，SAC-16P試件空隙率大于5%，并不滿足3%～5%的空隙率要求，試件成型方法對于瀝青混合料空隙率影響較大。雖然《公路工程瀝青及瀝青混合料試驗規程》(JTG E20-2011)考慮國內無法迅速普及SGC、GTM等先進成型方法的實際情況，保留了較為簡單實用的靜壓成型方法，但通過研究發現，靜壓法成型導致空隙率過大，不適用于SAC-16試件成型。

為了分析三維空隙計算與試驗實測空隙率的差異，采用表干法對試件的空隙率進行測試，比較結果見圖7。

圖7 三維圖像計算與表干法實測空隙率比較

從圖7可以看出，利用三維圖像計算得出的計算空隙率要高于表干法實測空隙率，這主要是由于試件表面開口空隙的影響。

表干法測量空隙率計算公式為：

VV=(1-γf/γt)×100

(3)

式中，γt為試件的理論最大相對密度。

試件的毛體積密度γf通過式(4)計算：

γf=ma/(mf-mw)

(4)

式中，ma為干燥試件空中質量；mf為試件表干質量；mw為試件的水中質量。

SAC級配中提高了碎石比例，較普通密級配要粗，雖然增加了礦粉含量以減少空隙率，但SAC與外界相連通的開口空隙體積和數量較密級配均有所提高。仍采用表干法進行毛體積密度測量時，試件從水中取出后，開口空隙中的水會自然流出，用毛巾擦拭時又會吸走一部分水，導致測量值mf小于實際值，由式(3)和式(4)推導可知，表干法測量空隙率也將小于實際值。而三維重構能夠真實模擬空隙的空間分布狀態，且計算空隙率考慮了瀝青混合料的開口空隙，因此空隙率計算值與真實值更為接近。

4 結論

采用CT設備對標準馬歇爾和靜壓兩種方法成型的SAC-16試件進行掃描，利用圖像處理技術對瀝青混合料的空隙率進行計算，得出如下結論：

(1)Mimics圖像處理系統可以實現空隙的三維重構，試件側向邊緣均存在較多連通而成的大空隙，試件內部空隙多為不連通且分布均勻的小空隙。

(2)靜壓法成型的SAC-16試件空隙率大于5%，并不滿足3%～5%的空隙率要求，靜壓法成型將會導致瀝青混合料空隙率過大，不適用于SAC-16試件成型。

(3)三維圖像計算SAC-16M與SAC-16P試件空隙率分別為4.5%、7.3%，均高于表干法實測值，產生差異主要是由于開口空隙的影響，利用三維圖像得到的計算空隙率與真實值更為接近。