RU2726275C1 - Method of determining liquid flow rate and vortex flow meter for its implementation - Google Patents

Method of determining liquid flow rate and vortex flow meter for its implementation Download PDF

Info

Publication number
RU2726275C1
RU2726275C1 RU2019145775A RU2019145775A RU2726275C1 RU 2726275 C1 RU2726275 C1 RU 2726275C1 RU 2019145775 A RU2019145775 A RU 2019145775A RU 2019145775 A RU2019145775 A RU 2019145775A RU 2726275 C1 RU2726275 C1 RU 2726275C1
Authority
RU
Russia
Prior art keywords
signal
spectrum
flow rate
vibration
vortex
Prior art date
Application number
RU2019145775A
Other languages
Russian (ru)
Other versions
RU2726275C9 (en
Inventor
Рустам Саинтурович Вахрушев
Павел Владимирович Изотов
Алексей Владимирович Кузнецов
Original Assignee
Акционерное общество "Ижевский радиозавод"
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Акционерное общество "Ижевский радиозавод" filed Critical Акционерное общество "Ижевский радиозавод"
Priority to RU2019145775A priority Critical patent/RU2726275C9/en
Application granted granted Critical
Publication of RU2726275C1 publication Critical patent/RU2726275C1/en
Publication of RU2726275C9 publication Critical patent/RU2726275C9/en

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01FMEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/05Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects
    • G01F1/20Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects by detection of dynamic effects of the flow
    • G01F1/32Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects by detection of dynamic effects of the flow using swirl flowmeters

Landscapes

  • Measuring Volume Flow (AREA)
  • Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • General Physics & Mathematics (AREA)

Abstract

FIELD: measuring equipment.SUBSTANCE: invention relates to measurement equipment, namely to methods and facilities intended for measurement of liquid flow rate, and can be used in various industries for purposes of control, regulation and accounting of medium flows. In particular, the invention can be used for determination of liquid flow rate in a well by a vortex flow meter under conditions of action of low-frequency vibration loads. Method of determining liquid flow rate consists in continuously receiving a signal which includes a useful frequency component which characterizes the velocity of the liquid flow and a noise frequency component which characterizes the intensity of vibration. PT signal obtained at the current moment of time is digitized and placed into the spectrum. Vibration sensor measures vibration, digitizes the measured VS signal and spreads it into a spectrum. Further, the useful signal spectrum is picked up by subtracting the signal spectrum VS from the spectrum of the PT signal. Derived spectrum of the useful signal is used to calculate the vortex generation frequency from which the fluid flow rate is determined, which determines the flow rate of the liquid at the current moment in order to transmit information on the determined liquid flow to the consumer. Vortex flow meter configured to supply voltage to power-requiring elements and comprising series-connected primary transducer (PT), processing device (hereinafter, PD), vibration sensor (further, VS). PD is configured to decompose in the spectra of the digitized PT signal and the VS signal, extracting the useful signal spectrum from the PT signal spectrum, application of a useful adaptive filtration signal in the form of a Hamming Window to the said selected spectrum and evaluation of weight functions, calculation of the vortex generation frequency, calculation of the liquid flow rate value at the current moment in time, automatic adjustment of the PT signal level, filtration of its low frequency.EFFECT: high accuracy of determining the frequency of vortex formation, especially during strong vibration action due to: higher vibration stability of the vortex flow meter, possibility of excluding false measurement of flow rate and increasing the minimum measured level of flow rate in conditions of influencing the flow meter of dynamically varying vibration loads.12 cl, 6 dwg

Description

Изобретение относится к измерительной технике, а именно к способам и средствам, предназначенным для измерения расхода жидкостей, и может быть использовано в различных отраслях народного хозяйства для целей контроля, регулирования и учета потоков среды. В частности, изобретение может быть использовано для определения расхода жидкости в скважине вихревым расходомером в условиях воздействия на него низкочастотных вибрационных нагрузок.The invention relates to measuring equipment, and in particular to methods and means intended for measuring the flow of liquids, and can be used in various sectors of the economy for the purpose of monitoring, regulating and accounting for environmental flows. In particular, the invention can be used to determine the flow rate of a fluid in a well by a vortex flow meter under conditions of exposure to low-frequency vibration loads.

Известно, что для определения расхода жидкости в скважине вихревой расходомер устанавливают в ней в непосредственной близости от погружного электродвигателя (ПЭД), например, на выкиде насоса. Частота оборотов работы ПЭД, как правило, находится в диапазоне от 10 до 80 Гц. При этом на данных и кратных им частотах возникает вибрация. Эта вибрационная нагрузка передается по колонне НКТ близлежащему оборудованию, в том числе вихревому расходомеру. В случаях, когда частота и амплитуда вибрации приближаются к частоте вихреобразования, в частности, при низкочастотных вибрационных нагрузках, частота вибрации и частота полезного сигнала сливаются друг с другом и перестают различаться. Вследствие этого частота вибрации может ошибочно восприниматься устройством обработки расходомера как частота полезного сигнала, характеризующая расход жидкости. В результате это приводит к низкой точности и достоверности определения расхода жидкости.It is known that to determine the flow rate of a fluid in a well, a vortex flowmeter is installed in it in the immediate vicinity of a submersible electric motor (SEM), for example, on a pump outlet. The speed of operation of the PEM, as a rule, is in the range from 10 to 80 Hz. At the same time, vibration occurs at the data and at its multiple frequencies. This vibration load is transmitted through the tubing string to nearby equipment, including a vortex flowmeter. In cases when the frequency and amplitude of vibration approach the frequency of vortex formation, in particular, at low-frequency vibration loads, the vibration frequency and the frequency of the useful signal merge with each other and cease to differ. As a result, the vibration frequency may be mistakenly perceived by the flowmeter processing device as the frequency of the useful signal characterizing the flow rate of the liquid. As a result, this leads to low accuracy and reliability of determining fluid flow.

Наиболее близким является способ определения расхода жидкости, реализованный в универсальном вихревом расходомере (патент РФ на изобретение №2192623, МПК G01F 1/32, 2000 г.), содержащем первичный преобразователь и электрическую схему обработки сигнала первичного преобразователя, которая содержит последовательно включенные первый нормирующий усилитель, фильтр, второй нормирующий усилитель и аналого-цифровой преобразователь, микроконтроллер, управляющие выходы которого подключены соответственно к управляющим входам первого и второго нормирующих усилителей, интерфейсы расходомера. Функцию адаптивно перестраиваемых фильтров выполняет алгоритм цифровой адаптивной фильтрации, реализованный в микроконтроллере. Способ осуществляется следующим образом. Получают сигнал с первичного преобразователя, производит выборку сигнала с частотой опроса, необходимой для восстановления сигнала в диапазоне измеряемых частот, оценивают уровень оцифрованного входного сигнала и выполняют автоматическую регулировку усиления, производит поиск и выделение сигнала вихревой частоты в соответствии с алгоритмом цифровой адаптивной фильтрации, вычисляют текущее значение вихревой частоты, вычисляют значение расхода. Способ позволяет определить скорость или расход жидкости жидких и газообразных сред в скважине при воздействии возмущений в виде вибрации трубопровода, пульсаций потока и акустических шумов.The closest is a method for determining fluid flow, implemented in a universal vortex flowmeter (RF patent for the invention No. 2192623, IPC G01F 1/32, 2000), containing a primary transducer and an electrical signal processing circuit of the primary transducer, which contains a first normalizing amplifier in series , a filter, a second normalizing amplifier and an analog-to-digital converter, a microcontroller, the control outputs of which are connected respectively to the control inputs of the first and second normalizing amplifiers, flowmeter interfaces. The function of adaptively tunable filters is performed by the digital adaptive filtering algorithm implemented in the microcontroller. The method is as follows. They receive a signal from the primary converter, select the signal with the polling frequency necessary to restore the signal in the range of measured frequencies, evaluate the level of the digitized input signal and perform automatic gain control, search and extract the vortex frequency signal in accordance with the digital adaptive filtering algorithm, calculate the current vortex frequency value, the flow value is calculated. The method allows to determine the speed or flow rate of liquid of liquid and gaseous media in the well when exposed to disturbances in the form of pipeline vibration, flow pulsations and acoustic noise.

Недостатками способа является низкая точность и достоверность определения расхода жидкости, в том числе при низкочастотных вибрационных нагрузках, поскольку способ не позволяет выделить вибрацию при смешении вибрации и полезного сигнала, характеризующего расход жидкости и исключить ее при определении расхода жидкости.The disadvantages of the method are the low accuracy and reliability of determining the flow rate of the liquid, including at low frequency vibration loads, since the method does not allow vibration to be distinguished by mixing the vibration and the useful signal characterizing the flow rate of the liquid and to exclude it when determining the flow rate of the liquid.

Наиболее близким является универсальный вихревой расходомер (патент РФ на изобретение №2192623, МПК G01F 1/32, 2000 г.), содержащий первичный преобразователь и электрическую схему обработки сигнала первичного преобразователя, которая содержит последовательно включенные первый нормирующий усилитель, фильтр, второй нормирующий усилитель и аналого-цифровой преобразователь, микроконтроллер, управляющие выходы которого подключены соответственно к управляющими входами первого и второго нормирующих усилителей, интерфейсы расходомера. Функцию адаптивно перестраиваемых фильтров выполняет алгоритм цифровой адаптивной фильтрации, реализованный в микроконтроллере. Устройство позволяет определить расход жидкости.The closest is a universal vortex flowmeter (RF patent for the invention No. 2192623, IPC G01F 1/32, 2000), containing a primary transducer and an electrical signal processing circuit of the primary transducer, which contains a first normalizing amplifier, a filter, a second normalizing amplifier and an analog-to-digital converter, a microcontroller, the control outputs of which are connected respectively to the control inputs of the first and second normalizing amplifiers, flowmeter interfaces. The function of adaptively tunable filters is performed by the digital adaptive filtering algorithm implemented in the microcontroller. The device allows you to determine the flow rate.

Недостатком является низкая точность определения расхода жидкости, поскольку отсутствует возможность измерения вибрации и, ее учета при определении расхода жидкости, в том числе при низкочастотных вибрационных нагрузках. Кроме того, при уменьшении скорости потока жидкости, происходит ослабление полезного сигнала расхода, при этом полезный сигнал становится соизмерим с нижними гармониками. Это ограничивает нижний предел измерения потока расходомера.The disadvantage is the low accuracy of determining the flow rate of the liquid, since there is no possibility of measuring vibration and, taking it into account when determining the flow rate of the liquid, including low-frequency vibration loads. In addition, with a decrease in the fluid flow rate, the useful flow signal is weakened, while the useful signal becomes comparable with the lower harmonics. This limits the lower limit of flow meter flow measurement.

Целью изобретения является устранение указанных недостатков.The aim of the invention is to remedy these disadvantages.

Технический результат заключается в повышении точности определения расхода жидкости за счет применения способа определения расхода жидкости в скважине, в том числе при низкочастотных вибрационных нагрузках.The technical result consists in increasing the accuracy of determining the fluid flow rate through the application of the method for determining the fluid flow rate in the well, including at low frequency vibration loads.

Кроме того, технический результат также заключается в повышении точности определения расхода жидкости за счет обеспечения возможности измерить (выделить) вибрацию и возможности ее исключить при определении расхода жидкости, а также в расширении функциональных возможностей, заключающемся в расширении рабочего диапазона вихревого расходомера в области малых расходов.In addition, the technical result also consists in increasing the accuracy of determining the liquid flow rate by providing the ability to measure (isolate) the vibration and the possibility of eliminating it when determining the liquid flow rate, as well as in expanding the functionality, which consists in expanding the working range of the vortex flowmeter in the field of low flow rates.

Технический результат достигается в способе определения расхода жидкости, заключающемся в том, что первичным преобразователем непрерывно принимают сигнал, включающий полезную частотную составляющую, характеризующую скорость потока жидкости и шумовую частотную составляющую, характеризующую интенсивность вибрации. Полученный в текущий момент времени сигнал ПП оцифровывают и раскладывают в спектр. Посредством датчика вибрации измеряют вибрацию, оцифровывают измеренный сигнал ДВ, раскладывают в спектр. Далее выделяют спектр полезного сигнала путем вычитания из спектра сигнала ПП спектр сигнала ДВ. По выделенному спектру полезного сигнала вычисляют значение частоты вихреобразования, по которой вычисляют значение скорости потока жидкости, определяющее расход жидкости в текущий момент времени для передачи информации об определенном расходе жидкости потребителю.The technical result is achieved in the method for determining the flow rate of the liquid, which consists in the fact that the primary converter continuously receives a signal including a useful frequency component characterizing the speed of the fluid flow and a noise frequency component characterizing the intensity of the vibration. The PP signal received at the current time is digitized and laid out in the spectrum. Using a vibration sensor, vibration is measured, the measured DW signal is digitized, and laid out in a spectrum. Next, the spectrum of the useful signal is isolated by subtracting the spectrum of the DW signal from the spectrum of the PP signal. From the selected spectrum of the useful signal, the value of the vortex formation frequency is calculated, from which the value of the fluid flow rate is calculated, which determines the fluid flow rate at the current time to transmit information about the specific fluid flow rate to the consumer.

Предпочтительно, чтобы измерение и оцифровка вибрации датчиком вибрации ДВ осуществлялись одновременно с оцифровкой сигнала ПП.It is preferable that the measurement and digitization of vibration by the vibration sensor of the LW are carried out simultaneously with the digitization of the PP signal.

Предпочтительно, чтобы после выделения спектра полезного сигнала к нему применяли адаптивную фильтрацию, после чего вычисляли бы значение частоты вихреобразования.It is preferable that, after the spectrum of the useful signal has been allocated, adaptive filtering is applied to it, after which the value of the vortex formation frequency would be calculated.

Предпочтительно, чтобы адаптивную фильтрацию осуществляли посредством последовательного применения к выделенному спектру полезного сигнала Окна Хемминга и оценки весовых функций.Adaptive filtering is preferably carried out by sequentially applying the Hamming Window to the selected spectrum of the useful signal and evaluating the weight functions.

Предпочтительно, чтобы перед оцифровкой уровень амплитуды полученного сигнала ПП автоматически регулировали.Preferably, before digitizing, the amplitude level of the received PP signal is automatically adjusted.

Предпочтительно, чтобы после автоматической регулировки уровня сигнала ПП фильтровали его нижние частоты.It is preferable that after automatic adjustment of the signal level the PP filtered its low frequencies.

Предпочтительно, чтобы вычисленное значение скорости потока жидкости корректировали посредством заранее заданных калибровочных коэффициентов.Preferably, the calculated value of the fluid flow rate is corrected by means of predetermined calibration factors.

Кроме того, указанная цель достигается в вихревом расходомере, выполненном с возможностью подачи питающего напряжения на требующие питания элементы, и содержащем соединенные последовательно первичный преобразователь (далее, ПП), устройство обработки (далее, УО), датчик вибрации (далее, ДВ). При этом УО выполнен с возможностью разложения в спектры оцифрованных сигнала ПП и сигнала ДВ, выделения спектра полезного сигнала из спектра сигнала ПП, применения к упомянутому выделенному спектру полезного сигнала адаптивной фильтрации в виде Окна Хемминга и оценки весовых функций, вычисления значения частоты вихреобразования, вычисления значения скорости потока жидкости в текущий момент времени, автоматической регулировки уровня сигнала ПП, фильтрации его нижних частот.In addition, this goal is achieved in a vortex flowmeter, configured to supply power to power-demanding elements, and containing a series-connected primary Converter (hereinafter, PP), a processing device (hereinafter, UO), a vibration sensor (hereinafter, DV). In this case, the RO is capable of decomposing into the spectra of the digitized PP signal and the DV signal, extracting the spectrum of the useful signal from the spectrum of the PP signal, applying adaptive filtering in the form of a Hamming window and estimating weight functions to the selected spectrum of the useful signal, computing the value of the vortex frequency, calculating the value fluid flow rate at the current time, automatic adjustment of the PP signal level, filtering its low frequencies.

Предпочтительно, чтобы УО включал контроллер, с которым соединены физический интерфейс, оперативное запоминающее устройство (далее, ОЗУ), с которым, в свою очередь, последовательно соединены первый аналогово-цифровой преобразователь (первый АЦП), соединенный с ПП, второй аналогово-цифровой преобразователь (второй АЦП), соединенный с ДВ.Preferably, the UO includes a controller with which a physical interface is connected, a random access memory (hereinafter RAM), which, in turn, is connected in series with a first analog-to-digital converter (first ADC) connected to the PC, a second analog-to-digital converter (second ADC) connected to the DV.

Предпочтительно, чтобы УО дополнительно включал узел автоматической регулировки уровня сигнала первичного преобразователя (далее, УАРУС ПП), содержащий последовательно соединенные друг с другом усилительный узел (далее, УУ), регулирующий усилитель (далее, РУ), при этом УУ подключен к ПП, РУ соединен с первым АЦП, причем контроллер УО через обратную связь связан с РУ с возможностью регулирования уровня амплитуды сигнала ПП.It is preferable that the UO additionally include a node for automatically adjusting the signal level of the primary converter (hereinafter, UARUS PP), containing an amplifier unit (hereinafter referred to as UU), a regulating amplifier (hereinafter RU) connected in series, while the UU is connected to the PP, RU connected to the first ADC, and the controller UO through feedback is connected to the RU with the ability to control the amplitude level of the signal of the PP.

Предпочтительно, чтобы УУ содержал усилитель заряда.Preferably, the SU contains a charge amplifier.

Предпочтительно, чтобы УУ содержал, по крайней мере, два усилителя заряда, соединенные параллельно друг с другом, подключенные последовательно с РУ.Preferably, the SU contains at least two charge amplifiers connected in parallel with each other, connected in series with the switchgear.

Предпочтительно, чтобы УАРУС ПП дополнительно содержал фильтр низких частот (далее, ФНЧ), включенный между РУ и первым АЦП.Preferably, UARUS PP additionally contains a low-pass filter (hereinafter, the low-pass filter), included between the switchgear and the first ADC.

Далее на фиг. 1 структурная схема предлагаемого вихревого расходомера, на фиг. 2 приведен пример предлагаемого вихревого расходомера, который реализует предлагаемый способ определения расхода жидкости в предпочтительном варианте, на фиг. 3 показан спектр сигнала ПП, включающий гармонику вибрации и гармонику вихреобразования, на фиг. 4 показан спектр сигнала ДВ, на фиг. 5 показан спектр полезного сигнала, фиг. 6 показан спектр выделенного полезного сигнала после адаптивной фильтрации.Next, in FIG. 1 is a structural diagram of the proposed vortex flowmeter; FIG. 2 shows an example of the proposed vortex flowmeter, which implements the proposed method for determining fluid flow in the preferred embodiment, FIG. 3 shows the spectrum of the PP signal, including the harmonic of vibration and the harmonic of vortex formation, in FIG. 4 shows the spectrum of the DW signal; FIG. 5 shows a spectrum of a useful signal; FIG. Figure 6 shows the spectrum of the extracted useful signal after adaptive filtering.

Вихревой расходомер на фиг. 1 включает последовательно соединенные ПП 1, УО 2, ДВ 3, при этом чувствительный элемент ПП 1 представляет собой пьезоэлемент, улавливающий вихри, образованные за телом обтекания. УО 2 (фиг. 2) содержит УАРУС 4 ПП (показан пунктирной линией), включающий последовательно соединенные друг с другом УУ 5, РУ 6, контроллер 7, физический интерфейс 8, ОЗУ 9, к которому последовательно подключены первый АЦП 10, второй АЦП 11. РУ 6 последовательно соединен первым АЦП 10. Для повышения точности, уменьшения помех между РУ 6 и первым АЦП 10 может быть дополнительно подключен фильтр низких частот (ФНЧ)12. Контроллер 7 охвачен обратной связью с РУ 6, обеспечивающей возможность автоматической регулировки коэффициентов усиления для согласования диапазонов выходного напряжения ПП 1 и входного напряжения первого АЦП 10. ДВ 3 соединен через второй АЦП 11 с контроллером 7.The vortex flowmeter of FIG. 1 includes series-connected PP 1, UO 2, DV 3, while the sensitive element of PP 1 is a piezoelectric element that traps the vortices formed behind the flow body. UO 2 (Fig. 2) contains UARUS 4 PP (shown by a dashed line), including UU 5, RU 6, controller 7, physical interface 8, RAM 9, to which the first ADC 10, the second ADC 11 are connected in series The RU 6 is connected in series with the first ADC 10. In order to increase accuracy and reduce interference between the RU 6 and the first ADC 10, a low-pass filter (LPF) 12 can be additionally connected. The controller 7 is covered by feedback from the RU 6, providing the ability to automatically adjust the gain to match the ranges of the output voltage of the PP 1 and the input voltage of the first ADC 10. DV 3 is connected through the second ADC 11 to the controller 7.

На фиг. 3 показан спектр сигнала ПП, в котором гармоники смешаны, а также в отдельности гармоника вибрации (вибрационная шумовая составляющая) и гармоника вихреобразования (составляющая полезного сигнала). На фигуре 4 показан спектр сигнала ДВ, в котором выделена максимальная гармоника, характеризующая интенсивность вибрации. На фигуре 5 показан спектр полезного сигнала. В частности, спектр полезного сигнала может быть получен программным исключением вибрационной гармоники из спектра сигнала ПП. Для более точного выделения гармоники, отвечающей за скорость потока жидкости, а также для улучшения соотношения сигнал-шум к выделенному спектру полезного сигнала применяют адаптивную фильтрация в виде последующих друг за другом окна Хемминга и оценки весовых функций. Результат применения адаптивной фильтрации представлен на фигуре 6 и отражает увеличенное соотношение сигнал-шум, благодаря чему появляется возможность более точно получить основную частоту вихреобразования (в спектре полезного сигнала), по которой определяют скорость потока (расход) жидкости.In FIG. Figure 3 shows the spectrum of the PP signal in which harmonics are mixed, as well as individually the harmonic of vibration (vibrational noise component) and harmonic of vortex formation (component of the useful signal). Figure 4 shows the spectrum of the DW signal, in which the maximum harmonic characterizing the intensity of vibration is highlighted. The figure 5 shows the spectrum of the useful signal. In particular, the spectrum of the useful signal can be obtained by software excluding vibrational harmonics from the spectrum of the PP signal. For a more accurate selection of the harmonics responsible for the fluid flow rate, as well as for improving the signal-to-noise ratio to the selected spectrum of the useful signal, adaptive filtering is applied in the form of Hamming windows that follow one after another and weight functions are estimated. The result of the use of adaptive filtering is shown in FIG. 6 and reflects an increased signal to noise ratio, which makes it possible to more accurately obtain the main vortex frequency (in the spectrum of the useful signal), which determines the flow rate (flow rate) of the liquid.

В качестве ПП 1 может быть использован, например, датчик изгибающего момента пьезоэлектрический 108м-23-10мм ТУ4213-108-24172160-08, располагаемый за телом обтекания в проточной части вихревого расходомера, установленного в скважине в разрыв труб НКТ. Чувствительный элемент ПП 1 улавливает вихри, образовывающиеся за телом обтекания и преобразует их в аналоговый сигнал, характеризующий интенсивность вихреобразования и, как следствие, скорость потока жидкости. Этот сигнал является комплексным и включает в себя полезную частотную составляющую, характеризующую скорость потока жидкости и шумовую частотную составляющую, характеризующую интенсивность вибрации в текущий момент времени. ПП 1 соединен УО 2.As PP 1, for example, a piezoelectric bending moment sensor 108m-23-10mm TU4213-108-24172160-08 located behind the flow body in the flow part of the vortex flowmeter installed in the borehole into the tubing rupture can be used. The sensitive element PP 1 captures the vortices formed behind the body of the flow and converts them into an analog signal characterizing the intensity of the vortex formation and, as a consequence, the fluid flow rate. This signal is complex and includes a useful frequency component that characterizes the fluid flow rate and a noise frequency component that characterizes the intensity of the vibration at the current time. PP 1 is connected to UO 2.

УО 2 осуществляет следующие основные функций:UO 2 carries out the following main functions:

- оцифровку сигнала ПП 1 и сигнала ДВ 3 и их сохранение в ОЗУ 9;- digitization of the signal PP 1 and the signal DV 3 and their storage in RAM 9;

- автоматическую регулировку уровня амплитуды сигнала ПП 1 под соответствующий диапазон первого АЦП 10;- automatic adjustment of the amplitude level of the signal PP 1 for the corresponding range of the first ADC 10;

- разложение сигнала ПП 1 и сигнала ДВ 3 в спектры;- decomposition of the signal PP 1 and the signal DV 3 into spectra;

- выделение спектра полезного сигнала (очищен от вибраций) (путем вычитания из спектра сигнала ПП 1 спектр сигнала ДВ 3;- selection of the spectrum of the useful signal (cleaned from vibrations) (by subtracting from the spectrum of the signal PP 1 spectrum signal 3;

- применение адаптивной фильтрации к спектру выделенного полезного сигнала посредством последовательного применения Окна Хемминга и оценки весовых функций (позволяет придать полезной гармонической составляющей вихреобразования больший вес в результирующем значении по сравнению с шумовой составляющей);- the use of adaptive filtering to the spectrum of the selected useful signal by successive application of the Hamming Window and the estimation of weight functions (allows you to give the useful harmonic component of the vortex formation more weight in the resulting value compared to the noise component);

- вычисления значения частоты вихреобразования,- calculating the value of the vortex frequency,

- вычисления значения скорости потока жидкости;- calculating the value of the fluid flow rate;

- вычисления значения частоты вибрации и пр.- calculating the value of the vibration frequency, etc.

В данной реализации УО 2 включает УАРУС 4 ПП, контроллер 7 (в частности, микроконтроллер), с которым соединены физический интерфейс 8, ОЗУ 9, с которым в свою очередь последовательно соединены первый АЦП 10, предназначенный для преобразования в цифровой вид сигнала ПП 1 и второй АЦП 11, предназначенный для преобразования в цифровой вид измеренного сигнала ДВ 3.In this implementation, UO 2 includes UARUS 4 PP, a controller 7 (in particular, a microcontroller), which is connected to a physical interface 8, RAM 9, which in turn is connected in series with the first ADC 10, intended for digitalization of the signal PP 1 and the second ADC 11, designed to convert into digital form the measured signal DV 3.

УАРУС 4 ПП предназначен для предварительной обработки сигнала ПП 1 под соответствие сигнала ПП 1 диапазону первого АЦП 10 и включает последовательно соединенные УУ 5, РУ 6. Если принятый сигнал ПП 1 сразу удовлетворяет диапазону первого АЦП 10, сигнал ПП 1 может быть подан сразу на первый АЦП 10, минуя УАРУС 4 ПП, как показано на фиг. 2.UARUS 4 PP is intended for preprocessing the signal of PP 1 to match the signal of PP 1 to the range of the first ADC 10 and includes series-connected UU 5, RU 6. If the received signal PP 1 immediately satisfies the range of the first ADC 10, the signal of PP 1 can be applied immediately to the first ADC 10, bypassing UARUS 4 PP, as shown in FIG. 2.

УУ 5 предназначен для предварительного усиления сигнала ПП 1 и включает, по крайней мере, один усилитель заряда, в частности, интегратор тока (построенный на основе метода интегрирования тока и выполненный на операционном усилителе), в частности, выполненный на микросхеме TL082, фирма Texas Instruments. В частности, для фильтрации симфазных помех, образующихся в результате работы электронных устройств вихревого расходомера, может быть использован двухканальный ПП 1 (на выходе ПП 1 имеем прямой и инверсный сигналы). В этом случае УУ 5 может включать первый усилитель заряда, обеспечивающий усиление прямого сигнала ПП 1, параллельно которому подключают второй усилитель заряда, обеспечивающий усиление инверсного сигнала ПП 1. С выхода УУ 5 усиленные прямой сигнал и инверсный сигнал поступают на вход РУ 6, где за счет их разницы на его выходе формируется один (общий) усиленный сигнал ПП 1, поступающий через первый АЦП 10 в ОЗУ 9, где он сохраняется.UU 5 is designed to pre-amplify the signal of PP 1 and includes at least one charge amplifier, in particular, a current integrator (built on the basis of the current integration method and executed on an operational amplifier), in particular, made on a TL082 chip, Texas Instruments . In particular, for filtering the common-mode interference resulting from the operation of electronic devices of the vortex flowmeter, two-channel PP 1 can be used (at the output of PP 1 we have direct and inverse signals). In this case, SU 5 can include a first charge amplifier that provides amplification of the direct signal of PP 1, in parallel with which a second charge amplifier is connected that provides amplification of the inverse signal of PP 1. From the output of SU 5, the amplified direct signal and inverse signal are fed to input RU 6, where due to their difference, at its output, one (common) amplified signal of PP 1 is generated, which comes through the first ADC 10 to RAM 9, where it is stored.

В качестве РУ 6 может быть использован, например, широкоизвестный измерительный усилитель (инструментальный), выполненный на трех операционных усилителях, в частности, выполненный на микросхеме TL082, производитель Texas Instruments. РУ 6 соединен через первый АЦП 10 с контроллером 7, охваченным с ним обратной связью с возможностью автоматического регулирования уровня амплитуды сигнала ПП 1 посредством коэффициентов усиления в случае несоответствия уровня амплитуды оцифрованного сигнала ПП 1 заданным значениям диапазона уровня амплитуды сигнала (хранящимся в ОЗУ 9). Таким образом, обратная связь позволяет автоматически адаптировать уровень амплитуды сигнала ПП 1 под требуемый аппаратной частью диапазон значений амплитуд с максимально возможной дискретностью, что дает возможность получить указанный сигнал в цифровом виде с высокой точностью.As RU 6 can be used, for example, a well-known measuring amplifier (instrumental), made on three operational amplifiers, in particular, made on a TL082 chip, manufactured by Texas Instruments. The switchgear 6 is connected through the first ADC 10 to the controller 7, which is covered by feedback with the possibility of automatically adjusting the amplitude level of the signal of the signal 1 by means of amplification factors if the amplitude level of the digitized signal signal 1 does not correspond to the specified values of the signal amplitude level range (stored in RAM 9). Thus, the feedback allows you to automatically adapt the amplitude level of the PP 1 signal to the range of amplitudes with the maximum possible discreteness required by the hardware, which makes it possible to obtain the specified signal in digital form with high accuracy.

Физический интерфейс 8 соединен с контроллером 7 и предназначен для передачи данных по каналам связи на внешние устройства потребителя. По физическому интерфейсу 8 передают следующие данные: идентификатор устройства, номер устройства, версия программного обеспечения, статус выполнения команд, накопленный объем, частота вибрации, коэффициенты усиления для регулирования амплитуды РУ 6, температура, время работы и адрес вихревого расходомера, скорость физического интерфейса 8, минимальная и максимальная измеряемые частоты сигнала. Передача данных, в частности, может быть осуществлена, например, посредством геофизического кабеля, подключаемого к физическому интерфейсу 8 на диспетчерский пункт и/или станцию управления и др. В качестве физического интерфейса 8 в данной реализации использован, в частности, цифровой физический интерфейс UART. В других вариантах исполнения в качестве физического интерфейса 8 может быть использован любой другой известный интерфейс, в частности CAN, RS-232, RS-485, USB.The physical interface 8 is connected to the controller 7 and is designed to transmit data via communication channels to external consumer devices. The following data is transmitted via physical interface 8: device identifier, device number, software version, command execution status, accumulated volume, vibration frequency, amplification factors for controlling the amplitude of RU 6, temperature, operating time and address of the vortex flowmeter, speed of the physical interface 8, minimum and maximum measured signal frequencies. Data transmission, in particular, can be carried out, for example, by means of a geophysical cable connected to a physical interface 8 to a control room and / or a control station, etc. As a physical interface 8 in this implementation, in particular, a digital physical UART interface is used. In other embodiments, any other known interface, in particular CAN, RS-232, RS-485, USB, can be used as the physical interface 8.

ОЗУ 9 служит для хранения оцифрованных сигнала ПП 1 и сигнала ДВ 3, хранения их спектров, характеристик, например, таких каких: амплитуда, частота, величина постоянной составляющей, хранения сигналов ПП 1 и ДВ 3, хранения вычисленных данных, хранения параметров работы расходомера в текущий момент времени, промежуточных вычисленных данных, не требующих длительного хранения, но необходимых для вычисления на следующих этапах вычисления, хранения заданных допустимых значений диапазонов уровня амплитуды сигнала.RAM 9 is used to store the digitized signal PP 1 and the signal DV 3, store their spectra, characteristics, for example, such as: amplitude, frequency, magnitude of the DC component, store signals PP 1 and DV 3, store the calculated data, store the parameters of the flowmeter in the current moment of time, intermediate calculated data that does not require long-term storage, but necessary for calculation at the next stages of calculation, storage of the set allowable values of signal amplitude level ranges.

В данном варианте исполнения ОЗУ 9, первый АЦП 10, второй АЦП 11 выполнены в составе контроллера 7 и могут быть реализованы, в частности, на микросхеме ATMEGA168-20AU, фирмы MICROCHIP. В других реализациях, второй АЦП 11 может быть выполнен, например, в составе ДВ 3, в частности, ДВ марки AIS1120SXTR, изготовитель STM. Кроме того, первый АЦП 10, второй АЦП 11, ОЗУ 9 могут быть выполнены в виде отдельных микросхем или устройств. В этом случае, в качестве первого АЦП 10 и второго АЦП 11 могут быть использованы любые широко известные аналогово-цифровые преобразователи, например марки MCP3201-CI/P, производитель Microchip, а в качестве ОЗУ - любой широкоизвестный ОЗУ, в частности ОЗУ марки AT25DF041A-SSH-В, производитель Atmel.In this embodiment, the RAM 9, the first ADC 10, the second ADC 11 are made as part of the controller 7 and can be implemented, in particular, on the microcircuit ATMEGA168-20AU, MICROCHIP. In other implementations, the second ADC 11 can be performed, for example, as part of the DV 3, in particular, the AIS1120SXTR DV brand manufactured by STM. In addition, the first ADC 10, the second ADC 11, RAM 9 can be made in the form of separate chips or devices. In this case, as the first ADC 10 and the second ADC 11, any widely known analog-to-digital converters can be used, for example, the MCP3201-CI / P brands manufactured by Microchip, and any widely known RAM, in particular the AT25DF041A- RAM, can be used as RAM. SSH-B, manufacturer of Atmel.

Для подавления спектральных компонент усиленный сигнал ПП 1 с выхода РУ 6, подаваемый на вход первого АЦП 10, пропускают через фильтр ФНЧ 12, чтобы частота усиленного сигнала ПП 1 превышала половину частоты дискретизации первого АЦП 10. Это позволяет исключить наложение спектров. В качестве ФНЧ 12 может быть использован широко известный anti-aliasing (антиалиасинговый) фильтр, в частности RC-фильтр.To suppress the spectral components, the amplified signal of PP 1 from the output of RU 6 supplied to the input of the first ADC 10 is passed through the low-pass filter 12 so that the frequency of the amplified signal of PP 1 exceeds half the sampling frequency of the first ADC 10. This eliminates the overlap of spectra. As the low-pass filter 12, a well-known anti-aliasing filter, in particular an RC filter, can be used.

ДВ 3 предназначен для непосредственного измерения вибрации. Присутствующими в сигнале ДВ 3 прочими шумовыми составляющими (например, помехами, наводками и пр.) пренебрегают, не учитывая их в расчете, вследствие того, что они имеет ничтожно малую амплитуду по сравнению с полезным сигналом. ДВ 3 может быть выполнен, в частности, на микросхеме ADIS 16006, производитель Analog Devices. Кроме того, в качестве ДВ 3 может быть использован, в частности, датчик вибрации, выполненный с возможностью осуществления функции оцифровки измеренного сигнала, например, датчик вибрации марки AIS1120SXTR, производитель STM. В этом случае на выходе ДВ 3 получают оцифрованный сигнал, поступающий в контроллер 7, минуя второй АЦП 11.DV 3 is intended for direct measurement of vibration. Other noise components present in the DV 3 signal (for example, interference, interference, etc.) are neglected, not taking them into account, because they have a negligible amplitude compared to the useful signal. DV 3 can be performed, in particular, on the chip ADIS 16006, manufacturer Analog Devices. In addition, as the DV 3 can be used, in particular, a vibration sensor configured to digitize the measured signal, for example, a vibration sensor brand AIS1120SXTR, manufacturer STM. In this case, at the output of the DV 3 receive a digitized signal supplied to the controller 7, bypassing the second ADC 11.

Вихревой расходомер выполнен с возможностью подачи питающего напряжения на требующие питания элементы. Питание вихревого расходомера осуществляется от внешнего источника «от +20 В до +60 В». На УАРУС 4 ПП подают два стабилизированных напряжения питания «+3.3 В», при этом первое стабилизированное напряжение питания по шине питания (не показана) поступает на УУ 5, РУ 6, ФНЧ 12, а второе - по шине питания (не показана) поступает на УО 2 и ДВ 3, а также на первый АЦП 10, ОЗУ 9, контроллер 7, физический цифровой интерфейс 8, второй АЦП 11.The vortex flowmeter is configured to supply voltage to the elements requiring power. The vortex flowmeter is powered from an external source “from +20 V to +60 V”. Two stabilized supply voltages “+3.3 V” are supplied to UARUS 4 PP, while the first stabilized supply voltage via the power bus (not shown) is supplied to UU 5, RU 6, low-pass filter 12, and the second is supplied via the power bus (not shown) on UO 2 and DV 3, as well as on the first ADC 10, RAM 9, controller 7, physical digital interface 8, and the second ADC 11.

Осуществление способа определения расхода жидкости будет показано на работе вихревого расходомера, установленного в скважине в непосредственной близости от погружного электродвигателя (ПЭД), например, на выкиде насоса. Расходомер работает следующим образом.The implementation of the method for determining the flow rate of the liquid will be shown on the work of the vortex flowmeter installed in the well in the immediate vicinity of the submersible electric motor (SEM), for example, on the pump discharge. The flow meter operates as follows.

Чувствительным элементом ПП 1 улавливают вихри, образовывающиеся за телом обтекания, а также вибрацию, воздействующую на расходомер. При одноканальном исполнении ПП 1 с его выхода в текущий момент времени получают сигнал, например, со следующими характеристиками: сигнал, в котором преобладают две частоты - 21 Гц и 19 Гц, амплитудой- 20 мВ, среднеквадратичным значением сигнала - 1600 мВ. Указанный сигнал поступает в УАРУС 4 ПП, где его амплитуда усиливается до значения 1300 мВ (на выходе УАРУС 4 ПП).The sensing element PP 1 capture the vortices formed behind the body flow around, as well as the vibration acting on the flow meter. In single-channel execution of software 1, a signal is received from its output at the current time, for example, with the following characteristics: a signal in which two frequencies prevail - 21 Hz and 19 Hz, amplitude - 20 mV, rms signal value - 1600 mV. The specified signal enters the UARUS 4 PP, where its amplitude is amplified to a value of 1300 mV (at the output of the UARUS 4 PP).

При двухканальном исполнении ПП 1 (исключает синфазные помехи) с выхода ПП 1 получают прямой и инверсный сигналы (отличаются между собой по фазе), каждый из которых имеет преобладающие частоты 21 Гц и 19 Гц, амплитудой 20 мВ, среднеквадратичным значением- 1600 мВ. Указанные сигналы поступают в УАРУС 4 ПП, где их амплитуда усиливается до значения 1300 мВ (остальные параметры не изменяются). При этом амплитуда прямого сигнала ПП 1 усиливается первым усилителем заряда, а амплитуда инверсного сигнала- вторым усилителем заряда УУ 5, подключенным параллельно первому усилителю заряда. Усиленные прямой сигнал ПП 1 и инверсный сигнал ПП 1 поступают на вход РУ 6, где за счет их разницы на его выходе формируется один (общий) усиленный сигнал ПП 1, в частности с преобладающими частотами - 21 Гц и 19 Гц, амплитудой - 1300 мВ, среднеквадратичным значением - 1600 мВ. Таким образом устраняют синфазную помеху.In the two-channel design of the software 1 (eliminates common mode interference), the direct and inverse signals (differ in phase) are received from the output of the software 1, each of which has the prevailing frequencies of 21 Hz and 19 Hz, with an amplitude of 20 mV, and a rms value of 1600 mV. These signals are fed to UARUS 4 PP, where their amplitude is amplified to a value of 1300 mV (other parameters do not change). The amplitude of the direct signal PP 1 is amplified by the first charge amplifier, and the amplitude of the inverse signal by the second charge amplifier UU 5, connected in parallel with the first charge amplifier. The amplified direct signal PP 1 and the inverse signal PP 1 are input to RU 6, where, due to their difference, one (common) amplified signal PP 1 is generated at its output, in particular with prevailing frequencies - 21 Hz and 19 Hz, amplitude - 1300 mV , the rms value is 1600 mV. This eliminates common mode interference.

С выхода УАРУС 4 ПП усиленный сигнал ПП 1 оцифровывают посредством первого АЦП 10, при этом в данном варианте первый АЦП 10 имеет следующие параметры: разрядность 10 бит, диапазон входных напряжений от 0 до 3300 мВ, частота дискретизации 1200 Гц, длительность оцифровки сигнала 430 мсек. На выходе первого АЦП 10 получают оцифрованный сигнал ПП 1 с количеством оцифрованных значений сигнала равным 1536 сэмплов. Усреднением каждых N соседних оцифрованных значений, в нашем случае N=3 соседних значений сигнала, получают оцифрованный сигнал ПП 1 с количеством оцифрованных значений сигнала равным 512 сэмплов, тем самым повышают его точность. При этом осуществляем передискретизацию первого АЦП 10 путем деления частоты дискретизации на N соседних оцифрованных значений сигнала. После передискретизации АЦП получаем частоту дискретизации АЦП равной 400 Гц.From the output of UARUS 4 PP, the amplified signal of PP 1 is digitized by the first ADC 10, and in this embodiment, the first ADC 10 has the following parameters: 10 bit capacity, input voltage range from 0 to 3300 mV, sampling frequency 1200 Hz, signal sampling duration 430 ms . The output of the first ADC 10 receives a digitized signal PP 1 with the number of digitized signal values equal to 1536 samples. By averaging every N neighboring digitized values, in our case N = 3 neighboring signal values, we get a digitized signal PP 1 with the number of digitized signal values equal to 512 samples, thereby increasing its accuracy. At the same time, we resample the first ADC 10 by dividing the sampling frequency by N adjacent digitized signal values. After oversampling the ADC, we obtain an ADC sampling frequency of 400 Hz.

После чего оцифрованный сигнал ПП 1 поступает в ОЗУ 9, где он сохраняется. Для исключения наложения спектров перед оцифровкой усиленный сигнал ПП 1 фильтруют посредством ФНЧ 12, в частности RC-цепочкой.After that, the digitized signal PP 1 enters RAM 9, where it is stored. To avoid overlapping spectra before digitization, the amplified signal PP 1 is filtered by the low-pass filter 12, in particular, an RC chain.

По запросу контроллера 7 сохраненный в ОЗУ 9 сигнал ПП 1 поступает в контроллер 7, где его уровень сравнивается (постоянно в каждый момент времени) с заранее заданными допустимыми диапазонами значений сигнала, при которых уровень сигнала будет находиться в пределах заданного диапазона входных напряжений первого АЦП 10. В частности, в данном варианте максимальное значение уровня сигнала ПП 1 должно находится в диапазоне от 2200 до 2750 мВ, минимальное значение уровня сигнала - в диапазоне от 550 до 1100 мВ.At the request of the controller 7, the signal of the software 1 stored in RAM 9 is sent to the controller 7, where its level is compared (constantly at each moment of time) with the predetermined allowable ranges of signal values at which the signal level will be within the specified input voltage range of the first ADC 10 In particular, in this embodiment, the maximum value of the signal level PP 1 should be in the range from 2200 to 2750 mV, the minimum value of the signal level should be in the range from 550 to 1100 mV.

Максимальное значение уровня сигнала ПП 1 можно вычислить как сумму среднеквадратичного значения сигнала и его амплитуды. Минимальное значение уровня сигнала ПП 1 можно вычислить как разность среднеквадратичного значения сигнала и его амплитуды. В данном случае, с учетом того, что текущая амплитуда сигнала ПП 1 равна 1300 мВ, а среднеквадратичное значение - 1600 мВ, минимальное значение уровня сигнала ПП 1 составит 300 мВ, что ниже, чем нижняя граница предела заданного диапазона, а максимальное значение уровня сигнала ПП 1 составит 2900 мВ, что превышает верхнюю границу предела заданного диапазона (т.е. имеет место несоответствие, которое может привести к некорректной оцифровке сигнала). В этом случае (при не соответствии) принимается решение о необходимости корректировки уровня сигнала ПП 1 до соответствия заданным допустимым диапазонам значений сигнала. Контроллер 7 посредством обратной связи РУ 6 автоматически регулирует уровень напряжения сигнала ПП 1 в текущий момент времени с помощью коэффициентов усиления, (ранее сохраненный сигнал ПП 1 при этом отбрасывается). Процесс повторяют вплоть до тех пор, пока уровень напряжения сигнала ПП 1 не установится в пределах заданных значений диапазона. В данном примере уровень амплитуды регулируют с помощью коэффициента усиления 0,8. В результате на выходе УАРУС 4 ПП получают сигнал с максимальным значением сигнала равным 2320 мВ. Это значение соответствует требуемому диапазону. Тем самым ослабляются наведенные помехи (шум) по отношению к полезному сигналу. Затем поступивший оптимизированный сигнал оцифровывается с помощью первого АЦП 10, и затем сохраняется в ОЗУ 9. Цикл повторяется до тех пор, пока вихревой расходомер включен и находится в режиме измерения.The maximum value of the signal level PP 1 can be calculated as the sum of the rms value of the signal and its amplitude. The minimum value of the signal level PP 1 can be calculated as the difference of the rms value of the signal and its amplitude. In this case, taking into account the fact that the current amplitude of the PP 1 signal is 1300 mV and the rms value is 1600 mV, the minimum value of the PP 1 signal level will be 300 mV, which is lower than the lower limit of the limit of the specified range, and the maximum value of the signal level PP 1 will be 2900 mV, which exceeds the upper limit of the limit of the specified range (i.e., there is a mismatch that can lead to incorrect digitization of the signal). In this case (if not compliance), a decision is made on the need to adjust the signal level PP 1 to match the specified allowable ranges of signal values. The controller 7 by means of feedback RU 6 automatically adjusts the voltage level of the signal PP 1 at the current time using the amplification factors (the previously saved signal PP 1 is discarded). The process is repeated until the voltage level of the signal PP 1 is not established within the specified values of the range. In this example, the amplitude level is adjusted using a gain of 0.8. As a result, at the output of UARUS 4 PP receive a signal with a maximum signal value of 2320 mV. This value corresponds to the required range. This reduces the induced noise (noise) with respect to the desired signal. Then the incoming optimized signal is digitized using the first ADC 10, and then stored in RAM 9. The cycle is repeated until the vortex flowmeter is turned on and is in measurement mode.

Посредством ДВ 3 измеряют вибрацию. Полученный сигнал ДВ 3 оцифровывают либо посредством второго АЦП 11 либо посредством ДВ 3 с функцией оцифровки, раскладывают в спектр и сохраняют в ОЗУ 9. Получение сигнала ДВ 3 может быть осуществлено как одновременно с приемом и оцифровкой сигнала ПП 1 (при условии, что значение величины вибрации постоянно изменяется во времени), так и отдельно от него в любой момент времени (при условии, что величина вибрации иметь неизменное/постоянное/статичное значение в любой момент времени). В данном примере измеренный сигнал ДВ имеет следующие характеристики: частота сигнала ДВ - 21 Гц, амплитуда - 500 мВ, среднеквадратичное значение - 1500 мВ.By means of DV 3 vibration is measured. The received signal of DV 3 is digitized either by means of the second ADC 11 or by DV 3 with the digitizing function, decomposed into a spectrum and stored in RAM 9. The signal of DV 3 can be received both simultaneously with the reception and digitization of signal PP 1 (provided that the value of vibration is constantly changing in time), and separately from it at any moment in time (provided that the magnitude of vibration has a constant / constant / static value at any time). In this example, the measured LF signal has the following characteristics: LF signal frequency - 21 Hz, amplitude - 500 mV, rms value - 1500 mV.

При параметрах второго АЦП И, в частности, разрядностью 10 бит, диапазон входных напряжений от 0 до 3300 мВ, частотой дискретизации - 1200 Гц, длительностью оцифровки сигнала 430 мсек., после оцифровки на выходе второго АЦП 11 получим оцифрованный сигнал ДВ с количеством оцифрованных значений сигнала равным 1536 сэмплов. Усреднением каждых N соседних оцифрованных значений сигнала, в нашем случае N=3 соседних значений, повышаем точность и получаем оцифрованный сигнал ДВ 3 с количеством оцифрованных значений сигнала равным 512 сэмплов, который затем сохраняем в ОЗУ 9. При этом осуществляем передискретизацию (понижают частоту дискретизации) второго АЦП 11 путем деления частоты дискретизации на N соседних оцифрованных значений сигнала. В результате получают частоту дискретизации второго АЦП 11 равную 400 Гц.With the parameters of the second ADC AND, in particular, with a resolution of 10 bits, the input voltage range is from 0 to 3300 mV, the sampling frequency is 1200 Hz, the signal digitization duration is 430 ms., After digitizing at the output of the second ADC 11, we obtain a digitized DV signal with the number of digitized values signal equal to 1536 samples. By averaging every N adjacent digitized signal values, in our case N = 3 neighboring values, we increase the accuracy and obtain a digitized signal DV 3 with the number of digitized signal values equal to 512 samples, which we then store in RAM 9. At the same time, we resample (lower the sampling frequency) the second ADC 11 by dividing the sampling frequency by N adjacent digitized signal values. The result is a sampling frequency of the second ADC 11 equal to 400 Hz.

Сохраненные в ОЗУ 9 сигнал ПП 1 и сигнал ДВ 3 считывают контроллером 7 (по его запросу) и раскладывают в спектры известными способами, например, с помощью быстрого преобразования Фурье (чтобы оценить наличие частотных составляющих в. сигнале, характеризующих наличие определенных частот в предлагаемых сигналах). В частности, оцифрованный сигнал ПП 1 раскладывают в спектр (фиг. 3), содержащий гармоники в количестве, зависящем от частоты дискретизации соответствующего АЦП. Например, в нашем примере при частоте дискретизации первого АЦП 10 равной 1200 Гц, спектр сигнала ПП 1 будет содержать 256 гармоник.The PP 1 signal and the DV 3 signal stored in RAM 9 are read by the controller 7 (upon request) and laid out in the spectra by known methods, for example, using the fast Fourier transform (to evaluate the presence of frequency components in the signal that characterize the presence of certain frequencies in the proposed signals ) In particular, the digitized signal PP 1 is laid out in a spectrum (Fig. 3) containing harmonics in an amount depending on the sampling frequency of the corresponding ADC. For example, in our example, at the sampling frequency of the first ADC 10 equal to 1200 Hz, the spectrum of the signal of PP 1 will contain 256 harmonics.

Аналогичным способом оцифрованный сигнал ДВ 3 раскладывают в спектр, содержащий от 1 до N гармоник, в частности, с помощью Быстрого преобразования Фурье. Полученный спектр представлен на фиг. 4 и в рассматриваемом примере, при частоте дискретизации второго АЦП 11 равной 1200 Гц, спектр имеет 256 гармоники. Сравнивают гармоники спектра сигнала ДВ 3 между собой, выделяют наибольшую по мощности гармонику, соответствующую преобладающей частоте в сигнале ДВ 3 (т.е. частоте вибрации, которую принимают (считают) равной частотной шумовой составляющей сигнала ПП 1). Определяют ее порядковый номер для дальнейшего вычисления частоты вибрации и исключения соответствующей частоты из спектра с ПП 1. В данном примере выделена наиболее мощная гармоника с порядковым номером 27.In a similar way, the digitized signal of DW 3 is decomposed into a spectrum containing from 1 to N harmonics, in particular, using the Fast Fourier transform. The resulting spectrum is shown in FIG. 4 and in this example, with a sampling frequency of the second ADC 11 equal to 1200 Hz, the spectrum has 256 harmonics. The harmonics of the spectrum of the DV 3 signal are compared with each other, the highest harmonic in power corresponding to the prevailing frequency in the DV 3 signal is selected (i.e., the vibration frequency that is assumed to be equal to the frequency noise component of the signal 1). Its serial number is determined for further calculation of the vibration frequency and the exclusion of the corresponding frequency from the spectrum with PP 1. In this example, the most powerful harmonic with serial number 27 is highlighted.

Выделяют спектр полезного сигнала путем вычитания наибольшей по мощности гармоники спектра сигнала ДВ 3 из спектра сигнала ПП 1 (при этом две одинаковые гармоники (вибрация) гасятся и, остается очищенный от вибрации спектр полезного сигнала (фиг. 5). Затем к выделенному спектру полезного сигнала (сигнал, очищенный от вибраций) применяют адаптивную фильтрацию, а именно, посредством последовательного использования окна Хэмминга и оценки весовой функции. При наличии шумовых составляющих, на данном этапе улучшится соотношение сигнал\шум, то есть гармоники, отвечающие за шум, станут менее мощными в сравнении с частотой вихреобразования. В результате получают спектр полезного сигнала (фиг. 6) с улучшенным соотношением частота-шум, тем самым повышается точность нахождения полезной частоты. Сравнивают гармоники выделенного спектра полезного сигнала между собой, определяют наибольшую по мощности гармонику (в данном случае максимальная по мощности гармоника имеет порядковый номер 24), по которой находят текущее значение частоты вихреобразования.The spectrum of the useful signal is isolated by subtracting the highest harmonic of the spectrum of the signal of the DW 3 signal from the spectrum of the signal PP 1 (in this case, two identical harmonics (vibration) are suppressed and the spectrum of the useful signal cleared from vibration remains (Fig. 5). Then, to the selected spectrum of the useful signal (a signal cleared of vibrations) adaptive filtering is used, namely, by sequentially using the Hamming window and estimating the weight function.When there are noise components, the signal-to-noise ratio will improve at this stage, that is, the harmonics responsible for the noise will become less powerful in compared with the vortex formation frequency. As a result, a useful signal spectrum is obtained (Fig. 6) with an improved frequency-to-noise ratio, thereby improving the accuracy of finding the useful frequency. The harmonics of the selected useful signal spectrum are compared with each other, and the maximum power harmonic is determined (in this case, the maximum in power, the harmonic has serial number 24), according to which find the current value of the vortex formation frequency.

Текущее значение частоты вихреобразования может быть вычислено, например, по формуле:The current value of the vortex formation frequency can be calculated, for example, by the formula:

Fтек=Xi*(H/N),Ftec = Xi * (H / N),

где Fтек - Текущее значение частоты вихреобразования, Гцwhere Ftek - Current value of the vortex frequency, Hz

Xi - порядковый номер бина спектра, при этом i принимает значение от 1 до N, где N-соответствует максимальному порядковому номеру гармоники (в нашем случае 256), которое определяется параметрами оцифровки сигнала.Xi is the serial number of the spectrum bin, and i takes a value from 1 to N, where N-corresponds to the maximum harmonic serial number (in our case 256), which is determined by the signal digitization parameters.

Н - Частота дискретизации АЦП (Гц);N - ADC sampling frequency (Hz);

N - Количество точек сигнала (отсчет сигнала - численное значение напряжения сигнала в определенный момент времени).N - The number of signal points (signal count - the numerical value of the signal voltage at a certain point in time).

Fтек=24*(400/512)=18.75 ГцFtec = 24 * (400/512) = 18.75 Hz

После чего определяют скорость потока жидкости (объем за текущий момент времени). В частности, скорость потока жидкости в текущий момент времени может быть определена по формуле:Then determine the fluid flow rate (volume for the current time). In particular, the fluid flow rate at the current time can be determined by the formula:

Vтек= Fтек*D/Sh, гдеVtek = Ftek * D / Sh, where

Vтек - скорость потока жидкости в текущий момент времени (л/час);Vtek - fluid flow rate at the current time (l / h);

Fтек - текущая частота вихреобразования (Гц);Ftek - current vortex frequency (Hz);

D - ширина тела обтекания (мм);D is the width of the body flow (mm);

Sh - число Струхаля.Sh is the Strouhal number.

При ширине тела обтекания равном 7 мм, и известным значением числа Струкаля равным 0,295 преобразовываем частоту вихреобразования в скорость потока жидкости следующим образом:With the width of the flow body equal to 7 mm, and the known value of the Strukal number equal to 0.295, we transform the vortex frequency into the fluid flow rate as follows:

Vтек=18.75 Гц * 7 мм/0.295=444.9 л/часVtec = 18.75 Hz * 7 mm / 0.295 = 444.9 l / h

Далее вычисленная скорость потока жидкости в текущий момент времени корректируется с помощью коррелирующего коэффициента (заранее заданного, учитывающего отклонение от эталона, хранится в ОЗУ 9), тем самым сглаживается погрешность в техническом исполнении модели устройства.Further, the calculated fluid flow rate at the current time is corrected using a correlation coefficient (predefined, taking into account the deviation from the standard, stored in RAM 9), thereby smoothing the error in the technical performance of the device model.

444.9 л/ч * 1.005=447.139 л/час444.9 l / h * 1.005 = 447.139 l / h

Вычисленный расход жидкости передается через интерфейс 8 потребителю.The calculated fluid flow rate is transmitted via the interface 8 to the consumer.

Способ и устройство обеспечивают повышение точности определения частоты вихреобразования, особенно при сильном вибрационном воздействии за счет: повышения вибрационной устойчивости вихревого расходомера, возможности исключения ложного измерения расхода и увеличения минимального измеряемого уровня расхода в условиях воздействия на расходомер динамически изменяющихся вибрационных нагрузок - повышение чувствительности расходомера с соотношения сигнал/шум (расход/вибрация) с 6 до -6 дБ, вызванных работой двигателей, насосов, клапанов, увеличения рабочего диапазона работы расходомера со значения 1:10 до 1:20 в сторону меньших значений расхода.The method and device provide improved accuracy in determining the frequency of vortex formation, especially with strong vibration due to: increased vibration stability of the vortex flowmeter, the possibility of eliminating false flow measurements and increasing the minimum measured flow rate under conditions of dynamically changing vibration loads acting on the flowmeter - increasing the sensitivity of the flowmeter from the ratio signal / noise (flow / vibration) from 6 to -6 dB caused by the operation of motors, pumps, valves, increasing the operating range of the flow meter from 1:10 to 1:20 to lower flow rates.

Понижение частоты дискретизации сигнала (программным способом) путем усреднения нескольких ближайших точек захвата позволяет повысить точность оцифровки сигнала и уменьшить погрешность дискретизации.Decreasing the sampling frequency of the signal (programmatically) by averaging several nearest capture points allows you to increase the accuracy of digitization of the signal and reduce the sampling error.

Таким образом, несмотря на то, что заявленное в качестве изобретения техническое решение показано и описано на примере его конкретной реализации, среднему специалисту в данной области должно быть понятно, что различные изменения по форме и содержанию могу быть сделаны без отступления от сущности и объема изобретения, определенных приведенной в заявке формуле изобретения.Thus, despite the fact that the technical solution claimed as an invention is shown and described by the example of its specific implementation, it should be clear to a person skilled in the art that various changes in form and content can be made without departing from the essence and scope of the invention, defined in the application claims.

Claims (12)

1. Способ определения расхода жидкости, характеризующийся тем, что с первичного преобразователя (ПП) получают сигнал, включающий полезную частотную составляющую, характеризующую скорость потока жидкости и шумовую частотную составляющую, характеризующую вибрацию в текущий момент времени, полученный сигнал ПП оцифровывают, раскладывают в спектр, при этом посредством датчика вибрации (ДВ) измеряют вибрацию в текущий момент времени, измеренный сигнал ДВ оцифровывают и раскладывают в спектр, выделяют спектр полезного сигнала путем вычитания из спектра сигнала ПП спектр сигнала ДВ, на основании выделенного спектра полезного сигнала вычисляют значение частоты вихреобразования, по которой определяют расход жидкости в текущий момент времени.1. The method of determining the flow rate of the liquid, characterized in that the primary Converter (PP) receive a signal that includes a useful frequency component that characterizes the fluid flow rate and a noise frequency component that characterizes the vibration at the current time, the received PP signal is digitized, laid out in the spectrum, while using a vibration sensor (VL), vibration is measured at the current moment of time, the measured LW signal is digitized and laid out in the spectrum, the spectrum of the useful signal is extracted by subtracting the spectrum of the LW signal from the spectrum of the PP signal, based on the extracted spectrum of the useful signal, calculate the vortex frequency, which determine the flow rate at the current time. 2. Способ по п. 1, отличающийся тем, что к выделенному спектру полезного сигнала применяют адаптивную фильтрацию, при этом значение частоты вихреобразования вычисляют после ее применения.2. The method according to p. 1, characterized in that adaptive filtering is applied to the selected spectrum of the useful signal, while the value of the vortex formation frequency is calculated after its application. 3. Способ по п. 2, отличающийся тем, что упомянутую адаптивную фильтрацию осуществляют посредством последовательного применения Окна Хемминга и оценки весовых функций.3. The method according to p. 2, characterized in that the said adaptive filtering is carried out by successive application of the Hamming Window and the evaluation of weight functions. 4. Способ по п. 1, отличающийся тем, что перед оцифровкой уровень полученного сигнала ПП автоматически регулируют.4. The method according to p. 1, characterized in that before digitizing the level of the received signal PP automatically adjust. 5. Способ по п. 4, отличающийся тем, что после автоматической регулировки уровня сигнала ПП осуществляют фильтрацию его нижних частот.5. The method according to p. 4, characterized in that after automatic adjustment of the signal level PP carry out filtering of its low frequencies. 6. Способ по любому из предшествующих пунктов, отличающийся тем, что сигнал ДВ оцифровывают одновременно с оцифровыванием сигнала ПП.6. The method according to any one of the preceding paragraphs, characterized in that the DV signal is digitized simultaneously with the digitization of the PP signal. 7. Вихревой расходомер, реализующий способ по п. 1, выполненный с возможностью подачи питающего напряжения на требующие питания элементы и содержащий соединенные последовательно первичный преобразователь (ПП), узел автоматической регулировки уровня сигнала ПП (УАРУС ПП), устройство обработки (УО), датчик вибрации (ДВ), при этом УО выполнен с возможностью разложения в спектры оцифрованных сигнала ПП и измеренного сигнала ДВ, выделения спектра полезного сигнала, применения к этому выделенному спектру адаптивной фильтрации посредством Окна Хемминга и оценки весовых функций.7. A vortex flowmeter that implements the method according to claim 1, configured to supply power to power-demanding elements and comprising a primary converter (PP) connected in series, an automatic signal level adjustment unit PP (UARUS PP), a processing device (UO), a sensor vibration (LW), while the RL is capable of expanding into the spectra of the digitized PP signal and the measured LF signal, extracting the spectrum of the useful signal, applying adaptive filtering to this selected spectrum by means of the Hamming Window and estimating weight functions. 8. Вихревой расходомер по п. 7, отличающийся тем, что УО содержит контроллер, физический интерфейс, ОЗУ, с которым последовательно соединены первый аналогово-цифровой преобразователь (первый АЦП), второй аналогово-цифровой преобразователь (второй АЦП), причем ДВ соединен последовательно через второй АЦП с контроллером.8. The vortex flow meter according to claim 7, characterized in that the UO comprises a controller, a physical interface, RAM, to which the first analog-to-digital converter (first ADC) is connected in series, the second analog-to-digital converter (second ADC), and the DV is connected in series through the second ADC with the controller. 9. Вихревой расходомер по п. 1, отличающийся тем, что УАРУС ПП включает последовательно соединенные друг с другом усилительный узел (УУ), регулирующий усилитель (РУ), с которым связан через обратную связь контроллер УО с возможностью регулирования уровня амплитуды сигнала ПП, при этом РУ, в свою очередь, соединен с первым АЦП УО.9. The vortex flowmeter according to claim 1, characterized in that the UARUS PP includes an amplifying unit (UI) connected in series with each other, a regulating amplifier (RU), to which the UO controller is connected via feedback with the ability to control the amplitude level of the PP signal, when this RU, in turn, is connected to the first ADC UO. 10. Вихревой расходомер по п. 9, отличающийся тем, что УУ содержит усилитель заряда.10. The vortex flowmeter according to claim 9, characterized in that the SU contains a charge amplifier. 11. Вихревой расходомер по п. 9, отличающийся тем, что УУ содержит, по крайней мере, два усилителя заряда, соединенные параллельно друг с другом, подключенные последовательно с РУ.11. The vortex flowmeter according to claim 9, characterized in that the control unit comprises at least two charge amplifiers connected in parallel with each other, connected in series with the switchgear. 12. Вихревой расходомер по п. 9, отличающийся тем, что УАРУС ПП дополнительно содержит фильтр низких частот (ФНЧ), включенный между РУ и первым АЦП.12. The vortex flow meter according to claim 9, characterized in that the UARUS PP further comprises a low-pass filter (LPF) included between the switchgear and the first ADC.
RU2019145775A 2019-12-31 2019-12-31 Method of determining liquid flow rate and vortex flow meter for its implementation RU2726275C9 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
RU2019145775A RU2726275C9 (en) 2019-12-31 2019-12-31 Method of determining liquid flow rate and vortex flow meter for its implementation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
RU2019145775A RU2726275C9 (en) 2019-12-31 2019-12-31 Method of determining liquid flow rate and vortex flow meter for its implementation

Publications (2)

Publication Number Publication Date
RU2726275C1 true RU2726275C1 (en) 2020-07-10
RU2726275C9 RU2726275C9 (en) 2020-10-05

Family

ID=71510042

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2019145775A RU2726275C9 (en) 2019-12-31 2019-12-31 Method of determining liquid flow rate and vortex flow meter for its implementation

Country Status (1)

Country Link
RU (1) RU2726275C9 (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2192623C2 (en) * 2000-11-13 2002-11-10 Мартынов Евгений Васильевич Multi-purpose vortex-type flowmeter
US20120192657A1 (en) * 2011-01-31 2012-08-02 Krohne Messtechnik Gmbh Vortex flowmeter
CN101701834B (en) * 2009-11-10 2013-01-23 合肥工业大学 Strong interference resistant digital vortex shedding flowmeter based on double sensors
CN105067048B (en) * 2015-07-17 2018-02-13 合肥科迈捷智能传感技术有限公司 A kind of measuring method of the digitlization vortex-shedding meter of anti-vibration interference

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2192623C2 (en) * 2000-11-13 2002-11-10 Мартынов Евгений Васильевич Multi-purpose vortex-type flowmeter
CN101701834B (en) * 2009-11-10 2013-01-23 合肥工业大学 Strong interference resistant digital vortex shedding flowmeter based on double sensors
US20120192657A1 (en) * 2011-01-31 2012-08-02 Krohne Messtechnik Gmbh Vortex flowmeter
CN105067048B (en) * 2015-07-17 2018-02-13 合肥科迈捷智能传感技术有限公司 A kind of measuring method of the digitlization vortex-shedding meter of anti-vibration interference

Also Published As

Publication number Publication date
RU2726275C9 (en) 2020-10-05

Similar Documents

Publication Publication Date Title
JP6740363B2 (en) Device, system, and method for measuring internal impedance of a test battery using frequency response
AU2011236558B2 (en) System and method of detecting cavitation in pumps
TW442655B (en) Laplace transform impedance spectrometer and its measurement method
WO2017143649A1 (en) Kalman filter-based vortex flowmeter anti-transient shock interference signal processing method and system
JP5300174B2 (en) Jitter measuring apparatus, jitter measuring method, test apparatus, and electronic device
US20100000879A1 (en) Electrochemical noise as a localized corrosion indicator
US10459013B2 (en) Apparatus for and method of providing measurements of uncertainty in respect of a transfer function
US20140303913A1 (en) Broken rotor bar detection based on current signature analysis of an electric machine
AU2010271044B2 (en) Determining the resonance parameters for mechanical oscillators
CN104296786B (en) Digital bridge capacitive measuring module
CN110763903B (en) Residual current detection method, device and circuit
CN107124179B (en) Phase-locked amplifier for detecting weak photocurrent signal
US11262429B2 (en) Method, apparatus and device for detecting abnormity of energy metering chip, and medium
CN102043084A (en) Method for detecting resistive current of lightning arrester
RU2726275C1 (en) Method of determining liquid flow rate and vortex flow meter for its implementation
RU2014153145A (en) CHECK THE CURRENT CONTROL OF THE PROCESS
US6973839B2 (en) Electromagnetic flow meter having processing means resolving output data into a non-flow waveform component
CN108680210B (en) Transient electromagnetic flow transmitter based on voltage and current differential
CN107270954B (en) A kind of vibrating string type sensor performance evaluation device and method
CN115932355A (en) Self-adaptive signal processing system and method of clamp-on ammeter
Kam et al. Simple schemes for measuring autocorrelation functions
CN109254251B (en) Battery impedance measuring device, method and chip
RU2653150C1 (en) Method of finding-out of information image of electrical signal
AU2020407219A1 (en) Gating energy consumption accumulation by detecting a fundamental component of a current
US9759751B1 (en) Line cycle correlated spectral analysis for power measurement systems

Legal Events

Date Code Title Description
TK4A Correction to the publication in the bulletin (patent)

Free format text: CORRECTION TO CHAPTER -FG4A- IN JOURNAL 19-2020 FOR INID CODE(S) (72)

TH4A Reissue of patent specification