流體速度的測量方法
2023-09-22 02:59:20 2
專利名稱：流體速度的測量方法
技術領域：
本發明涉及一種利用超聲波測量流體速度的方法，特別是涉及一種精確測量在管中流動的流體的速度地方法。
背景技術：
很多工業領域使用象氣體或液體的流體。例如，通過管道提供城市用氣和原油，建設露天溝渠提供農業用水，在鋼鐵、化工和石油行業等使用循環冷水，建立供水和排水系統設施以向住戶和工廠供水和排洩住戶和工廠產生的汙水。在這些領域，精確地調整和管理流體的流量尤為重要。
一般地，使用流量計來測量流動路徑中流體的流量。近年來，已廣泛使用安裝在流動路徑上的超聲傳感器，通過其接收和發送超聲波來測量流體速度的方法，然後應用測得的流體速度來計算流體流量。相應地，利用超聲傳感器精確測量流體速度顯得尤為重要。
圖1為描述利用超聲傳感器測量流體速度的方法的示意圖。請參閱圖1所示，流體的流量由下面的公式確定Q＝A×VQ流體流量A流體流動路徑的橫截面積V流體的平均速度已知流動路徑中流體的橫截面積和流體速度，可容易地計算流體的流量。假定流體充滿流動路徑，流體的橫截面積就等於流動路徑的橫截面積。
否則，若採用超聲波多束流量計(ultrasonic_multi-beam_flow_meter)來測量流體速度，一般則基於傳輸時間差。更加具體的講，在流動路徑的A點和B點相互對稱地安裝一對超聲傳感器，以使線段AB與流體的流動方向構成一夾角(θ)。假定在流動路徑中流體不流動時超聲傳感器發射的超聲波在流體中的速度為C，流體的平均速度為V，兩超聲傳感器之間的距離為L，從A發射的超聲波到達B所需時間為tAB，從B發射的超聲波到達A的所需時間為tBA，則tAB、tBA給定如下
tBA=LC-Vcos,tAB=LC+Vcos]]>超聲波相對於流體流動的方向向前發射(從A到B)時的傳輸時間tAB比超聲波相對於流體流動的方向向後發射(從B到A)時的傳輸時間tBA較短，其時間差Δt為t=tBA-tAB=2VLcosC2(1-V2cosC2)]]>由於 很小，可忽略不計，流體的平均速度為V=tC22Lcos]]>該速度為流體在對稱安裝的一對超聲傳感器之間直線流動的速度。由流體的橫截面積乘以測得的流體速度所得到的流體流量，產生很大的誤差，這是因為流體的速度是變化的。
圖2A描述一流體管道10，其半徑為r，構成流體流動路徑，其上裝有5對超聲傳感器21，22，23，24和25。圖2B為描述利用圖1中流體管道10和5對超聲傳感器所測量的流體速度分布曲線。參見圖2A，流經流動路徑P管壁的流體速度V很低，而流經流動路徑P中心位置的流體速度高。因此，圖2B所示的流體速度分布曲線為一拋物線。如果流體流動路徑P中的流體是均勻流動的，則形成的速度分布曲線對稱於管道10的中心。沿著管道10的管壁流動的流體由於與管壁發生摩擦，其流速V較低，而流經流動路徑P的中心位置的流體由於與管壁發生的摩擦較少，其流速V較高。因此，將流體的橫截面積乘於由任意一對超聲傳感器測得的流體速度，是不可能計算出流體的流量的。如圖2A所示，該組超聲傳感器沿著流動路徑半徑周界安裝，利用每對超聲傳感器測得的流體速度計算出流體的平均速度，再將流體的橫截面積乘於該平均速度計算出流體的流量。
按照慣例，在利用由每對超聲傳感器獲得的流體速度計算流體的平均速度時，使用加權係數法。該加權係數法將加權係數W與由超聲傳感器測得的流體速度的乘積加在一起。如圖2A所示，5對超聲傳感器21，22，23，24和25，間隔地安裝在流動路徑P的半徑為r的周界上，所測得流體速度V1到V5。將加權係數W1到W5分別與速度V1到V5的乘積求和，計算出流體平均速度Vmean如下Vmean＝W1V1+W2V2+W3V3+W4V4+W5V5於是，將該平均速度Vmean乘於流體的(流體流動路徑)橫截面積A得到流體的流量如下Q＝Vmean×A＝(W1V1+W2V2+W3V3+W4V4+W5V5)×A加權係數W1到W5，由影響流體速度的係數例如摩擦係數確定，因製造商的不同而不同。
如上所示，傳統的加權係數法將加權係數W乘以測得的速度V，以糾正由於流體與管道之間的摩擦、雷諾數(Reynolds_Number)引起的流體湍流(turbulent_flow)等因素所產生的速度誤差。換言之，傳統加權係數法是將加權係數與多對超聲傳感器測得的速度值相乘，計算出平均速度的代表值的近似值。相應地，這種方法不能精確地計算流動路徑中流動的流體的速度的分布曲線，因此，所計算出的流體的流量具有很大的誤差。
發明內容
本發明提供一測量流體速度的方法，該方法考慮流體的流動特徵，將流動路徑的橫截面分成三部分，並提出該三部分合適的速度函數，以反應流動路徑中流動的流體的實際速度分布曲線，從而提供流體速度的精確值。
根據本發明的特徵，本發明所提供的測量流體的流動速度的方法包括設定半徑為r的流體流動路徑的一半徑方向為X軸，在X軸(0≤x≤2r)上設定n個坐標(xi，i＝1，2，3...n)，分別對應於每一對在坐標xi處垂直於X軸的n條直線上的關於X軸上下對稱的接收和發射超聲波的超聲傳感器，並在流動路徑上平行地安裝n對超聲傳感器；每對安裝在坐標xi(i＝1，2，3…n)上下位置的超聲傳感器接收和發送超聲波來測量流體速度Vi(i＝1，2…n)；利用安裝在X軸上點x1和點x2上下位置的超聲傳感器測得的速度V1和V2，應用公式V0＝V1×s1×(s2×V1/V2)2(其中0.7≤s1≤0.8以及1.1≤s2≤1.2)，計算點x＝0處的流體速度V0；利用安裝在X軸上點xn-1和xn上下位置的超聲傳感器測得的速度Vn-1和Vn，應用公式Vf＝Vn×s1×(s2×Vn/Vn-1)2(其中0.7≤s1≤0.8以及1.1≤s2≤1.2)，計算點x＝2r處的流體速度Vf；得到一流體速度的第一函數Vi(x)，在三次樣條函數插值公式中(其中x是可變量)代入xi，Vi(i＝1，2，3，4，5)，以及V0，Vf，可得到如下關於區域xi≤x≤xi+1中的流體區域速度的區域多項式(polynomial)，Vi(x)=ai-1(xi+1-x)36hi+ai(x-xi)36hi+(yihi-ai+1hi6)(xi+1-x)+(yi-1hi-aihi6)(x-xi)]]>其中，hi6ai-1+hi+hi+13ai+hi+16ai+1=yi+2-yi+1hi+1-yi+1-yihi,]]>i＝0，1，2…n-1，a-1＝an(an-1)＝0，hi＝xi+1-xi，以及Vi代表區域xi≤x≤xi+1；中的流體速度；考慮流體流動路徑的管壁在區域0≤x＜xa(這裡，0＜xa＜x1)中影響流體的速度，利用V0，在區域0≤x＜xa中，把x設為可變量，得到流體速度的第二函數Vj(x)＝s3×V0×(1-exp(-x/s42r))，其中，0.2≤s3≤0.4，0.005≤s4≤0.015；考慮流體流動路徑的管壁在區域xb≤x＜2r(這裡，0＜xb＜x2r)中影響流體的速度，利用Vf，在區域xb≤x＜2r中，把x設為可變量，得到流體速度的第三函數Vk(x)＝s3×Vf×(1-exp(-x/s42r))，其中，0.2≤s3≤0.4，0.005≤s4≤0.015；以及在區域0≤x＜xa中，應用第一函數計算流體速度，在區域xa≤x≤xb中，應用第二函數計算流體速度，在區域xb≤x＜2r中，應用第三函數計算流體速度。
xa可位於從點x0(x＝0)開始，沿半徑朝著流動路徑中心方向距離為r/10的地方，而xb可位於從點xf(x＝2r)開始，沿半徑朝著流動路徑中心方向距離為r/10的地方。
i可以是1，2，3，4，和5，安裝5對超聲傳感器，其中，x1位於從點x0(x＝0)開始，沿半徑朝著流動路徑中心方向距離為r/3的地方，x2位於從點x0(x＝0)開始，沿半徑朝著流動路徑中心方向距離為2r/3的地方，x3位於從點x0(x＝0)開始，沿半徑朝著流動路徑中心方向距離為r的地方，x4位於從點xf(x＝2r)開始，沿半徑朝著流動路徑中心方向距離為2r/3的地方，x5位於從點xf(x＝2r)開始，沿半徑朝著流動路徑中心方向距離為r/3的地方。這裡，s1＝0.734959.
s2＝1.163523.
s3＝0.29903.
s4＝0.00958.
上述說明僅是本發明技術方案的概述，為了能夠更清楚了解本發明的技術手段，並可依照說明書的內容予以實施，以下以本發明的較佳實施例並配合附圖詳細說明如後。
參照以下所附圖式，詳細描述本發明的較佳實施例，則不難理解本發明以上和其它特徵和優點。
圖1是描述利用超聲傳感器測量流體速度的方法的示意圖；圖2A是說明流體流動路徑上安裝5對超聲傳感器的示意圖；圖2B是說明利用圖1中5對超聲傳感器測得的速度分布曲線的示意圖；圖3是說明本發明一實施例中流體流動路徑上安裝5對超聲傳感器以及由其測得的速度分布曲線的示意圖；圖4是說明本發明一實施例的速度分布曲線的示意圖；圖5是遵照冪定律的曲線示意圖；圖6是說明超聲傳感器的位置以及傳統加權係數法的加權係數的示意圖；圖7和圖8是說明由傳統加權法獲得的流體流量的誤差率的示意圖；以及圖9說明利用本發明實施例測得的流體速度而獲得的流體流量的誤差率的示意表。
具體實施例方式
舉本發明最佳實施例，結合所附圖式，詳細描述本發明如下。
圖2A是說明流體流動路徑上安裝5對超聲傳感器的示意圖，圖3是說明本發明一實施例中流體流動路徑上安裝5對超聲傳感器以及由其測得的速度分布曲線的示意圖，圖4是說明本發明一實施例的速度分布曲線的示意圖。
參照圖2A、3、和4，現描述本發明一實施例的流體速度測量方法(如下)。
超聲傳感器安裝在流體流動路徑P上，用以測量流體速度。該超聲傳感器安裝在構成流體流動的流動路徑P上，例如是輸油管道、供水和排水管道等的流體管道10。超聲傳感器可以配合敲打插入方式在流體管道10建成後安裝，或當流體管道10正在建設時預先安裝。流體流動管道10的半徑為r，將其一半徑方向設為X軸。由於流體流動管道10的半徑為r，所以X-坐標軸上的變量x的範圍為0≤x≤2r。將X-軸上的n個坐標設為xi，這裡i為1，2，3...n.本實施例設定5個坐標，因此n為1，2，3，4，5，X-坐標分別為x1，x2，x3，x4.x5。坐標x1位於從點x0(x＝0)開始，沿半徑朝著流動路徑中心方向距離為r/3的地方，坐標x2位於從點x0(x＝0)開始，沿半徑朝著流動路徑中心方向距離為2r/3的地方，坐標x3位於從點x0(x＝0)開始，沿半徑朝著流動路徑中心方向距離為r的地方，坐標x4位於從點Xf(x＝2r)開始，沿半徑朝著流動路徑中心方向距離為2r/3的地方，坐標x5位於從點xf(x＝2r)開始，沿半徑朝著流動路徑中心方向距離為r/3的地方。因此超聲傳感器相對於點x3對稱安裝。於是，設定垂直於X-軸並過點x1的直線。所設定的多條直線與水平面的夾角設為θ。每一對安裝在所設定的直線上的超聲傳感器關於X-軸上下對稱，因此所安裝的5對超聲傳感器平行於流體管道10的半徑方向。每一對超聲傳感器由兩個分別安裝在位於x-軸上、下位置的B點和A點的超聲傳感器組成。這兩個超聲傳感器之間的距離為Li，每對超聲傳感器之間的距離分別為L1，L2，L3，L4，L5。例如，超聲傳感器21b和21a分別安裝在對應於點x1的點B和點A，其之間距離為L1。
流體速度Vi(i＝1，2…n)由該5對超聲傳感器21，22，23，24，25測量。組成一對超聲傳感器的兩個超聲傳感器相互對稱，接收和發送穿過流體管道10中的流體的超聲波。更明確地講，安裝在X-軸下方點A的超聲傳感器發射超聲波，安裝在X-軸上方點B的超聲傳感器接收超聲波，反之亦然。測得的超聲波發送和接收之間的每一傳輸時間分別為tAB=LC+Vcos,tBA=LC-Vcos]]>這裡，C是超聲波在靜止的流體中傳播的速度，V是流體流動速度。
相對於流體向前流動的方向向前發射(從A到B)的超聲波的傳輸時間tAB比相對於流體向前流動的方向向後發射(從B到A)的超聲波的傳輸時間tBA較短，其時間差計算如下t=tBA-tAB=2VLcosC2(1-V2cosC2)]]>由於 較小，可忽略不計，流體的平均速度描述如下V=tC22Lcos]]>相應地，安裝在點x1的一對超聲傳感器之間的流體速度V1為V1=tC22L1cos]]>安裝在點x2的一對超聲傳感器之間的流體速度V2為V2=tC22L2cos,]]>安裝在點x3的一對超聲傳感器之間的流體速度V3為V3=tC22L3cos,]]>安裝在點x4的一對超聲傳感器之間的流體速度V4為V4=tC22L4cos,]]>安裝在點x5的一對超聲傳感器之間的流體速度V5為V5=tC22L5cos.]]>點x＝0的虛擬速度，即流過流體流動路徑P的管壁附近的流體虛擬速度V0，由安裝在X-軸上點x1和x2上下位置的每對超聲傳感器21和22測得的速度V1和V2確定。也就是說，速度V0不是一實際值，而是一隨機虛擬值，暫時用於設定樣條函數(稍後將作描述)。這裡，可通過如下公式獲得V0。
___V0＝V1×s1×(s2×V1/V2)2
這裡，s1依流動條件即流體管道10的尺寸或流體條件等而變化，其範圍為0.7≤s1≤0.8.本發明實施例中s1＝0.734959.然而，如果值s1小於0.7或大於0.8，由利用V0設定的樣條函數獲得的速度分布曲線就會與實際的速度分布曲線相差很大。S2依流動條件即流體管道10的尺寸或流體條件等而變化，其範圍為1.1≤s2≤1.2.本發明實施例中s2＝1.163523.然而，如果值s2小於1.1或大於1.2，由利用V0設定的樣條函數獲得的速度分布曲線就會與實際的速度分布曲線相差很大。
點x＝2r的虛擬速度，即流過流體流動路徑P的管壁附近的流體虛擬速度Vf，由安裝在X-軸上點xn-1和xn(本發明實施例中為x4和x5)上下位置的每對超聲傳感器24和25測得的速度Vn-1和Vn確定。也就是說，速度Vf不是一實際值，而是一隨機虛擬值，暫時用於設定樣條函數(稍後將作描述)。這裡，可通過如下公式獲得Vf。
Vf＝Vn×s1×(s2×Vn/Vn-1)2這裡，s1依流動條件即流體管道10的尺寸或流體條件等而變化，其範圍為0.7≤s1≤0.8.本發明實施例中s1＝0.734959.然而，如果值s1小於0.7或大於0.8，由利用Vf設定的樣條函數獲得的速度分布曲線就會與實際的速度分布曲線相差很大。S2依流動條件即流體管道10的尺寸或流體條件等而變化，其範圍為1.1≤s2≤1.2.本發明實施例中s2＝1.163523.然而，如果值s2小於1.1或大於1.2，由利用Vf設定的樣條函數獲得的速度分布曲線就會與實際的速度分布曲線相差很大。
利用流體速度V0到Vf計算流體速度的第一函數Vi(x)。使用三次樣條函數插值得到的曲線(函數)採用區域多項式，以彌補把所有區域作為單一區域的更高次多項式函數的缺陷。
在本發明實施例中，在三次樣條函數插值公式中(其中x是可變量)代入xi，Vi(i＝1，2，3，4，5)，以及V0，Vf，可得到如下關於區域xi≤x≤xi+1中的流體區域速度的區域多項式，Vi(x)=ai-1(xi+1-x)36hi+ai(x-xi)36hi+(yihi-ai+1hi6)(xi+1-x)+(yi+1hi-aihi6)(x-xi)]]>這裡，hi6ai-1+hi+hi+13ai+hi+16ai+1=yi+2-yi+1hi+1-yi+1-yihi]]>i＝0，1，2…n-1，a-1＝an(an-1)＝0，hi＝xi+1-xi，以及Vi代表區域xi≤x≤xi+1中的流體速度。
區域多項式描述如下。
將i＝0代入三次樣條函數插值公式，點x0(x＝0)和點x1(x＝x1)之間的區域多項式V0(x)為
V0(x)=a-1(x1-x)36h0+a0(x-x0)36h0+(y0h0-a1h06)(x1-x)+(y1h0-a0h06)(x-x0),]]>將i＝1代入三次樣條函數插值公式，點x1和點x2之間的區域多項式V1(x)為V1(x)=a0(x2-x)36h1+a1(x-x1)36h1+(y1h1-a2h16)(x2-x)+(y2h1-a1h16)(x-x1).]]>將i＝2代入三次樣條函數插值公式，點x2和點x3之間的區域多項式V2(x)為V2(x)=a1(x3-x)36h2+a2(x-x2)36h2+(y2h2-a3h26)(x3-x)+(y3h2-a2h26)(x-x2),]]>將i＝3代入三次樣條函數插值公式，點x3和點x4之間的區域多項式V3(x)為V3(x)=a2(x4-x)36h3+a3(x-x3)36h3+(y3h3-a4h36)(x4-x)+(y4h3-a3h36)(x-x3),]]>將i＝4代入三次樣條函數插值公式，點x4和點x5之間的區域多項式V4(x)為V4(x)=a3(x3-x)36h4+a4(x-x4)36h4+(y4h4-a5h46)(x5-x)+(y5h4-a4h46)(x-x4).]]>將i＝5代入三次樣條函數插值公式，點x5和點x6之間的區域多項式V5(x)為V5(x)=a4(x4-x)36h5+a5(x-x5)36h5+(y5h5-a6h56)(x6-x)+(y6h5-a5h56)(x-x5)]]>這裡x6＝xf.
如圖3所示，區域多項式具有相同的變量x，並設定為連續，構成第一函數Vi(x)。
考慮流體流動管道10影響流體速度，設定區域0≤x＜xa(這裡，xa為0＜xa＜x1)中的流體速度的第二函數Vj(x)。在區域0≤x＜xa之間，由於流體與流體流動管道的管壁之間的摩擦力大，該區域的流體速度明顯地低於管道中央的流體速度，因此需要新的速度函數。如圖4所示，在區域0≤x＜xa中，用實線表示的新速度函數代替用虛線表示的第一函數Vi(x)。本實施例中，xa位於從點x0(x＝0)開始，沿半徑朝著管道中心方向距離為r/10的地方，即流體流動管道10的直徑的5％處。一般地，流體管道影響的區域大於或小於管道直徑的8％。
利用V0得到的第二函數Vj(x)，其公式如下Vi(x)＝s3×V0×(1-exp(-x/s42r))這裡，exp代表指數函數，s3依流動條件即流體管道10的尺寸或流體條件等而變化，其範圍為0.2≤s3≤0.4.本發明實施例中，s3＝0.29903.然而，如果值s3小於0.2或大於0.4，由利用受流體管道影響的區域的第二函數Vj(x)獲得的速度分布曲線就會與實際的速度分布曲線相差很大。S4依流動條件即流體管道10的尺寸或流體條件等而變化，其範圍為0.005≤s4≤0.015.本發明實施例中，s4＝0.00958.然而，如果值s4小於0.005或大於0.015，由利用受流體管道影響的區域的第二函數Vj(x)獲得的速度分布曲線就會與實際的速度分布曲線相差很大。
考慮流體流動管道10影響流體速度，設定區域xb≤x＜x2r(這裡，xb為0＜xb＜x2r)中的流體速度的第三函數Vk(x)。在區域xb≤x＜x2r之間，由於流體與流體流動管道的管壁之間的摩擦力大，該區域的流體速度明顯地低於管道中央的流體速度，因此需要新的速度函數。如圖4所示，在區域xb≤x＜x2r中，用實線表示的新速度函數代替用虛線表示的第一函數Vi(x)。本實施例中，xb位於從點xf(x＝2r)開始，沿半徑朝著管道中心方向距離為r/10的地方，即流體流動管道10的直徑的5％處。一般地，流體管道影響的區域大於或小於管道直徑的8％。
利用Vf得到的第三函數Vk(x)，其公式如下Vk(x)＝s3×Vf×(1-exp(-x/s42r))這裡，exp代表指數函數，s3依流動條件即流體管道10的尺寸或流體條件等而變化，其範圍為0.2≤s3≤0.4.本發明實施例中，s3＝0.29903.然而，如果值s3小於0.2或大於0.4，由利用受流體管道影響的區域的第二函數Vj(x)獲得的速度分布曲線就會與實際的速度分布曲線相差很大。S4依流動條件即流體管道10的尺寸或流體條件等而變化，其範圍為0.005≤s4≤0.015.本發明實施例中，s4＝0.00958.然而，如果值s4小於0.005或大於0.015，由利用受流體管道影響的區域的第三函數Vk(x)獲得的速度分布曲線就會與實際的速度分布曲線相差很大。
獲得的第二函數Vj(x)、第一函數Vi(x)和第三函數Vk(x)是關於變量x的函數。第一函數Vi(x)由6個函數組成，所以共獲得8個區域函數。由於每一函數有一變量x，如果輸入流體管道10的一x坐標，就可以精確地測得流過在x坐標點處垂直於X-軸的直線的流體的速度。
另外，流體流動管道10的橫截面積沿X-軸(半徑方向)變化，可被表達為一函數。兩次求解該函數與速度函數的乘積的積分，得到流體流量的精確值。
由上述函數產生的流體速度所構成的速度分布曲線與實際速度分布曲線相比，誤差範圍很窄，因此所測得的速度精確，繼之計算出精確的流體流量。
現以應用測得的速度算出的流體流量的誤差來比較傳統加權係數法和本發明實施例的方法。
利用由Honggreung科學出版公司(Hongreung_Science_PublishingCompany)於2003年出版的「流體力學的基本原理」(「Fundamentals_ofFluid_Mechanics」)一書中第415和416頁所描述的冪定律進行比較。該冪定律被廣泛用作均勻對稱流動(這裡指流體速度不隨時間變化)的流體(尤其是湍流)的速度分布曲線的理論模型。已知該冪定律近似等同於在流體流動管道的大多數區域內估計的速度分布曲線。該冪定律表達為VmaxV=(1-rR)1/n]]>該冪定律的曲線圖如圖5所示。其水平坐標軸無量綱，代表離流體管道10的中心的距離x與流體管道半徑r的比率，垂直軸無量綱，代表流體的區域速度v與流體的最大速度Vmax的比率。在這點上，由Honggreung科學出版公司於2003年出版的「流體力學的基本原理」一書中第415和416頁所描述的冪指數n是管壁粗糙度和雷諾數的函數，由以下公式給出，n=1f.]]>該冪指數由以下Colebrook和White方程式計算1f=1.74-2log10(2kD+18.7Ref)]]>其中，Re代表雷諾數，在104到108之間變化，k/D代表相對粗糙度，在0到10-2之間變化。兩者均作為流動模型。
現描述由傳統加權係數法利用上面描述的流動模型所計算的流體流量和其誤差率。
在傳統加權係數法中，5對超聲傳感器21，22，23，24和25沿著流體流動路徑P的半徑方向，按照Gaussian，Chebyshev，和Tailor結構(configuration)配置。對應於該三種結構配置的特定的坐標如圖6所示。圖6的表格參見由韓國流體機械協會於2003年12月5日和6日舉辦的流體機械研究和發展研討會的文獻彙編第157頁。在該表格中，路徑數(thenumber_of_paths)代表多少對超聲傳感器安裝在流體流動管道10上。例如，在Gaussian結構中的路徑數5，表明該例共安裝了5對超聲傳感器。更加詳細地說，一對超聲傳感器安裝於在中心點位置垂直於流體流動路徑P的直線的上、下方，一對超聲傳感器安裝於在中心點的兩邊的任一邊、半徑為Rx_0.538469的點上垂直於流體流動路徑P的直線上，一對超聲傳感器安裝於在中心點的兩邊的任一邊、半徑為Rx_0.906179的點上垂直於流體流動路徑P的直線上，因此共有5對超聲傳感器平行地安裝在流動路徑P上。
將由超聲傳感器分別測得的流體速度V1到V5視為由冪定律獲得的速度值，加權係數W1到W5分別乘以速度V1到V5。因為流體是均勻和對稱流動的，所以與速度值V2和V4分別相乘的加權係數W2和W4均為0.362165，與速度值V1和V5分別相乘的加權係數W1和W5均為0.063786。把這些乘積加起來得到流體平均速度Vmean，在將其與橫截面積A相乘，得到流體的流量Q，如下所示。
Q＝Vmean×A＝(W1V1+W2V2+W3V3+W4V4+W5V5)×A對由傳統加權係數法計算的流體流量Q的誤差率的說明見圖7和圖8的曲線圖。水平軸代表雷諾數和相對粗糙度，垂直軸代表流體的實際流量與公稱流量Qref之間的誤差。改變雷諾數和相對粗糙度，計算流體流量。利用根據理論流動模型由冪定律得到的速度值計算流體的公稱流量Qref。這裡誤差被定義為Error＝{(Q-Qref)/Qref)×100.
圖7所示的流量Q，其路徑數為5，利用了5對超聲傳感器。在Chebyshev結構中，當雷諾數為104，相對粗糙度為10-2時，誤差率為0.2現在描述由本發明實施例測得的速度所計算的流體流量及其誤差。如圖3所示，共安裝5對超聲傳感器21，22，23，24和25，其中，一對安裝於點x1的上、下位置，一對安裝於點x2的上、下位置，一對安裝於點x3的上、下位置，一對安裝於點x4的上、下位置，一對安裝於點x5的上、下位置。由各對超聲傳感器分別測得的速度V1到V5被視為由冪定律獲得的速度值。利用速度V1到V5構建第一函數Vi(x)、第二函數Vj(x)和第三函數Vk(x)。設定點xa和xb對應於流體流動管道10的直徑的5％處，按照以上傳統加權係數法的相同方式變化雷諾數和相對粗糙度，用本發明實施例的方法將橫截面積與第一函數Vi(x)、第二函數Vj(x)和第三函數Vk(x)相乘，並將乘積二次積分得到流量。流量Q與公稱流量Qref之間的誤差率見圖9所示的曲線圖。利用根據理論流動模型由冪定律得到的速度值，來計算流體的公稱流量Qref。誤差率的計算方式與傳統加權係數法誤差率的計算方式相同。參見圖9，當雷諾數為104，相對粗糙度為10-2時，誤差率為0.015，當雷諾數為108，相對粗糙度為10-5時，誤差率為-0.047。這些誤差率低於傳統加權係數法的誤差率。
從圖7到圖9，可以看出，不管雷諾數和相對粗糙度，由本發明實施例利用所測得的速度計算的流量的誤差率明顯低於由傳統加權係數法計算的流量誤差率。
如上所述，根據本發明的一實施例所測量流體速度的方法，因考慮流體的流動特徵，將流動路徑的橫截面分成三部分，並為之構建反映流體在流動路徑中流動的實際速度分布曲線的合適的速度函數。因而提供精確的流體速度值，從而顯著降低流量的計算誤差。
以上所述，僅是本發明的較佳實施例而已，並非對本發明作任何形式上的限制，雖然本發明已以較佳實施例揭露如上，然而並非用以限定本發明，任何熟悉本專業的技術人員，在不脫離本發明技術方案範圍內，當可利用上述揭示的結構及技術內容作出些許的更動或修飾為等同變化的等效實施例，但是凡是未脫離本發明技術方案的內容，依據本發明的技術實質對以上實施例所作的任何簡單修改、等同變化與修飾，均仍屬於本發明技術方案的範圍內。
權利要求
1.一種測量流體流動速度的方法，其特徵在於其包括以下步驟設定半徑為r的流體流動路徑的一半徑方向為X-軸，在X-軸(0≤x≤2r)上設定n個坐標(xi，i＝1，2，3...n)，將每一對接收和發送超聲波的超聲傳感器對稱地配置於在xi坐標處垂直於X-軸的n條直線上的位於X-軸上方和下方的位置，在流動路徑上所安裝的n對超聲傳感器相互平行；每一對安裝在坐標xi(i＝1，2，3...n)處的上下位置的超聲傳感器，通過接收和發送超聲波來測量流體速度Vi，(i＝1，2，3...n)；利用位於X-坐標上點x1和點x2的上、下位置的每對超聲傳感器測得的速度V1和V2，由一公式V0＝V1×s1×(s2×V1/V2)2計算點x＝0的流體速度V0，其中0.7≤s1≤0.8以及1.1≤s2≤1.2；利用位於X-坐標上點xn-1和點xn的上、下位置的每對超聲傳感器測得的速度Vn-1和Vn，由一公式Vf＝Vn×s1×(s2×Vn/Vn-1)2計算點x＝2r的流體速度Vf，其中0.7≤s1≤0.8以及1.1≤s2≤1.2；在三次樣條函數插值公式中(其中x是可變量)代入xi，Vi(i＝1，2，3，4，5)，以及V0，Vf，得到關於在區域xi≤x≤xi+1中的流體區域速度的區域多項式，從而獲得流體速度的第一函數Vi(x)Vi(x)=ai-1(xi+1-x)36hi+ai(x-xi)36hi+(yihi-ai+1hi6)(xi+1-x)+(yi+1hi-aihi6)(x-xi)]]>這裡，hi6ai-1+hi+hi+13ai+hi+16ai+1=yi+2-yi-1hi+1-yi+1-yihi]]>i＝0，1，2...n-1，a-1＝an(an-1)＝0，hi＝xi+1-xi，以及Vi代表區域xi≤x≤xi+1中的流體速度；考慮流體流動路徑的管壁在區域0≤x＜xa(這裡，0＜xa＜x1)中影響流體的速度，利用V0，在區域0≤x＜xa中，把x設為可變量，得到流體速度的第二函數Vj(x)＝s3×V0×(1-exp(-x/s42r))，其中，0.2≤s3≤0.4，0.005≤s4≤0.015；考慮流體流動路徑的管壁在區域xb≤x＜2r(這裡，0＜xb＜x2r)中影響流體的速度，利用Vf，在區域xb≤x＜2r中，把x設為可變量，得到流體速度的第三函數Vk(x)＝s3×Vf×(1-exp(-x/s42r))，其中，0.2≤s3≤0.4，0.005≤s4≤0.015；以及在區域0≤x＜xa中，應用第一函數計算流體速度，在區域xa≤x≤xb中，應用第二函數計算流體速度，在區域xb≤x＜2r中，應用第三函數計算流體速度。
2.根據權利要求1所述的方法，其特徵在於，xa可位於從點x0(x＝0)開始，沿半徑朝著流動路徑中心方向距離為r/10的地方，而xb可位於從點xf(x＝2r)開始，沿半徑朝著流動路徑中心方向距離為r/10的地方。
3.根據權利要求1所述的方法，其特徵在於，i可以是1，2，3，4，和5，安裝5對超聲傳感器，其中，x1位於從點x0(x＝0)開始，沿半徑朝著流動路徑中心方向距離為r/3的地方，x2位於從點x0(x＝0)開始，沿半徑朝著流動路徑中心方向距離為2r/3的地方，x3位於從點x0(x＝0)開始，沿半徑朝著流動路徑中心方向距離為r的地方，x4位於從點xf(x＝2r)開始，沿半徑朝著流動路徑中心方向距離為2r/3的地方，x5位於從點xf(x＝2r)開始，沿半徑朝著流動路徑中心方向距離為r/3的地方。
4.根據權利要求1所述的方法，其特徵在於s1＝0.734959。
5.根據權利要求1所述的方法，其特徵在於s2＝1.163523。
6.根據權利要求1所述的方法，其特徵在於s3＝0.29903。
7.根據權利要求1所述的方法，其特徵在於s4＝0.00958。
全文摘要
本發明提供一種測量流體速度的方法，該方法考慮流體的流動特徵，將一流動路徑的橫截面分成三部分，並構建分別適合於該三部分的函數，即第一函數
文檔編號G01F1/66GK1734267SQ20051005997
公開日2006年2月15日 申請日期2005年4月4日 優先權日2004年8月11日
發明者金秉燦 申請人:昌民Tech株式會社