KR20080039498A - 유량계의 다상 유동 물질에 대한 센서 신호를 처리하기위한 방법 및 계측 전자장치 - Google Patents

유량계의 다상 유동 물질에 대한 센서 신호를 처리하기위한 방법 및 계측 전자장치 Download PDF

Info

Publication number
KR20080039498A
KR20080039498A KR1020087006616A KR20087006616A KR20080039498A KR 20080039498 A KR20080039498 A KR 20080039498A KR 1020087006616 A KR1020087006616 A KR 1020087006616A KR 20087006616 A KR20087006616 A KR 20087006616A KR 20080039498 A KR20080039498 A KR 20080039498A
Authority
KR
South Korea
Prior art keywords
flow
gas
density
frequency
component
Prior art date
Application number
KR1020087006616A
Other languages
English (en)
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 마이크로 모우션, 인코포레이티드
Publication of KR20080039498A publication Critical patent/KR20080039498A/ko

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/76—Devices for measuring mass flow of a fluid or a fluent solid material
    • G01F1/78—Direct mass flowmeters
    • G01F1/80—Direct mass flowmeters operating by measuring pressure, force, momentum, or frequency of a fluid flow to which a rotational movement has been imparted
    • G01F1/84—Coriolis or gyroscopic mass flowmeters
    • G01F1/845—Coriolis or gyroscopic mass flowmeters arrangements of measuring means, e.g., of measuring conduits
    • G01F1/8468—Coriolis or gyroscopic mass flowmeters arrangements of measuring means, e.g., of measuring conduits vibrating measuring conduits
    • G01F1/8472—Coriolis or gyroscopic mass flowmeters arrangements of measuring means, e.g., of measuring conduits vibrating measuring conduits having curved measuring conduits, i.e. whereby the measuring conduits' curved center line lies within a plane
    • G01F1/8477—Coriolis or gyroscopic mass flowmeters arrangements of measuring means, e.g., of measuring conduits vibrating measuring conduits having curved measuring conduits, i.e. whereby the measuring conduits' curved center line lies within a plane with multiple measuring conduits
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/76—Devices for measuring mass flow of a fluid or a fluent solid material
    • G01F1/78—Direct mass flowmeters
    • G01F1/80—Direct mass flowmeters operating by measuring pressure, force, momentum, or frequency of a fluid flow to which a rotational movement has been imparted
    • G01F1/84—Coriolis or gyroscopic mass flowmeters
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/76—Devices for measuring mass flow of a fluid or a fluent solid material
    • G01F1/78—Direct mass flowmeters
    • G01F1/80—Direct mass flowmeters operating by measuring pressure, force, momentum, or frequency of a fluid flow to which a rotational movement has been imparted
    • G01F1/84—Coriolis or gyroscopic mass flowmeters
    • G01F1/8409—Coriolis or gyroscopic mass flowmeters constructional details
    • G01F1/8413—Coriolis or gyroscopic mass flowmeters constructional details means for influencing the flowmeter's motional or vibrational behaviour, e.g., conduit support or fixing means, or conduit attachments
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/76—Devices for measuring mass flow of a fluid or a fluent solid material
    • G01F1/78—Direct mass flowmeters
    • G01F1/80—Direct mass flowmeters operating by measuring pressure, force, momentum, or frequency of a fluid flow to which a rotational movement has been imparted
    • G01F1/84—Coriolis or gyroscopic mass flowmeters
    • G01F1/8409—Coriolis or gyroscopic mass flowmeters constructional details
    • G01F1/8436—Coriolis or gyroscopic mass flowmeters constructional details signal processing
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F15/00—Details 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/02—Compensating or correcting for variations in pressure, density or temperature
    • G01F15/022—Compensating or correcting for variations in pressure, density or temperature using electrical means
    • G01F15/024—Compensating or correcting for variations in pressure, density or temperature using electrical means involving digital counting

Landscapes

  • Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Measuring Volume Flow (AREA)
  • Physical Or Chemical Processes And Apparatus (AREA)

Abstract

본원발명의 일 실시예에 따라 유량계(5)의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 계측 전자장치(20)가 제공된다. 이러한 계측 전자장치는 다상 유동 물질에 대한 제1 센서 신호(210) 및 제2 센서 신호(211)를 수신하기 위한 인터페이스(201), 및 상기 인터페이스(201)와 통신하는 처리 시스템(203)을 포함한다. 상기 처리 시스템(203)은 상기 인터페이스(201)로부터 제1 센서 신호(210) 및 제2 센서 신호(211)를 수신하고, 제1 센서 신호(210)로부터 제1 90도 위상 편이(213)를 생성하고 제2 센서 신호(211)로부터 제2 90도 위상 편이(214)를 생성하고, 제1 90도 위상 편이(213) 또는 제2 90도 위상 편이(214) 중 하나를 사용하여 주파수(221)를 계산하고, 제1 90도 위상 편이(213) 및 제2 90도 위상 편이(214) 중 하나 이상을 사용하여 위상차(220)를 계산하고, 다상 유동 물질에 대한 질량 유량(223), 밀도(224), 또는 용적 유량(225) 중 하나 이상을 계산하도록 구성된다.
코리올리, 유량계, 유동 특성

Description

유량계의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 방법 및 계측 전자장치 {METER ELECTRONICS AND METHODS FOR PROCESSING SENSOR SIGNALS FOR A MULTI-PHASE FLOW MATERIAL IN A FLOWMETER}
본원발명은 유량계의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 방법 및 계측 전자장치에 관한 것이다.
1985년 1월 1일 J.E. Smith 등에게 허여된 미국 특허 번호 4,491,025 및 1982년 2월 11일 J.E. Smith 에게 허여된 Re. 31,450 호에 개시된 것과 같은 파이프라인을 통해 유동하는 물질의 질량 유동, 밀도, 용적 유동 및 기타 정보를 측정하기 위한 코리올리 질량 유량계를 사용하는 것이 공지되어 있다. 이러한 유량계는 다른 구성의 유동관을 하나 이상 구비한다. 각각의 도관 구성은 예를 들어, 단수 굽힘, 비틀림, 방사상 및 결합된 모드를 포함하는 자연 진동 세트를 갖는 것으로 간주할 수 있다. 통상의 코리올리 질량 유동 측정 장치에서, 도관 구성은 물질이 도관을 통해 유동함에 따라 하나 이상의 진동 모드로 여기되며, 도관의 운동은 도관을 따라 이격되어 있는 위치에서 측정된다.
물질이 채워진 시스템의 진동 모드는 유동관 내의 물질 및 유동관의 결합된 질량에 의해 부분적으로 규정된다. 물질은 유량계의 유입구 측 상에 연결된 파이 프라인으로부터 유량계로 흐른다. 이후 물질은 유동관 또는 유동관들로 향하게 되고 유량계를 떠나 배출구 측 상에 연결된 파이프라인으로 흐르게 된다.
드라이버는 유동관에 하중을 가한다. 하중은 유동관이 진동하도록 한다. 유량계를 통해 흐르는 물질이 없다면, 유동관을 따른 모든 지점은 동일한 위상으로 진동한다. 물질이 유동관을 통해 유동함에 따라, 코리올리 가속도가 유동관을 따른 각각의 위치들이 유동관을 따른 다른 위치에 대하여 다른 위상을 갖도록 한다. 유동관의 유입구 측에서의 위상은 드라이버에 뒤처지는 반면, 배출구 측에서의 위상은 드라이버를 앞선다. 유동관의 다른 위치에 센서들이 배치되어 다른 위치에서의 유동관의 운동을 나타내는 사인함수 신호를 생성한다. 두 개의 센서 신호 간의 위상차는 유동관 또는 유동관을 통해 유동하는 물질의 질량 유량에 비례한다. 종래의 방식 중 하나에서는 이산 푸리에 변환(Discrete Fourier Transform; DFT) 또는 고속 푸리에 변환(Fast Fourier Transform; FFT)이 센서 신호 사이의 위상차를 결정하기 위해 사용된다. 유동관 조립체의 진동 주파수 응답 및 위상차가 질량 유량을 얻기 위해 사용된다.
종래의 방식 중 하나에서는 예를 들어 진동 드라이버 시스템으로 송신되는 주파수를 사용함으로써 독립된 기준 신호가 픽오프(pickoff) 신호 주파수를 결정하기 위해 사용된다. 다른 접근 방식에서는, 구동 주파수에 대한 기준 신호 또는 픽오프의 주파수를 잠그기(lock) 위하여 위상 잠금식 루프(phase locked loop; PLL)가 사용된다. 다른 종래의 방식에서는 픽오프 센서에 의해 생성된 진동 주파수 응답이 노치 필터(notch filter)에 있는 주파수에 중심을 맞춤(centering)으로써 결 정될 수 있으며, 여기서 종래 기술의 유량계는 노치 필터의 노치를 픽오프 센서 주파수에 유지하려고 시도한다. 이러한 종래 기술은 정지된 상태에서는 매우 잘 작동하며, 여기서 유량계 내의 유동 물질은 균일하며 최종 픽오프 신호 주파수는 비교적 안정적이다. 그러나 액체 유동 물질 내에 공기 기포가 있거나 유동 물질이 액체 및 고체를 포함하는 2-상 유동(two-phase flow)과 같이, 유동 물질이 균일하지 않은 경우에는 종래기술의 위상 측정은 실패하게 된다. 이러한 상황에서, 종래 기술에서 결정된 주파수는 빠르게 요동한다. 빠르고 큰 주파수 변이의 상태 동안에는, 픽오프 신호가 필터 대역폭 외부로 이동하여 잘못된 위상 및 주파수 측정이 이루어질 수 있다. 이는 또한, 유량계가 빈 상태(empty condition)와 채워진 상태(full condition)에서 번갈아 가며 반복적으로 작동하는 엠티-풀-엠티 배칭(empty-full-empty batching) 에서도 문제가 된다. 또한, 센서의 주파수가 빠르게 이동하면, 복조(demodulation) 과정은 실제 또는 측정된 주파수를 따라갈 수 없어서 잘못된 주파수에서의 복조가 이루어지게 된다. 결정된 주파수가 부정확하거나 잘못되면, 후속적으로 얻어지는 밀도, 용적 유량 등의 값도 부정확하거나 잘못된다는 것을 이해해야 한다. 더욱이, 이러한 오차는 계속되는 유동 특성 측정에서 증가하게 된다.
종래 기술에서, 픽오프 신호는 노치 필터를 실행하기 위하여 디지털화되어 디지털 방식으로 다루어진다. 노치 필터는 단지 주파수의 좁은 대역만을 수용한다. 그러므로 목적 주파수가 변경되면 노치 필터는 주기 시간에 대해 목적 신호를 추적할 수 없게 된다. 통상적으로, 디지털 노치 필터 실행은 요동하는 목적 신호를 추 적하는데 1-2 초가 소요된다. 주파수를 결정하기 위해 종래 기술에서 요구되는 시간으로 인하여, 주파수 및 위상 결정이 오차를 포함할 뿐만 아니라 오차 측정은 오차 및/또는 2-상 유동이 실제적으로 일어나는 시간 간격(time span)을 초과하는 시간 간격을 포함하게 되는 결과가 발생한다. 이는 노치 필터 실행의 비교적 늦은 응답에 기인한다.
결과적으로 종래 기술의 유량계는 유량계 내의 유동 물질의 2-상 유동 과정 동안에 픽오프 센서 주파수를 정확하거나, 빠르거나, 만족스럽게 추적하거나 결정할 수 없다. 결과적으로, 종래 기술은 결정된 픽오프 주파수를 사용하여 위상차를 유도하므로, 위상 결정도 마찬가지로 느리고 오차를 갖기 쉽다. 따라서 주파수 측정에서의 어떠한 오차도 위상 결정에서 증가하게 된다. 결과적으로 주파수 결정 및 위상 결정에서 오차가 증가하여 질량 유량의 결정에 있어 오차를 증가시키게 된다. 또한, 결정된 주파수 값은 밀도 값을 결정하기 위해 사용되므로(밀도는 제곱 주파수 분의 1과 거의 동일하다), 주파수 결정에서의 오차는 밀도 결정에서도 반복되거나 증가하게 된다. 이는 또한 용적 유량의 경우에도 마찬가지인데, 용적 유량은 질량 유량을 밀도로 나눈 것과 동일하기 때문이다.
통상적인 오일 유정에서, 유정 산출물은 일반적으로 오일뿐만 아니라 유정 산출물 흐름 내에 천연 가스와 물의 다양한 양도 포함한다. 따라서 오일 유정 산출물은 통상적으로 다상 유체 흐름(multi-phase fluid stream)을 포함한다. 결과적으로, 유정 흐름은 보통 분리 장치로 공급된다. 분리 장치는 유정 흐름으로부터 하나 이상의 성분을 추출한다. 분리 장치는 2-상 분리기 또는 3-상 분리기를 포함 할 수 있다. 2-상 분리기는 통상적으로 다상 유체로부터 혼입 가스를 추출한다. 2-상 분리기의 산출물은 혼입 가스가 없는 2-상 유체 및 가스 산출물을 포함할 수 있다. 예를 들어, 유정 흐름이 오일, 물, 및 천연 가스를 포함하는 곳에서는, 2-상 분리기 액체 산출물이 2-상 오일 및 물 유동 흐름을 포함할 수 있다. 3-상 분리기는 혼입 가스를 분리해낼 수 있으며 오일로부터 물을 분리할 수도 있다.
그러나 분리기는 유동 성분을 완전히 분리하지는 않는다. 예를 들어, 추출된 가스는 작은 량의 잔여 액체를 포함할 수 있다. 분리기의 가스 출구가 통상적으로 액체의 출구보다 물리적으로 높게 위치하므로, 이는 통상적으로 캐리오버(carryover) 또는 액체 캐리오버라고 불린다. 또한, 3-상 분리기의 오일 산출물(또는 2-상 분리기의 오일 및 물 산출물)은 여전히 작은 양의 혼입 가스를 포함할 수 있다. 분리기의 액체/물 출구가 통상적으로 가스 출구보다 물리적으로 낮게 위치하므로, 이는 캐리언더(carryunder) 또는 가스 캐리언더라고 불린다.
종래 기술에 따른 분리기 산출물의 측정은 일정한 유형의 유량계로 각각의 분리기 산출물 흐름을 측정하는 것을 포함한다. 이러한 측정은 분리기의 가스 산출물 내에 무시할만한 액체 캐리오버가 있다는 것과 오일/물 산출물 내에 무시할만한 가스 캐리언더가 있다는 것을 가정한다. 종래 기술에서는 혼입 가스 및/또는 액체가 실질적으로 순간적으로 결정될 수 없기 때문에, 따라서 정확하게 결정될 수 없기 때문에 이러한 가정이 만들어진다. 그러나 캐리오버 및 캐리언더는 통상적으로 존재하며, 종래 기술의 유정 산출물 측정에 있어 심각한 부정확성을 초래한다. 이러한 부정확성은 가스 및 액체 종래 기술 측정을 조합함으로써 증가하게 된다.
유량계의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 방법 및 계측 전자장치를 제공함으로써 상술한 문제점들 및 기타의 문제점들이 해결되며 본 기술 분야에서의 기술 진보가 이루어진다.
본원발명의 일 실시예에 따라 유량계의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 계측 전자장치가 제공된다. 이러한 계측 전자장치는 다상 유동 물질에 대한 제1 센서 신호 및 제2 센서 신호를 수신하기 위한 인터페이스, 및 상기 인터페이스와 통신하는 처리 시스템을 포함한다. 상기 처리 시스템은 상기 인터페이스로부터 제1 센서 신호 및 제2 센서 신호를 수신하고, 제1 센서 신호로부터 제1 90도 위상 편이를 생성하고 제2 센서 신호로부터 제2 90도 위상 편이를 생성하고, 제1 90도 위상 편이 또는 제2 90도 위상 편이 중 하나를 사용하여 주파수를 계산하고, 제1 90도 위상 편이 및 제2 90도 위상 편이 중 하나 이상을 사용하여 위상차를 계산하고, 다상 유동 물질에 대한 질량 유량, 밀도, 또는 용적 유량 중 하나 또는 둘 이상을 계산하도록 구성된다.
본원발명의 일 실시예에 따라, 유량계 내의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 방법이 제공된다. 상기 방법은 다상 유동 물질에 대한 제1 센서 신호 및 제2 센서 신호를 수신하는 단계 및 제1 센서 신호로부터 제1 90도 위상 편이를 생성하고 제2 센서 신호로부터 제2 90도 위상 편이를 생성하는 단계를 포함한다. 상기 방법은 제1 90도 위상 편이 및 제2 90도 위상 편이 중 하나를 이용하여 주파수를 계산하는 단계를 더 포함한다. 상기 방법은 제1 90도 위상 편이 및 제2 90도 위상 편이 중 하나를 이용하여 위상차를 계산하는 단계를 더 포함한다. 상기 방법은 다상 유동 물질에 대한 질량 유량, 밀도, 또는 용적 유량 중 하나 또는 둘 이상을 계산하는 단계를 더 포함한다.
본원발명의 일 실시예에 따라, 유량계 내의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 방법이 제공된다. 상기 방법은 다상 유동 물질에 대한 제1 센서 신호 및 제2 센서 신호를 수신하는 단계 및 제1 센서 신호로부터 제1 90도 위상 편이를 생성하고 제2 센서 신호로부터 제2 90도 위상 편이를 생성하는 단계를 포함한다. 상기 방법은 제1 90도 위상 편이 및 제2 90도 위상 편이 중 하나를 이용하여 주파수를 계산하는 단계를 더 포함한다. 상기 방법은 제1 90도 위상 편이 및 제2 90도 위상 편이 중 하나를 이용하여 위상차를 계산하는 단계를 더 포함한다. 상기 방법은 다상 유동 물질에 대한 질량 유량, 밀도, 또는 용적 유량 중 하나 이상을 계산하는 단계 및 다상 유동 물질 내의 가스 캐리언더 또는 액체 캐리오버 중 하나 또는 둘 이상을 계산하는 단계를 더 포함한다.
계측 전자장치의 일 측면에서, 상기 인터페이스가 센서 신호를 디지털화하도록 구성된 디지타이저(digitizer)를 포함한다.
계측 전자장치의 다른 측면에서, 상기 처리 시스템이 다상 유동 물질 내의 액체 캐리오버 또는 가스 캐리언더 중 하나 이상을 결정하도록 추가적으로 구성된다.
계측 전자장치의 또 다른 측면에서, 상기 생성이 제1 90도 위상 편이 및 제2 90도 위상 편이를 생성하기 위하여 힐버트 변환을 이용하는 것을 포함한다.
계측 전자장치의 또 다른 측면에서, 상기 주파수 계산이 제1 센서 신호 및 제1 90도 위상 편이로부터 주파수를 계산하는 것을 포함한다.
계측 전자장치의 또 다른 측면에서, 상기 위상차 계산이 제1 센서 신호, 제1 90도 위상 편이, 및 제2 센서 신호로부터 위상차를 계산하는 것을 포함한다.
계측 전자장치의 또 다른 측면에서, 상기 위상차 계산이 제1 센서 신호, 제1 90도 위상 편이, 제2 센서 신호, 및 제2 90도 위상 편이로부터 위상차를 계산하는 것을 포함한다.
계측 전자장치의 또 다른 측면에서, 상기 처리 시스템이 주파수를 적어도 가스 주파수 성분 및 유체 주파수 성분으로 분리하고, 가스 주파수 성분 및 유체 주파수 성분 중 하나 이상 및 주파수 응답으로부터 액체 분율 또는 가스의 보이드 분율 중 하나 이상을 결정하고, 가스의 보이드 분율을 이용하여 가스 유동 성분의 가스상 밀도 또는 다상 유동 물질의 액체 유동 성분의 액체상 밀도 중 하나 이상을 결정하고, 그리고 액체 분율 또는 가스의 보이드 분율 중 하나 이상을 이용하여 액체 캐리오버 또는 다상 유동 물질의 가스 캐리언더 중 하나 이상을 결정하도록 추가적으로 구성된다.
계측 전자장치의 또 다른 측면에서, 상기 분리가 주파수 응답을 가스 주파수 성분 및 유체 주파수 성분 중 하나를 실질적으로 차단하는 하나 또는 둘 이상의 필터로 처리하는 것을 포함한다.
계측 전자장치의 또 다른 측면에서, 상기 분리가 가스 주파수 성분을 실질적으로 차단하고 유체 주파수 성분을 실질적으로 통과시키는 제1 필터로 주파수 응답을 필터링하는 것과, 유체 주파수 성분을 실질적으로 차단하고 가스 주파수 성분을 실질적으로 통과시키는 제2 필터로 주파수 응답을 필터링하는 것을 포함하고, 여기서 제1 필터가 유체 주파수 성분을 출력하고 제2 필터가 가스 주파수 성분을 출력한다.
계측 전자장치의 또 다른 측면에서, 액체 분율 또는 가스의 보이드 분율 중 하나 이상을 결정하는 것이, 주파수 응답으로부터 전체 밀도를 계산하는 것, 유체 주파수 성분으로부터 유체 성분 밀도를 계산하는 것, 가스 주파수 성분으로부터 가스 성분 밀도를 계산하는 것, 및 유체 성분 밀도에서 전체 밀도를 뺀 값을 유체 성분 밀도에서 가스 성분 밀도를 뺀 값으로 나눈 비율로서 가스의 보이드 분율을 계산하는 것을 포함한다.
계측 전자장치의 또 다른 측면에서, 상기 처리 시스템이 주파수 응답을 적어도 가스 주파수 성분 및 유체 주파수 성분으로 분리하고, 주파수 응답으로부터 전체 밀도를 결정하고, 가스 주파수 성분으로부터 가스 밀도를 결정하고, 가스 주파수 성분 및 유체 주파수 성분 중 하나 이상과 주파수 응답으로부터 가스의 보이드 분율을 결정하고, 그리고 가스의 보이드 분율에 가스의 밀도를 전체 밀도로 나눈 비율을 곱한 값으로부터 질량 분율을 결정하도록 추가적으로 구성된다.
계측 전자장치의 또 다른 측면에서, 상기 처리 시스템이 주파수 응답으로부터 유동 물질의 질량 유량을 결정하고 질량 분율 및 질량 유량을 이용하여 제1 유동 성분 질량과 제2 유동 성분 질량 중 하나 이상을 결정하도록 추가적으로 구성된다.
계측 전자장치의 또 다른 측면에서, 상기 처리 시스템이 제곱된 주파수 응답을 생성하도록 주파수 응답을 제곱하고, 실질적으로 순간적인 유동 흐름 밀도를 생성하도록 제곱된 주파수 응답의 역수를 취하고, 실질적으로 순간적인 유동 흐름 밀도를 가스 유동 물질의 가스 유동 분율을 나타내는 소정의 가스 밀도 및 액체 유동 분율을 나타내는 소정의 액체 밀도 중 하나 이상과 비교하고, 이러한 비교로부터 액체 유동 분율 또는 가스 유동 분율 중 하나 이상을 결정하도록 추가적으로 구성된다.
방법의 일 측면에서, 상기 방법은 다상 유동 물질 내의 액체 캐리오버 또는 가스 캐리언더 중 하나 이상을 결정하는 단계를 포함한다.
방법의 다른 측면에서, 상기 생성 단계가 제1 90도 위상 편이 및 제2 90도 위상 편이를 생성하기 위하여 힐버트 변환을 이용하는 단계를 포함한다.
방법의 또 다른 측면에서, 상기 주파수 계산 단계가 제1 센서 신호 및 제1 90도 위상 편이로부터 주파수를 계산하는 단계를 포함한다.
방법의 또 다른 측면에서, 상기 위상차 계산 단계가 제1 센서 신호, 제1 90도 위상 편이, 및 제2 센서 신호로부터 위상차를 계산하는 단계를 포함한다.
방법의 또 다른 측면에서, 상기 위상차 계산 단계가 제1 센서 신호, 제1 90도 위상 편이, 제2 센서 신호, 및 제2 90도 위상 편이로부터 위상차를 계산하는 단계를 포함한다.
방법의 또 다른 측면에서, 상기 방법은 주파수를 적어도 가스 주파수 성분 및 유체 주파수 성분으로 분리하는 단계, 가스 주파수 성분 및 유체 주파수 성분 중 하나 이상 및 주파수 응답으로부터 액체 분율 또는 가스의 보이드 분율 중 하나 이상을 결정하는 단계, 가스의 보이드 분율을 이용하여 가스 유동 성분의 가스상 밀도 또는 다상 유동 물질의 액체 유동 성분의 액체상 밀도 중 하나 이상을 결정하는 단계, 및 액체 분율 또는 가스의 보이드 분율 중 하나 이상을 이용하여 액체 캐리오버 또는 다상 유동 물질의 가스 캐리언더 중 하나 이상을 결정하는 단계를 더 포함한다.
방법의 또 다른 측면에서, 상기 분리 단계가 주파수 응답을 가스 주파수 성분 및 유체 주파수 성분 중 하나를 실질적으로 차단하는 하나 또는 둘 이상의 필터로 처리하는 단계를 포함한다.
방법의 또 다른 측면에서, 상기 분리 단계가 가스 주파수 성분을 실질적으로 차단하고 유체 주파수 성분을 실질적으로 통과시키는 제1 필터로 주파수 응답을 필터링하는 단계, 및 유체 주파수 성분을 실질적으로 차단하고 가스 주파수 성분을 실질적으로 통과시키는 제2 필터로 주파수 응답을 필터링하는 단계를 포함하고, 상기 제1 필터가 유체 주파수 성분을 출력하고 상기 제2 필터가 가스 주파수 성분(1412)을 출력한다.
방법의 또 다른 측면에서, 액체 분율 또는 가스의 보이드 분율 중 하나 이상을 결정하는 단계가, 주파수 응답으로부터 전체 밀도를 계산하는 단계, 유체 주파수 성분으로부터 유체 성분 밀도를 계산하는 단계, 가스 주파수 성분으로부터 가스 성분 밀도를 계산하는 단계, 및 유체 성분 밀도에서 전체 밀도를 뺀 값을 유체 성분 밀도에서 가스 성분 밀도를 뺀 값으로 나눈 비율로서 가스의 보이드 분율을 계산하는 단계를 포함한다.
방법의 또 다른 측면에서, 상기 방법이 주파수 응답을 적어도 가스 주파수 성분 및 유체 주파수 성분으로 분리하는 단계, 주파수 응답으로부터 전체 밀도를 결정하는 단계, 가스 주파수 성분으로부터 가스 밀도를 결정하는 단계, 가스 주파수 성분 및 유체 주파수 성분 중 하나 이상과 주파수 응답으로부터 가스의 보이드 분율을 결정하는 단계, 및 가스의 보이드 분율에 가스의 밀도를 전체 밀도로 나눈 비율을 곱한 값으로부터 질량 분율을 결정하는 단계를 더 포함한다.
방법의 또 다른 측면에서, 상기 방법이 주파수 응답으로부터 유동 물질의 질량 유량을 결정하는 단계 및 질량 분율 및 질량 유량을 이용하여 제1 유동 성분 질량과 제2 유동 성분 질량 중 하나 이상을 결정하는 단계를 더 포함한다.
방법의 또 다른 측면에서, 상기 방법이 제곱된 주파수 응답을 생성하도록 주파수 응답을 제곱하는 단계, 실질적으로 순간적인 유동 흐름 밀도를 생성하도록 제곱된 주파수 응답의 역수를 취하는 단계, 실질적으로 순간적인 유동 흐름 밀도를 가스 유동 물질의 가스 유동 분율을 나타내는 소정의 가스 밀도 및 액체 유동 분율을 나타내는 소정의 액체 밀도 중 하나 이상과 비교하는 단계, 및 상기 비교로부터 액체 유동 분율 또는 가스 유동 분율 중 하나 이상을 결정하는 단계를 더 포함한다.
전체 도면에 있어서 동일한 도면 부호는 동일한 부재를 나타낸다.
도 1은 본원발명의 일 실시예의 코리올리 유량계를 도시한다.
도 2는 본원발명의 일 실시예에 따른 계측 전자장치를 도시한다.
도 3은 본원발명의 일 실시예에 따른, 유량계의 센서 신호를 처리하는 방법의 순서도이다.
도 4는 본원발명의 일 실시예에 따른 계측 전자장치를 도시한다.
도 5는 본원발명의 일 실시예에 따른, 유량계의 제1 및 제2 센서 신호를 처리하는 방법의 순서도를 도시한다.
도 6은 본원발명의 일 실시예에 따른 처리 시스템의 일부의 블록 다이어그램을 도시한다.
도 7은 본원발명의 일 실시예에 따른 힐버트 변환 블록을 상세하게 도시한다.
도 8 및 9는 본원발명의 일 실시예에 따른 분석 블록의 두 개의 독립 분기의 블록 다이어그램을 도시한다.
도 10은 정상 상태에 있는 유량계의 픽-오프 센서 신호의 파워 스펙트럼 밀도 그래프를 도시한다.
도 11은 단일 위상 편이 실시예에 따른 힐버트 변환 블록을 도시한다.
도 12는 단일 위상 편이 실시예에 대한 분석 블록을 도시한다.
도 13은 종래 기술과 비교하여 본원발명의 센서 처리를 도시하며, 여기서 각각의 시간차(Δt) 값이 비교된다.
도 14는 본원발명의 다른 실시예에 따른 계측 전자장치를 도시한다.
도 15는 공기, 유체, 및 공기/유체 혼합물(즉, 혼입 공기를 포함하는 유체) 에 대한 유량계 주파수 응답의 그래프를 도시한다.
도 16은 본원발명의 일 실시예에 따른 유량계를 통해 유동하는 유동 물질 내의 가스의 보이드 분율(void fraction)을 결정하기 위한 방법의 순서도를 도시한다.
도 17은 본원발명의 일 실시예에 따른 유량계를 통해 유동하는 유동 물질 내의 가스의 보이드 분율을 결정하기 위한 방법의 순서도를 도시한다.
도 18은 본원발명의 일 실시예에 따라 가스 주파수 성분 및 유체 주파수 성분의 분리를 위해 사용될 수 있는 로패스(low-pass) 및 하이패스(high-pass) 필터 응답을 도시하는 주파수 그래프이다.
도 19는 본원발명의 일 실시예에 따라 유량계를 통해 유동하는 유동 물질 내의 가스의 보이드 분율을 결정하기 위한 방법의 순서도이다.
도 20은 노치 필터 주파수 응답의 그래프이다.
도 21은 본원발명의 일 실시예에 따라 유량계를 통해 유동하는 유동 물질의 유동 성분의 질량 분율을 결정하기 위한 방법의 순서도를 도시한다.
도 22는 본원발명의 일 실시예에 따라 유량계를 통해 유동하는 유동 물질의 유동 성분의 질량 분율을 결정하기 위한 방법의 순서도를 도시한다.
도 23은 본원발명의 일 실시예에 따라 유량계를 통해 유동하는 유동 물질의 유동 성분의 질량 분율을 결정하기 위한 방법의 순서도를 도시한다.
도 24는 본원발명의 일 실시예에 따라 유량계를 통해 유동하는 가스 유동 물질의 액체 유동 분율을 결정하기 위한 방법의 순서도를 도시한다.
도 25는 천연 가스 밀도 대 글리콜 백분율(즉, 액체 유동 분율)의 그래프이다.
도 26은 본원발명의 일 실시예에 따라 유량계를 통해 유동하는 가스 유동 물질의 액체 유동 분율을 결정하기 위한 방법의 순서도를 도시한다.
도 27은 본원발명의 일 실시예에 따라 유량계의 센서 신호를 처리하기 위한 방법의 순서도를 도시한다.
도 1-27 및 및 이하의 설명은 본원발명이 속하는 기술분야에서 통상의 지식을 가진 자에게 본원발명의 최상의 모드를 사용하고 제작하는 방법에 대하여 설명하기 위한 구체적인 실시예를 설명한다. 본원발명의 원리를 설명하기 위하여, 몇몇의 통상적인 측면은 간략하게 설명되거나 생략되었다. 본원발명이 속하는 기술분야에서 통상의 지식을 가진 자는 이러한 실시예의 수정이 본원발명의 기술분야에 속한다는 것을 알 수 있다. 본원발명이 속하는 기술분야에서 통상의 지식을 가진 자는 이하에 설명된 특징들이 다양한 방식으로 결합되어 본원발명의 다양한 수정을 형성하게 된다는 것을 알 수 있다. 결과적으로, 본원발명은 이하에 기술된 특정 실시예에 의하여 제한되지 않으며, 단지 청구범위 및 그 균등물에 의해서만 제한된다.
도 1은 유량계 조립체(10) 및 계측 전자장치(20)를 포함하는 코리올리 유량계(5)를 도시한다. 유량계 조립체(10)는 프로세스 물질의 밀도 및 질량 유량에 응답한다. 계측 전자장치(20)는 리드(lead; 100)를 통해 유량계 조립체(10)에 연결 되어 경로(26)로 밀도, 질량 유량(mass flow rate), 및 온도에 관한 정보를 제공할 뿐만 아니라 본원발명과 관련없는 기타의 정보도 제공한다. 비록 본원발명이 속하는 기술분야에서 통상의 지식을 가진 자에게 본원발명이 코리올리 질량 유량계에 의해 제공되는 추가적인 측정 능력이 없이도 진동의 관 농도계(densitometer)로서 실시될 수 있다는 것이 자명하기는 하지만, 코리올리 유량계 구조체에 관해 기술한다.
유량계 조립체(10)는 한 쌍의 매니폴드(150, 150'), 플랜지 넥(flange neck)을 갖는 플랜지(103, 103'), 한 쌍의 평행한 유동관(130, 130'), 드라이버 메커니즘(180), 온도 센서(190), 및 한 쌍의 속도 센서(170L, 170R)를 포함한다. 유동관(130, 130')은 두 개의 실질적으로 곧은 유입구 다리(131, 131') 및 배출구 다리(134, 134')를 갖는데, 이는 유동관 장착 블록(120, 120')에서 서로를 향해 모인다. 유동관(130, 130')은 그 길이를 따라 두 개의 대칭 지점에서 굽혀지며 그 길이에 걸쳐 실질적으로 평행하다. 브레이스 바(brace bar; 140, 140')는 축(W, W')을 형성하는데 사용되며, 이 축 주위로 각각의 유동관이 진동한다.
유동관(130, 130')의 측면 다리(131, 131') 및 (134, 134')는 유동과 장착 블록(120, 120')에 고정적으로 부착되며, 이들 블록은, 계속해서, 매니폴드(150, 150')에 고정적으로 부착된다. 이는 코리올리 유량계 조립체(10)를 통하는 연속적인 밀폐 물질 경로를 제공한다.
구멍(102, 102')을 갖는 플랜지(103, 103')가 유입 단부(104) 및 배출 단부(104')를 통해 측정되는 프로세스 물질을 운반하는 프로세스 라인(도시되지 않 음)에 연결되면, 물질은 플랜지(103)의 오리피스(101)를 통해 유량계의 단부(104)로 유입되어 매니폴드(150)를 통해 표면(121)을 갖는 유동관 장착 블록(120)으로 안내된다. 매니폴드(150) 내에서 물질은 분리되어 유동관(130, 130')을 통해 전달된다. 유동관(130, 130')을 여기시키면, 프로세스 물질은 매니폴드(150')에서 단일한 스트림으로 재결합되고, 이후 볼트 구멍(102')을 갖는 플랜지(103')에 의하여 프로세스 라인(도시되지 않음)에 연결된 배출 단부(104')로 전달된다.
유동관(130, 130')은 각각, 굽힘축(W--W 및 W'--W')에 대해 실질적으로 동일한 질량 분포, 관성 모멘트, 및 영률(Young's modulus)을 갖도록, 유동관 장착 블록(120, 120')에 대해 선택되어 적절하게 장착된다. 이러한 굽힘축은 브레이스 바(140, 140')를 통해 지나간다. 유동관의 영률이 온도에 따라 변하므로, 또한 이러한 변화가 밀도 및 유동의 계산에 영향을 미치므로, 저항성 온도 검측기(resistive temperature detector; RTD)(190)가 유동관(130')에 장착되어 유동관의 온도를 계속적으로 측정한다. 유동관의 온도 및 주어진 통과 전류에 대해 RTD 에 걸쳐 나타나는 전압은 유동관을 통과하는 물질의 온도에 의해 지배된다. RTD에 걸쳐 나타나는 온도 종속성 전압은 계측 전자장치(20)에서 공지의 방법으로 사용되어 유동관의 임의의 온도 변화로 기인하는 유동관(130, 130')의 탄성계수의 변화를 보상한다. RTD 는 리드(195)에 의하여 계측 전자장치(20)에 연결된다.
양 유동관(130, 130')은 드라이버(180)에 의하여 그 각각의 굽힘축(W 및 W')에 대하여 반대 방향으로 소위 유량계의 제1 이상(out-of-phase) 굽힘 모드에서 구동된다. 드라이버 메커니즘(180)은 유동관(130')에 장착되는 자석 및 유동관(130) 에 장착되는 대향 코일과 같은 많은 공지의 장치 중 하나를 포함할 수 있으며, 양 도관을 진동시키기 위하여 교류가 이를 통과한다. 계측 전자장치(20)에 의하여 적당한 구동신호가 리드(185)를 통해 드라이버 메커니즘(180)으로 인가된다.
계측 전자장치(20)는 리드(195) 상에서 RTD 신호를 수신하며, 각각 리드(165L, 165R) 상에 나타나는 좌우 속도 신호를 수신한다. 계측 전자장치(20)는 리드(185) 상에 나타나는 구동 신호를 생성하여 부재(180)를 구동시키고 유동관(130, 130')을 진동시킨다. 계측 전자장치(20)는 유량계 조립체(10)를 통과하는 물질의 밀도 및 질량 유량을 계산하기 위하여 RTD 신호와 좌우 속도 신호를 처리한다. 이러한 정보는 다른 정보와 함께 계측 전자장치(20)에 의하여 경로(26)를 통해 이용 수단(29)에 적용된다.
도 2는 본원발명의 일 실시예에 따른 계측 전자장치(20)를 도시한다. 계측 전자장치(20)는 인터페이스(201) 및 처리 시스템(203)을 포함할 수 있다. 계측 전자장치(20)는 픽오프/속도 센서 신호와 같이 유량계 조립체(10)로부터의 제1 및 제2 센서 신호를 수신한다. 계측 전자장치(20)는 코리올리 유량계로서의 작동을 포함하여, 농도계로서 작동하거나 질량 유량계로서 작동할 수 있다. 계측 전자장치(20)는 유량계 조립체(10)를 통과하여 유동하는 유동 물질의 유동 특성을 얻기 위하여 제1 및 제2 센서 신호를 처리한다. 예를 들어, 계측 전자장치(20)는 센서 신호로부터 위상차, 주파수, 시간차(Δt), 밀도, 질량 유량, 용적 유량 중 하나 이상을 결정할 수 있다. 또한, 다른 유동 특성도 본원발명에 따라 결정될 수 있다. 이러한 결정에 대해서는 이하에서 기술된다.
위상차 결정 및 주파수 결정은 종래 기술에서의 그러한 결정보다 더 빠르고 정확하며 신뢰성이 있다. 일 실시예에서, 위상차 결정 및 주파수 결정은, 주파수 기준 신호에 대한 필요성 없이, 단 하나의 센서 신호의 위상 편이(phase shift)로부터 직접적으로 유도된다. 이는 유동 특성을 계산하기 위해 요구되는 처리 시간을 유리하게 감소시킨다. 다른 실시예에서, 주파수는 단지 하나의 위상 편이 신호로부터 유도되는 반면, 위상차는 양 센서 신호의 위상 편이로부터 유도된다. 이는 양 유동 특성의 정확성을 증가시키며, 이들 모두는 종래기술에서보다 더욱 빠르게 결정될 수 있다.
종래 기술의 주파수 결정 방법은 실행하는 데 있어 통상적으로 1~2 초가 소요된다. 대조적으로, 본원발명에 따른 주파수 결정은 50 미리 세컨드(ms) 정도의 적은 시간에 실행될 수 있다. 처리 시스템의 유형 및 구성, 진동 응답의 샘플링 율(sampling rate), 필터 크기, 데시메이션 비율(decimation rate) 등등에 따라 더욱 빠른 주파수 결정도 가능하다. 50 ms 의 주파수 결정 속도에서 본원발명에 따른 계측 전자장치(20)는 종래 기술보다 약 40배 빠를 수 있다.
인터페이스(201)는 도 1의 리드(100)를 통해서 속도 센서(170L, 170R) 중 하나로부터 센서 신호를 수신한다. 인터페이스(201)는 포매팅(formatting), 증폭(amplification), 버퍼링(buffering) 등의 방식과 같이, 필요하거나 요구되는 신호 조정(signal conditioning)을 수행할 수 있다. 대안적으로, 신호 조정 중 일부 또는 전부는 처리 시스템(203)에서 실행될 수 있다.
또한, 인터페이스(201)는 계측 전자장치(20)와 외부 장치 사이의 통신을 가 능하게 할 수 있다. 인터페이스(201)는 전자, 광학, 또는 무선 통신의 어떠한 방식도 가능할 수 있다.
일 실시예에서 인터페이스(201)는 디지타이저(digitizer; 202)와 연결되며, 여기서 센서 신호는 아날로그 센서 신호를 포함한다. 디지타이저(202)는 아날로그 센서 신호를 표본추출하고 디지털화하여 디지털 센서 신호를 생성한다. 디지타이저(202)는 또한 요구되는 어떠한 데시메이션도 실행할 수 있으며, 여기서 디지털 센서 신호는 요구되는 신호 처리의 양을 줄이고 처리 시간을 줄이기 위하여 데시메이트(decimate)된다. 데시메이션은 이하에서 보다 자세히 설명된다.
처리 시스템(203)은 계측 전자장치(20)의 가동을 수행하며 유량계 조립체(10)로부터 유동 측정치를 처리한다. 처리 시스템(203)은 하나 이상의 처리 과정(processing routine)을 실행하며 이로써 하나 이상의 유동 특성을 생성하기 위하여 유동 측정치를 처리한다.
처리 시스템(203)은 일반적 용도의 컴퓨터, 마이크로 프로세스 시스템, 논리 회로, 또는 기타의 일반 용도나 주문형 처리 장치를 포함할 수 있다. 처리 시스템(203)은 다수의 처리 장치 사이에 분배될 수 있다. 처리 시스템(202)은 저장 시스템(204)과 같이, 어떠한 형태의 일체형 또는 독립형 전자 저장 매체도 포함할 수 있다.
처리 시스템(203)은 센서 신호(210)로부터 하나 이상의 유동 특성을 결정하기 위하여 센서 신호(210)를 처리한다. 하나 이상의 유동 특성으로는 예를 들어, 유동 물질에 대한 위상차, 주파수, 시간차(Δt), 질량 유량, 및/또는 밀도가 있다.
도시된 실시예에서, 처리 시스템(203)은 하나의 센서 신호 위상 편이(213) 및 두 개의 센서 신호(210, 211)로부터 유동 특성을 결정한다. 처리 시스템(203)은 하나의 위상 편이(213) 및 두 개의 센서 신호(210, 211)로부터 적어도 위상 차 및 주파수를 결정할 수 있다. 결과적으로, (상류 및 하류 픽오프 신호 중 하나와 같은) 제1 또는 제2 위상 편이된 센서 신호가 본원발명에 따라 처리 시스템(203)에 의하여 처리되어 유동 물질에 대한 위상차, 주파수, 시간차(Δt), 및/또는 질량 유량을 결정할 수 있다.
저장 시스템(204)은 유량계 파라미터 및 데이터, 소프트웨어 루틴(software routines), 상수 값, 및 변수 값을 저장할 수 있다. 일 실시예에서, 저장 시스템(204)은 처리 시스템(203)에 의하여 실행되는 처리 과정을 포함한다. 저장 시스템(204)은 일 실시예에서 위상 편이 루틴(212), 위상차 루틴(215), 주파수 루틴(216), 시간차(Δt) 루틴(217), 및 유동 특성 루틴(218)을 저장한다.
일 실시예에서, 저장 시스템(204)은 코리올리 유량계(5)와 같은 유량계를 작동시키기 위해 사용되는 변수들을 저장한다. 저장 시스템(204)은 일 실시예에서 제1 센서 신호(210) 및 제2 센서 신호(211)와 같은 변수들을 저장하는데, 이는 속도/픽오프 센서(170L, 170R)로부터 수신된다. 또한, 저장 시스템(204)은 유동 특성을 결정하기 위해 생성되는 90도 위상 편이(213)를 저장할 수 있다.
일 실시예에서, 저장 시스템(204)은 유동 측정으로부터 얻어지는 하나 이상의 유동 특성을 저장한다. 저장 시스템(204)은 일 실시예에서 위상차(220), 주파수(221), 시간차(Δt)(222), 질량 유량(223), 밀도(224), 및 용적 유량(225)과 같 은 유동 특성을 저장하며, 이들은 모두 센서 신호(210)로부터 결정된다.
위상 편이 루틴(212)은 입력 신호에 대하여, 즉 센서 신호(210)에 대하여 90도 위상 편이를 실행한다. 위상 편이 루틴(212)은 일 실시예에서, 힐버트 변환(Hilbert transform)(이하에서 설명됨)을 실행한다.
위상차 루틴(215)은 하나의 90도 위상 편이(213)를 사용하여 위상차를 결정한다. 위상차를 결정하기 위하여 추가적인 정보가 사용될 수도 있다. 일 실시예에서 위상차는 제1 센서 신호(210), 제2 센서 신호(211), 및 90도 위상 편이(213)로부터 계산될 수 있다. 결정된 위상차는 저장 시스템(204)의 위상차(220)에 저장될 수 있다. 위상차는, 90도 위상 편이(213)로부터 결정될 때, 종래 기술에서보다 더욱 빠르게 계산되어 얻어질 수 있다. 이는 높은 유량 속도를 갖는 유량계 장치 또는 다중 위상 유동이 발생하는 경우에 있어 중요한 차이를 제공할 수 있다. 또한, 위상차는 센서 신호(210) 또는 센서 신호(211) 중 하나의 주파수와 무관하게 결정될 수 있다. 더욱이, 위상차가 주파수와 무관하게 결정되기 때문에, 위상차의 에러 성분은 주파수 결정의 에러 성분을 포함하지 않는다. 즉, 위상차 측정에는 복합된 에러가 없다. 결과적으로, 위상차 에러는 종래 기술의 위상차에 비해 감소된다.
주파수 루틴(216)은 90도 위상 편이(213)로부터 (제1 센서 신호(210) 또는 제2 위상 신호(211)에 의해 나타내지는 것과 같은) 주파수를 결정한다. 결정된 주파수는 저장 시스템(204)의 주파수(221) 내에 저장될 수 있다. 주파수는, 90도 위상 편이(213)로부터 결정될 때, 종래 기술에서보다 더욱 빠르게 계산되어 얻어질 수 있다. 이는 높은 유량 속도를 갖는 유량계 장치 또는 다중 위상 유동이 발생하는 경우에 있어 중요한 차이를 제공할 수 있다.
시간차(Δt) 루틴(217)은 제1 센서 신호(210) 및 제2 센서 신호(211) 사이의 시간차(Δt)를 결정한다. 시간차(Δt)는 저장 시스템(204)의 시간차(Δt)(222)에 저장될 수 있다. 시간차(Δt)는 결정된 주파수에 의해 나눠진 결정된 위상을 실질적으로 포함하며, 따라서 질량 유량을 결정하기 위해 사용된다.
유동 특성 루틴(218)은 하나 이상의 유동 특성을 결정할 수 있다. 유동 특성 루틴(218)은 이러한 추가적인 유동 특성을 얻기 위하여, 예를 들어, 결정된 주파수(221) 및 결정된 위상차(220)를 사용할 수 있다. 이러한 결정을 위하여, 예를 들어 질량 유량이나 밀도와 같은 추가적인 정보가 사용될 수 있다는 것을 이해해야 한다. 유동 특성 루틴(218)은 시간차(Δt)(222)로부터, 따라서 위상차(220) 및 주파수(221)로부터 질량유량을 결정할 수 있다. 질량 유량을 결정하기 위한 공식은, 본 명세서에 참조로서 병합되는 티트로(Titlow) 등의 미국 특허 번호 5,027,662호에 개시되어 있다. 질량 유량은 유량계 조립체(10)의 유동 물질의 질량 유동과 관련된다. 마찬가지로, 유동 특성 루틴(218)은 밀도(224) 및/또는 용적 유량(225)도 결정할 수 있다. 결정된 질량 유량, 밀도, 및 용적 유량은 저장 시스템(204)의 질량 유량(223), 밀도(224), 및 용적 유량(225)에 각각 저장될 수 있다. 또한, 유동 특성은 계측 전자장치(20)에 의하여 외부 장치로 전달될 수 있다.
도 3은 본원발명의 실시예에 따른 유량계의 센서 신호를 처리하는 방법의 순서도(300)이다. 단계(301)에서, 제1 및 제2 센서 신호가 수신된다. 제1 센서 신 호는 상류 또는 하류 픽오프 센서 신호를 포함할 수 있다.
단계(302)에서, 센서 신호는 조정될 수 있다. 일 실시예에서, 이러한 조정은 의도치 않은 신호 및 노이즈를 제거하기 위한 필터링을 포함할 수 있다. 일 실시예에서, 필터링은 유량계의 예상된 기본 주파수 부근에 센터링된 대역 통과 필터링(band-pass filtering)을 포함할 수 있다. 또한, 증폭이나 버퍼링 등과 같은 다른 조정 작업이 실행될 수도 있다. 센서 신호가 아날로그 신호를 포함하면, 상기 단계는 디지털 센서 신호를 생성하기 위해 실행되는 샘플링, 디지털화(digitization), 및 데시메이션의 방법을 더 포함할 수도 있다.
단계(303)에서, 단일의 90도 위상 편이가 생성된다. 90도 위상 편이는 센서 신호의 90도 위상 편이를 포함한다. 90도 위상 편이는 어떠한 방식의 위상 편이 메커니즘이나 작동에 의해서도 실행될 수 있다. 일 실시예에서, 90도 위상 편이는 디지털 센서 신호에 대해 작동하는 힐버트 변환을 사용하여 실행된다.
단계(304)에서, 위상차가 단일의 90도 위상 편이를 사용하여 계산된다. 위상차를 계산하기 위하여 추가 정보가 사용될 수도 있다. 일 실시예에서, 위상차가 제1 센서 신호, 제2 센서 신호, 및 단일의 90도 위상 편이로부터 결정된다. 위상차는 응답 신호, 즉 픽오프 신호에 있는 위상차를 포함하는데, 이는 진동하는 유량계 조립체(10)에서의 코리올리 효과에 기인하여 나타난다.
최종 위상차는 계산에 어떠한 주파수 값을 필요로 하지 않고 결정된다. 최종 위상차는 주파수를 사용하여 계산되는 위상차보다 더 빠르게 얻어질 수 있다. 최종 위상차는 주파수를 사용하여 계산된 위상차보다 더 큰 정확성을 갖는다.
단계(305)에서, 주파수가 계산된다. 본원발명에 따른 주파수는 90도 위상 편이로부터 바람직하게 계산된다. 일 실시예에서 주파수는 90도 위상 편이 및 90도 위상 편이가 유도되는 해당 센서 신호를 사용한다. 주파수는 제1 센서 신호 및 제2 센서 신호 중 하나의 진동 응답 주파수이다(이러한 두 센서의 주파수는 작동에 있어서 실질적으로 동일하다). 주파수는 드라이버(180)에 의해 생성된 진동에 대한 유동관 또는 유동관들의 진동 주파수 응답을 포함한다.
이렇게 유도되는 주파수는 어떠한 독립적인 주파수 기준 신호에 대한 필요 없이 얻어진다. 주파수는 종래 기술에서보다 더 빠른 작동에서 단일의 90도 위상 편이로부터 얻어진다. 최종 주파수는 종래기술에서 계산된 주파수보다 더 큰 정확성을 갖는다.
단계(306)에서, 유동 물질의 질량 유량이 계산된다. 질량 유량은 단계(304) 및 단계(305)에서 계산된 최종 주파수 및 최종 위상차로부터 계산된다. 또한, 질량 유량 계산은 주파수 및 위상차로부터 시간차(Δt)를 계산할 수 있으며, 시간차(Δt)는 최종적으로 질량 유량을 계산하는데 사용된다.
단계(307)에서, 밀도가 선택적으로 결정될 수 있다. 밀도는 유동 특성의 하나로서 결정될 수 있으며, 예를 들어 주파수로부터 결정될 수 있다.
단계(308)에서, 용적 유량이 선택적으로 결정될 수 있다. 용적 유량은 유동 특성의 하나로서 결정될 수 있으며, 예를 들어 질량 유량 및 밀도로부터 결정될 수 있다.
도 4는 본원발명의 일 실시예에 따른 계측 전자장치(20)를 도시한다. 도 2 와 공통되는 부재는 동일한 도면부호가 사용되었다.
이 실시예의 계측 전자장치(20)는 제1 센서 신호(210) 및 제2 센서 신호(211)를 포함한다. 처리 시스템(203)은 이러한 신호들로부터 하나 이상의 유동 특성을 결정하기 위하여 제1 및 제2 (디지털) 센서 신호(210, 211)를 처리한다. 앞서 기술한 바와 같이, 하나 이상의 유동 특성은 유동 물질에 대한 위상차, 주파수, 시간차(Δt), 질량 유량, 밀도, 및/또는 용적 유량을 포함할 수 있다.
도시된 실시예에서, 처리 시스템(203)은, 외부 주파수 측정에 대한 필요성 및 외부 주파수 기준 신호에 대한 어떠한 필요성 없이, 단지 두 개의 센서 신호(210, 211)로부터 유동 특성을 결정한다. 처리 시스템(203)은 두 개의 센서 신호(210, 211)로부터 적어도 위상차 및 주파수를 결정할 수 있다.
전술한 바와 같이, 저장 시스템(204)은 위상 편이 루틴(212), 위상차 루틴(215), 주파수 루틴(216), 시간차(Δt) 루틴(217), 및 유동 특성 루틴(218)을 저장한다. 저장 시스템(204)은 제1 센서 신호(210) 및 제2 센서 신호(211)를 저장한다. 저장 시스템(204)은 또한 유동 특성을 결정하기 위하여 센서 신호들로부터 발생되는 제1 90도 위상 편이(213) 및 제2 90도 위상 편이도 저장한다. 전술한 바와 같이, 저장 시스템(204)은 위상차(220), 주파수(221), 시간차(Δt)(222), 질량 유량(223), 밀도(224), 및 용적 유량(225)을 저장한다.
위상 편이 루틴(212)은 제1 센서 신호(210) 및 제2 센서 신호(211)에 대한 것을 포함하여 입력 신호에 대한 90도 위상 편이를 실행한다. 위상 편이 루틴(212)은 일 실시예에서 힐버트 변환(아래에 설명됨)을 실행한다.
위상차 루틴(215)은 제1 90도 위상 편이(213) 및 제2 90도 위상 편이(214)를 사용하여 위상차를 결정한다. 위상차를 계산하기 위하여 추가 정보가 사용될 수도 있다. 일 실시예에서 위상차는 제1 센서 신호(210), 제2 센서 신호(211), 제1 90도 위상 편이(213), 및 제2 90도 위상 편이(214)로부터 계산된다. 결정된 위상차는, 전술한 바와 같이, 저장 시스템(204)의 위상차(220)에 저장될 수 있다. 위상차는, 제1 및 제2 90도 위상 편이를 사용하여 결정될 때, 종래기술에서보다 더 빠르게 계산되어 얻어질 수 있다. 이는 높은 유량을 갖는 유량계 장치나 다중 위상 유동이 발생하는 경우에 중요한 차이를 제공할 수 있다. 또한, 위상차는 센서 신호(210, 211)의 주파수와 무관하게 결정될 수 있다. 더욱이, 위상차는 주파수와 무관하게 결정되므로, 위상차의 에러 성분은 주파수 결정의 에러 성분의 영향을 받지 않는다. 즉, 위상차 측정에는 복합된 에러가 없다. 결과적으로, 위상차 에러는 종래기술의 위상차에 비해 줄어들게 된다.
주파수 루틴(216)은 제1 90도 위상 편이(213) 및 제2 90도 위상 편이(214)로부터 (제1 센서 신호(210) 또는 제2 센서 신호(211)에 의해 나타내지는 것과 같은) 주파수를 결정한다. 결정된 주파수는 전술한 바와 같이 저장 시스템(204)의 주파수(221)에 저장될 수 있다. 주파수는, 제1 및 제2 90도 위상 편이로부터 결정될 때, 종래기술에서보다 더 빠르게 계산되어 얻어질 수 있다. 이는 높은 유량을 갖는 유량계 장치나 다중 위상 유동이 발생하는 경우에 중요한 차이를 제공할 수 있다.
시간차(Δt) 루틴(217)은 제1 센서 신호(210) 및 제2 센서 신호(211) 사이의 시간차(Δt)를 결정한다. 시간차(Δt)는 전술한 바와 같이 저장 시스템(204)의 시간차(Δt)(222)에 저장될 수 있다. 시간차(Δt)는 결정된 주파수에 의해 나눠진 결정된 위상을 실질적으로 포함하며, 따라서 질량 유량을 결정하기 위해 사용된다.
유동 특성 루틴(218)은 전술한 바와 같이, 질량 유량, 밀도, 및/또는 용적 유량 중 하나 이상을 결정할 수 있다.
도 5는 본원발명의 실시예에 따른 유량계의 제1 및 제2 센서 신호를 처리하는 방법의 순서도(500)이다. 단계(501)에서, 제1 및 제2 센서 신호가 수신된다. 일 실시예에서, 제1 센서 신호는 상류 또는 하류 픽오프 센서 신호를 포함할 수 있다.
단계(502)에서 제2 센서 신호가 수신된다. 일 실시예에서, 제2 센서 신호는 하류 또는 상류 픽오프 센서 신호(즉, 제1 센서 신호의 반대)를 포함할 수 있다.
단계(503)에서, 센서 신호는 조정될 수 있다. 일 실시예에서, 이러한 조정은 의도치 않은 신호 및 노이즈를 제거하기 위한 필터링을 포함할 수 있다. 일 실시예에서, 필터링은 전술한 바와 같이 대역 통과 필터링을 포함할 수 있다. 또한, 증폭이나 버퍼링 등과 같은 다른 조정 작업이 실행될 수도 있다. 센서 신호가 아날로그 신호를 포함하면, 상기 단계는 디지털 센서 신호를 생성하기 위해 실행되는 샘플링, 디지털화, 및 데시메이션의 방법을 더 포함할 수도 있다.
단계(504)에서, 제1 90도 위상 편이가 생성된다. 제1 90도 위상 편이는 제1 센서 신호의 90도 위상 편이를 포함한다. 90도 위상 편이는 어떠한 방식의 메커니즘이나 작동에 의해서도 실행될 수 있다. 일 실시예에서, 90도 위상 편이는 디지 털 센서 신호에 대해 작동하는 힐버트 변환을 사용하여 실행된다.
단계(505)에서, 제2 90도 위상 편이가 생성된다. 제2 90도 위상 편이는 제2 센서 신호의 90도 위상 편이를 포함한다. 제1 90도 위상 편이에서와 같이, 90도 위상 편이는 어떠한 방식의 메커니즘이나 작동에 의해서도 실행될 수 있다.
단계(506)에서, 위상차가 제1 90도 위상 편이 및 제2 90도 위상 편이를 사용하여 제1 센서 신호 및 제2 센서 신호 사이에서 계산된다. 위상차를 계산하기 위하여 추가 정보가 사용될 수도 있다. 일 실시예에서, 위상차가 제1 센서 신호, 제2 센서 신호, 제1 90도 위상 편이 및 제2 90도 위상 편이로부터 결정된다. 위상차는 응답 신호, 즉 픽오프 신호에 있는 위상차를 포함하는데, 이는 진동하는 유량계 조립체(10)에서의 코리올리 효과에 기인하여 나타난다.
최종 위상차는 계산에 어떠한 주파수 값을 필요로 하지 않고 결정된다. 최종 위상차는 주파수를 사용하여 계산되는 위상차보다 더 빠르게 얻어질 수 있다. 최종 위상차는 주파수를 사용하여 계산된 위상차보다 더 큰 정확성을 갖는다.
단계(507)에서, 주파수가 계산된다. 본원발명에 따른 주파수는 제1 90도 위상 편이 및 제2 90도 위상 편이로부터 바람직하게 계산된다. 일 실시예에서 주파수는 90도 위상 편이 및 90도 위상 편이가 유도되는 해당 센서 신호를 사용한다. 주파수는 제1 센서 신호 및 제2 센서 신호 중 하나의 진동 응답 주파수이다(이러한 두 센서의 주파수는 작동에 있어서 실질적으로 동일하다). 주파수는 드라이버(180)에 의해 생성된 진동에 대한 유동관 또는 유동관들의 진동 주파수 응답을 포함한다.
이렇게 유도되는 주파수는 어떠한 독립적인 주파수 기준 신호에 대한 필요 없이 얻어진다. 주파수는 종래 기술에서보다 더 빠른 작동에서 90도 위상 편이로부터 얻어진다. 최종 주파수는 종래기술에서 계산된 주파수보다 더 큰 정확성을 갖는다.
단계(508)에서, 유동 물질의 질량 유량이 계산된다. 질량 유량은 단계(506) 및 단계(507) 계산된 최종 주파수 및 최종 위상차로부터 계산된다. 또한, 질량 유량 계산은 주파수 및 위상차로부터 시간차(Δt)를 계산할 수 있으며, 시간차(Δt)는 최종적으로 질량 유량을 계산하는데 사용된다.
단계(509)에서, 전술한 바와 같이 밀도가 선택적으로 결정될 수 있다.
단계(510)에서, 전술한 바와 같이 용적 유량이 선택적으로 결정될 수 있다.
도 6은 본원발명의 일 실시예에 따른 처리 시스템(203)의 일부에 대한 블록 다이어그램(600)을 도시한다. 도면에서, 블록들은 처리 회로나 처리 작동/루틴을 나타낸다. 블록 다이어그램(600)은 스테이지 1 필터 블록(stage 1 filter block; 601), 스테이지 2 필터 블록(stage 2 filter block; 602), 힐버트 변환 블록(603), 및 분석 블록(604)을 포함한다. LPO 및 RPO 입력은 좌측 픽오프 신호 입력 및 우측 픽오프 신호 입력을 포함한다. LPO 나 RPO 중 하나는 제1 센서 신호를 포함할 수 있다.
일 실시예에서, 스테이지 1 필터 블록(601) 및 스테이지 2 필터 블록(602)은 처리 시스템(203)에서 실행되는 디지털 유한 임펄스 응답(Finite Impulse Response; FIR) 다상 데시메이션 필터(polyphase decimation filter)를 포함한다. 이러한 필터는 센서 신호 중 하나 또는 둘 모두에 대한 데시메이팅 및 필터링을 위한 최적의 방법을 제공하며, 이러한 필터링 및 데시메이팅은 시간 순서에서 동일한 데시메이션 비율로 실행된다. 대안적으로 스테이지 1 필터 블록(601) 및 스테이지 2 필터 블록(602)은 무한 임펄스 응답(Infinite Impulse Response; IIR) 필터 또는 기타 적당한 디지털 필터나 필터 프로세스를 포함할 수 있다. 그러나 다른 필터링 프로세스 및/또는 필터링 실시가 고려될 수 있고 이는 본 명세서 및 청구범위에 있다는 것을 이해하여야 한다.
도 7은 본원발명의 일 실시예에 따른 힐버트 변환 블록(603)을 상세히 도시한다. 도시된 실시예에서, 힐버트 변환 블록(603)은 LPO 분기(700) 및 RPO 분기(710)를 포함한다. LPO 분기(700)는 LPO 필터 블록(702)과 평행하게 LPO 지연 블록(701)을 포함한다. 마찬가지로, RPO 분기는 RPO 필터 블록(712)과 평행하게 RPO 지연 블록(711)을 포함한다. LPO 지연 블록(701) 및 RPO 지연 블록(711)은 샘플링 지연(sampling delay)을 도입한다. 따라서 LPO 지연 블록(701) 및 RPO 지연 블록(711)은 LPO 필터 블록(702) 및 RPO 필터 블록(712)에 의해 필터링되는 LPO 및 RPO 디지털 신호 샘플보다 시간에 있어 뒤에 있는 LPO 및 RPO 디지털 신호 샘플을 선택한다. LPO 필터 블록(702) 및 RPO 필터 블록(712)은 입력된 디지털 신호 샘플에 대해 90도 위상 편이를 실행한다.
힐버트 변환 블록(603)은 위상 측정의 제공에 대한 제1 단계이다. 힐버트 변환 블록(603)는 필터링되고 데시메이팅된 LPO 및 RPO 신호를 수신하여 힐버트 변환을 실행한다. 힐버트 변환은 LPO 및 RPO 신호의 90도 위상-편이된 버전을 생성 한다. 즉 힐버트 변환은 원래의 동상(in-phase)(I) 신호 성분의 이상(quadrature)(Q) 성분을 생성한다. 따라서 힐버트 변환 블록(603)의 출력은, 원래의 동상(I) 신호 성분 LOP I 및 RPO I와 함께, 새로운 이상(Q) 성분 LPO Q 및 RPO Q 를 제공한다.
힐버트 변환 블록(603)에 대한 입력은 다음과 같이 나타낼 수 있다:
Figure 112008019671264-PCT00001
힐버트 변환을 사용하면 출력은 다음과 같이 된다:
Figure 112008019671264-PCT00002
원래의 항을 힐버트 변환의 출력과 결합하면 다음과 같다:
Figure 112008019671264-PCT00003
도 8 및 9는 본원발명의 일 실시예에 따른 분석 블록(604)의 두 개의 독립된 분기부에 대한 블록 다이어그램이다. 분석 블록(604)은 주파수, 차등 위상, 및 델타 T(Δt) 측정의 마지막 단계이다. 도 8은 동상(I) 및 이상(Q) 성분으로부터 위상차를 결정하는 제1 분기를 포함하는 위상 부분(604a)이다. 도 9는 단일 센서 신호의 동상(I) 및 이상(Q) 성분으로부터 주파수를 결정하는 주파수 부분(604b)이다. 단일 센서 신호는 도시된 바와 같이 LPO 신호를 포함할 수 있거나, 대안적으로 RPO 신호를 포함할 수 있다.
도 8의 실시예에서, 분석 블록(604)의 위상 부분(604a)은 합 블록(join block; 801a, 801b), 켤레 블록(conjugate block; 802), 복소수 곱셈 블록(complex multiplication block; 803), 필터 블록(804), 및 위상 각 블록(phase angle block; 805)을 포함한다.
합 블록(801a, 801b)은 센서 신호의 동상(I) 및 이상(Q) 성분 모두를 수신하여 전달한다. 켤레 블록(802)은 센서 신호(여기서는 LPO 신호)에 대한 켤레 복소수를 실행하여 음의 허수 신호(imaginary signal)를 형성한다. 복소수 곱셈 블록(803)은 아래의 식(8)을 실행하여 RPO 신호와 LPO 신호를 곱한다. 필터 블록(804)은 상술한 FIR 필터와 같은 디지털 필터를 실행한다. 필터 블록(804)은 센서 신호의 동상((I) 및 이상(Q) 성분으로부터 고조파함유기(harmonic content)를 제거할 뿐만 아니라 신호를 데시메이트하기 위해 사용되는 다상 데시메이션 필터를 포함할 수 있다. 필터 계수는 예를 들어 10의 인자에 의한 데시메이션과 같이 입력된 신호의 데시메이션을 제공하도록 선택될 수 있다. 위상 각 블록(805)은 LPO 신호 및 RPO 신호의 동상(I) 및 이상(Q) 성분으로부터 위상 각을 결정한다. 위상 각 블록(805)은 아래의 식(11)을 실행한다.
도 8에 도시된 위상 부분(604a)은 아래의 식을 실행한다:
Figure 112008019671264-PCT00004
여기서
Figure 112008019671264-PCT00005
는 LPO의 켤레 복소수이다.
Figure 112008019671264-PCT00006
라고 가정하면:
Figure 112008019671264-PCT00007
최종 차등 위상 각(differential phase angle)은:
Figure 112008019671264-PCT00008
도 9는 본원발명에 따른 분석 블록(604)의 주파수 부분(604b)의 블록 다이어그램이다. 주파수 부분(604b)은 좌 또는 우측 픽오프 신호(LPO 또는 RPO) 중 하나에 대해 작동할 수 있다. 도시된 실시예에서 주파수 부분(604b)은 합 블록(901), 켤레 복소수 블록(902), 샘플링 블록(903), 복소수 곱셈 블록(904), 필터 블록(905), 위상 각 블록(906), 상수 블록(907), 및 나눗셈 블록(908)을 포함한다.
전술한 바와 같이, 합 블록(901)은 센서 신호의 동상(I) 및 이상(Q) 성분 모두를 수신하여 전달한다. 켤레 블록(902)은 센서 신호(여기서는 LPO 신호)에 대한 켤레 복소수를 실행하여 음의 허수 신호(imaginary signal)를 형성한다. 지연 블록(903)은 주파수 부분(604b)으로 샘플링 지연(sampling delay)을 도입하고, 따라서 시간순으로 오래된 디지털 신호 샘플을 선택한다. 이러한 오래된 디지털 신호 샘플은 복소수 곱셈 블록(904)에서 현재의 디지털 신호와 곱해진다. 복소수 곱셈 블록(904)은 아래의 식(12)을 실행하여 LPO 신호와 LPO 켤레 신호를 곱한다. 필터 블록(905)은 상술한 FIR 필터와 같은 디지털 필터를 실행한다. 필터 블록(905)은 센서 신호의 동상((I) 및 이상(Q) 성분으로부터 고조파함유기(harmonic content)를 제거할 뿐만 아니라 신호를 데시메이트하기 위해 사용되는 다상 데시메이션 필터를 포함할 수 있다. 필터 계수는 예를 들어 10의 인자에 의한 데시메이션과 같이 입력된 신호의 데시메이션을 제공하도록 선택될 수 있다. 위상 각 블록(906)은 LPO 신호의 동상(I) 및 이상(Q) 성분으로부터 위상 각을 결정한다. 위상 각 블록(906)은 아래의 식(13)을 실행한다. 상수 블록(907)은 식(14)에 도시된 바와 같이 2π로 나누어지는 샘플율(Fs)을 포함하는 인자를 공급한다. 나눗셈 블록(908)은 식(14)의 나눗셈 과정을 실행한다.
주파수 부분(604b)은 아래의 식을 실행한다:
Figure 112008019671264-PCT00009
따라서 두 개의 연속적이 샘플 사이의 각도는:
Figure 112008019671264-PCT00010
이는 좌측 픽오프의 라디안 주파수이다. Hz 로 변환하면:
Figure 112008019671264-PCT00011
여기서 "Fs" 는 힐버트 변환 블록(603)의 비율이다. 전술한 실시예에서 "Fs"는 약 2 kHz 이다.
도 10은 정상 조건에 있는 유량계의 픽오프 신호에 대한 파워 스펙트럼 밀도 그래프(power spectrum density plot)이다. 유량계의 기본 주파수는 그래프의 가장 높은 봉(spike)이며 약 135 Hz 근방에 위치한다. 이 도면은 또한 주파수 스펙트럼 내에 몇 개의 다른 높은 봉도 도시한다(제1 비-기본 모드는 기본 모드 주파수의 약 1.5배 주파수에 있는 비틀림 모드(twist mode)이다). 이러한 봉은 유량계의 고조파(harmonic frequency)를 포함하며, 또한 기타의 바람직하지 않은 센서 모드(즉, 비틀림 모드, 제2 굽힘 모드 등등)를 포함할 수도 있다.
도 11은 단일 위상 편이 실시예에 따른 대안적인 힐버트 변환 블록(603')을 도시한다. 이 실시예의 힐버트 변환 블록(603')은 LPO 분기(1100) 및 RPO 분기(1110)를 포함한다. LPO 분기(1100)는 필터 블록(702)과 평행하게 지연 블록(701)을 포함한다. 이 실시예의 RPO 분기(1110)는 단지 지연 블록(701)만을 포함한다. 앞서와 마찬가지로, 지연 블록(701)들은 샘플링 지연을 도입한다. 앞서와 마찬가지로, 필터 블록(702)은 입력된 디지털 신호 샘플에 대해 90도 위상 편이를 실행한다. 대안적으로 힐버트 변환 블록(603')이 단지 RPO 신호를 위상 편이할 수 있다는 것을 이해해야 한다.
이러한 처리 실시예는 주파수 및 위상차 모두를 얻기 위하여 단지 하나의 센서 신호에 대한 힐버트 변환/위상 편이를 사용한다. 이는 위상 측정을 실행하는데 필요한 계산의 수를 상당히 감소시키며, 또한 질량 유량을 얻기 위해 필요한 계산의 수를 상당히 감소시킨다.
이러한 실시예에서, 힐버트 변환(603')의 출력은 좌측 또는 우측 센서 신호 중 하나의 이상(Q) 성분을 제공하나 모두에 제공하지는 않는다. 아래의 예에서, LPO 신호는 위상 편이된다.
Figure 112008019671264-PCT00012
힐버트 변환을 사용하면, 출력은 다음과 같이 된다:
Figure 112008019671264-PCT00013
LPO 원래 항을 힐버트 변환의 출력(즉, 90도 위상 편이)과 결합하면 다음과 같으며:
Figure 112008019671264-PCT00014
RPO 는 동일하다:
Figure 112008019671264-PCT00015
도 12는 단일 위상 편이 실시예에 대한 분석 블록(604a')을 도시한다. 이 실시예의 분석 블록(604a')은 하나의 합 블록(801), 복소수 곱셈 블록(803), 저 로패스 필터(low-pass filter) 블록(1201), 및 위상 각 블록(805)을 포함한다. 이 실시예의 분석 블록(604a')은 아래의 식을 실행한다:
Figure 112008019671264-PCT00016
로패스 필터 블록(1201)은 복소수 곱셈 블록(803)에 의해 생성된 고주파수 성분을 제거하는 로패스 필터를 포함한다. 로패스 필터 블록(1201)은 어떠한 방식의 로패스 필터링 작동도 실행할 수 있다. 곱셈 과정은 두 개의 항을 생성한다. (-ωt + ωt + φ)항은 결합되어 위상만의 φ항(DC 결과)으로 간략해지는데, 이는 (-ωt)항과 (ωt) 항이 서로 상쇄되기 때문이다. (ωt + ωt + φ)항은 2배의 주파수인 (2ωt + φ)항으로 간략화된다. 결과는 2 항의 합이므로, 고주파수 (2ωt + φ) 항은 제거될 수 있다. 여기서 관심 있는 유일한 항은 DC 항이다. 고주파수 (2ωt + φ) 항은 로패스 필터를 사용하여 결과로부터 필터링될 수 있다. 로패스 필터의 컷오프(cut-off)는 0과 2ω 사이의 어떠한 곳에도 위치할 수 있다.
필터링 후의 결과는 다음과 같다:
Figure 112008019671264-PCT00017
따라서, 차등 위상 각은 다음과 같다.
Figure 112008019671264-PCT00018
두 개 대신에 하나의 픽오프 신호에 대해 힐버트 변환을 취함으로써, 코리올리 질량 유량계에서 위상 및 주파수 추정을 실행하기 위해 필요한 계산 부하는 바람직하게 감소된다. 따라서 위상 및 주파수는 두 개의 센서 신호를 사용하기는 하지만 단지 하나의 90도 위상 편이를 사용하여 결정될 수 있다.
도 13은 종래 기술과 비교한 본원발명의 센서 처리를 도시하며, 여기서는 각각의 시간차(Δt) 값이 비교된다. 표는 가스 유동(즉, 예를 들어 가스 기포)을 포 함하는 유동 물질을 도시한다. 이러한 조건에서, 유동 노이즈는 위상 및 주파수 계산의 속도로 인하여 새로운 알고리즘에서 실질적으로 감소된다. 본원발명에 의해 얻어진 결과는 종래 기술의 (Δt) 측정에서 반영되는 높은 마루 및 골(peak and valley)들을 나타내지 않는다는 것을 그래프로부터 알 수 있다.
도 14는 본원발명의 다른 실시예에 따른 계측 전자장치(20)를 도시한다. 이 실시예의 계측 전자장치(20)는 앞서 설명한 바와 같이 인터페이스(201), 디지타이저(202), 처리 시스템(203), 및 저장 시스템(204)을 포함할 수 있다. 다른 실시예와 공통되는 구성요소 및/또는 루틴에 대해서는 동일한 도면부호를 사용한다. 이러한 구성의 계측 전자장치(20)는 앞서 언급된 것들과 같이 다양한 구성요소 및/또는 루틴을 포함할 수 있다는 것을 이해해야 한다.
작동에 있어서, 계측 전자장치(20)는 유량계(5)를 통해 유동하는 유동 물질에 있는 유동 성분의 질량 분율을 결정하기 위하여 유량계 조립체(10)로부터의 제1 센서 신호(210) 및 제2 센서 신호(211)를 수신하여 처리한다. 질량 분율은 2상 유동 흐름에 있는 제1 유동 성분 및 제2 유동 성분 간의 질량 유동의 비율이다. 질량 분율은 다양한 유동 성분의 질량을 결정하는데 사용될 수 있다. 예를 들어, 유동은 유체 성분 및 기체 성분을 포함할 수 있다. 유체 성분 질량 유량 및 가스 성분 질량 유량 중 하나 이상을 유도하기 위하여 유동 물질의 전체 질량 유량에 질량 분율이 곱해질 수 있다. 유체는 어떠한 양식의 유체도 포함할 수 있으며, 가스는 어떠한 양식의 가스도 포함할 수 있다. 가스는 예를 들어 공기를 포함할 수 있다. 이하의 설명은 유체 내의 공기에 초점을 두고 있으나 본원발명은 어떠한 가스에도 적용될 수 있음을 이해해야 할 것이다.
계측 전자장치(20)는 제1 센서 신호(210) 및 제2 센서 신호(211)를 수신하여 처리한다. 하나의 처리 작용은, 전술한 바와 같이, 센서 신호들 중 하나 또는 이들 모두로부터 주파수 응답(1410)을 생성하는 것이다. 주파수 응답(1410)은 공급된 구동 진동에 응답하여 유량계 조립체(10)의 주파수를 포함한다. 계측 전자장치(20)는 주파수 응답(1410)을 가스 주파수 성분(1412) 및 유체 주파수 성분(1416)으로 분리(break out)한다. 계측 전자장치(20)는 주파수 응답(1410)m로부터 전체 밀도(ρmix)를 결정한다. 마찬가지로, 가스 성분 밀도(ρgas)(1421)는 가스 주파수 성분(1412)으로부터 결정된다. 계측 전자장치(20)는 주파수 응답(110)과 가스 주파수 성분(1412) 및 유체 주파수 성분(1416) 중 하나 이상을 이용하여 가스의 보이드(void) 분율(1418)을 결정한다. 계측 전자장치(20)는 질량 분율(1419)을 결정하기 위하여 보이드 분율(1418), 전체 밀도(1420), 및 가스 밀도(1421)를 추가로 이용한다. 질량 분율(mf)은 다음과 같이 정의된다:
Figure 112008019671264-PCT00019
일 실시예에서, 질량 분율을 가스의 질량 분율(mfgas)을 포함한다. 가스의 질량 분율은 다음과 같다:
Figure 112008019671264-PCT00020
그러나 본원발명은 대안적으로 유동 물질 내의 유체의 질량 분율(mffluid), 또는 임의의 모든 질량 분율을 결정할 수도 있다는 것을 이해해야 한다. 유체의 질량 분율(mffluid)은 가스의 질량 분율의 보완(complement)을 포함한다:
Figure 112008019671264-PCT00021
그러나 이러한 본 명세서는 설명의 명료성을 위하여 가스의 질량 유동(mfgas)에 초점을 맞출 것이다.
처리 시스템(203)은 제1 센서 신호(210) 및 제2 센서 신호(211)를 사용하여 가스 유동 물질의 전체 밀도(1420)를 결정하고, 전체 밀도(1420)를 가스 유동 물질의 가스 유동 분율을 나타내는 가스 밀도(1421) 또는 액체 유동 분율을 나타내는 유체 밀도(1422) 중 하나 이상에 비교하고, 가스 밀도(1421) 및 유체 밀도(1422) 중 하나 이상과 전체 밀도(1420)로부터 액체 분율(1427)을 결정하도록 추가적으로 구성될 수 있다. 몇몇 실시예에서, 전체 밀도(1420)는 가스 밀도(1421) 및 유체 밀도(1422) 모두에 비교된다. 또한, 처리 시스템(203)은 주파수 응답(1410)을 사용하여 가스 유동 물질의 전체 밀도(1420)를 결정하고, 전체 밀도(1420)로부터 액체 분율(1427) 및/또는 가스 분율(1428)을 결정하도록 구성될 수 있다.
가스 유동 물질은 어떠한 가스 또는 가스 혼합물도 포함할 수 있다. 한가지 형태의 가스 유동 물질은 오일 유정으로부터 추출된 가스를 포함한다. 예를 들어, 가스 유동 물질은 천연 가스를 포함할 수 있다. 그러나, 다른 가스나 가스 혼합물 도 고려될 수 있으며, 이는 상세한 설명 및 청구항의 범위 내에 놓인다.
액체 유동 분율은 가스 유동 물질 내에 혼입된 임의의 액체를 포함할 수 있다. 액체 유동 분율은 액체 흐름, 액체 증기, 또는 액체 방울(drop)을 포함할 수 있다. 예를 들어, 액체는 천연 가스 유동 내의 물을 포함할 수 있다. 대안적으로, 액체는 건조 공정으로부터의 글리콜(glycol) 캐리오버와 같이, 천연 가스 유동 내의 글리콜을 포함할 수 있다. 다른 대안으로서, 액체는 펌프, 조절기, 또는 다른 유동 취급 메커니즘에 의하여 가스 유동 물질 내로 유입되는 오일과 같이, 가스 유동 물질 내의 기름을 포함할 수 있다. 그러나 다른 액체나 액체 조합이 고려될 수도 있으며, 이러한 것 역시 본원발명의 상세한 설명 및 청구항의 범위 내에 속하게 된다.
유량계 조립체(10)의 주파수는 액체 유동 분율이 유량계 조립체(10)를 통과할 때 변하게 된다. 순수 가스 유동과는 대조적으로, 주파수는 액체 유동 분율이 증가함에 따라 감소한다. 이는 액체 유동 분율이 존재할 때 증가되는 밀도에 기인한다. 따라서, 유동 흐름 밀도는 액체 유동 분율을 결정하기 위해 사용될 수 있다.
제1 센서 신호(210) 및 제2 센서 신호(211)는 픽오프 센서(170L 및 170R)로부터의 신호와 같이 계측 전자장치(20)에 의해 실질적으로 연속적으로 수신되어 처리되는 시간-변화(time-varying) 전자 신호를 포함한다. 주파수 응답(1410)은 전술한 처리 블록(도 6-7 및 도 9 참조)을 사용하여 결정될 수 있다. 바람직하게, 전술한 고속 주파수 결정을 사용할 때, 본원발명은 빠르고, 정확하고, 신뢰성 있게 가스의 보이드 분율(1418)을 결정할 수 있다.
이 실시예의 처리 시스템(203)은 보이드 분율 루틴(1401), 노치 필터 루틴(1402), 질량 분율 루틴(1405), 가스 캐리언더 루틴(1406), 액체 캐리오버 루틴(1407), 밀도 루틴(1408), 및 액체 유동 분율 루틴(1409)을 포함할 수 있다. 처리 시스템(203)은 로패스 필터 루틴(1403) 및 하이 패스 필터 루틴(1404)과 같이, 하나 이상의 필터 또는 필터 루틴을 더 포함할 수 있다. 대안적으로, 하나 이상의 필터 또는 필터 루틴은 노치 필터 구성이나 좁은 주파수 대역을 차단하는 기타의 필터 구성을 포함할 수 있다. 처리 시스템(203)은 주파수 응답 측정, 보이드 분율 결정, 및 질량 분율 결정을 각각 저장할 수 있는 주파수 응답(1410) 및 보이드 분율(1418)을 더 포함할 수 있다. 처리 시스템(203)은 보이드 분율 및 질량 분율 결정에 대한 작동 주파수 값을 저장하는 가스 주파수 성분(1412) 및 유체 주파수 성분(1416)을 더 포함할 수 있다. 처리 시스템(203)은 보이드 분율 및 질량 분율 결정에 대한 작동 밀도 값들을 저장하는 전체 밀도(203), 가스 성분 밀도(1421), 및 유체 성분 밀도(1422)를 더 포함할 수 있다. 처리 시스템(203)은 각각의 캐리언더/캐리오버 양을 저장하는 가스 캐리언더(1425) 및 액체 캐리오버(1426)를 더 포함할 수 있다. 처리 시스템(203)은 가스 성분 양을 저장하는 가스 분율(1428) 및 액체 성분 양을 저장하는 액체 분율(1427)을 더 포함할 수 있다.
주파수 응답(1410)은 혼합 주파수(fmix)를 포함하는데, 여기서 주파수 응답(1410)은 가스 주파수 성분(fgas)(1412) 및 유체 주파수 성분(ffluid)(1416)을 포함 할 수 있다. 보이드 분율 및 질량 분율은 이러한 주파수 성분이 혼합 주파수(fmix)로부터 분리되어 결정된 이후에 결정될 수 있다. 언제라도, 주파수 응답(1410)은 임의의 양의 가스 주파수 성분(fgas)(1412), 즉 혼입 가스로 인한 주파수 성분을 포함할 수 있다.
도 15는 공기, 유체, 및 공기/유체 혼합물(즉, 혼입 공기를 포함하는 유체)에 대한 유량계 주파수 응답의 그래프를 도시한다. 가스의 밀도는 유량계를 통해 유동하는 유동 물질의 유체 밀도와 구별될 수 있다. 밀도는 측정 주파수로부터 유도될 수 있으므로, 공기와 관련된 주파수는 유체의 주파수로부터도 구별될 수 있다. 이는 다른 가스나 가스 혼합물에 대해서도 역시 적용된다.
주파수를 계산하기 위한 방정식은 다음과 같다:
Figure 112008019671264-PCT00022
여기서 ω는 코리올리 유량계의 라디안 주파수이다. ω-1 항은 앞의 또는 그 전의 샘플 기간으로부터 라디안 주파수 샘플을 나타낸다. 라디안 주파수 ω를 헤르쯔(Hz)의 주파수 f 로 변환하면:
Figure 112008019671264-PCT00023
이 방정식은 오직 하나의 주파수만 존재하는 것을 가정한다. 혼입 공기의 경우에서와 같이(공기의 주파수 및 유동 물질 유체의 주파수) 만약 두 개의 주파수 가 존재하게 되면, 새로운 방정식은 다음과 같다:
Figure 112008019671264-PCT00024
여기서 fmix 는 가스 주파수 성분(fgas) 및 유체 주파수 성분(ffluid)을 포함하는 전체 유동 물질의 주파수 응답을 나타낸다.
도 14를 다시 참조하면, 로패스 필터 루틴(1403)은 로패스 필터를 실행한다. 로패스 필터는 로패스 컷오프 주파수보다 실질적으로 작은 저 주파수를 통과시킨다. 따라서 로패스 필터는 고 주파수를 제거하기 위하여 사용될 수 있다.
하이패스 필터 루틴(1404)은 하이패스 필터를 실행한다. 하이패스 필터는 하이패스 컷오프 주파수보다 실질적으로 높은 고 주파수를 통과시킨다. 따라서 하이패스 필터는 저 주파수를 제거하기 위하여 사용될 수 있다.
노치 필터 루틴(1402)은 노치 필터를 실행한다. 노치 필터는 노치 필터의 주파수 응답에 있는 "노치(notch)"에 중심을 두는 좁은 대역의 주파수를 제거한다. 노치에 있는 주파수만이 노치 필터에 의해 차단된다. 따라서 노치 필터는 알려진 바람직하지 않은 주파수를 주파수 응답(1410)으로부터 제거하는데 매우 유용하다.
보이드 분율 루틴(1401)은 유동 물질 내의 (통상적으로 가스의) 보이드 분율을 결정한다. 보이드 분율은 유동 성분의 밀도로부터 결정될 수 있으며, 여기서 전체 밀도(ρmix)는 가스 성분 밀도(ρgas) 및 유체 성분 밀도(ρfluid)를 포함한다.
밀도(ρ)는 실질적으로 다음과 같다:
Figure 112008019671264-PCT00025
여기서 f 는 유체 주파수 성분(1416)(즉, fmix)의 주파수 측정치이다. 유체 성분 밀도(ρfluid)(1422)는 유체 주파수 성분(1416)을 사용하여 계산될 수 있다. 일 실시예에서, 유체 주파수 성분(1416)은 평균 혼합 주파수를 포함한다. 가스 성분 밀도(ρgas)(1421)는 가스 주파수 성분(1412)을 사용하여 계산될 수 있다. 결과적으로, 가스의 보이드 분율(1418)은 유체 성분 밀도(ρfluid)(1422)에서 전체 밀도(ρmix)(1420)를 뺀 값을 유체 성분 밀도(ρfluid)(1422)에서 가스 성분 밀도(ρgas)(1421)를 뺀 값으로 나눈 값으로서 계산된다. 이러한 계산은 다음과 같다:
Figure 112008019671264-PCT00026
결과적인 가스의 보이드 분율(1418)은 유동 물질 내의 유체에 대한 가스의 비율을 반영한다.
질량 분율 루틴(1405)은 주파수 응답(1410)으로부터 질량 분율(1419)을 결정한다. 일 실시예에서, 질량 분율 루틴(1405)은 질량 분율(1419)을 계산하기 위하여, 유도된 밀도값과 함께, 결정된 보이드 분율(VF)을 사용한다.
질량(m) 및 용적(V)은 밀도(ρ)에 의해 관련된다. 따라서 밀도는 다음과 같 다:
Figure 112008019671264-PCT00027
결과적으로, 질량 분율(mf)은 다음과 같이 간략화될 수 있다:
Figure 112008019671264-PCT00028
보이드 분율(VF)은 용적 비율과 같으므로:
Figure 112008019671264-PCT00029
질량 분율(mf)은 다음과 같다:
Figure 112008019671264-PCT00030
결과적으로, 질량 분율은 보이드 분율(VF), 가스 성분 밀도(ρgas)(1421), 및 전체 밀도(ρmix)(1422)로부터 결정될 수 있다. 가스 성분 밀도(ρgas)(1421) 및 전체 밀도(ρmix)(1422)는 각각 가스 주파수 성분(1412) 및 주파수 응답(1410)으로부터 결정될 수 있다.
가스나 유체 중 하나가 알려져 있다면 가스 주파수 성분(1412) 및 유체 주파수 성분(1416) 중 하나만이 필요하다는 것을 이해해야 한다. 예를 들어, 가스가 공기를 포함하면, 표준 공기 주파수 응답(및 밀도)은 가정될 수 있다. 결과적으로, 공지된 가스나 유체 주파수는 필터링될 수 있으며, 오직 하나의 필터링 단계만 이 요구된다.
가스 캐리언더 루틴(1406)은 다상 유동 물질 내의 가스 캐리언더 양을 결정한다. 일 실시예의 가스 캐리언더 루틴(1406)은 가스의 보이드 분율(VF)을 결정함으로써 가스 캐리언더 양을 결정한다. VF 값은 전체 유동 물질 용적의 용적 백분율을 포함한다. 따라서 VF는 가스 캐리언더 양을 포함하거나, 또는 다상 유동 물질 내의 가스 질량 유량(mfGAS)을 결정하도록 조작될 수 있다. 가스 성분의 질량 유량이 계산된다. 결정된 가스 캐리언더 양은 가스 캐리언더(1425)에 저장될 수 있다.
액체 캐리오버 루틴(1407)은 다상 유동 물질 내의 액체 캐리오버 양을 결정한다. 액체 캐리오버 루틴(1407)은 액체 분율을 결정함으로써 액체 캐리오버 양을 결정하며 다상 유동 물질 내에서 액체 성분 밀도, 액체 질량 유량 등등 중 하나 이상을 결정하기 위하여 액체 분율을 이용한다. 결정된 액체 캐리오버 양은 액체 캐리오버(1426) 내에 저장될 수 있다.
밀도 루틴(1408)은 제1 센서 신호(210) 및 제2 센서 신호(211)로부터 전체 밀도(1420)를 결정한다. 일 실시예에서, 밀도 루틴(1408)은 유동 흐름 밀도 결정의 주파수 응답(1410)을 이용한다. 전체 밀도(1420)는 가스 유동 물질의 밀도에 관련되며, 가스 유동 물질 내의 액체 유동 분율에 따라 변한다.
액체 유동 분율 루틴(1409)은 전체 밀도(1420)를 이용하여 액체 분율(1427)을 생성한다. 액체 분율(1427)은 가스 유동 물질 내의 액체의 양(또는 백분율)에 관련된다. 이러한 결정은 이하에서 도 15와 연관하여 설명된다. 또한, 액체 유동 분율 루틴(1409)은 가스 분율(1428)을 생성할 수도 있으며, 여기서 가스 분율(1428)은 가스 유동 물질 내의 가스의 양/백분율에 관련된다.
계측 전자장치(20)는 전체 질량 유량, 성분 질량 유동, 성분 용적 등과 같은 기타의 유동 특성을 추가적으로 결정할 수 있다. 계측 전자장치(20)는 유량계 조립체(10)와 통신할 수 있으며, 여기서 유량계 조립체(10)는 주파수 응답을 생성하는 어떠한 방식의 유량계도 포함할 수 있다. 일 실시예에서, 유량계 조립체(10)는 코리올리 유량계를 포함한다. 다른 실시예에서, 유량계 조립체(10)는 진동 농도계(vibrating densitometer)를 포함한다.
도 16은 본원발명의 일 실시예에 따라 유량계를 통과하여 유동하는 유동 물질 내에 있는 가스의 보이드 분율을 결정하기 위한 방법의 순서도(1600)이다. 단계(1601)에서, 주파수 응답이 수신된다. 주파수 응답은 예를 들어 계측 전자장치(20)에 수신될 수 있다. 주파수 응답은 유동 물질을 포함하는 진동 유량계 조립체(10)에 대한 주파수 응답을 포함한다. 유동 물질은 혼입 가스를 포함할 수 있다.
일 실시예에서, 주파수 응답은 제1 센서 신호 및 제2 센서 신호를 포함할 수 있다. 제1 센서 신호 및 제2 센서 신호는 예를 들어 유량계 조립체(10)의 픽오프 센서(170L 및 170R)로부터 수신될 수 있다. 센서 신호로부터 90도의 위상 편이가 생성될 수 있다. 90도 위상 편이 및 제1 및 제2 센서 신호는 주파수 응답을 결정하는데 사용될 수 있으며, 여기서 주파수 응답은 유동 물질의 질량 유량 및 혼입 가스의 양과 존재에 따라 시간에 따라 변한다.
단계(1602)에서, 주파수 응답은 가스 주파수 성분(1412)과 유체 주파수 성분(1416)으로 분리된다. 주파수 응답(1410)이 유동 물질 내의 가스 유량에 관련되는 가스 주파수 성분과 유체 유량에 관련되는 유체 주파수 성분을 포함하기 때문에 이러한 것이 가능하다. 유체는 임의 형태의 유체를 포함할 수 있다. 이러한 분리(breaking out)는 앞서 설명한 바와 같이 제1 및 제2 필터에 의해 실행될 수 있다. 또한, 분리는 앞서 설명한 바와 같이 로패스 필터 및 하이패스 필터에 의해 실행될 수 있다.
단계(1603)에서, 앞서 설명한 바와 같이, 가스의 보이드 분율(1418)은 주파수 응답(1410), 가스 주파수 성분(1412), 및 유체 주파수 성분(1416)을 사용하여 결정된다. 이러한 결정은 앞서 설명한 바와 같이, 주파수 응답(1410), 가스 주파수 성분(1412), 및 유체 주파수 성분(1416)으로부터 밀도값을 결정하는 과정을 포함할 수 있다. 결과적인 가스의 보이드 분율(1418)은 비율, 백분율, 또는 기타의 수단으로 표현될 수 있다.
도 17은 본원발명의 실시예에 따라 유량계를 통해 유동하는 유동 물질에 있는 가스의 보이드 분율을 결정하기 위한 방법의 순서도(1700)를 도시한다. 단계(1701)에서, 앞서 설명한 바와 같이 주파수 응답이 수신된다.
단계(1702)에서, 주파수 응답이 제1 필터로 처리된다. 제1 필터는 가스 주파수 성분을 실질적으로 차단하고 유체 주파수 성분을 실질적으로 통과시킨다(도 18참조). 일 실시예에서, 제1 필터는 로패스 필터를 포함하며, 여기서 로패스 필 터의 로패스 컷오프 주파수는 유체 주파수 성분보다 실질적으로 높다. 결과적으로, 로패스 필터는 유체 주파수 성분을 실질적으로 통과시키고 가스 주파수 성분을 실질적으로 차단한다.
단계(1703)에서, 주파수 응답은 제2 필터로 처리된다. 제2 필터는 유체 주파수 성분을 실질적으로 차단하며 가스 주파수 성분을 실질적으로 통과시킨다. 일 실시예에서, 제2 필터는 하이패스 필터를 포함하며, 여기서 하이패스 필터의 하이패스 컷오프 주파수는 가스 주파수 성분보다 실질적으로 낮다(그러나 유체 주파수 성분보다는 높다). 결과적으로, 하이패스 필터는 가스 주파수 성분을 실질적으로 통과시키고 유체 주파수 성분을 실질적으로 차단한다.
단계(1704)에서, 앞서 설명한 바와 같이, 가스의 보이드 분율(1418)은 주파수 응답(1410), 가스 주파수 성분(1412), 및 유체 주파수 성분(1416)을 사용하여 결정된다.
도 18은 본원발명의 일 실시예에 따라 가스 주파수 성분 및 유체 주파수 성분을 분리하는데 사용될 수 있는 로패스 및 하이패스 필터 응답을 도시하는 주파수 그래프이다. 그래프의 하부에 있는 라인은 유체 주파수 성분 로브(lobe) 및 가스 주파수 성분 로브를 포함하는 유량계 주파수 응답을 나타내는 그래프이다. 유체 주파수 성분 로브는 가스 주파수 성분 로브보다 주파수가 낮다. 상부의 라인은 컷오프 주파수와 함께 하이패스 주파수 응답 및 로패스 주파수 응답을 포함한다. 여기서, 로패스 및 하이패스 필터 모두에 대한 컷오프 주파수는 두 개의 로브 사이에 실질적으로 중심을 두고 있다. 로패스 및 하이패스 필터는 유체 및 가스 주파수 성분에 따라서 서로 다른 컷오프 주파수를 가지거나 공통의 컷오프 주파수를 가질 수 있다. 로패스 필터는 유체 주파수 성분을 출력하며 하이패스 필터는 가스 주파수 성분을 출력한다는 것을 알 수 있다. 따라서, 두 개의 필터는 주파수 응답(1410)을 가스 주파수 성분(1412) 및 유체 주파수 성분(1416)으로 분리할 수 있다.
유체 및 가스 주파수 성분을 분리하는 다른 방법은 하나의 공지된 주파수 성분을 필터링 하는 과정과 필터 작용에 의해 통과된 주파수 성분을 사용하여 유체 및 가스 주파수 밀도를 결정하는 과정을 포함한다. 예를 들어, 유동 물질 내의 가스가 공기인 경우, 필터링 작동은 통상적인 공기 주파수 응답에 중심을 두는 비교적 좁은 주파수 대역을 필터링하도록 구성될 수 있다. 후속적으로, 주파수 응답으로부터 유도된 전체 밀도와 잔존하는 유체 주파수 성분으로부터 유도된 유체 밀도 성분이 공기 밀도 항(term)을 결정하기 위해 사용될 수 있다. 예를 들어, 가스가 대기의 공기라고 알려진 경우, (예를 들어 노치 필터와 같은) 필터가 주파수 응답의 공기 주파수 성분을 실질적으로 차단하기 위해 사용될 수 있다. 결과적으로 전체 밀도(ρmix)(1420)는 주파수 응답(1410)으로부터 계산될 수 있으며, 유체 성분 밀도(ρfluid)(1422)는 유체 주파수 성분(1416)으로부터 계산될 수 있다. 따라서, 공기 성분 밀도(ρgas)(1421)는 다음과 같다:
Figure 112008019671264-PCT00031
상기 식은 다음과 같이 고쳐질 수 있다:
Figure 112008019671264-PCT00032
대안적으로, 유체 주파수 성분은 제거/필터링될 수 있으며 보이드 분율은 가스 주파수 성분을 사용하여 결정될 수 있다는 것을 이해해야 한다. 앞서와 같이, 이러한 단일 주파수 제거는 유체가 공지된 특성의 주파수 응답 및 밀도를 갖는 경우에 실행될 수 있다. 따라서 단일 주파수 제거 방법은 가스 주파수 성분 또는 유체 주파수 성분을 제거할 수 있다.
일 실시예에서, 단일 주파수 성분은 하나 이상의 필터에 의해 제거될 수 있는 반면, 다른 주파수 성분은 필터링 과정을 통해 통과된다. 일 실시예에서 하나 이상의 필터는 노치 필터를 포함한다. 노치 필터는 좁은 대역(즉, 주파수 응답의 노치) 내의 주파수를 제외한 모든 주파수를 통과시킨다. 대안적으로, 하나 이상의 필터는 임의의 적당한 필터 또는 필터의 조합을 포함할 수 있다.
도 19는 본원발명의 일 실시예에 따라 유량계를 통해 유동하는 유동 물질 내에 있는 가스의 보이드 분율을 결정하기 위한 방법의 순서도(1900)를 도시한다. 단계(1901)에서, 주파수 응답(1410)은 앞서 설명한 바와 같이 수신된다.
단계(1902)에서, 주파수 응답이 노치 필터로 처리된다. 노치 필터는 이 실시예에서 가스 주파수 성분 위 및 아래와 같이, 노치의 위 및 아래의 주파수를 통과시킨다. 따라서, 노치 필터는 가스 주파수 성분(1412)을 실질적으로 차단한다. 노치 필터는 유체 주파수 성분(1416)을 실질적으로 통과시킨다.
도 20은 노치 필터 주파수 응답의 그래프이다. 도시된 실시예에서, 노치는 가스 주파수에 중심이 맞춰진다. 노치 필터는 노치의 위 및 아래에 있는 모든 주파수를 실질적으로 통과시키며 오직 가스 주파수만이 노치 필터에 의해 실질적으로 차단된다.
다시 도 19를 참조하면, 단계(1903)에서, 앞서 설명한 바와 같이, 가스의 보이드 분율(1418)은 주파수 응답(1410), 가스 주파수 성분(1412), 및 유체 주파수 성분(1416)을 사용하여 결정된다.
도 21은 본원발명의 일 실시예에 따라 유량계를 통해 유동하는 유동 물질의 유동 성분의 질량 분율을 결정하기 위한 방법의 순서도(2100)를 도시한다. 단계(2101)에서, 전술한 바와 같이 주파수 응답이 수신된다.
단계(2102)에서, 전술한 바와 같이 주파수 응답이 가스 주파수 성분(1412) 및 유체 주파수 성분(1416)으로 분리된다.
단계(2103)에서, 전체 밀도(ρmix)가 주파수 응답으로부터 결정된다. 전체 밀도(ρmix)는 결합된 유체 및 가스 유동 성분의 밀도를 반영한다. 전술한 바와 같이, 전체 밀도(ρmix)는 실질적으로 1을 주파수 응답으로 나눈 값(즉, 주파수 응답의 역수)을 제곱한 값을 포함한다.
단계(2104)에서, 가스 성분 밀도(ρgas)가 가스 주파수 성분(fgas)으로부터 결정된다. 가스 성분 밀도(ρgas)는 정확히 가스 유동 성분의 밀도를 반영한다.
단계(2105)에서, 전술한 바와 같이, 가스의 보이드 분율(VF)(1418)이 주파수 응답(1410), 가스 주파수 성분(1412), 및 유체 주파수 성분(1416)을 사용하여 결정된다. 결과적인 가스의 보이드 분율(1418)은 비율, 백분율, 또는 다른 치수로서 표현될 수 있다.
단계(2106)에서, 식(46)에 표현된 바와 같이, 보이드 분율(VF)(1418) 및 전체 밀도(ρmix)에 대한 가스 밀도(ρgas)의 비율로부터 질량 분율이 결정된다.
도 22는 본원발명의 일 실시예에 따라 유량계를 통해 유동하는 유동 물질의 유동 성분의 질량 분율을 결정하기 위한 방법의 순서도(2200)를 도시한다. 주파수 응답으로부터 유체 및 가스 주파수 성분을 분리하는 한 방법은 두 개의 필터링 작업을 수행하는 단계를 포함한다. 제1 필터링 작업은 가스 주파수 성분은 실질적으로 차단하고 유체 주파수 성분은 실질적으로 통과시키는 제1 필터로 주파수 응답을 필터링하는 단계를 포함한다. 제2 필터링 작업은 유체 주파수 성분은 실질적으로 차단하고 가스 주파수 성분은 실질적으로 통과시키는 제2 필터로 주파수 응답을 필터링하는 단계를 포함한다. 결과적으로, 제1 필터는 유체 주파수 성분을 출력하는 반면에, 제2 필터는 가스 주파수 성분을 출력한다.
단계(2201)에서, 전술한 바와 같이 주파수 응답이 수신된다.
단계(2202)에서, 전술한 바와 같이 주파수 응답이 제1 필터로 필터링된다.
단계(2203)에서, 전술한 바와 같이 주파수 응답이 제2 필터로 필터링된다.
단계(2204)에서, 전술한 바와 같이 전체 밀도(ρmix)가 결정된다.
단계(2205)에서, 전술한 바와 같이 가스 밀도(ρgas)가 결정된다.
단계(2206)에서, 전술한 바와 같이, 가스의 보이드 분율(1418)이 주파수 응답(1410), 가스 주파수 성분(1412), 및 유체 주파수 성분(1416)을 사용하여 결정된다.
단계(2207)에서, 전술한 바와 같이 질량 분율(1419)이 결정된다.
도 23은 본원발명의 일 실시예에 따라 유량계를 통해 유동하는 유동 물질의 유동 성분의 질량 분율을 결정하기 위한 방법의 순서도(2300)를 도시한다. 단계(1901)에서, 전술한 바와 같이 주파수 응답(1410)이 수신된다.
단계(2302)에서, 전술한 바와 같이 주파수 응답이 노치 필터로 처리된다.
단계(2303)에서, 전술한 바와 같이 전체 밀도(ρmix)가 결정된다.
단계(2304)에서, 전술한 바와 같이 가스 밀도(ρgas)가 결정된다.
단계(2305)에서, 전술한 바와 같이, 가스의 보이드 분율(1418)이 결정된다.
단계(2306)에서, 전술한 바와 같이 질량 분율(1419)이 결정된다.
도 24는 본원발명의 일 실시예에 따라 유량계를 통해 유동하는 가스 유동 물질의 액체 유동 분율을 결정하기 위한 방법의 순서도(2400)를 도시한다. 단계(2401)에서, 전술한 바와 같이 유량계 조립체(10)로부터 제1 및 제2 센서 신호가 수신된다.
단계(2402)에서, 전술한 바와 같이, 센서 신호가 조정될 수 있다.
단계(2403)에서, 가스 유동 물질의 유동 흐름 밀도(flow stream density)가 결정된다. 유동 흐름 밀도는 제1 센서 신호 및 제2 센서 신호를 사용하여 결정된다. 다른 유동 특성, 변수, 및/또는 상수도, 필요하다면, 결정에 사용될 수 있다.
단계(2404)에서, 유동 흐름 밀도가 가스 밀도 및 액체 밀도와 비교된다. 가스 밀도는 가스 유동 흐름의 가스 유동 분율을 나타내며 액체 밀도는 액체 유동 분율을 나타낸다. 일 실시예에서, 양 밀도가 알려지고 사용된다. 다른 실시예에서는, 이들 두 밀도 중 하나만이 알려져서 사용된다.
도 25는 천연 가스 밀도 대 글리콜 백분율(즉, 액체 유동 분율)의 그래프이다. 표는 비교에 사용되는 데이터 세트를 나타낸다. 그러나 어떠한 형태의 데이터 구조도 사용될 수 있으며, 데이터가 표 형태일 필요는 없다. 사선의 그래프 선은 글리콜(glycol) 대(versus) 가스의 다양한 분율에 대한 가스 유동 물질 밀도를 나타낸다. 그래프로부터 전체 가스 유동 물질 밀도가 액체 유동 분율에 비례한다는 것을 알 수 있다. 따라서, 적어도 측정된 가스 유동 물질(즉, 전체) 밀도 및 알려진 가스 밀도를 사용하여 액체 유동 분율이 결정될 수 있다. 데이터 순람(lookup)은 예를 들어 압력 및 온도와 같은 다른 인자를 포함할 수 있다는 것을 이해해야 한다. 상기 그래프가 비록 천연 가스와 글리콜에 대한 것이지만, 다른 가스들이나 액체들도 사용되고 결정될 수 있다는 것을 이해해야 한다.
도 24를 다시 참조하면, 단계(2405)에서 액체 유동 분율이 결정된다. 액체 유동 분율은 가스 유동 물질에 있는 액체의 양 또는 백분율을 포함한다. 액체 유동 분율은 예를 들어 가스 질량 유량 및/또는 액체 질량 유량과 같은 다른 계산을 수행하기 위하여 후속적으로 사용될 수 있다. 액체 유동 분율은 적어도 알려진 가 스 밀도에 대한 유동 흐름 밀도의 비교로부터 결정될 수 있다. 대안적으로, 이러한 비교는 유동 흐름 밀도는 알려진 가스 밀도 및 알려진 액체 밀도 모두에 비교하는 것을 포함할 수 있다.
유동 흐름 밀도를 알려진 가스 밀도 및 알려진 액체 밀도에 비교함으로써, 액체 유동 분율 결정에 사용될 수 있는 비율이 형성될 수 있다. 이러한 비율은 다음과 같다:
Figure 112008019671264-PCT00033
여기서 Dmeasured 는 센서 신호로부터 결정된 유동 흐름 밀도이며, DG 는 알려진 가스 밀도, DL 은 알려진 액체 밀도, 그리고 X 는 액체 유동 분율이다. 결과적으로, 액체 유동 분율(X)은 다음과 같이 결정될 수 있다:
Figure 112008019671264-PCT00034
도 26은 본원발명의 일 실시예에 따라 유량계를 통해 유동하는 가스 유동 물질의 액체 유동 분율을 결정하기 위한 방법의 순서도(2600)를 도시한다. 단계(2601)에서, 전술한 바와 같이 유량계 조립체(10)로부터 제1 및 제2 센서 신호가 수신된다.
단계(2602)에서, 전술한 바와 같이 센서 신호가 조정될 수 있다.
단계(2603)에서, 제1 센서 신호로부터 90도 위상 편이가 생성된다. 90도 위상 편이는 전술한 바와 같이 생성될 수 있다. 비록 제1 센서 신호가 예로서 위상 편이되었지만, 양 센서 신호가 사용될 수 있다는 것을 이해해야 한다.
단계(2604)에서, 유량계 조립체(10)의 주파수 응답이 계산된다. 전술한 바와 같이, 주파수 응답은 90도 위상 편이 및 제1 센서 신호를 사용하여 계산될 수 있다.
단계(2605)에서, 가스 유동 물질의 유동 흐름 밀도는 주파수 응답을 사용하여 결정된다. 일 실시예에서, 유동 흐름 밀도는 주파수 응답을 제곱하고 제곱된 주파수 응답의 역수를 취함으로써 결정된다(밀도≒ 1/f2 이므로). 필요하다면, 다른 유동 특성, 변수, 및/또는 상수도 이러한 결정을 하는데 사용될 수 있다는 것을 이해해야 한다.
단계(2606)에서, 액체 유동 분율이 결정된다. 액체 유동 분율은 전술한 바와 같이, 알려진 가스 밀도에 대한 유동 흐름 밀도의 비교로부터 결정될 수 있다.
도 27은 본원발명의 일 실시예에 따라 유량계의 센서 신호를 처리하기 위한 방법의 순서도(2700)를 도시한다. 단계(2701)에서, 제1 센서 신호 및 제2 센서 신호가 유량계로부터 수신된다. 제1 센서 신호 및 제2 센서 신호는 유량계에서 유동하는 다상 유동 물질에 대한 것이다.
단계(2702)에서, 제1 센서 신호 및 제2 센서 신호 중 하나 또는 모두가 위상 편이된다. 위상 편이 단계는 실질적으로 90도 만큼 센서 신호를 편이시키는 단계를 포함한다. 이는 제1 90도 위상 편이 및 제2 90도 위상 편이를 초래한다.
단계(2703)에서, 유량계의 주파수(또는 주파수 응답)이 계산된다. 일 실시 예에서, 주파수는 전술한 바와 같이 제1 센서 신호 및 제1 90도 위상 편이로부터 생성된다.
단계(2704)에서, 센서 신호들 사이의 위상차가 계산된다. 일 실시예에서, 위상차는 제1 센서 신호, 제1 90도 위상 편이, 제2 센서 신호, 및 제2 90도 위상 편이로부터 결정된다(도 8 및 이에 관한 설명 참조). 다른 실시예에서, 위상차는 제1 센서 신호, 제1 90도 위상 편이, 제2 센서 신호로부터 결정된다(도 12 및 이에 관한 설명 참조).
단계(2705)에서, 질량 유량, 밀도, 또는 용적 유량 중 하나 이상이 계산된다. 전술한 바와 같이, 질량 유량, 밀도, 및/또는 용적 유량은 빠른 주파수 및/또는 빠른 위상 결정을 이용하여 다상 유동 물질에 대해 정확하게 계산된다.
단계(2706)에서, 액체 캐리오버(carryover) 또는 가스 캐리언더(carryunder) 중 하나 또는 모두가 계산된다. 액체 캐리오버 및/또는 가스 캐리언더는 전술한 바와 같이 빠른 주파수 및/또는 빠른 위상 결정을 이용하여 다상 유동 물질에 대해 정확하게 계산된다. 가스 캐리언더 및/또는 액체 캐리오버는 예를 들어 전술한 보이드 분율 방법을 사용함으로써 결정될 수 있다(도 16, 17, 및 19와 이에 관한 설명 참조). 대안적으로, 가스 캐리언더 및/또는 액체 캐리오버는 예를 들어 전술한 질량 분율 방법을 사용함으로써 결정될 수 있다(도 21-23 및 이에 관한 설명 참조). 또 다른 대안으로서, 가스 캐리언더 및/또는 액체 캐리오버는 예를 들어 전술한 액체 유동 분율을 사용하여 결정될 수 있다(도 24 및 26과 이에 관한 설명 참조). 더욱이, 이러한 방법들은 조합하여 사용될 수 있다. 이후 가스와 액체의 백 분율(즉, 가스의 보이드 분율 및 액체 분율)은 다상 유동 물질의 모든 성분의 양을 정확하고 완전하게 정하는데 사용될 수 있다.
본원발명에 따른 방법 및 계측 전자장치는 필요하다면 몇 가지 장점을 얻기 위하여 어떠한 실시예에 따라서도 실행될 수 있다. 본원발명은 다상 유동 물질에 대한 질량 유량, 밀도, 및 용적 유량 중 하나 이상을 계산할 수 있다. 본원발명은 다상 유동 물질에 대한 유동 특성을 신속하게 계산할 수 있다. 본원발명은 보다 높은 정확성과 신뢰성을 갖는 유동 특성의 계산을 제공할 수 있다. 본원발명은 종래기술에 비해 처리 시간을 덜 소요하면서도 더 빠른 유동 특성의 계산을 제공할 수 있다.

Claims (39)

  1. 유량계(5)의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 계측 전자장치(20)로서,
    다상 유동 물질에 대한 제1 센서 신호(210) 및 제2 센서 신호(211)를 수신하기 위한 인터페이스(201); 및
    상기 인터페이스(201)와 통신하는 처리 시스템(203); 을 포함하고,
    상기 처리 시스템이 상기 인터페이스(201)로부터 제1 센서 신호(210) 및 제2 센서 신호(211)를 수신하고, 제1 센서 신호(210)로부터 제1 90도 위상 편이(213)를 생성하고 제2 센서 신호(211)로부터 제2 90도 위상 편이(214)를 생성하고, 제1 90도 위상 편이(213) 또는 제2 90도 위상 편이(214) 중 하나를 사용하여 주파수(221)를 계산하고, 제1 90도 위상 편이(213) 및 제2 90도 위상 편이(214) 중 하나 이상을 사용하여 위상차(220)를 계산하고, 다상 유동 물질에 대한 질량 유량(223), 밀도(224), 또는 용적 유량(225) 중 하나 또는 둘 이상을 계산하도록 구성되는,
    유량계(5)의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 계측 전자장치(20).
  2. 제1항에 있어서,
    상기 인터페이스가 센서 신호를 디지털화하도록 구성된 디지타이저를 포함하 는,
    유량계(5)의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 계측 전자장치(20).
  3. 제1항에 있어서,
    상기 처리 시스템(203)이 다상 유동 물질 내의 액체 캐리오버(1426) 또는 가스 캐리언더(1425) 중 하나 이상을 결정하도록 추가적으로 구성되는,
    유량계(5)의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 계측 전자장치(20).
  4. 제1항에 있어서,
    상기 생성이 제1 90도 위상 편이(213) 및 제2 90도 위상 편이(214)를 생성하기 위하여 힐버트 변환을 이용하는 것을 포함하는,
    유량계(5)의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 계측 전자장치(20).
  5. 제1항에 있어서,
    상기 주파수(221) 계산이 제1 센서 신호(210) 및 제1 90도 위상 편이(213)로부터 주파수(221)를 계산하는 것을 포함하는,
    유량계(5)의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 계측 전자장치(20).
  6. 제1항에 있어서,
    상기 위상차(220) 계산이 제1 센서 신호(210), 제1 90도 위상 편이(213), 및 제2 센서 신호(211)로부터 위상차(220)를 계산하는 것을 포함하는,
    유량계(5)의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 계측 전자장치(20).
  7. 제1항에 있어서,
    상기 위상차(220) 계산이 제1 센서 신호(210), 제1 90도 위상 편이(213), 제2 센서 신호(211), 및 제2 90도 위상 편이(214)로부터 위상차(220)를 계산하는 것을 포함하는,
    유량계(5)의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 계측 전자장치(20).
  8. 제1항에 있어서,
    상기 처리 시스템(203)이 주파수(221)를 적어도 가스 주파수 성분(1412) 및 유체 주파수 성분(1416)으로 분리하고, 가스 주파수 성분(1412) 및 유체 주파수 성분(1416) 중 하나 이상 및 주파수 응답(1410)으로부터 액체 분율(1427) 또는 가스의 보이드 분율(1418) 중 하나 이상을 결정하고, 가스의 보이드 분율(1418)을 이용하여 가스 유동 성분의 가스상 밀도(1421) 또는 다상 유동 물질의 액체 유동 성분의 액체상 밀도(1422) 중 하나 이상을 결정하고, 그리고 액체 분율(1427) 또는 가스의 보이드 분율(1418) 중 하나 이상을 이용하여 액체 캐리오버(1426) 또는 다상 유동 물질의 가스 캐리언더(1425) 중 하나 이상을 결정하도록 추가적으로 구성되는,
    유량계(5)의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 계측 전자장치(20).
  9. 제8항에 있어서,
    상기 분리가 주파수 응답(1410)을 가스 주파수 성분(1412) 및 유체 주파수 성분(1416) 중 하나를 실질적으로 차단하는 하나 또는 둘 이상의 필터로 처리하는 것을 포함하는,
    유량계(5)의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 계측 전자장 치(20).
  10. 제8항에 있어서,
    상기 분리가 가스 주파수 성분(1412)을 실질적으로 차단하고 유체 주파수 성분(1416)을 실질적으로 통과시키는 제1 필터로 주파수 응답(1410)을 필터링하는 것과, 유체 주파수 성분(1416)을 실질적으로 차단하고 가스 주파수 성분(1412)을 실질적으로 통과시키는 제2 필터로 주파수 응답(1410)을 필터링하는 것을 포함하고, 여기서 제1 필터가 유체 주파수 성분(1416)을 출력하고 제2 필터가 가스 주파수 성분(1412)을 출력하는,
    유량계(5)의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 계측 전자장치(20).
  11. 제8항에 있어서,
    액체 분율(1427) 또는 가스의 보이드 분율(1418) 중 하나 이상을 결정하는 것이, 주파수 응답(1410)으로부터 전체 밀도(1420)를 계산하는 것, 유체 주파수 성분(1416)으로부터 유체 성분 밀도(1422)를 계산하는 것, 가스 주파수 성분(1412)으로부터 가스 성분 밀도(1421)를 계산하는 것, 및 유체 성분 밀도(1422)에서 전체 밀도(1420)를 뺀 값을 유체 성분 밀도(1422)에서 가스 성분 밀도(1421)를 뺀 값으 로 나눈 비율로서 가스의 보이드 분율(1420)을 계산하는 것을 포함하는,
    유량계(5)의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 계측 전자장치(20).
  12. 제1항에 있어서,
    상기 처리 시스템(203)이 주파수 응답(1410)을 적어도 가스 주파수 성분(1412) 및 유체 주파수 성분(1416)으로 분리하고, 주파수 응답(1410)으로부터 전체 밀도(1420)를 결정하고, 가스 주파수 성분(1412)으로부터 가스 밀도(1421)를 결정하고, 가스 주파수 성분(1412) 및 유체 주파수 성분(1416) 중 하나 이상과 주파수 응답(1410)으로부터 가스의 보이드 분율(1418)을 결정하고, 그리고 가스의 보이드 분율(1418)에 가스의 밀도(1421)를 전체 밀도(1420)로 나눈 비율을 곱한 값으로부터 질량 분율(1419)을 결정하도록 추가적으로 구성되는,
    유량계(5)의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 계측 전자장치(20).
  13. 제12항에 있어서,
    상기 처리 시스템(203)이 주파수 응답(1410)으로부터 유동 물질의 질량 유량(223)을 결정하고 질량 분율(1419) 및 질량 유량(223)을 이용하여 제1 유동 성분 질량과 제2 유동 성분 질량 중 하나 이상을 결정하도록 추가적으로 구성되는,
    유량계(5)의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 계측 전자장치(20).
  14. 제1항에 있어서,
    상기 처리 시스템(203)이 제곱된 주파수 응답을 생성하도록 주파수 응답(1410)을 제곱하고, 실질적으로 순간적인 유동 흐름 밀도를 생성하도록 제곱된 주파수 응답의 역수를 취하고, 실질적으로 순간적인 유동 흐름 밀도를 가스 유동 물질의 가스 유동 분율(1428)을 나타내는 소정의 가스 밀도(1421) 및 액체 유동 분율(1427)을 나타내는 소정의 액체 밀도(1422) 중 하나 이상과 비교하고, 이러한 비교로부터 액체 유동 분율(1427) 또는 가스 유동 분율(1428) 중 하나 이상을 결정하도록 추가적으로 구성되는,
    유량계(5)의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 계측 전자장치(20).
  15. 유량계(5) 내의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 방법으로서,
    다상 유동 물질에 대한 제1 센서 신호(210) 및 제2 센서 신호(211)를 수신하 는 단계;
    제1 센서 신호(210)로부터 제1 90도 위상 편이(213)를 생성하고 제2 센서 신호(211)로부터 제2 90도 위상 편이(214)를 생성하는 단계;
    제1 90도 위상 편이(213) 및 제2 90도 위상 편이(214) 중 하나를 이용하여 주파수(221)를 계산하는 단계;
    제1 90도 위상 편이(213) 및 제2 90도 위상 편이(214) 중 하나를 이용하여 위상차(220)를 계산하는 단계; 및
    다상 유동 물질에 대한 질량 유량(223), 밀도(224), 또는 용적 유량(225) 중 하나 또는 둘 이상을 계산하는 단계; 를 포함하는,
    유량계(5) 내의 다상 유동 물질에 대한 센서 신호 처리 방법.
  16. 제15항에 있어서,
    다상 유동 물질 내의 액체 캐리오버(1426) 또는 가스 캐리언더(1425) 중 하나 이상을 결정하는 단계를 포함하는,
    유량계(5) 내의 다상 유동 물질에 대한 센서 신호 처리 방법.
  17. 제15항에 있어서,
    상기 생성 단계가 제1 90도 위상 편이(213) 및 제2 90도 위상 편이(214)를 생성하기 위하여 힐버트 변환을 이용하는 단계를 포함하는,
    유량계(5) 내의 다상 유동 물질에 대한 센서 신호 처리 방법.
  18. 제15항에 있어서,
    상기 주파수(221) 계산 단계가 제1 센서 신호(210) 및 제1 90도 위상 편이(213)로부터 주파수(221)를 계산하는 단계를 포함하는,
    유량계(5) 내의 다상 유동 물질에 대한 센서 신호 처리 방법.
  19. 제15항에 있어서,
    상기 위상차(220) 계산 단계가 제1 센서 신호(210), 제1 90도 위상 편이(213), 및 제2 센서 신호(211)로부터 위상차(220)를 계산하는 단계를 포함하는,
    유량계(5) 내의 다상 유동 물질에 대한 센서 신호 처리 방법.
  20. 제15항에 있어서,
    상기 위상차(220) 계산 단계가 제1 센서 신호(210), 제1 90도 위상 편이(213), 제2 센서 신호(211), 및 제2 90도 위상 편이(214)로부터 위상차(220)를 계산하는 단계를 포함하는,
    유량계(5) 내의 다상 유동 물질에 대한 센서 신호 처리 방법.
  21. 제15항에 있어서,
    주파수(221)를 적어도 가스 주파수 성분(1412) 및 유체 주파수 성분(1416)으로 분리하는 단계;
    가스 주파수 성분(1412) 및 유체 주파수 성분(1416) 중 하나 이상 및 주파수 응답(1410)으로부터 액체 분율(1427) 또는 가스의 보이드 분율(1418) 중 하나 이상을 결정하는 단계;
    가스의 보이드 분율(1418)을 이용하여 가스 유동 성분의 가스상 밀도(1421) 또는 다상 유동 물질의 액체 유동 성분의 액체상 밀도(1422) 중 하나 이상을 결정하는 단계; 및
    액체 분율(1427) 또는 가스의 보이드 분율(1418) 중 하나 이상을 이용하여 액체 캐리오버(1426) 또는 다상 유동 물질의 가스 캐리언더(1425) 중 하나 이상을 결정하는 단계; 를 더 포함하는,
    유량계(5) 내의 다상 유동 물질에 대한 센서 신호 처리 방법.
  22. 제21항에 있어서,
    상기 분리 단계가 주파수 응답(1410)을 가스 주파수 성분(1412) 및 유체 주 파수 성분(1416) 중 하나를 실질적으로 차단하는 하나 또는 둘 이상의 필터로 처리하는 단계를 포함하는,
    유량계(5) 내의 다상 유동 물질에 대한 센서 신호 처리 방법.
  23. 제21항에 있어서,
    상기 분리 단계가 가스 주파수 성분(1412)을 실질적으로 차단하고 유체 주파수 성분(1416)을 실질적으로 통과시키는 제1 필터로 주파수 응답(1410)을 필터링하는 단계; 및
    유체 주파수 성분(1416)을 실질적으로 차단하고 가스 주파수 성분(1412)을 실질적으로 통과시키는 제2 필터로 주파수 응답(1410)을 필터링하는 단계; 를 포함하고,
    상기 제1 필터가 유체 주파수 성분(1416)을 출력하고 상기 제2 필터가 가스 주파수 성분(1412)을 출력하는,
    유량계(5) 내의 다상 유동 물질에 대한 센서 신호 처리 방법.
  24. 제21항에 있어서,
    액체 분율(1427) 또는 가스의 보이드 분율(1418) 중 하나 이상을 결정하는 단계가,
    주파수 응답(1410)으로부터 전체 밀도(1420)를 계산하는 단계;
    유체 주파수 성분(1416)으로부터 유체 성분 밀도(1422)를 계산하는 단계;
    가스 주파수 성분(1412)으로부터 가스 성분 밀도(1421)를 계산하는 단계; 및
    유체 성분 밀도(1422)에서 전체 밀도(1420)를 뺀 값을 유체 성분 밀도(1422)에서 가스 성분 밀도(1421)를 뺀 값으로 나눈 비율로서 가스의 보이드 분율(1420)을 계산하는 단계; 를 포함하는,
    유량계(5) 내의 다상 유동 물질에 대한 센서 신호 처리 방법.
  25. 제15항에 있어서,
    주파수 응답(1410)을 적어도 가스 주파수 성분(1412) 및 유체 주파수 성분(1416)으로 분리하는 단계;
    주파수 응답(1410)으로부터 전체 밀도(1420)를 결정하는 단계;
    가스 주파수 성분(1412)으로부터 가스 밀도(1421)를 결정하는 단계;
    가스 주파수 성분(1412) 및 유체 주파수 성분(1416) 중 하나 이상과 주파수 응답(1410)으로부터 가스의 보이드 분율(1418)을 결정하는 단계; 및
    가스의 보이드 분율(1418)에 가스의 밀도(1421)를 전체 밀도(1420)로 나눈 비율을 곱한 값으로부터 질량 분율(1419)을 결정하는 단계; 를 더 포함하는,
    유량계(5) 내의 다상 유동 물질에 대한 센서 신호 처리 방법.
  26. 제25항에 있어서,
    주파수 응답(1410)으로부터 유동 물질의 질량 유량(223)을 결정하는 단계; 및
    질량 분율(1419) 및 질량 유량(223)을 이용하여 제1 유동 성분 질량과 제2 유동 성분 질량 중 하나 이상을 결정하는 단계; 를 더 포함하는,
    유량계(5) 내의 다상 유동 물질에 대한 센서 신호 처리 방법.
  27. 제15항에 있어서,
    제곱된 주파수 응답을 생성하도록 주파수 응답(1410)을 제곱하는 단계;
    실질적으로 순간적인 유동 흐름 밀도를 생성하도록 제곱된 주파수 응답의 역수를 취하는 단계;
    실질적으로 순간적인 유동 흐름 밀도를 가스 유동 물질의 가스 유동 분율(1428)을 나타내는 소정의 가스 밀도(1421) 및 액체 유동 분율(1427)을 나타내는 소정의 액체 밀도(1422) 중 하나 이상과 비교하는 단계; 및
    상기 비교로부터 액체 유동 분율(1427) 또는 가스 유동 분율(1428) 중 하나 이상을 결정하는 단계; 를 더 포함하는,
    유량계(5) 내의 다상 유동 물질에 대한 센서 신호 처리 방법.
  28. 유량계(5) 내의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 방법으로서,
    다상 유동 물질에 대한 제1 센서 신호(210) 및 제2 센서 신호(211)를 수신하는 단계;
    제1 센서 신호(210)로부터 제1 90도 위상 편이(213)를 생성하고 제2 센서 신호(211)로부터 제2 90도 위상 편이(214)를 생성하는 단계;
    제1 90도 위상 편이(213) 또는 제2 90도 위상 편이(214) 중 하나를 이용하여 주파수(221)를 계산하는 단계;
    제1 90도 위상 편이(213) 및 제2 90도 위상 편이(214) 중 하나를 이용하여 위상차(220)를 계산하는 단계;
    다상 유동 물질에 대한 질량 유량(223), 밀도(224), 또는 용적 유량(225) 중 하나 또는 둘 이상을 계산하는 단계; 및
    다상 유동 물질 내의 가스 캐리언더(1425) 또는 액체 캐리오버(1426) 중 하나 이상을 계산하는 단계; 를 포함하는,
    유량계(5) 내의 다상 유동 물질에 대한 센서 신호 처리 방법.
  29. 제28항에 있어서,
    상기 생성 단계가 제1 90도 위상 편이(213) 및 제2 90도 위상 편이(214)를 생성하기 위하여 힐버트 변환을 이용하는 단계를 포함하는,
    유량계(5) 내의 다상 유동 물질에 대한 센서 신호 처리 방법.
  30. 제28항에 있어서,
    상기 주파수(221) 계산 단계가 제1 센서 신호(210) 및 제1 90도 위상 편이(213)로부터 주파수(221)를 계산하는 단계를 포함하는,
    유량계(5) 내의 다상 유동 물질에 대한 센서 신호 처리 방법.
  31. 제28항에 있어서,
    상기 위상차(220) 계산 단계가 제1 센서 신호(210), 제1 90도 위상 편이(213), 및 제2 센서 신호(211)로부터 위상차(220)를 계산하는 단계를 포함하는,
    유량계(5) 내의 다상 유동 물질에 대한 센서 신호 처리 방법.
  32. 제28항에 있어서,
    상기 위상차(220) 계산 단계가 제1 센서 신호(210), 제1 90도 위상 편이(213), 제2 센서 신호(211), 및 제2 90도 위상 편이(214)로부터 위상차(220)를 계산하는 단계를 포함하는,
    유량계(5) 내의 다상 유동 물질에 대한 센서 신호 처리 방법.
  33. 제28항에 있어서,
    주파수(221)를 적어도 가스 주파수 성분(1412) 및 유체 주파수 성분(1416)으로 분리하는 단계;
    가스 주파수 성분(1412) 및 유체 주파수 성분(1416) 중 하나 이상 및 주파수 응답(1410)으로부터 액체 분율(1427) 또는 가스의 보이드 분율(1418) 중 하나 이상을 결정하는 단계;
    가스의 보이드 분율(1418)을 이용하여 가스 유동 성분의 가스상 밀도(1421) 또는 다상 유동 물질의 액체 유동 성분의 액체상 밀도(1422) 중 하나 이상을 결정하는 단계; 및
    액체 분율(1427) 또는 가스의 보이드 분율(1418) 중 하나 이상을 이용하여 액체 캐리오버(1426) 또는 다상 유동 물질의 가스 캐리언더(1425) 중 하나 이상을 결정하는 단계; 를 더 포함하는,
    유량계(5) 내의 다상 유동 물질에 대한 센서 신호 처리 방법.
  34. 제33항에 있어서,
    상기 분리 단계가 주파수 응답(1410)을 가스 주파수 성분(1412) 및 유체 주 파수 성분(1416) 중 하나를 실질적으로 차단하는 하나 또는 둘 이상의 필터로 처리하는 단계를 포함하는,
    유량계(5) 내의 다상 유동 물질에 대한 센서 신호 처리 방법.
  35. 제33항에 있어서,
    상기 분리 단계가 가스 주파수 성분(1412)을 실질적으로 차단하고 유체 주파수 성분(1416)을 실질적으로 통과시키는 제1 필터로 주파수 응답(1410)을 필터링하는 단계; 및
    유체 주파수 성분(1416)을 실질적으로 차단하고 가스 주파수 성분(1412)을 실질적으로 통과시키는 제2 필터로 주파수 응답(1410)을 필터링하는 단계; 를 포함하고,
    상기 제1 필터가 유체 주파수 성분(1416)을 출력하고 상기 제2 필터가 가스 주파수 성분(1412)을 출력하는,
    유량계(5) 내의 다상 유동 물질에 대한 센서 신호 처리 방법.
  36. 제33항에 있어서,
    액체 분율(1427) 또는 가스의 보이드 분율(1418) 중 하나 이상을 결정하는 단계가,
    주파수 응답(1410)으로부터 전체 밀도(1420)를 계산하는 단계;
    유체 주파수 성분(1416)으로부터 유체 성분 밀도(1422)를 계산하는 단계;
    가스 주파수 성분(1412)으로부터 가스 성분 밀도(1421)를 계산하는 단계; 및
    유체 성분 밀도(1422)에서 전체 밀도(1420)를 뺀 값을 유체 성분 밀도(1422)에서 가스 성분 밀도(1421)를 뺀 값으로 나눈 비율로서 가스의 보이드 분율(1420)을 계산하는 단계; 를 포함하는,
    유량계(5) 내의 다상 유동 물질에 대한 센서 신호 처리 방법.
  37. 제28항에 있어서,
    주파수 응답(1410)을 적어도 가스 주파수 성분(1412) 및 유체 주파수 성분(1416)으로 분리하는 단계;
    주파수 응답(1410)으로부터 전체 밀도(1420)를 결정하는 단계;
    가스 주파수 성분(1412)으로부터 가스 밀도(1421)를 결정하는 단계;
    가스 주파수 성분(1412) 및 유체 주파수 성분(1416) 중 하나 이상과 주파수 응답(1410)으로부터 가스의 보이드 분율(1418)을 결정하는 단계; 및
    가스의 보이드 분율(1418)에 가스의 밀도(1421)를 전체 밀도(1420)로 나눈 비율을 곱한 값으로부터 질량 분율(1419)을 결정하는 단계; 를 더 포함하는,
    유량계(5) 내의 다상 유동 물질에 대한 센서 신호 처리 방법.
  38. 제37항에 있어서,
    주파수 응답(1410)으로부터 유동 물질의 질량 유량(223)을 결정하는 단계; 및
    질량 분율(1419) 및 질량 유량(223)을 이용하여 제1 유동 성분 질량과 제2 유동 성분 질량 중 하나 이상을 결정하는 단계; 를 더 포함하는,
    유량계(5) 내의 다상 유동 물질에 대한 센서 신호 처리 방법.
  39. 제28항에 있어서,
    제곱된 주파수 응답을 생성하도록 주파수 응답(1410)을 제곱하는 단계;
    실질적으로 순간적인 유동 흐름 밀도를 생성하도록 제곱된 주파수 응답의 역수를 취하는 단계;
    실질적으로 순간적인 유동 흐름 밀도를 가스 유동 물질의 가스 유동 분율(1428)을 나타내는 소정의 가스 밀도(1421) 및 액체 유동 분율(1427)을 나타내는 소정의 액체 밀도(1422) 중 하나 이상과 비교하는 단계; 및
    상기 비교로부터 액체 유동 분율(1427) 또는 가스 유동 분율(1428) 중 하나 이상을 결정하는 단계; 를 더 포함하는,
    유량계(5) 내의 다상 유동 물질에 대한 센서 신호 처리 방법.
KR1020087006616A 2005-08-18 2006-08-15 유량계의 다상 유동 물질에 대한 센서 신호를 처리하기위한 방법 및 계측 전자장치 KR20080039498A (ko)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US70927105P 2005-08-18 2005-08-18
US60/709,271 2005-08-18

Related Child Applications (1)

Application Number Title Priority Date Filing Date
KR1020117005940A Division KR20110032009A (ko) 2005-08-18 2006-08-15 유량계의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 방법 및 계측 전자장치

Publications (1)

Publication Number Publication Date
KR20080039498A true KR20080039498A (ko) 2008-05-07

Family

ID=37527015

Family Applications (3)

Application Number Title Priority Date Filing Date
KR1020117005940A KR20110032009A (ko) 2005-08-18 2006-08-15 유량계의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 방법 및 계측 전자장치
KR1020127033903A KR101410003B1 (ko) 2005-08-18 2006-08-15 유량계의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 방법 및 계측 전자장치
KR1020087006616A KR20080039498A (ko) 2005-08-18 2006-08-15 유량계의 다상 유동 물질에 대한 센서 신호를 처리하기위한 방법 및 계측 전자장치

Family Applications Before (2)

Application Number Title Priority Date Filing Date
KR1020117005940A KR20110032009A (ko) 2005-08-18 2006-08-15 유량계의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 방법 및 계측 전자장치
KR1020127033903A KR101410003B1 (ko) 2005-08-18 2006-08-15 유량계의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 방법 및 계측 전자장치

Country Status (12)

Country Link
US (1) US8781759B2 (ko)
EP (1) EP1938051B9 (ko)
JP (1) JP4966306B2 (ko)
KR (3) KR20110032009A (ko)
CN (1) CN101268341B (ko)
AR (1) AR054920A1 (ko)
AU (1) AU2006279675B2 (ko)
BR (1) BRPI0614823B1 (ko)
CA (1) CA2619689C (ko)
HK (1) HK1124658A1 (ko)
RU (1) RU2371680C1 (ko)
WO (1) WO2007022118A1 (ko)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10508997B1 (en) 2018-05-28 2019-12-17 Korea Institute Of Geoscience And Mineral Resources Method of analyzing fines migration in multiphase flow in sediment layer using X-ray computed tomography image

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2626245C (en) * 2005-10-18 2013-05-14 Micro Motion, Inc. Meter electronics and methods for determining a phase difference between a first sensor signal and a second sensor signal of a flow meter
US8751171B2 (en) * 2007-03-07 2014-06-10 Invensys Systems, Inc. Coriolis frequency tracking
EP2179254B1 (en) * 2007-07-30 2020-03-18 Micro Motion, Inc. Flow meter system and method for measuring flow characteristics of a three phase flow
MX2010011639A (es) * 2008-05-01 2010-11-30 Micro Motion Inc Flujometro vibratorio de muy alta frecuencia.
CN102625905B (zh) * 2009-05-04 2013-10-30 琼脂有限公司 多相流体测量装置和方法
US20140136125A1 (en) * 2010-05-04 2014-05-15 Agar Corporation Ltd. System and method for multi-phase fluid measurement
JP5510747B2 (ja) * 2010-11-05 2014-06-04 横河電機株式会社 コリオリ質量流量計
JP6437106B2 (ja) * 2014-09-18 2018-12-12 マイクロ モーション インコーポレイテッド 密度差を判定するための方法及び装置
DE102015102834A1 (de) * 2015-02-27 2016-09-01 Endress + Hauser Gmbh + Co. Kg Vibronischer Sensor
US9863798B2 (en) 2015-02-27 2018-01-09 Schneider Electric Systems Usa, Inc. Systems and methods for multiphase flow metering accounting for dissolved gas
US10876874B2 (en) * 2016-01-13 2020-12-29 Micro Motion, Inc. Multi-phase coriolis measurement device and method
EP3638989B1 (en) * 2017-06-14 2023-03-01 Micro Motion, Inc. A notch filter in a vibratory flow meter
DE102017115251A1 (de) 2017-07-07 2019-01-10 Endress+Hauser Flowtec Ag Die vorliegende Erfindung betrifft einen Messaufnehmer zum Bestimmen des Massedurchflusses einer Flüssigkeit
US11196359B2 (en) * 2018-05-11 2021-12-07 Micro Motion Inc. Supervisory monitor for energy measurement
CN111060168A (zh) * 2019-12-26 2020-04-24 深圳市佳运通电子有限公司 一种流量信号采样方法和装置

Family Cites Families (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2865201A (en) * 1954-08-26 1958-12-23 Roth Wilfred Gyroscopic mass flowmeter
GB1280997A (en) * 1968-10-29 1972-07-12 Solartron Electronic Group Improvements in or relating to fluid density transducers
US4812739A (en) * 1986-09-15 1989-03-14 Swanson Claude V Apparatus and method for using microwave radiation to measure water content of a fluid
US5069074A (en) * 1987-07-22 1991-12-03 Exac Corporation Apparatus and method for measuring the mass flow rate of material flowing through at least one vibrating conduit
US4876879A (en) * 1988-08-23 1989-10-31 Ruesch James R Apparatus and methods for measuring the density of an unknown fluid using a Coriolis meter
US4852395A (en) * 1988-12-08 1989-08-01 Atlantic Richfield Company Three phase fluid flow measuring system
US5224372A (en) * 1990-05-14 1993-07-06 Atlantic Richfield Company Multi-phase fluid flow measurement
US5295084A (en) * 1991-10-08 1994-03-15 Micromotion, Inc. Vibrating tube densimeter
JP3200827B2 (ja) * 1993-12-24 2001-08-20 横河電機株式会社 コリオリ質量流量計
JP3219122B2 (ja) * 1994-07-11 2001-10-15 横河電機株式会社 コリオリ質量流量計
US5594180A (en) * 1994-08-12 1997-01-14 Micro Motion, Inc. Method and apparatus for fault detection and correction in Coriolis effect mass flowmeters
EP0698783A1 (de) * 1994-08-16 1996-02-28 Endress + Hauser Flowtec AG Auswerte-Elektronik eines Coriolis-Massedurchflussaufnehmers
US5524475A (en) * 1994-11-10 1996-06-11 Atlantic Richfield Company Measuring vibration of a fluid stream to determine gas fraction
US6268727B1 (en) * 1997-06-24 2001-07-31 Southwest Research Institute Measurement of flow fractions flow velocities and flow rates of a multiphase fluid using ESR sensing
US7124646B2 (en) * 1997-11-26 2006-10-24 Invensys Systems, Inc. Correcting for two-phase flow in a digital flowmeter
ID29450A (id) 1998-12-08 2001-08-30 Micro Motion Inc Pengontrol aliran massa coriolis
US6748813B1 (en) * 1998-12-08 2004-06-15 Emerson Electric Company Coriolis mass flow controller
US6526839B1 (en) * 1998-12-08 2003-03-04 Emerson Electric Co. Coriolis mass flow controller and capacitive pick off sensor
US6505131B1 (en) 1999-06-28 2003-01-07 Micro Motion, Inc. Multi-rate digital signal processor for signals from pick-offs on a vibrating conduit
US6412354B1 (en) 1999-12-16 2002-07-02 Halliburton Energy Services, Inc. Vibrational forced mode fluid property monitor and method
GB0221782D0 (en) * 2002-09-19 2002-10-30 Univ Sussex Methods of measuring two-phase fluid flow using single-phase flowmeters
US6847898B1 (en) * 2003-08-21 2005-01-25 Appleton Papers Inc. Real time determination of gas solubility and related parameters in manufacturing processes
NO320172B1 (no) * 2004-02-27 2005-11-07 Roxar Flow Measurement As Stromningsmaler og fremgangsmate for maling av individuelle mengder av gass, hydrokarbonvaeske og vann i en fluidblanding
US6912918B1 (en) * 2004-03-10 2005-07-05 General Electric Company Mass flow sensor and methods of determining mass flow of a fluid
RU2371678C2 (ru) 2004-12-29 2009-10-27 Майкро Моушн, Инк. Высокоскоростная оценка частоты и фазы расходомеров
CA2608838C (en) * 2005-05-20 2013-05-14 Micro Motion, Inc. Meter electronics and methods for rapidly determining a mass fraction of a multi-phase fluid from a coriolis flow meter signal
AR054760A1 (es) * 2005-05-20 2007-07-18 Micro Motion Inc Equipo de medicion electronico y metodos para determinar una fraccion de vacio de gas
WO2008011587A2 (en) * 2006-07-21 2008-01-24 Invensys Systems, Inc. Multi-phase coriolis flowmeter

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10508997B1 (en) 2018-05-28 2019-12-17 Korea Institute Of Geoscience And Mineral Resources Method of analyzing fines migration in multiphase flow in sediment layer using X-ray computed tomography image

Also Published As

Publication number Publication date
KR20110032009A (ko) 2011-03-29
JP2009505100A (ja) 2009-02-05
EP1938051B9 (en) 2023-05-31
JP4966306B2 (ja) 2012-07-04
EP1938051B1 (en) 2022-12-21
CN101268341A (zh) 2008-09-17
RU2371680C1 (ru) 2009-10-27
AU2006279675B2 (en) 2010-12-16
BRPI0614823B1 (pt) 2018-02-14
BRPI0614823A2 (pt) 2011-04-19
CN101268341B (zh) 2012-07-11
AU2006279675A1 (en) 2007-02-22
HK1124658A1 (en) 2009-07-17
KR101410003B1 (ko) 2014-06-20
KR20130014625A (ko) 2013-02-07
US20080243400A1 (en) 2008-10-02
CA2619689A1 (en) 2007-02-22
US8781759B2 (en) 2014-07-15
EP1938051A1 (en) 2008-07-02
AR054920A1 (es) 2007-07-25
WO2007022118A1 (en) 2007-02-22
CA2619689C (en) 2013-08-13

Similar Documents

Publication Publication Date Title
KR101206381B1 (ko) 가스 유동 물질 내의 액체 유동 분율을 결정하기 위한 방법 및 계측 전자장치
KR101410003B1 (ko) 유량계의 다상 유동 물질에 대한 센서 신호를 처리하기 위한 방법 및 계측 전자장치
KR101132778B1 (ko) 코리올리 유량계 신호로부터 다상 유체의 질량 분율을신속히 결정하기 위한 방법 및 계측 전자장치
KR20110039399A (ko) 유량계의 제1 센서 신호와 제2 센서 신호 사이의 위상차를 결정하기 위한 방법 및 계측 전자장치
KR101153466B1 (ko) 가스의 보이드 분율을 결정하기 위한 방법 및 계측전자장치
MX2008002110A (es) Electronica de medicion y metodos para procesar señales de sensor para un material del flujo multi-fase en un flujometro

Legal Events

Date Code Title Description
A201 Request for examination
E902 Notification of reason for refusal
AMND Amendment
E601 Decision to refuse application
J201 Request for trial against refusal decision
A107 Divisional application of patent
AMND Amendment
E801 Decision on dismissal of amendment
B601 Maintenance of original decision after re-examination before a trial
J301 Trial decision

Free format text: TRIAL DECISION FOR APPEAL AGAINST DECISION TO DECLINE REFUSAL REQUESTED 20110211

Effective date: 20120328

WITB Written withdrawal of application