用于确定流体的流属性的设备 【技术领域】
本发明涉及用于确定流体的流属性的设备、 方法和计算机程序。背景技术 S. Sudo 等人的论文 “Detection of small particles in fluid flow using a self-mixing laser” ，Optics Express, Vol. 115, Issue 13, pp. 8135-8145 公开了一 种通过使用具有极高光学灵敏度的激光二极管泵浦的薄片固体激光器的自混合激光多普 勒测量来检测流体流中的小颗粒的实时方法。 由来自在稀样本流中移动通过小直径玻璃管 的小颗粒的重新注入散射光调制的激光输出的不对称功率谱被观察并且显示出反应流体 流的速度分布， 该速度分布遵从泊肃叶 （Poiseuille） 定律。流体流的速度分布与该不对称 功率谱的依存性被确定并且用于基于测量的不对称功率谱确定流体流的速度分布。
这种方法具有的缺点为只有在流体流具有小的光学厚度时， 可以确定流体流。对 于更大的光学厚度， 无法确定流体流。
发明内容 本发明的目的是提供用于确定流体的流属性的设备、 方法和计算机程序， 其允许 确定更大光学厚度的流体的流属性。
在本发明的一方面中， 提供了一种用于确定流体的流属性的设备， 其中该设备包 括： - 距离和速度确定单元， 其用于确定流体的元到该距离和速度确定单元的距离并且用 于同时确定元的速度， 该距离和速度确定单元包括具有激光器腔体的激光器， 其中该距离 和速度确定单元调适为， 通过将在激光器腔体内生成的激光辐射引导到流体用于被流体反 射并且通过将反射辐射与激光器腔体内的辐射混合而生成自混合干涉信号， 并且基于所生 成的自混合干涉信号确定距离和速度， - 流确定单元， 其用于基于所确定的距离和速度的至少其一确定流体的流属性。
由于流体的元到该距离和速度确定单元的距离以及元的速度同时被确定， 所述元 所在之处与该距离和速度确定单元的距离以及这些元在这些距离处具有的速度已知。 这允 许基于所确定的距离和速度的至少其一确定流体的流属性。 例如， 如果流体的光学厚度小， 使得不一定需要元的所确定的距离用于确定期望流属性， 例如因为激光辐射可以深深地穿 透到流体内使得足以用于确定期望流属性的速度被确定， 该流确定单元可以仅仅基于所确 定的速度确定流属性。 然而， 如果光学厚度使得所确定的速度不足以用于确定流属性， 例如 因为激光辐射无法穿透流体内足够深， 该流确定单元可以使用所确定的速度和所确定的距 离用于确定流属性。 因而， 即使光学厚度大， 通过使用同时确定的速度和距离而可以确定流 属性。
该距离和速度确定单元调适为， 通过将在激光器腔体内生成的激光辐射引导到流 体用于被流体反射以及通过将反射辐射与激光器腔体内的辐射混合而生成自混合干涉信
号。激光器腔体内的这种混合引起激光功率的波动， 该波动可以作为自混合干涉信号被检 测。这种自混合干涉信号依赖于流体的元的速度和距离， 其中该距离和速度确定单元调适 为根据所生成的自混合干涉信号确定流体的这些元的距离和速度。
同时确定流体的元的距离和速度是指相同的被元反射的反射辐射被用于确定这 个元的距离和速度， 即这个元的所确定的距离和速度同时描述这个元的距离和速度。
流体的元为例如流体本身的元和 / 或添加到流体的元。
流确定单元调适为基于所确定的距离和速度的至少其一确定流体的流属性， 这是 指距离、 速度、 或者距离和速度二者被用于确定流体的流属性。
优选的是， 该流确定单元调适为确定最大流速度和体积流量 （volume flow） 的至 少其一作为流体的属性。 体积流量优先地定义为在预定义时间间隔内流动经过流体的截面 的流体体积。因而， 体积流量也可以看作是体积流率。
另外优选的是， 该距离和速度确定单元调适为 ： - 根据该自混合干涉信号确定多普勒频率， - 确定所确定的多普勒频率的最大多普勒频率， - 根据所确定的最大多普勒频率确定流体的元的最大流速度， 其中该流确定单元调适为确定该最大流速度作为流属性。 这允许容易地且高精确度地 确定最大流速度。
另外优选的是， 该流确定单元调适为 ： - 提供定义最大流速度和体积流量之间的关系的体积流量函数， - 通过使用该体积流量函数和该最大流速度， 确定该体积流量作为流属性。
另外优选的是， 该设备还包括流宽度确定单元， 其用于根据元的所确定的距离确 定流的宽度。
如果激光辐射横穿该流， 则在流以外不存在反射或散射元， 或者如果流体置于管 内， 则在管边缘以外不存在反射或散射元， 并且因此将不会分别从流之外或管之外返回距 离信息。此外， 相对于发射的激光辐射的传播方向在流的前方， 或者如果流体置于管内， 在 管的前方不存在流体的散射或反射元， 并且因此从该位置将不返回距离信息。 因而， 通过确 定到该距离和速度确定单元的最近距离和最大距离， 可以确定流的宽度。这个所确定的宽 度可以例如被用于， 例如通过将所确定的宽度与流的已知宽度比较， 确定激光辐射是否完 全穿透该流。
另外优选的是， 该流确定单元调适为 ： - 提供流模型函数， 该流模型函数定义依赖于流体的元到该距离和速度确定单元的距 离的元的速度， - 将该流模型函数拟合到元的所确定的距离和速度， - 根据所拟合的流模型函数确定流属性。
该流模型优先地为层流流模型， 该模型假设最大流速度是在流的中间并且零流量 值 （flow value） 位于流的边缘。
即使仅仅流体的少许元的距离和速度被确定， 仍可以执行流模型函数到流体的元 的所确定的距离和速度的拟合。因而， 即使流体在光学上是厚的， 仍可以执行这种拟合。这 因此进一步提高了确定具有大光学厚度的流体的流属性的能力。另外优选的是， 该设备还包括 ： - 流宽度确定单元， 其用于根据元的所确定的距离确定流的宽度， 其中该距离和速度确定单元和该流确定单元调适为使得 a） 如果所确定的宽度等于或大于预定义最大速度宽度， 该距离和速度确定单元确定该 自混合干涉信号的最大频率并且根据所确定的最大频率确定流体的元的最大流速度， 并且 该流确定单元确定该最大流速度作为流属性， b） 如果所确定的宽度小于预定义最大速度宽度， 该流确定单元提供定义依赖于流体的 元到该距离和速度确定单元的距离的元的速度的流模型函数， 将该流模型函数拟合到元的 所确定的距离和速度， 并且根据所拟合的流模型函数确定流属性。
该预定义最大速度宽度定义流的宽度， 从而允许根据所生成的自混合干涉信号确 定最大流速度， 其中该流的宽度至少必须由流宽度确定单元确定。流的这个所确定的宽度 依赖于流体的光学厚度， 即在流体中辐射的穿透深度。 因而， 通过依赖于流的所确定的宽度 确定流体的流属性， 流属性的确定依赖于流体的光学厚度。 如果光学厚度小， 使得激光辐射 达到最大速度宽度， 则该距离和速度确定单元和该流确定单元调适为， 通过根据该自混合 干涉信号确定多普勒频率， 通过确定所确定的多普勒频率的最大多普勒频率， 并且通过根 据所确定的最大多普勒频率确定流体的元的最大流速度， 确定最大流速度。如果光学厚度 大得使得激光辐射无法达到最大速度宽度， 该流确定单元提供定义依赖于流体的元到该距 离和速度确定单元的距离的元的速度的流模型函数， 将该模型函数拟合到元的所确定的距 离和速度， 并且根据所拟合的流模型函数确定流属性。这允许确定依赖于流体的光学厚度 的流体的流属性。 另外优选的是， 该流确定单元调适为根据该自混合干涉信号确定流体的流是层流 还是湍流作为流属性。确定流体的流是层流还是湍流可以用于控制流体的流， 使得该流保 持或者变为层流。 例如， 如果流体被泵浦通过管， 泵压力可以被控制使得， 如果该流为湍流， 泵压力减小使得该流变为层流。这最优化了由于内部摩擦引起的流内的损耗。
优先地， 该流确定单元调适为， 在自混合干涉信号的频率谱图具有混沌行为时确 定该流是湍流， 以及其中该流确定单元调适为， 在自混合干涉信号的频率谱图不具有混沌 行为时确定该流是层流。
在本发明的另一方面中， 提供了一种用于确定流体的流属性的方法， 其中该方法 包括下述步骤 ： - 确定流体的元到距离和速度确定单元的距离以及用于同时确定元的速度， 其中通过 将在激光器腔体内生成的激光辐射引导到流体用于被流体反射而生成自混合干涉信号， 其 中反射辐射与激光器腔体内的辐射混合， 并且其中基于所生成的自混合干涉信号确定距离 和速度， - 基于所确定的距离和速度的至少其一确定流体的流属性。
在本发明的又一方面中， 提供了一种用于确定流体的流属性的计算机程序， 其中 该计算机程序包括程序代码装置， 其用于在该计算机程序在控制设备的计算机上运行时， 致使如权利要求 1 限定的该设备实施如权利要求 10 限定的方法的步骤。
应理解， 本发明的优选实施例也可以是从属权利要求与相应独立权利要求的任何 组合。
附图说明
参考下文描述的实施例， 本发明的这些和其它方面将清楚明白并得到阐释。在下 述附图中 ： 图 1 示意性和示例性示出用于确定流体的流属性的设备的实施例， 图 2 示意性和示例性示出用于确定流体的流属性的设备的距离和速度确定单元的实 施例， 图 3 示意性和示例性示出依赖于流体流内的归一化位置的归一化速度， 图 4 示意性和示例性示出依赖于流体流内的归一化速度的元数目， 图 5 示意性和示例性示出自混合干涉信号的谱图， 图 6 示意性和示例性示出用于确定流体的流属性的布置， 以及 图 7 示例性示出流程图， 该流程图说明用于确定流体的流属性的方法的实施例。具体实施方式
图 1 示意性和示例性示出用于确定流体 2 的流属性的设备 1。该设备包括距离和 速度确定单元 3， 其用于确定流体的元到距离和速度确定单元 3 的距离并且用于同时确定 元的速度。距离和速度确定单元 3 包括具有激光器腔体 6 的激光器 5（示于图 2） 。距离和 速度确定单元 3 调适为， 通过将在激光器腔体 6 内生成的激光辐射 7 引导到流体 2 用于被流 体 2 反射并且通过将反射辐射 8 与激光器腔体 6 内的辐射混合， 生成自混合干涉信号。距 离和速度确定单元 3 还调适为基于所生成的自混合干涉信号确定距离和速度。
用于确定流体 2 的流属性的设备 1 还包括流确定单元 4， 其用于基于所确定的距离 和速度的至少其一确定流体 2 的流属性。
在此实施例中， 流体 2 在管 10 内流动。
泵 11 仅仅示意性示于图 1， 该泵 11 可存在以用于控制管 10 内流体 2 的流。泵 11 可以与流确定单元 4 连接， 从而控制泵 11 使得获得预定流量值。
距离和速度确定单元 3 调适为， 通过将在激光器腔体 6 内生成的激光辐射 7 引导 到流体 2 用于被流体 2 反射并且通过将反射辐射 8 与激光器腔体 6 内的辐射混合， 生成自 混合干涉信号。激光器腔体 6 内的这种混合引起激光功率的波动， 该波动被示于图 2 的检 测器 12 检测。检测器 12 通过检测从激光器腔体 6 耦出的辐射 13 的强度而检测激光功率。 由检测器 12 检测到的自混合干涉信号依赖于流体 2 的元的速度和距离， 并且距离和速度确 定单元 3 包括分析单元 14， 其调适为根据所生成的自混合干涉信号确定流体 2 的这些元的 距离和速度。
该距离和速度确定单元首先确定速度在辐射 7 方向上的分量， 其中通过知晓流方 向和辐射 7 方向之间的角而用三角法确定在流方向上的速度， 例如通过将在辐射 7 方向上 的分量乘以此角的余弦。
流确定单元 4 调适为确定最大流速度和体积流量的至少其一作为流体的属性。体 积流量优先地定义为在预定义时间间隔内流动经过流体截面的流体 2 的体积。因而， 体积 流量也可以看作是体积流率。
流体 2 流优先地为层流流 （laminar flow） ， 层流流由抛物线速度分布表征。对于管 10 内液体的层流流， 速度在管边界为零并且在管 10 的中心具有最大值。这示例性且示 意性示于图 3。
图 3 示 出 依 赖 于 管 1 0 的 归 一 化 半 径 的归一化速度 。 归一化半径 0.0 表示管 10 的中心， 并且归一化半径 -1.0 和 1.0表示管边界。示于图 3 的速度分布可以用下述方程表达 ： （1） ， 其中
表示在半径处流体 2 的元的速度， 其中表示在管 10 的中心处的最大速度， 并且其中 P 表示管 10 的半径。 如果假设在流体 2 内散射元密度均匀， 可以获得在不同速度 处颗粒数目的分布 ， 如图 4 示意性和示例性所示。 在图 4 中可以看出， 由于层流流中速度分布图的抛物线性质， 大多数元以最大速 度行进。此图中的最大速度是由流体元的数目朝向最大速度的急剧增大来表征。在某一速 度的测量的自混合干涉信号的强度与在此速度的流体元的数目成比例。对于光学薄流体， 最大速度点将由自混合干涉信号中的强峰来标记。
优先地通过在管面积上对方程 （1）中定义的速度分布图积分来确定体积流量 。这得到下述方程 ： （2） 。在一实施例中， 距离和速度确定单元 3， 具体而言分析单元 14 调适为根据自混合 干涉信号确定多普勒频率。优先地， 距离和速度确定单元 3 使用下述方程用于根据多普勒 频率确定速度 ：
（3） ， 其中 表示多普勒频率， 表示沿着图 1 中的激光束 7 方向的速度分量， 并且 表示未扰动的激光器 5 的波长。来自在流体 2 内流动的元的反馈在激光器腔体 6 内部 生成具有此多普勒频率的变化的干涉信号， 该干涉信号为自混合干涉信号。 因此， 由检测器 12 检测到的激光输出功率以一频率被调制， 从该频率可以导出流体内散射元的速度。 因而， 通过使用方程 （3） ， 基于所确定的多普勒频率可以确定流体 2 的元的速度。
在一实施例中， 距离和速度确定单元 3， 具体而言分析单元 14 调适为确定所确定 的多普勒频率的最大多普勒频率， 并且通过使用方程 （3） 根据所确定的最大多普勒频率， 确 定流体 2 的元的最大流速度。在此实施例中， 流确定单元 4 调适为确定最大流速度为流属 性。再者， 流确定单元 4 优先地调适为通过使用所确定的最大流速度和方程 （2） 确定体积 流量。
对于透明流体， 单位为任意单位的自混合干涉信号的谱图的实例示于图 5。在此 实例中， 调制设备 15 不调制激光器 5 的频率， 即所示谱图也是多普勒频率的谱图。优先地 通过在自混合干涉信号的多个单独功率谱上求平均而获得该谱图。在示于图 5 的实例中，该频率谱图具有清晰的峰并且在观察到的频率迅速衰减之后， 最大频率是在大约 0.22MHz。 优先地， 放大器被用于测量自混合干涉信号。这种放大器的使用会在自混合干涉信号的谱 图中引入伪像， 其例如为在图 5 中可以看出的在 DC 处的峰。
优先地， 如果激光器 5 的辐射 7 可以达到具有最大流速度的流体的元， 具体而言， 如果激光器 5 的辐射 7 可以达到管 10 的中心， 则该流体被视为透明的。
为了确定依赖于最大流速度的体积流量， 流确定单元 4 优先地包括体积流量函 数， 该体积流量函数定义最大流速度和体积流量之间的关系。该函数优先地由方程 （2） 定 义。流确定单元 4 优先地调适为， 通过使用体积流量函数和最大流速度， 确定体积流量作为 流属性。
设备 1 还包括流宽度确定单元 9， 其用于根据元的所确定的距离确定流的宽度。 如 果激光辐射 7 横穿该流， 该流之外不存在反射或散射元， 具体而言， 在管 10 的边缘之外不存 在反射或散射元， 并且因此将不从管 10 之外返回距离信息。此外， 相对于发射的激光辐射 7 的传播方向， 在管前方不存在流体 2 的散射或反射元， 并且因此将不从该位置返回距离信 息。
因而， 通过确定到距离和速度确定单元 3， 具体而言到距离和速度确定单元 3 的激 光器 5 的最近距离和最大距离， 可以确定流的宽度。例如通过将所确定的宽度与流体 2 的 已知宽度比较， 这个所确定的宽度可以例如用于确定激光辐射 7 是否完全穿透流体 2。 在此 实施例中， 通过将所确定的宽度与管 10 的宽度比较， 可以确定激光辐射 7 是否完全穿透流 体 2。
此外， 所确定的流的宽度可以用于确定管 10 的直径， 该直径对于确定速度分布图 是至关重要的， 并且该直径在可用时使得校准是多余的。 此外， 所确定的流的宽度可以被用 作控制参数以确保所研究的流体是足够透明的， 使得从管 10 的整个截面获得散射功率。
距离和速度确定单元 3 和流确定 4 优先地调适为使用频率调制技术用于确定流体 2 的元的距离。为此， 距离和速度确定单元 3 优先地包括调制设备 15， 其用于调制激光器 5 的频率。调制设备 15 优先地为电流驱动单元， 其优先地施加三角形调制在激光器驱动电流 上。如果激光器 5 为半导体激光器， 如同在此优选实施例这样， 这种电流调制引起发射的辐 射 7 的波长中的相应调制。结果， 当改变注入电流 I 时， 每增加与从激光器 5 到流体 2 相应 元的往返长度适配的波长， 辐射的相位增加 360°。每个 360°相位旋转造成发射的辐射 7 的功率中的一个最小值和一个最大值。这些 “波动” 的数目 由下述方程定义 ： （4） ， 其中 d 表示从激光器到相应流体元的长度。
波长减小具有与散射元远离激光器 5 移动相似的效应， 而波长增大则模仿散射元 朝向激光器 5 移动。如果使用对激光器电流的三角形调制， 波长将周期性减小和增大， 从而 模仿周期性地远离和朝向激光器 5 移动。由检测器 12 测量的功率， 即从激光器腔体 6 耦出 的辐射 13 的强度， 在时间上按此三角形调制的频率变化， 但是现在在该频率上叠加了频率 为 的波动， 该频率 可以通过使用下述方程来确定 ：9与波长变动的函数可以102356297 A CN 102356307说明书7/12 页（5） 。
中的下标表示散射元不移动， 即散射元不具有在辐射 7 方向上的速度分量。这种情况下， 在三角形调制的上升段和下降段期间， 波动频率是相同的。考虑移动的元， 该频 率另外改变了多普勒频率。 当元移动离开时， 在波长减小期间多普勒频率将加到 长增大时从 减去多普勒频率。这可以由下述方程表达 ： （6） 以及 （7） 。
， 而当波在方程 （6） 中，表示在三角形调制的上升段自混合干涉信号的频率， 并且表示在三角形调制的下降段自混合干涉信号的频率。
通过根据下述方程计算频率： （8）并且利用所计算的频率 体而言到激光器 5 的距离。
通过使用方程 （5） ， 可以确定元到距离和速度确定单元 3， 具通过根据下述方程计算： （9）并且利用所计算的频率 度。
通过使用方程 （1） ， 可以确定沿着辐射 7 方向的元的速由于例如通过上述方程 （1） 至 （7） 已知流体的单个元的距离和速度与自混合干涉 信号的频率的依存性， 则例如通过将流体的元的贡献， 即距离和速度的贡献线性地结合到 自混合干涉信号的频率谱图， 也已知流体的若干元的距离和速度与自混合干涉信号的相应 频率谱图的依存性。 这种已知的依存性优先地被用于根据自混合干涉信号的频率谱图确定 距离分布和速度分布。例如， 知晓上述已知的自混合干涉信号谱图与距离和速度分布的依 存性， 可以执行类似蒙特 - 卡罗模拟的模拟， 其中利用不同距离和速度分布模拟自混合干 涉信号的不同谱图， 直到就比如相关性或平方差之和的相似性度量而言， 自混合干涉信号 的模拟的谱图与测量的谱图是相似的。还可以使用拟合过程， 其中距离分布和速度分布被 确定， 使得通过使用上述已知的自混合干涉信号的谱图与距离和速度分布的依存性确定的 自混合干涉信号的谱图被拟合到自混合干涉信号的测量的谱图。 还有可能通过使用上述已 知依存性， 根据自混合干涉信号的测量的谱图解析地计算距离和速度分布。
为了执行上述模拟过程或者上述拟合过程， 自混合干涉信号的谱图与距离和速度 分布的依存性可以看作一模型。除了将流体的单个元的速度和距离的贡献结合到， 具体而 言线性地结合到自混合干涉信号的谱图， 这个模型可以包含在指定速度对流体元的密度的 考虑和 / 或对在流体内的衰减的考虑。下面进一步将更详细描述这些附加考虑。
如果假设不调制电流， 谱图具有形状10， 则如果电流被调制， 谱图在向上102356297 A CN 102356307说明书8/12 页侧翼可以用表示并且在向下侧翼用表示， 其中 a 为比例常数并且 d 为相应元到距离和速度确定单元 3， 具体而言到激光器 5 的距离。散射颗粒生 成的频率是由到激光器的距离连同由于电流调制和颗粒速度引起的激光器频率变化来确 定。距离和速度确定单元 3 优先地调适为如上所述通过模拟、 拟合或解析计算从两个谱图 和
提取距离信息和速度信息。在下文中参考图 6 将更详细描述确定流体中的元的距离和速度的实施例。
图 6 仅仅示出距离和速度确定单元 3 的激光器 5。 激光器 5 的辐射 7 横穿透镜 16， 使得辐射 7 在传播通过管 10 时被聚焦在流体内部。该透镜具有最优化来自流体的自混合 信号的功能。
激光器 5 及其光学元件 16 置于管 10 外部。距离 l 0 为从激光器 5， 具体而言从激 光器腔体的耦出激光器镜到管 10 的中心的长度。字母 r 描述从管 10 的中心沿着辐射 7 方 向的相应元 17 的位置。流速度沿着管 10 的中心通常较高， 并且在壁处在 R max 减小到零。管 10 含有流体 2， 该流体具有散射元 17。
假设流体具有层流流， 其中使用一种校正来校正不垂直于管表面的入射角， 由方 程 （1） 给出流体 2 的速度与在管 10 中的位置的函数。
（10） 。如果流体 2 对于辐射 7 不是完全透明的， 来自管 10 中深度 r 的背散射光的数量由 下述给出 （11） ，表示光在流体 2 中的衰减系数。
在下文中假设元 17 的密度在管长度上恒定， 即在流方向上恒定， 并且自混合干涉 信号强度与距离 r 以及与管长度的依存性也恒定。这优先地是指透镜 16 的聚焦作用相当 弱， 即大的焦点深度。
在深度 r 处的元将散射光， 使得其速度沿着辐射 7 方向的其投影分量产生背散射 光中如方程 （3） 中定义的多普勒频率偏移。
如果也应用激光器 5 的电流的调制， 激光器 5 的波长被调制， 这引起背散射光中的 附加频率偏移， 其可以由下述方程描述 ： （12） ， 其中
其中表示附加频率偏移并且表示未扰动的激光器 5 的波长， 并且 l 0 为激光器到管的中心的距离。 在位置 r 处的元 17 形成功率谱， 该功率谱可以由下述方程描述 ： （13） ， 其中 S (f ) 表示功率谱， 具体而言由检测器 12 测量的强度谱图， 并且其中 g (…) 表示自 混合干涉信号的响应函数。注意， 功率谱不具有负频率。因此考虑所得到的频率的绝对值。多普勒频率偏移和调制频率偏移都依赖于在流体 2 中的位置。再者， 该信号在流体内被衰 减。
由检测器 12 测量的信号可以由下述方程描述， 该信号由将光背散射到激光器腔 体 6 内的所有元的贡献组成 ： （14） 。 通过在位置 r 上对方程 （13） 积分， 得到方程 （14） 。
如果使用三角形调制， 在向上和向下斜坡期间每次的波长变化是恒定的， 只是频 率偏移具有相反符号。
距离和速度确定单元 3 优先地调适为在三角形调制的向上侧翼并且单独地在三 角形调制的向下侧翼采集自混合干涉信号。根据这些自混合干涉信号， 可以针对向上侧翼 计算一功率谱并且可以针对向下侧翼计算一功率谱。优先地， 该距离和速度确定单元调适 为分别在若干向上侧翼上和若干向下侧翼上求平均， 从而提高功率谱的信噪比。
对于向上侧翼， 功率谱由下述方程定义 ：
（15） 。
对于向下侧翼， 功率谱可以由下述方程定义 ： （16） 。通过对响应函数 g (… ) 进行拟设， 例如利用在指定速度的流体元的密度通过进行 高斯函数的卷积， 可以将方程 （15） 和 （16） 拟合到针对向上侧翼和向下侧翼由检测器 12 测 量的功率谱。拟合过程含有调适拟设响应函数中的拟合参数。如果拟合参数为最大流速度 和 ， 这些参数由此拟合过程确定。最大流速度可以用于例如依据方程 （2） 确定体积流量。
优先地， 距离和速度确定单元 3 调适为选择调制单元 15 的调制频率， 使得功率谱 不通过零频率， 因为这会产生在检索速度和距离分布图中的一些不定性。 具体而言， 调制频 率优先地选择为使得多普勒频率和调制频率不具有相同值。这意味着， 如果已知多普勒频 率只能在某一频率范围内， 则调制频率优先地选择为使得其不在此频率范围内。
通过拟合可以将功率谱分离为第一部分和第二部分， 该第一部分仅仅依赖于
， 的距离分布。
对应于元的速度分布， 并且第二部分仅仅依赖于，对应于元对于流体在光学上是厚的情形， 光将没有足够深地穿透到流体内而达到具有最 大速度的元。这种情况下使用最大流速度作为拟合参数不再是有用的。然而， 根据向 上和向下侧翼的两个谱图， 可以确定流体元的相应速度和距离， 即位置。谱图的形式由 和 给出， 其中 仅仅由速度分布给出， 并且依赖于 l 0 的部分仅仅为流体元的位置的函数。使用拟合过程， 可以解开这 两种贡献。例如， 两个未知参数为吸收系数和最大流速度。吸收和速度分布图对 r 的依存 性已知并且在拟合过程中使用。拟合曲线被最优化， 使得它们与向上侧翼和向下侧翼的两 个测量的谱图最佳地对应。 通过该拟合， 速度分布与距离的函数， 即速度分布和距离分布同 时被获得。具体而言如方程 （1） 中定义并且示于图 3 的流模型函数优先地被拟合到所确定的 速度和距离分布， 其中根据所拟合的流模型函数可以确定期望流属性。
调制频率优先地选择为使得它不干涉对于自混合干涉信号是令人感兴趣的谱图 的部分， 即， 如已经在上文所述， 调制频率优先地选择为使得功率谱不通过零频率。 此外， 这 种调制的幅度应大得使得， 在三角形调制的向上或向下部分上的不移动对象的检测中可以 发现少许调制周期。
距离和速度确定单元 3， 具体而言分析单元 14 优先地分开对在三角形调制的向上 和向下部分上获得的功率谱积分。两个侧翼的功率谱优先地具有相同形状， 但是分别由于 位置依存性 同地缩放。 流确定单元 4 优先地调适为提供流模型函数， 其定义依赖于流体的元到距离和速 度确定单元 3， 具体而言到激光器 5 的距离的元的速度。 流确定单元 4 优先地还调适为将该 流模型函数拟合到元的所确定的距离和速度以及根据所拟合的流模型函数确定流属性。
流模型优先地为层流流模型， 该模型假设最大流速度位于流的中间并且零流量值 位于流的边缘。这种优选流模型函数示意性和示例性地示于图 3。可以执行流模型函数 到流体的元的所确定的距离和速度的这种拟合， 即使仅仅确定了流体的少许元的距离和速 度。因而， 即使流体在光学上是厚的。例如如果仅仅可以确定这样的元的距离和速度， 这些 元具有到距离和速度确定单元 3 的距离， 该距离对应于介于 -1.0 和 -0.5 的归一化半径， 则 该流模型函数可以被拟合到这些元的距离和速度， 并且基于所拟合的流模型函数例如可以 确定最大流速度并且因而也可以确定体积流量。
在一实施例中， 距离和速度确定单元 3 和流确定单元 4 调适为使得 a） 如果所确定的宽度等于或大于预定义最大速度宽度， 距离和速度确定单元 3 确定自 混合干涉信号的最大频率并且根据所确定的最大频率确定流体的元的最大流速度， 并且流 确定单元 4 确定该最大流速度作为流属性， b） 如果所确定的宽度小于预定义最大速度宽度， 该流确定单元提供定义依赖于流体的 元到该距离和速度确定单元的距离的元的速度的流模型函数， 将该流模型函数拟合到元的 所确定的距离和速度， 并且根据所拟合的流模型函数确定流属性。
预定义最大速度宽度定义流的宽度， 从而允许根据所生成的自混合干涉信号确定 最大流速度， 其中该流的宽度至少必须由流宽度确定单元 9 确定。这个所确定的流的宽度 依赖于流体 2 的光学厚度， 即在流体 2 中辐射 7 的穿透深度。因而， 通过确定依赖于所确定 的流的宽度的流体 2 的流属性， 流属性的确定依赖于流体 2 的光学厚度。如果光学厚度小， 使得激光辐射 7 达到最大速度宽度， 则距离和速度确定单元 3 与流确定单元 4 调适为， 通过 根据自混合干涉信号确定多普勒频率， 通过确定所确定的多普勒频率的最大多普勒频率， 并且通过根据所确定的最大多普勒频率确定流体 2 的元的最大流速度， 确定最大流速度。 如果光学厚度大得使得激光辐射 7 无法达到最大速度宽度， 则流确定单元 4 提供定义依赖 于流体 2 的元到距离和速度确定单元 3 的距离的元的速度的流模型函数， 将该模型函数拟 合到元的所确定的距离和速度， 并且根据所拟合的流模型函数确定流属性。这允许确定依 赖于流体 2 的光学厚度的流体 2 的流属性。
另外优选的是， 流确定单元 4 调适为， 根据自混合干涉信号确定流体 2 的流是层流
13、的原因， 频率轴被不102356297 A CN 102356307说明书11/12 页还是湍流作为流属性。优先地， 流确定单元 4 调适为， 在自混合干涉信号的频率谱图具有混 沌行为时确定该流是湍流， 以及在自混合干涉信号的频率谱图不具有混沌行为时， 具体而 言在自混合干涉信号的频率谱图具有完善 （well-established） 的形状时确定该流是层流。 流确定单元 4 耦合到泵 11， 使得在流确定单元 4 确定流体 2 的流是湍流时， 泵 11 被控制使 得流体 2 的流变为层流。因而， 设备 1 可以被调适使得， 通过使用包括距离和速度确定单元 3、 流确定单元 4 和泵 11 的控制环路， 获得和 / 或维持管 10 内的流体 2 的层流流。
在下文中， 将参考示于图 7 的流程图描述用于确定流体的流属性的方法。
在步骤 101， 距离和速度确定单元 3 确定流体 2 的元 17 到距离和速度确定单元 3 的距离， 并且同时确定元 17 的速度。通过将在激光器腔体 6 内生成的激光辐射 7 引导到流 体 2 用于被流体 2 反射， 生成自混合干涉信号。反射辐射 8 与激光器腔体 6 内的辐射混合， 并且基于所生成的自混合干涉信号确定流体 2 中元 17 的距离和速度。
在步骤 102， 流属性确定单元 4 基于所确定的距离和速度的至少其一确定流体的 流属性。优先地， 基于元 17 的所确定的距离和所确定的速度， 该流属性确定单元确定最大 流速度和体积流量。
流的精确测量在许多不同应用中是重要且关键的任务， 所述应用涵盖了从例如化 学或食品加工的工业工艺， 到类似例如汽车引擎的机器， 直至比如血液输入或注入的医学 应用。 在大多数情形下， 由于所使用的液体或气体的不同性质， 不同应用需要不同解决方案 来精确确定通过某个管的液体或气体的量 （体积流量） ， 该解决方案涵盖机械流量计， 涵盖 热检测器， 直至涵盖超声装置。 相比之下， 依据本发明的用于确定流体的流属性的设备能够 精确测量显著不同性质的介质中的流速度。 传感器， 即该距离和速度确定单元， 是基于激光 器 5 的激光器腔体 6 中的自混合干涉， 该激光器优先地为半导体激光器。依据本发明的用 于确定流体的流属性的设备能够确定对激光辐射 7 具有强衰减的介质中以及透明介质中 的流速度和体积流量。
依据本发明的用于确定流体的流属性的设备可以在非常多的应用， 具体而言在上 述各应用中使用。
如上文已经描述， 使用激光器腔体 6 内部的激光和外部反馈之间的干涉用于生成 自混合干涉信号。激光自混合干涉的原理允许非接触的速度和距离测量。如果激光器 5 瞄 向散射元 17， 小部分的散射光反射到激光器腔体 6 内， 在那里它与强激光场混合。当元 17 的移动具有沿着激光束 7 方向的分量时， 反射光的相位相对于原始激光连续地偏移， 从而 导致激光器腔体 6 中的反馈按照等于多普勒频率的频率而周期性变动。这在上文中被更详 细地解释并且也在通过引用结合于此的论文 “Laser diode self-mixing technique for sensing applications” ， G. Giuliani, M. Norgia, S. Donati 和 T. Bosch, J. Opt. A: Pure Appl. Opt. 4, 283–294, 2002 中解释。
用于确定流体的流属性的设备使用基于自混合干涉原理的激光传感器来优先地 测量液体或气体的层流流， 而不依赖于介质的性质。该设备可以用于测量透明以及散射或 吸收介质的流。 该设备通过将散射元的速度分布的测量与散射元的距离分布的测量结合而 实现这一点。对于具有大的衰减系数的介质， 即流体， 辐射 7 将不穿透流体 2 够深而达到最 大速度的区域。因此， 无法直接确定最大速度。为了解决这个问题， 该设备确定位置分布， 例如流体 2 内散射元 17 的距离分布连同这些元 17 的速度分布。这给出在流体 2 中辐射 7的穿透深度以及测量的自混合干涉信号中的最大速度的位置， 即辐射从其被背散射的元 17 的速度的最大速度的位置。通过使用流模型函数， 具体而言通过使用示于图 3 的流模型函 数， 通过将该流模型函数拟合到最大的所获得的速度和相应位置， 可以将这个最大的所获 得的速度和相应位置用于确定该流中的整体最大速度。
应指出， 相同设备， 具体而言相同激光传感器被用于确定元的速度和距离。
激光器 5 优先地为半导体激光器， 具体而言为垂直腔面发射激光器 （VCSEL） 。
尽管在上述实施例中已经描述了用于确定自混合干涉信号的某种配置， 在其它实 施例中可以使用用于确定自混合干涉信号的其它配置。
尽管在上述实施例中确定在管中的流体的流属性， 在其它实施例中可以确定不在 管中流动的流体的流属性。例如， 可以由依据本发明的设备确定自由地或者在不同于管的 通道或腔体中流动的流体的流属性。
本领域技术人员在实践所要求保护的发明时， 通过研究附图、 公开内容和所附权 利要求， 可以理解和达成对所公开实施例的其它变动。
在权利要求中， 词语 " 包括 " 不排除其它元或步骤， 并且不定冠词 " 一 " 或 " 一个 " 不排除多个。 单个单元或装置可以完成在权利要求中列出的若干项的功能。 在互不相同的从属 权利要求中提到某些措施的纯粹事实并不表示不能有利地使用这些措施的组合。
由一个或若干单元或装置执行的计算和 / 或确定， 比如确定流体中的元的距离和 速度或者比如确定最大流速度或体积流量， 可以由任何其它数目的单元或装置执行。所述 计算和 / 或确定和 / 或依据用于确定流体的流属性的上述方法对用于确定流体的流属性的 设备的控制可以被实施为计算机程序和 / 或专用硬件的程序代码装置。
计算机程序可以存储 / 发布在与其它硬件一起或者作为其它硬件的部分被供应 的合适介质上， 诸如光学存储介质或固体介质， 但是也可以在其它形式中发布， 诸如经由互 联网或者其它有线或无线的远距离通信系统。
权利要求中的任何附图标记不应解读为限制范围。