張超 李青 程一峰 劉章杰

摘要：巖土內(nèi)部的滲透（孔隙）水壓力采用滲壓計測量，進(jìn)而進(jìn)行滲流監(jiān)測，其測量的是標(biāo)量，存在動態(tài)響應(yīng)差、量程小、非直接測量等缺陷，難以滿足巖土中滲流監(jiān)測發(fā)展的需求。鑒于此，該文設(shè)計一種基于熱擴散原理的遠(yuǎn)程在線滲流計，主要由電源模塊、溫度檢測模塊、恒功率加熱模塊、顯示模塊組成;由三線制鉑電阻恒流源測溫電路檢測進(jìn)出水兩端溫度差，檢測結(jié)果返回顯示到STM32并經(jīng)過RS-485總線傳至主機進(jìn)行數(shù)據(jù)分析。實驗結(jié)果表明，該滲流計可實現(xiàn)在線監(jiān)測巖土中的滲流情況，能有效檢測巖土中滲流的速度，分辨率達(dá)到0.05mL/s，具有良好的靈敏度、準(zhǔn)確性和重復(fù)性，適用于大部分巖.土中的滲流測量。

關(guān)鍵詞：巖土;滲流;微小測量;熱傳遞

中圖分類號：TP23 文獻(xiàn)標(biāo)志碼：A 文章編號：1674-5124（2019）07-0092-05

收稿日期：2018-06-13;收到修改稿日期：2018-07-25

基金項目：國家重點研發(fā)計劃項目（2017YF00804604）;浙江省重點研發(fā)計劃項目（2018003040）

作者簡介：張超（1994-），男，內(nèi)蒙古包頭市人，碩士研究生，專業(yè)方向為檢測技術(shù)。

通信作者：李青（1955-），男，浙江杭州市人，教授，主要從事檢測技術(shù)方向研究。

0 引言

滲流形成的基坑塌陷[1]、滑坡、垮壩[2]等危害每年都會造成巨大的人員傷亡和經(jīng)濟損失。由于巖土自身特殊的物理性質(zhì)，巖土空隙中運動的地下水也就是滲流的測量[3]，很難通過傳統(tǒng)的聲、光、電等手段進(jìn)行檢測，巖土中的滲流具有流速低，流量小，組成復(fù)雜等特性[4]。目前大量使用的巖土滲流測量設(shè)備是孔隙水壓力計和測壓管[5]，其測量原理為測量巖土中高度不同的兩個點之間的水力坡降，通過達(dá)西定律Q=K·A·H₁-H₂/L來得出滲流量Q的大小進(jìn)而得到滲流速度的大小[6]。其中K是滲透系數(shù)，對于同一種土壤而言，滲透系數(shù)是不變的，其通過土壤滲透儀來測得，A為孔隙水壓力計的截面積。其最終測得的是標(biāo)量，不能反映出滲流的方向，測壓管只能測量滲流壓力變幅小的部位，且其滯后時間長，可靠性低[7]。滲壓計具有量程小，分辨率低，施工難度大等缺點[8]。目前尚無直接測量滲流速度的傳感器。

針對目前測量方法的缺點和不足，本文設(shè)計一種基于熱擴散的原理[9]在巖土中測量微小流量的液體流量計，采用微型陶瓷加熱棒，微型鉑電阻放置在細(xì)管中組成傳感器。利用熱擴散原理對巖土中的滲流實現(xiàn)較高精度的測量。設(shè)計采用RS-485總線通信的方式將采集到的電信號實時上傳到實驗主機，實驗主機對鉑電阻采集到的電信號和時間建立數(shù)學(xué)模型找出電信號變化和滲流速度之間的對應(yīng)關(guān)系。

本裝置是一個實時在線監(jiān)測設(shè)計裝置，安裝在野外采用太陽能電池板進(jìn)行供電，在日后還可以利用以太網(wǎng)實現(xiàn)遠(yuǎn)程數(shù)據(jù)共享，可以實現(xiàn)巖土中滲流的遠(yuǎn)程監(jiān)測，具有深遠(yuǎn)意義。

1 巖土中滲流的監(jiān)測方法與方案設(shè)計

該設(shè)計具體方法是采用兩個直徑為3～的PT100鉑電阻和一個直徑為3.3mm的陶瓷加熱棒，將鉑電阻和加熱棒通過膠水固定在直徑為6mm長度為200mm的聚碳酸酯細(xì)管內(nèi)，加熱棒固定在細(xì)管中央長度的位置，鉑電阻固定在細(xì)管兩端的位置。將細(xì)管和直徑為55mm的漏斗用膠水粘緊，將直徑50mm厚度5mm的透水石放置于漏斗內(nèi)用膠水粘緊，將此裝置放置于直徑和漏斗一致的60mm粗管內(nèi)，用膠水固定。滲流傳感器的設(shè)計方案如圖1所示。當(dāng)細(xì)管中有水流過時，水先經(jīng)過一端的鉑電阻，此時鉑電阻將溫度值轉(zhuǎn)化為電信號，然后水流過加熱棒進(jìn)行加熱，最后水流過另一端的鉑電阻并測量其溫度并轉(zhuǎn)化為電信號。加熱棒采用恒功率加熱[10]，兩端鉑電阻采集到的溫度差會隨著流速的變化而變化。

2 細(xì)管內(nèi)流體傳熱模型及仿真

2.1 層流原理及數(shù)學(xué)推導(dǎo)

水是一種粘性流體，當(dāng)水通過固定長度和內(nèi)徑的細(xì)管時，流體存在一定的阻力，細(xì)管越細(xì)，阻力越大。聚碳酸酯材質(zhì)的細(xì)管光滑的內(nèi)壁將導(dǎo)致穩(wěn)定的層流狀態(tài)。在這種條件下，最大的流速是在中心區(qū)，細(xì)管內(nèi)壁附近的流體保持幾乎靜止不動，根據(jù)層流原理對細(xì)管內(nèi)水的流動情況進(jìn)行數(shù)學(xué)分析[11]。

假設(shè)粘度為μ的流體在半徑為R的水平細(xì)管內(nèi)作層流運動，現(xiàn)取管軸中心處一半徑為r，長度為l的流體柱作為分析對象，作用于流體柱兩端面的壓差為△p，則作用在流體柱上的推動力為f=△pπr²。

設(shè)距離細(xì)管中心處r的流體速度為V_r（r+dr）處的相鄰流體層的速度為（V_r+dV_r），則流體速度沿半徑方向的變化率（即速度梯度）為dV_r/d_r，兩個相鄰流體層所產(chǎn)生的粘滯力為τ_r。層流時粘滯力遵循牛頓粘性定律，即τ_r=-μdV_r/d_r。

作用在流體柱上的阻力為τ_r=-μdV_r/d_r。流體作等速運動時，推動力與阻力大小相等，方向相反，所以△pπR²=-2πRlμdV_r/d_r。

當(dāng)r=r時，μ_r=μ_rr：當(dāng)r=R時，μ_r=0。故，將△△pπR²=-2πRlμdV_r/d_r積分后可得到

式（1）是流體在細(xì)管內(nèi)作層流運動時的速度分布表達(dá)式，表示在某一壓力降下，速度與半徑的關(guān)系是拋物線型，即在管路中心處的速度最大，面到管壁處的速度為零。流經(jīng)厚度為dr的流體柱的體積流量為dF=U_rdA ，其中dA=2πrdr。那么細(xì)管內(nèi)流體的流量就是這些流體柱的流量的積分F=

式（2）即為液體在細(xì)管內(nèi)層流狀態(tài)下的數(shù)學(xué)公式?？梢钥闯?，當(dāng)細(xì)管的長度一定，液體一定（粘度確定）時，流量與半徑的四次方和壓差成正比[12]。

2.2 細(xì)管中水的速度場和溫度場仿真

采用Ansys Fluent軟件來進(jìn)行建模并對細(xì)管中流體的流動狀態(tài)及溫度進(jìn)行仿真來驗證2.1中的數(shù)學(xué)推導(dǎo)。如圖2所示為細(xì)管中流體的溫度場，圖3為細(xì)管中流體的速度場，其中加熱棒為恒功率4W，水由右向左以0.05mL/s的速度流動。

從溫度場仿真圖中可得出，在加熱棒4W恒功率，水流速度0.05mL/s（一滴水）的情況下，上下游兩端溫度差大約為6.8℃。由速度場仿真圖可看出，速度與半徑的關(guān)系是呈拋物線型，即在細(xì)管中心處的速度最大，越靠近管壁流動的速度越慢，在緊貼管壁處速度幾乎為零。因此在滲流計制作過程中，要將加熱棒和PT100盡可能固定在細(xì)管中軸線上方可得到最好的測量效果。

3 硬件電路設(shè)計

3.1 鉑電阻測溫電路

測溫電路主要是對PT100鉑電阻上的電壓信號進(jìn)行采集和放大。通過STM32內(nèi)部A/D采樣端口對測溫電路輸出的電壓信號進(jìn)行采樣并顯示在屏幕上，運用RS-485總線將數(shù)據(jù)上傳到主機進(jìn)行處理。圖4為測溫電路的設(shè)計框圖。

測溫電路主要包括PT100鉑電阻、信號放大電路、濾波電路、STM32控制器以及電源電路等[73]，最終實現(xiàn)對鉑電阻電壓信號的實時采集和上傳。

1mA恒流源產(chǎn)生電路使用Howland運放電流源[14]，如圖5所示電路，有兩個電阻反饋網(wǎng)絡(luò)。在保持輸入電壓V_in不變的情況下，假設(shè)因負(fù)載電阻R_L減小而引起輸出電流i_L增大，則節(jié)點c和d間的電壓升高，則流過R₂和R₄的電流i_D和i_E增大，因R₂不變，則節(jié)點a的電壓升高，根據(jù)運算放大器“虛短”的概念，節(jié)點b的電壓也要升高，在相同輸入電壓的情況下，此時流過電阻R₁的電流減小，再根據(jù)運算放大器“虛斷”的概念，則流過R₃的電流也減小，而輸出電流為流過R₃和R₅的電流之和，所以此時輸出電流減小，通過閉環(huán)反饋從而抑制了輸出電流的增加，以達(dá)到恒流的作用。其恒流性能良好，最終可得到輸出電流與輸入電壓成正比，可以得出：

鉑電阻選用德國Heraeus高精度A級PT100薄膜型鉑電阻，0℃時阻值誤差為±0.06%，測溫范圍為-50～300℃，響應(yīng)時間t_0.5=0.05s，t_0.9=0.15s。采用不銹鋼套管進(jìn)行防水，總尺寸為φ3mm×15mm，多股鍍銀屏蔽線為引出線。

采用三線制PT100鉑電阻來減小導(dǎo)線電阻所帶來的附加誤差，使用op07c電壓運算放大器采集鉑電阻上的電壓信號并進(jìn)行100倍放大。如圖6所示為三線制鉑電阻測量原理，當(dāng)電橋平衡時R₁（R_T+r₆）=R₂（R₃+r₄），因為R₁=R₂，所以R_T+r₆=R₃+r₄，r₄=r₅=r₆為三線制鉑電阻導(dǎo)線上的電阻，所以R_T=R₃[15]。

3.2 測溫電路標(biāo)定

分別采用4個109.73，120，125，140Ω的RX70型萬分之一精度的高精密線繞電阻來代替電路中的鉑電阻對測溫電路進(jìn)行標(biāo)定，標(biāo)定結(jié)果放大倍數(shù)約為99.93，平均誤差為0.1655mV，對應(yīng)到PT100鉑電阻分度表中可得到測溫精度約為0.1℃，滿足對巖土中滲流監(jiān)測的需求。

4 實驗裝置的設(shè)計

4.1 流速與溫度的變化關(guān)系實驗

為驗證流速與溫度變化的關(guān)系，設(shè)計如圖7所示的實驗裝置，分為限流裝置、溫度信號檢測和放大電路、加熱棒電路，將水袋用支架掛起，水袋、滲流計細(xì)管、限流閥之間用軟管連接，滲流計細(xì)管用鐵甲臺固定并保持水平，調(diào)整限流閥旋鈕獲得不同的流速，加熱棒恒功率加熱。通過PT100采集到的滲流計細(xì)管兩端的溫度差來計算出流速和流量的大小。

4.2 滲流計標(biāo)定實驗裝置

為了得到巖土中滲流的真實情況，更好地模擬巖土中滲流的條件和環(huán)境，設(shè)計了如圖8所示的標(biāo)定實驗裝置，長度110cm直徑6cm的透明管和DN50口徑的彎頭以及水平的長度為40cm直徑為6cm的透明管進(jìn)行連接，并采用弧形支架固定在底板上保持穩(wěn)定。水平放置的透明管內(nèi)是滲流計細(xì)管和漏斗以及透水石的結(jié)合體，使用防水膠進(jìn)行連接和填充縫隙。在垂直透明管中加入一定量的土，然后加入水，水由于重力滲過土然后流經(jīng)彎頭接著流過透水石集聚在漏斗內(nèi)，當(dāng)漏斗中的水充滿一半時水位和滲流計細(xì)管平齊，此時水開始充滿細(xì)管，分別流過細(xì)管中的PT100，加熱棒、PT100，整套裝置連接處均采用防水膠進(jìn)行密封處理。目前采用控制限流閥來控制水的流速來模擬長管內(nèi)裝土的狀態(tài)。

5 實驗數(shù)據(jù)處理

為了模擬出土壤中的滲流速度，采用限流閥來對流速進(jìn)行限定，查閱各類土的滲透系數(shù)如表1所示。

由于巖土中的滲流是為了監(jiān)測垮壩、滑坡、基坑塌陷等一系列災(zāi)害事故，其巖土大多為顆粒直徑較大的粗砂或松散堆積物。故選用粉土質(zhì)砂（顆粒直徑約為0.05mm）滲透系數(shù)為0.6m/d，換算到長管截面為0.18mL/s，以每滴水0.05mL來計算，流速大約為3滴/s～4滴/s。使用限流閥將流速限制為1滴/s，2滴/s，3滴/s，4滴/s，在非同一天的同一時間、長管中水位一致、加熱棒功率恒定4W的情況下分別進(jìn)行多次實驗，實驗結(jié)果如圖9所示。其中橫坐標(biāo)代表時間，縱坐標(biāo)為下游PT100測溫的電壓值大小，電壓越大代表溫度越高。

可看出不同流速下，達(dá)到穩(wěn)定的時間，斜率都有明顯不同，其中0.05mL/s（1滴/s）大約需4500s達(dá)到穩(wěn)定，0.10mL/s（2滴/s）大約需3300s達(dá)到穩(wěn)定，0.15mL/s（3滴/s）大約需2100s達(dá)到穩(wěn)定，0.20mL/s（4滴/s）大約需390s達(dá)到穩(wěn)定。

根據(jù)圖9中的數(shù)據(jù)，將PT100上的電壓值轉(zhuǎn)化成溫度值。保持加熱棒恒功率4W加熱時，在流量為0.05mL/s（1滴/s）時，下游PT100上的最高溫度為51.1℃，此時上游PT100上的溫度值，也就是人水口的初始溫度為22.1℃，經(jīng)過4600s以后溫度差約為29℃，并保持不變。在流量為0.1mL/s（2滴/s）時，下游PT100的最高溫度為38.6℃，入水口溫度為22.1℃，經(jīng)過3000s以后溫度差約為16.5℃并保持不變。在流量為0.15mL/s（3滴/s），上下游溫度差約為2℃并保持不變。0.2mL/s及以上測出的溫差太小，難以檢測。故采用4W加熱棒恒功率加熱，細(xì)管內(nèi)PT100相距加熱棒都為10cm的情況下，適用于0.05～0.15mL/s的滲流速度，此時測量精度最高。

6 結(jié)束語

本文設(shè)計了一種基于熱擴散原理的液體微流量計用來測量巖土中的滲流。實驗結(jié)果表明這種檢測方式和傳統(tǒng)滲壓計相比，具有直接測量、測量范圍廣、重復(fù)性好等優(yōu)點，可以測量到最低0.05mL/s的滲流速度，將來可以實現(xiàn)多點測量，從而實現(xiàn)對巖土中滲流的在線監(jiān)測，是一種新的測量巖土中滲流的解決方案。

參考文獻(xiàn)

[1]朱海波.滲流對地下室基坑沉降的影響及預(yù)防[J].工程建設(shè)與設(shè)計，2017（16）：28-29.

[2]楊淑艷，趙玉坤.淺談水庫滲流安全問題[J].內(nèi)蒙古水利，2014（1）：152-152.

[3]楊裕云，楊紅剛，吳有才.與地下水作用有關(guān)的地質(zhì)災(zāi)害[J].水文地質(zhì)工程地質(zhì)，2003（S1）：1-7.

[4]劉泉聲，崔先澤，張程遠(yuǎn)，等.多孔介質(zhì)中沉積顆粒脫離特性試驗研究[J].巖土工程學(xué)報，2015，37（4）：747-754.

[5]孫嘉彬.張峰水庫大壩滲流監(jiān)測和儀器選型[J].山西水利，2008，24（5）：77-79.

[6]程冠初，巖土介質(zhì)滲流以及輸運從孔隙尺度到達(dá)西尺度的研究[D].杭州：浙江大學(xué)， 2012.

[7]劉玉峰，王消川.對土石壩滲流安全監(jiān)測儀器的幾點認(rèn)識[J].大壩觀測與土工測試，2000（5）：11-13.

[8]李剛，宋先海.滲壓計在滲流監(jiān)測中的誤差分析及對策[J].人民長江，2010，41（15）：59-62.

[9]楊世銘.傳熱學(xué)[M].4版.北京：高等教育出版社.2006：20-32.

[10]劉偉帥.基于熱擴散式恒功率流量測量研究[D].大慶：東北石油大學(xué)，2015.

[11]李雯，熱式質(zhì)量流量計的設(shè)計熱式質(zhì)量流量計的設(shè)計[D].杭州：浙江大學(xué)，2007.

[12]魏列江，液體微小流量的非定常流測量原理與方法的研究

[D].蘭州：蘭州理工大學(xué)，2009.

[13]MAHNAM A，YAZDANIAN H，SAMANI M M.Comprehensive study of Howland circuit with non-idealcomponents to design high performance current pumps[J].Measurement，2016（82）：94-104.

[14]NIE Y Z.Design of high precision temperature sensor basedon platinum resistance[J].Applied Mechanics&Materials，2014，539：177-180.

[15]MORRIS A S，LANGARI R.Measurement andInstrumentation（Second Edition）[M].Pittsburgh：AcademicPress，2016.

（編輯：徐柳）