СИСТЕМА И СПОСОБ ОБНАРУЖЕНИЯ НАРОСТА ОТЛОЖЕНИЙ В УЛЬТРАЗВУКОВОМ РАСХОДОМЕРЕ И МАШИНОЧИТАЕМЫЙ НОСИТЕЛЬ ИНФОРМАЦИИ
SYSTEM AND METHOD TO DETECT OUTGROWTH OF DEPOSITS IN ULTRASONIC FLOW METER AND MACHINE-READABLE INFORMATION MEDIUM

Область техники, к которой относится изобретение
Изобретение относится к способу и системе обнаружения и отслеживания отложений, формирующихся внутри ультразвукового расходомера.
Уровень техники
После добычи углеводородов из недр, поток текучей среды (например, сырая нефть или природный газ) перемещается из пункта добычи к потребителям по системе трубопроводов. При этом желательно знать с точностью количество протекающей в потоке текучей среды, а при транспортировки текучей среды потребителям, или при так называемом откачивании углеводородов с промысла потребителю по закрытой системе, требуется особая точность при определении количества текучей среды. Если последнее не имеет места, точность измерений тем не менее желательна.
В подобных случаях, таких как откачивание углеводородов с промысла потребителю по закрытой системе, возможно применение ультразвуковых расходомеров. В ультразвуковом расходомере имеют место обмен акустическими сигналами, проходящими через поток текучей среды, с последующим выполнением измерений между одной или более парами преобразователей. Каждая пара преобразователей заключена внутри корпуса расходомера или муфты так, что проходящий от одного преобразователя к другому акустический сигнал пересекает под некоторым углом протекающую через расходомер текучую среду. Группа электронных устройств, расположенных внутри корпуса расходомера выполнены с возможностью измерения разности между временем распространения сигнала, которое требуется акустическому сигналу для его прохождения от преобразователя, расположенного ниже по потоку, до преобразователя, расположенного выше по потоку, и временем распространения сигнала, которое требуется акустическому сигналу для его прохождения от преобразователя, расположенного выше по потоку, до преобразователя, расположенного ниже по потоку. Разность времен распространения акустического сигнала впоследствии используют для вычисления средней скорости и объемного расхода текучей среды, протекающей через расходомер.
При эксплуатации ультразвукового расходомера в трубопроводе на внутренней поверхности данного прибора может иметь место нарост отложений. Например, указанные ультразвуковые приборы для измерения расхода текучих сред могут применяться для измерения расхода сырой нефти в составе которой нередко присутствуют парафины. Таким образом, на внутренних поверхностях прибора со временем образуются наросты парафиновых отложений. Нарост таких отложений приводит к ухудшению точности измерений времени распространения акустических сигналов и вызывает уменьшение эффективной площади поперечного проходного сечения расходомера, тем самым внося погрешность в расчетное значение объемного расхода текучей среды, протекающей через расходомер.
Раскрытие изобретения
Указанные выше проблемы могут быть по крайней мере частично решены с помощью способа и системы для обнаружения нароста отложений внутри ультразвукового расходомера. По крайней мере некоторые из приведенных в качестве иллюстраций вариантов осуществления изобретения представляют собой системы, содержащие ультразвуковой расходомер и группу электронных устройств. Ультразвуковой расходомер содержит муфту, выполненную с возможностью взаимодействия с потоком текучей среды, первую пару преобразователей, механически закрепленных на муфте и расположенных так, что они находятся в контакте с потоком текучей среды, и группу электронных устройств. При этом первая пара преобразователей содержит преобразователь, находящийся выше по потоку, и преобразователь, находящийся ниже по потоку, связанный с преобразователем, находящимся выше по потоку, и определяет первую хорду между указанными преобразователями. Группа электронных устройств выполнена с возможностью получения диагностических данных, на основе информации об акустических сигналах, которыми обмениваются преобразователи первой пары, где диагностические данные в ключают скорость распространения звука через текучую среду. Указанная группа электронных устройств дополнительно выполнена с возможностью определения нароста отложений на внутренней поверхности муфты на основании линии основного направления скорости распространения звука через текучую среду.
Прочие приведенные в качестве иллюстраций варианты осуществления изобретения представляют собой способы получения диагностических данных, на основе информации об акустических сигналах, которыми обмениваются преобразователи пары преобразователей ультразвукового расходомера, в том числе скорости распространения звука через текучую среду, протекающую через ультразвуковой расходомер.
Иные варианты осуществления изобретения, дополнительные к вариантам, приведенным выше для иллюстрации, представляют собой машиночитаемые носители информации, которые включают блок группы машинных команд, связанный с процессором для определения нароста отложений внутри ультразвукового расходомера посредством диагностических данных и акустических сигналов, которые передают преобразователи пары преобразователей, при этом одним из указанных диагностических данных является скорость распространения звука через текучую среду, протекающую через ультразвуковой преобразователь.
В раскрытых в данном описании системах и способах могут сочетаться признаки и преимущества, позволяющие устранить недостатки, свойственные предшествующему уровню техники. Различные характеристики, указанные выше, а также другие признаки и особенности изобретения станут очевидными для специалистов в данной области техники после чтения нижеследующего подробного описания изобретения и просмотра сопровождающих чертежей.
Краткое описание чертежей
Подробное описание различных вариантов изобретения содержит ссылки на следующие сопровождающие чертежи:
На Фиг.1 представлен вид сверху в разрезе ультразвукового расходомера;
На Фиг.2 представлен вид сбоку ультразвукового расходомера, соответствующего по крайней мере нескольким вариантам изобретения; показаны муфта и хорды А-D;
На Фиг.3 представлен вид сверху ультразвукового расходомера, соответствующего по крайней мере нескольким вариантам изобретения; показаны корпус муфты и несколько пар преобразователей;
На Фиг.4 представлена блок-схема, иллюстрирующая вариант осуществления способа определения нароста отложений;
На Фиг.5 и 6 схематически представлены виды поперечных сечений в радиальном и осевом направлениях ультразвукового расходомера с наростом отложений; и
На Фиг.7 представлена блок-схема, иллюстрирующая вариант осуществления способа количественного определения нароста отложений и корректировки объемного расхода текучей среды с учетом нароста отложений.
Осуществление изобретения
Некоторые термины используются везде в следующих далее описании и формуле изобретения для ссылок на конкретные компоненты системы. В данном документе компоненты, отличные друг от друга по наименованию, но не по функциям, различаться не будут.
В следующем далее обсуждении и в формуле изобретения термин "включает" используется не в ограничительном смысле, и таким образом, должен пониматься как означающий "включая (но перечень не является ограничительным)…". Термин "соединять" и ему подобные означают как непосредственное, так и опосредованное соединение, подключение или связь. Таким образом, если говорится, что первое устройство соединено со вторым устройством, то соединение может быть как непосредственным, так и таким, которое опосредовано соединениями через другие устройства и соединители. Далее, термин "установлен" будет служить для обозначения как непосредственного, так и опосредованного соединения. Таким образом, если говорится, что первое устройство установлено на втором устройстве, то соединение установленного устройства может быть как непосредственным, так и опосредованным соединениями через другие устройства и соединители.
На Фиг.1 показан ультразвуковой расходомер, служащий для измерения расхода текучей среды, в соответствии с вариантами изобретения. Муфта 100, которая может быть установлена между двумя участками трубопровода, имеет предопределенные размеры, которыми определяется геометрия области, в которой проводятся измерения. Пара преобразователей 120 и 130 устанавливается в соответствующих держателях 125 и 135, разнесенных по длине муфты 100. Преобразователи 120 и 130 являются ультразвуковыми приемопередатчиками, то есть генерируют и принимают ультразвуковые сигналы. Термин "ультразвуковые сигналы" в данном контексте относится к акустическим сигналам, частота которых в некоторых вариантах превышают примерно 20 килогерц. В некоторых вариантах ультразвуковые сигналы могут иметь частоту приблизительно 125 килогерц (для газовых расходомеров) и 1 мегагерц (для жидкостных расходомеров). Ультразвуковые сигналы любой из этих частот могут генерироваться и приниматься пьезоэлектрическими элементами каждого преобразователя. Для генерации ультразвукового сигнала осуществляется воздействие электрического сигнала на пьезоэлектрический элемент, который отвечает вибрацией на это воздействие. При вибрации пьезоэлектрического элемента возбуждается ультразвуковой сигнал, который распространяется в продольно-поперечном направлении внутри муфты 100 через текучую среду, достигая соответствующего преобразователя, принадлежащего той же паре преобразователей. Ультразвуковой сигнал, достигший принимающего пьезоэлектрического элемента, заставляет его вибрировать, генерируя электрический сигнал, который детектируется, преобразуется в цифровую форму и анализируется группой электронных устройств, связанной с расходомером.
Пусть отрезок 110, который называют хордой, соединяет два преобразователя 120 и 130 и направлен под некоторым углом θ к оси 105. Длина хорды 110 это расстояние между передней поверхностью преобразователя 120 и передней поверхностью преобразователя 130. Точки 140 и 145 являются местами, где акустические сигналы, генерируемые преобразователями 120 и 130, входят в или выходят из текучей среды, протекающей через муфту 100. Положение преобразователей 120 и 130 может быть определено углом θ, первым расстоянием L, измеренным между преобразователями 120 и 130, вторым расстоянием X, представляющим собой расстояние между точками 140 и 145, измеренное вдоль оси, и третьим расстоянием D, соответствующим диаметру трубопровода или муфты. В большинстве случаев расстояния D, Х и L точно определяются при изготовлении расходомера. Кроме того, такие преобразователи как 120 и 130 обычно размещаются на определенном расстоянии от точек 140 и 145 соответственно; это расстояние не зависит от размеров расходомера, то есть от диаметра муфты.
Сначала нижний по потоку преобразователь 120 генерирует ультразвуковой сигнал, распространяющийся в направлении верхнего по потоку преобразователя 130 и достигающий его. Через некоторое время верхний по потоку преобразователь 130 генерирует ответный ультразвуковой сигнал, распространяющийся в направлении нижнего по потоку преобразователя 120 и достигающий его. Таким образом, обмен распространяющимися вдоль хорды 110 ультразвуковыми сигналами 115 между преобразователями 120 и 130 можно сравнить с подачей мяча и его захватом в ловушку. При работе расходомера каждая пара преобразователей может повторять такую последовательность действий тысячи раз в минуту.
Поток текучей среды в направлении 150 через муфту 100 характеризуется эпюрой 152 скоростей. Показанные на чертеже векторы скорости 153-158 демонстрируют тот факт, что скорость потока, проходящего через муфту 100, увеличивается при приближении к оси 105. Время прохождения ультразвукового сигнала 115 между преобразователями 120 и 130 зависит, в частности, от направления распространения (вверх или вниз по потоку) ультразвукового сигнала 115. Время прохождения ультразвукового сигнала 115 по потоку (то есть, в том же направлении, в котором движется поток) меньше, чем время прохождения сигнала в противоположном направлении (то есть, в направлении против потока). Значения времени прохождения сигнала в направлениях по потоку и против потока можно использовать для определения средней по всей хорде 110 скорости потока, а также для определения скорости звука в потоке текучей среды. При заданных размерах сечения расходомера, через который проходит текучая среда, и при известной средней скорости может быть рассчитан объем текучей среды, прошедшей через муфту 100 за заданное время.
Для более точного определения средней по всему сечению расходомера скорости в ультразвуковом расходомере могут использоваться несколько хорд; такие расходомеры называют многолучевыми. На Фиг.2 показан многолучевой ультразвуковой расходомер. В данных вариантах хорда А (обозначена как 225 на Фиг.2), хорда В (230), хорда С (235) и хорда D (240) в муфте 100 расположены в разных слоях потока. В альтернативных вариантах многолучевые расходомеры могут иметь разное число хорд. Каждая из хорд А-D соответствует двум преобразователям, работающим поочередно в качестве передатчика и в качестве приемника. Показана также группа электронных устройств 160 управления, которая получает и обрабатывает данные от пар преобразователей, соответствующих четырем хордам А-D. На Фиг.2 не показаны четыре пары преобразователей, соответствующие хордам А-D.
Схему расположения четырех пар преобразователей можно легко понять, обратившись к Фиг.3. Четыре пары входов/выходов преобразователей, или отверстий, расположены на муфте 100. Каждая пара отверстий преобразователей соответствует одной из хорд 110, показанных на Фиг.2. На муфте 100 расположена первая пара отверстий 125 и 135 преобразователей и соответствующие им преобразователи. Другая пара отверстий преобразователей, состоящая из отверстий 165 и 175 (видны частично), и соответствующие им преобразователи установлены таким образом, чтобы соединяющая их хорда была расположена, в общих чертах, Х-образно по отношению к хорде 110, соединяющей отверстия 125 и 135 преобразователей. Точно так же отверстия 185 и 195 преобразователей расположены параллельно отверстиям 165 и 175 преобразователей, но в других плоскостях (то есть, на разных радиальных расстояниях от оси трубопровода или муфты расходомера). На Фиг.3 не показаны явно четвертая пара отверстий преобразователей и соответствующие им преобразователи. Рассматривая совместно Фиг.2 и 3, можно видеть, что пары преобразователей размещены таким образом, что две верхние пары преобразователей, соответствующие хордам А и В, расположены Х-образно, и две нижних две пары преобразователей, соответствующих хордам С и D, также расположены Х-образно. На основе данных о времени прохождения сигналов может быть определена скорость потока текучей среды для каждой из хорд А-D, благодаря чему можно получить средние вдоль хорд скорости потока (далее - "скорости по хордам"); комбинируя скорости по хордам, можно определить скорость потока, среднюю по всему сечению трубопровода или муфты 100 расходомера. Объемный расход текучей среды, проходящей через муфту 100 расходомера, может быть далее определен как произведение средней скорости потока на площадь поперечного сечения муфты 100 расходомера.
i
1,i
2,i
Скорости потока, измеряемые вдоль хорд, получают по наборам значений времени распространения сигнала, получаемым от четырех пар преобразователей. Для хорды i набор значений времени распространения сигнала включает набор разностей ΔT между значениями времени распространения сигнала Т в направлении против потока и времени распространения сигнала T в направлении по потоку, получаемых, в основном, из следующего уравнения:
<mrow><mi>Δ</mi><mi>T</mi><mo>=</mo><msub><mi>T</mi><mrow><mn>1,</mn><mi>i</mi></mrow></msub><mo>−</mo><msub><mi>T</mi><mrow><mn>2,</mn><mi>i</mi></mrow></msub><mtext> </mtext><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>
i
i
i
Набор, состоящий из 20 значений ΔT, может использоваться для определения среднего значения ΔT. В альтернативных вариантах число значений ΔT может изменяться.
i
i
На основании среднего значения ΔT, может быть определена средняя вдоль хорды скорость V, потока, в основном, из следующего уравнения:
<mrow><msub><mi>V</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msubsup><mi>L</mi><mi>i</mi><mn>2</mn></msubsup></mrow><mrow><mn>2</mn><msub><mi>X</mi><mi>i</mi></msub></mrow></mfrac><mo>⋅</mo><mfrac><mrow><mi>Δ</mi><msub><mi>T</mi><mi>i</mi></msub></mrow><mrow><msub><mi>T</mi><mrow><mn>1,</mn><mi>i</mi></mrow></msub><msub><mi>T</mi><mrow><mn>2,</mn><mi>i</mi></mrow></msub></mrow></mfrac><mtext> </mtext><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>
i
A
B
C
D
i
i
i
где i - номер конкретной хорды, V - определяемая вдоль хорды скорость потока (то есть, V, V, V и V для хорд А-D соответственно), L - расстояние, или длина хорды, между преобразователями, и X - расстояние в потоке, измеряемое вдоль оси. Далее на основе средних определяемых вдоль хорд скоростей V может быть определена средняя скорость V потока текучей среды, проходящего через трубу или муфту 100 расходомера, в основном, с помощью следующего уравнения:
<mrow><mi>V</mi><mo>=</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder><mrow><msub><mi>W</mi><mi>i</mi></msub><msub><mi>V</mi><mi>i</mi></msub></mrow></mstyle><mtext> </mtext><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>
i
где W - весовой множитель, зависящий от номера i хорды.
Используя среднюю скорость V потока, следующим образом можно определить объемный расход Q текучей среды, проходящей через муфту 100 расходомера:
<mrow><mi>Q</mi><mo>=</mo><mi>V</mi><msub><mi>A</mi><mi>C</mi></msub><mtext> </mtext><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>
C
C
где A - площадь поперечного сечения муфты 100 расходомера. Площадь поперечного сечения A определяется:
<mrow><msub><mi>A</mi><mi>c</mi></msub><mo>=</mo><mfrac><mrow><mi>π</mi><msup><mi>D</mi><mn>4</mn></msup></mrow><mn>4</mn></mfrac><mtext> </mtext><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>
где D - внутренний диаметр муфты 100 расходомера, который обычно определяется при изготовлении муфты 100.
i
При работе ультразвукового расходомера группа электронных устройств, связанная с расходомером, выполняет группу функций, которые могут включать инициализацию работы преобразователя, получение выходного сигнала от одного или группы преобразователей и вычисление средней измеряемой вдоль хорды скорости V для каждой из хорд, средней скорости V потока и объемного расхода Q потока, проходящего через расходомер. Параметры, вычисленные группой электронных устройств могут затем передаваться электронной системе, удаленной от расходомера, например, системе диспетчерского управления и сбора данных (SCADA) и использоваться как входная информация для дальнейших вычислений. Альтернативно, все вычисления могут выполняться группой электронных устройств, связанной с расходомером, или электронными устройствами, удаленными от него.
i
1,i
2,i
C
1,i
2,i
i
C
При прохождении текучей среды через муфту расходомера на его внутренних поверхностях может начаться образование нароста отложений. Через некоторое время рост отложений может стать причиной возникновения погрешностей при вычислении средних вдоль хорд скоростей V, функций, зависящих от измеренных значений Т, T времени распространения сигнала, площади поперечного сечения A расходомера и, следовательно, к погрешностям в вычисленном значении объемного расхода Q потока текучей среды через расходомер. При образовании отложений акустический сигнал, генерируемый преобразователем, проходит через материал отложений, прежде чем быть принятым другим преобразователем. Таким образом, формирование нароста отложений влияет на измеряемые значения Т, Т времени распространения сигнала, приводя к погрешностям в вычисляемых значениях средней измеренной вдоль хорды скорости V. Отложения, далее, уменьшают эффективную площадь поперечного сечения A расходомера. Обе эти погрешности приводят к погрешностям при вычислении объемного расхода Q потока через расходомер.
Если нет возможности выявить и количественно оценивать рост отложений на внутренних поверхностях ультразвукового расходомера, то расходомер требует регулярного технического обслуживания, которое часто оказывается длительным и дорогостоящим дополнительным обслуживанием. Отметим также, что в данной области техники отсутствует какой-либо способ внесения достоверных поправок в вычисляемые значения объемного расхода, которые учитывали бы образование нароста отложений в периоды между техническим обслуживанием. Раскрываемое изобретение направлено на создание способов и систем обнаружения и количественной оценки формирования нароста отложений в ультразвуковых расходомерах. Раскрываемые способы и системы позволяют получать диагностические данные, как измеренные, так и вычисленные, позволяющие проводить техническое обслуживание, необходимое для устранения отложений, в основном, когда это действительно необходимо, а не в соответствии с жестким графиком. Если техническое обслуживание в данное время не выполняется, раскрытые здесь способы и системы позволят скорректировать вычисляемый объемный расход с учетом образования нароста отложений.
1,i
2,i
A
B
AB
В соответствии с некоторыми вариантами, погрешности в значениях времени распространения сигнала Т, Т и, следовательно, формирование нароста отложений, можно обнаружить, используя функцию η. Уровень техники и способ получения индикатора погрешности η представлены в патенте US No 6,816,808, который включается сюда по ссылке. Для двух хорд А и В, длины которых, равные соответственно L и L, значительно отличаются, индикатор погрешности η рассчитывается:
<mrow><msub><mi>η</mi><mrow><mi>A</mi><mi>B</mi></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>L</mi><mi>B</mi></msub><msub><mi>L</mi><mi>A</mi></msub></mrow><mrow><msub><mi>L</mi><mi>B</mi></msub><mo>−</mo><msub><mi>L</mi><mi>A</mi></msub></mrow></mfrac><mo>⋅</mo><mfrac><mrow><msub><mi>c</mi><mi>B</mi></msub><mo>−</mo><msub><mi>c</mi><mi>A</mi></msub></mrow><mrow><msub><mi>c</mi><mi>B</mi></msub><msub><mi>c</mi><mi>A</mi></msub></mrow></mfrac><mtext> </mtext><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>
A
B
где c - вычисленная скорость звука, проходящего через текучую среду, измеренная вдоль хорды A, c - вычисленная скорость звука, проходящего через текучую среду вдоль хорды В. Аналогичные выражения можно написать для других пар хорд со значительно отличающимися длинами, например, хорд В и D, хорд С и А и хорд С и D, показанных на Фиг.2.
A
Скорость звука c для хорды А определяется:
<mrow><msub><mi>c</mi><mi>A</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>L</mi><mi>A</mi></msub></mrow><mn>2</mn></mfrac><mo>⋅</mo><mfrac><mrow><msub><mi>T</mi><mrow><mn>1,</mn><mi>A</mi></mrow></msub><mo>+</mo><msub><mi>T</mi><mrow><mn>2,</mn><mi>A</mi></mrow></msub></mrow><mrow><msub><mi>T</mi><mrow><mn>1,</mn><mi>A</mi></mrow></msub><msub><mi>T</mi><mrow><mn>2,</mn><mi>A</mi></mrow></msub></mrow></mfrac><mtext> </mtext><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>
1,A
2,A
где Т - время распространения сигнала в направлении против потока вдоль хорды А и Т - время распространения сигнала в направлении потока вдоль хорды А. Скорость звука, относящаяся к хорде В, как и скорости, относящиеся к другим хордам, могут быть выражены аналогичным образом.
1,A
2,A
1,B
2,B
AB
AB
AB
AB
Если скорость распространения звука через текучую среду в расходомере является однородной, и отсутствуют погрешности при измерении значений Т, Т, Т и T времени распространения сигнала, то значение η равно нулю. Однако в случае погрешностей при измерении любого из значений времени распространения сигнала значение η будет отличным от нуля. Как обсуждалось выше, наросты отложения могут вызывать погрешности в измеренных значениях времени распространения сигнала, и, таким образом, вызывать отклонение η от нуля. Поэтому при отслеживании через некоторые промежутки времени значений η может быть обнаружен рост отложений внутри ультразвукового расходомера.
AB
AB
При вычислении η через некоторые промежутки время с целью обнаружения формирования отложений, следует учитывать режим, в котором происходит движение текучей среды. Могут образовываться слои, имеющие различные температуры, то есть, может возникать неоднородность поля температур внутри расходомера, когда скорость потока текучей среды мала или равна нулю. Например, текучая среда вблизи верхней части расходомера может иметь более высокую температуру, чем текучая среда вблизи нижней части расходомера. Так как скорость распространения звука через текучую среду зависит от температуры, скорость звука в объеме текучей среды не будет однородной функцией при неравномерном распределении температур. Различия в скорости прохождения звука через текучую среду вызовут, далее, отличия в значениях η при переходе от одной хорды к следующей.
AB
1,A
2,A
1,B
2,B
Но неоднородность поля температур может быть уменьшена, если текучая среда хорошо перемешивается, например, если при прохождении текучей среды через расходомер возникает турбулентность. Поэтому при вычислении η через некоторые промежутки времени с целью обнаружить сформировавшиеся отложения предпочтительно измерять значения Т, Т, T и Т времени распространения сигнала, когда текучая среда движется через расходомер.
1,i
2,i
1,i
2,i
1,i
2,i
На Фиг.4 показана блок-схема, иллюстрирующая способ обнаружения сформировавшихся наростов отложений внутри ультразвукового расходомера. Алгоритм 400 для этого способа начинает работу с периодического измерения группой электронных устройств значений Т, Т времени распространения сигнала по крайней мере между двумя парами преобразователей, установленных в расходомере (блок 405). В некоторых вариантах значения T, T времени распространения могут измеряться для всех пар преобразователей, а в других вариантах значения Т, T времени распространения могут измеряться лишь для некоторых выбранных пар преобразователей.
1,i
2,i
1,i
2,i
Индикатор погрешности η может быть рассчитан для любых двух пар преобразователей, для которых производятся измерения значений Т, Т времени распространения сигнала, у которых длины хорд значительно отличаются. Пусть, например, имеются четыре пары преобразователей, находящихся внутри ультразвукового расходомера, показанного на Фиг.2. Если значения Т, T времени распространения сигнала измеряются для каждой пары преобразователей, то значение η можно рассчитывать для пар преобразователей, связанных с хордами А и В, В и D, А и С, С и D.
1,i
2,i
1,i,avg
2,i,avg
1,i,avg
2,i,avg
1,i,avg
2,i,avg
i
i
avg
1,i,avg
2,i,avg
В некоторых вариантах значения Т, Т времени распространения сигнала измеряются многократно в течение некоторого промежутка времени, а затем усредняются по этому промежутку времени, и таким образом формируются средние значения T, T (блок 410). Средние значения T, T времени распространения сигнала могут далее использоваться для вычислений скорости звука, при которых в уравнении (7) значения переменных заменяются средними значениями Т, Т времени распространения сигнала и длины хорды L (блок 415). В свою очередь, значения скорости звука c и соответствующие значения длин хорд далее подставляются в уравнение (6) с целью вычисления η (блок 420). В некотором смысле вычисленное значение η  является средним значением η, или η, так как оно определяется с использованием средних значений T, T времени распространения сигнала.
i
1,i
2,i
i
i
avg
Альтернативно, скорость звука c может быть рассчитана при подстановке измеренных значений Т, T времени распространения сигнала и длин L хорд в уравнение (7) (блок 425). В свою очередь, скорость звука c и соответствующие длины хорд могут далее подставляться в уравнение (6) для вычисления η (блок 430). Эти вычисления могут повторяться в течение некоторого промежутка времени, а вычисленные значения η усредняться, формируя, таким образом, среднее значение η, или η (блок 435).
avg
avg
Затем каждое вычисленное значение η может быть оценено с целью определить, обнаруживаются ли внутри ультразвукового расходомера сформировавшиеся наросты отложения (блок 440). Одним из оснований для вывода о том, что происходит рост отложений, является предположение о том, что каждое из вычисленных значений η постепенно изменяется, в общем, на одну и ту же величину и превышает некоторое заранее заданное пороговое значение.
i,1
i,2
avg
Если не обнаружено формирования отложений, то техническое обслуживание в данное время не требуется. Группа электронных устройств, связанная с ультразвуковым расходомером, может продолжать измерять изменения значений T, T времени распространения сигнала и вычислять η для двух или большего числа пар преобразователей без какого-либо вмешательства, пока не будет обнаружен рост отложений. При обнаружении в некоторый момент времени роста отложений может потребоваться техническое обслуживание. В зависимости от ситуации, работа ультразвукового расходомера может быть приостановлена для выполнения технического обслуживания с целью уменьшения или удаления отложений (блок 445). Альтернативно, группа электронных устройств может выполнять дополнительные вычисления, позволяющие определять скорректированное значение объемного расхода текучей среды, проходящей через расходомер, учитывающее рост отложений (блок 450).
Возможность обнаруживать формирование нароста отложений, в соответствии с раскрытыми здесь способами, включая алгоритм 400, позволяет выполнять обслуживание ультразвукового расходомера, когда это необходимо, а не по жесткому графику. Кроме того, возможность количественно определять рост отложений и учитывать его путем корректировки вычисляемого значения объемного расхода потока через расходомер позволяет отложить техническое обслуживание на некоторое время, которое не задается заранее. По крайней мере в некоторых вариантах величина η дает основу для количественного определения роста отложений и корректировки значений объемного расхода текучей среды, проходящей через расходомер.
AB
A
B
A
1,A
1,B
A
A
Вернемся к уравнению (6), где η - функция от скорости звука c и скорости звука c. Предположим, что возникла погрешность ΔT в измеренном значении Т или Т времени распространения сигнала, обусловленная ростом отложений; тогда результирующая погрешность Δc вычисленной скорости звука c равна:
<mrow><mi>Δ</mi><msub><mi>c</mi><mi>A</mi></msub><mo>=</mo><mo>−</mo><mfrac><mrow><msup><mi>c</mi><mn>2</mn></msup></mrow><mrow><msub><mi>L</mi><mi>A</mi></msub></mrow></mfrac><mi>Δ</mi><msub><mi>T</mi><mi>A</mi></msub><mtext> </mtext><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>
A
Таким образом, скорость звука c для хорды А может быть записана в виде:
<mrow><msub><mi>c</mi><mi>A</mi></msub><mo>=</mo><mi>c</mi><mo>+</mo><mi>Δ</mi><msub><mi>c</mi><mi>A</mi></msub><mo>=</mo><mi>c</mi><mo>−</mo><mfrac><mrow><msup><mi>c</mi><mn>2</mn></msup></mrow><mrow><msub><mi>L</mi><mi>A</mi></msub></mrow></mfrac><mi>Δ</mi><msub><mi>T</mi><mi>A</mi></msub><mtext> </mtext><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>
Аналогичные выражения могут быть написаны для других хорд, например хорд В, С и D, показанных на Фиг.2.
A
B
AB
Подставляя выражение (9) и аналогичное выражение для хорды В в уравнение (6) и пренебрегая членами более высокого порядка, представляющими собой произведения ΔT и ΔT, можно представить η таким образом:
<mrow><msub><mi>η</mi><mrow><mi>A</mi><mi>B</mi></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>L</mi><mi>B</mi></msub><mi>Δ</mi><msub><mi>T</mi><mi>A</mi></msub><mo>−</mo><msub><mi>L</mi><mi>A</mi></msub><mi>Δ</mi><msub><mi>T</mi><mi>B</mi></msub></mrow><mrow><msub><mi>L</mi><mi>B</mi></msub><mo>−</mo><msub><mi>L</mi><mi>A</mi></msub></mrow></mfrac><mtext> </mtext><mrow><mo>(</mo><mrow><mn>10</mn></mrow><mo>)</mo></mrow></mrow>
deposit
deposit
Если образование нароста отложений происходит на поверхностях преобразователей, находящихся на противоположных концах хорды, то они отложения вызовут погрешности в значении скорости с распространения звука в текучей среде. Один из способов количественного определения этой погрешности состоит в предположении, что отложения вызывают погрешность в измеренном значении Т времени распространения сигнала. Предположив, что материал отложений характеризуется скоростью звука c, и что толщина нароста отложений на каждой из поверхностей преобразователей равна t, можно определить погрешность ΔТ для времени распространения сигнала, обусловленную отложениями, по формуле:
<mrow><mi>Δ</mi><mi>T</mi><mo>=</mo><mn>2</mn><msub><mi>t</mi><mrow><mi>d</mi><mi>e</mi><mi>p</mi><mi>o</mi><mi>s</mi><mi>i</mi><mi>t</mi></mrow></msub><mrow><mo>(</mo><mrow><mfrac><mn>1</mn><mrow><msub><mi>c</mi><mrow><mi>d</mi><mi>e</mi><mi>p</mi><mi>o</mi><mi>s</mi><mi>i</mi><mi>t</mi></mrow></msub></mrow></mfrac><mo>−</mo><mfrac><mn>1</mn><mi>c</mi></mfrac></mrow><mo>)</mo></mrow><mtext> </mtext><mrow><mo>(</mo><mrow><mn>11</mn></mrow><mo>)</mo></mrow></mrow>
deposit
deposit
A
B
Как показывает уравнение (11), погрешность ΔТ для времени распространения сигнала, обусловленная отложениями, зависит от толщины t нароста отложений, скорости звука с в текучей среде и скорости звука c в материале отложений, но не от длины хорды. Кроме того, если предположить, что рост отложений происходит внутри ультразвукового расходомера равномерно, то погрешность ΔT в значении времени распространения сигнала, обусловленная отложениями, не зависит от хорды. Таким образом, величины ΔT и ΔT в уравнении (10) не зависят от хорд А и В соответственно и могут быть обозначены просто через ΔT. Уравнение (10) тогда упрощается следующим образом:
<mrow><msub><mi>η</mi><mrow><mi>A</mi><mi>B</mi></mrow></msub><mo>=</mo><mi>Δ</mi><mi>T</mi><mtext> </mtext><mrow><mo>(</mo><mrow><mn>12</mn></mrow><mo>)</mo></mrow></mrow>
AB
deposit
Комбинируя уравнения (11) и (12), можно выразить η и погрешность ΔТ в значении времени распространения сигнала, обусловленную отложениями, через толщину t нароста отложений:
<mrow><msub><mi>η</mi><mrow><mi>A</mi><mi>B</mi></mrow></msub><mo>=</mo><mi>Δ</mi><mi>T</mi><mo>=</mo><mn>2</mn><msub><mi>t</mi><mrow><mi>d</mi><mi>e</mi><mi>p</mi><mi>o</mi><mi>s</mi><mi>i</mi><mi>t</mi></mrow></msub><mrow><mo>(</mo><mrow><mfrac><mn>1</mn><mrow><msub><mi>c</mi><mrow><mi>d</mi><mi>e</mi><mi>p</mi><mi>o</mi><mi>s</mi><mi>i</mi><mi>t</mi></mrow></msub></mrow></mfrac><mo>−</mo><mfrac><mn>1</mn><mi>c</mi></mfrac></mrow><mo>)</mo></mrow><mtext> </mtext><mrow><mo>(</mo><mrow><mn>13</mn></mrow><mo>)</mo></mrow></mrow>
deposit
Алгебраически преобразуя уравнение (13), можно выразить толщину t нароста отложений в виде:
<mrow><msub><mi>t</mi><mrow><mi>d</mi><mi>e</mi><mi>p</mi><mi>o</mi><mi>s</mi><mi>i</mi><mi>t</mi></mrow></msub><mo>=</mo><mfrac><mrow><mi>Δ</mi><mi>T</mi></mrow><mn>2</mn></mfrac><mo>⋅</mo><mfrac><mrow><mi>c</mi><msub><mi>c</mi><mrow><mi>d</mi><mi>e</mi><mi>p</mi><mi>o</mi><mi>s</mi><mi>i</mi><mi>t</mi></mrow></msub></mrow><mrow><mi>c</mi><mo>−</mo><msub><mi>c</mi><mrow><mi>d</mi><mi>e</mi><mi>p</mi><mi>o</mi><mi>s</mi><mi>i</mi><mi>t</mi></mrow></msub></mrow></mfrac><mtext> </mtext><mrow><mo>(</mo><mrow><mn>14</mn></mrow><mo>)</mo></mrow></mrow>
AB
1,i
2,i
AB
AB
deposit
Как описано выше со ссылками на Фиг.4, величина η может быть вычислена с использованием измеренных значений Т, Т времени распространения сигнала. Определив η и зная из уравнения (13), что величина η равна погрешности ΔT в значении времени распространения сигнала, можно определить далее толщину t нароста отложений, используя уравнение (14). Таким образом, уравнение (14) позволяет количественно определять рост отложений внутри ультразвукового расходомера.
deposit
deposit
В некоторых вариантах решение о выполнении технического обслуживания ультразвукового расходомера с целью удаления отложений может зависеть от толщины t нароста отложений. Например, может быть выбрано некоторое пороговое значение, такое что обслуживание должно производиться, когда толщина t нароста отложений превысит это значение. Если обслуживание с целью уменьшения или устранения отложений не выполняется, то значение объемного расхода текучей среды через расходомер может быть скорректировано с учетом роста отложений.
Объемный расход текучей среды через расходомер определяется по формуле:
<mrow><mi>Q</mi><mo>=</mo><mi>V</mi><msub><mi>A</mi><mi>c</mi></msub><mtext> </mtext><mrow><mo>(</mo><mrow><mn>15</mn></mrow><mo>)</mo></mrow></mrow>
c
c
где A - площадь поперечного сечения, V - средняя скорость потока. Как обсуждалось выше, формирование отложений ведет к погрешностям как в значении средней скорости V потока, так и в значении площади поперечного сечения A. Таким образом, относительная неопределенность ΔQ/Q объемного расхода может выражаться, как:
<mrow><mfrac><mrow><mi>Δ</mi><mi>Q</mi></mrow><mi>Q</mi></mfrac><mo>=</mo><mfrac><mrow><mi>Δ</mi><mi>V</mi></mrow><mi>V</mi></mfrac><mo>+</mo><mfrac><mrow><mi>Δ</mi><msub><mi>A</mi><mi>c</mi></msub></mrow><mrow><msub><mi>A</mi><mi>c</mi></msub></mrow></mfrac><mtext> </mtext><mrow><mo>(</mo><mrow><mn>16</mn></mrow><mo>)</mo></mrow></mrow>
c
c
c
c
c
Здесь ΔV/V и ΔA/A - относительные неопределенности значений средней скорости потока V и площади поперечного сечения A соответственно. Чтобы определить относительную неопределенность объемного расхода ΔQ/Q и, таким образом, погрешность ΔQ в значении объемного расхода Q, обусловленную отложениями, оцениваются относительные неопределенности средней скорости потока и площади поперечного сечения, то есть, соответственно, ΔV/V, ΔA/A.
c
Площадь поперечного сечения A:
<mrow><msub><mi>A</mi><mi>c</mi></msub><mo>=</mo><mfrac><mrow><mi>π</mi><msup><mi>D</mi><mn>2</mn></msup></mrow><mn>4</mn></mfrac><mtext> </mtext><mrow><mo>(</mo><mrow><mn>17</mn></mrow><mo>)</mo></mrow></mrow>
c
c
где D - внутренний диаметр муфты расходомера, измеренный при отсутствии отложений. Формирование отложений уменьшает эффективный внутренний диаметр муфты расходомера. Таким образом, относительная неопределенность ΔA/A площади поперечного сечения может быть выражена так:
<mrow><mfrac><mrow><mi>Δ</mi><msub><mi>A</mi><mi>c</mi></msub></mrow><mrow><msub><mi>A</mi><mi>c</mi></msub></mrow></mfrac><mo>=</mo><mfrac><mn>1</mn><mrow><msub><mi>A</mi><mi>c</mi></msub></mrow></mfrac><mo>⋅</mo><mfrac><mrow><mo>∂</mo><msub><mi>A</mi><mi>c</mi></msub></mrow><mrow><mo>∂</mo><mi>D</mi></mrow></mfrac><mi>Δ</mi><mi>D</mi><mtext> </mtext><mrow><mo>(</mo><mrow><mn>18</mn></mrow><mo>)</mo></mrow></mrow>
c
deposit
c
c
Выражая A через D с помощью уравнения (17) и замечая, что ΔD равно -2t , получим, что относительная неопределенность ΔA/A площади поперечного сечения равна:
<mrow><mfrac><mrow><mi>Δ</mi><msub><mi>A</mi><mi>c</mi></msub></mrow><mrow><msub><mi>A</mi><mi>c</mi></msub></mrow></mfrac><mo>=</mo><mo>−</mo><mfrac><mrow><mn>4</mn><msub><mi>t</mi><mrow><mi>d</mi><mi>e</mi><mi>p</mi><mi>o</mi><mi>s</mi><mi>i</mi><mi>t</mi></mrow></msub></mrow><mi>D</mi></mfrac><mtext> </mtext><mrow><mo>(</mo><mrow><mn>19</mn></mrow><mo>)</mo></mrow></mrow>
Средняя скорость V потока задается выражением:
<mrow><mi>V</mi><mo>=</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder><mrow><msub><mi>W</mi><mi>i</mi></msub><msub><mi>V</mi><mi>i</mi></msub></mrow></mstyle><mtext> </mtext><mrow><mo>(</mo><mrow><mn>20</mn></mrow><mo>)</mo></mrow></mrow>
i
где W - весовой множитель для хорды i. Таким образом, относительная неопределенность ΔV/V средней скорости потока может быть выражена в виде:
<mrow><mfrac><mrow><mi>Δ</mi><mi>V</mi></mrow><mi>V</mi></mfrac><mo>=</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder><mrow><msub><mi>W</mi><mi>i</mi></msub><mfrac><mrow><mi>Δ</mi><msub><mi>V</mi><mi>i</mi></msub></mrow><mrow><msub><mi>V</mi><mi>i</mi></msub></mrow></mfrac><mtext> </mtext><mrow><mo>(</mo><mrow><mn>21</mn></mrow><mo>)</mo></mrow></mrow></mstyle></mrow>
i
Средняя скорость V потока вдоль хорды i задается выражением:
<mrow><msub><mi>V</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msubsup><mi>L</mi><mi>i</mi><mn>2</mn></msubsup></mrow><mrow><mn>2</mn><msub><mi>X</mi><mi>i</mi></msub></mrow></mfrac><mo>⋅</mo><mfrac><mrow><msub><mi>T</mi><mrow><mn>1,</mn><mi>i</mi></mrow></msub><mo>−</mo><msub><mi>T</mi><mrow><mn>2,</mn><mi>i</mi></mrow></msub></mrow><mrow><msub><mi>T</mi><mrow><mn>1,</mn><mi>i</mi></mrow></msub><msub><mi>T</mi><mrow><mn>2,</mn><mi>i</mi></mrow></msub></mrow></mfrac><mtext> </mtext><mrow><mo>(</mo><mrow><mn>22</mn></mrow><mo>)</mo></mrow></mrow>
i
i
i
1,i
2,i
i
где L - длина хорды i, X - продольная составляющая L, Т - время распространения сигнала против потока (то есть, значение времени распространения акустического сигнала вдоль хорды i от преобразователя, находящегося ниже по потоку, к преобразователю, находящемуся выше по потоку), и T - время распространения сигнала в направлении потока. Таким образом, относительная неопределенность ΔV/V скорости потока для хорды i может быть выражена следующим образом:
<mrow><mfrac><mrow><mi>Δ</mi><msub><mi>V</mi><mi>i</mi></msub></mrow><mrow><msub><mi>V</mi><mi>i</mi></msub></mrow></mfrac><mo>=</mo><mfrac><mn>1</mn><mrow><msub><mi>V</mi><mi>i</mi></msub></mrow></mfrac><mrow><mo>(</mo><mrow><mfrac><mrow><mo>∂</mo><msub><mi>V</mi><mi>i</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>X</mi><mi>i</mi></msub></mrow></mfrac><mi>Δ</mi><msub><mi>X</mi><mi>i</mi></msub><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>V</mi><mi>i</mi></msub></mrow><mrow><msub><mi>T</mi><mrow><mn>1,</mn><mi>i</mi></mrow></msub></mrow></mfrac><mi>Δ</mi><mi>T</mi><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>V</mi><mi>i</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>T</mi><mrow><mn>2,</mn><mi>i</mi></mrow></msub></mrow></mfrac><mi>Δ</mi><mi>T</mi></mrow><mo>)</mo></mrow><mtext> </mtext><mrow><mo>(</mo><mrow><mn>23</mn></mrow><mo>)</mo></mrow></mrow>
1,i
2,i
где погрешности в значениях Т, Т времени распространения сигнала равны ДТ, как обсуждалось выше. С помощью уравнения (22) данное выше выражение может быть упрощено:
<mrow><mfrac><mrow><mi>Δ</mi><msub><mi>V</mi><mi>i</mi></msub></mrow><mrow><msub><mi>V</mi><mi>i</mi></msub></mrow></mfrac><mo>=</mo><mfrac><mrow><mi>Δ</mi><msub><mi>X</mi><mi>i</mi></msub></mrow><mrow><msub><mi>X</mi><mi>i</mi></msub></mrow></mfrac><mo>−</mo><mfrac><mrow><mn>2</mn><mi>c</mi></mrow><mrow><msub><mi>L</mi><mi>i</mi></msub></mrow></mfrac><mi>Δ</mi><mi>T</mi><mtext> </mtext><mrow><mo>(</mo><mrow><mn>24</mn></mrow><mo>)</mo></mrow></mrow>
где с - скорость звука для текучей среды.
i
Продольная составляющая X дается выражением:
<mrow><msub><mi>X</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><mi>D</mi><msqrt><mrow><mn>1</mn><mo>−</mo><msubsup><mi>f</mi><mi>i</mi><mn>2</mn></msubsup></mrow></msqrt></mrow><mrow><mi>tan</mi><mi>θ</mi></mrow></mfrac><mtext> </mtext><mrow><mo>(</mo><mrow><mn>25</mn></mrow><mo>)</mo></mrow></mrow>
i
i
i
i
i
i
i
где θ - угол, образуемый хордой i с направлением потока и f - относительное расстояние, на которое измерительная плоскость хорды i удалена от центральной плоскости расходомера. На Фиг.5 и 6 схематично показаны относительные расстояния f и протяженность X вдоль оси для хорды i, имеющей длину L. Рост отложений влечет к погрешностям в значениях среднего диаметра расходомера и относительного расстояния f, на котором измерительная плоскость хорды i находится от центральной плоскости расходомера. Таким образом, относительная неопределенность ΔX/X в протяженности хорды вдоль оси может быть выражена в виде:
<mrow><mfrac><mrow><mi>Δ</mi><msub><mi>X</mi><mi>i</mi></msub></mrow><mrow><msub><mi>X</mi><mi>i</mi></msub></mrow></mfrac><mo>=</mo><mfrac><mn>1</mn><mrow><msub><mi>X</mi><mi>i</mi></msub></mrow></mfrac><mrow><mo>(</mo><mrow><mfrac><mrow><mo>∂</mo><msub><mi>X</mi><mi>i</mi></msub></mrow><mrow><mo>∂</mo><mi>D</mi></mrow></mfrac><mi>Δ</mi><mi>D</mi><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>X</mi><mi>i</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>f</mi><mi>i</mi></msub></mrow></mfrac><mi>Δ</mi><msub><mi>f</mi><mi>i</mi></msub></mrow><mo>)</mo></mrow><mtext> </mtext><mrow><mo>(</mo><mrow><mn>26</mn></mrow><mo>)</mo></mrow></mrow>
deposit
Как и выше, ΔD можно выразить через толщину t нароста отложений:
<mrow><mi>Δ</mi><mi>D</mi><mo>=</mo><mo>−</mo><mn>2</mn><msub><mi>t</mi><mrow><mi>d</mi><mi>e</mi><mi>p</mi><mi>o</mi><mi>s</mi><mi>i</mi><mi>t</mi></mrow></msub><mtext> </mtext><mrow><mo>(</mo><mrow><mn>27</mn></mrow><mo>)</mo></mrow></mrow>
i
i
Кроме того, погрешность в значении относительного расстояния Δf/f равна:
<mrow><mfrac><mrow><mi>Δ</mi><msub><mi>f</mi><mi>i</mi></msub></mrow><mrow><msub><mi>f</mi><mi>i</mi></msub></mrow></mfrac><mo>=</mo><mfrac><mrow><mn>2</mn><msub><mi>t</mi><mrow><mi>d</mi><mi>e</mi><mi>p</mi><mi>o</mi><mi>s</mi><mi>i</mi><mi>t</mi></mrow></msub></mrow><mi>D</mi></mfrac><mtext> </mtext><mrow><mo>(</mo><mrow><mn>28</mn></mrow><mo>)</mo></mrow></mrow>
i
i
i
i
Подставляя выражения для ΔD и Δf/f из уравнений (27) и (28) в уравнение (26), получаем относительную неопределенность ΔX/X в протяженности хорды вдоль оси в таком виде:
<mrow><mfrac><mrow><mi>Δ</mi><msub><mi>X</mi><mi>i</mi></msub></mrow><mrow><msub><mi>X</mi><mi>i</mi></msub></mrow></mfrac><mo>=</mo><mo>−</mo><mfrac><mrow><mn>2</mn><msub><mi>t</mi><mrow><mi>d</mi><mi>e</mi><mi>p</mi><mi>o</mi><mi>s</mi><mi>i</mi><mi>t</mi></mrow></msub></mrow><mrow><mi>D</mi><mrow><mo>(</mo><mrow><mn>1</mn><mo>−</mo><msubsup><mi>f</mi><mi>i</mi><mn>2</mn></msubsup></mrow><mo>)</mo></mrow></mrow></mfrac><mtext> </mtext><mrow><mo>(</mo><mrow><mn>29</mn></mrow><mo>)</mo></mrow></mrow>
i
i
i
i
Заменяя далее ΔX/X в уравнении (24) на выражение, взятое из уравнения (29), можно упростить выражение для относительной неопределенности ΔV/V скорости потока для хорды i:
<mrow><mfrac><mrow><mi>Δ</mi><msub><mi>V</mi><mi>i</mi></msub></mrow><mrow><msub><mi>V</mi><mi>i</mi></msub></mrow></mfrac><mo>=</mo><mn>2</mn><msub><mi>t</mi><mrow><mi>d</mi><mi>e</mi><mi>p</mi><mi>o</mi><mi>s</mi><mi>i</mi><mi>t</mi></mrow></msub><mrow><mo>(</mo><mrow><mfrac><mn>1</mn><mrow><mi>D</mi><mrow><mo>(</mo><mrow><mn>1</mn><mo>−</mo><msubsup><mi>f</mi><mi>i</mi><mn>2</mn></msubsup></mrow><mo>)</mo></mrow></mrow></mfrac><mo>−</mo><mfrac><mn>2</mn><mrow><msub><mi>L</mi><mi>i</mi></msub></mrow></mfrac><mo>⋅</mo><mfrac><mrow><mi>c</mi><mo>−</mo><msub><mi>c</mi><mrow><mi>d</mi><mi>e</mi><mi>p</mi><mi>o</mi><mi>s</mi><mi>i</mi><mi>t</mi></mrow></msub></mrow><mrow><msub><mi>c</mi><mrow><mi>d</mi><mi>e</mi><mi>p</mi><mi>o</mi><mi>s</mi><mi>i</mi><mi>t</mi></mrow></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mtext> </mtext><mrow><mo>(</mo><mrow><mn>30</mn></mrow><mo>)</mo></mrow></mrow>
AB
1,i
2,i
deposit
AB
deposit
c
c
i
i
Как было описано выше со ссылками на Фиг.4, величина η может быть вычислена с использованием измеренных значений Т, Т времени распространения сигнала. Используя уравнение (13), можно определить толщину t нароста отложений как функцию от η. Толщина отложений t используется как известная величина ("вход") в уравнениях (19) и (30) для определения ΔA/A и ΔV/V соответственно. Наконец, относительная погрешность ΔQ/Q в значении объемного расхода может быть оценена при помощи уравнения (16).
c
С учетом сформировавшихся отложений скорректированное значение объемного расхода Q текучей среды, проходящей через расходомер, равно:
<mrow><msub><mi>Q</mi><mi>c</mi></msub><mo>=</mo><mi>Q</mi><mo>+</mo><mi>Δ</mi><mi>Q</mi><mtext> </mtext><mrow><mo>(</mo><mrow><mn>31</mn></mrow><mo>)</mo></mrow></mrow>
где Q - объемный расход текучей среды, проходящей через расходомер, при отсутствии отложений, величина которого может быть определена с использованием уравнения (15).
avg
На Фиг.7 приведена блок-схема, иллюстрирующая алгоритм 500 данного способа, служащий для количественного определения нароста отложений внутри ультразвукового расходомера и корректировки значений объемного расхода текучей среды, проходящей через расходомер, с учетом отложений. Алгоритм 500 начинает работу с выявления и отслеживания нароста отложений, выполняемых в соответствии с алгоритмом 400, который обсуждался выше со ссылками на Фиг.4 (блок 505). Одним из результатов алгоритма 400 является значение η которое может быть далее подставлено в уравнение (13), что позволяет вычислять толщину нароста отложений (блок 510).
c
deposit
i
i
i
i
i
i
c
c
Далее в алгоритме 500 выполняются шаги для вычисления скорректированного значения Q объемного расхода текучей среды, проходящей через расходомер, с использованием значения t толщины нароста отложений. Относительная неопределенность ΔV/V площади поперечного сечения определяется с помощью уравнения (19) (блок 515). Относительная неопределенность ΔV/V средней измеряемой вдоль хорды скорости потока для каждой хорды i вычисляется с использованием уравнения (30) (блок 520). Вычисленные значения ΔV/V для всех хорд i затем подставляются в уравнение (21), чтобы вычислить относительную неопределенность ΔV/V средней скорости потока (блок 525). Относительная неопределенность значения ΔQ/Q объемного расхода определяется по вычисленным значениям ΔA/A и ΔV/V с помощью уравнения (16) (блок 530).
c
Объемный расход О текучей среды, проходящей через расходомер при отсутствии отложений, может быть определен с помощью уравнения (15) (блок 535). Вычисленные значения ΔQ/Q могут быть далее скомбинированы так, чтобы получить величину ΔQ погрешности в значении объемного расхода, обусловленной ростом отложений (блок 540). Наконец, скорректированное значение объемного расхода Q текучей среды, проходящей через расходомер, может быть определено с использованием уравнения (31) (блок 545).
Специалист в данной области может, исходя из представленного здесь описания, комбинировать описанные способы для создания программного обеспечения, которое, в сочетании с соответствующими универсальными или специальными аппаратными средствами, может использоваться для создания компьютерной системы и/или компьютерных блоков, реализующих данное изобретение, для создания компьютерной системы и/или компьютерных блоков, выполняющих действия, предлагаемые в данном изобретении, и/или для создания машиночитаемых носителей информации с хранящейся на них информацией, связанных с процессором, реализующим некоторые из аспектов способа, предлагаемого в изобретении. При использовании указанного машиночитаемого носителя информации описанные здесь способы могут применяться периодически или же по мере надобности.
Приведенное выше обсуждение, как предполагается, лишь иллюстрирует принципы и различные варианты данного изобретения. Многочисленные варианты и модификации станут очевидными специалистам в данной области после оценки в полной мере раскрытия изобретения. Предполагается, что следующая далее формула изобретения будет интерпретироваться как охватывающая все такие варианты и модификации.