宋伯石

摘 要： 變形測量是砂箱模型試驗中的一個重要環(huán)節(jié)。通過對數(shù)字圖像測量中的標點法及相關(guān)法的介紹，給出了圖像獲取、畸變矯正及位移、應(yīng)變等參數(shù)計算的方法，并通過PIVLAB及自行設(shè)計的標點法程序?qū)Ρ粍訕赌Ｐ驮囼灲Y(jié)果進行了觀測，結(jié)果表明兩種方法均可完成試驗中土體位移、應(yīng)變發(fā)展的定性、定量測量，且兩者各具優(yōu)勢，其中標點法更適于大范圍觀測，相關(guān)法更適于局部變形測量，應(yīng)根據(jù)具體試驗規(guī)模選擇合適的方法。

關(guān)鍵詞： 數(shù)字圖像測量； 模型試驗； 標點法； 圖像相關(guān)法； PIV； 圖像畸變

1. 引言

砂箱模型試驗因為其邊界條件可控、相比原位實驗可以在更短時間內(nèi)獲得土體應(yīng)力、應(yīng)變、位移等數(shù)據(jù)，已經(jīng)是巖土工程領(lǐng)域中的重要研究手段。在模型試驗中，土體的位移破壞往往是局部變形累積并進一步發(fā)展的結(jié)果，傳統(tǒng)的位移測量手段如拉線法、位移計、百分尺等，因為大多屬于接觸式測量，為減小對模型的干擾，測量往往大而概括，難以細致的獲得土體局部變形的演化信息，因此不易深刻反映土體變形破壞的機理。

隨著計算機可視化和攝像技術(shù)的發(fā)展，數(shù)字圖像測量技術(shù)已經(jīng)被廣泛地運用于土體的位移測量上，例如李元海等[1-4]曾提出應(yīng)用標點法及相關(guān)性法對模型位移進行測量，并開發(fā)了相關(guān)圖像處理軟件；周擁軍等[5]基于圖像測量法對高速公路滑坡與路塹邊坡物理模型試驗的平面位移場進行了測量；毛靈濤等[6]通過數(shù)字標識點圖像相關(guān)法對模型箱內(nèi)土體位移進行了測量，王學(xué)濱等[7]通過數(shù)字圖像相關(guān)法對等應(yīng)變率下不同含水率砂樣的剪切帶進行了觀測。本文在一系列砂箱模型試驗的基礎(chǔ)上，對近年來發(fā)展的數(shù)字圖像測量法的原理、流程進行了詳細介紹，并對土體運動速度、位移、應(yīng)變等參數(shù)的計算進行了闡述。

2. 數(shù)字圖像位移測量法

2.1 圖像的獲取及畸變矯正

影響圖像獲取質(zhì)量的因素主要有相機參數(shù)、被拍攝對象表面處理、光照條件及硬件擺放布局。

在相機參數(shù)方面，高信噪比的傳感器和對高動態(tài)范圍（HDR）的支持是更加理想的，所以盡量選擇使用CCD而非CMOS的相機。在拍攝時應(yīng)在保證光照和景深的前提下盡量使用大光圈，可以有效減少鏡頭畸變[8]。在選取相機像素數(shù)時，應(yīng)將其保證在足夠且合理的范圍內(nèi)，在拍攝砂箱時，用兩個像素表示一個土顆粒是比較合適的[9]。

在被拍攝物體的表面處理方面，標識點法以標識點為基準測量位移，因此標識點應(yīng)與背景顏色對比鮮明。圖像相關(guān)性法依據(jù)模型表面斑紋特征選定位移測量參考點，因此模型表面的斑紋應(yīng)與模型背景有足夠的對比度，斑紋和孔隙應(yīng)該大于圖像中的噪聲。

在光照方面，應(yīng)盡量為被拍攝物體提供足夠且均勻的光照，光線不均或光抖動都會對圖像的后期處理造成干擾。在布局硬件時，相機鏡頭應(yīng)與被拍攝面平行來減少照片的變形。

在拍攝過程中，由于鏡頭與被攝面共面性不足、徑向和切向的鏡頭畸變、被拍攝面存在反光、圖像像素比不為1等因素[10]，會對圖像測量帶來誤差。為了矯正圖像畸變，需對相機進行標定以獲得相機內(nèi)參數(shù)。標定相機的方法有：傳統(tǒng)標定法、主動視覺標定法和自標定法。本文使用的是傳統(tǒng)標定法中的張正友標定法，目前該方法已經(jīng)集成到Matlab的相機標定工具箱中，可以完成相機內(nèi)參數(shù)的計算、圖像的畸變矯正及圖像坐標到世界坐標的轉(zhuǎn)換。

2.2 標點法

標點法即在被拍攝模型表面放置標識點作為物理測量點，通過閾值分割技術(shù)將標識點從圖像中識別出來，然后計算不同試驗階段各個標識點的位置變化，以標識點的位移代替其周圍土體的位移，最后通過插值方法獲得全場的位移場。

標識點的提取屬于圖像分割問題，本文使用的方法是基于顏色和飽和度的分割技術(shù)。在砂箱模型試驗中，標識點的顏色應(yīng)與土體有較高的對比，因此在HSV圖像模式下圖像中標識點的色調(diào)值（H）和飽和度值（S）會和土體的有很大不同，通過選擇合適的閾值，將圖像二值化便可以精確的提取到標識點并計算出質(zhì)心位置。

在得到質(zhì)心位置后，便要對土體變形進行解釋，標識點的位移可以通過測量圖像上標識點的坐標變化直接求得，應(yīng)變的計算則可借助有限元法中的四邊形等參數(shù)變換來計算[2]。

2.3 圖像相關(guān)法

圖像相關(guān)法又稱無標點法，是指利用巖土材料本身的紋理特征，或在其表面進行制斑，對變形前后巖土體進行拍攝，然后對一對圖像進行相關(guān)分析來進行變形的測量與解釋。本文使用的是PIV法即粒子圖像測速法。該方法的使用一般基于四個前提[11]，（1）紋理特征均布在被測表面。（2）被觀測的斑紋能夠很好地代表巖土體的位移。（3）目標圖片和參考圖片要足夠的相似。（4）在位移測量中，用于相關(guān)性計算的形函數(shù)應(yīng)該始終保持一致。

在進行PIV法測量時，一般分為兩個步驟，一是對參考圖像進行網(wǎng)格劃分，對每一網(wǎng)格子區(qū)域在目標圖像中進行相關(guān)性計算，然后根據(jù)相關(guān)性峰值確定該子區(qū)域在變形后圖像中的整數(shù)像素坐標。在參考圖像中劃分子區(qū)域網(wǎng)格，并在目標圖像中的搜索區(qū)域內(nèi)計算相關(guān)系數(shù)，圖中相關(guān)系數(shù)峰值點對應(yīng)的整數(shù)像素坐標即為子區(qū)域中心發(fā)生位移后的位置。二是通過亞像素插值函數(shù)對第一步得到的整像素位移進行優(yōu)化。通常使用雙三次樣條插值法或者高斯插值法，優(yōu)化后的位移值精度可以達到0.01倍像素尺寸。通過上面兩步得到位移場后，便可以進行應(yīng)變的計算，以四個相鄰像素點的位移作為四邊形單元的四個節(jié)點，則計算方法與標點法中應(yīng)變的計算方法相同。

3. 應(yīng)用實例驗證

為驗證上述方法在砂箱模型試驗中的可行性，筆者借助水平荷載作用下被動樁模型試驗為應(yīng)用實例來做簡單說明，其中標點法由筆者自行開發(fā)的Matlab程序?qū)崿F(xiàn)，相關(guān)法則借由PIVLAB實現(xiàn)[12]。

3.1 試驗概述

模型箱尺寸為長1.0m，寬0.4m，高0.3m，底部由鋼架支撐。模型箱左右兩面為鋼化玻璃，上部為有機玻璃蓋，以方便拍攝試驗圖像。水平荷載通過在推土板后放置配重來施加，10kg為一個配重單位。試驗用土為中粗砂，干密度為2.35g/cm3，孔隙度為0.36g/cm3。試驗時分層填土，多遍夯實，直至砂土填至設(shè)計高度，并放置標識點和石英砂來制斑，然后由數(shù)碼相機從上向下拍攝模型表面，如圖1所示。

3.2 試驗結(jié)果

在圖1中，選擇樁后布置紅色標識點的矩形范圍作為測量范圍，使用上述方法對實驗結(jié)果進行了測量。

圖2為在水平荷載下樁后土體的速度及位移場，其中a組圖為標點法結(jié)果，b組圖為相關(guān)法結(jié)果。圖中S為推土板位移，Q為所施加配重，即反映水平荷載大小。從圖6可以看出，在水平荷載作用下，土體向樁間及樁身處擠壓，形成相對位移并在樁間處形成明顯的位移拱。隨著配重增加，推土板繼續(xù)前進，土體最終會從樁間完全擠出。

為速度矢量圖，從上圖可以看出，標點法由于受到標識點數(shù)量限制，采樣測量點要少于相關(guān)法，但在標識點足夠的情況下，采用適當?shù)牟逯邓惴▽?shù)據(jù)進行處理后，也可以得到比較理想的試驗結(jié)果，并且由于測點較少，其測量速度要優(yōu)于相關(guān)法。試驗中由于條件限制使用的是單側(cè)光源，模型被拍攝表面布光有些不均，在使用相關(guān)法時反光處測量會精度下降，需要插值算法進行優(yōu)化結(jié)果。另外，樁后樁尖處是位移突變區(qū)，當配重加至110kg時此處土體產(chǎn)生拉裂縫，繼續(xù)使用相關(guān)法將產(chǎn)生測量誤差。而在使用標點法時，標識點則可不受光照條件被全部提取，并且土體大變形不會對標識點的測量產(chǎn)生影響。

以上試驗結(jié)果表明，數(shù)字圖像測量法可以獲得在不同荷載條件下土體的位移場及速率場，測量結(jié)果能保證足夠的精度且規(guī)律性良好，因此該方法適合在砂箱模型試驗中進行土體變形的定性及定量測量。

4. 結(jié)論

針對數(shù)字圖像測量法中的標點法及相關(guān)法，從在砂箱模型試驗中應(yīng)用的角度出發(fā)，對圖像的獲取、畸變矯正及兩方法的基本原理進行了闡述，并通過被動樁模型試驗進行了檢驗，結(jié)果表明兩種方法均能全面、細致的測量土體變形。

試驗中兩種方法各自體現(xiàn)出了不同的優(yōu)缺點，標點法在使用時不受觀測區(qū)域變形限制，圖像處理更快，但在布置標識點時有可能擾動土體，并且測點有限；相關(guān)法無需布置標識點，對土體無擾動，并且測點的數(shù)量及測量范圍設(shè)置更加靈活，但測量精度易受光線變化及測量區(qū)變形突變影響，且圖像處理時間較長。綜合來看，標點法更適宜大范圍變形測量，相關(guān)法更適宜局部化測量，因此在實際使用中，應(yīng)根據(jù)模型規(guī)模選用合適的方法。

目前數(shù)字圖像測量技術(shù)在使用中主要以模型表面的數(shù)碼照片為試驗數(shù)據(jù)，對于反映土體內(nèi)部變形情況方面還存在著局限性，隨著激光掃描技術(shù)及透明砂等試驗手段的不斷發(fā)展，研究人員現(xiàn)已逐漸可以獲得土體內(nèi)部變形的圖像數(shù)據(jù)，真正達到全場實時測量，因此數(shù)字圖像測量法在巖土模型試驗研究中必將擁有更加廣泛的應(yīng)用前景。

參考文獻：

[1] 李元海and朱合華，基于圖像分析的實驗?zāi)Ｐ妥冃螆隽繙y標點法. 同濟大學(xué)學(xué)報：自然科學(xué)版，2003.31（10）： p. 1141-1145.

[2] 李元海，靖洪文，and曾慶有，巖土工程數(shù)字照相量測軟件系統(tǒng)研發(fā)與應(yīng)用.巖石力學(xué)與工程學(xué)報，2006. 2.

[3] 李元海，et al.，基于數(shù)字照相的砂土剪切變形模式的試驗研究.同濟大學(xué)學(xué)報：自然科學(xué)版， 2007. 35（5）： p. 685-689.

[4] 李元海and 靖洪文， 基于數(shù)字散斑相關(guān)法的變形量測軟件研制及應(yīng)用. 中國礦業(yè)大學(xué)學(xué)報， 2008. 37（5）： p. 635-640.

[5] 周擁軍，任偉中， <一種僅需距離控制的模型試驗平面位移場單像視覺測量_周擁軍.pdf>. 巖石力學(xué)與工程學(xué)報， 2009. 28（4）： p. 6.

[6] 毛靈濤， et al.， 數(shù)字標記點圖像相關(guān)法在模型實驗中的應(yīng)用. 遼寧工程技術(shù)大學(xué)學(xué)報： 自然科學(xué)版， 2013（10）： p. 1367-1373.

[7] 王學(xué)濱， et al.， 基于數(shù)字圖像相關(guān)方法的等應(yīng)變率下不同含水率砂樣剪切帶觀測. 巖土力學(xué)， 2015. 36（3）.

[8] D， J. and A. S， Improved Image-Based Deformation Measurement in the Centrifuge Environment. Vol. 36. 2013. 1-14.

[9] Westerweel， J.， Effect of Sensor Geometry on the Performance of PIV Interrogation. 2000： Springer Berlin Heidelberg. 37-55.

[10] Nübel， K.， Experimental and numerical investigation of shear localization in granular material. 2002.

[11] Stanier， S.， et al.， Vermiculate artefacts in image analysis of granular materials. Computers & Geotechnics， 2016. 72： p. 100-113.

[12] 姚秋玲，et al.，堤基管涌微觀機理模型試驗研究.中國水利水電科學(xué)研究院學(xué)報，2014. 12（1）： p. 1-7.