RU2726275C9 - 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
RU2726275C9
RU2726275C9 RU2019145775A RU2019145775A RU2726275C9 RU 2726275 C9 RU2726275 C9 RU 2726275C9 RU 2019145775 A RU2019145775 A RU 2019145775A RU 2019145775 A RU2019145775 A RU 2019145775A RU 2726275 C9 RU2726275 C9 RU 2726275C9
Authority
RU
Russia
Prior art keywords
signal
spectrum
flow rate
vibration
vortex
Prior art date
Application number
RU2019145775A
Other languages
Russian (ru)
Other versions
RU2726275C1 (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, namely to methods and tools for measuring the flow of liquids, and can be used in various sectors of the national economy for the purpose of monitoring, regulating and accounting for medium flows. In particular, the invention can be used to determine the flow rate of a fluid in a well with 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 liquid in a well, a vortex flow meter is installed in it in the immediate vicinity of a submersible electric motor (SEM), for example, at the pump outlet. The speed of the SEM operation, as a rule, is in the range from 10 to 80 Hz. In this case, vibration occurs at these and multiple frequencies. This vibration load is transmitted through the tubing string to nearby equipment, including a vortex meter. In cases where 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 consequence, the vibration frequency can be mistakenly perceived by the flow meter 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 in determining the flow rate.

Наиболее близким является способ определения расхода жидкости, реализованный в универсальном вихревом расходомере (патент РФ на изобретение №2192623, МПК G01F 1/32, 2000 г.), содержащем первичный преобразователь и электрическую схему обработки сигнала первичного преобразователя, которая содержит последовательно включенные первый нормирующий усилитель, фильтр, второй нормирующий усилитель и аналого-цифровой преобразователь, микроконтроллер, управляющие выходы которого подключены соответственно к управляющим входам первого и второго нормирующих усилителей, интерфейсы расходомера. Функцию адаптивно перестраиваемых фильтров выполняет алгоритм цифровой адаптивной фильтрации, реализованный в микроконтроллере. Способ осуществляется следующим образом. Получают сигнал с первичного преобразователя, производит выборку сигнала с частотой опроса, необходимой для восстановления сигнала в диапазоне измеряемых частот, оценивают уровень оцифрованного входного сигнала и выполняют автоматическую регулировку усиления, производит поиск и выделение сигнала вихревой частоты в соответствии с алгоритмом цифровой адаптивной фильтрации, вычисляют текущее значение вихревой частоты, вычисляют значение расхода. Способ позволяет определить скорость или расход жидкости жидких и газообразных сред в скважине при воздействии возмущений в виде вибрации трубопровода, пульсаций потока и акустических шумов.The closest is a method for determining the flow rate of a liquid, implemented in a universal vortex flow meter (RF patent for invention No. 2192623, IPC G01F 1/32, 2000), containing a primary converter and an electrical circuit for processing the signal of the primary converter, which contains the first normalizing amplifier connected 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, the flow meter interfaces. The function of adaptively tunable filters is performed by the digital adaptive filtering algorithm implemented in the microcontroller. The method is carried out as follows. Receive a signal from the primary transducer, samples the signal with a sampling frequency required to restore the signal in the range of measured frequencies, estimate the level of the digitized input signal and perform automatic gain control, search for and extract the vortex frequency signal in accordance with the digital adaptive filtering algorithm, calculate the current vortex frequency value, the flow rate value is calculated. The method makes it possible to determine the speed or flow rate of liquid of liquid and gaseous media in a well under the influence of disturbances in the form of pipeline vibration, flow pulsations and acoustic noise.

Недостатками способа является низкая точность и достоверность определения расхода жидкости, в том числе при низкочастотных вибрационных нагрузках, поскольку способ не позволяет выделить вибрацию при смешении вибрации и полезного сигнала, характеризующего расход жидкости и исключить ее при определении расхода жидкости.The disadvantages of this 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 to isolate vibration when mixing vibration and a 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 invention No. 2192623, IPC G01F 1/32, 2000), containing a primary converter and an electrical circuit for processing the signal of the primary converter, 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 to the control inputs of the first and second normalizing amplifiers, respectively, the flow meter 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 at low-frequency vibration loads. In addition, as the fluid flow rate decreases, the useful flow signal is weakened, and the useful signal becomes commensurate with the lower harmonics. This limits the lower flow measurement of the flow meter.

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

Технический результат заключается в повышении точности определения расхода жидкости за счет применения способа определения расхода жидкости в скважине, в том числе при низкочастотных вибрационных нагрузках.The technical result consists in increasing the accuracy of determining the fluid flow rate by using 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 flow rate of a liquid by making it possible to measure (isolate) vibration and the ability to eliminate it when determining the flow rate of a liquid, as well as to expand the functionality, which consists in expanding the working range of the vortex flow meter in the area of low flow rates.

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

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

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

Предпочтительно, чтобы адаптивную фильтрацию осуществляли посредством последовательного применения к выделенному спектру полезного сигнала Окна Хемминга и оценки весовых функций.It is preferable that the adaptive filtering is carried out by sequentially applying a Hamming Window to the selected spectrum of the wanted signal and estimating the weighting functions.

Предпочтительно, чтобы перед оцифровкой уровень амплитуды полученного сигнала ПП автоматически регулировали.It is preferable that the amplitude level of the received PCB signal is automatically adjusted prior to digitization.

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

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

Кроме того, указанная цель достигается в вихревом расходомере, выполненном с возможностью подачи питающего напряжения на требующие питания элементы, и содержащем соединенные последовательно первичный преобразователь (далее, ПП), устройство обработки (далее, УО), датчик вибрации (далее, ДВ). При этом УО выполнен с возможностью разложения в спектры оцифрованных сигнала ПП и сигнала ДВ, выделения спектра полезного сигнала из спектра сигнала ПП, применения к упомянутому выделенному спектру полезного сигнала адаптивной фильтрации в виде Окна Хемминга и оценки весовых функций, вычисления значения частоты вихреобразования, вычисления значения скорости потока жидкости в текущий момент времени, автоматической регулировки уровня сигнала ПП, фильтрации его нижних частот.In addition, this goal is achieved in a vortex flow meter configured to supply a supply voltage to the elements requiring power, and containing a primary transducer (hereinafter referred to as PP), a processing device (hereinafter referred to as UO), and a vibration sensor (hereinafter referred to as DV) connected in series. In this case, the UO is configured to decompose the digitized PP signal and DW signal into spectra, extract the spectrum of the useful signal from the spectrum of the PP signal, apply adaptive filtering to the mentioned selected spectrum of the useful signal in the form of the Hamming Window and estimate the weight functions, calculate the value of the vortex formation frequency, calculate the value fluid flow rate at the current time, automatic adjustment of the PP signal level, filtering its low frequencies.

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

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

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

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

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

Далее на фиг. 1 структурная схема предлагаемого вихревого расходомера, на фиг. 2 приведен пример предлагаемого вихревого расходомера, который реализует предлагаемый способ определения расхода жидкости в предпочтительном варианте, на фиг. 3 показан спектр сигнала ПП, включающий гармонику вибрации и гармонику вихреобразования, на фиг. 4 показан спектр сигнала ДВ, на фиг. 5 показан спектр полезного сигнала, фиг. 6 показан спектр выделенного полезного сигнала после адаптивной фильтрации.Further, in FIG. 1 is a block diagram of the proposed vortex flow meter, in Fig. 2 shows an example of the proposed vortex flowmeter, which implements the proposed method for determining the flow rate of a liquid in a preferred embodiment, FIG. 3 shows the spectrum of the SP signal, including the vibration harmonic and the vortex harmonic; FIG. 4 shows the spectrum of the LW signal; FIG. 5 shows the spectrum of the useful signal, FIG. 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 serially connected PP 1, UO 2, DV 3, while the sensitive element PP 1 is a piezoelectric element that catches the vortices formed behind the streamline body. UO 2 (Fig. 2) contains UARUS 4 PP (shown by a dashed line), including serially connected 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 RU 6 is connected in series with the first ADC 10. To improve accuracy, reduce noise between RU 6 and the first ADC 10, a low-pass filter (LPF) 12 can be additionally connected. Controller 7 is covered by feedback with RU 6, which provides the ability to automatically adjust the gains to match the ranges of the output voltage PP 1 and the input voltage of the first ADC 10. DV 3 is connected through the second ADC 11 to controller 7.

На фиг. 3 показан спектр сигнала ПП, в котором гармоники смешаны, а также в отдельности гармоника вибрации (вибрационная шумовая составляющая) и гармоника вихреобразования (составляющая полезного сигнала). На фигуре 4 показан спектр сигнала ДВ, в котором выделена максимальная гармоника, характеризующая интенсивность вибрации. На фигуре 5 показан спектр полезного сигнала. В частности, спектр полезного сигнала может быть получен программным исключением вибрационной гармоники из спектра сигнала ПП. Для более точного выделения гармоники, отвечающей за скорость потока жидкости, а также для улучшения соотношения сигнал-шум к выделенному спектру полезного сигнала применяют адаптивную фильтрация в виде последующих друг за другом окна Хемминга и оценки весовых функций. Результат применения адаптивной фильтрации представлен на фигуре 6 и отражает увеличенное соотношение сигнал-шум, благодаря чему появляется возможность более точно получить основную частоту вихреобразования (в спектре полезного сигнала), по которой определяют скорость потока (расход) жидкости.FIG. 3 shows the spectrum of the SP signal, in which the harmonics are mixed, as well as separately the vibration harmonic (vibration noise component) and the vortex harmonic (useful signal component). Figure 4 shows the spectrum of the DV signal, in which the maximum harmonic characterizing the vibration intensity is highlighted. Figure 5 shows the spectrum of the wanted signal. In particular, the spectrum of the useful signal can be obtained by software exclusion of the vibration harmonic from the spectrum of the PP signal. To more accurately isolate the harmonic responsible for the fluid flow rate, as well as to improve the signal-to-noise ratio to the selected spectrum of the useful signal, adaptive filtering is used in the form of successive Hamming windows and weighting functions estimates. The result of applying adaptive filtering is shown in figure 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 is used to determine 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 can be used, located behind the flow body in the flow path of a vortex flowmeter installed in the well in the tubing rupture. The sensing element PP 1 catches the vortices formed behind the flow body and converts them into an analog signal characterizing the intensity of vortex formation and, as a consequence, the fluid flow rate. This signal is complex and includes a useful frequency component characterizing the fluid flow rate and a noise frequency component characterizing the vibration intensity at the current time. PP 1 is connected by UO 2.

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

- оцифровку сигнала ПП 1 и сигнала ДВ 3 и их сохранение в ОЗУ 9;- digitization of the signal PP 1 and 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 PP 1 signal and the DV 3 signal into spectra;

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

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

- вычисления значения частоты вихреобразования,- calculating the value of the vortex formation 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, controller 7 (in particular, microcontroller), to which the physical interface 8 is connected, RAM 9, to which, in turn, the first ADC 10 is connected in series, intended for converting the signal PP 1 into digital form and the second ADC 11, intended for converting the measured signal DV 3 into digital form.

УАРУС 4 ПП предназначен для предварительной обработки сигнала ПП 1 под соответствие сигнала ПП 1 диапазону первого АЦП 10 и включает последовательно соединенные УУ 5, РУ 6. Если принятый сигнал ПП 1 сразу удовлетворяет диапазону первого АЦП 10, сигнал ПП 1 может быть подан сразу на первый АЦП 10, минуя УАРУС 4 ПП, как показано на фиг. 2.UARUS 4 PP is designed for preprocessing the PP 1 signal for the correspondence of the PP 1 signal to the range of the first ADC 10 and includes serially connected UU 5, RU 6. If the received PP 1 signal immediately satisfies the range of the first ADC 10, the PP 1 signal can be applied directly 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 intended for preliminary amplification of the PP 1 signal and includes at least one charge amplifier, in particular, a current integrator (built on the basis of the current integration method and implemented on an operational amplifier), in particular, made on a TL082 microcircuit, Texas Instruments ... In particular, a two-channel PP 1 can be used to filter out the symphase interferences resulting from the operation of the electronic devices of the vortex flow meter (at the output of the PP 1 we have direct and inverse signals). In this case, UU 5 can turn on the first charge amplifier, providing amplification of the direct signal PP 1, in parallel to which the second charge amplifier is connected, providing amplification of the inverse signal of the PP 1. From the output of the UU 5, the amplified direct signal and the inverse signal are fed to the input of RU 6, where for due to their difference at its output, one (common) amplified signal PP 1 is formed, 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 instrumentation amplifier (instrumental), made on three operational amplifiers, in particular, made on the TL082 microcircuit, manufactured by Texas Instruments. RU 6 is connected through the first ADC 10 with the controller 7, which is covered with feedback with the possibility of automatic regulation of the amplitude level of the signal PP 1 by means of amplification factors in case of a discrepancy between the amplitude level of the digitized signal PP 1 and the specified values of the signal amplitude level range (stored in RAM 9). Thus, the feedback makes it possible to automatically adapt the amplitude level of the PP 1 signal to the range of amplitude values required by the hardware with the maximum possible discreteness, 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.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 RU 6 amplitude, temperature, operating time and address of the vortex flow meter, speed of 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 the physical interface 8 to the control room and / or control station, etc. As the physical interface 8 in this implementation, in particular, the digital physical interface UART is used. In other embodiments, any other known interface can be used as physical interface 8, in particular CAN, RS-232, RS-485, USB.

ОЗУ 9 служит для хранения оцифрованных сигнала ПП 1 и сигнала ДВ 3, хранения их спектров, характеристик, например, таких каких: амплитуда, частота, величина постоянной составляющей, хранения сигналов ПП 1 и ДВ 3, хранения вычисленных данных, хранения параметров работы расходомера в текущий момент времени, промежуточных вычисленных данных, не требующих длительного хранения, но необходимых для вычисления на следующих этапах вычисления, хранения заданных допустимых значений диапазонов уровня амплитуды сигнала.RAM 9 serves for storing the digitized signal PP 1 and signal DV 3, storing their spectra, characteristics, for example, such as: amplitude, frequency, value of the constant component, storing signals PP 1 and DV 3, storing calculated data, storing the operating parameters of the flow meter in the current moment of time, intermediate computed data that do not require long-term storage, but are necessary for computation at the next stages of computation, storage of the specified allowable values of the 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 ATMEGA168-20AU microcircuit, manufactured by MICROCHIP. In other implementations, the second ADC 11 can be made, for example, as part of the DV 3, in particular, the DV of the AIS1120SXTR brand, manufactured by STM. In addition, the first ADC 10, the second ADC 11, RAM 9 can be made in the form of separate microcircuits or devices. In this case, any widely known analog-to-digital converters can be used as the first ADC 10 and the second ADC 11, for example, MCP3201-CI / P, manufactured by Microchip, and as RAM, any well-known RAM, in particular, AT25DF041A- SSH-B, manufacturer Atmel.

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

ДВ 3 предназначен для непосредственного измерения вибрации. Присутствующими в сигнале ДВ 3 прочими шумовыми составляющими (например, помехами, наводками и пр.) пренебрегают, не учитывая их в расчете, вследствие того, что они имеет ничтожно малую амплитуду по сравнению с полезным сигналом. ДВ 3 может быть выполнен, в частности, на микросхеме ADIS 16006, производитель Analog Devices. Кроме того, в качестве ДВ 3 может быть использован, в частности, датчик вибрации, выполненный с возможностью осуществления функции оцифровки измеренного сигнала, например, датчик вибрации марки AIS1120SXTR, производитель STM. В этом случае на выходе ДВ 3 получают оцифрованный сигнал, поступающий в контроллер 7, минуя второй АЦП 11.DV 3 is designed for direct vibration measurement. Other noise components present in the DV 3 signal (for example, interference, pickups, etc.) are neglected without taking them into account, due to the fact that they have a negligible amplitude compared to the useful signal. DV 3 can be made, in particular, on the ADIS 16006 microcircuit, manufactured by Analog Devices. In addition, as DV 3 can be used, in particular, a vibration sensor configured to perform the function of digitizing the measured signal, for example, a vibration sensor of the AIS1120SXTR brand, manufactured by STM. In this case, at the output of DV 3, a digitized signal is received, which is fed to 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 a supply voltage to the elements requiring power. The vortex flowmeter is powered from an external source "from +20 V to +60 V". UARUS 4 PP is supplied with two stabilized supply voltages "+3.3 V", while the first stabilized supply voltage is supplied through the power bus (not shown) to UU 5, RU 6, LPF 12, and the second - through 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, second ADC 11.

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

Чувствительным элементом ПП 1 улавливают вихри, образовывающиеся за телом обтекания, а также вибрацию, воздействующую на расходомер. При одноканальном исполнении ПП 1 с его выхода в текущий момент времени получают сигнал, например, со следующими характеристиками: сигнал, в котором преобладают две частоты - 21 Гц и 19 Гц, амплитудой- 20 мВ, среднеквадратичным значением сигнала - 1600 мВ. Указанный сигнал поступает в УАРУС 4 ПП, где его амплитуда усиливается до значения 1300 мВ (на выходе УАРУС 4 ПП).The sensing element PP 1 captures the vortices formed behind the body of the streamline, as well as the vibration affecting the flow meter. With a single-channel version of the PP 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, an amplitude of 20 mV, and an rms signal value of 1600 mV. This 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 мВ. Таким образом устраняют синфазную помеху.With a two-channel version of PP 1 (excludes common-mode interference), direct and inverse signals are received from the output of PP 1 (they differ in phase), each of which has dominant frequencies of 21 Hz and 19 Hz, amplitude 20 mV, rms value - 1600 mV. These signals are sent to UARUS 4 PP, where their amplitude is amplified to a value of 1300 mV (other parameters do not change). In this case, the amplitude of the direct signal PP 1 is amplified by the first charge amplifier, and the amplitude of the inverse signal is amplified 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 fed to the input of RU 6, where, due to their difference, one (common) amplified signal PP 1 is formed at its output, in particular with the prevailing frequencies - 21 Hz and 19 Hz, amplitude - 1300 mV , rms value - 1600 mV. This eliminates common mode noise.

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

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

По запросу контроллера 7 сохраненный в ОЗУ 9 сигнал ПП 1 поступает в контроллер 7, где его уровень сравнивается (постоянно в каждый момент времени) с заранее заданными допустимыми диапазонами значений сигнала, при которых уровень сигнала будет находиться в пределах заданного диапазона входных напряжений первого АЦП 10. В частности, в данном варианте максимальное значение уровня сигнала ПП 1 должно находится в диапазоне от 2200 до 2750 мВ, минимальное значение уровня сигнала - в диапазоне от 550 до 1100 мВ.At the request of the controller 7, the PP signal 1 stored in the RAM 9 enters the controller 7, where its level is compared (constantly at each time point) with the predetermined allowable signal ranges at which the signal level will be within the specified range of input voltages 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 - 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 between the rms value of the signal and its amplitude. In this case, taking into account 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 specified range, and the maximum value of the signal level PP 1 will be 2900 mV, which exceeds the upper limit of the specified range (i.e., there is a discrepancy that can lead to incorrect digitization of the signal). In this case (in case of non-compliance), a decision is made on the need to adjust the signal level of the PP 1 to comply with the specified allowable signal ranges. The controller 7, through the feedback of RU 6, automatically regulates the voltage level of the PP 1 signal at the current time using the amplification factors, (the previously saved PP 1 signal is discarded in this case). The process is repeated up until the voltage level of the signal PP 1 is established within the specified range values. In this example, the amplitude level is adjusted using a gain of 0.8. As a result, a signal with a maximum signal value equal to 2320 mV is obtained at the output of the UARUS 4 PP. This value is within the required range. This reduces the induced interference (noise) in relation to the desired signal. Then, the received 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 the measurement mode.

Посредством ДВ 3 измеряют вибрацию. Полученный сигнал ДВ 3 оцифровывают либо посредством второго АЦП 11 либо посредством ДВ 3 с функцией оцифровки, раскладывают в спектр и сохраняют в ОЗУ 9. Получение сигнала ДВ 3 может быть осуществлено как одновременно с приемом и оцифровкой сигнала ПП 1 (при условии, что значение величины вибрации постоянно изменяется во времени), так и отдельно от него в любой момент времени (при условии, что величина вибрации иметь неизменное/постоянное/статичное значение в любой момент времени). В данном примере измеренный сигнал ДВ имеет следующие характеристики: частота сигнала ДВ - 21 Гц, амплитуда - 500 мВ, среднеквадратичное значение - 1500 мВ.Vibration is measured by DV 3. The received signal DV 3 is digitized either by means of the second ADC 11 or by means of DV 3 with the digitization function, decomposed into a spectrum and stored in RAM 9. Receiving the DV 3 signal can be carried out both simultaneously with the reception and digitization of the PP 1 signal (provided that the value vibration constantly changes in time), and separately from it at any moment of time (provided that the magnitude of vibration has a constant / constant / static value at any moment). In this example, the measured LW signal has the following characteristics: LW 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, the bit width is 10 bits, the input voltage range is from 0 to 3300 mV, the sampling frequency is 1200 Hz, the signal digitization duration is 430 m, after digitizing at the output of the second ADC 11, we will receive a digitized DV signal with the number of digitized signal values equal to 1536 samples. By averaging each N neighboring 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 is then stored in RAM 9. At the same time, we carry out oversampling (lower the sampling frequency) second ADC 11 by dividing the sampling rate 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 гармоник.Stored in RAM 9, the PP signal 1 and the DV signal 3 are read by the controller 7 (at its request) and decomposed into spectra using known methods, for example, using the fast Fourier transform (to assess the presence of frequency components in the signal, characterizing the presence of certain frequencies in the proposed signals ). In particular, the digitized signal PP 1 is decomposed into a spectrum (Fig. 3), containing harmonics in an amount that depends on the sampling frequency of the corresponding ADC. For example, in our example, with the sampling rate of the first ADC 10 equal to 1200 Hz, the spectrum of the PP 1 signal will contain 256 harmonics.

Аналогичным способом оцифрованный сигнал ДВ 3 раскладывают в спектр, содержащий от 1 до N гармоник, в частности, с помощью Быстрого преобразования Фурье. Полученный спектр представлен на фиг. 4 и в рассматриваемом примере, при частоте дискретизации второго АЦП 11 равной 1200 Гц, спектр имеет 256 гармоники. Сравнивают гармоники спектра сигнала ДВ 3 между собой, выделяют наибольшую по мощности гармонику, соответствующую преобладающей частоте в сигнале ДВ 3 (т.е. частоте вибрации, которую принимают (считают) равной частотной шумовой составляющей сигнала ПП 1). Определяют ее порядковый номер для дальнейшего вычисления частоты вибрации и исключения соответствующей частоты из спектра с ПП 1. В данном примере выделена наиболее мощная гармоника с порядковым номером 27.In a similar way, the digitized signal DV 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 the 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 power harmonic corresponding to the prevailing frequency in the DV 3 signal is isolated (i.e., the vibration frequency, which is taken (considered) equal to the frequency noise component of the PP 1 signal). Its serial number is determined to further calculate the vibration frequency and exclude 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 selected by subtracting the highest power harmonic of the spectrum of the signal DV 3 from the spectrum of the signal PP 1 (while two identical harmonics (vibration) are suppressed and the spectrum of the useful signal, free from vibration, remains (Fig. 5). Then, to the selected spectrum of the useful signal (signal cleaned of vibrations) adaptive filtering is applied, namely, through sequential use of the Hamming window and the estimation of the weighting function.In the presence of noise components, at this stage the signal-to-noise ratio will improve, that is, the harmonics responsible for the noise will become less powerful in as a result, a spectrum of the useful signal (Fig. 6) with an improved frequency-to-noise ratio is obtained, thereby increasing the accuracy of finding the useful frequency. The harmonics of the selected spectrum of the useful signal are compared with each other, the highest power harmonic is determined (in this case, the maximum in terms of power, the harmonic has a 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),Ftek = Xi * (H / N),

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

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

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

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

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

После чего определяют скорость потока жидкости (объем за текущий момент времени). В частности, скорость потока жидкости в текущий момент времени может быть определена по формуле:After that, the fluid flow rate is determined (volume at 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тек - скорость потока жидкости в текущий момент времени (л/час);Vtk is the fluid flow rate at the current time (l / h);

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

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

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

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

Vтек=18.75 Гц * 7 мм/0.295=444.9 л/часVtk = 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 (predetermined, taking into account the deviation from the standard, is stored in RAM 9), thereby smoothening 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 liquid flow rate is transmitted through the interface 8 to the consumer.

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

Понижение частоты дискретизации сигнала (программным способом) путем усреднения нескольких ближайших точек захвата позволяет повысить точность оцифровки сигнала и уменьшить погрешность дискретизации.Downsampling of a signal (programmatically) by averaging several nearby capture points improves the signal's digitization accuracy and reduces the sampling error.

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

Claims (12)

1. Способ определения расхода жидкости, характеризующийся тем, что с первичного преобразователя (ПП) получают сигнал, включающий полезную частотную составляющую, характеризующую скорость потока жидкости и шумовую частотную составляющую, характеризующую вибрацию в текущий момент времени, полученный сигнал ПП оцифровывают, раскладывают в спектр, при этом посредством датчика вибрации (ДВ) измеряют вибрацию в текущий момент времени, измеренный сигнал ДВ оцифровывают и раскладывают в спектр, выделяют спектр полезного сигнала путем вычитания из спектра сигнала ПП спектр сигнала ДВ, на основании выделенного спектра полезного сигнала вычисляют значение частоты вихреобразования, по которой определяют расход жидкости в текущий момент времени.1. A method for determining the flow rate of a liquid, characterized in that a signal is received from a primary transducer (PP), which includes a useful frequency component characterizing the fluid flow rate and a noise frequency component characterizing vibration at the current time, the received PP signal is digitized, decomposed into a spectrum, at the same time, by means of a vibration sensor (DV), the vibration is measured at the current time, the measured DV signal is digitized and decomposed into a spectrum, the spectrum of the useful signal is separated by subtracting the spectrum of the DV signal from the spectrum of the PP signal, based on the selected spectrum of the useful signal, the value of the vortex formation frequency is calculated, according to which determine the flow rate of the liquid at the current time. 2. Способ по п. 1, отличающийся тем, что к выделенному спектру полезного сигнала применяют адаптивную фильтрацию, при этом значение частоты вихреобразования вычисляют после ее применения.2. The method according to claim 1, characterized in that adaptive filtering is applied to the selected spectrum of the useful signal, and the vortex frequency is calculated after its application. 3. Способ по п. 2, отличающийся тем, что упомянутую адаптивную фильтрацию осуществляют посредством последовательного применения Окна Хемминга и оценки весовых функций.3. A method according to claim 2, characterized in that said adaptive filtering is performed by sequentially applying the Hamming Window and estimating the weighting functions. 4. Способ по п. 1, отличающийся тем, что перед оцифровкой уровень полученного сигнала ПП автоматически регулируют.4. The method according to claim 1, characterized in that the level of the received PP signal is automatically adjusted before digitizing. 5. Способ по п. 4, отличающийся тем, что после автоматической регулировки уровня сигнала ПП осуществляют фильтрацию его нижних частот.5. The method according to claim 4, characterized in that after automatic adjustment of the level of the PP signal, its low frequencies are filtered. 6. Способ по любому из предшествующих пунктов, отличающийся тем, что сигнал ДВ оцифровывают одновременно с оцифровыванием сигнала ПП.6. A method according to any of the preceding claims, characterized in that the LW signal is digitized simultaneously with the digitization of the PP signal. 7. Вихревой расходомер, реализующий способ по п. 1, выполненный с возможностью подачи питающего напряжения на требующие питания элементы и содержащий соединенные последовательно первичный преобразователь (ПП), узел автоматической регулировки уровня сигнала ПП (УАРУС ПП), устройство обработки (УО), датчик вибрации (ДВ), при этом УО выполнен с возможностью разложения в спектры оцифрованных сигнала ПП и измеренного сигнала ДВ, выделения спектра полезного сигнала, применения к этому выделенному спектру адаптивной фильтрации посредством Окна Хемминга и оценки весовых функций.7. Vortex flowmeter implementing the method according to claim 1, configured to supply a supply voltage to the elements requiring power and containing a primary transducer (PP) connected in series, a unit for automatic signal level control PP (UARUS PP), a processing device (UO), a sensor vibration (DV), while the UO is configured to decompose the digitized PP signal and the measured DV signal into spectra, extract the spectrum of the useful signal, apply adaptive filtering to this selected spectrum by means of the Hamming Window and estimate the weight functions. 8. Вихревой расходомер по п. 7, отличающийся тем, что УО содержит контроллер, физический интерфейс, ОЗУ, с которым последовательно соединены первый аналогово-цифровой преобразователь (первый АЦП), второй аналогово-цифровой преобразователь (второй АЦП), причем ДВ соединен последовательно через второй АЦП с контроллером.8. Vortex flowmeter according to claim 7, characterized in that the UO contains a controller, a physical interface, a RAM with which the first analog-to-digital converter (first ADC), the second analog-to-digital converter (second ADC) are connected in series, 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 amplifier unit (UU) connected in series with each other, a control amplifier (RU), with which the UO controller is connected through feedback with the possibility of adjusting the amplitude of the PP signal, when this RU, in turn, is connected to the first ADC of the UO. 10. Вихревой расходомер по п. 9, отличающийся тем, что УУ содержит усилитель заряда.10. Vortex flowmeter according to claim 9, characterized in that the CU contains a charge amplifier. 11. Вихревой расходомер по п. 9, отличающийся тем, что УУ содержит, по крайней мере, два усилителя заряда, соединенные параллельно друг с другом, подключенные последовательно с РУ.11. Vortex flowmeter according to claim 9, characterized in that the control unit contains at least two charge amplifiers connected in parallel with each other, connected in series with the RU. 12. Вихревой расходомер по п. 9, отличающийся тем, что УАРУС ПП дополнительно содержит фильтр низких частот (ФНЧ), включенный между РУ и первым АЦП.12. Vortex flowmeter according to claim 9, characterized in that the UARUS PP additionally contains a low-pass filter (LPF) connected between the RU 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 RU2726275C1 (en) 2020-07-10
RU2726275C9 true 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
RU2726275C1 (en) 2020-07-10

Similar Documents

Publication Publication Date Title
US7634382B2 (en) Diagnostic device for use in process control system
CN108534852B (en) Electromagnetic flowmeter based on dynamic detection and detection method
TW297876B (en)
US20100000879A1 (en) Electrochemical noise as a localized corrosion indicator
EP1660844B1 (en) Flow meter filter system and method
US5424636A (en) Process and device for detecting revolutions of a rotating element by twice fourier transform
JP2006500683A (en) Adaptive spatial notch filter
JP3105465B2 (en) Voice section detection method
US11262429B2 (en) Method, apparatus and device for detecting abnormity of energy metering chip, and medium
US11788875B2 (en) Measurement pickup for determining the mass flow rate of a liquid
RU2601186C2 (en) Process control current circuit verification
JP6085910B2 (en) Photoelectric sensor and method for controlling amplification of received light amount in photoelectric sensor
RU2726275C9 (en) Method of determining liquid flow rate and vortex flow meter for its implementation
US6973839B2 (en) Electromagnetic flow meter having processing means resolving output data into a non-flow waveform component
CN109142642A (en) A kind of gas concentration quick sensing device and cognitive method
JP2021526218A (en) Vortex flowmeter to detect flow instability
CN113340369B (en) Signal processing method and device for turbine fuel mass flowmeter
CN114838924A (en) Structural damping ratio identification method based on wind-induced vibration non-stationary response
CN106059684B (en) A kind of noise processed system and method based on Correlation Identification
RU2653150C1 (en) Method of finding-out of information image of electrical signal
US20150293047A1 (en) Method and apparatus for determining water content of oil and water mixtures by measurement of specific admittance
RU162515U1 (en) Borehole Spectral Sound Level Meter
SU1065789A1 (en) Electrical machine air gap eccentricity indirect determination method
CN204421980U (en) Numerical model analysis filtering low flow velocity electromagnetic flowmeter
CN112664177B (en) Oil well working fluid level measuring method based on normalized convolution and adaptive filtering

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