JP7175365B1 - Determination Method of Characteristic Time Reference Wave of Acoustic Signal of Ultrasonic Flowmeter - Google Patents
Determination Method of Characteristic Time Reference Wave of Acoustic Signal of Ultrasonic Flowmeter Download PDFInfo
- Publication number
- JP7175365B1 JP7175365B1 JP2021163641A JP2021163641A JP7175365B1 JP 7175365 B1 JP7175365 B1 JP 7175365B1 JP 2021163641 A JP2021163641 A JP 2021163641A JP 2021163641 A JP2021163641 A JP 2021163641A JP 7175365 B1 JP7175365 B1 JP 7175365B1
- Authority
- JP
- Japan
- Prior art keywords
- acoustic signal
- peak value
- characteristic
- time reference
- reference wave
- 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.)
- Active
Links
Images
Abstract
【課題】超音波流量計の音響信号の特性時間基準波の判断方法の提供。【解決手段】超音波流量計は、第１音響信号を送信し、第２音響信号を受信する第１音波送受信ユニットと、第２音響信号を送信し、第１音響信号を受信する第２音波送受信ユニットと、を備えている。方法は、（ａ）第１音響信号に対応する第１波形および第２音響信号に対応する第２波形を受信するステップと、（ｂ）第１波形および第２波形の複数のピーク値をサンプリングするステップと、（ｃ）これらのピーク値に基づき、サーチ範囲を設定するステップと、（ｄ）サーチ範囲内の第１ピーク値を特性ピーク値として設定するステップと、（ｅ）特性ピーク値後のゼロクロス点の前とゼロ点との間の第１時間およびゼロ点とゼロクロス点の後との間の第２時間を記録し、かつ第１時間および第２時間の平均時間を計算するステップと、（ｆ）平均時間に基づいて総飛行時間を計算するステップと、を含む。【選択図】図７A method for determining a characteristic time reference wave of an acoustic signal of an ultrasonic flowmeter. The ultrasonic flowmeter includes a first acoustic wave transmitting/receiving unit for transmitting a first acoustic signal and receiving a second acoustic signal, and a second acoustic wave transmitting and receiving unit for transmitting a second acoustic signal and receiving the first acoustic signal. a transmitting/receiving unit; The method includes the steps of (a) receiving a first waveform corresponding to the first acoustic signal and a second waveform corresponding to the second acoustic signal; and (b) sampling a plurality of peak values of the first waveform and the second waveform. (c) setting a search range based on these peak values; (d) setting a first peak value within the search range as a characteristic peak value; (e) after the characteristic peak value recording a first time before and between the zero-crossing points of and a second time between the zero-crossing points and after the zero-crossing points of and calculating an average time of the first and second times; , (f) calculating a total flight time based on the average time. [Selection drawing] Fig. 7
Description
本発明は、超音波流量計の音響信号の判断方法に関し、特に、超音波流量計の音響信号の特性時間基準波の判断方法に関する。 The present invention relates to a method for judging an acoustic signal of an ultrasonic flowmeter, and more particularly to a method for judging a characteristic time reference wave of the acoustic signal of an ultrasonic flowmeter.
流量計は、工業測定における重要な機器の一つであり、各産業の応用や科学研究において、密接な関連を有し、測定精度の要求もますます高くなっている。流量計は、各分野に広く応用され、半導体の製造過程のように、塗布設備、エッチング設備、洗浄設備、乾燥設備の製造過程でみな流量計の測定技術が使用されている。 A flowmeter is one of the important instruments in industrial measurement, and has a close relationship in the application and scientific research of each industry, and the demand for measurement accuracy is getting higher and higher. Flowmeters are widely applied in various fields, and flowmeter measurement technology is used in manufacturing processes such as coating equipment, etching equipment, cleaning equipment, and drying equipment, such as semiconductor manufacturing processes.
超音波技術は、従来はよく軍事や医療等の用途に使用されていたが、近年では多くの産業で応用されるようになっている。超音波流量計は、比較的最近出現した測定機器であるが、測定対象と非接触式であるため、流体に抵抗が生じず、センサの寿命が長く、汚染しない等の特性によって、近年はより多くの産業に注目され、使用されてきている。 Ultrasound technology has traditionally been used for military and medical applications, but in recent years it has found application in many industries. Ultrasonic flowmeters are a relatively recent measurement device, but since they do not come into contact with the object being measured, there is no resistance in the fluid, the sensor has a long service life, and it is free from contamination. It has been noticed and used in many industries.
超音波流量計の測定方式は、超音波の伝播時間差演算法を用いて管路内の流速を測定した後、流速により流量を換算するものである。超音波と流体とが同じ方向に移動する場合、流速が速いほど伝播の時間差は大きくなる。また、流体の状態が超音波の進む速度を変えるため、前後の検知器で得られた時間差によって流速および流量を概算することができる。 The measurement method of the ultrasonic flowmeter is to measure the flow velocity in the pipeline using the propagation time difference calculation method of ultrasonic waves, and then convert the flow rate from the flow velocity. When the ultrasonic wave and the fluid move in the same direction, the faster the flow velocity, the greater the time difference in propagation. Also, since the fluid conditions change the speed at which the ultrasonic wave travels, the time difference obtained by the front and rear detectors can be used to approximate the flow velocity and volume.
しかしながら、対応する前後の検知器が受信した音波（超音波）の波形が対応して一致しない場合、すなわち前の検知器と後の検知器とによって得られた音波特性時間基準波が一致しない場合、流速および流量の予測には極めて大きな誤差が生じ、精密測定の役目を果たすことができない。 However, if the waveforms of the sound waves (ultrasound) received by the corresponding front and rear detectors are correspondingly inconsistent, i.e., if the acoustic characteristic time reference waves obtained by the front and rear detectors are not in agreement , the prediction of flow velocity and flow rate has a very large error and cannot serve as a precise measurement.
このため、どのように超音波流量計の音響信号の特性時間基準波の判断方法を設計し、音波に対して前（予備）処理および音響信号の判断を行って、超音波流量計の精密測定を実現することは、本発明者が研究を行う大きな課題である。 Therefore, how to design the method of judging the characteristic time reference wave of the acoustic signal of the ultrasonic flowmeter, pre-process the sound wave and judge the acoustic signal, so that the precision measurement of the ultrasonic flowmeter is a major subject of research by the present inventors.
本発明は、従来技術の問題を解決できる、超音波流量計の音響信号の特性時間基準波の判断方法を提供することを目的とする。 SUMMARY OF THE INVENTION It is an object of the present invention to provide a method for determining a characteristic time reference wave of an acoustic signal of an ultrasonic flowmeter, which can solve the problems of the prior art.
前述の目的を実現するために、本発明が提案する超音波流量計の音響信号の特性時間基準波の判断方法では、前記超音波流量計は、第１音響信号を送信し、第２音響信号を受信する第１音波送受信ユニットと、第２音響信号を送信し、第１音響信号を受信する第２音波送受信ユニットと、を備えている。前記方法は、（ａ）第１音響信号に対応する第１波形および第２音響信号に対応する第２波形を受信するステップと、（ｂ）第１波形および第２波形の複数のピーク値をサンプリングするステップと、（ｃ）これらのピーク値に基づき、サーチ範囲を設定するステップと、（ｄ）サーチ範囲内の最初の前記ピーク値である第１ピーク値を特性ピーク値として設定するステップと、（ｅ）特性ピーク値後のゼロクロス点の前とゼロ点との間の第１時間およびゼロ点とゼロクロス点の後との間の第２時間を記録し、かつ第１時間および第２時間の平均時間を計算するステップと、（ｆ）平均時間に基づいて総飛行時間を計算するステップと、を含む。 To achieve the above object, the present invention proposes a method for judging a characteristic time reference wave of an acoustic signal of an ultrasonic flowmeter, wherein the ultrasonic flowmeter transmits a first acoustic signal and a second acoustic signal and a second acoustic wave transmitting/receiving unit for transmitting a second acoustic signal and receiving the first acoustic signal. The method comprises the steps of: (a) receiving a first waveform corresponding to a first acoustic signal and a second waveform corresponding to a second acoustic signal; and (b) determining a plurality of peak values of the first waveform and the second waveform. (c) setting a search range based on these peak values; and (d) setting a first peak value, which is the first peak value within the search range, as a characteristic peak value. , (e) recording a first time before and after the zero-crossing point after the characteristic peak value and a second time between the zero-crossing point and after the zero-crossing point, and and (f) calculating a total flight time based on the average time.
一実施例において、ステップ（ｂ）の前にさらに以下のステップを含む。（ｆ）標準差下限閾値と標準差上限閾値とを設定するステップ。ステップ（ｄ）はさらに以下のステップを含む。（ｄ１）第１ピーク値が標準差上限閾値より大きい場合、特性ピーク値として設定するステップと、（ｄ２）第１ピーク値が標準差下限閾値より小さい場合、特性ピーク値として排除するステップと、（ｄ３）続くピーク値が標準差上限閾値より大きいと判断される場合、特性ピーク値として設定するステップ。 In one embodiment, the following steps are further included before step (b). (f) setting a lower standard difference threshold and an upper standard difference threshold; Step (d) further includes the following steps. (d1) if the first peak value is larger than the standard difference upper limit threshold, set it as a characteristic peak value; (d2) if the first peak value is smaller than the standard difference lower limit threshold, exclude it as a characteristic peak value; (d3) If the subsequent peak value is determined to be greater than the standard difference upper threshold, setting it as a characteristic peak value.
一実施例において、ステップ（ｂ）の前にさらに以下のステップを含む。（ｆ）標準差下限閾値と標準差上限閾値とを設定するステップ。ステップ（ｄ）はさらに以下のステップを含む。（ｄ１）第１ピーク値が標準差上限閾値より大きい場合、特性ピーク値として設定するステップと、（ｄ２）第１ピーク値が標準差下限閾値より小さい場合、特性ピーク値として排除するステップと、（ｄ３）続くピーク値が標準差下限閾値以上であり、かつ標準差上限閾値以下であると判断される場合、後のピーク値が前のピーク値よりも大きいか否かを判断し、もし大きい場合は特性ピーク値として設定するステップ。 In one embodiment, the following steps are further included before step (b). (f) setting a lower standard difference threshold and an upper standard difference threshold; Step (d) further includes the following steps. (d1) if the first peak value is larger than the standard difference upper limit threshold, set it as a characteristic peak value; (d2) if the first peak value is smaller than the standard difference lower limit threshold, exclude it as a characteristic peak value; (d3) If it is determined that the subsequent peak value is greater than or equal to the standard difference lower threshold and is less than or equal to the standard difference upper threshold, determine whether the subsequent peak value is greater than the previous peak value; If so, the step of setting it as a characteristic peak value.
一実施例において、ステップ（ｄ３）は以下のステップを含む。（ｄ４）後のピーク値が前のピーク値より大きく、かつ後のピーク値が前のピーク値よりある変化量分大きい場合、特性ピーク値として設定するステップ。 In one embodiment, step (d3) includes the following steps. (d4) If the subsequent peak value is greater than the previous peak value, and the subsequent peak value is greater than the previous peak value by a certain amount of change, then setting it as a characteristic peak value.
一実施例において、変化量は任意の隣接する２つのピーク値の最大傾斜変化量の割合である。 In one embodiment, the variation is a percentage of the maximum slope variation of any two adjacent peak values.
一実施例において、特性時間基準波の判断方法はさらに以下のステップを含む。（ｅ１）第１音響信号の第１ピーク値と第２ピーク値を取得するステップと、（ｅ２）第２音響信号の第１ピーク値と第２ピーク値とを取得するステップと、（ｅ３）これら第１ピーク値とこれら第２ピーク値とを互いに比較し、第１音響信号と第２音響信号とが揃っているか否かを判断するステップと、（ｅ４）もし第１音響信号と第２音響信号とが揃っていない場合、互いに対応する第１ピーク値と第２ピーク値とを置き換えて、第１音響信号と第２音響信号とを揃うようにするステップ。 In one embodiment, the method for determining the characteristic time reference wave further includes the following steps. (e1) obtaining a first peak value and a second peak value of the first acoustic signal; (e2) obtaining a first peak value and a second peak value of the second acoustic signal; (e3) (e4) comparing the first peak values and the second peak values to determine whether the first acoustic signal and the second acoustic signal are aligned; If the acoustic signals are not aligned, replacing the first peak value and the second peak value corresponding to each other to align the first acoustic signal and the second acoustic signal.
一実施例において、ステップ（ｅ４）は以下を含む。第１音響信号が第２音響信号よりも前にある場合、第１音響信号の第２ピーク値を第１音響信号の第１ピーク値に置き換え、第２音響信号が第１音響信号よりも前にある場合、第２音響信号の第２ピーク値を第２音響信号の第１ピーク値に置き換える。 In one embodiment, step (e4) includes: If the first acoustic signal precedes the second acoustic signal, the second peak value of the first acoustic signal is replaced with the first peak value of the first acoustic signal, and the second acoustic signal precedes the first acoustic signal. , replace the second peak value of the second acoustic signal with the first peak value of the second acoustic signal.
一実施例において、標準差上限閾値は信号雑音比が１０に等しく、標準差下限閾値は信号雑音比が５に等しい。 In one embodiment, the standard difference upper threshold equals 10 for signal-to-noise ratio and the standard difference lower threshold equals 5 for signal-to-noise ratio.
一実施例において、第１音響信号の周波数および第２音響信号の周波数付近以外のピーク値は削除される。 In one embodiment, peak values other than around the frequency of the first acoustic signal and the frequency of the second acoustic signal are removed.
本発明が所定の目的を達成するために採用する技術、手段、効果をより一層明確に理解できるよう、以下に本発明の詳細な説明と図面を提供する。これにより本発明の目的、特徴、特性はより一層深く、かつ、具体的に理解できるようになるが、図面はあくまでも参考と説明のためであって、本発明を限定するものではないと理解されるべきである。 The detailed description and drawings of the present invention are provided below so that the techniques, means and effects adopted by the present invention to achieve its intended purpose can be more clearly understood. While the purpose, features and characteristics of the present invention may thereby be further and more specifically understood, it should be understood that the drawings are for reference and explanation only and are not limiting of the invention. should.
本発明によれば、流量の計算（予測）を実現することができる。また、超音波流量計の測定結果を向上させることができる。 According to the present invention, flow rate calculation (prediction) can be realized. Also, the measurement results of the ultrasonic flowmeter can be improved.
本発明の技術内容および詳細な説明について、以下、図面と合わせて説明する。 The technical content and detailed description of the present invention will be described below in conjunction with the drawings.
図１に示すように、図１は本発明の超音波流量計の操作概略図で縦断面を示す。本発明の超音波流量計１００は、第１音波送受信ユニット１１と第２音波送受信ユニット１２とを備えている。本実施例において、第１音波送受信ユニット１１と第２音波送受信ユニット１２とは、流通管路の外面の相対する位置に対となって設置されている。ただし、２つの送受信ユニットの対となった設置は、図１に示すものに限定されない。すなわち、他の実施例において、第１音波送受信ユニット１１と第２音波送受信ユニット１２とは、流通管路の外面の同一線上の位置に設置してもよい。第１音波送受信ユニット１１は第１音響信号Ｓ１を送信し、第２音響信号Ｓ２を受信するものであり、第２音波送受信ユニット１２は第２音響信号Ｓ２を送信し、第１音響信号Ｓ１を受信するものである。このうち、第１音響信号Ｓ１および第２音響信号Ｓ２は、超音波信号（ｕｌｔｒａｓｏｎｉｃ ｓｉｇｎａｌ）である。
As shown in FIG. 1, FIG. 1 is an operational schematic diagram of the ultrasonic flowmeter of the present invention and shows a longitudinal section. An
第１音波送受信ユニット１１は、流通管路内の流体の流動方向に対して斜め方向に超音波信号を伝送し、かつ第２音波送受信ユニット１２により受信するが、このとき、第１音波送受信ユニット１１および第２音波送受信ユニット１２の送受信を切り換え、第２音波送受信ユニット１２により流通管路内の流体の反流動方向に対して斜め方向に超音波信号を伝送し、かつ第１音波送受信ユニット１１により受信すると、超音波信号の流体内の伝達時間差から流量を測定できるようになる。詳細は後述する。
The first sound wave transmitting/receiving
図２に示すように、図２は本発明の超音波流量計の音響信号の波形概略図である。音響信号の波形は第１音波送受信ユニット１１が伝送し、第２音波送受信ユニット１２が受信する第１音響信号Ｓ１と、第２音波送受信ユニット１２が伝送し第１音波送受信ユニット１１が受信する第２音響信号Ｓ２とを含んでいる。２つの音響信号の間には時間差が存在するが、時間偏移の設定により２つの音響信号を、頭を揃えて比較することで、音響信号の特性時間基準波位置を正確に判断できるようになる。
As shown in FIG. 2, FIG. 2 is a waveform schematic diagram of the acoustic signal of the ultrasonic flowmeter of the present invention. The waveforms of the acoustic signals are a first acoustic signal S1 transmitted by the first sound wave transmitting/receiving
図３に示すように、図３は従来技術の超音波流量計の流量計算の概略図である。これらのうち、
ただし、Ｍは音波経路、Ｄは管路の内径、θは超音波の入射角、Ｃ０は音の媒体内での速度、Ｖは流体の平均速度（流速）、Ｔｕｐは上流の所要時間、Ｔｄｏｗｎは下流の所要時間、Δｔは２つの音波送受信ユニットの時間差である。したがって、式（１）～式（４）に基づき、流通管路内の流量を（予測）計算することができる。 where M is the sound wave path, D is the inner diameter of the pipe, θ is the incident angle of the ultrasonic wave, C0 is the velocity in the sound medium, V is the average velocity of the fluid (flow velocity), and T up is the required upstream time. , T down is the time required downstream, and Δ t is the time difference between the two acoustic wave transmitting/receiving units. Therefore, the flow rate in the flow pipe can be (predicted) calculated based on the equations (1) to (4).
本発明の主な技術的特徴は、音波送受信ユニットが受信した音響信号（例えば第２音波送受信ユニット１２が受信した第１音響信号Ｓ１と、第１音波送受信ユニット１１が受信した第２音響信号Ｓ２と）に対して前（予備）処理および音響信号の判断を行うことにある。
The main technical feature of the present invention is that the acoustic signals received by the sound wave transmitting/receiving unit (for example, the first sound signal S1 received by the second sound wave transmitting/receiving
図４に示すように、図４は音波送受信ユニットが受信した音響信号の波形概略図である。超音波流量計の音響信号の特性時間基準波を正確に判断（探知）する（すなわち流量計算における音波の１番目である初期波として用いる）ため、本発明の技術手段は以下のように説明される。 As shown in FIG. 4, FIG. 4 is a waveform diagram of the acoustic signal received by the acoustic wave transmitting/receiving unit. In order to accurately determine (detect) the characteristic time reference wave of the acoustic signal of the ultrasonic flowmeter (that is, use it as the initial wave, which is the first sound wave in the flow rate calculation), the technical means of the present invention are described as follows. be.
理想的な超音波信号は完全な包絡形態である。ただし測定の位置や測定の条件などの要素により、超音波信号はそのような完全な包絡形態を表せない。したがって、本発明が提案する音響信号の特性時間基準波の判断方法は、主に以下の幾つかの重要なステップを含むことができる。 An ideal ultrasound signal is a perfect envelope form. However, due to factors such as the position of measurement and the conditions of measurement, the ultrasonic signal cannot represent such a perfect envelope form. Therefore, the method for judging the characteristic time reference wave of the acoustic signal proposed by the present invention can mainly include the following key steps.
ステップ１：音波が受信されるタイミングを設定する。例を挙げて説明すると、第１音波送受信ユニット１１から発した第１音響信号Ｓ１が第２音波送受信ユニット１２に届く時間には遅延が生じる。したがって、第１音響信号Ｓ１が第２音波送受信ユニット１２に近づいたときに、例えば音波の速度、媒体の種類、・・・等の要素に基づき、音波が受信されるタイミングをどのように設定すれば、第２音波送受信ユニット１２の受信器をオン（有効）にして、所望の音響信号を受信し始めるようにするかは、例えば図４に示す時間ｔ１の箇所である。
Step 1: Set the timing at which sound waves are received. To explain with an example, there is a delay in the time when the first acoustic signal S1 emitted from the first sound wave transmitting/receiving
ステップ２：ノイズの標準差を計算する。音波を送信する周波数は既知のもの（例えば１ＭＨＺ）であるため、受信する音響信号の周波数に基づいてノイズか否かを判断し、ノイズ(大きさ)の標準差をさらに計算することができる。したがって、合理的なノイズの大きさの標準差に基づき、後に判断する閾値の大きさを動態的に設計することができる。 Step 2: Compute the standard difference of the noise. Since the frequency at which the sound waves are transmitted is known (e.g., 1 MH Z ), it is possible to determine whether it is noise or not based on the frequency of the received acoustic signal, and further calculate the standard difference of the noise (magnitude). . Therefore, based on a reasonable noise magnitude standard difference, the threshold magnitude to be determined later can be dynamically designed.
ステップ３：音波のピーク値を構築する。受信した音波のピーク値（バレー値を含む）をサンプリングすることにより、当該音波のピーク値、すなわち音波の波形特性を取得（構築）する。実際は、音波を送信する周波数は既知のものであるため、二次フィルタリング方式により正確な波形に属さない可能性のある正しくはないノイズまたはグリッチ（ｇｌｉｔｃｈｉｎｇ）を除去し、所望の音響信号として見なすことができ、すなわち所望の周波数上のピーク値だけを残すことができる。 Step 3: Construct the peak value of the sound wave. By sampling the peak value (including the valley value) of the received sound wave, the peak value of the sound wave, that is, the waveform characteristic of the sound wave is obtained (constructed). In practice, since the frequency at which the sound waves are transmitted is known, a secondary filtering scheme removes false noise or glitching that may not belong to the correct waveform and regards it as the desired acoustic signal. , i.e., leave only the peak values on the desired frequency.
ステップ４：サーチウィンドウを構築する。音響信号と音波のノイズおよびサンプリングするピーク値とを決定すると、予測された特性時間基準波位置に近づく前にサーチウィンドウ（ｓｅａｒｃｈ ｗｉｎｄｏｗ）を開くことができる。構築されたサーチウィンドウにより、サーチウィンドウ内のデータだけを計算でき、計算時間と計算するデータ量とを節約することができる。例を挙げて説明すると、通常、取得（サンプリング）するデータは４０００件余りであるが、４０００件余りのデータ全てを計算する場合、データ量の計算に負担が生じる。サーチウィンドウの範囲について言えば、特性時間基準波データをカバーできることを原則として設計を行う。例えば図４に示す時間ｔ２においてサーチウィンドウを開き、時間ｔ３においてサーチウィンドウを終了する。 Step 4: Construct a search window. Once the acoustic signal and sound wave noise and sampling peak values are determined, a search window can be opened prior to approaching the predicted characteristic time reference wave position. The constructed search window allows only the data within the search window to be calculated, saving computation time and the amount of data to compute. Taking an example, the data to be acquired (sampled) is usually more than 4000 cases, but when calculating all the data of more than 4000 cases, the calculation of the data amount becomes a burden. As for the search window range, the design is based on the principle that the characteristic time reference wave data can be covered. For example, the search window is opened at time t2 shown in FIG. 4, and the search window is terminated at time t3.
ステップ５：サーチウィンドウ内の第１ピーク値を取得する。サーチウィンドウとサンプリングされる複数のピーク値とに基づき、サーチウィンドウ内における第１ピーク値を取得することができる。例えば図４に示す時間ｔ２における当該ピーク値である。また、サーチウィンドウ内の第１ピーク値に基づいて音波の特性時間基準波であるか否かを判断し、もし特性時間基準波でなければ、次のピーク値を判断する。ただし、特性時間基準波であるか否かを判断する方式は以下のとおりである。 Step 5: Get the first peak value within the search window. Based on the search window and the sampled peak values, a first peak value within the search window can be obtained. For example, it is the peak value at time t2 shown in FIG. Also, based on the first peak value within the search window, it is determined whether or not it is the characteristic time reference wave of the sound wave, and if it is not the characteristic time reference wave, the next peak value is determined. However, the method for determining whether or not it is the characteristic time reference wave is as follows.
（１）信号雑音比（ＳＮＲ）が５より小さいかで判断する。もしピーク値がＳＮＲ＝５より小さい場合、当該ピーク値は特性時間基準波ではないと直接判断する。 (1) Determine whether the signal-to-noise ratio (SNR) is smaller than 5. If the peak value is less than SNR=5, it is directly determined that the peak value is not the characteristic time reference wave.
（２）信号雑音比（ＳＮＲ）が１０より大きいかで判断する。もしピーク値がＳＮＲ＝１０より大きい場合、当該ピーク値は特性時間基準波であると直接判断する。図４に示す時間ｔ４における当該ピーク値はＳＮＲ＝１０より大きいため、当該ピーク値は特性時間基準波であると判断される。 (2) Determine whether the signal-to-noise ratio (SNR) is greater than 10. If the peak value is greater than SNR=10, it is directly determined that the peak value is the characteristic time reference wave. Since the peak value at time t4 shown in FIG. 4 is greater than SNR=10, the peak value is determined to be the characteristic time reference wave.
（３）もし信号雑音比（ＳＮＲ）が５と１０との間にある場合、さらなる判断を行い、ピーク値が特性時間基準波であるか否かを決定する。詳細な説明は後述する。 (3) If the signal-to-noise ratio (SNR) is between 5 and 10, make a further decision to determine if the peak value is the characteristic time reference wave. A detailed description will be given later.
ただし、信号雑音比の値は前述の数値（ＳＮＲ＝５、ＳＮＲ＝１０）に限られるものではなく、あくまでもこれら２つの数値の大きさの関係によりこれらピーク値が特性時間基準波であるか否かを直接判断すること、または、２つの数値の間にある場合はより細かな判断を行うということを説明するものである。 However, the value of the signal-to-noise ratio is not limited to the above-mentioned values (SNR=5, SNR=10). or to make a finer judgment if it is between the two values.
ステップ６：ゼロクロス点時間を計算する。ピーク値の信号雑音比が１０より大きく特性時間基準波であると判断された後（例えば図４の時間ｔ４）、ゼロクロス点の前とゼロ点との間の時間差と、ゼロ点とゼロクロス点の後ｔｐとの間の時間差とにより算術平均を行って計算し、求められた時間差により音波の波形が有する偏移（ｏｆｆｓｅｔ）によって生じる誤差を解消することができる。 Step 6: Calculate the zero-crossing point time. After it is determined that the signal-to-noise ratio of the peak value is greater than 10 and that it is a characteristic time reference wave (for example, time t4 in FIG. 4), the time difference between before the zero crossing point and the zero crossing point, and between the zero point and the zero crossing point. An error caused by an offset of the waveform of the sound wave can be eliminated from the calculated time difference by arithmetically averaging the time difference between the post tp and the calculated time difference.
ステップ７：総飛行時間（ＴＯＦ、ｔｉｍｅ ｏｆ ｆｌｉｇｈｔ）を計算する。平均時間に基づいて総飛行時間を計算し、流量の計算（予測）を実現する。 Step 7: Calculate the total time of flight (TOF). Calculate the total flight time based on the average time to achieve flow rate calculation (prediction).
特性時間基準波を判断する方式のうち、（３）において、もし信号雑音比が５と１０との間にある場合は、比較的説明する意味がある。サンプリングしたピーク値が信号雑音比ＳＮＲ＝５または信号雑音比ＳＮＲ＝１０の近くにある場合、特性時間基準波であるか否かの誤判断が生じやすい。以下、図面と合わせて説明する。 Among the methods for determining the characteristic time reference wave, in (3), if the signal-to-noise ratio is between 5 and 10, it is relatively meaningful to explain. When the sampled peak value is near the signal-to-noise ratio SNR=5 or the signal-to-noise ratio SNR=10, it is easy to make an erroneous judgment as to whether it is the characteristic time reference wave. A description will be given below with reference to the drawings.
図５に示すように、図５は本発明の音響信号の特性時間基準波の判断に対して最適化を行った波形概略図である。以下の説明は、前述のステップ３（音波のピーク値を構築する）とステップ４（サーチウィンドウを構築する）とを前提としており、後述のピーク値はサーチウィンドウにおいて取得したピーク値である。図５に示すように、時間ｔ１のとき、可能性のある特性時間基準波を取得する（そのピーク値が信号雑音比５と１０との間にあるため）。したがって、サンプリングしたピーク値が時間ｔ１のときに取得したピーク値より大きいか否かを後に再判断する。もし大きくない場合、時間ｔ１のときに取得したピーク値は特性時間基準波ではないと判断する。また、新しいピーク値が特性時間基準波であるか否かを、正確な特性時間基準波が見つかるまで、後に判断し直す。 As shown in FIG. 5, FIG. 5 is a waveform schematic diagram optimized for determining the characteristic time reference wave of the acoustic signal of the present invention. The following description assumes steps 3 (constructing the peak values of the acoustic wave) and 4 (constructing the search window) above, and the peak values described below are the peak values obtained in the search window. As shown in FIG. 5, at time t1, we obtain a possible characteristic time reference wave (because its peak value is between the signal-to-noise ratio of 5 and 10). Therefore, it is later re-determined whether the sampled peak value is greater than the peak value obtained at time t1. If not, it is determined that the peak value obtained at time t1 is not the characteristic time reference wave. Further, whether or not the new peak value is the characteristic time reference wave is re-determined later until an accurate characteristic time reference wave is found.
図５に示すように、９０個のサンプリング点を経過した後、時間ｔ２のときに別のピーク値が見つかったと（仮定）する。このとき、第２ピーク値（Ｔｕｐ２）が第１ピーク値（Ｔｕｐ１）よりも大きいか否かを判断する。さらに、もし第２ピーク値（Ｔｕｐ２）が第１ピーク値（Ｔｕｐ１）より大きい場合、第２ピーク値（Ｔｕｐ２）が第１ピーク値（Ｔｕｐ１）よりも十分に大きいか否かを判断する。但し、判断の基準は隣接する２つのピーク値の最大傾斜を根拠として判断することができ、例えば、２つのピーク値の最大傾斜はＳｍａｘであり、もし第２ピーク値（Ｔｕｐ２）と第１ピーク値（Ｔｕｐ１）との差が最大傾斜の２０％（すなわち０．２＊Ｓｍａｘ）より大きい場合、このように限定するわけではないが、第２ピーク値（Ｔｕｐ２）は第１ピーク値（Ｔｕｐ１）より大きいと判断することができる。したがって、この第２ピーク値（Ｔｕｐ２）は新しい（可能性のある）特性時間基準波として認定（設定）される。そして、サンプリングピーク値に対してもさらに後に同様の判断を行い、ピーク値が信号雑音比５と１０との間にあるときの真の特性時間基準波を判断することができる。 As shown in FIG. 5, after 90 sampling points (assuming) another peak value is found at time t2. At this time, it is determined whether or not the second peak value (Tup2) is greater than the first peak value (Tup1). Further, if the second peak value (Tup2) is greater than the first peak value (Tup1), it is determined whether the second peak value (Tup2) is sufficiently greater than the first peak value (Tup1). However, the judgment criterion can be judged based on the maximum slope of two adjacent peak values. For example, the maximum slope of two peak values is Smax, and if the second peak value (Tup2) If the difference from the value (Tup1) is greater than 20% of the maximum slope (i.e. 0.2*Smax), the second peak value (Tup2) is replaced by the first peak value (Tup1), although not so limited. can be determined to be greater than Therefore, this second peak value (Tup2) is recognized (set) as a new (possible) characteristic time reference wave. Then, the sampling peak value is subjected to the same determination later, and the true characteristic time reference wave when the peak value is between the signal-to-noise ratio of 5 and 10 can be determined.
図６に示すように、図６は本発明の音響信号の特性時間基準波の判断に対して別の最適化を行った波形概略図である。前述のように、超音波流量計は、第１音響信号Ｓ１を送信し、第２音響信号Ｓ２を受信する第１音波送受信ユニット１１と、第２音響信号Ｓ２を送信し、第１音響信号Ｓ１を受信する第２音波送受信ユニット１２とを備えているため、音波の信号は似ている２組である。したがって、前述した判断ステップは、第１音波送受信ユニット１１にとって言えば第１ピーク値（Ｔｕｐ１）と第２ピーク値（Ｔｕｐ２）とで、第２音波送受信ユニット１２にとって言えば第１ピーク値（Ｔｄｎ１）と第２ピーク値（Ｔｄｎ２）とで、かつ同じ判断ステップが２回実行される。すなわち、第１音響信号Ｓ１に対して１回実行され、第２音響信号Ｓ２に対して１回実行される。
As shown in FIG. 6, FIG. 6 is a waveform diagram of another optimization for determining the characteristic time reference wave of the acoustic signal of the present invention. As mentioned above, the ultrasonic flow meter includes a first acoustic wave transmitting/receiving
より好ましくは、上流送受信ユニットが送信した音波と下流送受信ユニットが送信した音波とを比較する。まず、前述のステップに基づいて第１音波送受信ユニット１１の第１ピーク値（Ｔｕｐ１）および第２ピーク値（Ｔｕｐ２）と、第２音波送受信ユニット１２の第１ピーク値（Ｔｄｎ１）および第２ピーク値（Ｔｄｎ２）とを取得して、“大きさが近い”か否かを２つずつ比較する。すなわち、Ｔｕｐ１とＴｄｎ１とを比較し（（Ｔｄｎ１－Ｔｕｐ１）／Ｔｕｐ１により）、Ｔｕｐ１とＴｄｎ２とを比較し（（Ｔｄｎ２－Ｔｕｐ１）／Ｔｕｐ１により）、Ｔｕｐ２とＴｄｎ１とを比較し（（Ｔｄｎ１－Ｔｕｐ２）／Ｔｕｐ２により）、Ｔｕｐ２とＴｄｎ２とを比較し（（Ｔｄｎ２－Ｔｕｐ２）／Ｔｕｐ２により）、第１比較値Ｃｍｐ１、第２比較値Ｃｍｐ２、第３比較値Ｃｍｐ３、第４比較値Ｃｍｐ４をそれぞれ得る。
More preferably, the acoustic wave transmitted by the upstream transceiver unit and the acoustic wave transmitted by the downstream transceiver unit are compared. First, based on the above steps, the first peak value (Tup1) and the second peak value (Tup2) of the first sound wave transmitting/receiving
Ｔｕｐ１とＴｄｎ１とが近く、かつＴｕｐ２とＴｄｎ２とも近い場合は、第１音響信号Ｓ１と第２音響信号Ｓ２とは揃っていると見なすことができると判断する。逆に、Ｔｕｐ２とＴｄｎ１とが近い場合は、第１音響信号Ｓ１と第２音響信号Ｓ２とは揃っておらず、かつ第１音響信号Ｓ１は第２音響信号Ｓ２よりも前（先）にあることを表す。同様に、Ｔｕｐ１とＴｄｎ２とが近い場合は、第１音響信号Ｓ１と第２音響信号Ｓ２とは揃っておらず、かつ第２音響信号Ｓ２は第１音響信号Ｓ１よりも前（先）にあることを表す。 When Tup1 is close to Tdn1 and Tup2 is close to Tdn2, it is determined that the first acoustic signal S1 and the second acoustic signal S2 can be regarded as aligned. Conversely, when Tup2 and Tdn1 are close, the first acoustic signal S1 and the second acoustic signal S2 are not aligned, and the first acoustic signal S1 is ahead of the second acoustic signal S2. represents Similarly, when Tup1 and Tdn2 are close, the first acoustic signal S1 and the second acoustic signal S2 are not aligned, and the second acoustic signal S2 is ahead of the first acoustic signal S1. represents
２つの音波が揃っていないと一旦判断されると、ピーク値の順序を調整する。すなわち、第１音響信号Ｓ１が第２音響信号Ｓ２よりも前（先）にある時、Ｔｕｐ２はＴｕｐ１に置き換えられる。逆に、第２音響信号Ｓ２が第１音響信号Ｓ１よりも前（先）にある時、Ｔｄｎ２はＴｄｎ１に置き換えられる。このようにして、Ｔｕｐ１およびＴｄｎ１の計算と、Ｔｕｐ２およびＴｄｎ２の計算によって合致（マッチング）させることができ、すなわちＴｕｐ１はＴｄｎ１に対応して計算を行い、Ｔｕｐ２はＴｄｎ２に対応して計算を行う。信号を揃える調整（すなわち特性時間基準波を正確な位置に置き換えること）により、信号送信の際に消失が生じ（ノイズに入り）、特性時間基準波に偏移（ｓｈｉｆｔ）が生じる状況を除去（相殺）することができる。 Once it is determined that the two waves are not aligned, the order of the peak values is adjusted. That is, Tup2 is replaced with Tup1 when the first sound signal S1 precedes (begins) the second sound signal S2. Conversely, when the second acoustic signal S2 precedes the first acoustic signal S1, Tdn2 is replaced with Tdn1. In this way, the calculations of Tup1 and Tdn1 can be matched by the calculations of Tup2 and Tdn2, i.e. Tup1 performs the calculation corresponding to Tdn1 and Tup2 performs the calculation corresponding to Tdn2. Aligning the signal (i.e. replacing the characteristic time reference wave with the correct position) eliminates the situation where erasures (into the noise) occur in the signal transmission and the characteristic time reference wave shifts ( offset).
さらに好ましくは、微小な波形変動によって特性時間基準波の誤判断が起きないために、ピーク値を置換することによりすでに決定した特性時間基準波が動かないように確保することができる。例を挙げて説明すると、もしすでに特性時間基準波として決定されたＴｕｐ１、またはすでに特性時間基準波として決定されたＴｄｎ１の変動が大きくない場合、特性時間基準波をＴｕｐ１とＴｄｎ１とに維持し、これにより特性時間基準波が頻繁に変動することを低減することができる。ただし、特性時間基準波として決定されたＴｕｐ１または特性時間基準波として決定されたＴｄｎ１の変動が非常に大きい場合は、新しい特性時間基準波により特性時間基準波を置換する。 More preferably, in order to prevent erroneous determination of the characteristic time reference wave due to minute waveform fluctuations, it is possible to ensure that the already determined characteristic time reference wave does not move by replacing the peak value. For example, if Tup1 already determined as the characteristic time reference wave or Tdn1 already determined as the characteristic time reference wave does not fluctuate significantly, maintaining the characteristic time reference wave at Tup1 and Tdn1, This can reduce frequent fluctuations of the characteristic time reference wave. However, if Tup1 determined as the characteristic time reference wave or Tdn1 determined as the characteristic time reference wave fluctuates significantly, the new characteristic time reference wave replaces the characteristic time reference wave.
こうして、前述の基本的なステップ（ステップ１～ステップ７）により、特性時間基準波の位置を大体判断でき、かつ総飛行時間を算出することで、流量の計算（予測）を実現することができる。さらに、最適化された判断と調整とにより、見つけ出された特性時間基準波の位置をより正確なものとすることができ、それにより超音波流量計の測定結果を向上させ、本発明の技術的な効果を実現させることができる。
Thus, the position of the characteristic time reference wave can be roughly determined by the basic steps (
図７に示すように、図７は本発明の超音波流量計の音響信号の特性時間基準波の判断方法のフローチャートである。図１と合わせて参照すると、前記超音波流量計は、第１音響信号Ｓ１を送信し、第２音響信号Ｓ２を受信する第１音波送受信ユニット１１と、第２音響信号Ｓ２を送信し、第１音響信号Ｓ１を受信する第２音波送受信ユニット１２とを備えている。前記の方法は以下のステップを含む。まず、第１音響信号に対応する第１波形および第２音響信号に対応する第２波形を受信する（Ｓ１１）。その後、第１波形および第２波形の複数のピーク値をサンプリングする（Ｓ１２）。その後、これらのピーク値に基づき、サーチ範囲（ｓｅａｒｃｈ ｗｉｎｄｏｗ）を設定する（Ｓ１３）。その後、サーチ範囲内の第１ピーク値を特性ピーク値として設定する（Ｓ１４）。その後、特性ピーク値後のゼロクロス点の前とゼロ点との間の第１時間およびゼロ点とゼロクロス点の後との間の第２時間を記録し、かつ第１時間および第２時間の平均時間を計算する（Ｓ１５）。最後に、平均時間に基づいて総飛行時間（ｔｉｍｅ ｏｆ ｆｌｉｇｈｔ、ＴＯＦ）を計算する（Ｓ１６）。
As shown in FIG. 7, FIG. 7 is a flow chart of the method for determining the characteristic time reference wave of the acoustic signal of the ultrasonic flowmeter of the present invention. Referring in conjunction with FIG. 1, the ultrasonic flowmeter includes a first acoustic wave transmitting/receiving
ステップ（Ｓ１２）の前にさらに以下のステップを含む。標準差下限閾値と標準差上限閾値とを設定するステップ。ただし標準差下限閾値と標準差上限閾値とは信号雑音比（ｓｉｇｎａｌ－ｔｏ－ｎｏｉｓｅ ｒａｔｉｏ、ＳＮＲ）である。 The following steps are further included before step (S12). Setting a lower standard difference threshold and an upper standard difference threshold. However, the standard difference lower threshold and the standard difference upper threshold are the signal-to-noise ratio (SNR).
ステップ（Ｓ１４）はさらに以下のステップを含む。第１ピーク値が標準差上限閾値(例えばＳＮＲ＝１０)より大きい場合、前記特性ピーク値として設定する。第１ピーク値が標準差下限閾値(例えばＳＮＲ＝５)より小さい場合、前記特性ピーク値として排除する。または、続くピーク値が前記標準差上限閾値より大きいと判断される場合、前記特性ピーク値として設定する。 Step (S14) further includes the following steps. If the first peak value is greater than the standard difference upper limit threshold (for example, SNR=10), it is set as the characteristic peak value. If the first peak value is less than the lower standard difference threshold (eg, SNR=5), it is excluded as the characteristic peak value. Alternatively, if the subsequent peak value is determined to be greater than the standard difference upper threshold, it is set as the characteristic peak value.
ステップ（Ｓ１４）はさらに以下のステップを含む。第１ピーク値が前記標準差上限閾値より大きい場合、前記特性ピーク値として設定する。第１ピーク値が前記標準差下限閾値より小さい場合、前記特性ピーク値として排除する。または、続くピーク値が前記標準差下限閾値以上であり、かつ前記標準差上限閾値以下であると判断される場合、後のピーク値が前のピーク値よりも大きいか否かを判断し、もし大きい場合は前記特性ピーク値として設定する。ただし、後のピーク値が前のピーク値より大きく、かつ後のピーク値が前のピーク値よりある変化量分大きい場合、前記特性ピーク値として設定する。一実施例において、変化量は任意の隣接する２つのピーク値の最大傾斜変化量の割合である。 Step (S14) further includes the following steps. If the first peak value is greater than the standard difference upper threshold, it is set as the characteristic peak value. If the first peak value is smaller than the standard difference lower threshold, it is excluded as the characteristic peak value. Alternatively, if it is determined that the subsequent peak value is greater than or equal to the standard difference lower threshold and is less than or equal to the standard difference upper threshold, determine whether the subsequent peak value is greater than the previous peak value, and if If it is larger, it is set as the characteristic peak value. However, if the later peak value is greater than the previous peak value and the later peak value is greater than the previous peak value by a certain amount of change, it is set as the characteristic peak value. In one embodiment, the variation is a percentage of the maximum slope variation of any two adjacent peak values.
特性時間基準波の判断方法はさらに以下のステップを含む。第１音響信号の第１ピーク値と第２ピーク値とを取得する。第２音響信号の第１ピーク値と第２ピーク値とを取得する。前記第１ピーク値と前記第２ピーク値とを比較して、第１音響信号と第２音響信号とが揃っているか否かを判断する。および、もし第１音響信号と第２音響信号とが揃っていない場合、互いに対応する第１ピーク値と第２ピーク値とを置き換えて、第１音響信号と第２音響信号とを揃うようにする。ただし、第１音響信号が第２音響信号よりも前にある場合、第１音響信号の第２ピーク値を第１音響信号の第１ピーク値に置き換える。第２音響信号が第１音響信号よりも前にある場合、第２音響信号の第２ピーク値を第２音響信号の第１ピーク値に置き換える。 The method for determining the characteristic time reference wave further includes the following steps. A first peak value and a second peak value of the first acoustic signal are obtained. A first peak value and a second peak value of the second acoustic signal are obtained. By comparing the first peak value and the second peak value, it is determined whether or not the first acoustic signal and the second acoustic signal are aligned. And, if the first acoustic signal and the second acoustic signal are not aligned, the first peak value and the second peak value corresponding to each other are replaced so that the first acoustic signal and the second acoustic signal are aligned. do. However, if the first acoustic signal precedes the second acoustic signal, the second peak value of the first acoustic signal is replaced with the first peak value of the first acoustic signal. If the second acoustic signal precedes the first acoustic signal, the second peak value of the second acoustic signal is replaced with the first peak value of the second acoustic signal.
このようにして、本発明が提供する超音波流量計の音響信号の特性時間基準波の判断方法によれば、基本的なステップ（ステップ１～ステップ７）により特性時間基準波の位置を大体判断でき、かつ、総飛行時間を算出することで、それにより流量の計算（予測）を実現することができる。さらに、最適化された判断と調整とにより、見つけ出された特性時間基準波の位置をより正確なものとすることができ、それにより超音波流量計の測定結果を向上させ、本発明の技術的な効果を実現させることができる。
Thus, according to the method for judging the characteristic time-reference wave of the acoustic signal of the ultrasonic flowmeter provided by the present invention, the position of the characteristic time-reference wave is roughly judged through the basic steps (
以上は本発明の好ましい具体的な実施例についての詳細な説明および図であって、本発明の特徴はこれらに限られるものではなく、これにより本発明を限定するものではない。本発明の範囲は以下に記載の特許請求の範囲を基準とし、本発明の特許出願の範囲に合致する思想とそれに類似する変形例も全て本発明の範疇に含まれ、当業者が本発明の範囲において容易に行える変形や修正もみな以下の本特許請求の範囲に含まれると理解されるべきである。 The foregoing are detailed descriptions and illustrations of preferred and specific embodiments of the present invention, which are not intended to be limiting of the features of the present invention. The scope of the present invention is based on the claims set forth below, and all ideas conforming to the scope of the patent application of the present invention and similar modifications are also included in the scope of the present invention. It should be understood that all variations and modifications within the scope that can be readily made are included in the scope of the following claims.
１００ 超音波流量計
１１ 第１音波送受信ユニット
１２ 第２音波送受信ユニット
Ｓ１ 第１音響信号
Ｓ２ 第２音響信号
Ｓ１１～Ｓ１６ ステップ
100
Claims (9)
前記超音波流量計は、第１音響信号を送信し、第２音響信号を受信する第１の音波送受信ユニットと、前記第２音響信号を送信し、前記第１音響信号を受信する第２の音波送受信ユニットと、を備え、
前記方法は、
（ａ）前記第１音響信号に対応する第１波形および前記第２音響信号に対応する第２波形を受信するステップと、
（ｂ）前記第１波形および前記第２波形の複数のピーク値をサンプリングするステップと、
（ｃ）これらのピーク値に基づき、サーチ範囲を設定するステップと、
（ｄ）前記サーチ範囲内の最初の前記ピーク値である第１ピーク値を特性ピーク値として設定するステップと、
（ｅ）前記特性ピーク値後のゼロクロス点の前とゼロ点との間の第１時間およびゼロ点とゼロクロス点の後との間の第２時間を記録し、かつ前記第１時間および前記第２時間の平均時間を計算するステップと、
（ｆ）前記平均時間に基づいて総飛行時間を計算するステップと、を含む判断方法。 A method for determining a characteristic time reference wave of an acoustic signal of an ultrasonic flowmeter, comprising:
The ultrasonic flowmeter includes a first acoustic wave transmitting/receiving unit that transmits a first acoustic signal and receives a second acoustic signal, and a second acoustic wave transmitting/receiving unit that transmits the second acoustic signal and receives the first acoustic signal. a sound wave transmitting/receiving unit;
The method includes:
(a) receiving a first waveform corresponding to said first acoustic signal and a second waveform corresponding to said second acoustic signal;
(b) sampling a plurality of peak values of said first waveform and said second waveform;
(c) setting a search range based on these peak values;
(d) setting a first peak value, which is the first peak value within the search range, as a characteristic peak value;
(e) recording a first time before and after the zero-crossing point and a second time after the zero-crossing point after said characteristic peak value, and calculating an average time of 2 hours;
(f) calculating a total flight time based on said average time.
（ｆ）標準差下限閾値と標準差上限閾値とを設定するステップを含み、
ステップ（ｄ）は、さらに、
（ｄ１）第１ピーク値が前記標準差上限閾値より大きい場合、前記特性ピーク値として設定するステップと、
（ｄ２）第１ピーク値が前記標準差下限閾値より小さい場合、前記特性ピーク値として排除するステップと、
（ｄ３）続くピーク値が前記標準差上限閾値より大きいと判断される場合、前記特性ピーク値として設定するステップと、を含む、
請求項１に記載の超音波流量計の音響信号の特性時間基準波の判断方法。 prior to step (b), further
(f) setting a lower standard difference threshold and an upper standard difference threshold;
Step (d) further comprises:
(d1) if the first peak value is greater than the standard difference upper limit threshold, setting it as the characteristic peak value;
(d2) if the first peak value is smaller than the standard difference lower limit threshold, rejecting it as the characteristic peak value;
(d3) if the subsequent peak value is determined to be greater than the standard difference upper threshold, setting it as the characteristic peak value;
A method for determining a characteristic time reference wave of an acoustic signal of an ultrasonic flowmeter according to claim 1.
（ｆ）標準差下限閾値と標準差上限閾値とを設定するステップを含み、
ステップ（ｄ）は、さらに、
（ｄ１）第１ピーク値が前記標準差上限閾値より大きい場合、前記特性ピーク値として設定するステップと、
（ｄ２）第１ピーク値が前記標準差下限閾値より小さい場合、前記特性ピーク値として排除するステップと、
（ｄ３）続くピーク値が前記標準差下限閾値以上であり、かつ前記標準差上限閾値以下であると判断される場合、後のピーク値が前のピーク値よりも大きいか否かを判断し、もし大きい場合は前記特性ピーク値として設定するステップと、を含む、
請求項１に記載の超音波流量計の音響信号の特性時間基準波の判断方法。 prior to step (b), further
(f) setting a lower standard difference threshold and an upper standard difference threshold;
Step (d) further comprises:
(d1) if the first peak value is greater than the standard difference upper limit threshold, setting it as the characteristic peak value;
(d2) if the first peak value is smaller than the standard difference lower limit threshold, rejecting it as the characteristic peak value;
(d3) if it is determined that the subsequent peak value is greater than or equal to the standard difference lower threshold and is less than or equal to the standard difference upper threshold, determining whether the subsequent peak value is greater than the previous peak value; if larger, set as the characteristic peak value;
A method for determining a characteristic time reference wave of an acoustic signal of an ultrasonic flowmeter according to claim 1.
（ｄ４）後のピーク値が前のピーク値より大きく、かつ後のピーク値が前のピーク値よりある変化量分大きい場合、前記特性ピーク値として設定するステップを含む、
請求項３に記載の超音波流量計の音響信号の特性時間基準波の判断方法。 Step (d3) is
(d4) setting the characteristic peak value when the subsequent peak value is greater than the previous peak value and the subsequent peak value is greater than the previous peak value by a certain amount of change;
A method for determining a characteristic time reference wave of an acoustic signal of an ultrasonic flowmeter according to claim 3.
請求項４に記載の超音波流量計の音響信号の特性時間基準波の判断方法。 The amount of change is the ratio of the maximum slope change of any two adjacent peak values,
A method for determining a characteristic time reference wave of an acoustic signal of an ultrasonic flowmeter according to claim 4.
（ｅ１）前記第１音響信号の第１ピーク値と第２ピーク値とを取得するステップと、
（ｅ２）前記第２音響信号の第１ピーク値と第２ピーク値とを取得するステップと、
（ｅ３）これら第１ピーク値とこれら第２ピーク値とを互いに比較し、前記第１音響信号と前記第２音響信号とが揃っているか否かを判断するステップと、
（ｅ４）もし前記第１音響信号と前記第２音響信号とが揃っていない場合、互いに対応する前記第１ピーク値と前記第２ピーク値とを置き換えて、前記第１音響信号と前記第２音響信号とを揃うようにするステップと、を含む、
請求項４に記載の超音波流量計の音響信号の特性時間基準波の判断方法。 moreover,
(e1) obtaining a first peak value and a second peak value of the first acoustic signal;
(e2) obtaining a first peak value and a second peak value of the second acoustic signal;
(e3) comparing the first peak values and the second peak values to determine whether the first acoustic signal and the second acoustic signal are aligned;
(e4) If the first acoustic signal and the second acoustic signal are not aligned, the first peak value and the second peak value corresponding to each other are replaced to obtain the first acoustic signal and the second acoustic signal. aligning the acoustic signal;
A method for determining a characteristic time reference wave of an acoustic signal of an ultrasonic flowmeter according to claim 4.
前記第１音響信号が前記第２音響信号よりも前にある場合、前記第１音響信号の前記第２ピーク値を前記第１音響信号の前記第１ピーク値に置き換え、
前記第２音響信号が前記第１音響信号よりも前にある場合、前記第２音響信号の前記第２ピーク値を前記第２音響信号の前記第１ピーク値に置き換える、ことを含む、
請求項６に記載の超音波流量計の音響信号の特性時間基準波の判断方法。 Step (e4) is
replacing the second peak value of the first acoustic signal with the first peak value of the first acoustic signal if the first acoustic signal precedes the second acoustic signal;
replacing the second peak value of the second acoustic signal with the first peak value of the second acoustic signal if the second acoustic signal precedes the first acoustic signal;
A method for determining a characteristic time reference wave of an acoustic signal of an ultrasonic flowmeter according to claim 6.
請求項２または３に記載の超音波流量計の音響信号の特性時間基準波の判断方法。 the standard difference upper threshold equals a signal to noise ratio of 10 and the standard difference lower threshold equals a signal to noise ratio of 5;
4. A method for determining a characteristic time reference wave of an acoustic signal of an ultrasonic flowmeter according to claim 2 or 3.
請求項１に記載の超音波流量計の音響信号の特性時間基準波の判断方法。
In step (b), peak values other than near the frequency of the first acoustic signal and the frequency of the second acoustic signal are deleted;
A method for determining a characteristic time reference wave of an acoustic signal of an ultrasonic flowmeter according to claim 1.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2021163641A JP7175365B1 (en) | 2021-10-04 | 2021-10-04 | Determination Method of Characteristic Time Reference Wave of Acoustic Signal of Ultrasonic Flowmeter |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2021163641A JP7175365B1 (en) | 2021-10-04 | 2021-10-04 | Determination Method of Characteristic Time Reference Wave of Acoustic Signal of Ultrasonic Flowmeter |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP7175365B1 true JP7175365B1 (en) | 2022-11-18 |
| JP2023054661A JP2023054661A (en) | 2023-04-14 |
Family
ID=84101936
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2021163641A Active JP7175365B1 (en) | 2021-10-04 | 2021-10-04 | Determination Method of Characteristic Time Reference Wave of Acoustic Signal of Ultrasonic Flowmeter |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP7175365B1 (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5831175A (en) * | 1996-06-12 | 1998-11-03 | Welch Allyn, Inc. | Method and apparatus for correcting temperature variations in ultrasonic flowmeters |
| JP2005201836A (en) * | 2004-01-19 | 2005-07-28 | Tokyo Keiso Co Ltd | Propagation time difference type ultrasonic flowmeter |
| JP2014224685A (en) * | 2013-05-15 | 2014-12-04 | パナソニック株式会社 | Flow rate measuring device |
| JP2020034367A (en) * | 2018-08-29 | 2020-03-05 | アズビル株式会社 | Zero cross detection circuit, ultrasonic flowmeter, and zero cross detection method |
-
2021
- 2021-10-04 JP JP2021163641A patent/JP7175365B1/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5831175A (en) * | 1996-06-12 | 1998-11-03 | Welch Allyn, Inc. | Method and apparatus for correcting temperature variations in ultrasonic flowmeters |
| JP2005201836A (en) * | 2004-01-19 | 2005-07-28 | Tokyo Keiso Co Ltd | Propagation time difference type ultrasonic flowmeter |
| JP2014224685A (en) * | 2013-05-15 | 2014-12-04 | パナソニック株式会社 | Flow rate measuring device |
| JP2020034367A (en) * | 2018-08-29 | 2020-03-05 | アズビル株式会社 | Zero cross detection circuit, ultrasonic flowmeter, and zero cross detection method |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2023054661A (en) | 2023-04-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| WO2014205884A1 (en) | Time difference type ultrasonic flowmeter measurement method based on zero point analysis | |
| US7096135B2 (en) | Method and system for calculating the transit time of an ultrasonic pulse | |
| US6647805B2 (en) | Transit-time difference type ultrasonic flowmeter | |
| CN103542901B (en) | A kind of flowmeter | |
| RU2182335C2 (en) | Procedure measuring time of propagation of sound signal in fluid medium and method measuring velocity of flow of fluid medium | |
| US9140594B2 (en) | Ultrasonic, flow measuring device | |
| CN107860430B (en) | Time difference measurement method of ultrasonic gas flowmeter based on time difference method | |
| US20150160052A1 (en) | Ultrasonic, Flow Measuring Device | |
| CN107655533B (en) | A kind of Ultrasonic Wave Flowmeter signal processing method and system based on backward energy integral | |
| CN108548578B (en) | Ultrasonic echo signal characteristic peak identification method based on self-adaptive threshold | |
| CN110346004B (en) | Flow measurement data fusion method of dual-channel ultrasonic time difference method | |
| JP7175365B1 (en) | Determination Method of Characteristic Time Reference Wave of Acoustic Signal of Ultrasonic Flowmeter | |
| CN112304376B (en) | Ultrasonic flowmeter flow measuring method based on data fusion | |
| JP2007187506A (en) | Ultrasonic flowmeter | |
| US6595070B1 (en) | Acoustic flow meters | |
| TWI790714B (en) | Method of determining characteristic time reference wave of acoustic signal of ultrasonic flowmeter | |
| US20230243682A1 (en) | Ultrasonic flow measurement | |
| CN115727907A (en) | Method for judging characteristic time reference wave of sound wave signal of ultrasonic flowmeter | |
| CN114674384A (en) | Ultrasonic flowmeter error wave prevention detection method, device and equipment and flowmeter | |
| CN111473827B (en) | V-shaped sound channel zero drift elimination method | |
| CN110455360B (en) | Ultrasonic water meter | |
| Ma | Ultrasonic Gas Flowmeter Based on Dynamic Reference Waveform Cross-Correlation Algorithm | |
| Li et al. | A novel differential time-of-flight algorithm for high-precision ultrasonic gas flow measurement | |
| RU2768295C1 (en) | Method of measuring fluid velocity in a pipeline | |
| CN116338240B (en) | Ultrasonic liquid flow velocity measurement method and device based on parabolic fitting |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20211004 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20211029 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20220614 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20220830 |
|
| 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: 20221101 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20221108 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 7175365 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |