Ультразвуковая диагностика восходящих газожидкостных потоков с использованием распределенного электромеханического преобразователя

Владимиров, Илья Александрович

Ультразвуковая диагностика восходящих газожидкостных потоков с использованием распределенного электромеханического преобразователя тема автореферата и диссертации по физике, 01.04.06 ВАК РФ

Владимиров, Илья Александрович АВТОР

кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ

Нижний Новгород МЕСТО ЗАЩИТЫ

2012 ГОД ЗАЩИТЫ

01.04.06 КОД ВАК РФ

Диссертация по физике на тему «Ультразвуковая диагностика восходящих газожидкостных потоков с использованием распределенного электромеханического преобразователя»

Автореферат диссертации на тему "Ультразвуковая диагностика восходящих газожидкостных потоков с использованием распределенного электромеханического преобразователя"

На правах рукописи

Владимиров Илья Александрович

УЛЬТРАЗВУКОВАЯ ДИАГНОСТИКА ВОСХОДЯЩИХ ГАЗОЖИДКОСТНЫХ ПОТОКОВ С ИСПОЛЬЗОВАНИЕМ РАСПРЕДЕЛЕННОГО ЭЛЕКТРОМЕХАНИЧЕСКОГО ПРЕОБРАЗОВАТЕЛЯ

01.04.06 - Акустика

1 7 МЛ Г, 2012

Автореферат диссертации на соискание ученой степени кандидата физико-математических наук

Нижний Новгород - 2012

005044316

Работа выполнена в федеральном государственном бюджетном образовательном учреждении высшего профессионального образования "Нижегородский государственный университет им. Н.И. Лобачевского"

Научный руководитель: кандидат физико-математических наук, доцент

Муякшин Сергей Иванович

Официальные оппоненты:

доктор технических наук Салин Борис Михайлович ИПФ РАН

кандидат физико-математических наук, доцент Мартьянов Александр Иванович ННГУ им. Н.И. Лобачевского

Ведущая организация: Федеральное государственное научное учреждение

"Научно-исследовательский радиофизический

институт"

Защита состоится " 30 " мая 2012 г. в 15:00 на заседании диссертационного совета Д 212.166.07 при Нижегородском государственном университете им. Н.И. Лобачевского по адресу: 603950, Нижний Новгород, пр. Гагарина, 23, корп. 4, ауд. 203.

С диссертацией можно ознакомиться в фундаментальной библиотеке Нижегородского государственного университета им. Н.И. Лобачевского.

Автореферат разослан "28" ¿Ь \Л2012 г

Ученый секретарь диссертационного совета к.ф.-м.н., доцент I В.В. Черепенников

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность исследования. Двухфазные потоки в виде смеси жидкости и газа встречаются в разных областях техники: в добыче нефти и газа, в химической и пищевой промышленности, в энергетике, в коммунальном хозяйстве. Для учета добытого углеводородного сырья, управления технологическими установками, контроля стоков, необходимо измерение расхода жидкой и газовой фаз по отдельности. В нефтяной промышленности проблема измерения расходов таких потоков возникла еще в начале 1980-х годов. Необходимо, чтобы приборы, измеряющие расход многофазного потока (расходомеры) удовлетворяли следующим требованиям: обладали приемлемой точностью и надежностью, не создавали возмущений в потоке, не зависели от режима течения и были пригодными для применения во всем возможном диапазоне изменения расходов компонентов потока. Несмотря на большое количество технических решений [1], предложенных в последние годы, ни один из двухфазных расходомеров, появившихся на рынке, всем указанным требованиям все же не отвечает [2]. В настоящее время для определения расходов потока в основном используется предварительная дегазация потока [2], Такие установки дороги и сложны в эксплуатации. Поэтому создание надежных, точных и не слишком дорогих многофункциональных приборов для измерения расхода компонентов потока без их предварительного разделения является актуальной задачей. Отметим, что наиболее перспективными методами диагностики двухфазных потоков, являются активные акустические методы. В представленной работе для изучения и определения параметров двухфазных потоков используется ультразвуковая система, основанная на методе измерения скорости неоднородностей с помощью эффекта Доплера. Принцип работы системы приводиться в [3,4]. В работе [3] описан способ ввода ультразвукового пучка в поток, основанный на использовании эффекта трансформации изгибной волны в стенке трубы (волны Лэмба) в продольную звуковую волну в жидкости. Там же предложено использовать для формирования волны Лэмба накладную фазированную антенную решетку из пьезоэлементов. В патенте [4] описан способ измерения параметров двухфазного потока и система для его реализации, в которой используется описанный выше способ ввода ультразвука в поток и эффект Доплера. Однако в [4] раскрывается только лишь принципиальная возможность определения расходов фаз потока. Поэтому полная реализация возможностей системы, включая достижение в большом диапазоне изменения расходов двухфазного потока достаточной точности измерения, требует систематического и углубленного исследования процессов распространения и рассеяния ультразвука в газожидкостных смесях (ГЖС).

преобразователя, выполненного в виде фазированной антенной решетки (ФАР), и изучение возможности использования этих закономерностей для раздельного определения расходов газовой и жидкой фаз.

Для достижения этой цели в работе были поставлены следующие основные задачи:

1. Создание и испытание лабораторной экспериментальной установки и программного обеспечения для изучения процессов распространения и рассеяния в ГЖС ультразвуковых волн, излучаемых накладными ФАР.

2. Теоретические и экспериментальные исследования процессов распространения и рассеяния в ГЖС ультразвуковых волн, излучаемых накладными ФАР, при различных режимах течения.

3. Разработка на основе проведенных теоретических и экспериментальных исследований, способа определения расхода фаз газожидкостного потока в «пузырьковом» и «снарядном» режимах течения.

Положения, выносимые на защиту:

1. Использование разработанной модели формирования зависимости энергии рассеянных ультразвуковых сигналов от дальности (профилей), в «пузырьковом» режиме газожидкостного потока, позволяет определять газосодержание по форме профиля энергии.

2. В «снарядном» режиме течения, при прохождении мимо датчика жидкой перемычки, рассеяние происходит на содержащихся в ней пузырьках газа, а при прохождении газовой полости - на неоднородностях толщины стекающей пленки. Неоднородности толщины пленки, на которых происходит рассеяние, движутся относительно стенки со скоростью близкой к скорости стекания пленки.

3. Для определения расходов фаз потока, в «пузырьковом» режиме течения, наиболее информативными являются зависимости от дальности энергии и средней частоты рассеянного в потоке сигнала, а в «снарядном» режиме - временные реализации средней частоты сигнала, получаемые с различных удалений от датчика.

Научная новизна работы заключается в следующем:

1. Объяснен механизм формирования зависимостей энергии рассеянного сигнала от дальности в «пузырьковом» потоке. На его основе предложен способ измерения газосодержания (при фиксированном составе газовой и жидкой фаз), преимуществом которого является использование одного приемно-излучающего датчика.

2. При «снарядном» режиме течения исследованы механизм рассеяния ультразвуковых волн Лэмба в стенке трубы, по которой стекает пленка жидкости, а также поведение зависимостей различных характеристик рассеянного сигнала от времени (энергии, доплеровского сдвига частоты, ширины спектра).

3. Экспериментально проверена модель расчета интегральных расходных характеристик газожидкостного потока с использованием ультразвуковой системы с накладными ФАР.

Достоверность результатов. Научные положения и выводы, представленные в диссертационной работе, аргументированы, обоснованы и базируются на фундаментальных теоретических положениях общей акустики и теории распространения и рассеяния волн, а также результатах целенаправленных экспериментальных исследований.

Достоверность полученных в диссертационной работе результатов обосновывается использованием современных средств и методик проведения исследований, повторяемостью результатов многократных экспериментов и согласованностью экспериментальных и теоретических оценок.

Практическое значение. Полученные в работе теоретические и экспериментальные результаты могут быть положены в основу разрабатываемых ультразвуковых измерителей расходов двухфазных потоков, а созданная и опробованная экспериментальная установка может рассматриваться как прототип ультразвукового расходомера.

Апробация работы. Основные результаты работы докладывались на XII Всероссийском семинаре «Волновые явления в неоднородных средах», (Звенигород, МГУ, 2010), XII, XIII и XIV научных конференциях по радиофизике (Нижний Новгород, ННГУ, 2008-2010), XV и XVI Нижегородских сессиях молодых ученых (Нижний Новгород, 2009,2010).

По теме диссертационной работы опубликовано 5 статей в трудах конференций и 3 статьи в журналах из перечня, рекомендованного ВАК для публикации результатов кандидатских диссертаций. Работа выполнена при поддержке проекта ФЦП "Научные и научно-педагогические кадры инновационной России" (контракт № П2134), руководителем которого был автор диссертации, и гранта государственной поддержки ведущих научных школ НШ-3700.2010.2.

Личный вклад автора. На протяжении всей исследовательской работы диссертант принимал непосредственное участие в постановке задач, проведении экспериментов и расчетов, построении моделей, обсуждении и физической интерпретации результатов исследований. Исключением является раздел 3.1.3, в котором автору принадлежат только проведение экспериментов и обсуждение результатов.

Структура и объем работы. Диссертация состоит из введения, четырех глав, заключения и списка цитируемой литературы. Общий объем диссертации составляет 112 страниц, включая 79 рисунков, 3 таблицы и список литературы из 73 наименований.

СОДЕРЖАНИЕ РАБОТЫ

Во введении сформулирована проблема исследования, обоснована ее актуальность, определены цели работы и задачи исследований. Сформулированы основные научные результаты, выносимые на защиту, определена их научная новизна и практическая значимость. Приведены сведения по апробации результатов работы.

Первая глава посвящена общему описанию и классификации двухфазных потоков. Приводятся количественные характеристики двухфазных потоков, а так же их возможные режимы течения, формируемые в вертикальной трубе в зависимости от различных соотношений расходов газа и жидкости. Также представлен обзор методов измерения параметров газожидкостных потоков, основой которого является проведенный в работе патентный анализ. Результаты проведенного анализа позволили определить направления исследований и сформулировать основные задачи, решение которых обусловило достижение цели данной диссертационной работы.

Вторая глава посвящена описанию используемой в данной работе ультразвуковой измерительной системы. В системе применяется бесконтактный способ ввода ультразвука в поток, основанный на трансформации антисимметричной волны Лэмба, возбуждаемой в стенке трубы с помощью накладной ФАР, параметры которой для этого подобраны специальным образом, в продольную звуковую волну в жидкости. При этом в жидкости одновременно формируются два симметричных луча. Рис. 1 иллюстрирует этот способ (показан только один луч) в случае, когда в потоке содержится небольшое количество неоднородностей, например, в виде пузырьков газа малых размеров.

Датчик (1) возбуждает в стенке трубы импульс изгибной волны (2), бегущий (как вверх, так и вниз) со скоростью Су. Эксперимент показывает, что изгибная волна распространяется по стенке трубопровода на расстояние, в несколько раз превышающее длину преобразователя. Эта волна будет излучать в жидкость звуковую волну (3,4) (ее скорость Сх), фронт которой, образует со стенкой трубы угол 0=агся'и(С., /Су). Таким образом, датчик и возбужденный участок трубы играют роль распределенного электромеханического преобразователя. На приеме за счет специальных приемов обработки данных формируются два независимых наклонных луча. В зависимости от задержек приема < соответствующий импульсный объем (или иначе «элемент дальности») будет перекрывать сечение трубы либо частично (3), либо полностью (4). Доплеровское смещение частоты (ДСЧ) сигнала, рассеянного одиночной

Рис. 1 Бесконтактный ввод ультразвука в поток. Белые стрелки-путь звука в стенке, черные-в жидкости.

неоднородностью, определяется соотношением: А/=2-/и-У-в1п(0)/С5 (V - скорость рассеивателя относительно антенны, /0 - частота заполнения импульса). Поскольку и'л(в)=С/С/, то доплеровский сдвиг А/ =2-/1,-У/ Cf не зависит от скорости звука в среде.

Для выделения доплеровского сдвига частоты в системе используется широко распространенный в радио- и акустической локации когерентно-импульсный метод обработки эхосигналов [8]. Сигналы, относящиеся к определенным элементам дальности, подвергались спектральному анализу со скользящим окном. Перед вычислением спектра из анализируемого участка сигнала вычиталось среднее значение. Для спектров по нижепредставленным формулам вычислялись их интегральные характеристики: энергия Еп, средняя частота спектра и ширина АР.

Здесь: 5„- отсчеты спектра; Sf=FJN - частотное разрешение, Рг - частота квантования (в данном случае - частота повторения импульсов), .¡V - число отсчетов в спектральном окне. Экспериментально измеряя зависимости интегральных характеристик от номера «элемента дальности» или от времени (при выбранном «элементе дальности»), мы получим соответственно профили этих характеристик или их временные зависимости.

Измерительная система имела следующие параметры. Рабочая частота ФАР /</=962 КГц, расстояние между ее элементами /=0,78 мм, число элементов - М-16, их ширина - 8 мм. Использовалось импульсное излучение с частотой повторения /,'г=4 кГц, длительность импульса составляла 13 периодов, т.е. 13,5 мкс. Частота квантования /^сигналов составила 189 КГц.

Для правильной интерпретации результатов измерений характеристик потока необходимо знать структуру акустического поля в трубе и расположение импульсных объемов, соответствующих определенным временным стробам. В результате исследования акустического поля внутри трубы, заполненной жидкостью, экспериментально определено, что в ней существует пучек излучаемого сигнала, ширина которого примерно равна ширине излучающего датчика. Выяснено также, что зондирующий импульс несколько раз отражается от противоположных стенок трубы.

Третья глава посвящена теоретическому и экспериментальному исследованию процессов рассеяния ультразвуковых волн в «пузырьковом» потоке, а также в стекающей пленке жидкости, которая формируется при «снарядном» течении газожидкостной смеси. Все представленные в главе экспериментальные результаты получены на специально созданной в ИНГУ лабораторной установке.

Параграф 3.1 посвящен исследованию «пузырькового» потока.

4^ =

В разделе 3.1.1 представлена упрощенная модель формирования профилей энергии рассеянного сигнала, полученная с использованием первого приближения теории многократного рассеяния. Схема импульсной локации моностатическая, излучающий датчик работает и на прием. В модели использованы результаты экспериментального исследования акустического поля в трубе.

Основной результат, вытекающий из модели, заключается в том, что профиль энергии рассеянного сигнала имеет несколько максимумов (рис.2), и фактически представляет собой наложение (сумму) двух различных профилей. Первый формируется при обычном рассеянии «назад». Ему соответствует первый максимум.

сигнал проходят по одному и тому же пути. Этому случаю соответствует импульсный объем с протяженностью СДУ2 (рис.3), где ^-длительность излученного

импульса. Второй профиль рассеянного сигнала формируется иначе. В этом случае сигнал доходит до рассеивателей, например, после отражения от верхней стенки. На датчик же попадает сигнал, рассеянный уже не назад (по отношению к падающему на рассеиватель сигналу), а в переднюю полусферу («вперед»). Соответствующий импульсный объем /1У2(0 имеет уже другую ориентацию и продольные размеры, которые для угла 0=28° (случай, реализуемый в эксперименте) примерно вдвое больше ЛУ^) (рис.3). Другая форма второго импульсного объема АУ2(0 связана с геометрией задачи, вследствие которой время задержки рассеянного сигнала не зависит от координаты поперек трубы (вдоль диаметра). Поэтому, если в трубе есть несколько рассеивателей, находящихся на одной линии вдоль диаметра, сигналы, рассеянные от них «вперед», на приемник будут приходить одновременно.

В разделе 3.1.2 представлено экспериментальное исследование рассеяния ультразвуковых волн в потоке пузырьков и проведено сравнение эксперимента с теорией. На рис.4 показаны зависимости энергии от номера элемента дальности, полученные в опытах с водой при увеличении газосодержания в трубе диаметром 40

мм. Эти зависимости, действительно, имеют несколько максимумов, что согласуется с модельными расчетами.

На рис. 5 представлено отношение величин первого максимума энергии ко второму в зависимости от относительного газосодержания для пузырькового потока в

воде и гидросиле, полученное в эксперименте. Рядом, на рис.6 представлена зависимость той же величины для воды, рассчитанная с помощью модели. Эта зависимость имеет практически линейный рост, при этом поведение теоретической и экспериментальной зависимостей отношения первого и второго максимумов качественно совпадает. Наблюдаемое на графике аномальное поведение зависимости для гидросила при газосодержании более 3%, вероятно, объясняется

использованием для интерпретации результатов приближенной модели рассеяния.

В-да

Вомф ЭД

-П'О.ОО?

.....6=0 025

----(1=0 056

.....П«0 092

Рис.4 Профили энергии для воды. /2-газосодержание

®«> Вт '''' ■ Гвдроецт

"" .......:......■

2 5 1 Галосли.»] АХНИ? %

Рис.5 Отношение первого и второго максимумов энергии в зависимости от относительного газосодержания

Рис.6 Отношение первого максимумов энергии. Модель.

и второго

В разделе 3.1.3 экспериментально исследован метод определения размеров мелких газовых пузырьков в достаточно разреженном потоке с помощью набора полосовых фильтров настроенных на различные частоты. Эта задача достаточно важна, поскольку размер пузырьков главным образом определяет их рассеивающие свойства, и знание этого параметра позволяет строить для «пузырькового» потока ясную модель рассеяния, которая, в свою очередь, является основой для определения содержания газа в потоке.

Для получения информации о размерах пузырька, производилось численное моделирование процесса локации на основе решения уравнения Рэлея - Плессета для

радиальных колебаний пузырька [6], а затем частотная фильтрация принимаемого сигнала, содержащего собственные колебания пузырька, с помощью набора полосовых фильтров, настроенных на различные частоты. Методика разделения сигналов от пузырьков разного размера была проверена экспериментально. Для возбуждения собственных колебаний пузырьков в широком диапазоне создавался широкополосный зондирующий импульс. При этом регистрация сигналов, порождаемых колебаниями пузырька, производилась в полосе ниже центральной частоты зондирования. Генератор радиоимпульсов вырабатывал импульсные сигналы с длительностью около 2 мкс с частотой заполнения 1 МГц. Частотный диапазон, в котором осуществлялось возбуждение и регистрация собственных колебаний пузырька, составлял от нескольких килогерц до 750 кГц. В результате экспериментов была продемонстрирована возможность выделения сигналов от пузырьков с оценкой их размеров в широком диапазоне радиусов. Недостатком такого метода является то, что сигналы от пузырьков, необходимые для их идентификации, должны приходить достаточно редко, чтобы на осциллограмме существовала возможность разделения времен их приходов. Поэтому данный метод измерения работоспособен только в разреженных потоках пузырьков.

Параграф 3.2 посвящен исследованию процессов

рассеяния ультразвуковых

сигналов в стекающей пленке жидкости, возникающей при «снарядном» режиме течения. Ниже представлены результаты эксперимента при всплытии двух, следующих друг за другом пузырей. Первый пузырь всплывает в невозмущенной жидкости. На рис.7 показаны временные зависимости

доплеровского сдвига частоты (ДСЧ) и энергии рассеянного сигнала. Обе реализации относятся к лучу (и волне Лэмба), направленному вверх - т.е. в сторону движения пузырей.

Основной интерес представляет тот факт, что энергия рассеянного на пленке сигнала в момент всплытия второго пузыря превышает энергию, рассеянную на пленке, стекающей в момент прохождения первого пузыря. Первый пузырь, который всплывает в невозмущенной жидкости, имеет правильную, симметричную форму. При его движении деформации поверхности стекающей пленки незначительны. Пузырь, следующий за первым, имеет уже искаженную, не симметричную форму, а стекающая пленка имеет значительные деформации толщины. Вероятно, именно за

. 1Е+12

: '£+11 лг»::

¡«ш-

А II 1 А' Ц ДО (

и ) и № г л Д и /1'

3 0.5 1 15 2 2,5 3 3.5 4 4,5 5 5.5 6 5.5 7 7.5 8 8.5 9 9.5 10 10.5 И Время, сек

Г 1,, /У

ч, « 1,

ц \

О 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 Б 3.5 9 9.5 10 10.5 11

Рис. 7 ДСЧ и энергия рассеянного сигнапа для гидросича

счет этого энергия, рассеянная при подъеме второго пузыря, значительно больше, нежели первого. Искажение формы второго пузыря вызвано воздействием на него возмущенной жидкости, находящейся в следе первого.

В дальнейшем анализу подвергались именно те сигналы, которые получены при рассеянии на пленке, возмущенной следом предыдущего пузыря. Ниже, на одном графике (рис. 8, 9) представлен набор зависимостей ДСЧ от времени, полученных для пузырей различной протяженности. Эти зависимости совмещены так, что их начала (моменты времени, когда ДСЧ близко к нулю) совпадают. Для самого длинного пузыря, как видно из рисунка, временная зависимость (график) располагается в интервале времени от 0 до 1 сек. Видно также, что графики от более длинных пузырей продолжают графики от более коротких. Таким образом, при наложении друг на друга эти графики образуют зависимость, которая характеризует максимальную скорость стекания пленки как функцию ее длины.

Рис.8 Временные реализации ДСЧ для нескольких пузырей различной длительности в воде.

Рис.9 Временные реализации ДСЧ для нескольких пузырей различной длительности в гидросиле.

— 20011» —- 400

—■ 600II»

— 1000(Го

— 1400 ПК

Время. еж

-10 0 Ю

.....400 пй

"■■600 т ¡000 т —1400 ш;

|>1>еш, сек

Известно [7], что для стекающей пленки, средняя скорость жидкости относительно пузыря в направлении, противоположном его движению, на расстоянии Ь от его верхней точки, определяется только инерционными эффектами и равна д1. Тогда относительно стенки скорость стекания пленки будет определяться выражением: vm = — va¡, где скорость всплытия пузыря. Эта величина

определялась с помощью видеосъемки. В результате для нее получились экспериментальные значения, лежащие в интервале 20 - 26 см/сек. Существует теоретическая оценка скорости всплытия пузыря [7], которая определяется выражением = 0,345^дО, Т) - диаметр трубы. В нашем случае теория дает значение = 22 с.и/сех, что хорошо согласуется с экспериментом.

Теоретическая зависимость скорости стекаиия пленки, пересчитанная во временную зависимость доплеровского смещения частоты подстановкой L=tv,, также нанесена на рисунки. Доплеровское смещение частоты Fd рассчитывалось по формулe:Fd=(2fi/C)V^. Здесь f0= 1 МГц- несущая частота зондирующих импульсов, С= Cflex=3060 м/с - скорость излучаемой изгибной волны Лэмба в стенке трубы, VIL, -скорость стекающей пленки. Теоретическая зависимость и экспериментальные графики, как для воды, так и для гидросила, практически совпадают почти в течение всего процесса стекания пленки. Совпадение измеренного и оцененного по формулам доплеровского сдвига достигается при скорости С, равной скорости изгибной волны Cj7„. Отсюда следует, во-первых, что доплеровское смещение частоты возникает именно за счет рассеяния изгибной волны, во-вторых, что рассеивающие эту волну неоднородности движутся со скоростью (относительно стенки) близкой к средней скорости стекания пленки.

Для воды совпадение экспериментальных графиков с расчетным наблюдается в диапазоне доплеровских смещений от 0 до 1000 Гц. Такое смещение частоты соответствует довольно значительной по величине скорости стекания относительно стенки VOTaI=l,5 м/с, а относительно головной части всплывающего пузыря - даже 1,7 м/с. Для гидросила экспериментальные графики выходят на некоторый постоянный уровень, для которого ДСЧ составляет около 800 Гц. Этот экспериментальный результат свидетельствует, что, в конечном счете, достигается некая предельная скорость стекания, при которой вес пленки полностью уравновешивается касательными напряжениями на стенке. Чем больше вязкость жидкости, тем раньше достигается это предельное значение скорости. Вязкость гидросила больше чем у воды, поэтому процесс выхода стекания пленки на постоянную скорость для него уже заметен. Для воды же на интервалах времени, доступных в эксперименте, эта стадия процесса стекания просто еще не достигнута.

Четвертая глава посвящена методу определения расходов и поиску характеристик сигналов, пригодных для определения расходов компонентов потока.

В параграфе 4.1 представлен общий подход к определению расходных параметров двухфазного потока, который заключается в нахождении решения Qj,

(Pi= If.Qa) \P2= 4>2(Qr-Qg)'

где, Pi, Pi - набор характеристик сигналов, выбранных для определения параметров потока, 0i(...), <t>2(---) - эмпирические, полуэмпирические или теоретические функциональные зависимости характеристик сигналов от расходов газа Qe и жидкости Qf. Первая трудность при решении этой обратной задачи заключается в том, что правые части уравнений определяются в результате калибровочных измерений ультразвуковой системой и поэтому точно известны только в ограниченном числе точек. Вторая связана со структурой (видом) уравнений, которая может оказаться,

Qg системы уравнений вида:

во-первых, нелинейной, а во-вторых - допускающей неоднозначные решения. Таким образом, при выборе характеристик сигналов Р, необходимо отбирать те из них, которые монотонно убывают (или нарастают) при изменении расходов компонентов потока.

Решение обратной задачи можно получить различными методами. Основой метода, использованного в работе, является то, что хотя способ постановки опыта, т.е. задание параметров потока и измерение соответствующих им характеристик сигналов определяет причинно-следственную связь между ними, с математической точки зрения оба пространства равноправны. Поэтому для «измерения» параметров потока, т.е. - для определения точки в пространстве параметров потока по значениям характеристик сигналов, не содержащимся в калибровочной базе, может быть использован алгоритм 2-х мерной интерполяции точки в пространстве параметров потока по точке в пространстве характеристик сигналов, полученных в результате измерения. Данный метод был реализован с помощью проведения 2-х мерной интерполяции с использованием функции '^пскШа" из пакета МАТЬАВ.

Параграфы 4.2 и 4.3 посвящены определению расходов газа и жидкости для «пузырькового» и «снарядного» режимов течения. Для изучения работы ультразвуковой системы в условиях существования в трубе реального двухфазного газожидкостного потока и построения калибровочной базы, то есть набора функциональных зависимостей характеристик сигналов от расходов Р( = ФгСФ/'СдХ была проведена серия экспериментов на сертифицированном промышленном проливном стенде ООО «АЗМИП». На стенде существовала возможность проводить измерения в достаточно широких пределах изменения расходов жидкой и газовой фаз потока (0-800 м3/сут), и создавать как «пузырьковый» так и в «снарядный» режимы течения.

В результате исследований было определено, что наиболее характерные отличия между «пузырьковым» и «снарядным» режимами течения проявляются во временных зависимостях ДСЧ. На рис.10 и 11 приведены временные реализации ДСЧ для 5-ых элементов дальности двух противоположно направленных лучей антенны.

Рис.10 Зависимости ДСЧ (ГгО от времени для Рис.11 Зависимости ДСЧ [Гц] от времени для Q/JOm3 /сут. и Qg-]00 м1 /сут. «Снарядный» Qf=200 м! /сут.. газа Qg=15 м* /сут. режим течения. «Пузырьковый» режим течения.

Можно заметить, что для «снарядного» режима зависимость ДСЧ является квазипериодической функцией времени и отслеживает чередующуюся структуру «газовый пузырь - жидкая перемычка». В противоположность к «снарядному» режиму течения, в «пузырьковом» режиме временная зависимость ДСЧ не имеет квазипериодической структуры, и представляет собой шумовой процесс. Среднее значение ДСЧ Рр при этом почти линейно зависит от расхода жидкости (рис.12) и определяет некоторую среднюю скорость движения потока. Расход газа влияет на величину £)р, но весьма незначительно.

При «снарядном» режиме течения (рис.10) хорошо виден временной сдвиг между зависимостями ДСЧ противоположных лучей антенны, направленных по потоку и против него. Поэтому, вследствие коррелированное™ между собой зависимостей ДСЧ, существует возможность непосредственного определения скорости движения структуры «газовый пузырь - жидкая перемычка» с помощью соотношения и=иЛкор, где Лкор - временной сдвиг между реализациями ДСЧ, определяемый по сдвигу максимума взаимной корреляционной функции, -некоторая база, которую можно определить по скорости распространения изгибной волны по стенке и задержкам элементов дальности, из которых взяты последовательности ДСЧ. На рис.13 представлена экспериментальная зависимость параметра II от суммарного расхода компонентов потока, которая представляет собой почти линейно нарастающую функцию.

Таким образом, с помощью ультразвуковой системы, при использовании одного

Г— ' ]-гк"Т

я Ш КО 23 МО та

.-V

Я> 10С 15С 2£» 255 Ш 350 400

О'.иУч'

Рис. 12 Семейства зависимостей Ир(Ц/) при (2р=сопМ (м3/сут). «Пузырьковый» режим (дисперсия Орр ~ 4 Гц)

а гк ж ¡со -ж

Рис.13 Зависимость скорости структуры потока С/ от суммарного расхода газа и жидкости. «Снарядный» режим (дисперсия в и ~ 0.08 м/с)

датчика, для обоих режимов удалось получить характеристики сигналов, определяющие среднюю скорость течения потока, которая связанна непосредственным образом с расходами фаз потока. Однако, для определения расходов фаз в каждом из режимов, необходимо использование минимум двух характеристик сигнала. Поэтому в работе для «пузырькового» и «снарядного» режимов течения были найдены и исследованы другие характеристики, обладающие условием монотонности при изменении расходов или (2Г Для «снарядного» режима наиболее информативной является зависимость ДСЧ от времени, поэтому все

характеристики были получены из нее. А именно, это минимальное значение доплеровского сдвига частот /•',„,■„, среднее значение доплеровского сдвига Эр, и отношение длительностей участков ДСЧ с разными знаками ГУГ+.

Для «пузырькового» режима наиболее информативными оказались профили энергии. Характеристики профилей энергии тем или иным образом связаны с затуханием в потоке, которое, в свою очередь, связано с газосодержанием. Использовались следующие характеристики профилей энергии: Еп„2/Епа1-

_ /пЕав(п) (Зи

отношение величины второго максимума к первому, лет = ^ - координата

«центра масс» профиля, которая также описывает соотношение между его пиками.

В результате, имея в распоряжении несколько различных характеристик сигналов для каждого из режимов, для их сочетаний (Р;, Р2), были оценены погрешности определения расходов газа и жидкости. По оценкам погрешностей производилось сравнение используемых пар характеристик. Определение погрешности решения

у = 0 ¿Гч * 21

..............:......... | 6

у : :

> зс ■ /

1 280

^ 153 Ж

у ;

100

8 160 233 250 ЗСС 350 4 Расход жидкости. ЕыстаапвниыА из стенд». иУс .1 3 450 "о 30 4 50 60 70 8 ставленный на с!гиде. м&'с Э0 106

Рис. 14 Результаты восстановления расхода жидкости. «Пузырьковый» р-м. Рис. 15 Результаты восстановления газа. «Пузырьковый» р-м. расхода

обратной задачи производилось путем удаления из калибровочной базы одной точки и ее последующего восстановления по предварительно отобранным

2С0

У = <¡97^x-i /

у >

ж у

\ У ■У

Рк*од Живгсок. ?ыстзепе-«ый «э

«.л л / /

ч % /

Л?

Рис.16 Результаты восстановления расхода жидкости. «Снарядный» р-м.

Рис. 17 Результаты восстановления расхода газа. «Снарядный» р-м.

характеристикам сигналов с помощью алгоритма интерполяции. Оказалось, что в «пузырьковом» режиме течения, наименьшая погрешность одновременного определения расходов жидкости (±1,3%) и газа (±10%) получается при использовании для их определения среднего ДСЧ (Р]=1)р), и координаты «центра масс» профиля энергии (Р2-Хст). Результаты восстановления расходов жидкости и газа с использованием этих параметров показаны на рис.14 и рис.15 соответственно.

В «снарядном» режиме течения наименьшая погрешность одновременного определения расходов жидкости (±1,3%) и газа (±1,7%) получается при использовании для их определения скорости потока, измеряемой корреляционным методом (Р1=Ц), и среднего значения ДСЧ (Р2=Вр). Результаты восстановления расходов жидкости и газа с использованием этих параметров показаны на рис. 16 и рис.17 соответственно.

В заключении приводятся основные результаты работы.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ ДИССЕРТАЦИИ

1. Разработана упрощенная модель формирования профилей энергии рассеянных ультразвуковых сигналов в потоке, переносящем дискретные неоднородности. Результаты расчетов по этой модели находятся в согласии с данными опытов.

2. Экспериментально реализован метод определения размеров и количества всплывающих газовых пузырьков в разреженном потоке, основанный на фильтрации их собственных колебаний, возникающих в результате их возбуждения короткими импульсами.

3. В «снарядном» режиме течения, при прохождении мимо датчика жидкой перемычки, рассеяние происходит на содержащихся в ней пузырьках газа, а при прохождении газовой полости - на неоднородностях толщины стекающей пленки. Неоднородности толщины пленки движутся относительно стенки со скоростью близкой к скорости стекания пленки. Поэтому доплеровский сдвиг частоты рассеянного на пленке сигнала соответствует средней скорости ее стекания.

4. Выявлено, что уровень сигнала, рассеянного стекающей пленкой, увеличивается при росте интенсивности гидродинамических возмущений жидкости, в которой происходит всплытие газовой полости. В эксперименте гидродинамические возмущения жидкости возникали вследствие всплытия первого пузыря в невозмущенной жидкости, последующие пузыри всплывали уже в жидкости, возмущенной движением предыдущего.

5. На основании модели, связывающей параметры сигналов, рассеянных в газожидкостном потоке, с расходами фаз потока выяснено, что для определения расходов, в «пузырьковом» режиме течения, наиболее информативными являются зависимости от дальности энергии и средней частоты рассеянного в потоке сигнала. В

«снарядном» режиме наиболее информативными являются временные реализации средней частоты сигнала, получаемые с различных удалений от датчика.

6. В «пузырьковом» режиме течения, наименьшая погрешность одновременного определения расходов жидкости и газа, с применением одного приемно-излучающего датчика, получается при использовании для их определения среднего доплеровского сдвига частоты (ДСЧ), и координаты «центра масс» профиля энергии рассеянного сигнала. В «снарядном» режиме течения наименьшая погрешность получается при использовании для их определения скорости потока, измеряемой корреляционным методом, и среднего значения ДСЧ.

СПИСОК ЦИТИРУЕМОЙ ЛИТЕРАТУРЫ

1. Кремлевский, П.П. Расходомеры и счетчики количества веществ: Справочник: Кн. 2 / П.П. Кремлевский - СПб.: Политехника, 2004.

2. Фролов А.И. Проблемы измерения продукции скважин нефтяных месторождений // Датчики и системы, 2001, № 9, с. 37-40.

3. Муякшин С.И., Односевцев В.А., Орлов И.Я. О возможности использования изгибных волн при создании акустических преобразователей для систем неразрушающего контроля // Изв. вузов. Радиофизика, 2000, t.XLIII, № 12, с. И 00-1108

4. В.И.Косарев, Л.Н.Мухин,С.И.Муякшин и др. Способ определения расхода компонентов двухфазного потока и система для его осуществления // Патент на изобретение № RU 2339915 С1. Опубликовано 2008.11.27 Бюл. № 33

5. Исимару А. Распространение и рассеяние волн в случайно-неоднородных средах, т. 1 // А. Исимару - М„ Мир,1981.-280 с.

6. Keller J.В., Miksis М., Bubble oscillation of large amplitude. J. Acoust. Soc. Am., 1980, 68 pp. 628-633.

7. Уоллис, Г. Одномерные двухфазные течения / Г. Уоллис - М., Мир, 1972.-440 с.

8. Р.Довиак, Д.3рнич. Допплеровские радиолокаторы и метеорологические наблюдения -Л., Гидрометеоиздат, 1988 г.

СП11КОК ПУБЛИКАЦИЙ ПО ДИССЕРТАЦИИ

Список публикаций в журналах, рекомендованных ВАК:

1. Муякшин С.И., Владимиров П.А. Ультразвуковая доплеровская система для измерения расхода газожидкостных потоков // «Датчики и системы», 2009, №8, с. 26-30.

2. Мансфельд А.Д., Волков Г.П., Санин А. Г., Владимиров И.А. Импульсная ультразвуковая спектроскопия газовых пузырьков // «Акустический Журнал», 2010, том 56, № 3, с. 323-332

3. Владимиров И.А., Муякшин С.И. Экспериментальное исследование рассеяния ультразвуковых волн Лэмба в пластине, покрытой стекающей пленкой жидкости переменной толщины // Вестник Нижегородского государственного университета им. Н.И. Лобачевского, 2011, № 5(3), с. 67-74.

!1уГш!кация статей в трудах конференций:

4. Владимиров H.A., Волков Г.П., Мансфельд А.Д. Импульсные ультразвуковые методы обнаружения газовых пузырьков// Труды XII научной конференции по радиофизике, ИНГУ, 2008

5. Муякшин С. И., Владимиров H.A. Метод определения расхода фаз газожидкостного потока в вертикальной трубе при "снарядном" режиме течения // Труды XIII научной конференции по радиофизике, ННГУ, 2009

6. Муякшин С. И., Владимиров H.A. Ультразвуковая система диагностики газожидкостных потоков в закрытой трубе // Труды XIV научной конференции по радиофизике, ННГУ, 2010

7. Владимиров И.А., Муякшин С.И. Ультразвуковые методы диагностики газожидкостных потоков // Сборник трудов XII Всероссийской школы семинара «Волновые явления в неоднородных средах», 2010г., Звенигород, В электронном виде на CD - диске

8. Владимиров И.А., Муякшин С.И. Об особенностях рассеяния звука на потоке пузырьков в трубе, заполненной жидкостью // XV Нижегородская сессия молодых ученых. Труды молодых ученых по естественнонаучным дисциплинам, Нижний Новгород, 2010

Отчеты о НИР

9. Отчет о НИР № 02201153204 от 09.03.2011 «Разработка перспективных ультразвуковых методов определения параметров трехкомпонентных двухфазовых потоков», рук. И.А. Владимиров

СОДЕРЖАНИЕ ДИССЕРТАЦИИ Введение

ГЛАВА 1 Двухфазные потоки в технике: обзор гидродинамики и методов измерения

1.1 Классификация

1.2 Количественное описание

1.3 Режимы течения

1.4 Методы измерения расходных характеристик Основные результаты главы 1

ГЛАВА 2 Описание измерительной системы, экспериментальных установок и методов измерения

2.1 Принцип работы датчика с фазированной антенной решеткой

2.2 Блок-схема измерительной системы

2.3 Экспериментальная установка для изучения «пузырькового» потока

2.4 Экспериментальная установка для изучения рассеяния волн Лэмба в пластине, покрытой стекающей пленкой жидкости

2.5 Проливной стенд ООО «АЗМИП»

2.6 Когерентно-импульсный метод измерения скорости потока

2.7 Исследование акустического поля, создаваемого накладной ФАР в круглой трубе с однородной жидкостью

Основные результаты главы 2

ГЛАВА 3 Теоретическое и экспериментальное исследование процессов рассеяния ультразвука в «пузырьковом» потоке и стекающей пленке жидкости

3.1 Исследование процессов рассеяния ультразвуковых импульсов в потоке пузырьков

3.1.1 Задача определения интенсивности рассеянного сигнала, принимаемого ФАР

3.1.2 Экспериментальные исследования процессов рассеяния ультразвуковых импульсов в потоке пузырьков, сравнение эксперимента с теорией

3.1.3 Определение размеров мелких газовых пузырьков в разреженном потоке с помощью частотной полосовой фильтрации

3.2 Исследование процессов рассеяния ультразвуковых волн Лэмба в пластине, покрытой стекающей пленкой жидкости

3.2.1 Объяснение механизма рассеяния ультразвуковых волн Лэмба, распространяющихся в пластине, покрытой деформированной пленкой жидкости

3.2.2 Экспериментальные исследования рассеяния ультразвуковых волн Лэмба в пластинке, покрытой стекающей пленкой жидкости. Сравнение эксперимента с теорией

Основные результаты главы 3

ГЛАВА 4 Определение расходов газа и жидкости в двухфазном потоке

4.1 Общий подход к определению расходов фаз потока

4.2 Определение расходов при «пузырьковом» режиме течения

4.2.1 Анализ информативных характеристик рассеянных и прошедших сигналов для «пузырькового» потока

4.2.2 Решение обратной задачи - определение расходов фаз потока по характеристикам сигналов

4.3 Определение расходов при «снарядном» режиме течения

4.3.1 Особенности гидродинамики «снарядного» режима течения и временная изменчивость рассеянных сигналов

4.3.2 Анализ информативных характеристик рассеянных сигналов для «снарядного» потока

4.3.3 Решение обратной задачи - определение расходов фаз потока по характеристикам сигналов Основные результаты главы 4

Заключение

Список использованной литературы

Подписано в печать 28.04.2012. Формат 60x84 1/16. Бумага офсетная. Печать офсетная. Гарнитура Тайме. Усл. печ. л. 1. Заказ № 299. Тираж 100.

Отпечатано в Центре цифровой печати Нижегородского госуниверситета им. Н.И. Лобачевского 603950, г. Нижний Новгород, пр. Гагарина, 23

Содержание диссертации автор исследовательской работы: кандидата физико-математических наук, Владимиров, Илья Александрович

Введение.

ГЛАВА 1 Двухфазные потоки в технике: обзор гидродинамики и методов измерения

1.1 Классификация.

1.2 Количественное описание.

1.3 Режимы течения.

1.4 Методы измерения расходных характеристик

Основные результаты главы 1.

ГЛАВА 2 Описание измерительной системы, экспериментальных установок и методов измерения.

2.1 Принцип работы датчика с фазированной антенной решеткой.

2.2 Блок-схема измерительной системы.

2.3 Экспериментальная установка для изучения «пузырькового» потока.

2.4 Экспериментальная установка для изучения рассеяния волн Лэмба в пластине, покрытой стекающей пленкой жидкости.

2.5 Проливной стенд ООО «АЗМИП».

2.6 Когерентно-импульсный метод измерения скорости потока.

2.7 Исследование акустического поля, создаваемого накладной ФАР в круглой трубе с однородной жидкостью.

Основные результаты главы 2.

ГЛАВА 3 Теоретическое и экспериментальное исследование процессов рассеяния ультразвука в «пузырьковом» потоке и стекающей пленке жидкости.

3.1 Исследование процессов рассеяния ультразвуковых импульсов в потоке пузырьков.

3.1.1 Задача определения интенсивности рассеянного сигнала, принимаемого ФАР.

3.1.2 Экспериментальные исследования процессов рассеяния ультразвуковых импульсов в потоке пузырьков, сравнение эксперимента с теорией.

3.1.3 Определение размеров мелких газовых пузырьков в разреженном потоке с помощью частотной полосовой фильтрации.

3.2 Исследование процессов рассеяния ультразвуковых волн Лэмба в пластине, покрытой стекающей пленкой жидкости.

3.2.1 Объяснение механизма рассеяния ультразвуковых волн Лэмба, распространяющихся в пластине, покрытой деформированной пленкой жидкости.

3.2.2 Экспериментальные исследования рассеяния ультразвуковых волн Лэмба в пластинке, покрытой стекающей пленкой жидкости. Сравнение эксперимента с теорией.

Основные результаты главы

ГЛАВА 4 Определение расходов газа и жидкости в двухфазном потоке.

4.1 Общий подход к определению расходов фаз потока.

4.2 Определение расходов при «пузырьковом» режиме течения.

4.2.1 Анализ информативных характеристик рассеянных и прошедших сигналов для «пузырькового» потока.

4.2.2 Решение обратной задачи - определение расходов фаз потока по характеристикам сигналов.

4.3 Определение расходов при «снарядном» режиме течения.

4.3.1 Особенности гидродинамики «снарядного» режима течения и временная изменчивость рассеянных сигналов.

4.3.2 Анализ информативных характеристик рассеянных сигналов для «снарядного» потока.

4.3.3 Решение обратной задачи - определение расходов фаз потока по характеристикам сигналов.

Основные результаты главы 4.

Введение диссертация по физике, на тему "Ультразвуковая диагностика восходящих газожидкостных потоков с использованием распределенного электромеханического преобразователя"

Актуальность исследования

Двухфазные потоки в виде смеси жидкости и газа встречаются в разных областях техники: в добыче нефти и газа, в химической и пищевой промышленности, в энергетике, в коммунальном хозяйстве. Для учета добытого углеводородного сырья, управления технологическими установками, контроля стоков, необходимо измерение расхода жидкой и газовой фаз по отдельности. В нефтяной промышленности проблема измерения расходов таких потоков возникла еще в начале 1980-х годов. Необходимо, чтобы приборы, измеряющие расход многофазного потока (расходомеры) удовлетворяли следующим требованиям: обладали приемлемой точностью и надежностью, не создавали возмущений в потоке, не зависели от режима потока и были пригодными для применения во всем возможном диапазоне изменения расходов компонентов потока. Несмотря на большое количество технических решений [1-4,8-10], предложенных в последние годы, ни один из двухфазных расходомеров, появившихся на рынке, всем указанным требованиям все же не отвечает [2]. В настоящее время для определения расходов потока в основном используется предварительная дегазация потока [2,3]. Такие установки дороги и сложны в эксплуатации. Поэтому создание надежных, точных и не слишком дорогих многофункциональных приборов для измерения расхода компонентов потока без их предварительного разделения является актуальной задачей. Отметим, что наиболее перспективными методами диагностики двухфазных потоков, являются активные акустические методы. В представленной работе для изучения и определения параметров двухфазных потоков используется ультразвуковая система, основанная на методе измерения скорости неоднородностей с помощью эффекта Доплера. Принцип работы системы приводиться в [5,6]. В работе [5] описан способ ввода ультразвукового пучка в поток, основанный на использовании эффекта трансформации изгибной волны в стенке трубы (волны Лэмба) в продольную звуковую волну в жидкости. Там же предложено использовать для формирования волны Лэмба накладную фазированную антенную решетку из пьезоэлементов. В патенте [6] описан способ измерения параметров двухфазного потока и система для его реализации, в которой используется описанный выше способ ввода ультразвука в поток и эффект Доплера. Однако в [6] раскрывается только лишь принципиальная возможность определения расходов фаз потока. Поэтому полная реализация возможностей системы, включая достижение в большом диапазоне изменения расходов двухфазного потока достаточной точности измерения [7], требует систематического и углубленного исследования процессов распространения и рассеяния ультразвука в газожидкостных смесях (ГЖС).

Целью работы является детальное изучение закономерностей рассеяния ультразвуковых импульсов в восходящем газожидкостном потоке, в который они вводятся через стенку трубы с помощью накладного пьезоэлектрического преобразователя, выполненного в виде фазированной антенной решетки (ФАР), и изучение возможности использования этих закономерностей для раздельного определения расходов газовой и жидкой фаз.

Для достижения этой цели в работе были поставлены следующие основные задачи:

Положения, выносимые на защиту:

1. Использование разработанной модели формирования зависимости энергии рассеянных ультразвуковых сигналов от дальности (профилей), в пузырьковом» режиме газожидкостного потока, позволяет определять газосодержание по форме профиля энергии.

Научная новизна работы заключается в следующем:

Достоверность результатов

Научные положения и выводы, представленные в диссертационной работе, аргументированы, обоснованы и базируются на фундаментальных теоретических положениях общей акустики и теории распространения и рассеяния волн, а также результатах целенаправленных экспериментальных исследований.

Практическое значение

Полученные в работе теоретические и экспериментальные результаты могут быть положены в основу разрабатываемых ультразвуковых измерителей расходов двухфазных потоков, а созданная и опробованная экспериментальная установка может рассматриваться как прототип ультразвукового расходомера.

Апробация работы

Основные результаты работы докладывались на XII Всероссийском семинаре «Волновые явления в неоднородных средах», (Звенигород, МГУ, 2010), XII, XIII и XIV научных конференциях по радиофизике (Нижний Новгород, ННГУ, 2008-2010), XV и XVI Нижегородских сессиях молодых ученых (Нижний Новгород, 2009, 2010).

Личный вклад автора

На протяжении всей исследовательской работы диссертант принимал непосредственное участие в постановке задач, проведении экспериментов и расчетов, построении моделей, обсуждении и физической интерпретации результатов исследований. Исключением является раздел 3.1.3, в котором автору принадлежат только проведение экспериментов и обсуждение результатов.

Публикации

Список публикаций в журналах, рекомендованных ВАК: Муякшин С.И., Владимиров И.А. Ультразвуковая доплеровская система для измерения расхода газожидкостных потоков // «Датчики и системы», 2009, №8, с. 26-30.

Мансфельд А.Д., Волков Г.П., Санин А. Г., Владимиров И.А. Импульсная ультразвуковая спектроскопия газовых пузырьков // «Акустический Журнал», 2010, том 56, №3, с. 323-332.

Владимиров И.А., Муякшин С.И. Экспериментальное исследование рассеяния ультразвуковых волн Лэмба в пластине, покрытой стекающей пленкой жидкости переменной толщины // Вестник Нижегородского государственного университета им. Н.И. Лобачевского, 2011, № 5(3), с. 67-74.

Список статей в трудах конференций: Владимиров И.А., Волков Г.П., Мансфельд А.Д. Импульсные ультразвуковые методы обнаружения газовых пузырьков// Труды XII научной конференции по радиофизике, ННГУ, 2008

Муякшин С. И., Владимиров И.А. Метод определения расхода фаз газожидкостного потока в вертикальной трубе при "снарядном" режиме течения // Труды XIII научной конференции по радиофизике, ННГУ, 2009

Владимиров И.А., Муякшин С. И Ультразвуковая система диагностики газожидкостных потоков в закрытой трубе // Труды XIV научной конференции по радиофизике, ННГУ, 2010

Владимиров И.А., Муякшин С.И. Ультразвуковые методы диагностики газожидкостных потоков // Сборник трудов XII Всероссийской школы семинара «Волновые явления в неоднородных средах», 2010г., Звенигород, В электронном виде на СО - диске

Владимиров И.А., Муякшин С.И. Об особенностях рассеяния звука на потоке пузырьков в трубе, заполненной жидкостью // XV Нижегородская сессия молодых ученых. Труды молодых ученых по естественнонаучным дисциплинам, Нижний Новгород, 2010 Отчеты о НИР

Отчет о НИР № 02201153204 от 09.03.2011 «Разработка перспективных ультразвуковых методов определения параметров трехкомпонентных двухфазовых потоков», рук. И.А. Владимиров Структура и объем работы.

Диссертация состоит из введения, четырех глав, заключения и списка цитируемой литературы. Общий объем диссертации составляет 112 страниц, включая 79 рисунков, 3 таблицы и список литературы из 73 наименований.

Заключение диссертации по теме "Акустика"

Основные результаты главы 4

- В результате проведенных экспериментов на специализированном проливном стенде, с помощью УЗ системы получены зависимости энергии Еп и доплеровского сдвига частоты 7ч/ от номера ЭД и от времени, а так же набор параметров, определенных по этим зависимостям, при различных сочетаниях расходов газа и жидкости.

- На основании модели, связывающей параметры сигналов, рассеянных в газожидкостном потоке, с его интегральными расходными характеристиками получены следующие результаты:

1. Использование способа определения содержания газа в «пузырьковом» потоке, с использованием профилей энергии только рассеянного сигнала, принимаемого одним излучающим датчиком, дает приблизительно такую же погрешность измерения, как и метод определения количества газа по затуханию прошедшего через поток сигнала с использованием двух, противоположно расположенных датчиков.

2. В «пузырьковом» режиме течения, наименьшая погрешность одновременного определения расходов жидкости и газа, с применением одного приемно-излучающего датчика, получается при использовании для их определения среднего доплеровского сдвига частоты (ДСЧ), и координаты «центра масс» профиля энергии рассеянного сигнала.

3. В «снарядном» режиме течения наименьшая погрешность одновременного определения расходов жидкости и газа получается при использовании для их определения скорости потока, измеряемой корреляционным методом, и среднего значения ДСЧ.

4. Для определения расходных параметров потока, в «пузырьковом» режиме течения, наиболее информативными являются зависимости энергии и средней частоты от номера ЭД (профили), а в «снарядном» режиме - временные зависимости энергии и средней частоты из разных элементов дальности, относящихся к разным лучам.

Заключение

В заключении отметим следующие основные результаты, полученные в диссертационном исследовании:

- Испытана ультразвуковая система, содержащая: 1) датчик с ФАР; 2) набор модулей в форм-факторе NI SXCI для генерации, усиления мощности ультразвуковых импульсов, их приема и преобразования в НЧ квадратурные компоненты; 3) шасси для их установки; 4) универсальную карту ввода-вывода NI PXI-6251 (USB-6251); 5) персональный компьютер.

- С использованием среды графического программирования Lab VIEW разработано программное обеспечение для сбора, отображения и записи данных с этой системы в реальном времени. Данные записываются в виде последовательностей моментов скользящего спектра для нескольких задержек и двух лучей.

- Создана модификация лабораторной установки для моделирования пузырькового потока в воде и гидросиле (при неподвижном столбе жидкости). В установке использована прозрачная стеклянная трубка, что обеспечило возможность визуального наблюдения потока. Для этой установки разработан метод определения истинного относительного газосодержания основанный на измерении высоты подъема уровня (для воды) или отслеживания динамики изменения высоты пены (для гидросила).

- Создана модификация лабораторной установки для изучения рассеяния ультразвуковых волн на стекающей пленке жидкости, формируемых датчиком с ФАР. Пленка возникает при всплытии в трубке с неподвижной жидкостью больших воздушных пузырей, перекрывающих все ее сечение. Пузыри генерируются с помощью программно-управляемого клапана.

- Исследование структуры акустического поля, создаваемого накладной ФАР в круглой трубе с однородной жидкостью, показало, что в жидкости формируется ультразвуковой пучок большой протяженности (вдоль трубки), который несколько раз отражается от стенки трубы. Ширина его примерно равна ширине ФАР.

- Разработана модель формирования профилей энергии рассеянных ультразвуковых сигналов в потоке, переносящем дискретные неоднородности. Результаты расчетов по этой модели находятся в согласии с данными опытов.

- При различных величинах газосодержания измерены профили энергии рассеянных сигналов для газожидкостной смеси, представляющей собой всплывающие пузырьки газа в воде и гидросиле.

- Получены экспериментальные зависимости между амплитудами пиков профиля энергии эхосигналов и газосодержанием. Найдено, что для воды и гидросила отношение пиков монотонно растет при увеличении газосодержания (в определенном диапазоне).

- Экспериментально реализован метод определения размеров и количества всплывающих газовых пузырьков в достаточно разреженном потоке на основе фильтрации принимаемого сигнала.

- Получены экспериментальные зависимости средней частоты, энергии и ширины скользящего спектра эхосигнала от времени, при прохождении мимо УЗ датчика всплывающих газовых полостей.

- Для «снарядного» режима течения объяснен процесс рассеяния волн Лэмба. Из объяснений следует, что обратное рассеяние поверхностных волн происходит на деформированной (переменной по толщине) пленке.

Проведенные экспериментальные исследования показали, что доплеровское смещение частоты сигнала, рассеянного на стекающей пленке, соответствует оценке средней скорости стекания пленки, основанной на гидродинамике всплытия газовой полости в вертикальной тубе. Отсюда следует, что рассеивающие звук неоднородности движутся со средней скоростью стекающей пленки.

- Выявлено, что уровень сигнала, рассеянного пленкой, увеличивается при росте интенсивности гидродинамических возмущений жидкости, в которой происходит всплытие газовой полости (В эксперименте гидродинамические возмущения жидкости возникали вследствие всплытия первого пузыря в невозмущенной жидкости, второй пузырь всплывал уже в жидкости, возмущенной движением первого)

- В результате проведенных экспериментов на специализированном проливном стенде, для «пузырькового» и «снарядного» режимов течения, с помощью УЗ системы получены зависимости энергии Еп и доплеровского сдвига частоты Г(1 от номера ЭД и от времени, а так же набор параметров, определенных по этим зависимостям, при различных сочетаниях расходов газа и жидкости.

Список источников диссертации и автореферата по физике, кандидата физико-математических наук, Владимиров, Илья Александрович, Нижний Новгород

1. Кремлевский, П.П. Расходомеры и счетчики количества веществ: Справочник: Кн. 2 // П.П. Кремлевский СПб.: Политехника, 2004.

2. Фролов А.И. Проблемы измерения продукции скважин нефтяных месторождений// Датчики и системы, 2001, №9, с. 37-40.

3. Дробков В.П. Разработка и исследование ультразвуковых методов и информационно-измерительной системы измерения расхода нефтеводогазового потока: автореф. Диссертации доктора технических наук : 05.11.16 / Дробков Владимир Петрович; Москва, 2007, 43 с.

4. Отчет о НИР № 02201153204 от 09.03.2011 «Разработка перспективных ультразвуковых методов определения параметров трехкомпонентных двухфазовых потоков», рук. И.А. Владимиров

5. Муякшин С.И., Односевцев В. А., Орлов И.Я. О возможности использования изгибных волн при создании акустических преобразователей для систем неразрушающего контроля // Изв. вузов. Радиофизика, 2000, T.XLIII, № 12, с.1100-1108

6. В.И.Косарев, Л.Н.Мухин,С.И.Муякшин и др. Способ определения расхода компонентов двухфазного потока и система для его осуществления // Патент на изобретение № RU 2339915 С1. Опубликовано 2008.11.27 Бюл. № 33

7. ГОСТ Р 8.615-2005. Измерения количества извлекаемой из недр нефти и нефтяного газа // Москва, Стандартинформ, 2006, 3 с.

8. Мельников В.И. Акустические методы диагностики газожидкостных потоков. / В.И. Мельников, В.П. Дробков, В.В. Контелев М.: Энергоатомиздат, 2006. 382 с.

9. Цейтлин В.Г. Расходоизмерительная техника. М.:Стандарты, 1977. -239 с.

10. Абрамов Г.С. Практическая расходометрия в промышленности. / Г.С. Абрамов, A.B. Барычев, М.И. Зимин М.: ОАО «ВНИИОЭНГ», 2000. -472 с.

11. Лабунцов, Д.А. Механика двухфазных систем. Учебное пособие для вузов / Д.А. Лабунцов, В.В. Ягов Москва: Издательство МЭИ, 2000. - 374 с: ил. ISBN 5-7046-0758-6

12. Уоллис, Г. Одномерные двухфазные течения / Г. Уоллис М., Мир, 1972.440 с.

13. Мамаев В. А. Гидродинамика газо-жидкостных смесей в трубах / Мамаев В. А., Одишария Г. Э., Семенов Н. И., Точигин А. А. // М., «Недра», 1969. 208 с.

14. Нигматулин Р. И. Динамика многофазных сред. Т. 1,2. М.:Наука, 1987.

15. Кутателадзе С. С., Стырикович М. А. Гидродинамика газожидкостных систем. М.:Энергия, 1976.

16. Мамаев В.А. и др. Движение газожидкостных смесей в трубах. М., "Недра", 1978, 344 с.

17. Танеев Ф.К. «Автоматизированные групповые установки типа «Спутник» для покомпонентного измерения продукции нефтяных скважин» // ЭИ серия «Автоматизация и телемеханизация в нефтяной промышленности», ВНИИОЭНГ, М. вып. 10, 1991.

18. Ковшов В.Д. Автоматизация технологических процессов. 4.2 Уфа.: Иэд-во УГНТУ, 1994, 132 с.

19. Установка измерительная групповая «Спутник-М». Руководство по эксплуатации. Тюмень, 2002.

20. Технологический учет жидкости в продукции нефтяных скважин. Групповая замерная установка «Дельта-М», НПО «НТЭС». www.nponts.ru.

21. Дробков В. П., Мельников В.И., Лабутин С.А. Методы и средства измерений влажности нефти. // Датчики и системы, 2002, №11, с. 24-27.

22. US Patent №6601458, United States, 73/861.04, G01F 1/74, 05.08.2003 Distributed sound speed measurements for multiphase flow measurement.

23. Международная заявка № WO 2003/073047, G01F 1/66, 04.09.2003, Appratus and method for measuring parameters of a mixture having liquid droplets suspendet in vapor flowing in a pipe.

24. Патент РФ № 2183012, G01F 1/74, 27.05.2002 / Томпсон Лэрд Берри (US), Способ измерения многофазного потока и устройство для его осуществления

25. Патент РФ № 2243510, G01F 1/74, 27.12.2004 / Аспелунн Эудун (NO), Видерее Top (US), Измерение многофазного потока в трубопроводе

26. Den Boer. Measurement of gas and liquid flowrates and watercut of multiphase mixtures of oil, water and gas // US patent № 5287752, 73/861.04, 22.02.1994

27. Патент 2138023 РФ, МКИ G01F1/74. Способ определения расхода компонентов многофазной среды. /В.П. Дробков (РФ), В.И. Мельников (РФ). №98103920/28; Заявлено 02.03.98; Опубликовано 20.09.99; Бюлл.-1999. № 26.

28. Патент 2126143 РФ, МКИ G01F1/74. Ультразвуковой расходомер компонентов многофазной среды. /В.И. Мельников (РФ), В.П. Дробков (РФ), А.В. Шустов (РФ). № 98104060/28; Заявлено 02.03.98; Опубликовано 10.02.99; Бюлл.- 1999. -№4.

29. Международная заявка № W02001/067050 Al, G01F 1/66, 13.09.2001 Simultaneous determination of multiphase flowrates and concentration/ Melnikov V. (RU) Drobkov V. (RU)

30. Roach G.I. et al. Muitphase Flowmeter for Oil Water a Gas by Dual Energy Gamma-Ray Transmission. Nucl. Geophys. vol. 8, N 3, p. 225, 1994

31. Д.Д.Булкин, Г.А.Соколов. Бесконтактный тепловой расходомер для измерения газожидкостных потоков // Датчики и системы, №12, 2008.-с.5-9.

32. Ю.И.Орехов, И.Н.Москалев, В.Е.Костюков и др. Способ определения покомпонентного расхода потока газожидкостной смеси продуктов газонефтедобычи в трубопроводе и устройство для его реализации //Патент на изобретение № RU 2164340 С2.

33. Патент США № 6732595, 73/861.27, G01F 1/66, 11.05.2004 Method of and system for determining the mass flow rate of a fluid flowing in a conduit

34. Международная заявка № WO 2002/039804, G01F 1/66, 23.05.2002 Antiparallel tag flow measurement system / Liu Yi; (US)

35. Патент США № 6293156 Bl, 73/861.26, G01F1/66. Coherent multi-path flow measurement system. / Chang Shen, Saul A. Jacobson 25.09.2001

36. Международная заявка № WO 2000/043736, G01F 1/66, 27.07.2000 Clamp-on gas flow meter / Shen Chang; (US), Jacobson Saul, A.; (US)

37. Бургасов B.C., Дунцев A.B., Лабутин С.А., Мельников В.И. Ультразвуковой измеритель волновых характеристик стекающих пленок жидкости // Приборы и техника эксперимента, 1992, № 2, с. 229-231.

38. Лабутин С.А., Мельников В.И. Ультразвуковые волноводные измерители параметров однофазных и двухфазных сред и потоков // Датчики и системы, 2001, №2, с. 54-62.

39. Вьюгин А. Б., Дунцев A.B., Лабутин С.А., Мельников В.И. Акустозондовая система исследования распределения газовой фазы в двухфазной среде // Датчики и системы, 2001, №11, с. 14-18.

40. Мельников В. И., Лабутин С.А., Шаронов Д. А. Анализ ультразвуковых импульсных методов измерений влажности нефти. // Датчики и системы, 2006, №1, с. 21-25.

41. Мельников В. П., Дробков В. П., Лабутин С. А., Шаронов Д. А. Исследования скорости распространения ультразвуковых волн в нефтеводяной эмульсии. // Датчики и системы, 2006, №2, с. 7-10.

42. Дунцев А.Н., Лабутин С.А., Мельников В.И., Хохлов В.Н. Допплеровсий ультразвуковой импульсный измеритель скорости и расхода двухфазных потоков жидкости // Приборы и техника эксперимента, 1994, № 2, с. 210.

43. Лабутин С.А., Пугин М.В. Исследование погрешностей ультразвукового измерителя волновых характеристик стекающих пленок жидкости // Межвуз. сб.: Радиоэлектронные и телекоммуникационные системы и устройства. -Н.Новгород, 1998, вып. 4, с. 66-74

44. Лабутин С.А., Мельников В.И. Ультразвуковые волноводные датчики и системы. Н.Новгород: НГТУ, 2001. - 86 с.

45. Кольцова, И.С. Распространение ультразвуковых волн в гетерогенных средах. / И.С. Кольцова СПб.: Изд-во С.-Петерб. ун-та, 2007. 247 с.

46. Муякшин, С.И. Ультразвуковая доплеровская система для измерения расхода газожидкостных потоков/ С.И. Муякшин, И.А. Владимиров // Датчики и системы, 2009, №8, с. 26-30.

47. И.А.Викторов. Физические основы применения волн Рэлея и Лэмба в технике М., Наука, 1966 -168 с.

48. Викторов И.А. Физические Звуковые поверхностные волны в твердых телах // М., Наука, 1981. -287 с.

49. Дж. Трэвис, Дж.Кринг. Lab VIEW для всех М.: ДМК Пресс, 2008-880 с.

50. Федосов В.П., Нестеренко А.К. Цифровая обработка сигналов в LabVIEW -М.: ДМК Пресс, 2007

51. Р.Довиак, Д.Зрнич. Допплеровские радиолокаторы и метеорологические наблюдения Л., Гидрометеоиздат, 1988 г.

52. Лабунцов, Д.А. Механика двухфазных систем. Учебное пособие для вузов //Д.А. Лабунцов, В.В. Ягов Москва: Издательство МЭИ, 2000. — 374 с: ил. ISBN 5-7046-0758-6

53. Исимару А. Распространение и рассеяние волн в случайно-неоднородных средах, т. 1 // А. Исимару М., Мир, 1981.-280 с.

54. Татарский В. И. Распространение волн в турбулентной атмосфере //В. И. Татарский М., Наука, 1967. - 548 с.

55. Рытов С. М. Введение в статистическую радиофизику, ч. 2 / С. М. Рытов, Ю. А. Кравцов, В. И. Татарский // М., Наука, 1978. 464 с.

56. Ануфриев И. Е. Самоучитель MatLab 5.3/б.х. // СПб., БХВ-Петербург, 2004.-736 с

57. Иглин С.П. Математические расчеты на базе Matlab. // "BHV-Санкт-Петербург" 2005г. 640 с.

58. Алексеев Е.Р., Чеснокова О.В. Решение задач вычислительной математики в пакетах Mathcad 12, MATLAB 7, Maple 9. // НТ Пресс, 2006,- 496 с.

59. Knowles J.K. J. Geophys. Res. // 1966. V. 71. N 22. P. 5480.

60. Coldren L.A. Appl. Phys. Lett. //1974. V. 25. N 7. P. 367.

61. Plesset M., Prosperetti A., Bubble dynamics and cavitation, Ann.Rev.Fluid Mech. 1977, 9, p.145-185

62. Keller J.B., Miksis M., Bubble oscillation of large amplitude. J. Acoust. Soc. Am., 1980, 68 pp. 628-633.

63. Мансфельд А.Д., Волков Г.П., Санин А. Г., Владимиров И.А. Импульсная ультразвуковая спектроскопия газовых пузырьков // «Акустический Журнал», 2010, том 56, № 3, с. 323-332

64. Слододник А., мл. Поверхностные акустические волны // Под ред. А. Олинера. М.: Мир, 1981. Гл. 6. 305 с. (Slobodnik A.J., jr. II Acoustic Surface Waves / Ed. A.A. Oliner. Topics in Appl. Phys. V. 24. Berlin et al.: Springer, 1978. Ch. 6).

65. Auld B.A. Acoustic Fields and Waves in Solids // New York: Wiley, 1973. V. 2. P. 278.

66. Шкадов В.Я. Гидродинамика и тепломассообмен с поверхностью раздела //Л.П.Холпанов, В.Я. Шкадов.-М.: Наука.- 1990.-271 c.-ISBN 5-02-001320-Х.

67. Зырынова А. В. Можаев В. Г. Резкое усиление локализации волн Релея под следом капли на анизотропной подложке // Письма в ЖТФ, Т. 35, вып. 10, 2009. стр.34-40.

68. Черданцев Андрей Викторович. Исследование свойств волн возмущения и капиллярной ряби в дисперсно-кольцевом течении : диссертация . кандидата физико-математических наук : 01.02.05.- Новосибирск, 2006.- 96 с.

69. Патент России №2375707, G01N29/02 10.12.2009 Способ контроля наличия газа в потоке жидкости (варианты) // Васильев Александр Алексеевич (RU), Шарипов Рафаил Кимович (RU), Краузе Александр Сергеевич (RU)

70. Урик Д.Р. Основы гидроакустики. Д.: Судостроение. 1978. 448 с.

71. Лямшев Л.М. Отражение звука тонкими пластинками и оболочками в жидкости. М.: Издво АН СССР, 1955.

72. Вучков И., Бояджиева Л., Солаков Е. Прикладной линейный регрессионный анализ. М.: Финансы и статистика, 1987. 239 с.