用于气体和液体应用的流量测量系统和方法
阅读：1018发布：2020-09-24
专利汇可以提供用于气体和液体应用的流量测量系统和方法专利检索，专利查询，专利分析的服务。并且本 发明 公开一种用于测量在含有 流体 的 导管 内传播的多个 声波 的 叠加 的系统，所述系统具有沿着所述导管的壁大致平行于流动方向 定位 的多个换能器。所述系统包括用于对多个声波在它们在所述导管内传播时的叠加进行建模的装置。，下面是用于气体和液体应用的流量测量系统和方法专利的具体信息内容。
1.一种用于测量在含有流体的导管内传播的多个声波的叠加的系统，所述导管具有纵向长度，其中所述流体在流动方向上流动穿过所述导管，并且其中所述多个声波中的至少一个声波由流动穿过所述导管的所述流体生成，所述系统包括：
多个换能器，其沿着所述导管的所述纵向长度的至少一部分大致平行于所述流动方向定位，其中所述多个换能器中的每个换能器包括用于感测在所述流动方向上流动穿过所述导管的流体的速率的装置，其中所述多个换能器中的每个换能器还包括用于感测所述多个声波传播穿过所述流体的速度的装置，并且其中所述多个换能器中的每个换能器被定位在相对于所述导管上的预定参考点的间隔位置，所述间隔位置对应于所述换能器与所述参考点之间的纵向距离；
用于对所述多个声波在它们在所述导管内传播时的所述叠加进行建模的装置，其中用于对所述叠加进行建模的所述装置与所述多个换能器中的每个换能器通信，并且其中用于对所述叠加进行建模的所述装置被配置来：
从每个换能器接收时域信号；
给所述时域信号分配位置值，所述位置值指示生成所述时域信号的所述换能器的所述间隔位置；并且
存储时域信号阵列及其对应的位置值，
其中，在存储选定数量的时域信号后，用于对所述叠加进行建模的所述装置被配置来处理所述时域信号阵列，以便产生所述多个声波在它们在所述导管内传播时的相位的模型，并且
其中使用以下波动方程来对所述多个声波的所述叠加进行建模
和
其中：
Ψ＝从右到左(RL)或从左到右(LR)移动穿过所述导管的所述波的相位，rw＝所述导管的半径，
z＝相对于所述导管内的所述流体流动方向测量的所述波的水平位置，t＝自生成所述波以来经过的时间，
N＝所述声波的振幅，
J0＝第一类的圆柱形贝塞尔函数，
Mcal＝校准系数，
f＝所述声波的频率，
c＝所述导管内的所述声波的速度，并且
k＝与所述声波相关的波数。
2.如权利要求1所述的系统，其中用于对所述叠加进行建模的所述装置被配置来:
生成与每个换能器相关的频域信号；
形成用于所述多个换能器的二维时域和空间域矩阵，所述矩阵包括与每个相应的换能器相关的时域和频域信号；
对所述二维时域和空间域矩阵执行二维快速傅立叶变换(2DFFT)；并且使所述2DFFT的结果与由所述波动方程建模的结果关联。
3.如权利要求2所述的系统，其中用于对所述叠加进行建模的所述装置被配置来识别所述2DFFT的所述结果与由所述波动方程建模的所述结果之间的最大关联。
4.如权利要求3所述的系统，其中，基于所述2DFFT的所述结果与由所述波动方程建模的所述结果之间的所述最大关联，用于对所述叠加进行建模的所述装置被配置来计算在所述导管内在所述流体流动方向上的声音传播以及逆着所述流体流动方向的声音传播。
5.如权利要求4所述的系统，其中，基于所述导管内的所述计算的声音传播，用于对所述叠加进行建模的所述装置被配置来计算所述导管内的流体流动速率和音速中的至少一个。
6.如权利要求1所述的系统，其中与每个换能器相关的所述时域信号指示在所述流动方向上流动穿过所述导管的所述流体的感测速率以及所述声波传播穿过所述流体的感测速度。
7.如权利要求2所述的系统，其中用于对所述叠加进行建模的所述装置还被配置来：
对从每个相应的换能器获取的所述时域信号执行一维快速傅立叶变换(1DFFT)；并且生成针对每个相应的换能器的归一化的时域信号，
其中所述归一化的时域信号被包括在所述二维时域和空间域矩阵中。
8.如权利要求1所述的系统，其中多个声波包括第一声波和第二声波。
9.如权利要求1所述的系统，其中所述多个换能器包括2个换能器。
10.如权利要求1所述的系统，其中所述多个换能器包括至少64个换能器。
11.如权利要求1所述的系统，其中所述多个换能器包括至少128个换能器。
12.如权利要求1所述的系统，其中所述多个换能器中的每个换能器与相邻换能器纵向间隔开分离长度。
13.如权利要求12所述的系统，其中所述分离长度范围是约0.5英寸至约3.5英寸。
14.如权利要求12所述的系统，其中所述分离长度范围是约1.0英寸至约3.0英寸。
15.如权利要求12所述的系统，其中所述分离长度范围是约1.5英寸至约2.5英寸。
16.如权利要求1所述的系统，其中用于对所述多个声波的所述叠加进行建模的所述装置包括具有处理器的计算机。
17.如权利要求16所述的系统，其中所述处理器包括用于处理从所述多个换能器接收的数据集的模数转换器，并且其中所述模数转换器以采样率处理所述数据集。
18.如权利要求1所述的系统，其中所述导管包括外表面，并且其中所述多个换能器能够可拆卸地安装在所述导管的所述外表面上。
19.如权利要求18所述的系统，其中所述多个换能器互相连接以形成大致线性的换能器阵列。
20.一种用于测量在含有流体的导管内传播的多个声波的叠加的方法，其中所述流体在流动方向上流动穿过所述导管，并且其中所述多个声波中的至少一个声波由流动穿过所述导管的所述流体生成，所述方法包括：
提供多个换能器，所述多个换能器沿着所述导管的纵向长度的至少一部分大致平行于所述流动方向定位，其中所述多个换能器中的每个换能器包括用于感测在所述流动方向上流动穿过所述导管的流体的速率的装置，其中所述多个换能器中的每个换能器还包括用于感测所述多个声波传播穿过所述流体的速度的装置，并且其中所述多个换能器中的每个换能器定位在相对于所述导管上的预定参考点的间隔位置，所述间隔位置对应于所述换能器与所述参考点之间的纵向距离；
从每个相应的换能器接收时域信号；
给所述时域信号分配位置值，所述位置值指示生成所述时域信号的所述换能器的所述间隔位置；
存储时域信号阵列及其对应的位置值；以及
处理所述数据集阵列，以便产生所述多个声波在它们在所述导管中传播时的所述叠加的模型，
其中使用以下波动方程来对所述多个声波的所述叠加进行建模
和
其中：
Ψ＝从右到左(RL)或从左到右(LR)移动穿过所述导管的所述波的相位，rw＝所述导管的半径，
z＝相对于所述导管内的所述流体流动方向测量的所述波的水平位置，t＝自生成所述波以来经过的时间，
N＝所述声波的振幅，
J0＝第一类的圆柱形贝塞尔函数，
Mcal＝校准系数，
f＝所述声波的频率，
c＝所述导管内的所述声波的速度，并且
k＝与所述声波相关的波数。
用于气体和液体应用的流量测量系统和方法
[0001] 相关申请的交叉引用
[0002] 本申请要求2013年3月14日提交的标题为“Flow Measurement Systems and Methods for Gas and Liquid Applications,”的美国专利申请号13/828,942的提交日期的利益，所述申请以引用的方式整体并入本文。
[0003] 领域
背景技术
[0005] 存在用于通过测量传播声波的上场速度和下场速度来测量穿过导管的流量的各种已知技术。然而，除非上场声波可与下场声波分离，否则这些已知方法是无效的。
[0006] 最近，已提出稍微不同的方法来解决常规流量测量技术的限制，以对导管内的声波进行建模。美国专利号7,725,270中描述了一种这样的方法。根据所述方法，假设测量的声波是以线性方式传播的自由的(未束缚的)声波。因此，所述方法忽视了这样的事实，即声波实际上受导管束缚并且声波在三维方向上进行传播。通过简化对声波的建模，仅可产生声波流动的近似模型。
[0007] 因此，相关领域中需要能够精确地测量导管内流体流动速率和导管内声音速度的系统和方法。另外，相关领域中需要能够在考虑到导管的存在的同时，精确地对传播穿过导管内的流体的声波的相位进行建模的系统和方法。此外，相关领域中需要能够精确地对多个声波在它们传播穿过导管内的流体时的叠加进行建模的系统和方法。
[0008] 概述
[0009] 在一个实施方案中，本发明包括一种用于测量在含有流体的导管内传播的多个声波的叠加的系统。在所述实施方案中，导管具有外表面，并且流体在流动方向上流动穿过导管。在一个方面，多个声波包括至少两个声波。
[0010] 在一个方面，用于测量多个声波的叠加的系统包括沿着导管的外表面大致平行于流动方向定位的多个换能器。在另外的方面，用于测量多个声波的叠加的系统包括用于对多个声波在它们在导管内传播时的叠加进行建模的装置。
[0011] 在另一方面，用于对多个声波的叠加进行建模的装置被配置来从每个换能器接收指示流体流动穿过导管的速率以及声波传播穿过流体的速度的数据集。在存储选定数量的数据集后，用于对多个声波的叠加进行建模的装置被配置来处理数据集阵列，以便产生多个声波在它们在导管内传播时的叠加的模型。在一个方面，用于对多个声波的叠加进行建模的装置通过对数据集阵列执行一维快速傅立叶变换来处理数据集阵列。在这方面，用于对多个声波的叠加进行建模的装置可使用一维快速傅立叶变换的结果来生成针对每个相应换能器的归一化的时域信号。随后，用于对多个声波的叠加进行建模的装置可基于数据集阵列和每个换能器的归一化的时域信号产生二维时域和空间域矩阵。随后，用于对多个声波的叠加进行建模的装置可对二维时域和空间域矩阵执行二维快速傅立叶变换。随后，可使二维快速傅立叶变换的结果与对多个声波的叠加进行建模的波动方程关联。基于所述关联，可确定顺着和逆着流体流动的声波的速度。基于顺着和逆着流体流动的声波的已确定的速度，可确定导管内的流体流动速率和音速。
[0013] 在参考附图的详细描述中，本发明的优选实施方案的这些和其他特征将变得更加显而易见，其中：
[0014] 图1是如本文所述的沿着导管的外表面定位的多个换能器的透视图。
[0015] 图2是如本文所述的沿着导管的外表面定位的多个换能器的放大透视图。
[0016] 图3是导管的横截面侧视图，其示出多个换能器中的相应的换能器之间的间距。图3还示意性示出如本文所述的处理器与多个换能器之间的通信。
[0017] 图4示出如本文所述的确定在多个换能器中的相应的换能器处的归一化的时域信号的示例性方法。
[0019] 图6示出如本文所述的使二维快速傅立叶变换数据与用于从左到右和从右到左传播的波的流动模型方程关联的示例性方法。图6还示出使用关联的结果来确定导管内顺着和逆着流体流动的声音传播、流体流动速率和流体内的音速。
[0020] 图7-9示出本文所述的对从左到右的声波传播进行建模的方法的实验应用。图7示出声波的二维快速傅立叶变换。图8示出声波的理论曲面图。图9示出图7的等值线图。
[0021] 图10-12示出本文所述的对从右到左的声波传播进行建模的方法的实验应用。图10示出声波的二维快速傅立叶变换。图11示出声波的理论曲面图。图12示出图10的等值线图。
[0022] 图13-15示出本文所述的对在两个方向上(从左到右和从右到左)传播的声波传播进行建模的方法的实验应用。图13示出声波的二维快速傅立叶变换。图14示出声波的理论曲面图。图15示出图13的等值线图。
[0023] 图16-17使用本文所述的方法示出对应于图13的二维快速傅立叶变换数据的关联信息。图16是曲面图，而图17是等值线图。
[0024] 图18-19使用本文所述的方法示出与具有5.5m/s流动速率和347m/s音速的流体相关的关联信息。图18是曲面图，而图19是等值线图。
[0025] 图20-21使用本文所述的方法示出与具有9.0m/s的流动速率和347m/s音速的流体相关的关联信息。图20是曲面图，而图21是等值线图。
[0026] 详述
[0027] 可通过参考以下详细描述、实例、图和权利要求及其之前和之后的描述而更容易地理解本发明。然而，在公开和描述本设备、系统和/或方法之前，应理解，除非另外说明，否则本发明不限于所公开的具体设备、系统和/或方法，因而当然可以变化。还应理解，本文所使用的术语仅出于描述特定方面的目的，且并不意图加以限制。
[0028] 本发明的以下描述被提供为本发明的当前已知的最好的实施方案的能实现的教导。为此，相关领域的技术人员将认识并了解到可对本文中描述的本发明的各种方面做出许多改变，而仍然获得本发明的有益结果。还将显而易见的是，本发明的一些希望的益处可通过选择本发明的一些特征而不利用其他特征来获得。因此，本领域的从业人员将认识到对本发明进行许多修改和改编是可能的，并且在某些情况下甚至可能是希望的并且是本发明的一部分。因此，以下描述是作为说明本发明的原理而提供的并不在其限制之中。
[0029] 除非上下文另外明确规定，否则如通篇所使用的单数形式“一个(种)”和“所述”包括多个指示物。因此，例如，除非上下文另外说明，否则提及“换能器”可包括两个或更多个此类换能器。
[0030] 在本文中，范围可表达为从“约”一个特定值，和/或至“约”另一个特定值。当表达这种范围时，另一方面包括从一个特定值和/或至另一个特定值。类似地，当通过使用先行词“约”将值表达为近似值时，将理解的是特定值形成了另一方面。应进一步理解，所述范围的每个端点无论是相对于另一个端点还是独立于另一个端点都是有意义的。
[0031] 如本文所使用的，术语“任选的”或“任选地”意指随后描述的事件或情形可能发生也可能不发生，并且所述描述包括所述事件或情形发生的情况和不发生的情况。
[0032] 如本文所使用的，术语“流体”可指气体、液体及其混合物。
[0033] 在一个实施方案中并参考图1-3，描述了用于测量含有流体的导管10内的至少一个声波的传播的系统100。在示例性方面，至少一个声波可包括多个声波，并且系统100可被配置来测量在导管10内传播的多个声波的叠加。在所述实施方案中并如图3中示出的，导管10可具有内表面12和外表面16。可以设想，导管10的内表面12可限定被配置来接收流体的孔14。还可以设想，流体可在流动方向F上流动穿过导管10的孔14。在一个方面，多个声波可包括两个声波(第一声波和第二声波)。
[0034] 在一个方面，用于测量至少一个声波的传播的系统100可包括用于生成至少一个声波的装置。在这方面，用于生成至少一个声波的装置可包括本领域所公知的用于在流体内产生声波的任何装置。例如但不限于，可以设想，用于生成至少一个声波的装置可包括至少一个扬声器或其他振动设备。然而，如本领域普通技术人员将了解的，可以设想，至少一个声波可由流动穿过导管的流体生成。还可以设想，至少一个声波可包括由随机声音或环境声音产生的声波。
[0035] 在另一方面，用于测量至少一个声波的传播的系统100可包括联接至导管10的多个换能器20，使得多个换能器沿着导管10的外表面12的至少一部分大致平行于流动方向F。在示例性方面，多个换能器20可固定联接至导管10的外表面12。然而，如图3所示，在其他示例性方面，可以设想，多个换能器20可被固定在导管10的壁15内，使得多个换能器定位成与导管的孔14流体连通。在再其他示例性方面，可以设想，多个换能器20可固定联接至导管10的内表面12，使得多个换能器与导管的孔14流体连通。可以设想，多个换能器20沿其定位的导管10的部分可具有用于特定应用的任何合适的纵向长度。
[0036] 在示例性方面，多个换能器中的每个换能器20可被配置来获取指示正常听觉区域(范围从约0到约20kHz)内的声波的数据。在这些方面，多个换能器中的每个换能器20可包括用于感测在流动方向F上流动穿过导管10的流体的速率的装置。在示例性方面，可以设想，用于感测流动穿过导管10的流体速率的装置可包括常规流量传感器。
[0037] 另外，多个换能器中的每个换能器20还可包括用于感测至少一个声波传播穿过流体的速度的装置。在示例性方面，可以设想，每个换能器20可包括用于检测在导管10内接近换能器的声音的麦克风。在其他示例性方面，可以设想，每个换能器20可包括被配置来检测在导管内接近换能器的声音的光声设备。在另外的示例性方面，可以设想，多个换能器20可被配置来检测在导管内穿过导管的壁15的一个或多个声波。
[0038] 在另一方面，多个换能器中的每个换能器20可相对于流动方向F与相邻换能器间隔开。在这方面，多个换能器中的每个换能器20可定位在相对于导管10的外表面16上的预定参考点18的间隔位置。在这方面，间隔位置可对应于换能器20与参考点18之间的纵向距离。
[0039] 可以设想，多个换能器中的每个换能器20可被配置来检测相对于穿过导管的流体流动以线性方式在导管10内传播的声波。还可以设想，多个换能器中的每个换能器20还可被配置来检测除了在导管10内传播的声波外的声波，包括例如但不限于环境声波。
[0040] 如本文讨论的，还可以设想，系统100可用于分析至少一个声波来识别测量导管10内的流动速率和音速所需要的声波。在一个方面，多个换能器20可包括例如但不限于，驻极体换能器、压电换能器、光纤换能器、激光换能器、液体换能器、微电机械系统(MEMS)换能器等等。然而，可以设想，可在本文描述的系统100中使用能够检测声波并将声波转换成电信号的任何换能器。
[0041] 在另外的方面，用于测量至少一个声波的传播的系统100可包括用于对多个声波在它们在导管10内传播时的叠加进行建模的装置。在这方面，用于对叠加进行建模的装置可与多个换能器中的每个换能器20通信。尽管本文关于对多个声波的叠加进行建模进行了描述，但可以设想，用于对叠加进行建模的装置还可用于对导管10内的单个声波的流动进行建模。
[0042] 在另一方面，用于对多个声波的叠加进行建模的装置可被配置来从每个换能器20接收呈时域信号形式的数据集，所述数据集指示流体在流动方向F上流动穿过导管10的感测到的速率以及声波传播穿过流体的感测到的速度。如本文所使用的，术语“时域信号”指由换能器20随着时间的推移而获取的指示传播穿过导管的声波的信号。可以了解的是，时域信号可与空间域信号区别开来。例如，正弦波可根据方程sin(2πft)(时域)在时间上传播或者根据方程sin(2πkx)在位置上传播。
[0043] 在另外的方面，用于对多个声波的叠加进行建模的装置可被配置来将位置值分配给数据集，所述位置值指示生成数据集的换能器20的间隔位置。在再一方面，用于对多个声波的叠加进行建模的装置可被配置来存储数据集阵列及其对应的位置值。可以设想，数据集阵列可用于确定在给定位置处(对应于特定换能器的位置)的多个声波的速度和流体的速率。还可以设想，在存储了选定数量的数据集后，用于对叠加进行建模的装置可被配置来处理数据集阵列，以便产生多个声波在它们在导管10内传播时的叠加的模型。
[0044] 在一个方面，可以设想，可在导管10内的任何位置处生成多个声波中的每个声波，而不影响本文所述的系统100的功能。还可以设想，多个声波中的每个声波的相位、频率和振幅不影响本文所述的系统100的功能。
[0045] 应理解，传播穿过导管的波受导管束缚。因此，如本文所公开的，为了精确地对这种波进行建模，有必要应用精确的束缚条件。因为不能去除环境声音的干扰，常规实践已要求将传播穿过导管的波形的两个不同相位分量分离。然而，如下面更详细解释的，本文所公开的方法和模型不要求从左到右与从右到左传播的波分离。
[0046] 所公开的方法和模型部分基于这样的假定，即传播穿过导管10的声波根据以下方程在三个维度上不受束缚地向内和向外传播：
[0047]
[0048] 其中 表示波动方程Ψ的拉普拉斯算子，
[0049] 表示波动方程关于时间的二阶导数，并且
[0050] c表示传播波的速度。
[0051] 假设波被限制在圆柱形导管中，方程可用圆柱坐标表示。因为从左到右和从右到左的流动仅可在z轴线(平行于流动方向F)上分离，以下方程用于表示其他两个坐标(r和φ)：
[0052]
[0053] 以上方程可被分离成如下：
[0054]
[0055] 如果假定分离的解，其中R(r)表示波的束缚径向分量，Φ(φ)表示圆柱形导管中的波的角旋转，并且P(z,t)表示未束缚的传播波，那么可使用以下方程：
[0056] Ψ(r，z，t)＝R(r)Φ(φ)P(z，t)。
[0057] 现在，波动方程可减少至如下：
[0058]
[0059] 其中m表示在分离变量φ中采用的任意常数。
[0060] 可基于以下方程发现角旋转波的解：imφ
[0061] Φ(φ)＝e 。
[0062] 边界条件命令连续波，因此m＝0、±1、±2等。如以下进一步所述的，对于真实数据可假定m＝0，因为波传播大致没有旋转分量。因此，波动方程可被减少至如下：
[0063]
[0064] 具有以下形式的新常数A可被减小：
[0065]
[0066] 这可被重新排列来进一步分离变量如下：
[0067]
[0068] 如果假定传播波不具有角分量，并且m＝0，那么可以以圆柱形贝塞尔函数的形式重新排列方程：
[0069]
[0070] 所述方程具有第一类的第零阶(由于m＝0)圆柱形贝塞尔函数(J0)的解：
[0071] R(r)＝J0(Ar)。
[0072] 可利用波动方程的剩余部分来处理可从左到右和从右到左传播的波的未束缚部分：
[0073]
[0074] 如前所述，波动方程是通过假定波同时从左到右和从右到左传播(在三个维度上向内和向外)而推导出的。然而，可能将方程分离成从左到右的形式和从右到左的形式以便处理两个波动分量。将从左到右和从右到左的分量组合前之前的方程的基本形式可表示如下：
[0075]
[0076] 因此，偏微分算子可被分成两项，从左到右的项 和从右到左的项随后，可推导出从右到左和从左到右传播的波的并行解。从右到左传播的波可由以下方程表示：
[0077]
[0078] 因此，对于从右到左传播的波，存在以下解：
[0079] PRL(z，t)＝e2πi(ft+kx)，
[0081] 如果扩展从右到左的偏微分方程并且代入从右到左的解，那么得到以下方程：
[0082]
[0083] 在偏微分算子的运算之后，得到以下方程：
[0084]
[0085] 所述方程可被减少至如下：
[0086]
[0087] 进一步的减少揭示A的解：
[0088]
[0089] 其中
[0090] 可以设想，因为贝塞尔函数是偶函数，方程的形式对结果不会产生影响。例如，如下所示，从右到左偏微分算子的A的解可与从左到右偏微分算子的A的解相同。
[0091]
[0092]
[0093] 将单个解代入组合的波动方程揭示从左到右传播的束缚波的解为如下：
[0094]
[0095] 其中N是振幅分量。
[0096] 类似地，从右到左传播的束缚波的解的形式可表示为如下：
[0097]
[0098] 实验测试已确认所述波动方程基本精确地对导管中通过二维声学阵列(如本文所述的多个换能器)获取的声波的时域和空间域进行建模。可添加校准项Mcal以更好地配合真实数据：
[0099]
[0100]
[0101] 随后，可处理导管的壁处的边界条件。假设波动速率可被定义为 并且壁处的径向速率必须为零，可使用以下边界条件：
[0102]
[0103] 可在数字上找到所述零阶贝塞尔函数的第一导数的根。前六个根(ζi)在表格1中，如下。
[0104] 表格1
[0105]i 根(ζi)
0 0
0 0
[0106]1 3.831630
2 7.015578
3 10.173635
4 13.323847
6 16.470579
2 7.015578
3 10.173635
4 13.323847
6 16.470579
[0107] 先前描述表明存在二维快速傅立叶变换(2DFFT)产生的超过一个的对角线。这与已知系统的教导相反，如美国专利号7,725,270所教导的系统，所述系统教导在2DFFT中仅存在与 有关的单个对角线谱。事实上，存在多个对角线谱的可能性，其可根据以下方程得以确定：
[0108] 或
[0109] 假设rw是圆柱形导管的固定半径，以下方程可用于对圆柱形导管中的声音传播进行建模：
[0110]
[0111]
[0112] 在一个方面，多个声波可包括第一声波和第二声波。在这方面，第一声波可在流动方向上传播，并且第二声波可在流动方向相反的方向上传播。可以设想，可使用以上波动方程来对第一声波和第二声波进行建模，所述方程在本文中通常被称为“导管束缚传播分离模型(CBPSM)”方程。
[0113] 在示例性方面，用于对多个声波的叠加进行建模的装置可包括具有处理器的计算机。在这方面，处理器可包括用于处理从多个换能器20接收的数据集的模数转换器。如图3所示，可以设想，计算机的处理器可与多个换能器中的每个换能器20可操作地通信，使得处理器(以及模数转换器)被配置来接收由每个相应的换能器收集的数据集。处理器可通过任何常规通信装置，包括有线和无线电子连接与多个换能器20可操作地通信。在另一方面，处理器(以及，任选地，模数转换器)可以以采样率处理数据集。
[0114] 在一个示例性方面，并且并不意图限制，为了充分使用波动方程来对声波进行建N模，多个换能器20可包括2个换能器，其中N大于或等于3、大于或等于4并且优选地大于或等于5。然而，可以设想，系统100可如本文所公开地仅利用少至两个的换能器而起作用。
在一个方面，多个换能器20可包括32个换能器(N＝5)。在另一方面，多个换能器20可包括64个换能器(N＝6)。在再一方面，多个换能器20可包括128个换能器(N＝7)。在另一方面，多个换能器20可包括256个换能器(N＝8)。可以设想，可通过增加系统中换能器的数量来提高模型的特殊分量的精确度。还可以设想，可通过从多个换能器20获取另外的数据来提高模型的时间分量的精确度。
在一个方面，多个换能器20可包括32个换能器(N＝5)。在另一方面，多个换能器20可包括64个换能器(N＝6)。在再一方面，多个换能器20可包括128个换能器(N＝7)。在另一方面，多个换能器20可包括256个换能器(N＝8)。可以设想，可通过增加系统中换能器的数量来提高模型的特殊分量的精确度。还可以设想，可通过从多个换能器20获取另外的数据来提高模型的时间分量的精确度。
[0115] 在示例性方面，可以设想，可鉴于若干因素确定换能器阵列的大小和尺寸，所述因素包括例如但不限于：处理器针对时域的采样率(SR)，其确定最大频率(fmax＝SR/2)，在所述最大频率下可获取数据集；相邻换能器之间的分离距离(zmin)；以及至少一个波动的传播速度(c)。可以设想，最大波数(kmax＝1/(2zmin))可为相邻换能器之间的距离的倒数。还可以设想，2d-FFT中的主对角线可以在整个域上与流体的音速(cSOS)完美配合。
[0116] 在另外的方面，多个换能器中的每个换能器20可与相邻换能器纵向隔开分离长度(zmin)。在这方面，换能器20的分离长度(zmin)可与处理器(模数转换器)的采样率(SR)和多个声波传播穿过流体的速度(cSOS)相关。采样率、分离长度和声波速度之间的关系可由以下方程表示：
[0117]
[0118] 使用这些关系并参考图3，系统100的分量可用于计算换能器20的适当分离长度22。在一个方面，分离长度22可在约0.5英寸与3.5英寸之间，更优选地在约1.0英寸与
3.0英寸之间，并且最优选地在约1.5与2.5英寸之间。在示例性方面，分离长度22可为约
2英寸。然而，可以设想，系统中存在越多的换能器20，分离长度22可越小。还可以设想，换能器20沿着其定位的导管10的外表面16的部分越短，分离长度22可越小。具体地，分离长度22可被选择性地减小，以容纳另外的换能器20并提高由系统100产生的模型的精确度。在另一方面，可以设想，分离长度22可大致与采样率成反比。因此，在这方面，分离长度22将随着采样率增加而成比例地减小。特别地，如果多个声波包括多个高频声波，那么采样率可增加，从而导致分离长度的减小。另外，随着与导管10内的流体相关的音速降低，分离长度22可成比例地减小。
3.0英寸之间，并且最优选地在约1.5与2.5英寸之间。在示例性方面，分离长度22可为约
2英寸。然而，可以设想，系统中存在越多的换能器20，分离长度22可越小。还可以设想，换能器20沿着其定位的导管10的外表面16的部分越短，分离长度22可越小。具体地，分离长度22可被选择性地减小，以容纳另外的换能器20并提高由系统100产生的模型的精确度。在另一方面，可以设想，分离长度22可大致与采样率成反比。因此，在这方面，分离长度22将随着采样率增加而成比例地减小。特别地，如果多个声波包括多个高频声波，那么采样率可增加，从而导致分离长度的减小。另外，随着与导管10内的流体相关的音速降低，分离长度22可成比例地减小。
[0119] 然而，可以设想，采样率可取决于系统100的部件的性能以及换能器20的数量而更高或更低。另外，还可以设想，可参考导管10内的可辨别的声音的最大频率来确定采样率。在示例性方面，可以设想，采样率可小于每秒约40,000个样本。
[0120] 任选地，在另一方面，多个换能器20能够可拆卸地安装在导管10的外表面16上。可以设想，可使用本领域公知的任何常规装置来将多个换能器20安装在外表面16上，所述常规装置例如像但不限于，抽吸装置、临时粘合剂、夹紧装置、固定装置等等。在这方面，多个换能器20可互相连接以形成大致线性的换能器阵列。可以设想，大致线性的换能器20阵列可允许多个换能器的移动，同时维持相邻换能器之间的所需分离长度22。
[0121] 在使用中，使用本文所述系统的方法可包括将多个换能器20定位成大致平行于如本文所述的流动方向F。在示例性方面，多个换能器20可沿着导管10的外表面16的至少一部分定位。在另外的方面，多个换能器中的每个换能器20可定位在相对于导管10的外表面上的预定参考点的间隔位置。
[0122] 在另一方面，所述方法可包括从每个换能器20接收呈时域信号形式的至少一个数据集，所述数据集指示流体在流动方向F上流动穿过导管10的感测速率，以及声波传播穿过流体的感测速度。在另外的方面，所述方法可包括给每个数据集分配指示生成数据集和时域信号的换能器的间隔位置的位置值。在另一方面，所述方法可包括存储数据集阵列及其对应的位置值。在再一方面，所述方法可包括处理数据集阵列，以便产生多个声波在它们在导管10内传播时的叠加的模型。
[0123] 可以设想，至少一个声波穿过导管的流动可独立于振幅。因此，可以设想，电子/机械振幅校准不需要非常精确。还可以设想，因为每个时域信号(和数据集)由换能器获取，可发生振幅的数值归一化，使得换能器之间的任何差异被去除。
[0124] 在另一方面，在从每个相应的换能器收集时域信号(数据集)之后，处理数据集阵列的步骤可包括对每个相应的时域信号执行一维快速傅立叶变换(1DFFT)，以确定每个频real imag率的真实系数和虚拟系数(C 和C )。在另一方面，使用以下方程，可生成每个换能器的归一化的时域信号：
[0125]
[0126] 其中：
[0127] Tinorm＝在位置(i)处的换能器的归一化的时域信号
[0128] Ti＝在导管内的位置(i)处的给定换能器处测量的时域信号，
[0129] Cjreal＝基于在位置(j)处的给定换能器处测量的时域信号的真实频率系数，[0130] Cjimag＝基于在位置(j)处的给定换能器处测量的时域信号的虚拟频率系数，并且[0131] fj＝在位置(j)处的给定换能器处的频率。
[0132] 在另外的方面，处理数据集(时域信号)阵列的步骤可包括基于生产针对每个相应的换能器的标准时域信号和每个相应的换能器的水平位置形成二维时域和空间域矩阵。在再一方面，处理数据集(时域信号)阵列的步骤可包括对二维时域和空间域矩阵执行二维快速傅立叶变换(2DFFT)。
[0133] 在另一方面，处理数据集阵列的步骤可包括使用波动方程来对多个声波的叠加进行建模
[0134] 和
[0135]
[0136] 其中：
[0137] Ψ＝从右到左(RL)或从左到右(LR)移动穿过导管的波的相位，
[0138] rw＝导管的半径，
[0139] z＝相对于导管内的流体流动方向测量的波的水平位置，
[0140] t＝自生成波动以来经过的时间，
[0141] N＝声波的振幅，
[0142] J0＝第一类的圆柱形贝塞尔函数，
[0143] Mcal＝校准系数，
[0144] f＝声波的频率，
[0145] c＝导管内的声波的速度，并且
[0146] k＝与声波相关的波数。
[0147] 在另一方面，处理数据集阵列的步骤可包括对2DFFT数据和波动方程数据分类，以识别2DFFT数据与波动方程数据之间的最大关联。在这方面，可通过制备2DFFT数据和波动方程数据的三维关联曲面图和/或等值线图来识别最大关联，其中在相对的轴线上提供从左到右和从右到左方向上的声波速度。可以设想，2DFFT数据与波动方程数据之间的最大关联的位置可对应于声波在从左到右方向上的速度等于声波在从右到左方向上的速度的位置。还可以设想，所述最大关联的位置可对应于导管内测量的音速(c)。
[0148] 在示例性方面，可从CBPSM方程生成针对范围为从300m/s到400m/s的从左到右的传播速度，并且类似地，用于具有约1m/s精度的从右到左的传播速度的一组曲面。所述组的曲面可表示约10,000(100x100)个曲面。
[0149] 可以设想，曲面图的每个曲面可通过如本领域所已知的最小平方比较测试来与从多个换能器20获取的数据相关联。随后，可如图16、18和20所示将所得的关联系数对从右到左和从左到右的速度绘图。最大关联可表示数据到CBPSM方程的最佳拟合。可以设想，最大关联可揭示顺着和逆着导管内流体流动的传播速度。
[0150] 当存在导管内流体的非零速率时，那么可以设想，可确定用于从右到左的流动和从左到右的流动的最大关联，从而提供顺着流体流动方向的声音传播(cW)和逆着流体流动方向的声音传播(cA)。使用所述信息，可以设想，随后可使用以下方程来计算导管内流体流动的速率：
[0151] 流率(Flow)＝(cW-cA)/2
[0152] 还可以设想，随后可通过使用以下方程来计算导管内的音速(cSOS)：
[0153] cSOS＝(cW+cA)/2
[0154] 实例
[0155] 在一个实验性实例中，例用以下参数测试如本文所描述的示例性系统：采样率(SR)＝10000Hz；速度(cSOS)＝350m/s；以及导管直径(rw)＝1英寸(0.02540m)。给定这些参数，相邻换能器之间的理想间距根据以下计算方法被计算为0.0350m：zmin＝(350m/s)/(10000Hz)＝0.0350m。然而，为了构造方便，将换能器阵列构建成相邻换能器之间的间距为2.000”(0.05080m)。给定采样率，假设高于5000Hz的频率将不会在导管内含有的流体中传播。
[0156] 考虑到上述设计参数，测试了示例性系统。在范围为500Hz到5000Hz的频率下的声音被从左到右地引入到导管中，并且应用本文所述的方法来获取图7中的以下2DFFT数据。流体的音速是347m/s。使用本文所述的CBPSM波动方程，对在相同频率范围内的347m/s模拟理论曲面。所述理论曲面显示在图2中。如在美国专利号7,725,270中描述的那些，图9中的根产生的方程完美地建模。
[0157] 如图10-12所示，针对在347m/s下移动的从右到左的声波执行类似的实验。在两种情形(从左到右和从右到左)下，明显存在反射波。如图所示，波数轴线超过尼奎斯特极限并进入可被认为是负波数的区域中。
[0158] 图13-15中示出的下组数据，用于从左到右和从右到左传播的声波。图13中的真实2DFFT数据揭示保持数据超过尼奎斯特极限的原因—贝塞尔函数影响正和负波数象限。
[0159] 对图13中的数据执行本文所公开的处理步骤。为了执行这些步骤，使用本文所公开的CBPSM波动方程生成针对从左到右的传播速度和从右到左的传播速度的一组理论曲面。为了所述实践，以1.0m/s为级针对两个方向建立300m/s到400m/s的传播流动范围，表示一组基本的10,000(100x100)曲面。图16和17中示出这些曲面中的每一个相对于收集用于图13中的曲面的2DFFT数据的关联。
[0160] 图13中的数据被收集用于流体的零速率流动。因此，预计两个传播的声波将具有相同的传播速度。图17中可清晰地看到，最大关联发生在两个轴线的347m/s处。
[0162] 从图19，顺着流体流动的传播速度被发现为352m/s，而逆着流体流动的传播速度被发现为341m/s。使用本文所公开的公式，这些参数对应于5.5m/s流动速率和346.5m/s音速，其与输入参数极为一致。
[0163] 在图20和21中，高流量被建模。从图21，顺着流体流动的传播速度被发现为355m/s，而逆着流体流动的传播速度被发现为338m/s。使用本文所公开的公式，这些参数对应于9.0m/s流动速率和346.5m/s音速，其与输入参数极为一致。
[0164] 虽然已经在上述说明书中公开了本发明的若干实施方案，但本领域技术人员理解，受益于上述描述和相关附图中所呈现的教导，本发明所属领域技术人员将能想出本发明的许多修改和其他实施方案。因此应理解本发明不限于上文公开的具体实施方案，并且许多修改和其他实施方案旨在包括在所附权利要求书的范围之内。此外，虽然本文以及所附的权利要求书中采用特定的术语，但是它们仅以一般和描述性意义使用，并且不是为了限制所述发明，或所附的权利要求书。
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 一种基于随机偏微分方程的风速波动特性建模方法 | 2020-05-13 | 13 |
| 基于二阶偏微分波动方程的PML吸收边界条件的方法 | 2020-05-14 | 84 |
| 一种基于波动方程正演的近地表噪音压制方法 | 2020-05-12 | 156 |
| 一种基于波动方程正演的近地表噪音压制方法 | 2020-05-12 | 909 |
| 一种海量数据GPU波动方程逆时偏移成像方法 | 2020-05-12 | 794 |
| 一种基于泰勒展开的波动方程鬼波压制方法 | 2020-05-12 | 599 |
| 一种利用波动方程获取地震波路径的方法 | 2020-05-13 | 603 |
| 一种海量数据GPU波动方程逆时偏移成像方法 | 2020-05-13 | 71 |
| 波动方程偏移角集合的有效计算 | 2020-05-11 | 187 |
| 波动方程偏移直接产生角道集方法 | 2020-05-11 | 833 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
分析报告
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
QQ群二维码
意见反馈
意见反馈