JP2022502675A - パイプ内の軸方向の流速分布と流量を音波法で測定する方法、およびシステム - Google Patents

パイプ内の軸方向の流速分布と流量を音波法で測定する方法、およびシステム Download PDF

Info

Publication number
JP2022502675A
JP2022502675A JP2021531155A JP2021531155A JP2022502675A JP 2022502675 A JP2022502675 A JP 2022502675A JP 2021531155 A JP2021531155 A JP 2021531155A JP 2021531155 A JP2021531155 A JP 2021531155A JP 2022502675 A JP2022502675 A JP 2022502675A
Authority
JP
Japan
Prior art keywords
sound wave
flow velocity
pipe
section
measurement
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2021531155A
Other languages
English (en)
Other versions
JP7032842B2 (ja
Inventor
言欽 李
帥杰 孫
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Publication of JP2022502675A publication Critical patent/JP2022502675A/ja
Application granted granted Critical
Publication of JP7032842B2 publication Critical patent/JP7032842B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G01P5/24Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave
    • G01P5/245Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave by measuring transit time of acoustical waves
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01FMEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/66Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by measuring frequency, phase shift or propagation time of electromagnetic or other waves, e.g. using ultrasonic flowmeters
    • G01F1/662Constructional details
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01FMEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/66Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by measuring frequency, phase shift or propagation time of electromagnetic or other waves, e.g. using ultrasonic flowmeters
    • G01F1/667Arrangements of transducers for ultrasonic flowmeters; Circuits for operating ultrasonic flowmeters
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01FMEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F7/00Volume-flow measuring devices with two or more measuring ranges; Compound meters
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H5/00Measuring propagation velocity of ultrasonic, sonic or infrasonic waves, e.g. of pressure waves
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G01P5/001Full-field flow measurement, e.g. determining flow velocity and direction in a whole region at the same time, flow visualisation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G01P5/24Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G01P5/24Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave
    • G01P5/241Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave by using reflection of acoustical waves, i.e. Doppler-effect
    • G01P5/242Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave by using reflection of acoustical waves, i.e. Doppler-effect involving continuous, e.g. modulated or unmodulated, waves

Abstract

【課題】本発明は、パイプ内の軸方向の流速分布と流量を音波法で測定する方法、およびシステムを提供する。【解決手段】流速分布と流量を測定する方法は、パイプまたはコンテナの安定した区間または関心のある区間で一意の横断面を選択し、横断面の周りの壁に複数の音波センサを設置するステップであって、複数の音波センサ同士の間は、複数の有効な音波経路を形成し、ここで、各々のペアの音波センサ同士の間は、1つの有効な音波経路としてのみカウントされるステップと、各々の有効な音波経路での音波伝播時間をそれぞれ測定して、音波経路の屈折原理に基づいて、各々の音波伝播時間を本発明のモデル式に代入して再構築計算を実行することによって、測定待ちのパイプまたはコンテナの横断面の軸方向の流速分布u(x,y)を得ることができ、得られた軸方向の流速分布u(x,y)を当該横断面に沿って積分して、パイプ内の流量を得る。【選択図】図2

Description

本発明は、流体測定の分野に関し、具体的には、パイプ内の軸方向の流速分布と流量を音波法で測定する方法、およびシステムに関する。
パイプまたはコンテナ内の軸方向の流速分布と流量は、多くの産業および科学研究での応用において重要またはキーポイントの技術パラメータであるとともに、一般的に軸方向の流速は不均一に分布するため、効果的な測定方法が必要である。音波法は、適用可能かつ高度な測定方法である。
音の速度と音波伝播方向に平行な流速成分とが直接重畳されて音の伝播時間に影響を与えるという原理に基づいて、従来の音波法の測定ではパイプの軸方向に沿って互いに異なる位置の2つの断面に一定の数の音波センサをそれぞれ設置する必要があり、このようにして得られるものは、2つの断面間の軸方向の流速分布情報の平均値である。従来の音波法の測定は、次のような欠点を有する、すなわち、一面では、必要とされる音波センサの数が多く、システムが複雑であり、他の一面では、測定するパイプ条件に対する要件が比較的に厳しく、また、軸方向の流速によって引き起こされる音の経路の屈折がさまざまな程度で測定の正確性および精度に影響を与える。
以上の存在する問題を解決するために、人々は、常に、理想的な技術的解決策を模索してきた。
本発明の目的は、従来技術の欠点を克服して、パイプ内の軸方向の流速分布と流量を音波法で測定する方法、およびシステムを提供することにある。
上述した目的を実現するために、本発明が採用する技術的解決策のパイプ内の軸方向の流速分布と流量を音波法で測定する方法は、
パイプまたはコンテナの安定した区間または関心のある区間において一意の横断面を選択し、前記横断面の周りの壁に複数の音波センサを設置するステップであって、複数の音波センサ同士の間は、複数の有効な音波経路を形成し、ここで、各々のペアの音波センサ同士の間は、1つの有効な音波経路としてのみカウントされるステップと、
各々の有効な音波経路での音波伝播時間をそれぞれ測定して、以下のモデル式に代入して、軸方向の流速分布u(x,y)を再構築して得るステップであって、
Figure 2022502675
ここで、lは、i番目の音波経路であり、Lは、i番目の有効な音波経路での2つの音波センサ同士の間の直線距離であり、Δtは、i番目の有効な音波経路での音波伝播時間であり、Nは、有効な音波経路の数であり、cは、媒体の物理的特性の条件下の音波伝播速度であるステップと、
前のステップで得られた軸方向の流速分布u(x,y)を当該横断面に沿って積分して、パイプ内の流量を得るステップと、を備える。
上述した方法において、各々の有効な音波経路での2つの音波センサ同士の間の向かい合っている方向の2つの音波伝播時間をそれぞれ測定し、この2つの音波伝播時間の平均値を当該有効な音波経路での音波伝播時間とする。
上述した方法において、各々の音波センサは、順番に音を発し、任意の1つの音波センサによって音波が発射されるときに、他のすべての音波センサは、いずれも当該音波を記録し、すべての音波センサによって発射される音波の周波数は、いずれも同一である。
上述した方法において、2〜3個ごとの音波センサは、同時に音を発し、ここで、各々の音波センサによって発射される音波の周波数は、互いに異なっており、フィルタリングによってこれら互いに異なっている音波の周波数を効果的に識別する。
上述した方法において、前記音波センサは、発受信一体型の音波センサである。
上述した方法において、前記音波センサは、それぞれ別体である音波発射センサと音波受信センサとを組み合わせたものである。
上述した方法において、前記軸方向の流速分布u(x,y)は、Taylor級数多項式を採用してフィッティング計算される。
上述した方法において、前記軸方向の流速分布u(x,y)は、自由射流の特性に一致する場合、以下のガウス多項式を採用して近似フィッティング計算され、
Figure 2022502675
ここで、(x,y)は、自由射流の流速の極値点の座標であり、Uは、自由射流の流速の極値であり、uは、自由射流中心から離れた端部での流速値であり、σは、自由射流の膨張特性の幅である。
上述した方法において、選択された横断面は、パイプ軸またはコンテナ軸と直角またはほぼ直角に交差する。
本発明は、上述したパイプ内の軸方向の流速分布と流量を音波法で測定する方法のシステムをさらに提供し、当該システムは、パイプまたはコンテナの特定の横断面の周りの壁に沿って分布する複数の音波センサと、デジタル/アナログ変換カードと、アナログ/デジタル変換カードと、測定用のコンピューターと、を備え、前記測定用のコンピューターには、上述したパイプ内の軸方向の流速分布と流量を音波法で測定する方法によって構築された測定ソフトウェアが事前設定されており、
前記デジタル/アナログ変換カードは、それぞれ前記測定用のコンピューターおよび前記音波センサに接続され、前記測定用のコンピューターの測定ソフトウェアが発したオーディオデジタル信号を音波アナログ信号に変換し、前記音波センサを介して発信するために使用され、
前記アナログ/デジタル変換カードは、それぞれ前記測定用のコンピューターおよび前記音波センサに接続され、前記音波センサが収集したパイプ内の音波情報をデジタル信号に変換して前記測定用のコンピューターに入力するために使用され、
前記測定用のコンピューターは、前記測定ソフトウェアを利用してすべての音波センサが音を発するように制御し、各々の音波センサからその他の音波センサへの有効な音波経路での音波伝播時間を測定し、各々の有効な音波経路での音波伝播時間をパイプまたはコンテナの横断面における軸方向の流速分布の再構築式
Figure 2022502675
に代入することによって、軸方向の流速分布u(x,y)を得ることができ、ここで、lは、i番目の有効な音波経路であり、Lは、i番目の有効な音波経路での2つの音波センサ同士の間の直線距離であり、Δtは、i番目の有効な音波経路での音波伝播時間であり、Nは、有効な音波経路の数であり、cは、媒体の物理的特性の条件下の音波伝播速度であり、
前記測定用のコンピューターは、前記測定ソフトウェアを利用して、得られた軸方向の流速分布u(x,y)を当該測定横断面に沿って積分して、パイプ内の流量を得る。
本発明は、従来技術と比較して、卓越した実質的特徴と顕著な進歩を有している。具体的には、本発明は、単一の横断面に配置した各音波センサ同士の間の音波伝播時間データを根拠として、音波経路の屈折理論に基づいてパイプまたはコンテナ内の軸方向の流速分布を再構築する方法によって、パイプまたはコンテナ内の軸方向の流速分布と流量を測定する。従来の音波法の測定と比較すると、本発明は、同じ精度の測定結果の要件を満たす上で、上述した音波センサの数が従来よりも2分の1と少なく、また、誤差が小さいし、精度が高く、信頼性が高い。さらに、本発明の測定精度は、横断面がパイプの軸に厳密に直角であるかどうか、音波センサが横断面の周りに沿って均一に配置されたかどうか、および、軸方向の流速場を再構築するモデル式の適切な変形などによって影響を受けない。
本発明に開示される方法は、気体、液体、および二相または多相が流れているパイプ内の軸方向の流速分布と流量に対する有効な測定に適用することができ、燃焼室、自由射流場、流動床、化学反応器などの内部流れ場や、開放空間での自由射流などの対象の測定に拡張することもできる。
1次元のx方向の軸方向の流速分布および該当する音線の分布の模式図である。 本発明の方法における配置の模式図である。 流れ場中心が幾何学的中心に位置する円形断面パイプ内の軸方向の流れ場のプリセットシミュレーションフィールドを示す。 流れ場中心が幾何学的中心に位置する、Taylor級数フィッティング3次精度の場合の6組の音波センサの円形断面パイプ内の軸方向の流れ場に対するシミュレーション測定による再構築の相対誤差REを示す。 流れ場中心が幾何学的中心に位置する、Taylor級数フィッティング4次精度の場合の6組の音波センサの円形断面パイプ内の軸方向の流れ場に対するシミュレーション測定による再構築の相対誤差REを示す。 流れ場中心が幾何学的中心から偏向する円形断面パイプ内の軸方向の流れ場のプリセットシミュレーションフィールドを示す。 流れ場中心が幾何学的中心から偏向する、Taylor級数フィッティング3次精度の場合の6組の音波センサの円形断面パイプ内の軸方向の流れ場に対するシミュレーション測定による再構築の相対誤差REである。 流れ場中心が幾何学的中心から偏向する、Taylor級数フィッティング4次精度の場合の6組の音波センサの円形断面パイプ内の軸方向の流れ場に対するシミュレーション測定による再構築の相対誤差REを示す。 流れ場中心が幾何学的中心から偏向する、Taylor級数フィッティング3次精度の場合の8組の音波センサの円形断面パイプ内の軸方向の流れ場に対するシミュレーション測定による再構築の相対誤差REを示す。 流れ場中心が幾何学的中心から偏向する、Taylor級数フィッティング4次精度の場合の8組の音波センサの円形断面パイプ内の軸方向の流れ場に対するシミュレーション測定による再構築の相対誤差REを示す。
以下、具体的な実施形態によって、本発明の技術的解決策をいっそう詳細に説明する。
図2に示すように、パイプ内の軸方向の流速分布と流量を音波法で測定する方法は、以下のステップを含む。
パイプまたはコンテナの安定した区間または関心のある区間において一意の横断面を選択し、前記横断面の周りの壁に複数の音波センサを設置する。前記音波センサの数は、流れ場の複雑さに基づいて設置することができ、好ましくは4〜8個の音波センサを設置する。前記音波センサは、発受信一体型の音波センサでもよいし、また、それぞれ別体である音波発射センサと音波受信センサとを組み合わせたものであってもよい。
複数の音波センサ同士の間は、複数の有効な音波経路を形成し、ここで、各々のペアの音波センサ同士の間は、1つの有効な音波経路としてのみカウントされる。各々の有効な音波経路での音波伝播時間をそれぞれ測定し、具体的には、各々の有効な音波経路での2つの音波センサ同士の間の向かい合っている方向の2つの音波伝播時間をそれぞれ測定し、この2つの音波伝播時間の平均値を当該有効な音波経路での音波伝播時間とする。
各々の有効な音波経路での音波伝播時間を以下のモデル式に代入して測定待ちのパイプまたはコンテナの横断面における軸方向の流速分布を再構築し、軸方向の流速分布u(x,y)を得る。
Figure 2022502675
ここで、lは、i番目の有効な音波経路であり、Lは、i番目の有効な音波経路での2つの音波センサ同士の間の直線距離であり、Δtは、i番目の有効な音波経路での音波伝播時間であり、Nは、有効な音波経路の数であり、cは、媒体の物理的特性の条件下の音波伝播速度である。
前のステップで得られた軸方向の流速分布u(x,y)を当該横断面に沿って積分して、パイプ内の流量を得る。
好ましくは、各々の音波センサは、順番に音を発し、また、すべての音波センサによって発射される音波の周波数は、いずれも同一である。任意の1つの音波センサによって音波が発射されるときに、他のすべての音波センサは、いずれも当該音波を記録する。一般的な測定対象は、約2秒で1回の完全な測定を完了でき、測定対象内の流れ場が一般的に比較的安定しているため、当該音を発する方法の測定精度に対する影響は十分に小さい。
好ましくは、2〜3個の音波センサごとに、同時に音を発することもでき、ここで、各々の音波センサは、互いに異なる音波の周波数を採用して発射する。各々の音波センサが順番に音を発することと比較すると、より短い時間でパイプ内の軸方向の流速分布を測定することができるが、互いに異なる周波数の音波信号を効果的に識別することはより難しい。
具体的に使用するときには、測定に求められる時間の長さに応じて制御することで、各々の音波センサが順番に音を発するか、2〜3個の音波センサごとに同時に音を発するようにすることができる。
測定にあたっては、測定対象のサイズおよび流体の媒体属性に応じて音波の周波数を合理的に選択することができることに注意すべきである。気相または多相の媒体の場合、音波の減衰率が周波数の2乗に比例するため、周波数が高いほど減衰が大きい。しかし、音波の周波数が高いほど、高精度の音波伝播時間データを得やすくなり、さらに、高精度の測定結果を得ることができるため、音波経路両端で十分な強さの音波信号を得ることができるように、可能な限り高い測定音波の周波数を採用して十分な測定精度を保証すべきである。
横断面内の半径方向の流速成分の影響を排除するために、2つの音波センサごとの向かい合っている方向の2つの音波伝播時間の平均値を、当該音波センサ同士の間の音波伝播時間とする。具体的な実施例において、横断面の周りに4、5、6、7、8個の音波センサを均一に配置することによって、それぞれ6、10、15、21、28個の有効な音伝播経路を形成することができる(図2参照)。
具体的に、前記軸方向の流速分布u(x,y)は、Taylor級数多項式を採用してフィッティング計算されるが、以下に3次精度の例を示す。
Figure 2022502675
ここで、o(x,y)は、3次の高次極小項であり、C、C、…、C10は、未定の多項式係数である。軸方向の流れ場の複雑さに基づいて構築する多項式の精度を決定し、その中の項の数を決定する。これに応じて、前記音波センサの数によって確定される有効な音波経路の数は、多項式係数の数以上である。
軸方向の流速分布が自由射流の特性に一致すると、前記軸方向の流速分布u(x,y)は、以下のガウス多項式を採用して近似フィッティング計算される。
Figure 2022502675
ここで、(x,y)は、自由射流の流速の極値点の座標であり、Uは、自由射流の流速の極値であり、uは、自由射流中心から離れた端部での流速値であり、σは、自由射流の膨張特性の幅である。対応する音波センサの数によって確定される有効な音波経路の数は、5個以上である。
具体的に、本発明で選択された横断面は、パイプ軸またはコンテナ軸と交差しており、選択された横断面は、パイプ軸またはコンテナ軸に直角であってもよいし、パイプ軸またはコンテナ軸にほぼ直角であってもよい。前記音波センサは、選択された横断面の周りの壁に沿って均一に配置されてもよいし、軸方向の流れ場の分布特徴に応じて選択された横断面の周りの壁に沿って不均一に配置されてもよく、軸方向の流速の変化が比較的に複雑な領域に対していっそう多い音波経路を得るようにする。
また、測定対象の横断面の形状、軸方向の流速場を再構築するモデル式の適切な変形などは、いずれも、本発明の效果に影響を及ぼさない。
本発明に係る測定方法は、円形、楕円形または長方形などのような、異なる形状の断面のパイプまたはコンテナに適用されるため、ここでは円形断面を有するパイプの軸方向の流速分布を例として、本発明の音波法による測定方法を検証する。
図2に示すように、当該円形断面のパイプの、パイプ内の安定した区間の1つの横断面を選択して、選択された横断面の周りの壁に沿って6個の音波センサを配置する。ここで、パイプ径は、1mであり、媒体は、空気であり、温度は、室温条件であり、本発明の方法を採用して、事前設定されたパイプ内の軸方向の流速分布に対して、音波法シミュレーション測定による再構築を実行する。
パイプ断面の幾何学的中心を原点として座標系を確立し、流れ場中心を幾何学的中心とする円形断面パイプ内の軸方向の流れ場のプリセットシミュレーションフィールド
Figure 2022502675
を、図3に示すように、構築する。
経路積分によって、表1に示すように、小数点以下6桁で、図2においてのすべての音波センサ同士の間の各々の音波経路の音波伝播時間が得られた。
表1:図3のシミュレーションパイプ内の軸方向の流れ場によって確定される各経路の伝播時間(単位:ms)
Figure 2022502675
前述したパイプ内の軸方向の流速分布と流量を音波法で測定する方法に基づき、Taylor級数フィッティング3次精度を採用して、前記軸方向の流速分布u(x,y)をフィッティング計算し、最終的に得られた、音波法シミュレーション測定再構築による流れ場の結果u(x,y)、およびその再構築の相対誤差REは、図4に示したとおりである。Taylor級数フィッティング4次精度を採用して、前記軸方向の流速分布u(x,y)をフィッティング計算し、最終的に得られた、音波法シミュレーション測定再構築による流れ場の結果u(x,y)、およびその再構築の相対誤差REは、図5に示したとおりである。
看取されるように、音波法シミュレーション測定の結果は、十分に正確で信頼性が高く、また、Taylor級数フィッティングの精度が向上すると、シミュレーション測定の精度も対応して向上することがわかる。
また、パイプ断面の幾何学的中心を原点として座標系を確立し、流れ場中心が幾何学的中心から偏向する円形断面パイプ内の軸方向の流れ場のプリセットシミュレーションフィールド
Figure 2022502675
を、図6に示すように、構築することもできる。
経路積分によって、表2に示すように、小数点以下6桁で、図2においてのすべての音波センサ同士の間の各々の音波経路の音波伝播時間が得られた。
表2:図6のシミュレーションパイプ内の軸方向の流れ場によって確定される各経路の伝播時間(単位:ms)
Figure 2022502675
前述したパイプ内の軸方向の流速分布と流量を音波法で測定する方法に基づき、Taylor級数フィッティング3次精度、および、Taylor級数フィッティング4次精度をそれぞれ採用して、前記軸方向の流速分布u(x,y)をフィッティング計算し、最終的に得られた、音波法シミュレーション測定再構築による流れ場の結果u(x,y)、およびその再構築の相対誤差REは、それぞれ図7および図8に示したとおりである。同様に看取されるように、音波法シミュレーション測定の結果は、十分に正確で信頼性が高く、また、Taylor級数フィッティングの精度が向上すると、シミュレーション測定の精度も対応して向上し、すなわち、音波法の測定精度がTaylor級数フィッティングの精度の影響を受けることがわかる。
音波法の測定精度に影響を与える要因をさらに説明するために、Taylor級数フィッティング3次精度、および、Taylor級数フィッティング4次精度の場合の、8組の音波センサの円形断面パイプ内の軸方向の流れ場に対するシミュレーション測定を実行した。得られた再構築の相対誤差REは、図9および図10に示したとおりである。
看取されるように、音波法の測定精度がTaylor級数フィッティングの精度の影響を受けるという結論に加えて、一般的で比較的に簡単な流れ場の場合、8組の音波センサのシミュレーション測定結果は、6組の音波センサのシミュレーション測定結果と比較して大幅に改善されなかったことがわかる。
上述したパイプ内の軸方向の流れ場の音波法シミュレーション測定結果の2つの代表例によれば、音波伝播経路の屈折原理に基づく音波法により軸方向の流れ場を測定する方法は、実現可能かつ信頼性があり、また、音波法の測定精度は、級数フィッティングの精度の影響を受け、音波センサの数とは関係がない。ただし、級数フィッティングの精度を向上させた場合、得られる基底関数係数の増加により、それに応じて音波センサの数を増やして、方程式を閉じて解く必要がある。
本発明は、上述したパイプ内の軸方向の流速分布と流量を音波法で測定する方法によるシステムをさらに提供し、当該システムは、パイプまたはコンテナの特定の横断面の周りの壁に沿って分布する複数の音波センサと、デジタル/アナログ変換カードと、アナログ/デジタル変換カードと、測定用のコンピューターと、を備え、前記測定用のコンピューターには、上述したパイプ内の軸方向の流速分布と流量を音波法で測定する方法によって構築された測定ソフトウェアが事前設定されている。
前記デジタル/アナログ変換カードは、それぞれ前記測定用のコンピューターおよび前記音波センサに接続され、前記測定用のコンピューターの測定ソフトウェアが発したオーディオデジタル信号を音波アナログ信号に変換し、前記音波センサを介して発信するために使用される。
前記アナログ/デジタル変換カードは、それぞれ前記測定用のコンピューターおよび前記音波センサに接続され、前記音波センサが収集したパイプ内の音波情報をデジタル信号に変換して前記測定用のコンピューターに入力するために使用される。
前記測定用のコンピューターは、前記測定ソフトウェアを利用してすべての音波センサが音を発するように制御し、各々の音波センサからその他の音波センサへの音波経路の音波伝播時間を測定し、各音波経路での音波伝播時間をパイプまたはコンテナの横断面における軸方向の流速分布の再構築式
Figure 2022502675
に代入することによって、軸方向の流速分布u(x,y)を得ることができ、ここで、lは、i番目の音波経路であり、Lは、i番目の音波経路に対応する音波センサペア間の距離であり、Δtは、当該音波経路での音波伝播時間、すなわち同一の音波センサペア間の向かい合っている方向の2つの音波伝播時間の平均であり、Nは、有効な音波経路の数であり、cは、音波伝播速度である。
前記測定用のコンピューターは、前記測定ソフトウェアを利用して、得られた軸方向の流速分布u(x,y)を当該測定横断面に沿って積分して、パイプ内の流量を得る。
発明に開示される方法は、気体、液体、および二相または多相が流れているパイプ内の軸方向の流速分布と流量に対する有効な測定に適用することができ、燃焼室、自由射流場、流動床、化学反応器などの内部流れ場や、開放空間での自由射流などの対象の測定に拡張することもできる。
最後に説明すべきことは、以上の実施例は、ただ、本発明の技術的解決策を説明するためのものであり、本発明に対する限定ではない。好ましい実施例を参照して本発明を詳細に説明したが、当業者は、依然として、本発明の具体的な実施形態に対して修正するか、または、一部の技術的特徴を同等に置き換えることができるが、本発明の技術的解決策の精神から逸脱することなく、それらのすべては、本発明の請求された技術的解決策の範囲によってカバーされなければならないと理解すべきである。

Claims (10)

  1. パイプ内の軸方向の流速分布と流量を音波法で測定する方法であって、
    パイプまたはコンテナの安定した区間または関心のある区間において一意の横断面を選択し、前記横断面の周りの壁に複数の音波センサを設置するステップであって、複数の音波センサ同士の間は、複数の有効な音波経路を形成し、ここで、各々のペアの音波センサ同士の間は、1つの有効な音波経路としてのみカウントされるステップと、
    各々の有効な音波経路での音波伝播時間をそれぞれ測定して、以下のモデル式に代入して測定待ちのパイプまたはコンテナの横断面の軸方向の流速分布u(x,y)を再構築するステップであって、
    Figure 2022502675
    ここで、lは、i番目の有効の音波経路であり、Lは、i番目の有効の音波経路での2つの音波センサ同士の間の直線距離であり、Δtは、i番目の有効の音波経路での音波伝播時間であり、Nは、有効の音波経路の数であり、cは、媒体の物理的特性の条件下の音波伝播速度であるステップと、
    前記軸方向の流速分布u(x,y)を当該横断面に沿って積分して、パイプ内の流量を得るステップと、を含む
    ことを特徴とする方法。
  2. 各々の有効な音波経路での2つの音波センサ同士の間の向かい合っている方向の2つの音波伝播時間をそれぞれ測定し、この2つの音波伝播時間の平均値を当該有効な音波経路での音波伝播時間とする
    ことを特徴とする請求項1に記載の方法。
  3. 各々の音波センサは、順番に音を発し、任意の1つの音波センサによって音波が発射されるときに、他のすべての音波センサは、いずれも当該音波を記録し、すべての音波センサによって発射される音波の周波数は、いずれも同一である
    ことを特徴とする請求項1または2に記載の方法。
  4. 2〜3個ごとの音波センサは、同時に音を発し、ここで、各々の音波センサによって発射される音波の周波数は、互いに異なっており、フィルタリングによってこれら互いに異なっている音波の周波数を効果的に識別する
    ことを特徴とする請求項1または2に記載の方法。
  5. 前記音波センサは、発受信一体型の音波センサである
    ことを特徴とする請求項1または2に記載の方法。
  6. 前記音波センサは、それぞれ別体である音波発射センサと音波受信センサとを組み合わせたものである
    ことを特徴とする請求項1または2に記載の方法。
  7. 前記軸方向の流速分布u(x,y)は、Taylor級数多項式を採用してフィッティング計算される
    ことを特徴とする請求項1または2に記載の方法。
  8. 前記軸方向の流速分布u(x,y)は、自由射流の特性に一致する場合、以下のガウス多項式を採用して近似フィッティング計算され、
    Figure 2022502675
    ここで、(x,y)は、自由射流の流速の極値点の座標であり、Uは、自由射流の流速の極値であり、uは、自由射流中心から離れた端部での流速値であり、σは、自由射流の膨張特性の幅である
    ことを特徴とする請求項1または2に記載の方法。
  9. 選択された横断面は、パイプ軸またはコンテナ軸と直角またはほぼ直角に交差する
    ことを特徴とする請求項1または2に記載の方法。
  10. 請求項1乃至9の中のいずれか1項に記載のパイプ内の軸方向の流速分布と流量を音波法で測定する方法のシステムであって、
    パイプまたはコンテナの特定の横断面の周りの壁に沿って分布する複数の音波センサと、デジタル/アナログ変換カードと、アナログ/デジタル変換カードと、測定用のコンピューターと、を備え、前記測定用のコンピューターには、前記したパイプ内の軸方向の流速分布と流量を音波法で測定する方法によって構築された測定ソフトウェアが事前設定されており、
    前記デジタル/アナログ変換カードは、それぞれ前記測定用のコンピューターおよび前記音波センサに接続され、前記測定用のコンピューターの測定ソフトウェアが発したオーディオデジタル信号を音波アナログ信号に変換し、前記音波センサを介して発信するために使用され、
    前記アナログ/デジタル変換カードは、それぞれ前記測定用のコンピューターおよび前記音波センサに接続され、前記音波センサが収集したパイプ内の音波情報をデジタル信号に変換して前記測定用のコンピューターに入力するために使用され、
    前記測定用のコンピューターは、前記測定ソフトウェアを利用してすべての音波センサが音を発するように制御し、各々の音波センサからその他の音波センサへの有効な音波経路での音波伝播時間を測定し、各々の有効な音波経路での音波伝播時間をパイプまたはコンテナの横断面における軸方向の流速分布の再構築式
    Figure 2022502675
    に代入することによって、軸方向の流速分布u(x,y)を得ることができ、ここで、lは、i番目の有効な音波経路であり、Lは、i番目の有効な音波経路での2つの音波センサ同士の間の直線距離であり、Δtは、i番目の有効な音波経路での音波伝播時間であり、Nは、有効な音波経路の数であり、cは、媒体の物理的特性の条件下の音波伝播速度であり、
    前記測定用のコンピューターは、前記測定ソフトウェアを利用して、得られた軸方向の流速分布u(x,y)を当該測定横断面に沿って積分して、パイプ内の流量を得る
    ことを特徴とするシステム。

JP2021531155A 2018-08-11 2018-09-29 パイプ内の軸方向の流速分布と流量を音波法で測定する方法、およびシステム Active JP7032842B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
CN201810910508.3 2018-08-11
CN201810910508 2018-08-11
PCT/CN2018/108837 WO2020034345A1 (zh) 2018-08-11 2018-09-29 声波法测量管内轴向流速分布、流量的方法及系统

Publications (2)

Publication Number Publication Date
JP2022502675A true JP2022502675A (ja) 2022-01-11
JP7032842B2 JP7032842B2 (ja) 2022-03-09

Family

ID=69525014

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2021531155A Active JP7032842B2 (ja) 2018-08-11 2018-09-29 パイプ内の軸方向の流速分布と流量を音波法で測定する方法、およびシステム

Country Status (5)

Country Link
US (1) US11454642B2 (ja)
EP (1) EP3835731A4 (ja)
JP (1) JP7032842B2 (ja)
CN (1) CN111051823B (ja)
WO (1) WO2020034345A1 (ja)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111051823B (zh) * 2018-08-11 2021-07-20 李言钦 声波法测量管内轴向流速分布、流量的方法及系统
CN113188746B (zh) * 2021-03-12 2023-03-14 同济大学 一种非接触式区域流体涡量测量方法

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2676321B2 (ja) * 1993-08-25 1997-11-12 株式会社昌民物産 超音波流量測定方法と装置
JP2935833B2 (ja) * 1996-05-27 1999-08-16 チャングミン シーオウ.エルティーディー. 多回線流量測定装置
JP3649028B2 (ja) * 1999-03-16 2005-05-18 富士電機システムズ株式会社 超音波流量計、及びその測定方法
CN102072966A (zh) * 2009-11-20 2011-05-25 湖北省电力试验研究院 一种声学法燃烧锅炉炉膛速度场的测量装置及测量方法
CN102174887A (zh) * 2011-01-05 2011-09-07 中国海洋石油总公司 一种超声波测量海底隔水管与钻柱间环空流量的装置
CN103336145A (zh) * 2013-06-18 2013-10-02 清华大学 一种针对管道内流体轴向流场的超声成像方法及设备
CN107607158A (zh) * 2017-09-18 2018-01-19 李言钦 声学法测量管内轴向流速分布及流量的系统
JP6582368B2 (ja) * 2014-07-23 2019-10-02 東京電力ホールディングス株式会社 流量計測装置および流量計測方法
US20200141776A1 (en) * 2018-08-11 2020-05-07 Yanqin Li Method and system of acoustic wave measurement of axial velocity distribution and flow rate
US20210223076A1 (en) * 2018-08-11 2021-07-22 Yanqin Li Method and system of acoustic wave measurement of axial velocity distribution and flow rate

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100401022C (zh) 2004-02-26 2008-07-09 富士电机系统株式会社 超声波流量计和超声波流量测量方法
CA2557432C (en) * 2004-02-26 2012-09-04 Fuji Electric Systems Co., Ltd. Ultrasonic flowmeter and ultrasonic flow rate measurement method
EP2282178B1 (en) 2009-07-23 2016-10-19 Elster NV/SA Ultrasonic flowmeter for measuring a flow characteristic of a fluid in a conduit
EP2383550A1 (en) 2010-04-28 2011-11-02 Miitors ApS Ultrasonic flow meter
CN102607653B (zh) 2012-04-06 2013-09-04 清华大学 一种基于超声波的管道流速成像和流量测量方法
CN202770855U (zh) 2012-08-15 2013-03-06 华北电力大学 一种基于声学的锅炉炉膛烟气流速监测系统
JP2014178202A (ja) 2013-03-14 2014-09-25 Osaka Gas Co Ltd 超音波流量計及び超音波流量測定方法
CN104395703A (zh) 2013-04-25 2015-03-04 株式会社又进 超声波流动测量系统
JP6212762B2 (ja) * 2013-10-18 2017-10-18 国立研究開発法人宇宙航空研究開発機構 流量計測方法及び装置
US9424674B2 (en) 2014-04-01 2016-08-23 Saudi Arabian Oil Company Tomographic imaging of multiphase flows
WO2019152041A1 (en) 2018-02-01 2019-08-08 Reliance Worldwide Corporation Sensor mount

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2676321B2 (ja) * 1993-08-25 1997-11-12 株式会社昌民物産 超音波流量測定方法と装置
JP2935833B2 (ja) * 1996-05-27 1999-08-16 チャングミン シーオウ.エルティーディー. 多回線流量測定装置
JP3649028B2 (ja) * 1999-03-16 2005-05-18 富士電機システムズ株式会社 超音波流量計、及びその測定方法
CN102072966A (zh) * 2009-11-20 2011-05-25 湖北省电力试验研究院 一种声学法燃烧锅炉炉膛速度场的测量装置及测量方法
CN102174887A (zh) * 2011-01-05 2011-09-07 中国海洋石油总公司 一种超声波测量海底隔水管与钻柱间环空流量的装置
CN103336145A (zh) * 2013-06-18 2013-10-02 清华大学 一种针对管道内流体轴向流场的超声成像方法及设备
JP6582368B2 (ja) * 2014-07-23 2019-10-02 東京電力ホールディングス株式会社 流量計測装置および流量計測方法
CN107607158A (zh) * 2017-09-18 2018-01-19 李言钦 声学法测量管内轴向流速分布及流量的系统
US20200141776A1 (en) * 2018-08-11 2020-05-07 Yanqin Li Method and system of acoustic wave measurement of axial velocity distribution and flow rate
US20210223076A1 (en) * 2018-08-11 2021-07-22 Yanqin Li Method and system of acoustic wave measurement of axial velocity distribution and flow rate

Also Published As

Publication number Publication date
US11454642B2 (en) 2022-09-27
EP3835731A4 (en) 2022-10-19
CN111051823A (zh) 2020-04-21
EP3835731A1 (en) 2021-06-16
WO2020034345A1 (zh) 2020-02-20
CN111051823B (zh) 2021-07-20
US20210223076A1 (en) 2021-07-22
JP7032842B2 (ja) 2022-03-09

Similar Documents

Publication Publication Date Title
CN106679741B (zh) 基于涡街流量计抗干扰信号的处理方法及系统
da Silva Hack et al. Measurement uncertainty: Literature review and research trends
CN103238065B (zh) 用超声波进行管线壁性能测试的系统和方法
JP7032842B2 (ja) パイプ内の軸方向の流速分布と流量を音波法で測定する方法、およびシステム
CN109375171B (zh) 一种基于正交匹配追踪算法的声源定位方法
CN108444653A (zh) 一种面阵型超声波气体泄漏点定位检测装置及方法
Fontana et al. Signal-theoretic characterization of waveguide mesh geometries for models of two-dimensional wave propagation in elastic media
Bridges et al. Direct evaluation of aeroacoustic theory in a jet
Hélie et al. One-dimensional acoustic models of horns and comparison with measurements
Jiang et al. A model-based hybrid ultrasonic gas flowmeter
Chen et al. Cobra probe measurements of mean velocities, Reynolds stresses and higher-order velocity correlations in pipe flow
CN109764956A (zh) 一种基于组合优化正则化方法的近场声全息方法
CN106707234B (zh) 一种联合时延差与角度测量的传感器网络目标定位方法
Liu et al. High temporal resolution pyrometry and velocimetry based on acoustic frequency division multiplexing
Zhang et al. Co-inversion of a scattering cavity and its internal sources: uniqueness, decoupling and imaging
Yu et al. A method for ultrasound thermal image distribution reconstruction
CN105115552A (zh) 基于侧音测相的超声波流量测量方法及装置
CN104840218B (zh) 一种呼吸速率的测量方法及电子设备
US10739174B2 (en) Method and system of acoustic wave measurement of axial velocity distribution and flow rate
Li et al. Acoustic reconstruction of the velocity field in a furnace using a characteristic flow model
RU2434206C1 (ru) Способ идентификации расхода теплоносителя по характеристической точке мощности и спектральной плотности эдс магнитоиндукционного преобразователя шарикового расходомера
CN106227962B (zh) 一种采用概率分布计算的水下管道周向模态识别方法
Zheng et al. Design of ultrasonic flowmeter for wet gas measurement
CN110261472A (zh) 小波基函数用于炉膛断面流场声波法测量重建的方法
CN115356372B (zh) 一种新型材料在飞行试验中的时变热响应测试方法及系统

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210222

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20210212

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210210

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20211106

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20211125

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20211129

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20220209

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220217

R150 Certificate of patent or registration of utility model

Ref document number: 7032842

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150