KR102183549B1 - Real-Time 3D Flow Meter with Ultrasonic Multiple Sensors - Google Patents

Real-Time 3D Flow Meter with Ultrasonic Multiple Sensors Download PDF

Info

Publication number
KR102183549B1
KR102183549B1 KR1020190112206A KR20190112206A KR102183549B1 KR 102183549 B1 KR102183549 B1 KR 102183549B1 KR 1020190112206 A KR1020190112206 A KR 1020190112206A KR 20190112206 A KR20190112206 A KR 20190112206A KR 102183549 B1 KR102183549 B1 KR 102183549B1
Authority
KR
South Korea
Prior art keywords
ultrasonic
flow
pipe
sensor
fluid
Prior art date
Application number
KR1020190112206A
Other languages
Korean (ko)
Inventor
장국진
신종훈
Original Assignee
(주)플로트론
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 (주)플로트론 filed Critical (주)플로트론
Priority to KR1020190112206A priority Critical patent/KR102183549B1/en
Application granted granted Critical
Publication of KR102183549B1 publication Critical patent/KR102183549B1/en

Links

Images

Classifications

    • 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
    • G01F15/00Details of, or accessories for, apparatus of groups G01F1/00 - G01F13/00 insofar as such details or appliances are not adapted to particular types of such apparatus
    • G01F15/06Indicating or recording devices

Landscapes

  • Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Measuring Volume Flow (AREA)

Abstract

The present invention relates to a real-time three-dimensional (3D) flowmeter using an ultrasonic multi-sensor and, more specifically, to a real-time three-dimensional (3D) flowmeter using an ultrasonic multi-sensor, capable of providing a timer obtaining ultrasonic speed and waveform data through a high-speed ADC with a PC board by enabling three-dimensional simulation measurement on a flow of fluids through the calculation of a route with an ultrasonic multi-sensor, and a high-precision ultrasonic flowmeter deriving an equation based on the collected data by carrying out a fluid test. According to the present invention, a flange for connection to a flow path of a pipe is formed at both ends formed into a hollow shape such that fluids can flow, a saddle is placed in a hole formed on the pipe to be inclined to an axial direction of the pipe such that ultrasonic sensors are inserted to be installed opposite each other, and the ultrasonic sensors are installed to be controlled in electric connection with a controller. The ultrasonic sensors are combined with the saddle to enable the two ultrasonic sensors to have an ultrasonic stroke by having an irregular grid shape such that the sensors face each other, and the ultrasonic sensors provide a structure for three-dimensionally measuring the flow velocity distribution of fluids in the hollow hole. Meanwhile, an error is led to be ±0.3% by using a multi-factor in a section of 500m^3/h-4000m^3/h through the grid stroke of the ultrasonic sensors.

Description

초음파 다중센서를 활용한 실시간 쓰리디 유량계{Real-Time 3D Flow Meter with Ultrasonic Multiple Sensors}Real-Time 3D Flow Meter with Ultrasonic Multiple Sensors}

본 발명은 초음파 다중센서를 활용한 실시간 쓰리디(3D) 유량계에 관한 것으로, 더욱 상세하게는 초음파 다중센서로 경로를 계산하여 유체 흐름을 3D로 모사 계량계측이 가능해 PC보드로 고속의 ADC를 사용해 초음파 속도와 파형 데이터를 획득할 수 있는 타이머, 유체 테스트를 진행하여 수집된 데이터를 기준으로 방정식을 도출해 고정밀도 초음파 유량계를 제공하는 초음파 다중센서를 활용한 실시간 쓰리디(3D) 유량계에 관한 것이다.The present invention relates to a real-time 3D flow meter using multiple ultrasonic sensors, and more specifically, it is possible to measure and measure the fluid flow in 3D by calculating a path with a multiple ultrasonic sensor. It relates to a real-time 3D flow meter using a timer that can acquire ultrasonic velocity and waveform data, and an ultrasonic multi-sensor that provides a high-precision ultrasonic flow meter by deriving an equation based on the data collected through a fluid test.

수도를 각 가정과 산업현장에 공급하기 위하여 상하수도를 설치하는 경우를 비롯하여, 도시가스의 공급, 송유관을 통한 원유의 이송, 농업용수를 위한 개수로의 설치 및 철강, 화학, 석유화학분야의 냉각수의 순환 등 오늘날 가스나 액체 등 유체를 이용하고 있는 산업분야는 수없이 많으며, 이러한 산업분야에서 유체의 유량을 적절히 조절하여 관리하는 것은 매우 중요한 문제이다.Including the installation of water and sewage to supply water to each home and industrial site, supply of city gas, transport of crude oil through pipelines, installation of open channels for agricultural water, and circulation of cooling water in the fields of steel, chemicals and petrochemicals There are numerous industrial fields that use fluids such as gases and liquids today, and it is very important to properly control and manage the fluid flow rate in these industrial fields.

이렇게 유체의 유량관리가 중요한 산업현장에서는 일반적으로 유량계를 사용하고 있으며, 유량계들 중 특히 초음파 유량계가 널리 사용되고 있다. 초음파 유량계는 유로에 설치된 초음파 센서에서 초음파를 수발신하여 유속을 측정하고, 이 유속을 이용하여 유량을 측정하도록 구성되어 있다.In industrial sites where fluid flow rate management is important, a flow meter is generally used, and among flow meters, an ultrasonic flow meter is particularly widely used. The ultrasonic flow meter is configured to measure a flow rate by receiving and transmitting ultrasonic waves from an ultrasonic sensor installed in a flow path, and measuring a flow rate using this flow rate.

종래의 초음파 유량계와 유속측정원리를 이용한 초음파 유량계(10)는 내측을 통해 유체가 흐를 수 있도록 중공형으로 이루어진 측정관을 구비한다. 이 측정관의 양단에는 배관(11) 등 유로와의 연결을 위한 플랜지(14)가 형성되어 있다. 또한, 배관(11)에는 새들(20)이 마련된다. 이 새들(20)은 배관(11)의 축방향에 대하여 경사지게 배치된다. 각 새들(20)에는 초음파 센서(30)가 삽입되어 설치되며, 두 개의 초음파 센서(30)는 상호 마주하게된다. 초음파 센서(30)는 컨트롤러(40)와 전기적으로 연결되어 컨트롤러(40)에 의하여 제어된다.The conventional ultrasonic flow meter and the ultrasonic flow meter 10 using the flow rate measurement principle has a measuring tube made of a hollow type so that the fluid flows through the inside. Flanges 14 for connection with a flow path such as a pipe 11 are formed at both ends of the measuring tube. In addition, a saddle 20 is provided in the pipe 11. The saddle 20 is disposed to be inclined with respect to the axial direction of the pipe 11. An ultrasonic sensor 30 is inserted into each saddle 20 and installed, and the two ultrasonic sensors 30 face each other. The ultrasonic sensor 30 is electrically connected to the controller 40 and is controlled by the controller 40.

상기한 구성의 초음파 유량계에서 유체의 유속측정은 일반적으로 전파시간차 방법에 의하여 얻어진다. 즉, 한쌍의 초음파 센서가 유체의 진행방향(측정관의 축방향)에 대하여 일정 각도(θ)를 형성하며 각각 유로의 한 지점과, 유체의 유동방향상 상기 한 지점의 하류측에 위치하는 하측지점에 서로 대면하도록 설치한다.In the ultrasonic flow meter of the above configuration, the flow velocity measurement of the fluid is generally obtained by the propagation time difference method. That is, a pair of ultrasonic sensors form a certain angle (θ) with respect to the moving direction of the fluid (the axial direction of the measuring tube), and each one point of the flow path and the lower side located downstream of the one point in the flow direction of the fluid. Install them so that they face each other at the points.

유체가 움직이지 않는 조건에서 일측 초음파 센서에서 발사된 초음파가 타측 초음파 센서에 도달한 전달시간과 역방향으로 발사된 초음파의 전달시간은 동일하여야 한다. 그러나 측정관을 지나는 유체의 속도로 인하여, 한 지점에서 발사된 초음파가 하측지점까지 도달하는 시간과 하측지점에서 발사된 초음파가 한 지점까지 도달하는 시간은 차이가 발생한다. 즉, 초음파가 유체의 진행방향에 대해 순방향으로 발사되는 경우의 전파시간은 초음파가 유체의 진행방향에 대해 역방향으로 발사된 경우의 전파시간에 비해서 짧으므로, 시간의 차이가 발생한다. 이러한 시간차를 이용하여 유체의 속도를 구할 수 있으며, 유체의 속도에 유로의 단면적을 곱해 유량을 산출하게 된다.Under the condition that the fluid does not move, the delivery time of the ultrasonic waves emitted from one ultrasonic sensor to the other ultrasonic sensor must be the same as the delivery time of the ultrasonic waves emitted in the reverse direction. However, due to the velocity of the fluid passing through the measuring tube, there is a difference between the time for the ultrasonic wave emitted from one point to reach the lower point and the time for the ultrasonic wave emitted from the lower point to reach one point. That is, the propagation time when the ultrasonic wave is emitted in the forward direction with respect to the moving direction of the fluid is shorter than the propagation time in the case where the ultrasonic wave is emitted in the reverse direction with respect to the moving direction of the fluid, a time difference occurs. Using this time difference, the velocity of the fluid can be obtained, and the flow rate is calculated by multiplying the velocity of the fluid by the cross-sectional area of the flow path.

상기한 구성의 초음파 유량계는 두 가지 방식으로 설치된다. 즉, 측정관을 미리 제작하여 수도관 등에 플랜지 접합시켜 사용하는 방식과, 별도로 제작된 측정관을 사용하는 것이 아니라 수도관 등 기존의 관을 측정관으로 사용하도록 수도관을 천공하여 초음파 센서를 설치하는 직접 시공방식이 있다.The ultrasonic flow meter of the above configuration is installed in two ways. In other words, a method of using a measuring pipe by flange bonding to a water pipe, etc., and a direct construction where an ultrasonic sensor is installed by drilling a water pipe to use an existing pipe such as a water pipe as a measuring pipe instead of using a separately manufactured measuring pipe. There is a way.

직접 시공방식에서는 유체의 흐름을 중단시키지 않은 상태에서(수도관의 경우 단수) 기존의 관로를 천공해야 하In the direct construction method, the existing pipe must be drilled without stopping the flow of fluid (single number in the case of water pipes).

므로, 천공된 단면에 부식방지용 도장을 하는 것이 불가능하였다. 이에 기존 관로의 천공된 단면이 부식될 뿐만 아니라, 부식면이 점차 확장되면 초음파 센서의 전면을 덮게 되어 초음파의 수발신이 정상적으로 이루어지지 않는 문제점이 있었다.Therefore, it was impossible to paint the perforated section for corrosion protection. Accordingly, not only the perforated cross section of the existing pipe is corroded, but when the corroded surface is gradually expanded, the entire surface of the ultrasonic sensor is covered, so that the reception and transmission of ultrasonic waves is not performed normally.

또한, 천공된 단면의 부식이 장기간 지속될 경우, 관의 두께가 얇아지게 되어 강도가 약화됨으로써 누수 및 관의 파손을 일으킬 위험도 있다.In addition, if the corrosion of the perforated cross-section continues for a long period of time, the thickness of the pipe becomes thin and the strength is weakened, and there is a risk of leakage and damage to the pipe.

[문헌 1] 특허등록번호 제1218780호(2012. 12. 28. 등록)[Document 1] Patent Registration No. 1218780 (registered on December 28, 2012) [문헌 2] 특허등록번호 제0560364호(2006. 03. 07. 등록)[Document 2] Patent Registration No. 0560364 (registered on March 07, 2006) [문헌 3] 특허등록번호 제0784415호(2007. 12. 04. 등록)[Document 3] Patent Registration No. 0784415 (registered on December 04, 2007) [문헌 4] 특허등록번호 제1897601호(2018. 09. 05. 등록)[Document 4] Patent registration number 1897601 (2018. 09. 05. registration)

따라서 본 발명은 이러한 종래의 결점을 해소하기 위하여 안출된 것으로, 본 발명의 해결과제는, 초음파 다중센서로 관로 내경에 공급되는 유체의 흐름을 3D로 모사 유량계측을 수행하도록 하는 것을 목적으로 한다.Therefore, the present invention was conceived to solve such a conventional shortcoming, and an object of the present invention is to simulate the flow of a fluid supplied to the inner diameter of a pipe by a multi-sonic ultrasonic sensor in 3D.

본 발명의 다른 해결과제는, PC 보드로 고속의 ADC를 사용해 초음파의 속도와 파형 데이터를 획득할 수 있는 타이머를 얻을 수 있도록 하는 것을 목적으로 한다.Another object of the present invention is to obtain a timer capable of acquiring ultrasonic speed and waveform data using a high-speed ADC with a PC board.

본 발명의 또 다른 해결과제는, 유체 테스트를 진행하여 수집된 데이터를 기준으로 방정식을 도출, 알고리즘을 이식하는 것을 목적으로 한다.Another object of the present invention is to derive an equation based on the collected data through a fluid test and implant an algorithm.

본 발명은 내측을 통해 유체가 흐를 수 있도록 중공형으로 이루어진 양단에는 배관의 유로와 연결을 위한 플랜지가 형성되어 있고, 상기 배관에 형성한 구멍에 새들이 배관의 축방향에 대하여 경사지게 배치되어 초음파 센서가 삽입되어 상호 마주하게 설치되며, 초음파 센서는 컨트롤러와 전기적으로 연결되어 제어되도록 설치하는 것을 포함하고,
상기 새들에 초음파 센서를 결합해 두 개의 초음파 센서가 상호 마주보도록 비규칙적인 격자 형태가 되어 초음파 행정을 갖도록 설치하며, 상기 초음파 센서는 중공의 유체 유속 분포를 입체적으로 측정하기 위한 구조를 제공하되;
상기 초음파 센서의 격자 행정을 통하여 500㎥/h∼4000㎥/h 구간에서 멀티 팩터를 사용해 오차가 ±0.3% 사이가 되도록 하는 것을 포함하는 것이다.
In the present invention, flanges for connection with the flow path of the pipe are formed at both ends of the pipe to allow fluid to flow through the inside, and the saddle is arranged inclined with respect to the axial direction of the pipe in the hole formed in the pipe, so that the ultrasonic sensor Is inserted and installed facing each other, and the ultrasonic sensor is electrically connected to the controller and installed to be controlled,
An ultrasonic sensor is coupled to the saddle to form an irregular grid so that the two ultrasonic sensors face each other to have an ultrasonic stroke, and the ultrasonic sensor provides a structure for three-dimensionally measuring a flow velocity distribution of a hollow fluid;
It includes making the error between ±0.3% by using a multi factor in the interval of 500㎥/h to 4000㎥/h through the grating stroke of the ultrasonic sensor.

삭제delete

삭제delete

본 발명은 초음파 다중센서로 관로 내경에 공급되는 유체의 흐름을 다양한 각도에서 측정하고 측정되는 유속분포를 3D로 모사 유량계측을 수행하는 효과를 제공하는 것이다.The present invention provides an effect of measuring the flow of fluid supplied to the inner diameter of a pipe at various angles with an ultrasonic multi-sensor and performing a flow rate measurement that simulates the measured flow velocity distribution in 3D.

본 발명은 PC 보드로 고속의 ADC를 사용해 초음파의 속도와 파형 데이터를 획득하고, 이러한 파형 데이터를 통해 적정한 보정계수를 찾아 보정할 수 있는 효과를 제공하는 것이다.The present invention provides an effect of acquiring speed and waveform data of ultrasonic waves using a high-speed ADC with a PC board, and finding and correcting an appropriate correction factor through the waveform data.

본 발명은 유체 테스트를 진행하여 수집된 데이터를 기준으로 방정식을 도출, 오차없는 데이터를 제공할 수 있도록 알고리즘을 이식하는 효과를 제공하는 것이다.The present invention provides an effect of implanting an algorithm so as to derive an equation based on data collected by performing a fluid test and provide data without errors.

도 1 은 본 발명의 바람직한 실시예를 나타낸 설치상태 정면도
도 2 는 본 발명의 설치상태를 나타낸 평면도
도 3 은 본 발명의 배관을 절단해 센서를 설치한 상태의 측면 단면도
도 4 는 본 발명의 배관에 대한 새들을 설치하기 이전에 대한 정면도
도 5 는 본 발명의 다른 실시예를 나타낸 정면도
도 6 은 본 발명의 다른 실시예를 나타낸 평면도
도 7 은 본 발명의 초음파 유량계에 대한 외형을 나타낸 사시도
도 8 은 본 발명의 배관 내부에 대한 초음파 진행방향을 나타낸 측면도
도 9 는 본 발명의 배관 내부에 대한 초음파 진행방향을 나타낸 다른 상태의 측면도
도 10 은 본 발명의 배관 내부에 대한 초음파 진행방향을 나타낸 또 다른 상태의 측면도
도 11 은 본 발명의 살라미 유속 함수 적분 영역을 나타낸 면적도
도 12 는 본 발명의 축대칭 유동 유속분포 보정계수를 나타낸 그림
도 13 내지 도 16 는 본 발명의 단일피크 유동 유속분포 보정계수를 나타낸 그림
도 16 과 도 17 은 본 발명의 이중피크 유동 유속분포 보정계수를 나타낸 그림
1 is a front view of an installation state showing a preferred embodiment of the present invention
Figure 2 is a plan view showing the installation state of the present invention
Figure 3 is a side cross-sectional view of a state in which a sensor is installed by cutting the pipe of the present invention
Figure 4 is a front view before installing the saddle for the piping of the present invention
5 is a front view showing another embodiment of the present invention
6 is a plan view showing another embodiment of the present invention
7 is a perspective view showing the appearance of the ultrasonic flow meter of the present invention
8 is a side view showing the direction of ultrasonic waves in the pipe of the present invention
9 is a side view of another state showing the direction of ultrasonic waves inside the pipe of the present invention
10 is a side view of another state showing the direction of ultrasonic waves in the pipe of the present invention
11 is an area diagram showing an integral region of the salami flow rate function of the present invention
12 is a diagram showing the correction factor of the axisymmetric flow velocity distribution of the present invention
13 to 16 are diagrams showing the correction factor for a single peak flow velocity distribution according to the present invention
16 and 17 are diagrams showing the correction factor of the double peak flow velocity distribution of the present invention

본 명세서 및 청구범위에 사용된 용어나 단어는 통상적이거나 사전적인 의미로 한정해서 해석되어서는 안 되며, 출원인이 그 자신의 고안을 가장 최선의 방법으로 설명하기 위해 용어의 개념을 적절하게 정의할 수 있다는 원칙에 입각하여 본 발명의 기술적 사상에 부합하는 의미와 개념으로 기재한 것으로 해석되어야만 한다.Terms or words used in this specification and claims should not be construed as limited to their usual or dictionary meanings, and the applicant may appropriately define the concept of terms in order to explain his or her design in the best way. It should be interpreted as having a meaning and concept consistent with the technical idea of the present invention based on the principle that there is.

따라서 본 명세서에 기재된 실시 예와 도면에 도시된 구성은 본 발명의 가장 바람직한 실시 예에 불과할 뿐이고 본 고안의 기술적 사상을 모두 대변하는 것은 아니므로, 본 출원시점에 이들을 대체할 수 있는 다양한 균등물과 변형 예들이 있을 수 있음을 이해하여야 한다.Therefore, the embodiments described in the present specification and the configurations shown in the drawings are only the most preferred embodiments of the present invention, and do not represent all the technical ideas of the present invention, and various equivalents that can replace them at the time of application It should be understood that there may be variations.

이하 도면을 참조하여 설명하기에 앞서, 본 발명의 요지를 드러내기 위해서 필요하지 않은 사항 즉 통상의 지식을 가진 당업자가 자명하게 부가할 수 있는 공지 구성에 대해서는 도시하지 않거나, 구체적으로 기술하지 않았음을 밝혀둔다.Prior to the following description with reference to the drawings, matters that are not necessary to reveal the gist of the present invention, that is, known configurations that can be obviously added by those skilled in the art are not shown or specifically described. Reveal.

이하에서 본 발명의 바람직한 실시예를 첨부시킨 도면에 따라 보다 상세하게 설명하기로 한다.Hereinafter, preferred embodiments of the present invention will be described in more detail with reference to the accompanying drawings.

도 1은 본 발명의 바람직한 실시예를 나타낸 설치상태 정면도이고, 도 2는 본 발명의 설치상태를 나타낸 평면도, 도 3은 본 발명의 배관을 절단해 센서를 설치한 상태의 측면 단면도, 도 4는 본 발명의 배관에 대한 새들을 설치하기 이전에 대한 정면도를 나타낸 것이다.1 is a front view of an installation state showing a preferred embodiment of the present invention, Figure 2 is a plan view showing the installation state of the present invention, Figure 3 is a side sectional view of a state in which a sensor is installed by cutting the pipe of the present invention, and Figure 4 is It shows a front view before installing the saddle for the piping of the present invention.

초음파 유량계와 유속측정원리를 이용한 초음파 유량계(10)는 내측을 통해 유체가 흐를 수 있도록 중공(12)을 갖는 배관(11)의 양쪽으로 유로와의 연결을 위한 플랜지(14)가 형성되어 있다. The ultrasonic flow meter 10 using the ultrasonic flow meter and the flow velocity measurement principle has flanges 14 for connection with a flow path on both sides of a pipe 11 having a hollow 12 so that a fluid can flow through the inside.

상기 배관(11)에는 초음파 센서(30)를 설치할 개수와 위치를 설정하여 8개의 구멍(13)을 형성하고, 이 구멍(13)에 새들(20)을 용접 고정하되;Eight holes 13 are formed in the pipe 11 by setting the number and location of the ultrasonic sensors 30 to be installed, and the saddle 20 is welded to the holes 13;

상기 새들(20)은 배관(11)의 축방향에 대하여 경사지게 배치되도록 하고, 상기 각 새들(20)에는 초음파 센서(30)가 삽입되어 설치되며, 두 개의 초음파 센서(30)는 상호 마주보도록 비규칙적인 격자 형태가 되도록 설치하고, 상기 초음파 센서(30)는 미도시된 컨트롤러(미도시)와 전기적으로 연결되어 컨트롤러에 연결된다.The saddle 20 is disposed to be inclined with respect to the axial direction of the pipe 11, and an ultrasonic sensor 30 is inserted and installed in each saddle 20, and the two ultrasonic sensors 30 are non-facing so as to face each other. It is installed so as to form a regular grid, and the ultrasonic sensor 30 is electrically connected to a controller (not shown) and is connected to the controller.

도 5는 본 발명의 다른 실시예를 나타낸 정면도이고, 도 6은 본 발명의 다른 실시예를 나타낸 평면도, 도 7은 본 발명의 초음파 유량계에 대한 외형을 나타낸 사시도를 나타낸 것이다.Figure 5 is a front view showing another embodiment of the present invention, Figure 6 is a plan view showing another embodiment of the present invention, Figure 7 is a perspective view showing the appearance of the ultrasonic flow meter of the present invention.

초음파 유량계(10)는 중공(12)을 갖는 의 양쪽으로 플랜지(14)가 설치되어 있고, 상기 배관(11)에 구멍(13)을 6개씩 12개 형성한 후 새들(20)을 설치해 상기 새들(20)에 초음파 센서(30)를 결합해 두 개의 초음파 센서(30)가 상호 마주보도록 비규칙적인 격자 형태가 되도록 설치하는 것이다.The ultrasonic flow meter 10 has flanges 14 installed on both sides of the pipe 11 having a hollow 12, and after forming 12 holes 13 in each of the pipe 11, a saddle 20 is installed. The ultrasonic sensor 30 is coupled to 20 and installed so that the two ultrasonic sensors 30 face each other in an irregular grid shape.

도 8은 본 발명의 배관 내부에 대한 초음파 진행방향을 나타낸 측면도이고, 도 9는 본 발명의 배관 내부에 대한 초음파 진행방향을 나타낸 다른 상태의 측면도, 도 10은 본 발명의 배관 내부에 대한 초음파 진행방향을 나타낸 또 다른 상태의 측면도를 나타낸 것이다.FIG. 8 is a side view showing the direction of ultrasonic propagation inside the pipe of the present invention, FIG. 9 is a side view showing the direction of ultrasonic propagation inside the pipe of the present invention, and FIG. 10 is a side view showing the ultrasonic propagation direction inside the pipe of the present invention. It shows a side view of another state showing the direction.

상측에 컨트롤러(40)가 설치되는 초음파 유량계(10)는 중공(12)을 갖는 의 양쪽으로 플랜지(14)가 설치되어 있고, 상기 배관(11)에 구멍(13)을 형성한 후 새들(20)을 설치해 상기 새들(20)에 초음파 센서(30)를 결합해 두 개의 초음파 센서(30)가 상호 마주보도록 비규칙적인 격자 형태가 되어 초음파 행정(50)을 갖도록 설치하는 것이다.The ultrasonic flow meter 10 on which the controller 40 is installed on the upper side has flanges 14 installed on both sides of it having a hollow 12, and after forming a hole 13 in the pipe 11, the saddle 20 ) Is installed to couple the ultrasonic sensor 30 to the saddle 20 to form an irregular grid so that the two ultrasonic sensors 30 face each other to have an ultrasonic stroke 50.

상기 초음파 센서(30)는 중공(12)의 유체 유속 분포를 입체적으로 측정하기 위한 구조를 제공하는 것이다. The ultrasonic sensor 30 provides a structure for three-dimensionally measuring the fluid flow velocity distribution in the hollow 12.

관내 유속은 위치에 따라 차이가 있으며, 단일 센서로 계측시 정밀도가 낮아 유체의 레이노즐 수 계산이 필요하다.The flow velocity in the pipe varies depending on the location, and the precision is low when measuring with a single sensor, so it is necessary to calculate the number of ray nozzles of the fluid.

초음파는 인간이 들을 수 있는 가청 최대 한계 범위를 넘어서는 주파수대를 의미하며, 압전효과가 발견된 이후 초음파 소나(항해 목적에 음파를 이용하는 것의 총칭)를 통해 발전하게 되었다. Ultrasound refers to a frequency band that exceeds the maximum audible range that humans can hear, and has developed through ultrasonic sonar (a generic term for using sound waves for navigation purposes) after the discovery of the piezoelectric effect.

초음파 기술은 유량측정에도 활용되어 1970년대에 산업현장에서 유량측정에 응용되기 시작하여 현재 0.1% 고정밀도의 초음파 유량계가 상용화 되어 있다.Ultrasonic technology was also used for flow measurement and began to be applied to flow measurement in industrial sites in the 1970s, and now, 0.1% high-precision ultrasonic flow meters are commercially available.

초음파 유량계는 유체가 흐르는 관내에 음파를 보내면 유량이 흐르는 같은 방향으로는 음속에 유속이 더해지고, 반대 방향으로는 음속에 유속이 빠진 속도로 전해진다. 두 방향의 음속을 측정하여 유속을 구하고, 배관경을 통해 체적유량을 측정할 수 있다. When sound waves are sent into a pipe through which a fluid flows, the ultrasonic flow meter is transmitted at a speed in which the speed of sound is added to the speed of sound in the same direction in which the flow rate flows, and the speed is subtracted from the speed of sound in the opposite direction. The flow velocity can be obtained by measuring the sound velocity in two directions, and the volumetric flow rate can be measured through the pipe diameter.

측정방법은 크게 전반속도차법과 도플러법 2가지로 분류할 수 있고, 전반속도차법은 전반시간차법, 주파수차법, 위상차법 3가지로 나누어 진다.Measurement methods can be largely classified into two types, the propagation speed difference method and the Doppler method, and the propagation speed difference method is divided into three types: the total time difference method, the frequency difference method, and the phase difference method.

초음파 센서의 부착위치에 따라 건식(Clamp type, Handy type)과 습식(Wet type)으로 나누어 지고, 센서 수에 따라 단일빔, 다중빔으로 나누어 진다.It is divided into dry (Clamp type, Handy type) and wet type (Wet type) according to the attachment position of the ultrasonic sensor, and is divided into single beam and multi-beam according to the number of sensors.

전반 시간차법은 유체 종류에 따른 유체 통과 음속을 알고 있어야 하며 그러기 위해서는 관내를 유체로 채우고 유량이 정지된 상태에서 유체를 통과하는 음속을 측정하여 프로세서는 알고 있어야 한다. 또한 온도에 따른 음속의 변화 특성을 알고 있어야 정밀도 높게 유속 측정이 가능하다.The propagation time difference method must know the speed of sound passing through the fluid according to the type of fluid, and to do so, the processor must know by measuring the speed of sound passing through the fluid when the tube is filled with fluid and the flow rate is stopped. In addition, it is possible to measure the flow rate with high precision only when you know the characteristics of the change of the sound speed according to temperature.

주파수차법은 유체 통과 음속이 수식에서 빠지게 되어 온도에 따른 음속의 변화로 야기되는 수식적 정밀성이 낮아지는 전반시간차법을 보완할 수 있는 방식으로 고정밀도 유량계에서 사용하기 적합하다.The frequency difference method is suitable for use in a high-precision flow meter as a method that can compensate for the total time difference method in which the mathematical precision caused by the change in the sound speed according to temperature due to the sound velocity passing through the fluid is omitted from the equation.

위상차법은 전반시잔차법과 원리는 동일하며 전반시간차법 보다 구현상 미세신호의 위상 검출이 가능하다면 전반시간차법 보다 정밀한 유속을 구할 수 있으나 유체통과 음속의 상수는 동일하게 수식에 들어가기 때문에 온도에 따른 음속 변화 특성은 알고 있어야 한다.The phase difference method has the same principle as the total time residual method, and if it is possible to detect the phase of a fine signal in an implementation rather than the total time difference method, a more precise flow velocity can be obtained than the total time difference method. It is important to be aware of the characteristics of the speed change.

도플러법은 엠블런스의 싸이렌 소리가 다가올 때와 지나갈 때의 소리가 틀리게 들리는 현상과도 같은 도플러효과를 활용한 것으로 유체의 흐름과 반대방향 및 같은 방향의 주파수를 측정하여 주파수의 변동으로 유량을 계측할 수 있다. 도플러법은 유체중에 특정 크기 이상의 부유물이 혼합되어 있어야 하므로 하수나 슬러지등의 유량계측에 적합하다.The Doppler method utilizes the Doppler effect, which is a phenomenon in which the siren sound of the ambulance approaches and passes incorrectly.It measures the frequency in the opposite direction and the same direction as the flow of the fluid, and measures the flow rate by fluctuation of the frequency. I can. The Doppler method is suitable for measuring the flow rate of sewage or sludge, as the fluid must contain suspended solids of a certain size or more.

초음파 유량계(10)는 주파수차법, 습식, 다중빔 구조로 제작하여 고정밀도 초음파 가스용 유량계 구현이 가능하도록 하였다.The ultrasonic flow meter 10 was manufactured in a frequency difference method, wet type, and a multi-beam structure to enable a high-precision ultrasonic gas flow meter.

v : 유속v: flow rate

C : 유체 통과 음속C: speed of sound through fluid

Figure 112019093202213-pat00001
(식1)
Figure 112019093202213-pat00001
(Equation 1)

Figure 112019093202213-pat00002
(식2)
Figure 112019093202213-pat00002
(Equation 2)

Figure 112019093202213-pat00003
(식3)
Figure 112019093202213-pat00003
(Equation 3)

Figure 112019093202213-pat00004
(식4)
Figure 112019093202213-pat00004
(Equation 4)

주파수차법은 수식에 음속 상수가 들어가지 않는다.The frequency difference method does not include a sound velocity constant in the equation.

다중빔을 사용하면 배관(11)내의 높이에 따른 유속을 측정할 수 있고 초음파 센서(30) 사이의 중간값은 스플라인 삽입(spline interpolation) 기법을 활용하여 높은 정도의 평균유속을 구할 수 있고, 초음파 센서(30) 일부 고장 발생 시에도 유속 산출이 가능하다는 장점이 있다.If multiple beams are used, the flow velocity according to the height of the pipe 11 can be measured, and the intermediate value between the ultrasonic sensors 30 can be obtained by using a spline interpolation technique to obtain a high degree of average flow velocity. There is an advantage in that the flow rate can be calculated even when some failures of the sensor 30 occur.

주파수 계산은 각각의 초음파 센서(30)로부터 데이터가 들어오면 제로-크로스(zero-crossing) 위치의 시간을 구하여 t(k-1), t(k), t(k+1), t(k+2) 라고 한다면 주파수 f 는 아래 방적식을 통해 계산한다.Frequency calculation is performed by calculating the time at the zero-crossing position when data is received from each ultrasonic sensor 30, t(k-1), t(k), t(k+1), t(k). If it is +2), then the frequency f is calculated using the equation below.

Figure 112019093202213-pat00005
(식5)
Figure 112019093202213-pat00005
(Equation 5)

Figure 112019093202213-pat00006
: A or B, 센서채널 frequency (
Figure 112019093202213-pat00007
or
Figure 112019093202213-pat00008
).
Figure 112019093202213-pat00006
: A or B, sensor channel frequency (
Figure 112019093202213-pat00007
or
Figure 112019093202213-pat00008
).

Figure 112019093202213-pat00009
: A or B, 센서채널 zero-crossing time (
Figure 112019093202213-pat00010
or
Figure 112019093202213-pat00011
).
Figure 112019093202213-pat00009
: A or B, sensor channel zero-crossing time (
Figure 112019093202213-pat00010
or
Figure 112019093202213-pat00011
).

제로-크로스(Zero-crossing) 시간(time) 계산은, 그래프의 마크(o) 포인트가 ADC 결과이고 마크간의 간격을 샘플링 타임(sampling time)으로 본다면 제로-크로스 타임(zero-crossing time)은 0을 기준으로 ADC결과가 음/양이 변화되는 구간으로 두 마크의 값과 샘플링 타임을 통해 제로-크로스를 계산한다. 계산은 정밀성을 높이기 위해 표준의 ‘큐브 스플라인 삽입(cube spline interpolation)’ 방식을 사용하거나 샘플링 타임을 감안하여 성능 검증 후 간단하게 ‘선형 보간법(linear interpolation)’방식을 사용한다.In the calculation of zero-crossing time, if the mark (o) point of the graph is the ADC result and the interval between marks is considered as the sampling time, the zero-crossing time is 0. Based on the ADC result, the zero-cross is calculated through the values of the two marks and the sampling time as the section in which the negative/positive changes. For the calculation, the standard'cube spline interpolation' method is used to increase precision, or the'linear interpolation' method is simply used after performance verification in consideration of the sampling time.

진폭, 위상 계산은, 진폭과 위상은 아래의 방정식을 통해 구한다.To calculate amplitude and phase, obtain the amplitude and phase through the equation below.

Figure 112019093202213-pat00012
(식6)
Figure 112019093202213-pat00012
(Equation 6)

Figure 112019093202213-pat00013
(식7)
Figure 112019093202213-pat00013
(Equation 7)

Figure 112019093202213-pat00014
(식8)
Figure 112019093202213-pat00014
(Equation 8)

Figure 112019093202213-pat00015
(식9)
Figure 112019093202213-pat00015
(Equation 9)

x : A or B, 채널x: A or B, channel

a : amplitudea: amplitude

D : data value (ADC result)D: data value (ADC result)

Φ : phaseΦ: phase

t' : t' = t-tx(k-1) or 0t': t'= t-tx(k-1) or 0

초음파 회선을 배열하는 방법으로는, 초음파 센서(30)를 설치하는 위치 계산 방법으로 쳬비셰프 적분에 적합하도록 초음파 회선을 배열하는 방법이 적당하다. As a method of arranging the ultrasound lines, a method of arranging the ultrasound lines so as to be suitable for the Chebbyshev integration is suitable as a method of calculating the position in which the ultrasound sensor 30 is installed.

초음파 거리 측정 방법론은, 대부분 초음파 유량계(10)는 초음파 센서(30) 사이의 설치 거리를 정확하게 알고 있어야 유속과 음속을 측정할 수 있다. 만약 초음파 센서(30) 설치 거리가 정확하게 입력되지 않으면, 측정 오차가 증가한다. 측정 오차를 줄이기 위해, 외벽부착식 초음파 유량계(10)를 설정할 때는, 배관(11) 두께를 입력하는 단계가 추가되기도 한다. 유량 환산 공식은 초음파 센서(30) 설치 거리를 가지고 유속과 음속을 계산한다. 음속을 가지고 유속을 계산하도록 이 공식을 변형하면 센서 설치 거리를 측정하려 애쓰지 않아도 된다. 그 대신, 음속을 계산하기 위해 온도와 압력을 측정하는 단계가 필요하다. 온도계, 압력계, 유량계는 배관을 흐르는 유체를 모니터링하기 위해 항상 사용하므로, 음속을 가지고 유속을 계산하도록 유량 환산 공식을 변형해야 한다.In the ultrasonic distance measurement methodology, most of the ultrasonic flowmeters 10 can measure the flow velocity and sound velocity only when the installation distance between the ultrasonic sensors 30 is accurately known. If the ultrasonic sensor 30 installation distance is not correctly input, the measurement error increases. In order to reduce measurement errors, when setting the external wall-mounted ultrasonic flow meter 10, a step of inputting the thickness of the pipe 11 may be added. The flow rate conversion formula calculates the flow velocity and sound velocity with the installation distance of the ultrasonic sensor 30. By modifying this formula to calculate the flow velocity from the sound velocity, you don't have to struggle to measure the distance between the sensor installation. Instead, it takes steps to measure temperature and pressure to calculate the speed of sound. Thermometers, pressure gauges, and flow meters are always used to monitor the fluid flowing through the pipe, so the flow conversion formula must be modified to calculate the flow rate with the speed of sound.

초음파 유량계(10)가 측정하는 유속은 선-적분된 유속으로 유량 측정에 필요한 평균 유속(면적-적분된 유속)과는 차이가 있다. 유속 보정계수는 선-적분된 유속을 평균 유속으로 환산하기 위한 계수로서 평균 유속을 선-적분된 유속으로 나눈 값이다. 유속 보정계수는 여러 가지가 있는데, 유속 보정계수가 가장 잘 알려져 있다.The flow rate measured by the ultrasonic flow meter 10 is a pre-integrated flow rate, which is different from the average flow rate (area-integrated flow rate) required for flow rate measurement. The flow rate correction factor is a coefficient for converting the pre-integrated flow rate into an average flow rate, and is a value obtained by dividing the average flow rate by the pre-integrated flow rate. There are several correction factors for flow velocity, and the correction coefficient for flow velocity is the best known.

Figure 112019093202213-pat00016
(1)
Figure 112019093202213-pat00016
(One)

Figure 112019093202213-pat00017
는 유속 보정계수,
Figure 112019093202213-pat00018
는 배관 직경과 평균 유속에 대한 레이놀즈 수(
Figure 112019093202213-pat00019
)이다.
Figure 112019093202213-pat00017
Is the flow rate correction factor,
Figure 112019093202213-pat00018
Is the Reynolds number for the pipe diameter and average flow rate (
Figure 112019093202213-pat00019
)to be.

유속분포 보정계수는 유속 보정계수보다 정밀한 유속 보정방법을 제시한다. K Ioose는 배관 단면 전체에 적용되는 선-적분된 유속을 가지고 평균 유속으로 보정한다. 이와 반대로, 유속분포 보정계수는 배관(11) 단면에서 초음파 경로가 차지하는 면적만큼 적분된 선-적분된 유속을 평균 유속으로 환산한다.

Figure 112019093202213-pat00020
는 초음파 센서 설치위치에 상관없이 똑같은 보정계수가 적용되지만, 유속분포 보정계수는 설치위치를 고려하여 서로 다른 보정계수를 적용한다.The flow velocity distribution correction factor suggests a more accurate flow velocity correction method than the flow velocity correction factor. K Ioose is corrected to the average flow rate with the pre-integrated flow rate applied to the entire pipe section. Conversely, the flow velocity distribution correction factor converts the line-integrated flow velocity integrated by the area occupied by the ultrasonic path in the cross section of the pipe 11 to the average flow velocity.
Figure 112019093202213-pat00020
The same correction factor is applied regardless of the ultrasonic sensor installation location, but different correction factors are applied for the flow velocity distribution correction factor considering the installation location.

유속분포는 축대칭 유동, 단일 피크 유동, 이중 피크 유동을 정의한 Salami 유속 함수를 수학적 모델로 사용하여 접근할 수 있다. 살라미(Salami) 유속 함수는 직관부 길이가 확보된 축대칭 유동이나 단일 곡관이 설치된 배관의 단일 피크 유동, 이중 곡관이 설치된 배관의 이중 피크 유동에 적용 가능하다. 살라미 유속 함수는 아래와 같다.The flow velocity distribution can be approached using the Salami flow velocity function, which defines axisymmetric flow, single peak flow, and double peak flow, as a mathematical model. The salami flow rate function can be applied to an axisymmetric flow with a straight pipe length secured, a single peak flow of a pipe with a single curved pipe, and a double peak flow of a pipe with a double curved pipe. The salami flow rate function is as follows.

Figure 112019093202213-pat00021
(2)
Figure 112019093202213-pat00021
(2)

단,

Figure 112019093202213-pat00022
는 유속 (m/s),
Figure 112019093202213-pat00023
은 무차원화된 반지름 방향 좌표,
Figure 112019093202213-pat00024
는 각도 방향 좌표 (rad),
Figure 112019093202213-pat00025
Figure 112019093202213-pat00026
는 유속분포 지수,
Figure 112019093202213-pat00027
은 임의의 값,
Figure 112019093202213-pat00028
는 임의의 함수이다. 각도 적용 범위는
Figure 112019093202213-pat00029
이다. Salami 유속 함수에는 6가지 모델이 있는데 이를 정리하면 [표 1]과 같다.only,
Figure 112019093202213-pat00022
Is the flow rate (m/s),
Figure 112019093202213-pat00023
Is the dimensionless radial coordinate,
Figure 112019093202213-pat00024
Is the angular coordinate (rad),
Figure 112019093202213-pat00025
Wow
Figure 112019093202213-pat00026
Is the flow velocity distribution index,
Figure 112019093202213-pat00027
Is any value,
Figure 112019093202213-pat00028
Is an arbitrary function. The angle coverage is
Figure 112019093202213-pat00029
to be. There are six models of the Salami flow rate function, which are summarized in [Table 1].

살라미 유속 함수 모델 Salami flow rate function model ProfileProfile ModelModel

Figure 112019093202213-pat00030
Figure 112019093202213-pat00030
Figure 112019093202213-pat00031
Figure 112019093202213-pat00031
Figure 112019093202213-pat00032
Figure 112019093202213-pat00032
Figure 112019093202213-pat00033
Figure 112019093202213-pat00033
AxisymmetricAxisymmetric AA 1/13 ~ 1/51/13 ~ 1/5 00 00 00 One peakOne peak BB 1/91/9 -0.5/
Figure 112019093202213-pat00034
-0.5/
Figure 112019093202213-pat00034
1/41/4
Figure 112019093202213-pat00035
Figure 112019093202213-pat00035
CC 1/91/9 0.04/
Figure 112019093202213-pat00036
0.04/
Figure 112019093202213-pat00036
1/41/4
Figure 112019093202213-pat00037
Figure 112019093202213-pat00037
DD 1/91/9 3.31703.3170 22
Figure 112019093202213-pat00038
Figure 112019093202213-pat00038
Two peaksTwo peaks EE 1/91/9 1/
Figure 112019093202213-pat00039
One/
Figure 112019093202213-pat00039
1/41/4
Figure 112019093202213-pat00040
Figure 112019093202213-pat00040
FF 1/91/9 0.68130.6813 1/91/9
Figure 112019093202213-pat00041
Figure 112019093202213-pat00041

살라미 유속 함수를 가지고 유속분포 보정계수를 계산하려면 도 11과 같이 적분 영역을 설정해야 한다. In order to calculate the correction factor for the flow velocity distribution using the salami flow velocity function, an integral region must be set as shown in FIG. 11.

빗금친 면적은 초음파 센서(30) 설치 위치에서 주고받는 초음파가 차지하는 면적이다. 초음파 센서(30)가 빗금친 면적이 배관(11)의 표면과 닿는 지점에 설치된다면, 원형(

Figure 112019093202213-pat00042
-
Figure 112019093202213-pat00043
) 좌표에서 식 (3)처럼 빗금친 면적에 해당하는 유속 분포를 적분하는 식을 나타낼 수 있다. 만약 빗금친 면적이 배관(11) 단면에 걸쳐 4개가 있다면, 이는 습식 4회선 초음파 유량계(10)에서 초음파가 주고받는 면적을 표현하는 셈이다.The hatched area is an area occupied by ultrasonic waves transmitted and received at the installation location of the ultrasonic sensor 30. If the ultrasonic sensor 30 is installed at a point where the hatched area contacts the surface of the pipe 11, a circular (
Figure 112019093202213-pat00042
-
Figure 112019093202213-pat00043
) In the coordinates, the equation for integrating the distribution of the flow velocity corresponding to the hatched area can be expressed as in equation (3). If there are four hatched areas across the cross-section of the pipe 11, this represents the area of ultrasonic waves exchanged in the wet four-wire ultrasonic flow meter 10.

Figure 112019093202213-pat00044
Figure 112019093202213-pat00044

Figure 112019093202213-pat00045
(3)
Figure 112019093202213-pat00045
(3)

Figure 112019093202213-pat00046
Figure 112019093202213-pat00046

단,

Figure 112019093202213-pat00047
은 선-적분된 유속 (m/s),
Figure 112019093202213-pat00048
은 적분 영역 (m2),
Figure 112019093202213-pat00049
은 배관 반지름 (m),
Figure 112019093202213-pat00050
는 최대 유속 (m/s),
Figure 112019093202213-pat00051
은 유속분포 보정계수이다.only,
Figure 112019093202213-pat00047
Is the line-integrated flow rate (m/s),
Figure 112019093202213-pat00048
Is the integral area (m 2 ),
Figure 112019093202213-pat00049
Is the pipe radius (m),
Figure 112019093202213-pat00050
Is the maximum flow rate (m/s),
Figure 112019093202213-pat00051
Is the correction factor for the flow velocity distribution.

식 (3)에서 주의할 점은, 유속분포 보정계수는 선-적분된 유속을 평균 유속(면적-적분된 유속)으로 보정해야 하므로,

Figure 112019093202213-pat00052
대신
Figure 112019093202213-pat00053
이 적용되어야 한다. 이 부분은 아직 검토 중으로 향후
Figure 112019093202213-pat00054
가 적용된 유속분포 보정계수를 제공하고자 한다. 식 (3)을 살라미 유속분포에 적용하려면 가능한 한 손으로 풀 수 있는 데까지 식을 전개하고 나서 수치해석을 시도하면 좋다. 식 (3)의 오른쪽 첫 번째 항을 손으로 풀 수 있는 데까지 풀면 아래와 같다.Note that in Equation (3), the correction factor for the flow velocity distribution must correct the pre-integrated flow velocity to the average velocity (area-integrated flow velocity).
Figure 112019093202213-pat00052
instead
Figure 112019093202213-pat00053
Should be applied. This part is still under review.
Figure 112019093202213-pat00054
It is intended to provide a correction factor for the flow velocity distribution to which is applied. To apply Equation (3) to the salami flow velocity distribution, it is good to develop the equation as far as possible to solve it by hand and then try numerical analysis. Solving the first term on the right of Equation (3) as far as it can be solved by hand is as follows.

Figure 112019093202213-pat00055
(4)
Figure 112019093202213-pat00055
(4)

단, 식 (4)의 오른쪽 첫 번째 적분식은 식 (2)의 오른쪽 첫 번째 항에 대응되고, 식 (4)의 오른쪽 두 번째 적분식은 식 (2)의 오른쪽 두 번째 항에 대응된다. 식 (3)과 비교하여 식 (4)가 가지는 장점은, 이중 적분을 단일 적분으로 변경하므로 일반적인 수치 적분 함수를 적용하기 쉽다는 점이다. 잘 알려진 매트랩(MATLAB)이나 C 언어 수치 적분 알고리즘을 적용하면 식 (4)를 쉽게 적분할 수 있다. 매트랩에서 상기 언급한 방법대로 프로그램으로 만든 다음, [표 1]에 제시된 순서대로

Figure 112019093202213-pat00056
를 계산하면 도 12 내지 도 17에 도시한 바와 같다.However, the first integral expression on the right of Equation (4) corresponds to the first term on the right of Equation (2), and the second integral expression on the right of Equation (4) corresponds to the second term on the right of Equation (2). The advantage of Equation (4) compared to Equation (3) is that it is easy to apply a general numerical integral function because it converts a double integral into a single integral. Equation (4) can be easily integrated by applying the well-known MATLAB or C language numerical integration algorithm. After creating the program in the manner mentioned above in MATLAB, in the order shown in [Table 1]
Figure 112019093202213-pat00056
Calculating is as shown in FIGS. 12 to 17.

도 12 내지 도 17에 제시된 유속분포 보정계수는, 초음파 센서(10)를 12시 방향으로 설치한다고 가정하고, 살라미(Salami) 유속 함수를 적분한 결과이다. 만약 초음파 센서를 1시나 3시, 5시 방향에 설치한다고 가정하면, 식 (2)에 나타낸

Figure 112019093202213-pat00057
Figure 112019093202213-pat00058
로 변경하고 수치 적분하면 계산할 수 있다. 이 때,
Figure 112019093202213-pat00059
는 초음파 센서(30)가 설치되는 각도이다. 만약 초음파 센서(30)를 12시 방향에 설치하면,
Figure 112019093202213-pat00060
는 0이 된다. 만약 초음파 센서(30)를 3시 방향에 설치하면,
Figure 112019093202213-pat00061
는 90도(
Figure 112019093202213-pat00062
)가 된다. 살라미 유속 함수를 시계 반대방향으로 회전하면, 유속분포 보정계수는 초음파 센서(30)를 시계 방향으로 회전시켜 수치 적분하는 것과 같은 효과를 낸다. 단, 유체가 흐르는 방향은 배관(11) 단면을 정면에서 바라볼 때, 유체가 뒷면에서 앞면으로 흐르는 방향을 기준으로 한다.The correction factor for the flow velocity distribution shown in FIGS. 12 to 17 is the result of integrating the Salami flow velocity function, assuming that the ultrasonic sensor 10 is installed in the 12 o'clock direction. Assuming that the ultrasonic sensor is installed at 1 o'clock, 3 o'clock, or 5 o'clock, the equation (2)
Figure 112019093202213-pat00057
To
Figure 112019093202213-pat00058
It can be calculated by changing to and numerical integration. At this time,
Figure 112019093202213-pat00059
Is the angle at which the ultrasonic sensor 30 is installed. If the ultrasonic sensor 30 is installed at 12 o'clock,
Figure 112019093202213-pat00060
Becomes 0. If the ultrasonic sensor 30 is installed at 3 o'clock,
Figure 112019093202213-pat00061
Is 90 degrees (
Figure 112019093202213-pat00062
). When the salami flow rate function is rotated counterclockwise, the flow rate distribution correction factor has the same effect as numerical integration by rotating the ultrasonic sensor 30 clockwise. However, the direction in which the fluid flows is based on the direction in which the fluid flows from the rear to the front when the cross section of the pipe 11 is viewed from the front.

초음파 센서(30) 설치위치는, 초음파 유량계(10)는 직관부 길이가 충분하게 확보된 상태에서 완전히 발달한 난류 유속분포(fully-developed turbulent flow velocity profile)를 가정하고 제작한다. 그러므로 초음파 유량계(10)가 필요로 하는 초음파 센서(30) 설치 위치는 거의 대부분 완전히 발달한 난류 유속분포를 기준으로 한다. 만약 초음파 회선 수가 매우 많으면, 완전히 발달한 난류 유속분포를 수치 적분하듯이 계산하여 유량을 계산할 수 있다. 그렇지만 초음파 회선 수는 4회선이나 5회선 등으로 한정되기 때문에, 효과적인 수치 적분 방법을 찾아야 한다.The ultrasonic sensor 30 is installed on the assumption that the ultrasonic flow meter 10 has a fully-developed turbulent flow velocity profile in a state in which the straight pipe length is sufficiently secured. Therefore, the installation position of the ultrasonic sensor 30 required by the ultrasonic flow meter 10 is almost entirely based on a fully developed turbulent flow velocity distribution. If the number of ultrasonic lines is very large, the flow rate can be calculated by calculating the fully developed turbulent flow velocity distribution as a numerical integration. However, since the number of ultrasonic lines is limited to 4 or 5 lines, it is necessary to find an effective numerical integration method.

가우스 적분(Gauss quadrature)이나 쳬비셰프 적분(Chebyshev quadrature)는 몇 개의 측정값을 가지고 유속분포를 적분하기 적합하다. 초음파 유량계(10)에서 가장 많이 사용하는 적분 방법은 가우스 적분이리라 예상된다. 이번 기술자문에서는 가중치가 1로 일정한 쳬비셰프 적분을 추천하고자 한다.The Gauss quadrature or Chebyshev quadrature is suitable for integrating the velocity distribution from several measurements. It is expected that the most commonly used integration method in the ultrasonic flow meter 10 is Gaussian integration. In this technical advisory, I would like to recommend the Chebbyshev integral with a constant weight of 1.

가우스 적분을 먼저 설명하자면, 가우스 적분은 배관 반지름이 1로 무차원화된 형상에 적용할 수 있다. 만약 4개 지점을 선택하려면, [-1, 1]을 기준으로 [-0.861136, -0.339981, 0.339981, 0.861136]인 무차원화된 위치에 초음파 회선을 배열하면 된다. 상기 제시된 값에 배관(11) 반지름을 곱하면, 초음파 회선이 중심축으로부터 떨어진 거리를 계산할 수 있다. 여기에는 가중치가 필요하다. 가중치란, 각 초음파 회선이 측정한 유속에 곱하는 0보다 크고 1보다 작은 계수이다. 4개 지점에 각각 대응하는 가중치는 [0.347855, 0.652145, 0.652145, 0.347855]이다. 이를 수학 공식으로 표현하면 아래와 같다.To explain the Gaussian integral first, the Gaussian integral can be applied to a dimensionless shape with a pipe radius of 1. If you want to select 4 points, you can arrange the ultrasonic lines in the dimensionless position of [-0.861136, -0.339981, 0.339981, 0.861136] based on [-1, 1]. When the suggested value is multiplied by the radius of the pipe 11, the distance of the ultrasonic line away from the central axis can be calculated. Here we need weights. The weight is a coefficient greater than 0 and less than 1 multiplied by the flow velocity measured by each ultrasonic line. The weights corresponding to the four points are [0.347855, 0.652145, 0.652145, 0.347855]. This can be expressed as a mathematical formula as follows.

Figure 112019093202213-pat00063
(5)
Figure 112019093202213-pat00063
(5)

단,

Figure 112019093202213-pat00064
은 배관 단면을 지나는 평균 유속 (m/s)이고,
Figure 112019093202213-pat00065
는 가중치,
Figure 112019093202213-pat00066
는 초음파 회선이 측정하는 선-적분된 유속 (m/s)이다.only,
Figure 112019093202213-pat00064
Is the average flow velocity (m/s) through the pipe section,
Figure 112019093202213-pat00065
Is the weight,
Figure 112019093202213-pat00066
Is the line-integrated flow velocity (m/s) measured by the ultrasonic convolution.

쳬비셰프 적분은 가우스 적분과 비슷하나 초음파 회선에 상관없이 가중치가 일정하다는 점이 다르다. 비교하자면, 이번 기술자문에서 설명하려는 가우스 적분의 정확한 명칭은 르장드-가우스 적분(Legendre-Gauss quadrature)로서 [-1, 1]을 구간으로 하는 임의의 함수를 적분할 수 있다. 이와는 반대로, 쳬비셰프 적분은 [-1, 1] 구간에서 양 끝점이 0인 임의의 함수를 적분할 수 있다. 배관 벽면에서 유속이 0이 되는 no-slip 조건은 [-1, 1] 구간의 양 끝점이 0이라는 것을 가리키므로, 쳬비셰프 적분을 적용하는 편이 더 유리하다.The Chebbyshev integral is similar to the Gaussian integral, except that the weight is constant regardless of the ultrasound line. For comparison, the exact name of the Gaussian integral to be described in this article is Legend-Gauss quadrature, which can integrate arbitrary functions with intervals [-1, 1]. Conversely, the Chebbyshev integration can integrate an arbitrary function whose end points are 0 in the interval [-1, 1]. Since the no-slip condition in which the flow velocity is 0 at the pipe wall indicates that both ends of the [-1, 1] section are 0, it is more advantageous to apply the Chebbyshev integral.

만약 가우스 적분과 마찬가지로 4개 지점을 선택하려면, [-1, 1]을 기준으로 [-0.794654, -0.187592, 0.187592, 0.794654]인 무차원화된 위치에 초음파 회선을 배열하면 된다. 가중치는 2/(초음파 회선 개수) = 0.5이다. 이를 정리하면 [표 2]와 같다. 그리고 초음파 회선이 5개일 때, 유속분포 적분값이 얼마나 정확한가 확인하기 위해 [표 2]와 [표 3]에 초음파 회선과 유속분포 적분값을 정리하였다.If you want to select 4 points like Gaussian integration, you can arrange the ultrasonic lines at the dimensionless position of [-0.794654, -0.187592, 0.187592, 0.794654] based on [-1, 1]. The weight is 2/(number of ultrasonic lines) = 0.5. This is summarized in [Table 2]. In addition, in order to check how accurate the flow velocity distribution integral value is when there are 5 ultrasonic lines, the ultrasound line and velocity distribution integral values are summarized in [Table 2] and [Table 3].

가우스 적분과 쳬비셰프 적분 비교 (4회선)Comparison of Gaussian Integral and Chebbyshev Integral (4 lines) 적분 방법Integration method 회선 1Line 1 회선 2Line 2 회선 3Line 3 회선 4Line 4 가우스 적분Gaussian integral 위치location -0.861136-0.861136 -0.339981-0.339981 0.3399810.339981 0.8611360.861136 가중치weight 0.3478550.347855 0.6521450.652145 0.6521450.652145 0.3478550.347855 쳬비셰프 적분Chembyshev integral 위치location -0.794654-0.794654 -0.187592-0.187592 0.1875920.187592 0.7946540.794654 가중치weight 0.50.5 0.50.5 0.50.5 0.50.5

가우스 적분과 쳬비셰프 적분 비교 (5회선) Comparison of Gaussian Integral and Chebbyshev Integral (5 lines) 적분 방법Integration method 회선 1Line 1 회선 2Line 2 회선 3Line 3 회선 4Line 4 회선 5Line 5 가우스 적분Gaussian integral 위치location -0.906180-0.906180 -0.538469-0.538469 00 0.5384690.538469 0.9061800.906180 가중치weight 0.2369270.236927 0.4786290.478629 0.5688890.568889 0.4786290.478629 0.2369270.236927 쳬비셰프 적분Chembyshev integral 위치location -0.832497-0.832497 -0.374541-0.374541 00 0.3745410.374541 0.8324970.832497 가중치weight 0.40.4 0.40.4 0.40.4 0.40.4 0.40.4

유속분포 적분값 비교 (5회선) Flow velocity distribution integral value comparison (5 lines) 속도 적분값 (이론식)Velocity integral value (theoretical formula) 가우스 적분 (5회선)Gaussian integral (5 lines) 쳬비셰프 적분 (5회선)Cheambishev integration (5 lines) 쳬비셰프 적분 (6회선)Cheambishev integral (6 lines) 1.751.75 1.7639731.763973 1.7679051.767905 1.7542921.754292

[표 3]에서 제시한 속도 적분값은 이론식에 따른 것으로 식 (2)의 오른쪽 첫째 항만 고려하고

Figure 112019093202213-pat00067
= 1/7이면, 속도 적분값은 1.75이다. 식 (2)의 오른쪽 첫째 항을 표 3에 나타낸 것처럼 5회선으로 배열하여 가우스 적분을 계산하면, 1.763943으로 상대 오차는 0.8 %이다. 5회선으로 배열하여 쳬비셰프 적분을 계산하면, 1.767905로 상대 오차는 1.02 %이다. 만약 6회선으로 배열하여 쳬비셰프 적분을 계산하면, 1.754292로 상대 오차는 0.25 %가 된다. 결국 예상과 달리, 초음파 회선 개수가 똑같으면, 쳬비셰프 적분보다 가우스 적분이 더 유리하다. 만약 초음파 회선 개수를 1회선 더 추가하면, 쳬비셰프 적분의 정확도는 향상될 수 있다. 그럼에도, 쳬비셰프 적분이 제안하는 장점은, 배관 벽면 속도가 0이라는 no-slip 조건이 물리적으로 맞고 가중치가 동일하므로 선-적분된 유속을 평균유속으로 환산하는 알고리즘이 쉬워진다는 점이다.The velocity integral value presented in [Table 3] is according to the theoretical equation, considering only the first term on the right of equation (2)
Figure 112019093202213-pat00067
= 1/7, the velocity integral is 1.75. If the Gaussian integral is calculated by arranging the first term on the right of Equation (2) in 5 lines as shown in Table 3, the relative error is 1.763943, which is 0.8%. If 5 lines are arranged and the Chebbyshev integral is calculated, the relative error is 1.767905, which is 1.02%. If the Chebbyshev integral is calculated by arranging it in 6 lines, the relative error is 0.25% at 1.754292. In the end, contrary to expectations, if the number of ultrasound lines is the same, the Gaussian integral is more advantageous than the Chebbyshev integral. If the number of ultrasonic lines is added by one more, the accuracy of the Chebbyshev integration can be improved. Nevertheless, the advantage proposed by the Chebbyshev integral is that the no-slip condition that the pipe wall velocity is 0 is physically met and the weights are the same, so the algorithm for converting the pre-integrated flow velocity into the average flow velocity becomes easier.

유량 환산공식은, 초음파 유량계(30)를 규정하는 ISO 규격은 액체용 초음파 유량계에 관한 ISO 12242:2012, 기체용 초음파 유량계에 관한 ISO 17089-1:2010, ISO 17089-2:2012가 있다. ISO 12242에 제시된 유량 공식을 정리하면 아래와 같다.As for the flow rate conversion formula, the ISO standards for specifying the ultrasonic flow meter 30 include ISO 12242:2012 for liquid ultrasonic flow meters, ISO 17089-1:2010 and ISO 17089-2:2012 for ultrasonic flow meters for gas. The flow rate formula presented in ISO 12242 is summarized below.

Figure 112019093202213-pat00068
(6)
Figure 112019093202213-pat00068
(6)

Figure 112019093202213-pat00069
Figure 112019093202213-pat00069

단,

Figure 112019093202213-pat00070
는 초음파 경로 길이 (m),
Figure 112019093202213-pat00071
는 초음파 경로와 유체흐름 방향 사이의 각도 (rad),
Figure 112019093202213-pat00072
는 초음파 도달시간 차이 (s),
Figure 112019093202213-pat00073
는 초음파 유량계에서 신호 처리하는 데 필요한 시간이나 초음파가 센서(30) 재질을 통과하는 시간을 모두 포함한 측정시간 (s),
Figure 112019093202213-pat00074
은 상류측 초음파 도달시간 (s),
Figure 112019093202213-pat00075
은 하류측 초음파 도달시간 (s),
Figure 112019093202213-pat00076
는 초음파 진행 속도 (m/s) 이다. only,
Figure 112019093202213-pat00070
Is the ultrasonic path length (m),
Figure 112019093202213-pat00071
Is the angle (rad) between the ultrasonic path and the direction of fluid flow,
Figure 112019093202213-pat00072
Is the ultrasonic arrival time difference (s),
Figure 112019093202213-pat00073
Is the measurement time (s) including the time required to process the signal in the ultrasonic flow meter or the time that the ultrasonic wave passes through the material of the sensor 30,
Figure 112019093202213-pat00074
Is the upstream ultrasonic arrival time (s),
Figure 112019093202213-pat00075
Is the downstream ultrasonic arrival time (s),
Figure 112019093202213-pat00076
Is the ultrasonic propagation speed (m/s).

ISO 17089-1과 ISO 17089-2에 제시된 유량 공식을 정리하면 아래와 같다. The flow formulas presented in ISO 17089-1 and ISO 17089-2 are summarized below.

Figure 112019093202213-pat00077
(7)
Figure 112019093202213-pat00077
(7)

Figure 112019093202213-pat00078
Figure 112019093202213-pat00078

식 (6)과 식 (7)은 같은 식을 다르게 변형했을 뿐으로, ISO 12242와 ISO 17089-1, ISO 17089-2는 초음파 도달시간과 초음파 경로 길이를 정확하게 측정한다는 가정에서 출발한다. 만약 초음파 경로 길이를 정확하게 측정하기 어렵다면 초음파 진행 속도와 유속을 정확하게 측정할 수 없다. 이와는 반대로, 온도와 압력을 가지고 초음파 진행 속도를 정확하게 측정하면, 유속을 정확하게 측정할 수 있다. 이를 식으로 표현하면 아래와 같다.Equations (6) and (7) are only modified differently from the same equation, and ISO 12242, ISO 17089-1, and ISO 17089-2 start from the assumption that the ultrasonic arrival time and the ultrasonic path length are accurately measured. If it is difficult to accurately measure the ultrasonic path length, it is not possible to accurately measure the ultrasonic propagation speed and flow rate. On the contrary, by accurately measuring the speed of ultrasonic waves with temperature and pressure, the flow rate can be accurately measured. This can be expressed as an equation as follows.

Figure 112019093202213-pat00079
(8)
Figure 112019093202213-pat00079
(8)

식 (8)은 초음파 진행 속도, 초음파 도달시간, 초음파 경로 각도가 정확하면 유속을 정확하게 계산할 수 있다고 알려준다. 초음파 진행 속도는 온도계와 압력계를 배관 내부를 지나는 유체 온도와 압력을 측정하면, 참조표준 데이터로부터 계산될 수 있다. NIST REFPROP (SRD 23)은 미국 NIST에서 개발한 유체 열물성 데이터베이스로서 온도와 압력을 입력하면 특정한 유체의 열물성(밀도, 점도, 음속 등)을 계산할 수 있다. 사실, 유량 측정이 필요한 배관(11)에는 온도계와 압력계가 공정 모니터링용으로 설치되기 때문에 온도계-압력계-유량계가 통합된 측정기기가 있다면, 식 (8)을 구현하기가 편리하다.Equation (8) tells us that the flow velocity can be accurately calculated if the ultrasonic propagation speed, ultrasonic arrival time, and ultrasonic path angle are correct. The ultrasonic propagation speed can be calculated from reference standard data by measuring the temperature and pressure of the fluid passing through the pipe with a thermometer and a manometer. NIST REFPROP (SRD 23) is a fluid thermal property database developed by NIST in the United States. By entering temperature and pressure, the thermal properties (density, viscosity, sound velocity, etc.) of a specific fluid can be calculated. In fact, since a thermometer and a pressure gauge are installed for process monitoring in the pipe 11 requiring flow rate measurement, if there is a measuring device in which a thermometer-pressure gauge-flowmeter is integrated, it is convenient to implement Equation (8).

초음파 도달시간은 초음파 유량계 컨트롤러(40)가 출력하는 측정값이다. 초음파 경로 각도는 유량계 배관(11)이 얼마나 정밀하게 제작되는가에 정확도가 달라진다. 따라서 식 (8)은 초음파 유량계 컨트롤러(40), 유량계 배관(11), 유체 특성을 모두 고려해야 유속을 측정할 수 있다는 점이다.The ultrasonic arrival time is a measured value output by the ultrasonic flow meter controller 40. The accuracy of the ultrasonic path angle varies depending on how precisely the flow meter pipe 11 is manufactured. Therefore, Equation (8) is that the flow rate can be measured considering all of the ultrasonic flow meter controller 40, the flow meter pipe 11, and fluid characteristics.

유체의 압력, 온도에 따른 초음파 이동속도 계산 방정식은, 음향 임피던스는 유체 밀도와 음속을 곱한 값이다. 음향 임피던스는 온도와 압력의 함수로 배관 내부를 지나는 유체 특성에 따라 달라진다. In the equation for calculating the ultrasonic movement velocity according to the pressure and temperature of the fluid, the acoustic impedance is the product of the fluid density and the sound velocity. Acoustic impedance depends on the characteristics of the fluid passing through the pipe as a function of temperature and pressure.

CO2, N2, air에 대한 음향 임피던스를 계산하기 위해서는 온도와 압력 범위를 정해야 한다. 온도 범위를 0 ℃ ~ 200 ℃, 압력 범위를 0.1 MPa (1 기압) ~ 3 MPa (30 기압)로 정했다. 단, 압력은 절대 압력으로 대기압은 0.1 MPa과 같다. 온도 범위를 21등분, 압력 범위를 30등분으로 나누어 계산했다. In order to calculate the acoustic impedance for CO 2 , N 2 , and air, the temperature and pressure ranges must be defined. The temperature range was set at 0°C to 200°C, and the pressure range was set at 0.1 MPa (1 atm) to 3 MPa (30 atm). However, the pressure is the absolute pressure and the atmospheric pressure is equal to 0.1 MPa. It was calculated by dividing the temperature range by 21 equal parts and the pressure range by 30 equal parts.

CO2, N2, air의 음향 임피던스는 온도와 압력이 증가할수록 증가한다. 상온, 대기압 조건에서 초음파 유량계 성능을 실험한 결과보다는 고온 고압 조건에서 초음파 유량계 성능을 실험한 결과가 더 좋을 것으로 예측할 수 있다.The acoustic impedance of CO 2 , N 2 , and air increases with increasing temperature and pressure. It can be predicted that the results of testing the performance of the ultrasonic flowmeter under high temperature and high pressure conditions will be better than the results of testing the performance of the ultrasonic flowmeter under room temperature and atmospheric pressure conditions.

만약 음향 임피던스가 수백 kg/ms 정도로 낮으면, 초음파 주파수가 낮을수록 초음파가 더 잘 전달된다. 그 이유는 초음파를 구동하는 데 필요한 에너지는 초음파 진폭과 주파수를 곱한 값에 비례하기 때문이다. 만약 에너지가 일정하다면, 초음파 주파수가 낮을수록 진폭이 증가한다. 초음파 진폭은 유체를 밀어내는 힘의 크기에 비례한다. 그러므로 초음파 주파수가 낮아질수록 유체를 밀거나 당기는 힘의 크기가 증가하고 이는 초음파 신호가 더 잘 전달됨을 의미한다.If the acoustic impedance is as low as several hundred kg/m s, the lower the ultrasonic frequency, the better the ultrasonic waves are transmitted. This is because the energy required to drive ultrasonic waves is proportional to the product of ultrasonic amplitude times frequency. If the energy is constant, the amplitude increases as the ultrasonic frequency decreases. The ultrasonic amplitude is proportional to the magnitude of the force pushing the fluid. Therefore, as the ultrasonic frequency decreases, the force pushing or pulling the fluid increases, which means that the ultrasonic signal is better transmitted.

결과만으로는 초음파 주파수와 음향 임피던스의 정량적인 인과 관계를 알 수 없다. 그 이유는 초음파 주파수와 음향 임피던스 사이의 인과 관계를 제시하는 수학적 모델은 실험에 의존하기 때문이다. 예를 들어, 한 쌍의 초음파 센서가 배관(11) 내부를 지나는 유체 흐름을 사이에 두고 초음파 신호를 주고받을 때, 초음파 주파수에 따라 진폭이 얼마나 변화하는가를 오실로스코프로 측정하면, 초음파 주파수와 음향 임피던스를 연관 짓는 수학적 모델을 구할 수 있다.The results alone do not reveal the quantitative causal relationship between ultrasonic frequency and acoustic impedance. The reason is that the mathematical model that suggests the causal relationship between the ultrasonic frequency and the acoustic impedance is experimental. For example, when a pair of ultrasonic sensors exchanges ultrasonic signals with a fluid flow passing through the inside of the pipe 11, measuring how much the amplitude changes according to the ultrasonic frequency with an oscilloscope, the ultrasonic frequency and acoustic impedance A mathematical model can be obtained that relates to.

1차 시험 결과 500㎥/h∼4000㎥/h 구간에서 오차가 ±0.6% 사이로 측정되었다. 1차 시험과 동일한 시료를 사용하고 1차 테스트(Test)를 바탕으로 멀티 팩터(Multi factor)를 사용하여 보정하였다.As a result of the first test, the error was measured within ±0.6% in the 500㎥/h∼4000㎥/h section. The same sample as the first test was used and corrected using a multi factor based on the first test.

2차 시험 결과 500㎥/h∼4000㎥/h 구간에서 오차가 ±0.3% 사이로 1차 시험시 보다 오차가 2배 개선 되었다.As a result of the second test, the error was within ±0.3% in the 500㎥/h∼4000㎥/h section, and the error was improved by 2 times compared to the first test.

이러한 결과를 바탕으로 초음파 유량계(10)의 배관(11)에 구멍(13)을 성형해 새들(20)을 설치한 후 초음파 센서(30)를 설치하는 위치와 각도 및 채널을 통한 초음파 행정(50)이 다양하게 설정하여 적용함으로써 배관(11) 내 유체 유속 분포 입체적 측정을 위한 초음파 행정(50)이 격자 행정의 다중센서 구조를 적용함으로써 유체 흐름을 3D로 모사하여 유량 계측이 가능하도록 하였다.Based on these results, after installing the saddle 20 by forming a hole 13 in the pipe 11 of the ultrasonic flow meter 10, the ultrasonic stroke 50 through the location and angle and the channel where the ultrasonic sensor 30 is installed. ) Is applied in various ways, so that the ultrasonic stroke 50 for three-dimensional measurement of the distribution of fluid flow velocity in the pipe 11 applies a multi-sensor structure of a grid stroke, thereby simulating the fluid flow in 3D and measuring the flow rate.

이러한 기술은 PC 보드로 고속의 ADC를 사용해 초음파 속도와 파형 데이터를 획득할 수 있는 타이머를 제공하고, 유체 테스트를 진행하여 수집된 데이터를 기준으로 방정식을 도출, 하드웨어에 알고리즘을 이식할 수 있어서 초음파 유량계측 기술을 발전시키게 되는 것이다.This technology provides a timer to acquire ultrasonic speed and waveform data using a high-speed ADC on a PC board, and derives an equation based on the collected data through a fluid test, and can implant an algorithm into the hardware. It is to develop flow measurement technology.

본 발명은 초음파 다중센서로 경로를 계산하여 유체 흐름을 3D로 모사 계량계측이 가능하도록 하고, PC보드로 고속의 ADC를 사용해 초음파 속도와 파형 데이터를 획득할 수 있는 타이머 및 유체 테스트를 진행하여 수집된 데이터를 기준으로 방정식을 도출해 고정밀도 초음파 유량계를 제공하는 매우 유용한 발명이다.The present invention calculates a path with an ultrasonic multi-sensor to simulate the flow of fluid in 3D, measure and measure the fluid flow, and conduct a timer and fluid test to acquire ultrasonic velocity and waveform data using a high-speed ADC with a PC board. It is a very useful invention to provide a high-precision ultrasonic flow meter by deriving an equation based on the obtained data.

10 : 초음파 유량계 11 : 배관
12 : 중공 13 : 구멍
14 : 플랜지 20 : 새들
30 : 초음파 센서 40 : 컨트롤러
20 : 센서부 50 : 초음파 행정
10: ultrasonic flow meter 11: piping
12: hollow 13: hole
14: flange 20: saddle
30: ultrasonic sensor 40: controller
20: sensor unit 50: ultrasonic stroke

Claims (2)

내측을 통해 유체가 흐를 수 있도록 중공(12)형으로 이루어진 양단에는 배관(11)의 유로와 연결을 위한 플랜지(14)가 형성되어 있고, 상기 배관(11)에 형성한 구멍(13)에 새들(20)이 배관(11)의 축방향에 대하여 경사지게 배치되어 초음파 센서(30)가 삽입되어 상호 마주하게 설치되며, 초음파 센서(30)는 컨트롤러(40)와 전기적으로 연결되어 제어되도록 설치하는 것을 포함하고,
상기 새들(20)에 초음파 센서(30)를 결합해 두 개의 초음파 센서(30)가 상호 마주보도록 비규칙적인 격자 형태가 되어 초음파 행정(50)을 갖도록 설치하며, 상기 초음파 센서(30)는 중공(12)의 유체 유속 분포를 입체적으로 측정하기 위한 구조를 제공하되;
상기 초음파 센서(30)의 격자 행정을 통하여 500㎥/h∼4000㎥/h 구간에서 멀티 팩터를 사용해 오차가 ±0.3% 사이가 되도록 하는 것을 포함하는 초음파 다중센서를 활용한 실시간 쓰리디 유량계.
Flanges 14 for connection with the flow path of the pipe 11 are formed at both ends of the hollow 12 shape so that fluid can flow through the inside, and a saddle in the hole 13 formed in the pipe 11 20 is disposed inclined with respect to the axial direction of the pipe 11 so that the ultrasonic sensor 30 is inserted and installed to face each other, and the ultrasonic sensor 30 is electrically connected to and controlled by the controller 40. Including,
An ultrasonic sensor 30 is coupled to the saddle 20 to form an irregular grid so that the two ultrasonic sensors 30 face each other and are installed to have an ultrasonic stroke 50, and the ultrasonic sensor 30 is hollow Provides a structure for three-dimensionally measuring the fluid flow velocity distribution of (12);
A real-time 3D flow meter using an ultrasonic multi-sensor comprising using a multi-factor to make an error between ±0.3% in the interval of 500㎥/h to 4000㎥/h through the grid stroke of the ultrasonic sensor 30.
삭제delete
KR1020190112206A 2019-09-10 2019-09-10 Real-Time 3D Flow Meter with Ultrasonic Multiple Sensors KR102183549B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020190112206A KR102183549B1 (en) 2019-09-10 2019-09-10 Real-Time 3D Flow Meter with Ultrasonic Multiple Sensors

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020190112206A KR102183549B1 (en) 2019-09-10 2019-09-10 Real-Time 3D Flow Meter with Ultrasonic Multiple Sensors

Publications (1)

Publication Number Publication Date
KR102183549B1 true KR102183549B1 (en) 2020-11-26

Family

ID=73679326

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020190112206A KR102183549B1 (en) 2019-09-10 2019-09-10 Real-Time 3D Flow Meter with Ultrasonic Multiple Sensors

Country Status (1)

Country Link
KR (1) KR102183549B1 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112683354A (en) * 2021-01-21 2021-04-20 沈阳佳德联益能源科技股份有限公司 NB-IoT (network B-IoT) valve control ultrasonic water meter

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100560364B1 (en) 2000-09-22 2006-03-21 신민철 apparatus for measuring the amount of flowing water using ultrasonic sensor
KR100784415B1 (en) 2006-04-14 2007-12-11 전기찬 flow meter using ultrasonic-sensor
KR101218780B1 (en) 2010-12-07 2013-01-04 (주)씨엠엔텍 Corrosion-resistant sensor and Ultrasonic flowmeter using the same
KR101267253B1 (en) * 2013-01-18 2013-05-27 가부시키가이샤 소닉 Ultrasonic flow measurement system
KR101897601B1 (en) 2017-02-21 2018-09-12 (주)씨엠엔텍 Ultrasonic transducer and ultrasonic flowmeter using the same

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100560364B1 (en) 2000-09-22 2006-03-21 신민철 apparatus for measuring the amount of flowing water using ultrasonic sensor
KR100784415B1 (en) 2006-04-14 2007-12-11 전기찬 flow meter using ultrasonic-sensor
KR101218780B1 (en) 2010-12-07 2013-01-04 (주)씨엠엔텍 Corrosion-resistant sensor and Ultrasonic flowmeter using the same
KR101267253B1 (en) * 2013-01-18 2013-05-27 가부시키가이샤 소닉 Ultrasonic flow measurement system
KR101897601B1 (en) 2017-02-21 2018-09-12 (주)씨엠엔텍 Ultrasonic transducer and ultrasonic flowmeter using the same

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112683354A (en) * 2021-01-21 2021-04-20 沈阳佳德联益能源科技股份有限公司 NB-IoT (network B-IoT) valve control ultrasonic water meter

Similar Documents

Publication Publication Date Title
RU2478919C2 (en) Calibration method under operating conditions of ultrasonic flow metres - counters of flow rate and volume of liquid single-phase media
US6209388B1 (en) Ultrasonic 2-phase flow apparatus and method
CN103808381B (en) A kind of temperature influence eliminating method of transit-time ultrasonic flow meter
JP4979838B2 (en) Method and apparatus for calibrating a measurement transducer of an ultrasonic flow measurement unit
US7506532B2 (en) Method for calibrating ultrasound clamp-on flowmeters
KR20020021558A (en) Ultrasonic wave flow measuring method
US10488237B2 (en) Measurement apparatus for measuring a flow rate of a fluid
CA2316588A1 (en) Ultrasonic 2-phase flow apparatus and stratified level detector
CN114088151B (en) External clamping type multichannel ultrasonic flow detection device and detection method
Mahadeva et al. Further studies of the accuracy of clamp-on transit-time ultrasonic flowmeters for liquids
CN106443059A (en) Fluid velocity measurement method, fluid metering method and flowmeter
Haugwitz et al. Multipath flow metering of high-velocity gas using ultrasonic phased-arrays
KR102183549B1 (en) Real-Time 3D Flow Meter with Ultrasonic Multiple Sensors
JP2002520583A (en) Multi-code flow meter
EP3063508A1 (en) A flow meter for ultrasonically measuring the flow velocity of fluids
JP2010256075A (en) Flowmeter and method of measuring flow rate
Waluś Mathematical modelling of an ultrasonic flowmeter primary device
JP2011530072A (en) Method and apparatus for determining fluid flow rate
CN106768104A (en) A kind of ultrasonic wave mass flowmenter
KR20110060536A (en) Flux measurement method
CN205861137U (en) Reducing two is popped one's head in time difference ultrasonic flow rate measurement apparatus
Sârbu Modern water flowmeters: Differential pressure flowmeters
KR100993617B1 (en) Clamp on typed multi-path ultrasonic flowmeter
CN105890684A (en) Novel setting method for determining sound channel positions by adopting Gauss-Jacobi polynomial
Gerasimov et al. Analysis of calibration and verification indirect methods of ultrasonic flowmeters

Legal Events

Date Code Title Description
GRNT Written decision to grant