РЕНТГЕНОВСКИЙ АНАЛИЗАТОР РАСХОДА И СОСТАВА КОМПОНЕНТОВ ГАЗОЖИДКОСТНОГО ПОТОКА

Предлагаемая полезная модель относится к измерительной технике и может быть использована в нефтедобывающей промышленности для контроля дебита нефтяных скважин.
Известны анализаторы параметров газожидкостного потока, основанные на облучении контролируемого потока пучком гамма-излучения (см. патент RU 2301985, МПК G01N 9/24 и патенты РФ на полезную модель 35892, МПК G01N 9/24 и 37222, МПК G01N 23/00).
Известные анализаторы содержат источник гамма-излучения, например, радиоизотопный источник гамма-квантов, первичный коллиматор, предназначенный для формирования пучка гамма-излучения на выходе радиоизотопного источника, сцинтилляционный детектор гамма-излучения и фотоэлектронный преобразователь, предназначенные для преобразования пучка гамма-излучения, прошедшего через контролируемую среду, в электрический сигнал, вторичный коллиматор, предназначенный для формирования пучка гамма-излучения на входе сцинтилляционного детектора, и вычислитель, предназначенный для определения параметров газожидкостного потока по информации о степени поглощения гамма-излучения контролируемой средой.
Недостатками известных устройств являются низкая точность измерения и необходимость непрерывного экологического мониторинга.
Первый из указанных недостатков связан с высокой энергией гамма-квантов и становится особенно существенным при контроле потока смеси нефть-вода-газ, поскольку атомные номера наиболее тяжелых элементов, входящих в состав нефти и воды - углерода (12) и кислорода (16) - мало отличаются между собой по величине, а коэффициент поглощения гамма-излучения отдельными компонентами контролируемой среды зависит, в основном, от атомных номеров элементов, входящих в состав этих компонентов. В связи с этим различия в поглощении высокоэнергетического излучения водой и нефтью незначительно отличаются между собой, что затрудняет точное определение компонентного состава
контролируемой среды.
Второй недостаток известных устройств объясняется, во-первых, высокой энергией гамма-излучения и, во-вторых, - невозможностью приостановить излучение радиоизотопного источника в нерабочие периоды жизненного цикла устройства: при хранении, транспортировке и утилизации. Это обстоятельство существенно затрудняет эксплуатацию и утилизацию известных устройств и требует непрерывного экологического мониторинга.
От указанных недостатков свободен известный рентгеновский анализатор состава газожидкостного потока, основанный на облучении контролируемой среды пучком низкоэнергетического рентгеновского излучения (см. патент США 5689540, МПК G01 №23/22, G01 №23/06, G01 №23/087).
Этот анализатор содержит корпус, источник рентгеновского излучения - рентгеновскую трубку, источник питания рентгеновской трубки, первичные коллиматоры, предназначенные для формирования пучков рентгеновского излучения на выходе рентгеновской трубки, сцинтилляционные детекторы рентгеновского излучения, вторичные коллиматоры, предназначенные для формирования пучков рентгеновского излучения на входе сцинтилляционных детекторов, фотоэлектронные преобразователи сигналов сцинтилляционных детекторов, а также вычислитель, предназначенный для определения параметров газожидкостного потока по информации от фотоэлектронных преобразователей и для управления источником питания рентгеновской трубки.
Недостаток известного устройства заключается в невозможности определения скорости и расхода газожидкостного потока.
Этот недостаток вызван тем, что информативные сигналы о состоянии газожидкостного потока, вырабатываемые детекторами рентгеновского излучения, не связаны со скоростными параметрами потока, а зависят только от относительного содержания компонентов потока и состава каждого компонента.
Также известен рентгеновский анализатор компонентного состава газожидкостного потока, основанный на последовательном облучении трехкомпонентной среды, содержащей нефть, воду и газ двумя уровнями рентгеновского излучения: излучением высокого уровня и излучением низкого уровня (см. патент США №2007/0291898 А1, МПК G01 №23/06, G01F 1/66).
Этот анализатор содержит корпус, источник рентгеновского излучения - рентгеновскую трубку, управляемый источник питания рентгеновской трубки, измерительные и контрольные сцинтилляционные детекторы рентгеновского излучения, рентгенопрозрачные вставки, установленные в корпусе, первичные и вторичные коллиматоры, а также вычислитель объемного и массового содержания компонентов контролируемой среды.
Недостатком известного анализатора является невозможность определения покомпонентного расхода газожидкостного потока.
Указанные недостатки отсутствует у наиболее близкого к предлагаемой полезной модели по технической сущности и достигаемому результату известного рентгеновского анализатора расхода и состава компонентов газожидкостного потока (см. патент США 6,097,786, МПК G01 №23/06).
Данный анализатор принят за ближайший аналог (прототип) предлагаемой полезной модели.
В известном анализаторе использованы два различных метода измерения.
Для измерения компонентного состава газожидкостного потока использован метод рентгеновского зондирования контролируемой среды с помощью источника низкоэнергетического рентгеновского излучения - рентгеновской трубки. Информация о компонентном объемном составе контролируемой среды формируется в известном анализаторе путем измерения степени ослабления рентгеновского излучения при его прохождении через контролируемую среду, где оно поглощается и рассеивается. Измерение степени ослабления производится сцинтилляционными детекторами, преобразующими рентгеновское излучение в видимый свет. Оптические сигналы с выходов детекторов преобразуются в электрические сигналы с помощью фотоэлектронных преобразователей и поступают в электронный блок, в составе которого предусмотрены модули обработки упомянутых электрических сигналов, вычислитель параметров газожидкостного потока и модуль управления источником питания рентгеновской трубки.
В данном методе выходной сигнал каждого сцинтилляционного детектора зависит от энергии рентгеновского излучения, прошедшего через контролируемую среду, и является функцией объемного компонентного состава этой среды, что
позволяет вычислить относительное объемное содержание каждого компонента контролируемого потока.
Для измерения скорости газожидкостного потока использован кросс-корелляционный метод, в котором данные для вычисления скорости формируются в результате рентгеновского облучения контролируемого потока в двух его различных зонах, последовательно расположенных по направлению потока, и контроля рентгеновских излучений, прошедших через контролируемый поток в каждой из упомянутых зон, с помощью вторых и первых сцинтилляционных детекторов рентгеновского излучения, последовательно установленных по направлению газожидкостного потока. Метод позволяет вычислить скорости одного или нескольких компонентов газожидкостного потока на основе информации, вырабатываемой вторыми и первыми сцинтилляционными детекторами, путем обнаружения движущейся со скоростью потока локальной неоднородности его состава и измерения времени перемещения локальной неоднородности от вторых до первых сцинтилляционных детекторов.
В состав известного анализатора входят генератор низкоэнергетического рентгеновского излучения - рентгеновская трубка, источник питания рентгеновской трубки, первичные коллиматоры, предназначенные для формирования пучков рентгеновского излучения на выходе рентгеновской трубки, корпус, представляющий собой отрезок трубопровода, предназначенного для протекания контролируемой среды, рентгенопрозрачные вставки, установленные в стенке корпуса, первые и вторые сцинтилляционные детекторы, предназначенные для приема рентгеновского излучения, прошедшего через рентгенопрозрачные вставки и контролируемую среду, одноканальные фотоэлектронные преобразователи, например, ионизационные камеры, каждый из которых предназначен для преобразования оптического сигнала с выхода соответствующего сцинтилляционного детектора в электрический сигнал, вторичные коллиматоры, предназначенные для формирования пучков рентгеновского излучения на входах сцинтилляционных детекторов, и электронный блок, в составе которого предусмотрены вычислитель, модуль управления и модули обработки электрических сигналов, предназначенные для преобразования электрических сигналов с выходов фотоэлектронных преобразователей в измерительную
информацию о компонентном объемном составе и скорости контролируемого потока, а также для управления источником питания рентгеновской трубки. Сцинтилляционные детекторы и рентгеновская трубка известного анализатора установлены на противоположных сторонах корпуса, причем первые сцинтилляционные детекторы установлены в плоскости поперечного сечения корпуса, проходящей через ось излучения рентгеновской трубки, а вторые сцинтилляционные детекторы установлены перед первыми детекторами по направлению потока.
Известный анализатор позволяет вычислять относительное объемное содержание компонентов газожидкостного потока на основе измеренных параметров поглощения рентгеновского излучения контролируемой средой; кроме того, известный анализатор позволяет вычислять скоростные параметры газожидкостного потока с использованием кросс-корелляционного метода путем измерения скорости перемещения локальной неоднородности состава контролируемой среды.
К недостаткам известного анализатора относятся:
- существенная погрешность измерения скоростных параметров газожидкостного потока;
- низкая точность измерения массового содержания компонентов контролируемой среды;
- высокая инструментальная погрешность измерения параметров контролируемого потока, вызванная разбросом по величине и дрейфом во времени коэффициентов передачи отдельных одноканальных сцинтилляционных детекторов и отдельных одноканальных фотоэлектронных преобразователей.
Первый недостаток известного анализатора - существенная погрешность измерения скоростных параметров газожидкостного потока - вызван различиями в способе определения местоположения движущейся локальной неоднородности состава контролируемой среды первыми и вторыми сцинтилляционными детекторами. Первые детекторы воспринимают рентгеновское излучение, пересекающее корпус в плоскости, содержащей ось излучения рентгеновской трубки, и определяют местоположение локальной неоднородности потока по отношению именно к этой плоскости, ортогональной продольной оси корпуса.
Однако вторые детекторы, расположенные по потоку перед первыми, воспринимают рентгеновское излучение, пересекающее корпус под некоторым углом к плоскости его поперечного сечения, и, в силу этого, определяют местоположение той же локальной неоднородности по отношению не к ортогональной, а к некоторой наклонной плоскости. Это приводит к неоднозначности в определении как конфигурации, так и местоположения локальной неоднородности потока и не дает возможности однозначным образом фиксировать момент времени прохождения неоднородности мимо первых и мимо вторых сцинтилляционных детекторов, что приводит к существенной погрешности вычисления скорости потока кросс-корелляционным методом.
Следует упомянуть, что для снижения указанной погрешности в прототипе предусмотрен вариант исполнения с использованием дополнительной рентгеновской трубки, установленной в плоскости расположения вторых сцинтилляционных детекторов. Однако это решение, несмотря на существенное усложнение конструкции известного устройства, не позволяет существенно снизить упомянутую погрешность из-за неизбежного различия между собой спектров, интенсивностей и временных дрейфов излучения двух рентгеновских трубок.
Второй недостаток известного анализатора - низкая точность измерения массового содержания компонентов контролируемой среды - вызван тем, что степень поглощения контролируемой средой рентгеновского излучения неоднозначно связана с плотностью среды и не может достоверно использоваться в качестве информативного параметра для определения массового содержания ее компонентов.
Третий недостаток известного анализатора - высокая инструментальная погрешность измерения параметров контролируемого потока - вызван разбросом коэффициентов передачи отдельных одноканальных сцинтилляционных детекторов и отдельных одноканальных фотоэлектронных преобразователей как по номинальным значениям, так и по величинам временных и температурных дрейфов. Этот разброс приводит к неконтролируемым различиям между собой амплитуд сигналов на выходе как отдельных одноканальных сцинтилляционных детекторов, так и отдельных одноканальных фотоэлектронных преобразователей
известного анализатора и создает неустранимую дополнительную инструментальную погрешность измерения.
Задачей предлагаемой полезной модели и ее техническим результатом является повышение точности и достоверности измерения параметров газожидкостного потока, в том числе - массового содержания компонентов потока и скорости однородных потоков.
Для решения поставленной задачи изменены конструкция и состав элементов рентгеновского анализатора расхода и состава компонентов газожидкостного потока.
В состав анализатора входят корпус, в стенке которого последовательно по направлению потока установлены вторая и первая рентгенопрозрачные вставки, рентгеновская трубка с источником питания, первый и второй первичные коллиматоры с несколькими коллимирующими отверстиями в каждом, первый и второй вторичные коллиматоры с несколькими коллимирующими отверстиями в каждом, первые и вторые детекторы, первые и вторые фотоэлектронные преобразователи и электронный блок. В состав электронного блока входят вычислитель, первый и второй модули обработки и модуль управления, причем каждый из первых детекторов оптически соединен с соответствующим первым фотоэлектронным преобразователем, каждый из вторых детекторов оптически соединен с соответствующим вторым фотоэлектронным преобразователем, а вычислитель соединен со входом модуля управления, подключенного к источнику питания.
В заявленном устройстве новым по отношению к прототипу является то, что, согласно полезной модели, в его состав дополнительно введены первое и второе рентгеновские зеркала, третья рентгенопрозрачная вставка, установленная в стенке корпуса после первой рентгенопрозрачной вставки по направлению потока, третий первичный и третий вторичный коллиматоры с несколькими коллимирующими отверстиями в каждом, третьи детекторы, первые, вторые, третьи и четвертые многоканальные световоды, контрольные детекторы и вторичные излучатели, четвертые детекторы и ортогональные коллиматоры.
Ортогональные коллиматоры установлены под прямым углом к оси излучения рентгеновской трубки, а контрольные детекторы, установлены таким образом,
чтобы прямая линия, соединяющая центр каждого из них с центром излучения рентгеновской трубки, не пересекала корпус.
Кроме того, в состав анализатора введены датчик давления и измерительный преобразователь, контроллер режимов, а также первый и второй таймеры. При этом первые детекторы объединены в первый многоканальный детектор, а в качестве первых и вторых фотоэлектронных преобразователей применены, соответственно, первый и второй многоканальные фотоэлектронные преобразователи. Вычислитель дополнительно снабжен двумя многоканальными входами, первый и второй модули обработки снабжены каждый дополнительным выходом, причем первый модуль обработки снабжен также и дополнительным входом, перед каждым из детекторов первого многоканального детектора установлен соответствующий вторичный излучатель, перед каждым четвертым детектором установлен соответствующий ортогональный коллиматор, датчик давления установлен в стенке корпуса. Выход каждого из детекторов первого многоканального детектора и выход каждого из контрольных детекторов соединен с соответствующим каналом первого многоканального световода, выход каждого из вторых и выход каждого из третьих детекторов соединен, соответственно, с соответствующим каналом второго и соответствующим каналом третьего многоканальных световодов, выход каждого из четвертых детекторов соединен с соответствующим каналом четвертого многоканального световода. Выходы первого и выходы четвертого многоканальных световодов подключены каждый к соответствующим входам первого многоканального фотоэлектронного преобразователя, подсоединенного своими выходами к соответствующим входам первого модуля обработки с помощью многоканальной информационной связи, выходы второго и выходы третьего многоканальных световодов подключены каждый к соответствующим входам второго многоканального фотоэлектронного преобразователя, подключенного своими выходами к соответствующим входам второго модуля обработки с помощью многоканальной информационной связи. При этом контроллер режимов и измерительный преобразователь подключены каждый к соответствующему многоканальному входу вычислителя с помощью соответствующей многоканальной информационной связи. Дополнительный выход первого модуля обработки соединен со входом контроллера режимов,
дополнительный вход первого модуля обработки - с выходом первого таймера, а дополнительный вход второго модуля обработки - с выходом второго таймера, датчик давления подключен ко входу измерительного преобразователя. При этом вход первого и вход второго таймеров предназначены для подсоединения к внешним системам, многоканальный вход-выход которых предназначен для обмена информацией с вычислителем по двусторонней информационной связи. Ортогональный коллиматор содержит коллимирующие отверстия, глубина каждого из которых существенно больше диаметра, а оси параллельны между собой и ортогональны оси излучения рентгеновской трубки; вторичный излучатель представляет собой выполненную из тяжелого металла, например, гадолиния, трубку, установленную в отверстии первого вторичного коллиматора.
Устройство и работа предложенного анализатора поясняются Фиг.1-Фиг.6.
На Фиг.1 представлена функциональная схема анализатора, на Фиг.2 - поперечное сечение корпуса в плоскости, содержащей ось излучения рентгеновской трубки, на Фиг.3 - функциональная схема, поясняющая работу анализатора при измерении скорости, на Фиг.4 - графическая зависимость коэффициентов ослабления рентгеновского излучения от энергии излучения, на Фиг.5 - временная последовательность оптических импульсов, вызванных рентгеновскими фотонами с различной энергией, на Фиг.6 - графическая зависимость плотности потока рентгеновских фотонов от их энергии при нескольких значениях напряжения питания рентгеновской трубки.
На Фигурах введены следующие обозначения: 1 - корпус, 2, 3 и 4 - первая, вторая и третья рентгенопрозрачные вставки, соответственно, 5 - рентгеновская трубка, 6 и 7 - первое и второе зеркала, соответственно, 8, 9 и 10 - первый, второй и третий первичные коллиматоры, соответственно, в каждом из которых выполнены коллимирующие отверстия, 11 - первый многоканальный детектор, 12 и 13 - вторые и третие сцинтилляционные детекторы, соответственно, 14, 15 и 16 - первый, второй и третий вторичные коллиматоры, соответственно, в каждом из которых выполнены коллимирующие отверстия, 17 - четвертые сцинтилляционные детекторы, 18 - ортогональные коллиматоры, 19, 20, 21 и 22 - первый, второй, третий и четвертый световоды, соответственно, 23 - контрольные сцинтилляционные детекторы, 24 - вторичные излучатели, 25 - первый
фотоэлектронный преобразователь, 26 - первый модуль обработки, 27 - второй фотоэлектронный преобразователь, 28 - второй модуль обработки, 29 - вычислитель, 30 - модуль управления, 31 - электронный блок, 32 - контроллер режимов, 33 - источник питания, 34 - внешние системы, 35 - датчик давления, 36 - измерительный преобразователь, 37 - первый таймер, 38 - второй таймер.
Предложенный анализатор содержит корпус 1, представляющий собой отрезок трубопровода с фланцами на его концах, предназначенными для подсоединения корпуса 1 к внешней магистрали.
В стенке корпуса 1 установлены рентгенопрозрачные вставки: первая, вторая и третья вставки 2, 3, и 4, например, кольцеобразные, выполненные из бериллия.
Низкоэнергетическая рентгеновская трубка 5, установлена таким образом, что ее ось излучения направлена в сторону первой вставки 2 вдоль ее диаметра. Напротив второй и напротив третьей вставок 3 и 4, расположены, соответственно, второе и первое зеркала 7 и 6, каждое из которых установлено по отношению к плоскости поперечного сечения корпуса под некоторым углом.
На выходе рентгеновской трубки 5 установлены первый, второй и третий первичные коллиматоры 8, 9 и 10, соответственно, с несколькими веером расходящимися коллимирующими отверстиями в каждом. При этом второй первичный коллиматор 9 и третий первичный коллиматор 10 расположены перед первым и перед вторым зеркалами 6 и 7, соответственно. Коллимирующие отверстия в каждом из коллиматоров 9, 10 выполнены таким образом, что оси каждого из них направлены по отношению к плоскости поперечного сечения корпуса под некоторым углом. Угол ориентации зеркал 6 и 7 и угол ориентации коллимирующих отверстий коллиматоров 9 и 10 выбираются таким образом, чтобы отраженные зеркалами пучки рентгеновского излучения пересекали корпус в плоскостях расположения вставок 3 и 4, соответственно.
В состав анализатора также входят первый многоканальный детектор 11 и вторые и третие детекторы 12 и 13, соответственно, на входе каждого из которых установлены соответствующие первые, вторые и третие вторичные коллиматоры 14, 15 и 16, соответственно, в каждом из которых содержится несколько коллимирующих отверстий, ось каждого из которых направлена от центра излучения рентгеновской трубки 5 к центру соответствующего детектора, при этом
число упомянутых отверстий равно числу соответствующих детекторов. Кроме того, в состав анализатора входят два четвертых детектора 17, два ортогональных коллиматора 18, каждый из которых установлен перед одним из четвертых детекторов 17, а также первый, второй, третий и четвертый световоды 19, 20, 21 и 22, соответственно. В плоскости установки первого многоканального детектора 11 расположены также контрольные детекторы 23, а в каждом, кроме двух крайних, отверстии первого вторичного коллиматора 14 установлен один из вторичных излучателей 24, представляющий собой коллимирующую трубку, выполненную из тяжелого металла, характеристическая линия которого расположена в низкоэнергетической части жесткого рентгеновского диапазона, например, из гадолиния или золота.
Каждый контрольный детектор 23 установлен таким образом, чтобы прямая линия, соединяющая его центр с центром излучения рентгеновской трубки 5, не пересекала первую вставку 2, т.е. проходила снаружи корпуса 1, а каждый из ортогональных коллиматоров 18 расположен в плоскости установки детекторов первого многоканального детектора 11 таким образом, что оси его коллимирующих отверстий расположены под прямым углом к оси излучения рентгеновской трубки 5.
Выход каждого из детекторов первого многоканального детектора 11, выход каждого из контрольных детекторов 23 и выход каждого из четвертых детекторов 17 подсоединен к соответствующему каналу первого световода 19, выход каждого из вторых детекторов 12 подсоединен к соответствующему каналу второго световода 20, а выход каждого из третьих детекторов 13 подсоединен к соответствующему каналу третьего световода 21.
Первый световод 19 и четвертый световод 22 подключены каждый своими выходами к соответствующим входам первого фотоэлектронного преобразователя 25, выход которого соединен многоканальной информационной связью со входом первого модуля обработки 26, а второй и третий световоды 20 и 21, соответственно, подключены каждый своими выходами к соответствующим входам второго фотоэлектронного преобразователя 27, выход которого соединен многоканальной информационной связью со входом второго модуля обработки 28.
Выход первого и выход второго модулей обработки 26 и 28, соответственно, подключены каждый к соответствующему входу вычислителя 29, выход которого
соединен с модулем управления 30. Вычислитель 29, первый и второй модули обработки 26 и 28, соответственно, и модуль управления 30 входят в состав электронного блока 31.
Предложенный анализатор также содержит контроллер режимов 32, выход которого подключен с помощью многоканальной информационной связи к соответствующему многоканальному входу вычислителя 29, а вход соединен с дополнительным выходом первого модуля обработки 26, и источник питания 33 рентгеновской трубки 5, вход которого подключен к выходу модуля управления 30.
Вычислитель 29 содержит выход, подсоединенный ко входу модуля управления 30, и может быть снабжен двусторонней информационной связью для обмена информацией с внешними системами 34.
В стенке корпуса 1 установлен датчик давления 35 подсоединенный ко входу измерительного преобразователя 36, подключенного к соответствующему многоканальному входу вычислителя 29 с помощью многоканальной информационной связи.
В состав анализатора также входят первый таймер 37, выход которого подключен к дополнительному входу первого модуля обработки 26, и второй таймер 38, выход которого соединен с дополнительным входом второго модуля обработки 28; вход первого таймера 37 и вход второго таймера 38 предназначены каждый для подключения к внешним системам 34.
Предложенный рентгеновский анализатор расхода состава компонентов газожидкостного потока работает следующим образом.
При подаче в вычислитель 29 сигнала запуска, например, из внешних систем 34 по двусторонней информационной связи, с выхода вычислителя 29 на вход модуля управления 30 поступает команда включения источника питания 33, который включается и запитывает рентгеновскую трубку 5 напряжением, соответствующим номинальному режиму питания, заданному модулем управления 30. В результате на выходе рентгеновской трубки 5 возбуждается и пересекает контролируемую среду, находящуюся в корпусе 1, низкоэнергетическое рентгеновское излучение, пучки которого формируются в каждом из коллимирующих отверстий первого, второго и третьего первичных коллиматоров 8, 9 и 10 соответственно.
Упомянутые первичные коллиматоры формируют три группы рентгеновских пучков, веерообразно расходящихся от центра излучения рентгеновской трубки 5 по направлениям к соответствующим отверстиям соответствующего вторичного коллиматора.
Первая группа веерообразно расходящихся пучков формируется первым первичным коллиматором 8 и лежит в плоскости поперечного сечения В-В корпуса 1, содержащем первую вставку 2 и ось излучения рентгеновской трубки 5; основная часть пучков этой группы, за исключением двух крайних пучков, пересекает корпус 1 через первую вставку 2, а крайние пучки не пересекают корпус 1 и проходят по воздуху у противоположных сторон корпуса 1 (см. Фиг.2).
Вторая группа веерообразно расходящихся пучков формируется вторым первичным коллиматором 9 и направлена в сторону зеркала 7. Отраженные зеркалом 7 пучки пересекают корпус 1 в сечении А-А, содержащем вторую вставку 3.
Третья группа веерообразно расходящихся пучков   формируется третьим первичным коллиматором 10 и направлена в сторону зеркала 6. Отраженные зеркалом 6 пучки пересекают корпус 1 в сечении С-С, содержащем третью вставку 4 (см. Фиг.3).
Каждый из двух крайних рентгеновских пучков первой группы, сформированных в крайних коллимирующих отверстиях первого первичного коллиматора 8 в сечении В-В корпуса 1, попадает, не пересекая корпус 1, в одно из крайних отверстий первого вторичного коллиматора 14, и, после формирования в этом отверстии, падает на соответствующий ему контрольный детектор 23. Каждый из сформированных первым первичным коллиматором 8 рентгеновских пучков, лежащих в сечении В-В корпуса 1, кроме двух крайних пучков, пересекает первую вставку 2 и контролируемую среду, заполняющую корпус 1, достигает соответствующего ему вторичного излучателя 24, установленного в отверстии первого вторичного коллиматора 14 и, после формирования в коллимирующем отверстии этого излучателя, падает на один из детекторов многоканального первого детектора 11.
В каждом из контрольных детекторов 23 и в каждом из детекторов многоканального первого детектора 11 формируется оптический сигнал, соответствующий интенсивности принятого рентгеновского излучения. Каждый из
сформированных оптических сигналов подается в соответствующий канал первого световода 19 и с соответствующего выхода последнего поступает на соответствующий вход первого фотоэлектронного преобразователя 25.
Одновременная передача по каналам первого световода 19 нескольких оптических сигналов, выработанных несколькими контрольными детекторами 23 и несколькими детекторами многоканального первого детектора 11, дает возможность использовать для фотоэлектронного преобразования этих сигналов единый многоканальный фотоэлектронный преобразователь - первый фотоэлектронный преобразователь 25.
Такое техническое решение позволяет свести к минимуму разброс коэффициентов преобразования и величин температурных и временных дрейфов нескольких отдельных детекторов и нескольких отдельных фотоэлектронных преобразователей и тем самым минимизировать соответствующую инструментальную погрешность.
Каждый из рентгеновских пучков, пересекающих контролируемую среду в плоскости расположения первой вставки 2, частично поглощается и рассеивается этой средой, что приводит к уменьшению плотности потока рентгеновских фотонов, падающих на каждый из детекторов первого многоканального детектора 11, и дает возможность оценивать компонентный состав поглощающей среды вдоль линии пересечения по степени уменьшения плотности потока фотонов в сравнении с исходной плотностью этого потока.
1
2
3
1
max1
2
mах2
3
mах3
1
2
3
1
2
Для получения более детальных данных об объемном компонентном составе следует последовательно облучать контролируемую среду пучками низкоэнергетического рентгеновского излучения, различающимися между собой уровнем энергии, при этом каждый уровень энергии рентгеновского излучения задается соответствующим ему напряжением питания рентгеновской трубки 5. Если, например, необходимо получить, как это показано на Фиг.6, три уровня энергии рентгеновского излучения, следует использовать три значения напряжения питания U, U и U, обеспечивающих уровни энергии E(I), Е(I) и Е(I), соответственно, отвечающие максимальным значениям плотностей потоков фотонов Imax, Imax и Imах на выходе рентгеновской трубки 5. Для этого рентгеновская трубка 5 последовательно запитывается напряжениями U, U
3
и U от источника питания 33 в соответствии с переключающими сигналами, поступающими на вход этого источника из модуля управления 30 при получении им соответствующих команд из вычислителя 29.
При облучении контролируемой среды рентгеновским излучением различной энергии в жестком рентгеновском диапазоне 30-100 КэВ изменяется характер поглощения рентгеновских фотонов каждым из компонентов контролируемой среды.
Это объясняется тем, что степень ослабления рентгеновского излучения компонентом контролируемой среды зависит не только от состава этого компонента, но и от энергии излучения, то есть является функцией энергии рентгеновских фотонов, генерируемых рентгеновской трубкой 5 (см. Фиг.4):
μ - линейный коэффициент ослабления компонентом контролируемой среды рентгеновского излучения с энергией Е;
Е - энергия рентгеновского излучения.
В случае фиксированного уровня энергии рентгеновского излучения линейный коэффициент ослабления падающего пучка конкретным компонентом контролируемой среды является постоянной величиной. Так, например, если выбрать конкретный компонент и зафиксировать определенный уровень энергии излучения, то значение линейного коэффициента ослабления μ будет постоянным.
Таким образом, зависимость степени ослабления плотности потока фотонов заданной энергии конкретным компонентом определяется эффективным линейным размером этого компонента вдоль направления соответствующего рентгеновского пучка (далее по тексту - эффективной толщиной компонента):
Y(x) - плотность потока фотонов на входе детектора многоканального первого детектора 11 при заполнении корпуса 1 однокомпонентной контролируемой средой;
I - плотность потока фотонов на входе детектора многоканального первого детектора 11 при отсутствии в корпусе 1 контролируемой среды;
e - основание натуральных логарифмов;
μ - линейный коэффициент ослабления рентгеновского излучения контролируемой средой;
x - эффективная толщина контролируемой среды.
1
2
3
Обозначим компоненты контролируемой среды следующим образом: K - нефть, К - вода и К - газ.
При наличии на пути падающего рентгеновского пучка всех трех компонентов: нефти, воды и газа выражение (2) принимает вид функционального уравнения:
m1
1
Y(t) - текущее значение плотности потока фотонов на входе детектора с порядковым номером m первого многоканального детектора 11 при напряжении U питания рентгеновской трубки 5 и при наличии в корпусе 1 трехкомпонентной контролируемой среды;
1
2
3
1
2
3
x, x и x - эффективные толщины компонентов K, К и К, контролируемой среды, соответственно;
1
q
2
q
3
q
q
1
2
3
μ(E), μ(E) и μ(E) - линейные коэффициенты ослабления рентгеновских излучений с фиксированной энергией E каждым из компонентов K, К и К,
m1
1
I(t) - текущее значение плотности потока фотонов на входе детектора с порядковым номером m, входящего в состав первого многоканального детектора 11, при напряжении U питания рентгеновской трубки 5 и при отсутствии в корпусе 1 контролируемой среды;
t - текущее время.
q
1
q
2
q
3
i
1
i
2
i
3
1
2
3
1
q
2
q
3
q
Последовательно подставляя в уравнение (3) значения линейных коэффициентов ослабления, соответствующие нескольким, например, трем, значениям энергии излучения E=E, E=Е, E=Е, а именно, значения μ(E), μ(Е) и μ(Е), где i=1, 2, 3, получим для каждого рентгеновского пучка, лежащего в плоскости расположения первой вставки 2, три системы функциональных уравнений вида (3) с тремя численными неизвестными x, x, x и тремя функциональными неизвестными μ(E), μ(E), μ(E).
m1
В случае, когда текущее значение параметра I(t) известно, то есть когда оно непрерывно измеряется, когда конкретное значение энергии излучения
q
1
q
2
q
3
1
q
2
q
3
q
1
2
3
1
2
3
1
2
3
E=E, E=E, E=Е зафиксировано, а также когда определены и занесены в память первого модуля обработки 26 экспериментальные зависимости линейных коэффициентов ослабления от энергии излучения μ(E), μ(Е), μ(Е), решение систем функциональных уравнений вида (3) в вычислителе 29 позволяет получить значения эффективных толщин x, x, и x каждого из компонентов K, К и К, соответственно, вдоль линии пересечения контролируемой среды одним из падающих рентгеновских пучков и тем самым определить объемную долю каждого из этих компонентов в поперечном сечении контролируемого потока вдоль направления упомянутого пучка. После определения эффективных толщин x, x, и x нефти, воды и газа, соответственно, вдоль каждого из рентгеновских пучков, падающих на один из m детекторов первого многоканального детектора 11, в вычислителе 29 производится контроль достоверности полученных значений эффективных толщин, для чего используется контрольное условие:
ℓ - длина хорды, стягивающей дугу окружности радиуса R вдоль соответствующего рентгеновского пучка;
R - внутренний радиус корпуса 1.
Для центрального рентгеновского пучка, направленного вдоль оси излучения рентгеновской трубки 5, длина ℓ упомянутой хорды составляет ℓ=2R, для всех других пучков она меньше указанного значения: ℓ<2R.
1
2
3
Операции решения системы функциональных уравнений вида (3) относительно искомых величин x, х и х выполняются в вычислителе 29.
m1
Параметр I(t) упомянутых функциональных уравнений характеризует текущее значение плотности потока фотонов на входе соответствующего детектора первого многоканального детектора 11 при отсутствии в корпусе 1 контролируемой среды.
m1
Поскольку при эксплуатации предложенного анализатора контролируемая среда всегда заполняет корпус 1, обеспечить состояние «при отсутствии контролируемой среды» весьма затруднительно. В тоже время, в процессе эксплуатации предложенного анализатора текущее значение параметра I(t) может существенно изменяться по величине относительно своего номинального
значения, измеренного и введенного в память первого модуля обработки 26 при юстировке предложенного анализатора.
m1
m1
Изменение текущего значения параметра I(t) может происходить, например, в связи с износом катода рентгеновской трубки 5, изменением напряжения ее питания, смещением пространственного положения рентгеновской трубки 5 относительно начального положения и т.п. Для того, чтобы уменьшить влияние указанных факторов на точность измерения компонентного состава, определение текущего значения параметра I(t) в процессе эксплуатации производится не прямым методом - по информации, формируемой детекторами первого многоканального детектора 11, а косвенно - на основании информации, вырабатываемой контрольными детекторами 23. При этом принимается во внимание, что коэффициенты передачи отдельных детекторов, входящих в состав первого многоканального детектора 11, и коэффициенты передачи контрольных детекторов 23 в случае, когда их сцинтиллирующие кристаллы принадлежат одной партии изготовления, практически не отличаются между собою по величине и одинаково изменяются как во времени, так и под действием внешних дестабилизирующих факторов.
В этом случае
, откуда
m1
1
I(t) - текущее значение сигнала на выходе детектора с порядковым номером m, входящего в состав первого многоканального детектора 11, при напряжении U питания рентгеновской трубки 5, пропорциональное текущему значению плотности потока фотонов на входе этого детектора при отсутствии в корпусе 1 контролируемой среды;
m1
I - начальное значение упомянутого сигнала;
n1
1
I(t) - текущее значение сигнала на выходе одного из двух контрольных детекторов 23 при напряжении U питания рентгеновской трубки 5, пропорциональное текущему значению плотности потока фотонов на входе этого детектора;
n1
I - начальное значение упомянутого сигнала;
t - текущее время.
n1
n1
Численное значение каждого из сигналов I(t) и I в выражении (5) определяется в первом модуле обработки 26 как среднеарифметическое значение двух оптических сигналов, каждый из которых вырабатывается одним из контрольных детекторов 23, передается по соответствующему каналу световода 19 на вход первого фотоэлектронного преобразователя 25 и, после оптоэлектронного преобразования в этом преобразователе, поступает в электронной форме по многоканальной информационной связи на вход первого модуля обработки 26.
n1
n1
Поскольку рентгеновское излучение, падающее на каждый из контрольных детекторов 23, пересекает только воздушную среду, поглощением и рассеянием в которой при энергиях излучения более 30 КэВ можно пренебречь, оптические сигналы, формируемые этими детекторами, не зависят от состава пропускающей среды, определяются только характеристиками излучения рентгеновской трубки 5 и могут использоваться в качестве контрольных сигналов I и I(t).
n1
Начальное значение I контрольного сигнала измеряется в процессе юстировки предложенного анализатора и вносится в память первого модуля обработки 26.
11
21
m1
1
12
22
m2
13
23
m3
2
3
Кроме того, в процессе юстировки измеряются и вносятся в память первого модуля обработки 26 начальные значения сигналов I, I, ... I, где m=1, 2 ..., 5, каждое из которых соответствует одному из детекторов первого многоканального детектора 11 при напряжении U питания рентгеновской трубки 5. В память упомянутого модуля также вносятся значения сигналов I, I, ..., I и I, I, ..., I, соответствующие напряжениям питания U и U рентгеновской трубки 5.
n1
При работе анализатора в первый модуль обработки 26 непрерывно поступает информация о текущем значении I(t) контрольного сигнала на выходе каждого из контрольных детекторов 23.
n1
m1
m1
На основе этой информации с использованием хранящихся в памяти первого модуля обработки 26 данных о значениях I и I вычисляется, в соответствии с выражением (5), и передается в вычислитель 29 текущее значение параметра I(t).
Из (3) и (5) имеем:
каждое из обозначений соответствует обозначениям, принятым в выражениях (3) и
(5).
1
q
2
q
3
q
q
Параметры μ(E), μ(E) и μ(Е) уравнения (6) представляют собой функциональные неизвестные, значения которых определяются в первом модуле обработки 26 на основе введенных в его память экспериментальных данных о взаимосвязи значения линейного коэффициента ослабления μ рентгеновского излучения с энергией E, отмасштабированной по реперному значению энергии.
1
2
3
1
2
3
Экспериментально полученные графические зависимости линейных коэффициентов ослабления рентгеновского излучения μ, μ, и μ каждым из компонентов контролируемой среды K, К, и К - нефтью, водой и газом, соответственно, от энергии Е рентгеновского излучения представлены на Фиг.4.
Как видно из этих графиков, упомянутые зависимости существенно нелинейны, поэтому для корректного решения функционального уравнения (6) необходимо учитывать не только интегральную характеристику рентгеновского излучения - плотность потока его фотонов I, но и дифференциальную характеристику - распределение рентгеновских фотонов по энергиям Е.
2
1
Диаграмма распределения во времени рентгеновских фотонов с различными энергиями, принятых одним из детекторов первого многоканального детектора 11 за некоторый фиксированный промежуток времени Δt=t-t, приведена на Фиг.5.
Данная диаграмма представляет собой временную последовательность немасштабированных по амплитуде оптических импульсов различной интенсивности J, каждый из которых сформирован одним из детекторов первого многоканального детектора 11 при падении на его вход соответствующего рентгеновского фотона; интенсивность J такого импульса пропорциональна энергии соответствующего фотона в жестком диапазоне рентгеновского излучения.
2
1
1
2
3
При достаточно большом числе S импульсов, выработанных каждым из детекторов первого многоканального детектора 11 за промежуток времени Δt=t-t, в первом модуле обработки 26 моделируется представленная на Фиг.6 дифференциальная характеристика рентгеновского излучения - зависимость плотности потока фотонов I от энергии фотонов Е при каждом из напряжений питания U, U и U из рентгеновской трубки 5. Промежуток времени Δt, в течение которого формируется количество S упомянутых импульсов, обычно называется временем экспозиции.
Для количественного масштабирования качественных зависимостей, представленных на Фиг.5 и на Фиг.6, необходимо сформировать реперное значение энергии, используя для этого достоверно определяемые параметры рентгеновских фотонов. Реперное значение энергии в предложенном анализаторе формируется с помощью вторичных излучателей 24.
Каждый вторичный излучатель 24 представляет собой установленную в отверстии первого вторичного коллиматора 14 коллимирующую трубку, изготовленную из тяжелого металла с ярко выраженной характеристической линией спектра излучения, лежащей в диапазоне энергий 30-100 КэВ, например, из гадолиния или платины.
При прохождении рентгеновского пучка сквозь упомянутую трубку фотоны пучка могут возбуждать флуоресценцию атомов гадолиния, расположенных в приповерхностном слое коллимирующего отверстия трубки, сопровождаемую возникновением вторичного характеристического излучения. Возбуждение флуоресценции происходит только в том случае, когда энергии ℏω возбуждающих фотонов совпадают с энергией характеристического излучения или превышают ее:
ℏ - постоянная Планка;
ω - круговая частота возбуждающих фотонов;
о
E - энергия характеристического излучения.
Спектр характеристического излучения является линейчатым; характеристической линии гадолиния в этом спектре соответствует энергия излучения
Gd
E=42,996 КэВ.
Если энергия возбуждающих фотонов существенно превосходит энергию характеристического возбуждения атомов гадолиния:
о
то плотность потока характеристических фотонов с уровнем энергии E существенно превосходит плотности потока фотонов с другими значениями энергии.
Поэтому при наложении на излучение рентгеновской трубки 5, характеризуемое непрерывным спектром, характеристического излучения, имеющего линейчатый спектр, существенно изменяется характер распределения I(E) плотности потока
рентгеновских фотонов по энергиям на входе каждого из детекторов первого многоканального детектора 11 (см. Фиг.6).
При отсутствии характеристического излучения распределение I(E) представляет собой в жестком рентгеновском диапазоне убывающую кривую, представленную на Фиг.6 сплошной линией.
Gd
Gd
Gd
Gd
max1
o
При наличии характеристического излучения с энергией фотонов Е=Е в упомянутом распределении возникает представленный на Фиг.6 штрих-пунктирной линией локальный максимум плотности потока фотонов I=I(E), где E - энергия возбуждения атомов гадолиния. Уровень этого максимума существенно превосходит наибольший уровень плотности потока фотонов I при всех других значениях энергии, отличающихся от значения Е:
Gd
max1
I>I.
Gd
o
o
mах
о
Gd
Gd
Возникновение локального максимума плотности потока фотонов I при характеристическом значении Е энергии объясняется тем, что значительное число первичных рентгеновских фотонов, распределенных в широком диапазоне энергии от Е до Е, поглощается атомами гадолиния, порождающими поток такого же числа вторичных флуоресцентных фотонов, сосредоточенных в весьма узкой энергетической области вблизи характеристического значения энергии E=Е. Поскольку все вторичные фотоны пересекают только воздушную среду, поглощением и рассеянием в которой можно пренебречь, практически все они достигают детекторов первого многоканального детектора 11 и обеспечивают на их входах высокую локальную плотность потока фотонов с характеристическим значением энергии E.
Gd
Gd
Каждому фотону с энергией Е, поглощенному детектором первого многоканального детектора 11, отвечает формируемый на выходе этого детектора световой импульс интенсивностью J (см. Фиг.5), передаваемый по соответствующему каналу первого световода 19 на вход первого преобразователя 25, где он преобразуется в электрический импульс, поступающий по многоканальной информационной связи в первый модуль обработки 26.
Поскольку интенсивность оптического импульса и, соответственно, амплитуда соответствующего ему электрического импульса пропорциональны энергии поглощенного фотона, возникает возможность обнаружения в первом
Gd
модуле обработки 26 наличия фотонов с характеристическим значением энергии E по частоте появления электрических импульсов одной и той же амплитуды. Способ обнаружения упомянутых фотонов поясняется временной последовательностью J(t) оптических импульсов, представленной на Фиг.5.
Gd
Gd
Gd
Gd
Анализ этой временной последовательности показывает, что наиболее часто, по сравнению с другими, встречаются импульсы интенсивности J, заштрихованные на Фиг.5. Эти импульсы соответствуют характеристическим фотонам с энергией Е и создают на кривой I(Е) (Фиг.6) локальный максимум плотности потока фотонов I в точке Е=Е.
Gd
Gd
Gd
Наличие локального максимума I позволяет объективно выделить из общего числа фотонов различной энергии группу фотонов с энергией Е и тем самым обозначить на оси энергий Е графика I(Е), Фиг.6 точную масштабную отметку E=42,996 КэВ, позволяющую достоверно отмасштабировать весь график I(Е) по значениям энергий.
1
q
2
q
3
q
1
2
q
Исходя из изложенного, для решения в вычислителе 29 уравнения (6) предварительно вычисляются в первом модуле обработки 26 численные значения каждого из параметров μ(E), μ(E) и μ(E) для нескольких фиксированных значений энергии Е=E, E, ..., E, ..., каждое из которых соответствует узкой энергетической полосе шириной 2ΔЕ, содержащей не менее числа s фотонов, а именно, численные значения
1
2
3
μ, μ и μ - линейные коэффициенты ослабления излучения, соответственно, нефтью, водой и газом;
q
E - фиксированное значение энергии излучения;
ΔE - разность двух последовательных фиксированных значений энергии излучения.
q
1
2
3
Количественные значения параметров μ(E) вычисляются в первом модуле обработки 26 на основе заложенных в его память численных данных, соответствующих приведенным на Фиг.4 экспериментальным зависимостям каждого из коэффициентов μ, μ и μ от энергии излучения Е. При этом, для уменьшения погрешности вычисления, при выборе численных данных выбираются возможно более узкие энергетические полосы шириной 2ΔЕ.
Для надежного формирования каждого значения плотности потока фотонов, усредненного на энергетической полосе шириной 2ΔЕ, необходимо получить за время экспозиции Δt не менее числа s=100 фотонов.
Численное значение времени экспозиции Δt выбирается и вводится в память первого таймера 49, исходя из двух противоречивых условий.
q
Во-первых, оно должно быть, по возможности, большим, для того чтобы обеспечить поглощение каждым детектором многоканального детектора 11 достаточно большого числа S рентгеновских фотонов, необходимого для достоверного вычисления в первом модуле обработки 26 усредненных значений плотности потоков фотонов в каждой из энергетических полос E±ΔЕ в диапазоне энергий от 30 до 100 КэВ.
Во-вторых, оно должно быть, по возможности, малым, для того, чтобы компонентный состав контролируемой среды, непрерывно протекающей сквозь подконтрольную область корпуса 1, не успел существенно измениться за время экспозиции.
В соответствии с указанными условиями значение времени экспозиции устанавливается, исходя из максимальной скорости контролируемого потока и эффективного размера подконтрольной области:
Δt - время экспозиции;
а - эффективный размер подконтрольной области, равный размеру первой вставки 2 по направлению контролируемого потока;
m
W - наибольшее значение продольной составляющей скорости контролируемого потока.
1
q
2
q
3
q
Для обеспечения требуемой точности измерения параметров контролируемого потока необходимо определить минимально возможное время экспозиции, достаточное для получения надежной информации о компонентном составе контролируемого потока. В свою очередь, это время зависит от числа S рентгеновских фотонов с энергией от 30 до 100 КэВ, принятых каждым из детекторов первого многоканального детектора 11. Поскольку, в соответствии с (9), численные значения коэффициентов μ(E), μ(E) и μ(E), характеризующих компонентный состав потока, вычисляются в узких энергетических полосах
q
q
шириной 2ΔЕ, число S определяется как сумма чисел s рентгеновских фотонов с энергией от (E-ΔЕ) до (E+ΔЕ), принятых каждым из детекторов первого многоканального детектора 11 за минимально возможное время экспозиции, т.е.:
S - число рентгеновских фотонов, падающих на каждый детектор первого многоканального детектора 11 за время экспозиции Δt;
q=1, 2, 3, ... - натуральное число, равное порядковому номеру энергетической полосы;
q
s - число фотонов, энергии которых находятся в энергетической полосе шириной 2ΔЕ.
Практически, с учетом числа s=100,
Указанное в условии (12) число S рентгеновских фотонов должно приниматься каждым детектором первого многоканального детектора 11 за время экспозиции Δt, наименьшее численное значение которого при скорости контролируемой среды Wm≤10 м/с и продольном размере первой вставки а≤5 мм, в соответствии с (10), задается условием
min
5
Δt - минимально возможное время экспозиции, обеспечивающее прием каждым из детекторов первого многоканального детектора 11 числа S≥5·10 рентгеновских фотонов.
min
Приведенные в неравенствах (12) и (13) численные значения параметров S и Δt наряду с числом m детекторов первого многоканального детектора 11, полностью определяют параметры генерации рентгеновской трубки 5, и, кроме того, характеризуют требуемую информационную производительность процессоров, входящих в состав первого модуля обработки 26 и вычислителя 29.
3
0
0
3
1
2
При определении состава контролируемой среды с высоким газосодержанием в вычислителе 29 необходимо учитывать текущие значения P(t) давления и T(t) температуры, поскольку зависимость от этих параметров линейного коэффициента μ ослабления рентгеновского излучения газом при высоких давлениях в корпусе 1 весьма существенна (см. Фиг.4). Так, например, если при нормальных климатических условиях: Р=1 ат, Т=20°С величина линейного коэффициента μ для газа существенно меньше величин линейных коэффициентов μ и μ для нефти и воды, соответственно, то
3
1
2
при значениях давления P(t)=10-100 ат величина коэффициента μ становится соизмеримой с величинами коэффициентов μ и μ.
3
q
3
q
Поэтому при давлениях P(t)>10 ат в выражении (6) вместо коэффициента μ(E), соответствующего нормальным климатическим условиям, следует использовать коэффициент μ (E,P,T), соответствующий фактическим значениям давления и температуры контролируемой среды в корпусе 1:
3
q
все обозначения, кроме μ(E,P,T), соответствуют обозначениям, использованным в выражении (6);
3
q
3
q
μ(E,P,T) - линейный коэффициент ослабления рентгеновского излучения компонентом К (газом) при энергии излучения E, давлении контролируемой среды Р и температуре контролируемой среды Т.
3
q
min
2
1
1
2
min
3
q
3
Значения линейного коэффициента μ(E,P,T) вычисляются в вычислителе 29 за промежуток времени Δt=t-t, причем сигналы о моменте времени t начала вычислений и моменте времени t окончания вычислений поступают на дополнительный вход первого модуля обработки 26 с выхода первого таймера 38, в память которого из внешних систем 34 введено численное значение параметра Δt, соответствующее неравенству (13), и передаются в вычислитель 29 с выхода первого модуля обработки 26. При этих вычислениях используются введенные в память вычислителя 29 из внешних систем 34 экспериментально полученные численные данные о зависимости величины линейного коэффициента μ от энергии излучения E при нескольких фиксированных значениях давления Р и температуры Т (см. Фиг.4), справочные данные о значениях коэффициента сжимаемости β, а также измерительная информация о текущих значениях давления P(t) и температуры P(t) контролируемой среды, сформированная датчиком давления 35, преобразованная в цифровую форму в   измерительном преобразователе 36 и переданная в вычислитель 29 по многоканальной информационной связи с выхода упомянутого измерительного преобразователя.
1
2
3
Решение функционального уравнения (14) для каждого из трех значений U, U, U напряжения питания рентгеновской трубки 5 и для каждого из m=1, 2, ... 5 рентгеновских пучков, принимаемых детекторами первого многоканального
min
детектора 11, производится в вычислителе 29 за время экспозиции, Δt, соответствующее неравенству (13).
1
2
3
Совокупность решений функциональных уравнений (14) при различных значениях напряжения U=U, U, U питания рентгеновской трубки 5 для каждого из m=1, 2, ... 5 рентгеновских пучков позволяет получить в вычислителе 29 и передать по соответствующей двусторонней многоканальной информационной связи во внешние системы 34 детальную количественную информацию об относительном содержании в единицах объема каждого из трех компонентов контролируемой среды: нефти, воды и газа по всему поперечному сечению контролируемого потока:
1
1
V(t) - текущее значение объемной доли компонента K (нефти);
2
2
V(t) - текущее значение объемной доли компонента К (воды);
3
3
V(t) - текущее значение объемной доли компонента К (газа);
t - текущее время.
Для получения данных о компонентном составе контролируемой среды в единицах массы необходимо использовать дополнительную информацию, которая вырабатывается четвертыми детекторами 17.
Оптические сигналы, снимаемые с выходов четвертых детекторов 17, вырабатываются этими детекторами под воздействием рассеянного рентгеновского излучения, представляющего собой вторичное излучение, возникающее в результате упругого рассеяния контролируемой средой первичного излучения рентгеновской трубки 5. Рассеянное рентгеновское излучение направлено, в основном, под большими углами к направлению первичного излучения, близкими к прямым углам. Поэтому на каждый из четвертых детекторов 17, установленных в плоскости расположения первой вставки 2 на линии, нормальной оси излучения рентгеновской трубки 5, падает, в основном, рассеянное излучение. Для полного исключения возможности попадания на вход четвертых детекторов 17 первичного излучения рентгеновской трубки 5 перед каждым из них установлен ортогональный коллиматор 18, содержащей несколько коллимирующих отверстий, оси которых параллельны между собой и направлены под прямым углом к оси излучения рентгеновской трубки 5, а глубина каждого отверстия существенно
превышает его диаметр:
b/d≫1, где
b - глубина коллимирующего отверстия;
d - диаметр коллимирующего отверстия.
Поскольку интенсивность рассеянного излучения в жестком рентгеновском диапазоне энергий до 100 КэВ определяется, в основном, плотностью рассеивающей среды и слабо зависит от интенсивности первичного излучения рентгеновской трубки 5, оптический сигнал, формируемый каждым из двух четвертых детекторов 17, пропорционален текущему значению средней плотности ρ(t) контролируемой среды.
После преобразования в электрический цифровой код, упомянутый сигнал поступает в вычислитель 29, где на основании принятой информации о текущем значении средней плотности ρ(t) контролируемой среды и ранее вычисленных данных (15) о текущих  значениях объемных долей каждого из ее компонентов определяются текущие значения относительного содержания каждого из компонентов газожидкостного потока в единицах массы:
1
1
M(t) - текущее значение массовой доли компонента K (нефти);
2
2
М(t) - текущее значение массовой доли компонента К (воды);
3
3
М(t) - текущее значение массовой доли компонента К (газа);
t - текущее время.
M1
M2
M3
Текущие значения (16) массового содержания компонентов необходимы для определения в вычислителе 29 текущих значений покомпонентного массового расхода Q(t), Q(t) и Q(t) каждого из трех компонентов газожидкостного потока - нефти, воды и газа, соответственно.
1
Так как, помимо текущего значения массовой доли нефти M(t), дебит нефти традиционно принято оценивать также и в объемных единицах (баррелях), в вычислителе 29 определяется номинальная объемная доля нефти, т.е. объемная доля, приведенная к нормальным климатическим условиям.
Для вычисления текущего значения номинальной объемной доли нефти в вычислителе 29 используется дополнительная информация о текущих значениях давления P(t) контролируемой среды, сформированная датчиком давления 35,
переданная с его выхода на вход измерительного преобразователя 36 и с выхода последнего по многоканальной информационной связи - на соответствующий многоканальный вход вычислителя 29.
1
0
0
1
В вычислителе 29 на основе принятых сигналов о текущих значениях P(t) производится приведение ранее вычисленной объемной доли нефти V(t) к нормальным климатическим условиям: давлению Р=1 ат и температуре Т=20°С с учетом хранящихся в памяти вычислителя 29 справочных данных о значениях коэффициента сжимаемости β нефти. В результате этой операции в вычислителе 29 определяется текущее значение номинальной объемной доли нефти:
01
1
V(t) - текущее значение объемной доли компонента K (нефти), приведенное к нормальным климатическим условиям;
t - текущее время.
Для измерения скорости W газожидкостного потока в предложенном анализаторе выбран автокорреляционный метод. При этом, в зависимости от режима контролируемого потока, определяемого контроллером режимов 32, для измерения скорости используется либо информация о движении локальной неоднородности поля скоростей контролируемого потока, либо информация о движении локальной неоднородности компонентного состава контролируемой среды.
Выявление упомянутых локальных неоднородностей производится во втором модуле обработки 28 путем обработки электрических сигналов, поступающих на его вход по многоканальной информационной связи с выхода второго фотоэлектронного преобразователя 27, на соответствующий вход которого по каналам второго световода 20 поступают оптические сигналы, вырабатываемые вторыми детекторами 12, а на другой вход, по каналам третьего световода 21, - оптические сигналы, вырабатываемые третьими детекторами 13 (см. Фиг.3). Упомянутые оптические сигналы формируются в результате воздействия излучения рентгеновской трубки 5, пучки которого, сформированные вторым первичным коллиматором 9 и третьим первичным коллиматором 10 и отраженные, соответственно, первым зеркалом 6 и вторым зеркалом 7, пересекают контролируемую среду, находящуюся в корпусе 1, в плоскостях расположения
второй и третьей вставок 3 и 4, соответственно, после чего дополнительно формируются вторым и третьим вторичными коллиматорами 15 и 16, соответственно, и принимаются вторыми и третьими детекторами 12 и 13, соответственно.
Так как при прохождении через находящуюся в корпусе 1 контролируемую среду каждый из упомянутых пучков рентгеновского излучения, взаимодействуя с этой средой, частично поглощается и рассеивается ею, он несет в себе информацию о наличии или отсутствии локальных неоднородностей контролируемого потока.
1
В качестве информативных параметров, характеризующих наличие локальной неоднородности поля скоростей или содержания компонента, например, условно, компонента K, в предложенном анализаторе используются плотности потока фотонов заданной энергии на входах вторых и третьих детекторов 12 и 13, соответственно:
21
31
Y(t), Y(t) - текущие значения плотности потока фотонов на входах вторых и третьих детекторов 12 и 13, соответственно, при наличии в корпусе 1 однокомпонентной контролируемой среды;
2
3
I, I - плотности потока фотонов на входах вторых и третьих детекторов 12 и 13, соответственно, при отсутствии в корпусе 1 контролируемой среды;
e - основание натуральных логарифмов;
1
1
μ - линейный коэффициент ослабления рентгеновского излучения компонентом K контролируемой среды при фиксированном значении энергии излучения;
1
1
x - эффективная толщина компонента K контролируемой среды;
t - текущее время.
1
21
При пересечении движущейся локальной неоднородностью поля скоростей или содержания компонента в момент времени τ сечения А-А корпуса 1 (см. Фиг.1) каждый из вторых детекторов 12 фиксирует наличие существенных изменений значения параметра Y(t) во времени и формирует первую временную последовательность изменений этого значения.
o
При дальнейшем перемещении локальной неоднородности на расстояние L между второй вставкой 3 и третьей вставкой 4, упомянутая неоднородность в момент
2
31
времени τ достигает сечения С-С корпуса 1 (см. Фиг.1) и инициирует в каждом из третьих детекторов 13 существенные изменения текущего значения плотности потока принимаемых им рентгеновских фотонов Y(t), формируя вторую временную последовательность этих изменений, коррелирующую с вышеупомянутой первой временной последовательностью.
21
22
23
31
32
33
1
2
3
При течении вдоль корпуса 1 трехкомпонентной контролируемой среды вторые и третьи детекторы 12 и 13, соответственно, непрерывно вырабатывают первые и вторые временные последовательности информативных сигналов Y(t), Y(t), Y(t) и Y(t), Y(t), Y(t), соответственно, текущие значения амплитуд которых зависят от эффективных толщин x, x и x нефти, воды и газа, соответственно. Первые временные последовательности поступают с выходов вторых детекторов 12 по соответствующим каналам вторых световодов 20 на один из входов второго фотоэлектронного преобразователя 27, вторые временные последовательности поступают с выходов третьих детекторов 13 по соответствующим каналам третьих световодов 21 на другой вход этого преобразователя, где преобразуются из оптической формы в электрическую и передаются по многоканальной информационной связи на вход второго модуля обработки 28, а с его выхода - в вычислитель 29. Алгоритм преобразования сигналов, принятых вычислителем 29, зависит от режима контролируемого потока. Режим потока определяется в контроллере режимов 32 на основании измерительной информации, поступающей на его вход с дополнительного выхода первого модуля обработки 26, где она формируется на основе оптических сигналов, выработанных детекторами первого многоканального детектора 11 и переданных с их выходов по каналам первого световода 19 на соответствующий вход первого преобразователя 25. В первом фотоэлектронном преобразователе 25 упомянутые сигналы преобразуются в электрическую форму и поступают на вход первого модуля обработки 26 по многоканальной информационной связи.
В первом модуле обработки 26 сравниваются поступающие значения плотностей потока фотонов с ранее принятыми и запомненными в памяти упомянутого модуля значениями этих плотностей и, при отсутствии существенных расхождений сравниваемых значений, вырабатывается сигнал о практически однородном режиме течения, а при их наличии - сигнал о существенно неустановившемся
режиме течения газожидкостного потока.
Сформированный в контроллере режимов 32 сигнал о режиме потока передается с его выхода по многоканальной информационной связи на соответствующий многоканальный вход вычислителя 29, в котором выбирается отвечающий принятому сигналу алгоритм преобразования.
Когда контроллером режимов 32 определен режим потока как режим существенно неустановившегося течения, а в вычислителе 29 выбран отвечающий этому режиму алгоритм преобразования измерительной информации, поступающие на вход вычислителя 29 с выхода второго модуля обработки 28 ранее описанные первые и вторые временные последовательности информативных сигналов непрерывно фиксируются в памяти вычислителя 29 в форме первых и вторых временных реализации, причем первые временные реализации соответствуют сечению А-А корпуса 1, а вторые временные реализации - сечению С-С корпуса 1(см. Фиг.1).
min
min
min
С учетом режима существенно неустановившегося движения, определенного контроллером режимов 32, в вычислителе 29 из группы алгоритмов «Вычисление скорости» выбирается алгоритм, отвечающий коду данного режима и, в соответствии с этим алгоритмом, производится обработка вышеуказанных первых и вторых временных реализации за промежуток времени Δτ, минимально необходимый для определения наличия корреляции. Численное значение промежутка времени Δτ вводится в память второго таймера 38 из внешних систем 34, передается с выхода этого таймера на дополнительный вход второго модуля обработки 28 и, с выхода последнего, - на вход вычислителя 29. Промежуток времени Δτ должен быть достаточным для приема вторыми детекторами 12 и третьими детекторами 13 экспериментально установленного числа рентгеновских фотонов
k
4
S≥10, где
k
k
min
S - число фотонов необходимое для корректного определения наличия локальной неоднородности состава контролируемой среды в сечениях А-А и С-С корпуса 1 (см. Фиг.1). Поскольку число S может быть выбрано существенно меньшим числа S фотонов, принимаемых за минимально возможное время экспозиции Δt, то и
min
min
время определения наличия корреляции Δτ составляет величину, существенно меньшую минимально возможного времени экспозиции Δt в выражении (13):
min
min
Δτ≪Δt.
min
1
2
После обработки за промежуток времени Δτ первых и вторых временных реализации, отвечающих моментам времени τ и τ, соответственно, в вычислителе 29 определяется их взаимная корреляционная функция и производится смещение второй реализации относительно первой во времени τ вплоть до получения максимума взаимной корреляционной функции.
к
2
1
o
При получении максимума взаимной корреляционной функции в вычислителе 29 определяется промежуток времени корреляционного смещения Δτ=τ-τ и, поскольку этот промежуток равен времени пробегания Δτ естественной меткой потока - его устойчивой флуктуацией - некоторой длины L корпуса 1, принятой за базовую длину:
к
Δτ=Δτ,
вычисляется скорость W контролируемого потока в соответствии с выражением
W - скорость потока вдоль продольной оси корпуса 1;
o
L - базовая длина, равная осевому расстоянию между геометрическими центрами второй и третьей вставок 3 и 4, соответственно;
Δτ - время пробегания естественной меткой потока базовой длины.
1
2
3
1
2
3
1
2
3
M1
M2
M3
На основании полученного значения скорости W в вычислителе 29 с использованием ранее найденных значений V(t), V(t) и V(t) объемных долей каждого из компонентов контролируемой среды, соответствующих выражению (15), вычисляются мгновенные значения покомпонентных объемных расходов Q(t), Q(t) и Q(t) каждого из трех компонентов газожидкостного потока: нефти, воды и газа, соответственно. Кроме того, в вычислителе 29 с использованием ранее найденных значений M(t), М(t) и М(t) массовых долей каждого из компонентов контролируемой среды, соответствующих выражению (16), вычисляются мгновенные значения покомпонентных массовых расходов Q(t), Q(t) и Q(t) каждого из трех компонентов газожидкостного потока: нефти, воды и газа,
соответственно.
Таким образом, в предложенном анализаторе существенно повышены точность измерения массового содержания компонентов, скорости и расхода контролируемого потока, снижена инструментальная погрешность измерения, вызванная разбросом по величине и дрейфом во времени параметров одноканальных сцинтилляционных детекторов и одноканальных фотоэлектронных преобразователей.