专利名称：用于确定气体空隙度的电子计量器和方法
技术领域：
本发明涉及用于确定流过流量计的流体材料中气体空隙度的电子计 量器和方法。
背景技术：
采用科里奥利质量流量计来测量流过管线的材料的质量流量、密度与
体积流量和其它信息是公知常识，在1985年1月1日授权给J. E. Smith 等人的美国专利4， 491， 025中以及1982年2月11日授权给J. E. Smith的 Re.31，450中是已知的。该流量计具有一个或多个不同结构的流管。每一 个管道结构可被视为具有一组自然振动模式，该自然振动模式例如包括简 单弯曲、扭转、径向和耦合模式。在一般的科里奥利流量计质量流量测量 应用中，当材料流经该管道时，以一个或多个振动模式激励管道结构，并
且沿着该管道在间隔点处测量该管道的运动。
材料填充系统的振动模式部分地由流管和该流管中材料的组合质量 限定。材料从流量计入口侧上连接的管线流到流量计。然后，将该材料引 入流管或多个流管，并从流量计流出而进入连接到出口侧上的管线中。
驱动器施加力给流管。该力引起该流管振荡。当没有材料流经该流量 计时，沿着流管的所有点以相同的相位振荡。当材料开始流经该流管时， 科里奥利加速度引起沿着该流管的每一点相对于沿着该流管的其它点来 讲具有不同的相位。该流管入口侧上的相位滞后该驱动器，而出口侧上的
相位领先该驱动器。该流管上不同点处放置了传感器，以产生表示不同点 处流管运动的正弦信号。这两个传感器信号之间的相差与流经该流管或多 个流管的材料质量流量成比例。在一种现有技术的方法中，利用离散傅里 叶变换(DFT)或快速傅里叶变换(FFT)来确定传感器信号之间的相差。 利用该相差和流体管组件的振动频率响应来获得质量流量。
在一种现有技术方法中，利用独立参考信号来确定拾取信号频率。例 如通过利用发送到振动驱动器系统的频率。在另一种现有技术方法中，可
以通过中心位于陷波滤波器中的频率来确定拾取传感器产生的振动响应
频率，其中该现有技术流量计试图将陷波滤波器的槽口 (notch)保持在
拾取传感器频率。在静止条件下，该现有技术工作性能相当好，其中流量 计中的流体材料是均匀的并且其中最终的拾取信号频率相对稳定。然而， 当流体材料不均匀例如在两相流体中时，该现有技术的相位测量性能低， 其中在上述两相流体中，流体材料包括液体和固体或者在液体流体材料中 包含气泡。在这种情况下，该现有技术确定的频率会快速波动。在快速且 大幅频率转换条件下，拾取信号可以变化到滤波器带宽外部，导致相位和
频率测量不准确。这在空一满一空定量(batching)中也是一个问题，其 中在交替式空和满条件下反复操作流量计。此外，如果传感器的频率快速 变化，则解调过程不能跟上实际或测量频率的步伐，导致在不准确的频率 下进行解调。应当理解的是，如果确定的频率不准确或不精确，则随后获 得的密度、体积流量等数值也会不准确和不精确。此外，会在随后的流动 特性确定中混合有该误差。
在现有技术中，拾取信号可以数字化以及数字化操控，从而实现陷波 滤波器。该陷波滤波器仅仅接收窄带频率。因此，当目标频率变化时，该 陷波滤波器在一定时间内不能追踪目标信号。通常，数字陷波滤波器追踪 波动目标信号需要l一2秒钟。因为现有技术确定该频率所需的时间，所 以结果是不仅频率和相位确定包含误差，而且误差测量所包含的时距 (time span)超过误差和/或两相流体实际发生期间的时距。这是因为陷波 滤波器方案的响应相对缓慢。
结果是，当流量计中的流体材料为两相流体时，该现有技术流量计不 能精确地、快速地或者令人满意地追踪或确定拾取传感器频率。因此，由 于该现有技术利用确定的拾取频率来获得相差，所以相位确定同样也是缓 慢的并且易出现误差。因此，在相位确定中混合有频率确定中的任何误 差。结果是，在频率确定以及在相位确定中误差增大，导致在确定质量流 量中误差增大。此外，因为确定的频率值用于确定密度值(密度约等于l 除以频率平方)，在密度确定中重复或混合有频率确定误差。对于体积流 量确定也是如此，其中在体积流量确定中体积流量等于质量流量除以密 度。
因此，在两相流体条件下该现有技术的准确度低或者损失了测量能 力。在许多流体应用中，流体材料中会带有空气(或者其它气体)。 一种 示例是油田生产，其中，从油井喷出的原油很可能在原油中混有空气。其 它示例有在食品制造过程中在液体中带有空气。所带的空气通常在流动液 体中以气泡排出。空气气泡会导致流量计测得的质量流量产生误差。非常 希望的是，即使由任意量空气混入到流动液体中，流量计也能够准确地测 量出流动液体的质量流量。非常希望的是，即使混入的空气产生波动，流 量计也能够准确地测量出流动液体的质量流量。
发明内容
通过提供用于确定流体材料中气体空隙度的电子计量器和方法，可以 解决上述和其它问题，并且实现本领域的技术进步。
根据本发明的实施例，提供用于确定流过流量计的流体材料中气体空 隙度的电子计量器。该电子计量器包括用于接收流体材料的频率响应的接 口以及与该接口通信的处理系统。该处理系统配置成从该接口接收频率响 应，将该频率响应分解成至少气体频率分量和流体频率分量，并且从该频
率响应以及一个或多个气体频率分量和流体频率分量中来确定气体空隙度。
根据本发明的实施例，提供一种用于确定流过流量计的流体材料中气 体空隙度的方法。该方法包括接收流体材料的频率响应；将该频率响应 分解成至少气体频率分量和流体频率分量；以及从该频率响应以及一个或 多个气体频率分量和流体频率分量中来确定气体空隙度。
根据本发明的实施例，提供一种用于确定流过流量计的流体材料中气 体空隙度的方法。该方法包括接收流体材料的频率响应；利用陷波滤波 器处理该频率响应，该陷波滤波器基本滤去气体频率分量和流体频率分量 之一；从该频率响应以及一个或多个气体频率分量和流体频率分量中来确
定气体空隙度。
根据本发明的实施例，提供一种用于确定流过流量计的流体材料中气
体空隙度的方法。该方法包括接收流体材料的频率响应以及利用第一滤 波器对频率响应进行滤波，该第一滤波器基本滤去气体频率分量并且基本
通过流体频率分量。该第一滤波器输出该流体频率分量。该方法还包括利
用第二滤波器对频率响应进行滤波，该第二滤波器基本滤去流体频率分量 并且基本通过气体频率分量。该第二滤波器输出该气体频率分量。该方法 还包括从该频率响应以及一个或多个气体频率分量和流体频率分量中来 确定气体空隙度。 发明方面
在电子计量器的一个方面中，气体空隙度包括空气空隙度。 在电子计量器的一个方面中，该接口包括数字化该频率响应的数字转 换器。
在电子计量器的另一个方面中，该处理系统还配置成用于从流量计接 收第一传感器信号和第二传感器信号，从第一传感器信号产生90度相 移，并且利用该90度相移来计算该频率响应。
在电子计量器的另一个方面中，该处理系统还配置成用于从流量计接 收第一传感器信号和第二传感器信号，从第一传感器信号产生90度相
移，并且利用该第一传感器信号和该90度相移来计算该频率响应。
在电子计量器的另一个方面中，该处理系统还配置成利用一个或多个 滤波器处理该频率响应从而分解该频率响应，其中该一个或多个滤波器基 本滤去气体频率分量和流体频率分量之一。
在电子计量器的另一个方面中，该处理系统还配置成利用陷波滤波器
处理该频率响应从而分解该频率响应，该陷波滤波器基本滤去气体频率分
量和流体频率分量之一。
在电子计量器的另一个方面中，该处理系统还配置成利用第一滤波器 对该频率响应进行滤波以及利用第二滤波器对该频率响应进行滤波，从而
分解该频率响应，该第一滤波器基本滤去气体频率分量以及基本通过流体 频率分量，该第二滤波器基本滤去流体频率分量以及基本通过气体频率分 量。该第一滤波器输出该流体频率分量，该第二滤波器输出该气体频率分
在电子计量器的另一个方面中，该处理系统还配置成利用低通滤波 器处理频率响应从而分解频率响应，该低通滤波器的低通截止频率基本高
于流体频率分量，该低通滤波器基本通过流体频率分量并且基本滤去气体 频率分量；以及利用高通滤波器处理频率响应从而分解频率响应，该高通
滤波器的高通截止频率基本低于气体频率分量，该高通滤波器基本通过气 体频率分量并且基本滤去流体频率分量。该低通滤波器输出该流体频率分
量，该高通滤波器输出该气体频率分量。
在电子计量器的另一个方面中，该处理系统还配置成从频率响应计 算整个密度从而确定气体空隙度；从流体频率分量计算流体分量密度；从
气体频率分量计算气体分量密度；以及计算气体空隙度，该气体空隙度等
于该流体分量密度减去该整个密度的差除以该流体分量密度减去该气体
分量密度的差的比率。
在电子计量器的另一个方面中，该流量计包括科里奥利流量计。 在该方法的一个方面中，气体空隙度包括空气空隙度。 在该方法的另一个方面中，接收频率响应的步骤还包括从流量计接
收第一传感器信号和第二传感器信号；从第一传感器信号产生90度相
移；以及利用该90度相移计算该频率响应。
在该方法的另一个方面中，接收频率响应的步骤还包括从流量计接
收第一传感器信号和第二传感器信号；从第一传感器信号产生90度相 移；以及利用该第一传感器信号和该90度相移来计算该频率响应。
在该方法的另一个方面中，该分解步骤包括利用一个或多个滤波器 处理该频率响应从而分解该频率响应，该一个或多个滤波器基本滤去气体
频率分量和流体频率分量之一。
在该方法的另一个方面中，该分解步骤包括利用陷波滤波器处理该 频率响应从而分解该频率响应，该陷波滤波器基本滤去气体频率分量和流 体频率分量之一。
在该方法的另一个方面中，该分解步骤包括利用第一滤波器对该频 率响应进行滤波以及利用第二滤波器对该频率响应进行滤波，从而分解该 频率响应，该第一滤波器基本滤去气体频率分量以及基本通过流体频率分 量，该第二滤波器基本滤去流体频率分量以及基本通过气体频率分量。该 第一滤波器输出该流体频率分量，该第二滤波器输出该气体频率分量。
在该方法的另一个方面中，该分解步骤包括利用低通滤波器处理频 率响应从而分解频率响应，该低通滤波器的低通截止频率基本高于流体频 率分量，其中该低通滤波器基本通过流体频率分量并且基本滤去气体频率
分量；以及利用高通滤波器处理频率响应从而分解频率响应，该高通滤波 器的高通截止频率基本低于气体频率分量，其中该高通滤波器基本通过气 体频率分量并且基本滤去流体频率分量。该低通滤波器输出该流体频率分
量，该高通滤波器输出该气体频率分量。
在该方法的另一个方面中，该确定步骤包括从频率响应计算整个密 度；从流体频率分量计算流体分量密度；从气体频率分量计算气体分量密
度；以及计算气体空隙度，该气体空隙度等于该流体分量密度减去该整个
密度的差除以该流体分量密度减去该气体分量密度的差的比率。 在该方法的另一个方面中，该流量计包括科里奥利流量计。
所有附图中相同的附图标记表示相同的部件。
图1是本发明实施例中的科里奥利流量计； 图2示出了根据本发明实施例的电子计量器；
图3是根据本发明实施例的科里奥利流量计中处理传感器信号的方 法流程图4示出了根据本发明实施例的电子计量器；
图5是根据本发明实施例的科里奥利流量计中处理第一和第二传感
器信号的方法流程图6是根据本发明实施例的处理系统一部分的结构图7示出了根据本发明实施例的希尔伯特变换(Hilbert transform)
块细节；
图8和9是根据本发明实施例的分析块的两个独立分支结构图； 图10是在通常情况下流量计的拾取传感器信号的功率频谱密度图； 图11示出了根据单个相移实施例的希尔伯特变换块； 图12示出了用于单个相移实施例的分析块；
图13示出了与现有技术相比的本发明的传感器处理，其中比较每一
个时差(At)值；
图14示出了根据本发明另一实施例的电子计量器；
图15是响应空气、流体和组合的空气/流体混合物(即包括夹带空气
的流体)的流量计频率图16是根据本发明的实施例，用于确定流过流量计的流体材料中气 体空隙度的方法的流程图。
图17是根据本发明的实施例，用于确定流过流量计的流体材料中气 体空隙度的方法的流程图。
图18是示出根据本发明的实施例的低通和高通滤波器响应的频率曲 线图，其中低通和高通滤波器响应用于分解流体频率分量和气体频率分
图19是根据本发明的实施例，用于确定流过流量计的流体材料中气 体空隙度的方法的流程图。
图20是陷波滤波器频率响应的曲线图。
具体实施例方式
图1_20和下面详细描述的具体示例教导本领域技术人员如何进行
和使用本发明的最佳模式。为了教导有创造性的原理，己简化或省略了一 些传统方面。本领域技术人员可以理解，这些示例的变形落入了本发明的 范围。本领域技术人员可以理解，下面描述的特征可以以各种方式组合， 以形成本发明的多种变形。因此，本发明不限于下面描述的具体示例，而 是通过权利要求书和它们的等价物来限定。
图1示出了科里奥利流量计5，其包括计量组件10和电子计量器20。 计量组件10响应于处理材料(process material)的质量流量和密度。 通过导线100，电子计量器20连接到计量组件10，以在通道26上方提供 密度、质量流量和温度信息，以及与本发明不相关的信息。尽管对于本领 域技术人员来说，本发明可实施为振动管比重计是明显的，其中该振动管 比重计没有科里奥利质量流量计提供的额外的测量能力，但是仍然描述了 科里奧利流量计的结构。
计量组件10包括一对支管150和150'、具有凸缘颈部110和110' 的凸缘103和103'、 一对平行流管130和130'、驱动机构180、温度传 感器190和一对速度传感器170L和170R。流管130和130'具有两个基本 上直的入口支路131和131'和出口支路134和134'，该入口支路131和 131'与出口支路134和134'在流管安装块120和120'处朝彼此会聚。流 管130和130'沿着它们的长度在两个对称位置处弯曲，并且在它们的整
个长度上基本上平行。撑杆140和140'用来限定轴W和W'，每一根流管
绕该轴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'中的单一流 中，并且之后被引入出口端104'，该出口端通过具有螺纹孔102'的凸缘 103'连接到处理管线(未示出)。
选择流管130和130'并且合适地安装到流管安装块120和120'，以 分别绕弯曲轴W—W和F—W'具有基本相同的质量分布、惯性矩(moment of inertia)和杨氏模量(Young's modulus)。这些弯曲轴穿过撑杆140 和140'。因为流管的杨氏模量随温度改变，并且该变化影响流量和密度 的计算，所以将抗温探测器(RTD) 190安装到流管130'上，以持续测量 流管温度。因此，通过流经流管的材料温度来控制该流管温度和给定电流 经过RTD出现的电压。在RTD上出现的温度相关电压通过电子计量器20 以已知的方法使用来补偿流管130和130'的弹性系数变化，该弹性系数 变化由流管温度的任何变化引起。该RTD通过导线195连接到电子计量器 20。
在流量计的第一异相弯曲模式下，通过驱动器180绕它们各自的弯曲 轴W和W'在相反方向上驱动两个流管130和130'。该驱动机构180可包 括任何一种已知的布置，例如磁铁安装到流管130'上且相对的线圈安装 到流管130上，并通过该线圈为振荡两根流管传递交流电流。通过电子计 量器20，合适的驱动信号经由导线185施加到驱动机构180上。电子计 量器20接收导线195上的RTD温度信号，并且分别在导线165L和165R 上出现左和右速度信号。电子计量器20产生出现在导线185上的驱动信 号以驱动部件180和振动管130和130'。电子计量器20处理左和右速度
信号以及RTD信号，以计算流经计量组件10的材料的质量流量和密度。
该信息与其它信息一起通过电子计量器20经由通路26施加到应用装置 29。
图2示出了根据本发明实施例的电子计量器20。该电子计量器20可 包括接口 201和处理系统203。该电子计量器20从计量组件10接收第一 和第二传感器信号，诸如拾取/速度传感器信号。该电子计量器20处理第 一和第二传感器信号，以获得流经计量组件10的流体材料的流动特性。 例如，电子计量器20可从这些传感器信号确定例如一个或多个相差、频 率、时差(At)、密度、质量流量和体积流量。此外，其它流动特性可根 据本发明确定。下面将讨论这些确定。
与现有技术的这些确定相比，相差确定和频率确定大大加快并且更加 精确和可靠。在一个具体实施例中，不需要任何频率参考信号，可直接从 仅仅一个传感器信号的相移获得该相差确定和频率确定。这有利地减少了 用以计算流动特性所需的处理时间。在另一个具体实施例中，从两个传感 器信号的相移获得相差，而仅从一个相移信号获得频率。这样增加了这两 个流动特性的精确度，并且与现有技术相比，可更快地确定这两个流动特 性。
现有技术的频率确定方法通常花费1一2秒来进行。相反，根据本发 明的频率确定可在如50毫秒(ms) —样少的时间内进行。取决于该处理 系统的类型和结构、振动响应的采样速率、滤波器尺寸、抽取
(decimation)速率等等，甚至可预期更快的频率确定。在50ms的频率 确定速率下，根据本发明的电子计量器20可比现有技术快约40倍。
接口 201经由图1的导线100接收来自速度传感器170L和170R之一 的传感器信号。该接口 201可进行任何需要的或期望的信号调节，例如任 何方式的格式化、放大、缓存等。可选择地，在该处理系统203中可执行 一些或全部的信号调节。
此外，接口 201可使电子计量器20和外部装置之间进行通信。接口 201可进行任何方式的电子、光学或无线通信。
在一个具体实施例中的接口 201耦合到数字转换器202，其中该传感 器信号包括模拟传感器信号。数字转换器202采样并数字化模拟传感器信
号并产生数字传感器信号。该数字转换器202还可进行任何需要的抽取， 其中该数字传感器信号可被抽取以减小需要处理的信号量，并减小处理时 间。该抽取将在下面更详细地讨论。
处理系统203管理电子计量器20的操作，并处理来自流体计量组件 IO的流体测量结果。该处理系统203执行一个或多个处理程序，因此处 理这些流体测量结果，以产生一个或多个流动特性。
处理系统203可包括通用计算机、微处理器系统、逻辑电路或一些其 它通用目的和惯用处理装置。处理系统203可分布在多个处理装置之间。 处理系统203可包括任何类型的集成或独立的电子存储介质例如存储系 统204。
处理系统203处理传感器信号210，以从该传感器信号210确定一个 或多个流动特性。例如，该一个或多个流动特性可包括流体材料的相差、 频率、时差(At)、质量流量和/或密度。
在示出的具体实施例中，处理系统203从两个传感器信号210和211 以及单个传感器信号相移213确定流动特性。处理系统203可从两个传感 器信号210和211以及单个传感器信号相移213至少确定相差和频率。因 此，第一或者第二相移传感器信号(例如一个上游或下游拾取信号)可通 过根据本发明的处理系统203而被处理，以确定流体材料的相差、频率、 时差(At)和/或质量流量。
存储系统204可存储流量计参数和数据、软件程序、常数值和变量 值。在一个具体实施例中，存储系统204包括通过处理系统203执行的程 序。在一个具体实施例中，该存储系统204存储相移程序212、相差程序 215、频率程序216、时差(At)程序217和流动特性程序218。
在一个具体实施例中，存储系统204存储用于操作科里奥利流量计5 的变量。 一个具体实施例中的存储系统204存储变量，例如第一传感器信 号210和第二传感器信号211，这些传感器信号从速度/拾取传感器170L 和170R接收。此外，存储系统204可存储产生的90度相移213，从而确 定流动特性。
在一个具体实施例中，存储系统204存储从流体测量结果获得的一个 或多个流动特性。 一个具体实施例中的存储系统204存储流动特性，例如 相差220、频率221、时差(At) 222、质量流量223、密度224和体积 流量225，这些流动特性都是从传感器信号210确定的。
相移程序212对输入信号即传感器信号210执行90度相移。在一个 具体实施例中，该相移程序212执行希尔伯特变换(下面讨论)。
相差程序215使用单个90度相移213来确定相差。还可使用额外的 信息来计算该相差。在一个具体实施例中，从第一传感器信号210、第二 传感器信号211和90度相移213计算该相差。确定的相差可存储在存储 系统204的相差220中。当从90度相移213确定相差时，可比现有技术 更快地计算和获得该相差。这可在具有高流速的流量计应用中或者在多相 流动发生的应用中提供临界差。此外，该相差可独立于传感器信号210 或211的频率来确定。而且，因为该相差独立于频率确定，因此相差中的 误微分量就不包括频率确定的误微分量，即在相差测量中没有混合误差。 因此，在现有技术的相差基础上，相差误差减小。
频率程序216从90度相移213确定频率(例如通过第一传感器信号 210或第二传感器信号211表示的)。该确定的频率可存储在存储系统204 的频率221中。当从单个90度相移213确定频率时，与现有技术相比， 可更快地计算和获得该频率。这可在具有高流速的流量计应用中或者在多 相流动发生的应用中提供临界差。
时差(A t)程序217确定第一传感器信号210和第二传感器信号211 之间的时差(At)。该时差(At)可存储在存储系统204的时差(At) 222中。该时差(At)基本上包括确定的相位除以确定的频率，因此用 于确定质量流量。
流动特性程序218可确定一个或多个流动特性。例如，该流动特性程 序218可使用确定的相差220和确定的频率221，以完成这些附加的流动 特性。可以理解，对于这些确定来说，需要附加信息，例如质量流量或密 度。流动特性程序218可从时差(A t) 222并因此从相差220和频率221 来确定质量流量。确定质量流量的公式在Titlow等人的美国专利 5， 027， 662中给出了，这里通过参考引入。该质量流量与计量组件10中 流体材料的质量流量相关。同样，流动特性程序218还可确定密度224 和/或体积流量225。确定的质量流量、密度和体积流量可分别存储在存
储系统204的质量流量223、密度224和体积225中。此外，通过电子计 量器20可将流动特性传送到外部装置。
图3是根据本发明具体实施例的科里奧利流量计中处理传感器信号 的方法的流程图300。在步骤301中，接收第一和第二传感器信号。第一 传感器信号可以包括上游或下游拾取传感器信号。
步骤302中，可调节这些传感器信号。在一个具体实施例中，该调节 可包括滤波以去除噪声和不期望的信号。在一个具体实施例中，该滤波可 包括带通滤波，该带通滤波中心在期望的科里奥利流量计5的基础频率周 围。此外，也可执行其它调节操作，例如放大、缓存等等。如果传感器信 号包括模拟信号，则该步骤还可包括任何类型的采样、数字化和抽取，执 行这些采样、数字化和抽取以产生数字传感器信号。
在步骤303中，产生单个90度相移。该90度相移包括传感器信号的 90度相移。该90度相移可通过任何类型的相移机构或操作进行。在一个 具体实施例中，对数字传感器信号进行操作，使用希尔伯特变换来进行该 90度相移。
在步骤304中，使用单个90度相移计算相差。还可使用附加信息以 计算该相差。在一个具体实施例中，从第一传感器信号、第二传感器信号 和单个90度相移来确定该相差。该相差包括响应信号中的相差，§卩，在 拾取传感器中，由于振动的计量组件10中的科里奥利效应，所以这种情 况是可理解的。
在计算中不需要任何频率值的情况下，确定最终的相差。与使用频率 计算的相差相比，该最终的相差可更快地获得。与使用频率计算的相差相 比，该最终的相差具有更高的精确度。
在步骤305中，计算频率。根据本发明的频率有利地从90度相移计 算。 一个具体实施例中的频率利用90度相移和对应的传感器信号，该90 度相移从该对应的传感器信号获得。该频率是第一传感器信号和第二传感 器信号之一的振动响应频率(这两个传感器信号的频率在操作中基本相 等)。该频率包括该流管或这些流管响应于驱动器180产生的振动的振动频率。
不需要任何独立的频率参考信号，而获得这样导出的频率。在比现有
技术更快的操作中，从单个90度相移获得该频率。最终的频率与现有技 术计算的频率相比具有更高的精确度。
在步骤306中，计算流体材料的质量流量。从步骤304和305中计算 的最终相差和最终频率中计算该质量流量。此外，质量流量计算可从该相 差和频率计算时差(At)，最终使用该时差(At)计算质量流量。
在步骤307中，可以选择性确定密度。该密度可以确定为流动特性之 一，并且可以通过例如频率来确定。
在步骤308中，可以选择性确定体积流量。该体积流量可以确定为流
动特性之一，并且可以通过例如质量流量和密度来确定。
图4示出了根据本发明具体实施例的电子计量器20。与图2中相同
的部件用相同附图标记表示。
该具体实施例中的电子计量器20包括第一传感器信号210和第二传 感器信号211。处理系统203处理第一和第二 (数字)传感器信号210和 211，以从这些信号确定一个或多个流动特性。如前所述，该一个或多个 流动特性可以包括流体材料的相差、频率、时差(At)、质量流量、密度 和/或体积流量。
在示出的该具体实施例中，处理系统203仅仅从该两个传感器信号 210和211确定流动特性，而不需要任何外部频率测量结果也不需要外部 频率参考信号。处理系统203可从这两个传感器信号210和211至少确定 相差和频率。
如前所述，存储系统204存储相移程序212、相差程序215、频率程 序216、时差(At)程序217和流动特性程序218。存储系统204存储第 一传感器信号210和第二传感器信号211。存储系统204还存储第一 90 度相移213和第二 90度相移以确定流动特性，该第一 90度相移213和第 二90度相移从这些传感器信号产生。如前所述，存储系统204存储相差 220、频率221、时差(At) 222、质量流量223、密度224和体积流量 225。
相移程序212对输入信号，包括对第一传感器信号210和对第二传感 器信号211执行90度相移。 一个具体实施例中的相移程序212进行希尔 伯特变换(下面讨论)。
相差程序215使用第一 90度相移213和第二 90度相移214确定相
差。还可使用额外的信息来计算该相差。在一个具体实施例中，从第一传 感器信号210、第二传感器信号211、第一90度相移212和第二90度相 移213计算该相差。如前所述，确定的相差可存储在存储系统204的相差 220中。当利用第一和第二90度相移确定相差时，可比现有技术更快地 计算和获得该相差。这可在具有高流速的流量计应用中或者在多相流动发 生的应用中提供临界差。此外，该相差可独立于传感器信号210和211 的频率来确定。而且，因为该相差独立于频率确定，因此相差中的误微分 量就不受频率确定的误微分量影响，即，在相差测量中没有混合误差。因 此，在现有技术的相差基础上，相差误差减小。
频率程序216从第一 90度相移213和第二 90度相移214确定频率(例 如通过第一传感器信号210或第二传感器信号211表示的)。如前所述， 该确定的频率可存储在存储系统204的频率221中。当从第一和第二 90 度相移确定频率时，与现有技术相比，可更快地计算和获得该频率。这可 在具有高流速的流量计应用中或者在多相流动发生的应用中提供临界 差。
时差(A t)程序217确定第一传感器信号210和第二传感器信号211 之间的时差(At)。如前所述，该时差(At)可存储在存储系统204的 时差(At) 222中。该时差(At)基本上包括确定的相位除以确定的频 率，并因此用于确定质量流量。
如前所述，流动特性程序218可确定质量流量、密度和/或体积流量
中的一个或多个。
图5是根据本发明具体实施例的科里奥利流量计中处理第一和第二 传感器信号的方法流程图500。在步骤501中，接收第一传感器信号。在 一个具体实施例中，第一传感器信号包括上游或下游拾取传感器信号。
在步骤502中，接收第二传感器信号。在一个具体实施例中，第二传 感器信号包括下游或上游拾取传感器信号(即与第一传感器信号相反)。
在步骤503中，可调节这些传感器信号。在一个具体实施例中，该调 节可包括滤波以去除噪声和不期望的信号。在一个具体实施例中，如前所 述，该滤波可包括带通滤波。此外，也可执行其它调节操作，例如放大、
缓存等等。如果传感器信号包括模拟信号，则该步骤还可包括任何类型的 采样、数字化和抽取，执行这些采样、数字化和抽取以产生数字传感器信号。
在步骤504中，产生第一 90度相移。该第一 90度相移包括第一传感 器信号的90度相移。该90度相移可通过任何类型的机构或操作来执行。 在一个具体实施例中，对数字传感器信号进行操作，使用希尔伯特变换来 进行该90度相移。
在步骤505中，产生第二90度相移。该第二90度相移包括第二传感 器信号的90度相移。与第一90度相移一样，该90度相移可通过任何类
型的机构或操作进行。
在步骤506中，使用第一 90度相移和第二 90度相移计算第一传感器 信号和第二传感器信号之间的相差。还可使用附加信息以计算该相差。在 一个具体实施例中，从第一传感器信号、第二传感器信号、第一90度相 移和第二90度相移来确定该相差。该相差包括响应信号中的相差，艮P， 在两个拾取传感器中，由于振动的计量组件10中的科里奧利效应，所以 这种情况是可以理解的。
在计算中不需要任何频率值的情况下，确定最终的相差。与使用频率 计算的相差相比，可更快地获得该最终的相差。与使用频率计算的相差相 比，该最终的相差具有更高的精确度。
在步骤507中，计算频率。根据本发明的频率有利地从第一90度相 移和第二 90度相移计算。 一个具体实施例中的频率利用这些90度相移和 对应的传感器信号获得，从该对应的传感器信号获得这些90度相移。该 频率是第一传感器信号和第二传感器信号之一的振动响应频率(这两个传 感器信号的频率在操作中基本相等)。该频率包括该流管或这些流管响应 于驱动器180产生的振动的振动频率。
不需要任何独立的频率参考信号，而获得这样得到的频率。在比现有 技术更快的操作中，从这些90度相移获得该频率。最终的频率与现有技 术计算的频率相比具有更高的精确度。
在步骤508中，计算流体材料的质量流量。该质量流量从步骤506 和507中计算的最终相差和最终频率中计算。此外，质量流量计算可从该相差和频率中计算时差(At),最终使用该时差(At)计算质量流量。 在步骤509中，如前所述，可选地确定密度。
在步骤510中，如前所述，可选地确定体积流量。
图6是根据本发明具体实施例的处理系统203 —部分的结构图。在该 图中，这些方块表示处理电路或处理动作/程序。该结构图600包括1阶 滤波器块601、 2阶滤波器块602、希尔伯特变换块603和分析块604。该 LPO和RPO输入包括左拾取信号输入和右拾取信号输入。该LPO或RPO 可包括第一传感器信号。
在一个具体实施例中，1阶滤波器块601和2阶滤波器块602包括在
处理系统203中实施的数字有限脉冲响应(RIR)多相抽取滤波器。这些
滤波器提供滤波和抽取一个或两个传感器信号的最优方法，在相同的时序
下并且在相同的抽取速率下进行该滤波和抽取。可替换地，l阶滤波器块
601和2阶滤波器块602可包括有限脉冲响应(IIR)滤波器或其它合适
的数字滤波器或滤波器处理。然而，可以理解在说明书和权利要求书的范
围内，可预期其它滤波处理和/或滤波的实施例。
图7示出了根据本发明具体实施例的希尔伯特变换块603的细节。在
示出的该具体实施例中，希尔伯特变换块603包括LPO支路700和RPO支
路710。该LPO支路700包括与LPO滤波器块702并联的LPO延迟块701。
同样，RP0支路包括与RP0滤波器块712并联的RP0延迟块711。 LPO延
迟块701和RP0延迟块711引导采样延迟。因此，LPO延迟块701和RPO
延迟块711选择LPO和RPO数字信号样本，该LPO和RPO数字信号样本在
时序上晚于LPO滤波器块702和RPO滤波器块712滤波的LPO和RPO数字
信号样本。LPO滤波器块702和RPO滤波器块712对输入的数字信号样本
进行90度相移。
希尔伯特变换块603是用以提供相位测量的第一步骤。希尔伯特变换 块603接收被滤波、抽取的LPO和RPO信号，并执行希尔伯特变换。希 尔伯特变换产生LPO和RP0信号的90度相移版本，即，它产生原始的、 同相(I)信号分量的90相移(Q)分量。因此，与原始的同相(I)信号 分量LPO I和RPO I —起，希尔伯特变换块603的输出提供新的90度相 移(Q)分量LPO Q和RPO Q。
希尔伯特变换块603的输入可表示为 LP6> = 4。 cos(— (2) 肿0 =《。咖(纽+ 0) (3) 利用希耳伯特变换，输出变成
LPOMbert=Alp。 sin( w t) (4)
RPOMbert二Arp。 sin(cot + (|) ) (5)
结合原始项与希尔伯特变换的输出产生
= 4,。[c。s—) + ,sinO,)] = ^。e/(ffl° (6)
WO =《。[cos(^ + + sin( W +州=《。—(°^) (7)
图8和9是根据本发明具体实施例的分析块604两个独立分支的结构 图。分析块604是最后阶的频率、微分相位和德尔塔T (At)测量结果。 图8是包括第一分支的相位部分604，该第一分支从同相(I)和90度相 移(Q)分量确定相差。图9是频率部分604b，频率部分604b从单个传 感器信号的同相(I)和90度相移(Q)分量确定频率。该单个传感器信 号可包括LPO信号，如示出的，或者可选择地包括RPO信号。
在图8的具体实施例中，分析块604的相位部分604a包括连接块801a 和801b、共轭块802、复数乘法块803、滤波器块804和相位角度块805。
连接块801a和801b接收传感器信号的同相(I)和90度相移(Q) 分量并传递它们。共轭块802对传感器信号(这里为LPO信号)进行复共 轭，并形成负的虚信号。复数乘法块803将RP0信号与LP0信号相乘，执 行下面的方程(8)。滤波器块804进行数字滤波，例如前述的FIR滤波。 滤波器块804可包括多相抽取滤波，利用该多相抽取滤波以从传感器信号 的同相(I)和90度相移(Q)分量中去除谐波分量，并且抽取该信号。 可选择滤波系数以提供输入信号的抽取，例如抽取因子为10。相位角度 块805从LPO信号和RPO信号的同相(I)和90度相移(Q)分量中确定 相位角度。该相位角度块805执行下面示出的方程(11)。
图8中示出的相位部分604a执行下面的等式
<formula>formula see original document page 24</formula>
然后
<formula>formula see original document page 25</formula>(10)
最终的微分相位角为 sin⑨ cos⑨
图9是根据本发明分析块604的频率部分604b的结构图。频率部分 604b可工作在左或右拾取信号(LP0或RP0)上。示出的该具体实施例中 的频率部分604b包括连接块901、复共轭块902、采样块903、复数乘法 块904、滤波器块905、相位角度块906、常数块907和除法块908。
如前所述，连接块901接收传感器信号的同相(I)和90度相移(Q) 分量并传递它们。共轭块902对传感器信号(这里为LPO信号)进行复共 轭，并形成负的虚信号。延迟块903引入采样延迟到频率部分604b，因 此，选择在时序上较老的数字信号采样。该时序上较老的数字信号采样乘 以复数乘法块904中的当前数字信号。复数乘法块904将LPO信号与LPO 共轭信号相乘，执行下面的方程(12)。滤波器块905进行数字滤波，例 如前述的FIR滤波。滤波器块905可包括多相抽取滤波，利用该多相抽取 滤波以从传感器信号的同相(I)和90度相移(Q)分量中去除谐波分量， 并且抽取该信号。可选择滤波系数以提供输入信号的抽取，例如抽取因子 为IO。相位角度块906从LP0信号的同相(I)和90度相移(Q)分量中 确定相位角度。该相位角度块906执行下面的方程(13)的一部分。该常 数块907提供因子，如等式(14)所示，该因子包括采样速率FJ余以两 倍pi。除法块908执行等式(14)的除法操作。
，部分604b执行下面的方程
<formula>formula see original document page 25</formula>(12)
两个连续样本之间的角度因此为
i[sm(K)] (13)
其是左拾取的角频率，转换成赫兹<formula>formula see original document page 25</formula>(14)
其中，Fs是希尔伯特转换块603的速率。在前述例子中，"Fs"约为
2kHz。
图10是在一般条件下的流量计拾取传感器信号的功率频谱密度图。 该流量计的基础频率为图的最高峰值，且位于约135Hz。该图还示出了该 频谱中的几个其它的大峰值(第一非基础模式是频率在约1. 5倍基础模式 频率下的扭曲模式)。这些峰值包括流量计的谐波频率，并且还包括其它 的、不期望的传感器模式(即扭曲模式、第二弯曲模式等等)。
图11示出了根据单相移具体实施例的可替换希尔伯特变换块603'。 该具体实施例中的希尔伯特变换块603'包括LP0支路1100和RP0支路 1110。该LP0支路1100包括与滤波器块702并联的延迟块701。该具体 实施例中的RPO支路1110仅包括延迟块701。如前所述，该延迟块701 引入釆样延迟。如前所述，滤波器块702对输入的数字信号样本进行90 度相移。可以理解，可替换的希尔伯特变换块603'可正好对RP0信号进 行相移。
该处理实施例利用仅仅一个传感器信号的希尔伯特变换/相移来获得 频率和相差(见图2—3)。这大大地减少了进行相位测量所需要的计算数 量，并且大大地减少了获得质量流量所需的计算数量。
该实施例中，希尔伯特变换块603'的输出将提供左或右传感器信 号、但不是两者的90度相移(Q)分量。在下面的示例中，对该LPO信号 进行了相移。
<formula>formula see original document page 26</formula> (27)
利用希耳伯特变换，该输出变成 <formula>formula see original document page 26</formula> (29)
结合LPO原始项与希尔伯特变换的输出(即通过90度相移)，产生
<formula>formula see original document page 26</formula> (30)
而RPO保持相同
一一 +e-)
<formula>formula see original document page 26</formula> (31)
图12示出了单个相移实施例的分析块604a'。该具体实施例中的分析块604a'包括一个连接块801、复数乘法块803、低通滤波器块1201和 相位角度块805。该具体实施例中的分析块604a'执行下面的等式
<formula>formula see original document page 27</formula>(32)
低通滤波器块1201包括低通滤波器，该低通滤波器去除复数乘法块 803产生的高频分量。低通滤波器块1201可进行任何类型的低通滤波操 作。乘法操作的结果产生两个项。由于(一cot)和(cot)项可互相抵消， 所以项(一cot+"t+dO结合并简化为仅仅相位小项(DC结果)。在两 倍频率下，该("t+"t+4))项简化为(2cot+4))。因为该结果是两 项和的结果，所以高频项(2"t+小)可去除。这里关注的唯一信号是 DC项。高频项(2"t+(tO可利用低通滤波器从结果中滤波掉。低通滤 波器的截止可定位在0到2co之间的任何位置。
滤波后，结果为
<formula>formula see original document page 27</formula>因此，微分相位角是
—tan-(34)
COS(0)
通过对一个拾取信号而不是对两个拾取信号进行希尔伯特变换，有利 地减小了在科里奥利质量流量计中进行相位和频率估计所需的计算负
荷。因此，该相位和频率可利用两个传感器信号、但只利用一个90度相
移来确定。
图13示出了与现有技术相比的本发明传感器处理，其中比较每一个 时差(At)值。该图表示出了包括气流(例如气泡)的流体材料。在该 条件下，因为相位和频率计算速率，所以流体噪声基本上在新的算法中减 小。从图可见，本发明获得的结果没有显示现有技术(At)测量中反映 出的那些大峰值和谷值。
图14示出根据本发明的另一个实施例的电子计量器20。该实施例的 电子计量器20可以确定流过相关流量计的流体材料的气体空隙度。该电 子计量器20可以接收流体材料的频率响应，将该频率响应分解成至少气 体频率分量和流体频率分量，并且从该频率响应以及一个或多个气体频率
分量和流体频率分量中来确定气体空隙度。该空隙度可以包括基本瞬时的 空隙度。该电子计量器20因此可以确定流体的流动特性，而不管含有的 气体的存在与否或量的多少。
该流体可以包括流体的任何形式，并且该气体可以包括气体的任何形 式。该气体可以包括例如空气。下面的说明集中在流体中的空气，但是应 当理解的是本发明可以应用到任何气体。
该实施例中的处理系统203可以包括空隙度程序1401、陷波滤波器 程序1402以及一个或多个滤波器或滤波器程序，例如低通滤波器程序 1403和高通滤波器程序1404。可替换地，该一个或多个滤波器或滤波器 程序可以包括陷波滤波器结构或可以滤去窄带频率的其它滤波器结构。该 处理系统203还包括频率响应1410和空隙度1418，它们能够分别存储频 率响应测量值和空隙度确定值。该处理系统203还包括流体频率分量1416 和气体频率分量1412，它们存储用于空隙度确定的工作频率值。该处理 系统203还包括整个密度1420、气体分量密度1421以及流体分量密度 1422，它们存储用于空隙度确定的工作密度值。
一个实施例中的频率响应1410包括前述从计量组件10接收到的频率 221 (参见图2和4)。该频率响应1410利用前述处理块来确定(参见图6 一7和9)。有利地，当利用前述高速频率确定时，本发明可以快速地、准 确地以及可靠地确定气体1418的空隙度。
在一个实施例中，频率响应1410从第一和第二传感器信号中导出， 其中这两个传感器信号被处理系统203处理，以便获得频率响应1410。 该处理可以包括相移一个传感器信号，如前所述，然后基于相移信号和这 两个传感器信号来产生频率值。
频率响应1410包括混合频率(fmix),其中频率响应1410可以包括气 体频率分量(fgas) 1412和流体频率分量(fflmd) 1416。在任何时候，频率 响应1410可以包括任意量的气体频率分量(fgas) 1412 (也就是，夹带气 体)。
图15是用于空气、用于流体以及用于组合的空气/流体混合物(也就 是，包括夹带空气的流体)的流量计频率响应的曲线图。本发明的基本理 念是，气体的密度可以区分于流过流量计的流体材料中流体的密度。由于
可以从通过流量计获得的测量频率中导出密度，所以空气相关的频率也可 以区分于流体的频率。这对于其它气体或气体混合物也是如此。 用于计算频率的等式为
<formula>formula see original document page 29</formula>(35)
其中C0是科里奥利流量计的角频率。该(Dd项表示从之前的或者较老
的样本周期获得的角频率样本。将该角频率co转换成单位赫兹(Hz)的
频率f得到
<formula>formula see original document page 29</formula>(36)
该等式假设仅存在一个频率。如果存在两个频率，如夹带空气的情况
下(空气的频率和流体材料流体的频率)，则新等式变成
<formula>formula see original document page 29</formula>(37)
其中fmix是整个流体材料的频率响应，包括气体频率分量(fgas)和流 体频率分量(fflmd)。
低通滤波器程序1403实现低通滤波器。低通滤波器通过基本低于低
通截止频率的低频。因此低通滤波器用于去除高频。
高通滤波器程序1404实现高通滤波器。高通滤波器通过基本高于高 通截止频率的高频。因此，高通滤波器用于去除低频。
陷波滤波器程序1402实现陷波滤波器。陷波滤波器滤去窄范围的频 率，该窄范围的频率中心位于陷波滤波器的频率响应中的"槽口"上。仅仅 槽口中的频率由陷波滤波器滤去。因此，陷波滤波器对于从频率响应1410 中去除已知的、不期望的频率非常有用。
在操作中，电子计量器20接收频率响应1410。电子计量器20将频 率响应1410分解成气体频率分量1412和流体频率分量1416。电子计量 器20采用频率响应1410和一个或多个气体频率分量1412和流体频率分 量1416，从而确定气体空隙度。电子计量器20从频率响应1410中确定 整个密度(pmix) 1420。类似的，从气体频率分量1412中确定体分量密度 (pai》1421，以及从流体频率分量1416中确定流体分量密度(pfluid) 1422 (参见下面的等式38和39)。这些密度值用于确定空隙度。因此，对于流
体材料中的基本任何气体和基本任何流体确定空隙度。 密度(P)包括
<formula>formula see original document page 30</formula>(38)
其中m是质量流量测量值，f是频率测量值，以及b是某些形式的校 准因数。因此，气体空隙度1418计算为等于流体分量密度(pfluid) 1422 减去整个密度(pmix) 1420的差除以流体分量密度(pfluid) 1422减去气体 分量密度(pgas) 1421的差得到的比率。计算具有以下形式
<formula>formula see original document page 30</formula>3 9)
最终的气体空隙度1418反映了流体材料中气体和流体之间的比率。
应当理解的是，如果气体或流体是已知的，则仅需要气体频率分量 1412和流体频率分量1416之一。例如，如果气体包括空气，则可以假设 标准空气频率响应(和密度)。因此，可以滤出已知的气体或流体频率， 并且仅需要一个滤波步骤。
电子计量器20还可以确定其它流动特性，例如整个质量流量、仅仅 流体分量的质量流量、流体密度等等。电子计量器20可以和计量组件10 进行通信，其中计量组件10包括产生频率响应的任何形式的流量计。在 一个实施例中，计量组件10包括科里奥利流量计。在另一个实施例中， 计量组件10包括振动密度计。
应当理解的是，该附图中的电子计量器20可以包括前述的各种其它 部件和/或程序。和其它实施例相同的部件和/或程序采用相同的附图标记 来表示。电子计量器20如前所述可以确定流体材料的频率响应和密度。
从频率响应中分解流体和气体频率分量的一种方法包括执行两个滤 波操作。 一个滤波操作包括利用第一滤波器对频率响应进行滤波，该第一 滤波器基本滤去气体频率分量并且基本通过流体频率分量。第二滤波操作 包括利用第二滤波器对频率响应进行滤波，该第二滤波器基本滤去流体频 率分量并且基本通过气体频率分量。因此，第一滤波器输出流体频率分 量，而第二滤波器输出气体频率分量。
图16是根据本发明的实施例用于确定流过流量计的流体材料中气体 空隙度的方法的流程图1600。在步骤1601中，接收频率响应。例如可以
在电子计量器20中接收该频率响应。该频率响应包括响应于振动计量组 件IO (包括流体材料)的频率响应。该流体材料可以包括气体。
在一个实施例中，该频率响应可以包括第一传感器信号和第二传感器
信号。可以从例如计量组件10的拾取传感器170L和170R中接收第一传 感器信号和第二传感器信号。从任一传感器信号中可以产生90度相移。 该90度相移和第一和第二传感器信号可以用于计算频率响应，其中该频 率响应根据流体材料的质量流量以及根据所夹带气体的存在与否和量大 小而随时间变换。
在步骤1602中，频率响应分解成气体频率分量1412和流体频率分量 1416。这是可以实现的，因为频率响应1410包括和流体材料中气体流量 相关的气体频率分量以及和流体流量相关的流体频率分量。该流体包括任 何形式的流体。该分解步骤可以通过第一滤波器和第二滤波器来执行，如 前所述。此外，该分解步骤可以通过低通滤波器和高通滤波器来执行，如 前所述。
在步骤1603中，如前所述，利用频率响应1410、气体频率分量1412 和流体频率分量1416来确定气体空隙度1418。该确定步骤可以包括从频 率响应1410、气体频率分量1412和流体频率分量1416中来确定密度值， 如前所述。最终的气体空隙度1418可以表示为比率、百分比或其它测量值。
图17是根据本发明的实施例用于确定流过流量计的流体材料中气体 空隙度的方法的流程图1700。在步骤1701中，接收频率响应，如前所述。
在步骤1702中，通过第一滤波器滤波该频率响应。第一滤波器基本 滤去气体频率分量并基本通过流体频率分量(见图18)。在一个具体实施 例中，第一滤波器包括低通滤波器，其中该低通滤波器的低通截止频率基 本上在流体频率分量以上。因此，低通滤波器基本通过流体频率分量并基 本滤去气体频率分量。
在步骤1703中，通过第二滤波器滤波该频率响应。第二滤波器基本 滤去流体频率分量并基本通过气体频率分量。在一个具体实施例中，第二 滤波器包括高通滤波器，其中该高通滤波器的高通截止频率基本上在气体 频率分量以下(但在流体频率分量以上)。因此，高通滤波器基本通过气
体频率分量并基本滤去流体频率分量。
在步骤1704中，如前所述，利用频率响应1410、气体频率分量1412 和流体频率分量1416来确定空气空隙度1418。
图18是频率图表，其示出了根据本发明实施例的低通和高通滤波器
响应，该低通和高通滤波器响应用于分解成流体频率分量和气体频率分
量。该图的下部线条表示包括流体频率分量波瓣和气体频率分量波瓣的流 量计频率响应。上部线条包括沿着截止频率的低通滤波器响应和高通滤波
器响应。这里，对于低通滤波器和高通滤波器的截止频率的中心基本上都 位于上述两个波瓣之间。低通滤波器和高通滤波器可具有共同的截止频 率，或者可以具有不同的截止频率，这取决于流体和气体频率分量而定。 可以看到，低通滤波器将输出流体频率分量，高通滤波器将输出气体频率 分量。因此，这两个滤波器可将频率响应1410分解成气体频率分量1412 和流体频率分量1416。
分解成流体和气体频率分量的另一方法包括滤波单个的、已知的频率 分量并利用滤波操作通过的频率分量来确定流体和气体分量密度。例如， 在流体材料中的气体为空气的情况下，那么可配置该滤波操作以滤波相对 窄的频带，该相对窄的频带以一般的空气频率响应为中心。随后，从频率 响应导出的整个密度以及从剩下的流体频率分量导出的流体密度分量可 用于确定空气密度项。例如，在气体己知为大气压下空气的情况下，滤波 器(例如陷波滤波器)可用于基本上滤去该频率响应的空气频率分量。因 此，整个密度(P mix) 1420可从频率响应1410中计算，流体分量密度(P fluid) 1422可从流体频率分量1416中计算。因此，空气分量密度(Pgas) 1421为整个密度(pmix) 1420减去流体分量密度(pfMd) 1422。
可替换地，可以理解的是，流体频率分量可去除/滤波掉，并且空隙 度可利用气体频率分量确定。如前所述，在该流体拥有已知的特征频率响 应和密度的情况下，可进行单个频率去除。因此，单个频率去除的方法可 去除流体频率分量或气体频率分量。
在一个具体实施例中，可通过一个或多个滤波器去除单个频率分量，
而其它频率分量可通过该滤波操作。 一个具体实施例中的一个或多个滤波 器包括陷波滤波器。除了窄带范围内的频率(即频率响应中的槽口)之外，
陷波滤波器通过所有频率。可替换地， 一个或多个滤波器可包括任何令人 满意的滤波器或滤波器组合。
图19是根据本发明的实施例用于确定流过流量计的流体材料中气体
空隙度的方法的流程图1900。在步骤1901中，接收频率响应1410，如前所述。
步骤1902中，通过陷波滤波器处理该频率响应。陷波滤波器通过槽 口以下或以上的频率，例如在该具体实施例中通过该气体频率响应以上或 以下的频率。因此，陷波滤波器基本上滤去气体频率分量1412。陷波滤 波器基本通过流体频率分量1416。
图20是陷波滤波器频率响应图。在示出的示例中，该槽口以气体频 率为中心。该陷波滤波器基本通过该槽口以上和以下的所有频率，并且该 陷波滤波器基本上仅仅滤去气体频率。
在步骤1903中，利用频率响应、气体频率分量1412和流体频率分量 1416来确定气体空隙度1418，如前所述。
根据本发明的电子计量器和方法可以根据任何实施例来实现，从而获 得所需的若干好处。本发明能够确定两相流体中的气体空隙度。本发明能 够确定多相流体中的气体空隙度。本发明能够确定空气的空隙度。本发明 能够提供更高准确度和可靠度的气体空隙度确定。本发明能够提供比现有 技术更快同时消耗更少处理时间的气体空隙度确定。
权利要求
1.用于确定流过流量计(5)的流体材料中气体空隙度的电子计量器(20)，该电子计量器(20)包括用于接收流体材料的频率响应的接口(201)；以及与该接口(201)通信的处理系统(203)，并且该处理系统配置成从该接口(201)接收频率响应，将该频率响应分解成至少气体频率分量和流体频率分量，并且通过该频率响应以及一个或多个气体频率分量和流体频率分量来确定气体空隙度。
2. 根据权利要求1所述的电子计量器(20)，其中该气体空隙度包括 空气空隙度。
3. 根据权利要求l所述的电子计量器(20)，其中该接口 (201)包 括数字化该频率响应的数字转换器(202)。
4. 根据权利要求1所述的电子计量器(20)，其中该处理系统(203) 还配置成用于从流量计(5)接收第一传感器信号和第二传感器信号，从 第一传感器信号中产生90度相移，并且利用该90度相移来计算该频率响 应。
5. 根据权利要求l所述的电子计量器(20)，其中该处理系统(203) 还配置成用于从流量计(5)接收第一传感器信号和第二传感器信号，从 第一传感器信号产生90度相移，并且利用该第一传感器信号和该90度相 移来计算该频率响应。
6. 根据权利要求1所述的电子计量器(20)，其中该处理系统(203) 还配置成通过利用一个或多个滤波器处理该频率响应，从而分解该频率响 应，该一个或多个滤波器基本滤去气体频率分量和流体频率分量之一 。
7. 根据权利要求l所述的电子计量器(20)，其中该处理系统(203)还配置成用于利用陷波滤波器处理该频率响应，从而分解该频率响应，该 陷波滤波器基本滤去气体频率分量和流体频率分量之一。
8. 根据权利要求1所述的电子计量器(20)，其中该处理系统(203) 还配置成用于利用第一滤波器对该频率响应进行滤波以及利用第二滤波 器对该频率响应进行滤波，从而分解该频率响应，其中该第一滤波器基本滤去气体频率分量以及基本通过流体频率分量，该第二滤波器基本滤去流 体频率分量以及基本通过气体频率分量，其中该第一滤波器输出该流体频 率分量，该第二滤波器输出该气体频率分量。
9. 根据权利要求l所述的电子计量器(20)，其中该处理系统(203) 还配置成用于利用低通滤波器处理频率响应从而分解频率响应，该低通 滤波器的低通截止频率基本高于流体频率分量，该低通滤波器基本通过流 体频率分量并且基本滤去气体频率分量；以及利用高通滤波器处理频率响应从而分解频率响应，该高通滤波器的高通截止频率基本低于气体频率分 量，该高通滤波器基本通过气体频率分量并且基本滤去流体频率分量，其 中该低通滤波器输出该流体频率分量，该高通滤波器输出该气体频率分
10. 根据权利要求1所述的电子计量器(20)，其中该处理系统(203) 还配置成用于从频率响应计算整个密度，从而确定气体空隙度；从流体 频率分量计算流体分量密度；从气体频率分量计算气体分量密度；以及计 算气体空隙度，该气体空隙度等于该流体分量密度减去该整个密度的差除 以该流体分量密度减去该气体分量密度的差的比率。
11. 根据权利要求l所述的电子计量器(20)，其中该流量计(5)包 括科里奧利流量计(5)。
12. —种用于确定流过流量计的流体材料中气体空隙度的方法，该方 法包括接收流体材料的频率响应；将该频率响应分解成至少气体频率分量和流体频率分量；以及 从该频率响应以及一个或多个气体频率分量和流体频率分量中来确定气体空隙度。
13. 根据权利要求12所述的方法，其中该气体空隙度包括空气空隙度。
14. 根据权利要求12所述的方法，其中该接收频率响应的步骤还包括从流量计接收第一传感器信号和第二传感器信号； 从第一传感器信号中产生90度相移；以及利用该90度相移计算该频率响应。
15. 根据权利要求12所述的方法，其中该接收频率响应的步骤还包括从流量计接收第一传感器信号和第二传感器信号； 从第一传感器信号中产生90度相移；以及利用该第一传感器信号和该90度相移来计算该频率响应。
16. 根据权利要求12所述的方法，其中该分解步骤包括利用一个或多个滤波器处理该频率响应，该一个或多个滤波器基本滤去气体频率分 量和流体频率分量之一。
17. 根据权利要求12所述的方法，其中该分解步骤包括利用陷波滤波器处理该频率响应，该陷波滤波器基本滤去气体频率分量和流体频率 分量之一。
18. 根据权利要求12所述的方法，其中该分解步骤包括利用第一滤波器对该频率响应进行滤波，该第一滤波器基本滤去气体频率分量以及基本通过流体频率分量；以及利用第二滤波器对该频率响应进行滤波，该第二滤波器基本滤去流体频率分量以及基本通过气体频率分量；其中该第一滤波器输出该流体频率分量，该第二滤波器输出该气体频 率分量。
19. 根据权利要求12所述的方法，其中该分解步骤包括利用低通滤波器处理频率响应，该低通滤波器的低通截止频率基本高 于流体频率分量，其中该低通滤波器基本通过流体频率分量并且基本滤去气体频率分量；以及利用高通滤波器处理频率响应，该高通滤波器的高通截止频率基本低于气体频率分量，其中该高通滤波器基本通过气体频率分量并且基本滤去流体频率分量；其中该低通滤波器输出该流体频率分量，该高通滤波器输出该气体频率分量。
20. 根据权利要求12所述的方法，其中该确定步骤包括从频率响应中计算整个密度；从流体频率分量中计算流体分量密度； 从气体频率分量中计算气体分量密度；以及计算气体空隙度，该气体空隙度等于该流体分量密度减去整个密度的 差除以该流体分量密度减去气体分量密度的差的比率。
21. 根据权利要求12所述的方法，其中该流量计包括科里奧利流量计。
22. —种用于确定流过流量计的流体材料中气体空隙度的方法，该方法包括接收流体材料的频率响应；利用陷波滤波器处理该频率响应，该陷波滤波器基本滤去气体频率分量和流体频率分量之一；以及从该频率响应以及一个或多个气体频率分量和流体频率分量中来确定气体空隙度。
23. 根据权利要求22所述的方法，该气体空隙度包括空气空隙度。
24. 根据权利要求22所述的方法，该接收频率响应的步骤还包括 从流量计接收第一传感器信号和第二传感器信号； 从第一传感器信号产生90度相移；以及利用该90度相移计算该频率响应。
25. 根据权利要求22所述的方法，该接收频率响应的步骤还包括 从流量计接收第一传感器信号和第二传感器信号； 从第一传感器信号产生90度相移；以及 利用该第一传感器信号和该90度相移来计算该频率响应。
26. 根据权利要求22所述的方法，该确定步骤包括从频率响应计算整个密度；如果通过陷波滤波器滤去气体频率分量，则从流体频率分量中计算流体分量密度，或者如果滤去流体频率分量则计算气体分量密度；如果通过陷波滤波器滤去气体频率分量，则估测气体分量密度，或者如果滤去流体频率分量则估测流体分量密度；计算气体空隙度，该气体空隙度等于该流体分量密度减去该整个密度的差除以该流体分量密度减去该气体分量密度的差的比率。
27. 根据权利要求22所述的方法，其中该流量计包括科里奥利流量计。
28. —种用于确定流过流量计的流体材料中气体空隙度的方法，该方法包括接收流体材料的频率响应；利用第一滤波器对频率响应进行滤波，该第一滤波器基本滤去气体频 率分量并且基本通过流体频率分量，其中该第一滤波器输出该流体频率分利用第二滤波器对频率响应进行滤波，该第二滤波器基本滤去流体频 率分量并且基本通过气体频率分量，其中该第二滤波器输出该气体频率分从该频率响应以及一个或多个气体频率分量和流体频率分量中来确 定气体空隙度。
29. 根据权利要求28所述的方法，该气体空隙度包括空气空隙度。
30. 根据权利要求28所述的方法，该接收频率响应的步骤还包括 从流量计接收第一传感器信号和第二传感器信号； 从第一传感器信号产生90度相移；以及利用该90度相移计算该频率响应。
31. 根据权利要求28所述的方法，该接收频率响应的步骤还包括 从流量计接收第一传感器信号和第二传感器信号； 从第一传感器信号产生90度相移；以及 利用该第一传感器信号和该90度相移来计算该频率响应。
32. 根据权利要求28所述的方法，该确定步骤包括 从频率响应计算整个密度；从流体频率分量计算流体分量密度； 从气体频率分量计算气体分量密度；以及计算气体空隙度，该气体空隙度等于该流体分量密度减去该整个密度的差除以该流体分量密度减去该气体分量密度的差的比率。
33. 根据权利要求28所述的方法，该流量计包括科里奥利流量计。
34. 根据权利要求28所述的方法，该第一滤波器包括低通滤波器， 该第二滤波器包括高通滤波器。
全文摘要
根据本发明的实施例提供用于确定流过流量计(5)的流体材料中气体空隙度的电子计量器(20)。该电子计量器(20)包括用于接收流体材料的频率响应的接口(201)；以及与该接口(201)通信的处理系统(203)。该处理系统(203)配置成从该接口(201)接收频率响应，将该频率响应分解成至少气体频率分量和流体频率分量，并且从该频率响应以及一个或多个气体频率分量和流体频率分量中来确定气体空隙度。
文档编号G01F1/84GK101194148SQ200680017445
公开日2008年6月4日 申请日期2006年5月19日 优先权日2005年5月20日
发明者C·B·麦卡纳利, M·J·贝尔 申请人:微动公司