CN106537136A - 虚拟多相流计量和砂检测 - Google Patents
虚拟多相流计量和砂检测 Download PDFInfo
- Publication number
- CN106537136A CN106537136A CN201580036015.6A CN201580036015A CN106537136A CN 106537136 A CN106537136 A CN 106537136A CN 201580036015 A CN201580036015 A CN 201580036015A CN 106537136 A CN106537136 A CN 106537136A
- Authority
- CN
- China
- Prior art keywords
- model
- characteristic vector
- flow
- stream parameter
- multiphase flow
- Prior art date
- Legal status (The legal status 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 status listed.)
- Pending
Links
- 239000004576 sand Substances 0.000 title abstract description 21
- 238000001514 detection method Methods 0.000 title description 7
- 238000000034 method Methods 0.000 claims abstract description 84
- 230000008569 process Effects 0.000 claims abstract description 47
- 238000013528 artificial neural network Methods 0.000 claims abstract description 40
- 238000012545 processing Methods 0.000 claims abstract description 32
- 239000012530 fluid Substances 0.000 claims abstract description 22
- 239000013598 vector Substances 0.000 claims description 76
- 238000012549 training Methods 0.000 claims description 37
- 238000013500 data storage Methods 0.000 claims description 25
- 238000003860 storage Methods 0.000 claims description 9
- 241000208340 Araliaceae Species 0.000 claims 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 2
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 2
- 235000008434 ginseng Nutrition 0.000 claims 2
- 230000005055 memory storage Effects 0.000 claims 1
- 229930195733 hydrocarbon Natural products 0.000 abstract description 9
- 150000002430 hydrocarbons Chemical class 0.000 abstract description 9
- 238000013473 artificial intelligence Methods 0.000 abstract 1
- 238000004422 calculation algorithm Methods 0.000 description 24
- 238000005259 measurement Methods 0.000 description 21
- 238000004458 analytical method Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 10
- 239000004215 Carbon black (E152) Substances 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 7
- 230000007704 transition Effects 0.000 description 7
- 230000008859 change Effects 0.000 description 6
- 239000007787 solid Substances 0.000 description 6
- 238000001228 spectrum Methods 0.000 description 6
- 238000009826 distribution Methods 0.000 description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 5
- 230000000875 corresponding effect Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 239000007789 gas Substances 0.000 description 4
- 239000000203 mixture Substances 0.000 description 4
- 239000003921 oil Substances 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 230000006399 behavior Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 239000000470 constituent Substances 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 241001269238 Data Species 0.000 description 2
- 230000003321 amplification Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000009530 blood pressure measurement Methods 0.000 description 2
- 238000004581 coalescence Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 238000005315 distribution function Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 238000010801 machine learning Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 210000002569 neuron Anatomy 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 238000013179 statistical model Methods 0.000 description 2
- 238000000536 acoustic emission spectrum Methods 0.000 description 1
- 230000005534 acoustic noise Effects 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000019771 cognition Effects 0.000 description 1
- 230000001149 cognitive effect Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000007701 flash-distillation Methods 0.000 description 1
- 125000003976 glyceryl group Chemical group [H]C([*])([H])C(O[H])([H])C(O[H])([H])[H] 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 230000009022 nonlinear effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000011896 sensitive detection Methods 0.000 description 1
- 230000009329 sexual behaviour Effects 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000005309 stochastic process Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/14—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object using acoustic emission techniques
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/10—Locating fluid leaks, intrusions or movements
- E21B47/107—Locating fluid leaks, intrusions or movements using acoustic means
-
- 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/74—Devices for measuring flow of a fluid or flow of a fluent solid material in suspension in another fluid
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Acoustics & Sound (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Environmental & Geological Engineering (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Measuring Volume Flow (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Electromagnetism (AREA)
Abstract
执行虚拟和非侵入性多相计量,用以识别碳氢化合物和其它流体的多相流态、流率、砂存在和其它多相流参数。无源声学检测器系统以声发射信号的形式接收声学流信息，并且数据处理器处理和分类声学模式。使用了统计信号处理方法。使用包括隐马尔可夫模型和人工神经网络的人工智能方法以及自动学习过程来为各种流态和流模式提供声学模型。该计量可以用于井下、顶侧和地面应用。
Description
技术领域
本发明涉及利用声发射测量来自动识别井下、顶侧和地表应用中的碳氢化合物和其它流体的多相流态的虚拟及非侵入性多相计量，并且涉及计量这种流态的流率、砂的存在以及其它参数。
背景技术
两相或更多相的同时流动被称为多相流。多相流的流动行为比单相流的流动行为复杂得多，并且多相流中的流态或流型取决于许多因素，包括：一种流体与另一种流体的相对密度比、流体之间的粘度差以及每种流体的速度(滑移)。术语“流体流动”可以包括油、水、气体和固体(砂)。为优化生产并确定井眼中是否产生了砂，对碳氢化合物流态中的多相流参数和流体中砂的存在的测量是非常重要的。
已经提出了用于非侵入性地测量多相流参数的许多方法。这些参数包括流态、流率、固体成分的存在、各个相的体积比和质量比。已经通过在流体中发送声学/超声频率并分析接收到的声学响应的有源系统实现了一种这样的方法，比如，美国专利No.6,672,131和美国专利No.7,775,125。
美国专利No.5,415,048实现了非侵入性振动响应和流动耦合压力测量的组合以确定流。此外，美国专利No.5,415,048结合压差测量使用了管的特征声频和振幅变化，以获得总质量流率和每个相的质量流率。美国专利No.6,575,043通过在管道的壁中产生声波来描述流动的特性。测量并分析完全在壁内传播的各种声波模式的衰减以确定流中的相的分布。美国专利No.6,412,352使用附接到承载多相流体的管的加速度计。分析由加速度计产生的信号以实现多相流的质量流率的非侵入式测量。
美国专利No.7,562,584涉及一种基于对1Hz至15KHz的频率范围内的多相流的声学特征的机械放大和分析的非侵入性和无源的流体流动测量系统。此外，可以确定定量流动数据和定性数据(比如，警报和状态的变化)，并将其无线地传输至远程位置。美国专利No.5,353,627还利用完全无源的声学检测器装置来确定封闭管线系统中的流态。检测到的声学模式被放大并与已知模式进行比较，以根据其声学指纹来识别流态。分析的频率范围是小于25KHz。
利用声发射来识别各种流态和固体成分的存在的非侵入性方法使用了来自流声学数据的各种参数(比如，信号中的信号振幅、rms值、能量和基频成分)，并且使用了阈值和/或模板匹配技术。这种方法所面临的挑战之一是系统中存在连续和随机的背景声电噪声以及声发射信号的非常低的信噪比(SNR)和随机性。因此，这些方法中的大多数不能在实际情况中，特别是在许多相关因素可以以复杂的方式影响多相流的声学特性的井下环境中，提供对碳氢化合物流态的精确的测量。这些方法也不考虑声发射信号的声学可变性和非稳态性。其它缺陷包括入侵性、高功耗、使用放射源、高成本、高复杂性和对于井下应用的大的物理尺寸。
发明内容
简要地，本发明提供了一种用于确定流管道中的流体的多相流的流参数的新的且改进的设备。该设备包括从流管道中的多相流感测声发射的换能器，以及将所感测的声发射转换成数字声发射信号的转换器。该设备还包括具有存储流管道中的流态数据的声学模型的数据库的数据存储器的计算机。该计算机还包括处理器，所述处理器形成来自数字声发射信号和流态数据的声学模型的流参数的测量值，以确定多相流的流参数的声学模型。该处理器执行计算机实现的以下步骤：将所感测的声发射信号分割成数字声发射段序列；确定数字声发射段序列的数字声发射段的特征向量，并且处理特征向量以确定多相流的流参数的模型。
本发明还提供了一种新的且改进的计算机实现的方法，其利用所述计算机的处理器基于来自多相流的声发射和存储在计算机的数据库中的管道中的流态数据的声学模型来确定流管道中的流体的多相流的流参数。该计算机实现的方法通过将声发射信号分割成数字声发射段序列并确定数字声发射段序列中的每一个的特征向量来完成。随后，处理特征向量以确定多相流的流参数的模型。
本发明还提供了一种新的且改进的数据处理系统，其基于来自多相流的声发射来确定流管道中的流体的多相流的流参数。数据处理系统包括存储流管道中的流态数据的声学模型的数据库的数据存储器。该数据处理系统还包括将声发射信号分割成数字声发射段序列的处理器。该处理器还确定数字声发射段序列中的每一个的特征向量，并且处理特征向量以确定多相流的流参数的模型。
本发明还提供了一种新的且改进的数据存储装置，其具有存储在非暂时性计算机可读介质中的计算机可操作指令，该计算机可操作指令用于使数据处理系统在计算机的处理器中基于来自多相流的声发射和流管道中的流态数据的声学模型来确定流管道中的流体的多相流的流参数。存储在数据存储装置中的指令使得数据处理系统中的处理器将声发射信号分割成数字声发射段序列，并且确定数字声发射段序列中的每一个的特征向量。该指令还使得处理器根据特征向量确定多相流的流参数的模型。
附图说明
图1是根据本发明的用于多相流计量和砂检测的设备的示意图。
图2是根据本发明的多相计量模块的示意图。
图3是根据本发明的一系列事件(例如，声发射)的隐马尔可夫建模拓扑和相关联的高斯分布因子的示意图。
图4是根据本发明的利用隐马尔可夫建模来构建声发射模型的方法的功能框图。
图5是根据本发明的用于基于声发射信号来构建声学模型的数据处理步骤的流程图的功能框图。
图6是根据本发明的用于在多相流计量和砂检测期间优化参数的数据处理步骤的流程图的功能框图。
图7是根据本发明的用于基于声发射信号的多相流计量和砂检测的数据处理步骤的流程图的功能框图。
图8是根据本发明的利用人工神经网络的虚拟多相流计量的处理模块的示意图。
图9是图5的处理模块的人工神经网络架构的示意图。
图10是根据本发明的用于基于声发射信号的利用人工神经网路的多相流计量和砂检测的数据处理步骤的流程图的功能框图。
具体实施方式
声发射
根据本发明，多相流中的声发射被定义为在碳氢化合物油气、水和砂的混合物的流体中和/或表面上发生的声学能量的物理现象。声发射由1KHz至100MHz的宽频范围内的弹性能量的自发性释放而产生，但是大部分释放的能量在1kHz到1MHz的频率范围内。声发射可以从很多源产生，包括：
(a)气泡形成、破裂和聚结；
(b)由涡流和旋涡产生的湍流噪声；
(c)多相流中的液体、气体和固体的相互作用；
(d)由高流旋涡造成的宽带湍流能量；以及
(e)由空化、闪蒸和再循环引起的间歇和瞬态能量变化和波动
在从各种流态的作为气泡尺寸和气泡数量的函数的多相流中发射的声音或较高频弹性能量的研究中已经进行了大量的调查。已经确定，从多相流发射的声能量是多相流内的气泡形成、破裂和聚结以及各相的相互作用的直接结果。此外，声能由于不同多相流态、流率以及流体中固体成分的量而变化。本发明通过感测声能数据和处理这种数据来提供虚拟流建模。
根据本发明，术语“虚拟计量”或“虚拟多相计量”是指如本文所述的在没有任何有源和直接的流测量方法或设备的情况下测量各种流参数的计量方法/技术。如将在下面描述的那样，根据本发明，根据一组无源测量(比如，压力、温度和声发射)来确定流参数。
硬件架构
如图1中所示，根据本发明的用于多相流虚拟建模的设备M包括安装在管、管子或其它流管道22上的作为换能器的声发射传感器或扩音器，如24处所示,通过该管道会发生要被计量的碳氢化合物/水混合物的多相流。流管道22可以是例如用于所谓的井下、顶侧和地表应用的在碳氢化合物的探测、生产或运输中利用的管道。
换能器20接收来源于多相流的具有1kHz到1MHz的频率的多相流声发射信号，以便捕获来自这种流的声发射。可以使用市售的声发射传感器/换能器(比如，来自VallenSystems的AE1045S或来自Mistras集团有限公司的宽带型AE传感器WSA)。应当理解，如果需要的话，也可使用其它市售的换能器。可利用1KHz至2-3MHz的频率范围的传感器。
由声发射传感器20捕获的信号是多相流声学信息和随机背景声学噪声的组合。施加耦合剂以将传感器20与流管道22耦合。通常，使用甘油基耦合剂或油基耦合剂。声发射传感器20将声学信号转换为电信号，使用前端预处理器P的前置放大器26来放大该电信号。来自前置放大器26的表示所感测的声发射的放大的信号在合适的滤波器28中被过滤，并使用高分辨率Σ-Δ模数转换器或A/D转换器30来将其转换成数字化的声发射信号。
来自A/D转换器30的被转换的数字声发射信号由数据处理器D接收。数据处理器D可为已编程的个人计算机(或PC)，或者其可为专用的专用数字信号处理器(或DSP)。数据处理器D可为具有适合的处理和存储能力的任何常规类型的处理器，比如膝上型计算机、平板计算机或其它合适的数字数据处理设备。
数据处理器D包括用于处理、分析和分类声学信号并提供计量结果的多相计量(MPM)软件模块32。原始数据、声学模型和测量结果可以作为数据库被存储在存储器中，用于以后的分析。如34处所示，来自模块32的处理的结果可用于在合适的显示器或绘图仪上的分析，或如将要阐述的那样，可以被传送到无线通信模块以便传送到另一计算机用以研究和分析。
如上所述，如果需要的话，设备M可以包括输出显示器和/或用户界面。例如，在地表或井口安装的情况下，设备M可以包括用于显示计量结果的输出显示器。
如果系统中需要任何修改的话，则设备M还可包括用户界面。示例包括用以访问数据处理器D的用户界面，以便在需要时更新多相计量软件模块32或前端模块40。用户界面还可用于访问和/或更新存储器38中的声学模型42、程序代码36或训练数据90。
处理器D访问来自转换器28的数字化的声发射信号，以执行本发明的可以作为一系列计算机可执行指令来被执行的处理逻辑和方法。指令可以作为程序代码36包含在作为计算机可读介质的存储器38(图2)中。该指令还可以以常规硬盘驱动器、电子只读存储器、计算机磁盘或磁带、光学存储装置或其它适当的数据存储装置的形式进行存储。
本文中的图5、图6、图7和图10中的流程图示出了根据本发明的如在计算机程序软件中实施的处理逻辑的结构。本领域技术人员应当理解，流程图示出了包括根据本发明来运行的集成电路上的逻辑电路的计算机程序代码元件的结构。因此，显而易见的是，本发明在其优选实施例中通过机器组件实现，所述机器组件以指导数字数据处理设备执行与所示的那些功能步骤相对应的功能步骤的序列的形式来呈现程序代码元件。
如将所阐述的那样，多相计量软件模块32是可训练的。根据本发明的一个实施例，隐马尔科夫模型(HMM)用作流态的统计模型的框架。隐马尔科夫模型是输出一系列符号或量的统计模型。隐马尔科夫模型可以自动训练并且是计算可执行的。
虽然多相声学信号是具有连续随机噪声的非线性信号，但是根据本发明，将它们视为短时稳定信号或线性信号。因此，根据本发明，在由传感器20感测的声发射信号被转换成数字格式之后，将其划分或分割成具体持续时间(通常为50ms至200ms)的小段。出于本发明的目的，假设声发射信号在每个这种时间段内是线性的。为每个稳定段计算特征向量，并且使用这些特征向量来创建和训练隐马尔科夫模型。为了解决由于存在连续随机噪声而引起的可变性，使用从实验室流动回路设置和/或现场数据收集的大量声学训练数据来创建和训练每种流型的隐马尔科夫模型。
多相计量模块
多相计量软件模块32(图2)以软件的形式实现，并且接收来自前端模块40的流数据用于预处理和特征计算，并且还接收存储的流态数据的来自综合模型库或数据库42的数据。模型库42中存储的数据采用从实验室和现场数据获得的各种流型的声学模型的形式，包括如在44处示意性示出的多相流态模型、如在46处示出的含砂量以及如在48处示出的流率。多相计量软件模块32被配置为辨识声学模式并识别各种参数(包括流态、流率和固体含量)的解码器。
前端模块40提供特定流段的声学观测值。这些观测值是被视为表示流的声学特征的特征向量。出于训练和检测的目的，从声学信号中提取特征向量。在预处理阶段中，连续声发射声学信号被分割，并被分成50ms至200ms的重叠段，并且在每个段上执行开窗。执行开窗是为了减少边缘处的能量，并且减少每个段的边缘处的不连续性。出于此目的，例如可以使用汉明(Hamming)窗。还可以使用包括海宁(Hanning)窗和布莱克曼(Blackman)窗的其它窗。在预处理之后，为每个段执行特征计算。
重要的声学特征可以包括频带中的频率分布、主频带和主能量。来自多相流的声发射是非稳定过程，因此，标准信号处理技术不适合于分析。倒谱系数也可以用作特征向量。倒谱系数是信号的对数振幅频谱的傅里叶分析的结果。如果对数振幅频谱包含许多规则间隔的谐波，则频谱的傅里叶分析显示与作为基本频率的谐波之间的间隔相对应的峰值。表示重要声学特征的其它参数(比如，小波变换)也可用于计算特征向量。
为了使多相流虚拟建模设备M测量各种流参数和砂存在，为将要遇到的各种流型中的每一个形成声学模型。为了为每种流型建立精确的声学模型(隐马尔可夫模型)，从实验室流动回路获取训练数据，并且使用实际现场。有益且优选的是，获取和累积了大量的这种训练数据。
隐马尔科夫模型本质上是对由一组状态s1、s2、s3、s4和s5(图3)和这些状态之间的跃迁概率a11、a12、a22、a23、a33、a34、a44、a45和a55定义的随机过程的建模，其中，每个状态描述了稳定随机过程，并且从一个状态至另一个状态的跃迁描述了该过程如何随时间改变其特征。每个状态s1、s2、s3、s4和s5具有通常由高斯概率分布函数表示的统计输出分布(如分别在b1至b5处示意性示出)，统计输出分布提供每个观测的特征向量的似然性和分布。优选使用所谓的Baum-Welch算法根据训练数据估算特定流态的隐马尔可夫模型的各种参数。
隐马尔科夫模型被正式定义为：
λ＝(A，B，π)
其中：
A是跃迁阵列，其存储由一个状态至另一个状态的跃迁概率；
B是观测值概率阵列，其存储从状态产生的观测值的概率；并且
π是状态的初始概率阵列。
使用隐马尔可夫模型为每种流型构建声学模型
图5中的流程图H中示意性的示出了根据本发明使用隐马尔可夫模型记录来为每种流型构建声学模型的方法。如步骤50处所示，记录表示特定流型的至少50个(每一个的时长为5-10秒)声发射信号的训练集。可以增加训练数据集的大小和各个声发射信号的长度以提高声学模型的精确度。声发射信号可以从实验室流动回路和/或实际现场获取。根据目标应用，也可以将从地面或井下状况获取的数据存储在存储器中并使用。
在步骤52期间，定义了隐马尔可夫模型的架构。使用通常具有5至20个状态的从左到右型隐马尔可夫模型。如果需要的话，也可以使用其它数量的状态。状态的数量取决于单个声发射信号中的每段的长度和段的数量。也可以使用具有不同拓扑结构(比如，遍历拓扑或并联拓扑)的隐马尔可夫模型。
针对每个训练数据，在步骤54期间使用上述过程计算特征向量。在步骤56期间，初始化各状态的概率。第一状态的初始概率是1，而将剩余状态的初始概率视为0。除了最后一个跃迁被初始化为1以外，跃迁概率中的每一个被初始化为0。
每个状态的统计输出分布是高斯分布。作为步骤58的一部分，计算完整训练声发射训练数据集中的所有特征向量的全局均值和方差。然后在步骤58期间用这个均值和方差来初始化所有高斯概率分布函数(即，所有状态输出分布)。
在初始化隐马尔可夫模型的所有参数之后，优选地使用如下所述的Baum-Welch算法来优化参数(A,B,λ)，并使用训练数据训练隐马尔可夫模型。
参数优化算法
根据本发明，优选地使用Baum-Welch算法来优化未知的隐马尔可夫模型参数(A,B,π)。Baum-Welch算法是广义期望最大化(GEM)算法的特定情况。Baum-Welch算法是基于前向-后向算法的估算方法。下面提供了该算法的简要描述，而详细的数学描述可参见“Amaximization technique occurring in the statistical analysis of probabilisticfunctions of Markov chains”,(Baum1970)(“马尔科夫链的概率函数的统计分析中出现的最大化技术”(Baum，1970年))和“Hidden Markov Models and the Baum–WelchAlgorithm”,IEEE Information Theory Society Newsletter,Dec.2003(“隐马尔科夫模型和Baum-Welch算法”(IEEE信息理论学会通信，2003年11月))。
为了使用训练数据集Otraining(其中，Otraining包含特定流型的大量(≥50个)的声发射特征向量(观测值序列))来构建和训练隐马尔可夫模型参数λa，模型参数λa将使得观测值最接近：
在由流程图C(图6)示意性地示出的序列中执行优化。如步骤60处所示，以上述方式初始化模型参数(A,B,π)。在步骤62期间初始化状态跃迁概率，并且在步骤64期间初始化状态概率。在步骤66期间，初始化各状态的平均值和方差。随后，在步骤68期间，通过模型运行每个观测值序列Otraining(n)，以估算每个模型参数的期望值。
在步骤70中，调整或改变模型参数，使得概率P(Otraining(n)|λa)最大化。随后，针对所有训练数据重复该过程，直到步骤72表示所有模型参数收敛到最佳值为止，并且重复训练过程直到步骤74表示处理了数据集中的所有信号为止。
一旦由此为所有流型创建了声学模型(隐马尔可夫模型)，则该模型可以由自动解码器模块32(图2)使用，以使用统计搜索算法实时地检测和识别多相流声发射信号。
多相流计量和砂检测
通过流程图F(图7)示意性地示出了用于多相流计量和砂检测的方法。多相流声发射信号由传感器20感测。感测的信号在被预放大之后被过滤，并且在A/D转换器30中被转换成数字信号。如步骤80中所示，通过多相计量解码器模块32获取该数字信号。在步骤82中，软件前端40将信号分割，并且在步骤84期间计算每段的特征向量。在步骤86期间，将声发射信号的代表性特征向量的序列传递至解码器32。解码器32也可使用各种流型的声学模型(隐马尔可夫模型)训练数据的综合库90(图2)。
在步骤92期间，使用维特比算法来寻找当前观测的声发射信号最可能归属的(多个)声学模型。在维特比处理中，给定各流型的声学模型库(λ1,λ2…)和接收到的声发射特征向量Oreceived，如下计算针对所有流型给定隐马尔可夫模型时所观测到的特征向量的似然性：
P(Oreceived|λ1)，P(Oreceived|λ2)，P(Oreceived|λ3),…
维特比算法是一种动态规划算法，用于寻找导致观测到的事件序列的被称为维特比路径的最有可能的隐藏状态的序列。下面提供了维特比算法的简要描述。该算法的细节可以在G.D.Forney,Jr.,“The Viterbi Algorithm,”Proc.IEEE,vol.61,pp.268–278,March 1973(G.D.Forney,Jr..“维特比算法”.《Proc.IEEE》，1973年3月第61卷第268–278页)和L.R.Rabiner,“A Tutorial on Hidden Markov Models and SelectedApplications in Speech Recognition,”Proc.IEEE,vol.77,pp.257-286,Feb.1989.(L.R.Rabiner.“隐马尔可夫模型和语音识别中的所选应用的教程”，《Proc.IEEE》，1989年2月第77卷第257-286页)中寻找到。
给定特征向量(观测值序列)Oreceived和隐马尔可夫模型λ1，维特比算法递归地确定模型中的最佳的状态序列Q＝{q1,q2,q3…}。最佳状态序列是最好地解释数据的状态序列，并且可以计算该状态序列的概率，由此计算P(Oreceived|λ1)。该算法基于如下假设：给定观测值序列，最佳路径(或模型中的最佳状态序列)比任何其它路径都更有可能，因此，P(Oreceived|λ1)≈P(Oreceived,δ|λ1),其中δ是最佳路径。在每个步骤中，仅存储一个路径，而丢弃其余路径。通过回溯可以重建最佳状态序列。
对于一个或一个以上声学模型而言，观测的声发射信号的可能值可以是高的。例如，弹状流和低含砂量的声学模型的高概率值可以表示砂含量低的多相流正在缓慢流动(slug)。因此，可以根据解码器输出来确定包括流态、流率范围和砂存在的重要的流参数。
为了减少随机噪声和可变性对实际井下或地表状况下的计量的影响，需要降低训练状况和识别状况之间的失配。这可以通过使用实际现场数据作为训练数据来实现。另外，可以针对特定计量状况来优化声学模型(隐马尔可夫模型)。此外，期望噪声频谱特性比声发射的频谱特性更稳定，因此，可以通过在频谱域中应用频谱减法来补偿噪声。
在步骤92期间以这种方式确定了最可能的声学模型之后，在步骤94期间，所确定的结果可用作为计量结果34(图1和图2)。
使用人工神经网络的声学建模
图8示出了根据本发明的用于多相流建模的另一个实施例，该实施例中使用人工神经网络或ANN方法来开发多相流模型。在如对于图2所述的模块P中的信号预处理和模数转换以及前端模块40中的特征计算之后，人工神经网络或ANN模块100接收数字输入信号。
如下面将描述的那样，ANN模块100包含适当大量的简单处理神经元样处理元件(节点)、以及各元件或各节点之间的适当大量的加权连接。ANN模块100还根据人工神经网络方法获得关于连接的知识的分布式表示，该知识由模块100中的网络通过学习处理而获得。
图8中所示的根据本发明的基于人工神经网络的多相流计量优选地是与用于训练的反向传播算法一起使用的多层前馈神经网络。
图9中提供了人工神经网络模块100的更详细的结构。提供了输入层中的适当数量的节点102，这种节点的数量取决于特征向量的数量和每个向量的长度。使用上述和图5中示出的过程来计算信号段的特征向量104。人工神经网络模块100中还提供了适当数量的输出节点106(图9)。输出节点的数量由流型的数量或要为多相流的模型开发的参数的数量来确定。还提供了隐藏层中的适当数量的节点108，该数量取决于输入节点102和输出节点106的数量。隐藏层节点采用非线性S型(sigmoidal)激活函数神经元的形式。
使用存储在存储器38中的大训练数据集和反向传播算法来训练人工神经网络模块100。训练数据集应包括每种流型的50个以上的声发射信号。一旦人工神经网络模块100被训练，其可用于使用与针对隐马尔可夫模型所述的相同过程的多相流计量。
构建人工神经网络
通过流程图N(图10)示意性地示出了用于为多相流计量和砂检测构建人工神经网络的方法。如110处所示，记录代表特定流型的每一个为5-10秒的至少50个声发射信号的训练集。所记录的声发射信号可以从实验室流动回路和/或实际现场数据获取。可以根据目标应用来使用从地面或井下状况获取的数据。
如步骤112处所示，将信号分割为具有合适的时间长度的段。针对每个训练信号，在步骤114中，使用如上针对步骤54所描述的过程来计算特征向量。在步骤116中，将从步骤114处得到的特征向量提供给人工神经网络模块100。
一旦创建了人工神经网络，作为步骤118的一部分，优选地使用合适的反向传播算法来训练人工神经网络。这种用于模块100的反向传播算法的示例包含在“ParallelDistributed Processing:Explorations in the Microstructure of Cognition”(“并行分布式处理：认知微观结构的探索”)(Rumelhart和McClelland,1986)中。在步骤118期间，从训练信号计算的特征向量104还被提供为训练输入，并且其相应的流型作为输出呈现给人工神经网络(监督式学习)。利用每种流型的训练集来重复该处理。更详细的描述可在Ethem,“Introduction to Machine Learning”(2nded.).Cambridge,Mass.:MIT Press,2010(Ethem.“机器学习简介”(第二版).剑桥，马萨诸塞州：MIT出版社，2010)中找到。
如其名称所示，在反向传播算法中，误差从输出节点反向传播至输入节点。从技术上讲，反向传播计算网络的误差关于网络的可修改权重的梯度。
利用人工神经网络建模进行计量
在利用人工神经网络建模进行的计量中，声发射换能器被激活，并且从管或管道22中的要被计量的碳氢化合物/水混合物的多相流中接收多相流声发射信号。如上所述，声发射信号在硬件前端模块P中被预防大、过滤并被转换成数字信号。
随后，通过多相计量软件模块32接收或收集来自模块P的数字信号。数字信号的处理在图10中示意性地示出。如步骤112处所示，软件前端40划分或分割该信号，并且计算每段的特征向量。
如步骤116处所示，声发射信号的代表性特征向量的序列从数据存储器或存储器作为输入被传递至人工神经网络模块100。在步骤118期间，执行人工神经网络处理并确定计量结果和流参数。如步骤120处所示，形成计量结果的输出，并可从输出节点106处获得计量结果，用以显示分析和评估。还将计量结果存储在数据存储器中以供进一步使用和分析。
相敏检测前端
虽然利用滤波器28进行了一定量的预滤波，不过图1所示的设备A利用了所描述的完整多相流声发射的宽带测量。在声发射频谱具有高水平的信噪比的情况下，可能需要包括高Q值测量。这可以通过以下方式之一执行：
(a)具有可调谐的高Q值滤波器的单锁定放大器信道，其中，所测量的电压水平表示在滤波器的中点频率处的噪声的振幅；或者
(b)每个锁定放大器被调谐到不同的频率的锁定放大器阵列，每个锁定放大器输出与滤波器的中点频率处的声音的振幅相对应的DC电压。
在这两种情况下，来自锁定放大器前端的输出提供了时间依赖性多相流声发射的傅里叶变换的量化的近似值。这可被内插并且被直接处理为声音频谱，并且可以执行逆傅里叶变换以重建时间依赖性行为或仅仅直接进行处理。
在训练阶段期间，锁定放大器滤波器的固定频率的选择可以被调谐到声音频谱内的相关频率，以最大化信号恢复并最小化噪声和其他机制(比如，混响)。
可以使用标准集成电路部件来组装锁定放大器，并且可以通过使用有限脉冲响应滤波器或无限脉冲响应滤波器的前端数字信号处理，或者通过可调谐模拟滤波器来实现可调谐频率滤波器。
应当注意，由于在用作为处理器D的PC/DSP中实现了信号处理和模式识别技术，所以基于隐马尔可夫模型或人工神经网络的计量的使用不影响系统的硬件架构。因此，根据本发明的多相流计量可以具有三种可能的配置：
(a)以基于隐马尔可夫模型的声学模型为基础的多相计量系统；
(b)以基于人工神经网络的声学模型为基础的多相计量系统；或者
(c)以基于隐马尔可夫模型和人工神经网络二者的声学模型为基础的多相计量系统。
在系统的另一个实施例中，设备M可以与一个或多个传感器测量(包括差压测量)相集成，以提高数据的精确度。设备M还可以与无线通信模块相集成，用以实现以下功能：
(a)在地面或井口附近应用的情况下，将数据发送至中央位置，并且实现远程操作；或者
(b)在井下应用的情况下，将数据发送至地面或井口。
因此，可以看出本发明提供了一种多相流计量解决方案，其能够在作为顶侧离岸位置或者岸上位置的地面应用中进行安装，并且还可部署在井下作为永久或可回收系统的一部分。该系统特性允许计量系统的紧凑型封装。
因此，本发明利用不同的方法，使用先进的信号处理技术和自动学习过程对多相流中的声学变化进行统计建模，以实现精确的多相计量。本发明以从多相流产生的声发射、在碳氢化合物油气、水和砂的混合物的流体的内部和/或表面发生的弹性能量的释放为基础。
本发明成本低廉、非放射性、超低功耗且尺寸小。本发明是非侵入式的，因为其不需要穿透管壁，因此不阻碍流动。本发明也不需要产生声学或射频(RF)激励信号来穿过流体和/或管壁。本发明已经得到充分的描述，使得具有对该主题的普通知识的人可以再现并获得本文的发明中提及的结果。然而，本发明主题的技术领域的任何技术人员可以进行未在本文请求中描述的修改，以将这些修改应用于确定的结构或其制造过程中(其需要所附的权利要求中所要求保护的主题)；这种结构应当包括在本发明的范围内。
应当注意和理解的是，在不脱离如所附权利要求中阐述的本发明的精神或范围的情况下，可以对上面详细描述的本发明进行改进和修改。
Claims (46)
1.一种用于确定流管道中的流体的多相流的流参数的设备，包括：
(a)换能器，其从所述流管道中的所述多相流感测声发射；
(b)转换器，其将所感测的声发射转换成数字声发射信号；
(c)计算机，其包括存储所述流管道中的流态数据的声学模型的数据库的数据存储器；并且
(d)所述计算机还包括处理器，所述处理器根据所述数字声发射信号和所述流态数据的声学模型来形成流参数的测量值，以确定所述多相流的流参数的声学模型，所述处理器执行计算机实现的以下步骤：
(1)将所感测的声发射信号分割成数字声发射段序列；
(2)确定所述数字声发射段序列的各数字声发射段的特征向量；
(3)处理所述特征向量以确定所述多相流的流参数的模型。
2.根据权利要求1所述的设备，还包括在所述声学模型的数据库中存储实际多相流状况的数据库的所述数据存储器。
3.根据权利要求2所述的设备，其中，所述处理器在处理所述特征向量以确定流参数的模型时接收来自所述数据库的实际多相流状况数据作为输入。
4.根据权利要求1所述的设备，其中，所述处理器在处理所述特征向量以确定流参数的模型时执行隐马尔可夫建模。
5.根据权利要求1所述的设备，其中，所述处理器在处理所述特征向量以确定流参数的模型时执行以下步骤：
基于所确定的所述数字声段序列的流向量和所存储的所述流管道中的流态数据的声学模型来确定所述多相流的流参数的模型。
6.根据权利要求5所述的设备，其中，所述处理器在确定流参数的模型时基于所确定的所述数字声段序列的流向量和所存储的所述流管道中的流态数据的声学模型来确定所述多相流的流参数的最可能的模型。
7.根据权利要求1所述的设备，其中，所述处理器在处理所述特征向量以确定流参数的模型时执行人工神经网络建模。
8.根据权利要求1所述的设备，其中，所述处理器在处理所述特征向量以确定流参数的模型时执行以下步骤：
(a)接收所述特征向量来作为人工神经网络处理的输入状态；
(b)基于所述输入状态执行所述人工神经网络处理以确定所述多相流的模型的流参数；以及
(c)提供所确定的所述多相流的模型的流参数来作为输出状态。
9.根据权利要求1所述的设备，其中，所述处理器在处理所述特征向量以确定流参数的模型时执行隐马尔可夫建模和人工神经网络建模。
10.根据权利要求1所述的设备，其中，所述处理器还形成用于执行处理所述特征向量的步骤的训练模型。
11.根据权利要求10所述的设备，其中，所述处理器还将所形成的训练模型存储在所述存储器中。
12.根据权利要求1所述的设备，其中，所述处理器还提供所确定的所述多相流的流参数的模型用于显示。
13.一种计算机实现的方法，其利用所述计算机的处理器基于来自多相流的声发射和存储在所述计算机的数据库中的管道中的流态数据的声学模型来确定流管道中的流体的多相流的流参数，所述方法包括以下计算机处理步骤：
(a)将所述声发射信号分割成数字声发射段序列；
(b)确定所述数字声发射段序列中的每一个数字声发射段的特征向量；
(c)处理所述特征向量以确定所述多相流的流参数的模型。
14.根据权利要求13所述的计算机实现的方法，其中，计算机数据存储器存储实际多相流状况的数据库，并且其中，基于来自所述数据库的实际多相流状况数据来执行处理所述特征向量以确定流参数的模型的步骤。
15.根据权利要求13所述的计算机实现的方法，其中，处理所述特征向量以确定流参数的模型的步骤包括隐马尔可夫建模。
16.根据权利要求13所述的计算机实现的方法，其中，处理所述特征向量以确定流参数的模型的步骤包括以下步骤：
基于所确定的所述数字声段序列的流向量和所存储的所述流管道中的流态数据的声学模型来确定所述多相流的流参数的模型。
17.根据权利要求13所述的计算机实现的方法，其中，处理所述特征向量以确定流参数的模型的步骤包括以下步骤：
基于所确定的所述数字声段序列的流向量和所存储的所述流管道中的流态数据的声学模型来确定所述多相流的流参数的最可能的模型。
18.根据权利要求13所述的计算机实现的方法，其中，处理所述特征向量以确定流参数的模型的步骤包括人工神经网络建模。
19.根据权利要求13所述的计算机实现的方法，其中，处理所述特征向量以确定流参数的模型的步骤包括以下步骤：
(a)接收所述特征向量作为人工神经网络处理的输入状态；
(b)基于所述输入状态执行所述人工神经网络处理，以确定所述多相流的模型的流参数；以及
(c)提供所确定的所述多相流的模型的流参数作为输出状态。
20.根据权利要求13所述的计算机实现的方法，其中，处理所述特征向量以确定流参数的模型的步骤包括隐马尔可夫建模和人工神经网络建模。
21.根据权利要求13所述的计算机实现的方法，还包括形成用于执行处理所述特征向量的步骤的训练模型的步骤。
22.根据权利要求21所述的计算机实现的方法，还包括将所形成的训练模型存储在所述计算机的存储器中的步骤。
23.根据权利要求13所述的计算机实现的方法，还包括提供所确定的所述多相流的流参数的模型以用于显示的步骤。
24.一种数据处理系统，其基于来自多相流的声发射来确定流管道中的流体的多相流的流参数，所述数据处理包括：
(a)数据存储器，其存储所述流管道中的流态数据的声学模型的数据库；以及
(b)处理器，其执行以下步骤：
(1)将声发射信号分割成数字声发射段序列；
(2)确定所述数字声发射段序列中的每一个数字声发射段的特征向量；以及
(3)处理所述特征向量以确定所述多相流的流参数的模型。
25.根据权利要求24所述的数据处理系统，还包括存储实际多相流状况的数据库的所述数据存储器。
26.根据权利要求24所述的数据处理系统，还包括在处理所述特征向量以确定流参数的模型时接收来自所述数据库的实际多相流状况数据作为输入的所述处理器。
27.根据权利要求24所述的数据处理系统，其中，所述处理器在处理所述特征向量以确定流参数的模型时执行隐马尔可夫建模。
28.根据权利要求24所述的数据处理系统，其中，所述处理器在处理所述特征向量以确定流参数的模型时执行以下步骤：
基于所确定的所述数字声段序列的流向量和所存储的所述流管道中的流态数据的声学模型来确定所述多相流的流参数的模型。
29.根据权利要求24所述的数据处理系统，其中，所述处理器在确定流参数的模型时，基于所确定的所述数字声段序列的流向量和所存储的所述流管道中的流态数据的声学模型来确定所述多相流的流参数的最可能的模型。
30.根据权利要求24所述的数据处理系统，其中，所述处理器在处理所述特征向量以确定流参数的模型时执行人工神经网络建模。
31.根据权利要求24所述的数据处理系统，其中,所述处理器在处理所述特征向量以确定流参数的模型时执行以下步骤：
(a)接收所述特征向量作为人工神经网络处理的输入状态；
(b)基于所述输入状态执行所述人工神经网络处理，以确定所述多相流的模型的流参数；以及
(c)提供所确定的所述多相流的模型的流参数作为输出状态。
32.根据权利要求24所述的数据处理系统，其中，所述处理器在处理所述特征向量以确定流参数的模型时执行隐马尔可夫建模和人工神经网络建模。
33.根据权利要求24所述的数据处理系统，其中，所述处理器还形成用于执行处理所述特征向量的步骤的训练模型。
34.根据权利要求24所述的数据处理系统，其中，所述处理器还将所形成的训练模型存储在所述存储器中。
35.根据权利要求24所述的数据处理系统，其中，所述处理器还提供所确定的所述多相流的流参数的模型用以显示。
36.一种数据存储装置，其在非暂时性计算机可读介质中存储有计算机可操作指令，该计算机可操作指令用于使数据处理系统在所述数据处理系统的处理器中基于来自多相流的声发射和存储在计算机的数据库中的管道中的流态数据的声学模型来确定流管道中的流体的多相流的流参数，存储在所述数据存储装置中的所述指令使得所述数据处理系统中的处理器执行以下步骤：
(a)将声发射信号分割成数字声发射段序列；
(b)确定所述数字声发射段序列中的每一个数字声发射段的特征向量；
(c)处理所述特征向量以确定所述多相流的流参数的模型。
37.根据权利要求36所述的数据存储装置，其中，所述数据存储器存储实际多相流状况的数据库，并且其中，所述指令还包括使所述处理器执行基于来自所述数据库的实际多相流状况数据来处理所述特征向量以确定流参数的模型的步骤的指令。
38.根据权利要求36所述的数据存储装置，其中，用于处理所述特征向量以确定流参数的模型的指令包括用于执行隐马尔可夫建模的指令。
39.根据权利要求36所述的数据存储装置，其中，用于处理所述特征向量以确定流参数的模型的指令包括用于执行以下步骤的指令：
基于所确定的所述数字声段序列的流向量和所存储的所述流管道中的流态数据的声学模型来确定所述多相流的流参数的模型。
40.根据权利要求36所述的数据存储装置，其中，用于处理所述特征向量以确定流参数的模型的指令包括用以执行以下步骤的指令：
基于所确定的所述数字声段序列的流向量和所存储的所述流管道中的流态数据的声学模型来确定所述多相流的流参数的最可能的模型。
41.根据权利要求36所述的数据存储装置，其中，用于处理所述特征向量以确定流参数的模型的指令包括用于执行人工神经网络建模的指令。
42.根据权利要求36所述的数据存储装置，其中，用于处理所述特征向量以确定流参数的模型的指令包括用于执行以下步骤的指令：
(a)接收所述特征向量作为人工神经网络处理的输入状态；
(b)基于所述输入状态执行所述人工神经网络处理，以确定所述多相流模型的流参数；以及
(c)提供所确定的所述多相流模型的流参数作为输出状态。
43.根据权利要求36所述的数据存储装置，其中，用于处理所述特征向量以确定流参数的模型的指令包括用于执行隐马尔可夫建模和人工神经网络建模的指令。
44.根据权利要求36所述的数据存储装置，其中，所述指令还包括用于执行以下步骤的指令：形成用于执行处理所述特征向量的步骤的训练模型。
45.根据权利要求36所述的数据存储装置，其中，所述指令还包括用于执行将所形成的训练模型存储在所述计算机的存储器中的步骤的指令。
46.根据权利要求36所述的数据存储装置，其中，所述指令还包括用于执行提供所确定的所述多相流的流参数的模型用以显示的步骤的指令。
Applications Claiming Priority (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US201462018727P | 2014-06-30 | 2014-06-30 | |
| US62/018,727 | 2014-06-30 | ||
| US14/735,217 US20150377667A1 (en) | 2014-06-30 | 2015-06-10 | Virtual multiphase flow metering and sand detection |
| US14/735,217 | 2015-06-10 | ||
| PCT/US2015/038457 WO2016003985A1 (en) | 2014-06-30 | 2015-06-30 | Virtual multiphase flow metering and sand detection |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106537136A true CN106537136A (zh) | 2017-03-22 |
Family
ID=54930151
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201580036015.6A Pending CN106537136A (zh) | 2014-06-30 | 2015-06-30 | 虚拟多相流计量和砂检测 |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US20150377667A1 (zh) |
| EP (1) | EP3161424A1 (zh) |
| JP (1) | JP6335390B2 (zh) |
| CN (1) | CN106537136A (zh) |
| CA (1) | CA2952295A1 (zh) |
| WO (1) | WO2016003985A1 (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111344710A (zh) * | 2017-09-26 | 2020-06-26 | 沙特阿拉伯石油公司 | 使用基于机器学习的模型进行成本有效的热力学流体特性预测的方法 |
| CN111771042A (zh) * | 2017-10-11 | 2020-10-13 | 英国石油勘探运作有限公司 | 使用声学频域特征来检测事件 |
| CN113272518A (zh) * | 2018-11-29 | 2021-08-17 | Bp探索操作有限公司 | 识别流体流入位置和流体类型的das数据处理 |
Families Citing this family (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP3311124B1 (en) * | 2015-06-22 | 2019-09-11 | Saudi Arabian Oil Company | System and method to provide entropy based characterization of multiphase flow |
| US10983238B2 (en) | 2016-09-26 | 2021-04-20 | Halliburton Energy Services, Inc. | Wellbore sand detection using passive acoustic array |
| US10883491B2 (en) * | 2016-10-29 | 2021-01-05 | Kelvin Inc. | Plunger lift state estimation and optimization using acoustic data |
| US10704935B2 (en) | 2016-12-04 | 2020-07-07 | Buoy Labs, Inc. | Fluid flow detector with tethered drag block |
| USD866375S1 (en) | 2017-08-02 | 2019-11-12 | Buoy Labs, Inc. | Water flow monitoring device |
| GB201713028D0 (en) * | 2017-08-14 | 2017-09-27 | Unify Group (Emea) Ltd | Detecting fluid flow |
| GB201714069D0 (en) * | 2017-09-01 | 2017-10-18 | Exnics Ltd | Apparatus and method |
| US10444196B2 (en) | 2017-09-20 | 2019-10-15 | Fisher Controls International Llc | Bandwidth-selectable acoustic emission apparatus and methods for transmitting time-averaged signal data |
| US10373598B2 (en) * | 2017-09-20 | 2019-08-06 | Fisher Controls International Llc | Integrated acoustic emission transducer apparatus and methods |
| US11781895B2 (en) | 2018-02-23 | 2023-10-10 | Buoy Labs, Inc. | Fluid flow analysis and management |
| US11567024B2 (en) * | 2018-12-14 | 2023-01-31 | Tata Consultancy Services Limited | System and method for detection of concentration of micro and nano particles in a fluid environment |
| US11079260B2 (en) * | 2019-06-25 | 2021-08-03 | Saudi Arabian Oil Company | Pipe cross-sectional fiber optic multiphase apparatus |
| CN111310900A (zh) * | 2020-02-24 | 2020-06-19 | 姜景元 | 一种多相流检测方法、装置、设备、介质 |
| CN111474243B (zh) * | 2020-04-27 | 2023-11-07 | 上海化工研究院有限公司 | 一种淤浆环管反应器内浆液浓度的测量系统 |
| US12044562B2 (en) | 2020-06-26 | 2024-07-23 | Schlumberger Technology Corporation | Multiphase flowmeter and related methods |
| US11428557B2 (en) | 2020-08-31 | 2022-08-30 | Saudi Arabian Oil Company | Determining fluid properties |
| US11525723B2 (en) | 2020-08-31 | 2022-12-13 | Saudi Arabian Oil Company | Determining fluid properties |
| US11592984B2 (en) * | 2020-09-11 | 2023-02-28 | Seagate Technology Llc | Onboard machine learning for storage device |
| US11748251B2 (en) * | 2021-01-08 | 2023-09-05 | Microsoft Technology Licensing, Llc | Storing tensors in memory based on depth |
| US12085687B2 (en) | 2022-01-10 | 2024-09-10 | Saudi Arabian Oil Company | Model-constrained multi-phase virtual flow metering and forecasting with machine learning |
| CN114609410B (zh) * | 2022-03-25 | 2022-11-18 | 西南交通大学 | 一种基于声学信号的便携式风特性测量设备及智能算法 |
| WO2023234955A1 (en) * | 2022-06-03 | 2023-12-07 | Halliburton Energy Services, Inc. | Memory tool for a retrievable flow meter device |
| US20240118118A1 (en) * | 2022-10-06 | 2024-04-11 | Halliburton Energy Services, Inc. | Virtual flow metering using acoustics |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020029883A1 (en) * | 2000-01-24 | 2002-03-14 | Vinegar Harold J. | System and method for fluid flow optimization |
| CN1912612A (zh) * | 2006-08-15 | 2007-02-14 | 天津工业大学 | 一种多相流物质的检测方法及检测装置 |
| KR20090000858A (ko) * | 2007-06-28 | 2009-01-08 | 주식회사 케이티 | 멀티모달 기반의 정보 검색 장치 및 방법 |
| EP2199755A1 (en) * | 2008-12-16 | 2010-06-23 | Nederlandse Organisatie voor toegepast-natuurwetenschappelijk Onderzoek TNO | An apparatus, a method and a computer program for recognition of flow regimes in a multiphase fluid flowing in a conduit |
| CN103733061A (zh) * | 2010-09-03 | 2014-04-16 | 洛斯阿拉莫斯国家安全股份有限公司 | 多相流体特性系统 |
| CN103852525A (zh) * | 2012-11-29 | 2014-06-11 | 沈阳工业大学 | 基于ar-hmm的声发射信号识别方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5741980A (en) * | 1994-11-02 | 1998-04-21 | Foster-Miller, Inc. | Flow analysis system and method |
| US5600073A (en) * | 1994-11-02 | 1997-02-04 | Foster-Miller, Inc. | Method and system for analyzing a two phase flow |
| BR0107821B1 (pt) * | 2000-01-24 | 2010-09-08 | método para controlar o escoamento de fluido em um conduto, e, poço de óleo de elevação de gás. | |
| AU3334800A (en) * | 2000-03-09 | 2001-09-17 | Vladimir Drobkov | Simultaneous determination of multiphase flowrates and concentrations |
| FR2812389B1 (fr) * | 2000-07-27 | 2002-09-13 | Inst Francais Du Petrole | Methode et systeme pour estimer en temps reel le mode d'ecoulement d'une veine fluide polyphasique, en tous points d'une conduite |
| JP5208124B2 (ja) * | 2007-11-15 | 2013-06-12 | 国立大学法人北海道大学 | 超音波混相流量計、超音波混相流量計測プログラム、および超音波を用いた混相流量計測方法 |
-
2015
- 2015-06-10 US US14/735,217 patent/US20150377667A1/en not_active Abandoned
- 2015-06-30 CA CA2952295A patent/CA2952295A1/en not_active Abandoned
- 2015-06-30 JP JP2017521020A patent/JP6335390B2/ja not_active Expired - Fee Related
- 2015-06-30 CN CN201580036015.6A patent/CN106537136A/zh active Pending
- 2015-06-30 EP EP15738544.4A patent/EP3161424A1/en not_active Withdrawn
- 2015-06-30 WO PCT/US2015/038457 patent/WO2016003985A1/en active Application Filing
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020029883A1 (en) * | 2000-01-24 | 2002-03-14 | Vinegar Harold J. | System and method for fluid flow optimization |
| CN1912612A (zh) * | 2006-08-15 | 2007-02-14 | 天津工业大学 | 一种多相流物质的检测方法及检测装置 |
| KR20090000858A (ko) * | 2007-06-28 | 2009-01-08 | 주식회사 케이티 | 멀티모달 기반의 정보 검색 장치 및 방법 |
| EP2199755A1 (en) * | 2008-12-16 | 2010-06-23 | Nederlandse Organisatie voor toegepast-natuurwetenschappelijk Onderzoek TNO | An apparatus, a method and a computer program for recognition of flow regimes in a multiphase fluid flowing in a conduit |
| CN103733061A (zh) * | 2010-09-03 | 2014-04-16 | 洛斯阿拉莫斯国家安全股份有限公司 | 多相流体特性系统 |
| CN103852525A (zh) * | 2012-11-29 | 2014-06-11 | 沈阳工业大学 | 基于ar-hmm的声发射信号识别方法 |
Non-Patent Citations (1)
| Title |
|---|
| ALI MAHVASH ET AL.: "Two-phase flow pattern identification using continuous hidden Markov model", 《INTERNATIONAL JOURNAL OF MULTIPHASE FLOW》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111344710A (zh) * | 2017-09-26 | 2020-06-26 | 沙特阿拉伯石油公司 | 使用基于机器学习的模型进行成本有效的热力学流体特性预测的方法 |
| CN111771042A (zh) * | 2017-10-11 | 2020-10-13 | 英国石油勘探运作有限公司 | 使用声学频域特征来检测事件 |
| CN113272518A (zh) * | 2018-11-29 | 2021-08-17 | Bp探索操作有限公司 | 识别流体流入位置和流体类型的das数据处理 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2017527823A (ja) | 2017-09-21 |
| WO2016003985A1 (en) | 2016-01-07 |
| JP6335390B2 (ja) | 2018-05-30 |
| EP3161424A1 (en) | 2017-05-03 |
| CA2952295A1 (en) | 2016-01-07 |
| US20150377667A1 (en) | 2015-12-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106537136A (zh) | 虚拟多相流计量和砂检测 | |
| US10472957B2 (en) | Entropy based multiphase flow detection | |
| US10380483B2 (en) | Method and apparatus for training language model, and method and apparatus for recognizing language | |
| Heittola et al. | The machine learning approach for analysis of sound scenes and events | |
| Lindgren et al. | Speech recognition using reconstructed phase space features | |
| CN111811617B (zh) | 一种基于短时傅里叶变换和卷积神经网络的液位预测方法 | |
| CN109192221A (zh) | 一种基于聚类的使用语音判断帕金森严重程度检测方法 | |
| US20100169085A1 (en) | Model based real time pitch tracking system and singer evaluation method | |
| CN115273904A (zh) | 一种基于多特征融合的愤怒情绪识别方法及装置 | |
| Mercado et al. | Classification of humpback whale vocalizations using a self-organizing neural network | |
| Balamurali et al. | Acoustic prediction of flowrate: varying liquid jet stream onto a free surface | |
| Perez-Carrillo | Statistical models for the indirect acquisition of violin bowing controls from audio analysis | |
| Ekiz et al. | Covid-19 Detection from Cough, Breath, And Speech Sounds with Short-Time Fourier Transform and a CNN Model | |
| Wiśniewski et al. | Improved approach to automatic detection of speech disorders based on the Hidden Markov Models approach | |
| Emmanuel et al. | Machine Learning Based Rain Gauge Using Acoustic Data | |
| Wache Ngateu | Acoustic detection of the short pulse call of Brydes whales using time domain features and hidden Marcov models | |
| Ham et al. | Discrimination of volcano activity and mountain-associated waves using infrasonic data and a backpropagation neural network | |
| Aslim et al. | Acoustic prediction of flowrate: varying liquid jet stream onto a free surface | |
| HONG | 1-D AND 2-D CONVOLUTION NEURAL NETWORK FOR BIRD SOUND DETECTION | |
| May | Nonlinear dynamic invariants for continuous speech recognition | |
| Ti et al. | On-site wind speed recovery from smartphone audio: Time domain deep learning approach, laboratory validation and outdoor field test | |
| Bastas | Nocturnal bird call recognition system for wind farm applications | |
| Curless | Classification of wolf call types using remote sensor technology | |
| Dheeraj | Evaluation Of Unsupervised Models Using Minimal Pair ABX Measure | |
| Stammers | Audio Event Classification for Urban Soundscape Analysis |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170322 |
|
| WD01 | Invention patent application deemed withdrawn after publication |