用於通過遠程測量徑向速度確定流體移動的方法和裝置製造方法

2023-09-23 06:55:50

用於通過遠程測量徑向速度確定流體移動的方法和裝置製造方法
【專利摘要】本發明涉及用於確定目標體積中流體的流動的方法，其包括在沿至少三個測量軸線分布的多個測量點處遠程測量該測量點附近的流體的徑向速度的遠程測量的步驟，該測量軸線穿過目標體積並具有不同的空間定向，並且該方法還用於計算分布為目標體積中的網格的多個計算點處的流體速度，其中流體速度的計算包括利用流體力學特性模型。本發明還涉及一種實施該方法的裝置。
【專利說明】用於通過遠程測量徑向速度確定流體移動的方法和裝置
【技術領域】
[0001]本發明涉及用於通過遠程測量流體的徑向移動速度來確定流體的移動的方法。本發明還涉及通過實施該方法遠程感測流體移動的裝置。更具體地，本發明所屬的領域涉及但不限於位於大氣底層的風的特性的遠程感測。
【背景技術】
[0002]對於許多應用，測量大氣或風的移動相當重要，尤其是在氣象學和用於監控諸如機場和風力發電廠的特徵的應用中。
[0003]通常有用的是，在很寬的高度範圍內或與大體積相對應的區域中測量空氣團的位移，而這是利用傳統的風速計如轉杯風速計無法實現的，並且這樣的測量需要能夠進行遠程測量的遠程感測儀器。
[0004]這些儀器具體包括雷達、雷射雷達和聲雷達。雷達和雷射雷達系統分別使用超高頻和光頻範圍中的電磁波。聲雷達系統使用聲波。
[0005]具體地，在由A.Dolf1-Bouteyre, M.Valla, B.Augere, J.-P.Cariou, D.Goular,D.Fleury, G.Canat, C.Planchat, T.Gaudo, L.Lombard, 0.Petilon, J.Lawson-Daku 於2007年7月8日至13日在美國科羅拉多州，Snowmass舉行的第14屆相干雷射雷達會議(CLRCXIV)上發表的「用於尾渦流監控的1.5 μ m全纖維脈衝雷射雷達」中描述了實施雷射雷達的示例。
[0006]利用這種儀器測量空氣團的移動通常以如下的方式執行:
[0007]-設備在待測區域中沿發射軸線連續地或以脈衝的形式發射一束或多束波(聲學波和/或電磁波)。沿不同發射軸線的發射可以同時或者順序進行。
[0008]-波束在大氣中受到散射效應,尤其由於所遇到的不均勻性(大氣微粒、顆粒、對於電磁波折射率的變化或對於聲波聲阻的變化)而受到散射效應。當波束在空氣團或移動粒子中散射時，這些波束由於都卜勒效應還會發生頻移。
[0009]-背散射波束通過根據測量軸線定向的一個或多個接收器來檢測。這些接收器在沿其測量軸線的方向上檢測由大氣散射的波。
[0010]-因此可以沿檢測器的測量軸線計算發生散射的距離，例如通過測量傳播時間的方法或通過幹涉測量法測量相移的方法。沿測量軸線的空氣團或粒子的徑向速度還可以通過測量由於都卜勒效應的波的頻移而獲得。所測量的徑向速度與散射區域的速度矢量在檢測器的測量軸線上的投影相對應。
[0011]遠程感測儀器、尤其包括適於測量位於大氣底層的風的特徵的雷射雷達系統的遠程感測儀器常為單站類型。這表示相同的光學器件或相同的天線(聲學的或電磁的)用於發送和接收信號。所探測的體積通常沿錐體分布，該錐體的頂點位於儀器的光學器件或天線的水平上。儀器沿錐體的每束脈衝測量粒子沿與發射軸線重合的測量軸線移動的徑向速度。因此獲得風的徑向速度的測量值，其表示風矢量在射束傳播軸線上的投影。
[0012]然後基於所獲得的徑向速度的測量值，計算整個目標體積的風矢量。[0013]在現有的裝置中，通常利用純幾何模型來執行這種計算。這些模型具有這樣的缺點，即這些模型基於有時相當不現實的假設，尤其是對於整個樣品的測量過程中風的空間和時間均勻性。根據該假設，在給定高度上，每個由儀器探測的大氣點處的風矢量都是相同的。
[0014]已知，例如可以通過求解具有三個未知數的具有至少三個方程的方程組來根據在同一高度處的至少三個方向上測量的至少三個徑向速度測量值計算給定高度處的風矢量的分量，其中所述方程組描述風矢量與其沿測量軸線的投影的幾何關係，該測量軸線通過徑向速度的測量而構成。
[0015]一種用於幾何計算的方法為人們所知，其在文獻中被稱為「速度方位顯示(VAD)」，該方法也是基於給定高度處風的空間均勻性的假設。如果裝置允許在覆蓋整個測量錐體的方向上進行測量，則該方法也是適用的。該方法包括用於豎直速度、方向以及水平速度的參數優化，充分利用在測量錐體的一整周(360度)上測量的所有徑向速度呈現調和函數形式的事實。該方法特別適用於處理利用雷射雷達系統獲得的測量值。
[0016]已知Werner等人的文件US4, 735, 503描述了來源於VAD方法的一種方法,根據該方法採集可被限制於成角度的扇形中。然而，該方法基於給定高度處風的空間均勻性的相同假設。
[0017]當在基本平坦的地形(具有很小起伏或沒有起伏的地形或者海面)上進行測量時，用於利用再建風矢量的幾何技術測量風的大氣遠程感測儀器允許精確測量風的平均速度。例如，利用雷射雷達系統，在10分鐘內取平均的測量值相對於通過校準的轉杯風速計構成的參考值的相對誤差低於2%。
[0018]另一方面，當在更複雜的地形如起伏的地形、山地、或覆蓋有森林的地形等上進行測量時，確定水平和豎直速度的準確性顯著降低。對於利用相同的雷射雷達系統執行的測量並且相對於校準的轉杯風速計，在複雜地點上觀察到的在10分鐘內計算的平均值的相對誤差約為5%至10%。
[0019]因此現有的實施幾何模型的裝置不能十分精確地測量水平和豎直速度，以及複雜地形上的風向。事實上，在這種特定的情況下，在由儀器所探測的大氣體積中的給定高度處的風不能再認為是均勻的。現在，在這些情況下、尤其是在發展風力發電廠的背景下，精確測量風的特徵是必需的。
[0020]當研究與高度有關的風速輪廓線(B卩，風的速度矢量與直接在儀器之上的高度有關)時也會出現該限制。
[0021]本發明的目的是提出適於在複雜環境中測量氣團移動的用於確定流體移動的方法。

【發明內容】

[0022]本發明的目的通過利用用於在目標體積中確定流體流動的方法來實現，該方法包括以下步驟:
[0023]-在沿至少三個測量軸線分布的多個測量點處，遠程測量測量點附近的流體的徑向速度，該測量軸線穿過目標體積且具有不同的空間定向；
[0024]-計算多個計算點處的流體的速度，多個計算點分布為目標體積中的網格；[0025]其中，計算流體的速度包括使用流體的力學特性模型(mechanical behaviormodel)。
[0026]有利地，測量流體的徑向速度可包括測量之前發射且在流體中散射的波的(由於都卜勒效應產生的)頻移。
[0027]徑向速度的測量可利用具有各種配置的儀器來執行，尤其是單站儀器(發送器和接收器在空間上一致)、雙站儀器(包括在空間上獨立的發送器和接收器)或在一個場所分布的多個儀器的組合。
[0028]當然根據本發明的方法適用於任何類型的流體的測量，尤其是液體或氣體，在流體中可獲得徑向速度的測量值，並且流體可以通過力學特性模型來描述。
[0029]該方法使得可以通過在多個測量軸線上得到的徑向速度的測量值來確定目標體積中流體的移動，該方法比已知的幾何方法更精確，尤其是在當體積中流體的行為不滿足現有技術的幾何模型所必需的空間均勻性的假設的情況下。
[0030]例如該假設可描述如下:表示流體局部移動的速度矢量對於目標體積的所有點大體上相同，這些點位於相對於參考物(如傳感器)大體上相等的距離處。而該假設實際上經常被證明是相當不切實際的。
[0031]使用流體行為模型使得可以描述更複雜的行為。因此可以大大提高流體速度矢量的分量的測量，尤其是在複雜的和/或有幹擾的環境中。
[0032]有利地，在根據本發明的方法直接提供一定體積中的風的速度矢量的三維表示的意義上講，根據本發明的方法本質上是三維的。此外，假定沿具有相對於目標體積合理分布的空間定向的測量軸線執行徑向速度的測量，流體速度矢量的分量都以相當的準確性來確定，例如水平和豎直分量之間沒有特別的區別。
[0033]目標體積可由連接計算點 的網格確定邊界或限定。該網格可包括結構化的(例如具有矩形或曲線截面的笛卡爾網格)或非結構化的網。測量點可包括在目標體積中並包括在網格中。
[0034]根據本發明的方法還可包括計算初始條件的步驟，該步驟包括:利用基於假設(目標體積中的流體速度在測量軸線通過的基本上平行定向的層中基本均勻)的幾何模型，通過徑向速度的測量值來計算位於計算點處的流體速度，計算初始條件包括以下步驟中的至少之一:
[0035]-邊界條件的計算，包括計算限制位於目標體積邊緣的計算點處的流體速度(或速度矢量)的限定條件，以及
[0036]-起始條件的計算,包括計算目標體積中計算點處的流體速度。
[0037]在目標體積中定義笛卡兒坐標系(X，Y，Z)，其中Z軸被定義為高度，均勻層可例如定向成基本平行於平面(Χ，γ)。
[0038]根據本發明的方法還可包括使用存在於目標體積的邊緣中或邊緣處材料表面的拓撲結構計算初始條件並且限制流體的延伸的步驟，該步驟包括:
[0039]-確定位於流體之外和/或材料表面附近的目標體積的計算點，以及
[0040]-歸因分析限制所述被限定計算點的速度(或速度矢量)的條件。
[0041]該初始條件可考慮到出現並且影響流體流動的表面的起伏。材料表面可以是例如位於目標體積底部處的地形表面，其拓撲結構已經預先確定。為了能夠應用該初始條件，目標體積及其網格被限定以使得全部或部分地包括該表面。
[0042]根據本發明的方法還可包括通過使用先前計算的初始條件並且通過求解流體的力學特性模型的方程來計算目標體積中的流體速度(或速度矢量)的步驟。
[0043]僅初始條件(包括邊界條件和起始條件)根據空間均勻性的假設進行計算。初始條件構成計算目標體積中流體速度的起始點:
[0044]-邊界條件使得可以對目標體積邊緣處的流體速度設立約束，從而求解流體力學特性模型的方程以在目標體積中得到確定的解；
[0045]-起始條件使得可以初始化以通過迭代法求解流體的力學特性模型的方程，從而使方程朝向確定的解收斂。
[0046]在本發明的背景中，可以實施流體力學特性的所有模型。取決於背景，模型的選擇可以有利地改變，以優化計算的複雜性。
[0047]根據本發明的實施方式，流體力學特性模型可包括這樣的假設，S卩，流體包括不可壓縮的牛頓流體並且其流動可近似地由納維-斯託克斯方程(Navier-Stokes equation)描述。
[0048]納維-斯託克斯方程(根據其一種通常的表現形式)可包括連續性方程和動量守恆方程。
[0049]根據本發明的實施方式，流體力學特性模型可包括流體的流動是靜止的假設。
[0050]根據其他實施方式，流體力學特性模型可包括這樣的假設，S卩，流體包括理想流體並且其流動可近似地由流體的歐拉方程來描述。
[0051]在這種情況下，和不可壓縮的牛頓流體相比，熱梯度和粘性的效果也可以忽略。
[0052]牛頓流體為粘性應力張量是應變張量的線性函數的流體。該模型很好地適用於多種常見流體，具體包括空氣和水。定義以下變量，採用以國際單位制表示的單位:
[0053]V ；流體的速度矢量(通常為m.S4);
[0054]V:流體的運動粘度(m2.S-1)；
[0055]P:流體的密度(kg.m_3);
[0056]/:作用在流體上的質量力的合成(N.kg—1)；
[0057]P:流體的壓強(Pa)。
[0058]因此可以通過包括質量守恆方程和動量守恆方程的方程組來描述流動:
【權利要求】
1.一種用於確定目標體積中流體的流動的方法，包括以下步驟: -在沿至少三個測量軸線(3)分布的多個測量點(21)處，遠程測量位於所述測量點(21)附近的所述流體的徑向速度(11)，所述測量軸線穿過所述目標體積並且具有不同的空間定向， -計算多個計算點(22、23)處的所述流體的速度(10)，其中所述多個計算點分布為所述目標體積中的網格(20)； 其中，計算所述流體的速度(10)包括利用所述流體的力學特性模型。
2.根據權利要求1所述的方法，其中，測量所述流體的徑向速度(11)包括:測量先前被發送並且在所述流體中被散射的波的都卜勒頻移。
3.根據權利要求1或2所述的方法，還包括計算初始條件的步驟，所述計算初始條件的步驟包括:利用基於假設的幾何模型，通過徑向速度(11)的測量值計算位於計算點(22，23)處的所述流體的速度(10)，所述假設為所述目標體積中所述流體的速度(10)在所述測量軸線(3)穿過的多個層中大體上是均勻的，所述多個層具有大體上平行的定向，所述初始條件的計算包括以下計算中的至少之一: -邊界條件的計算，包括對限制位於所述目標體積邊緣的計算點(23)處的所述流體的速度(10)的條件的計算，以及 -起始條件的計算，包括位於所述目標體積中計算點(22)處的所述流體的速度(10)的計算。
4.根據權利要求3所述的方法，還包括利用材料表面的拓撲結構計算初始條件的步驟，所述材料表面存在於所述目 標體積的邊緣中或存在於所述目標體積的邊緣處，並且限定所述流體的延伸，所述利用材料表面的拓撲結構計算初始條件的步驟包括: -確定位於所述流體之外和/或位於所述材料表面附近的所述目標體積的計算點(22、23)，以及 -歸因分析所限定的計算點(22、23)的速度(10)的限制條件。
5.根據權利要求3或4所述的方法，還包括通過利用先前計算的初始條件並且通過求解所述流體的力學特性模型的方程來計算所述目標體積中所述流體的速度(10)的步驟。
6.根據權利要求5所述的方法，其中所述流體的力學特性模型包括以下假設中任何一個: -所述流體包括不可壓縮的牛頓流體並且其流動近似由納維-斯託克斯方程描述，或者 -所述流體包括理想流體並且其流動近似由流體的歐拉方程描述。
7.根據權利要求6所述的方法，其中f表示所述流體的速度矢量(10)，V表示所述流體的運動粘度，P表示所述流體的密度，/表示作用在所述流體上的質量力的合力以及P表示所述流體的壓強，所述流體的流動由以下方程組中任何一個描述:
V-F = O
---1 -.V * I ' + (?, -Vl1- -1- -Vp 二 tOtP
8.根據權利要求7所述的方法，其中所述流體的力學特性模型的方程在所述目標體積的計算點處通過迭代法數值地求解。
9.根據權利要求8所述的方法，其中求解所述力學特性模型的方程包括利用條件矩陣。
10.根據權利要求5所述的方法，其中所述流體的所述力學特性模型包括一種假設，假設所述流體包括理想流體，並且所述流體的流動在所述目標體積中是近似無旋的。
11.根據權利要求10所述的方法，其中P表示所述流體的速度矢量(10)，P表示速度勢，所述流體的流動由以下方程描述:
12.根據權利要求11所述的方法，其中所述流體的所述力學特性模型的方程通過在所述目標體積的計算點處對拉普拉斯算子的矩陣進行求逆來數值地求解。
13.根據權利要求5-12中任一項所述的方法，還包括以下步驟: 通過求解所述流體的所述力學特性模型的方程從先前在所述目標體積中計算的所述流體的速度(10)來計算新的初始條件， 利用所述新的初始條件，通過求解所述流體的所述力學特性模型的方程來計算所述目標體積中的所述流體的速度(10)。
14.根據上述權利要求中任一項所述的方法，還包括計算所述目標體積中所述流體的速度(10)和測量徑向速度(11)的至少一個初步序列，並且計算所述目標體積中的所述流體的速度(10)還包括利用流體力學特性的動態模型，以及利用在一個或多個所述初步序列期間計算的所述流體速度(10)。
15.根據上述權利要求中任一項所述的方法，所述方法被實施以用於測量位於大氣底層中的風。
16.—種通過實施根據上述權利要求中任一項所述的方法來確定目標體積中流體的流動的裝置，包括: -用於遠程測量的部件(1)，在沿至少三個測量軸線(3)分布的多個測量點(21)處遠程測量位於所述測量點(21)附近的所述流體的徑向速度(11)，所述測量軸線(3)穿過所述目標體積並且具有不同的空間定向， -用於計算的部件(9)，用於計算多個計算點(22，23)處的所述流體的速度(10)，所述多個計算點分布為所述目標體積中的網格(20)；其中，所述用於計算所述流體的速度(10)的部件(9)被設置成用於實施所述流體的力學特性模型。
17.根據權利要求16所述的裝置，還包括以下裝置中的任何一個:雷射雷達、雷達、聲雷達。`
【文檔編號】G01P5/26GK103430030SQ201180064405
【公開日】2013年12月4日 申請日期:2011年12月28日 優先權日:2011年1月5日
【發明者】馬希歐·博奎特, 馬克西米·尼巴特, 阿蒙德·阿爾伯格爾 申請人:利奧斯費爾公司, 阿里亞技術公司