СПОСОБ И ИЗМЕРИТЕЛЬНЫЙ ПРИБОР ДЛЯ ОПРЕДЕЛЕНИЯ СОСТАВА МНОГОФАЗНОЙ ЖИДКОСТИ Российский патент 2002 года по МПК G01N23/08 G01N23/12 

Описание патента на изобретение RU2184367C2

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

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

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

Все эти способы исходят из того, что поглощение пучка фотонов в любом веществе выражается формулой: 1n I0/I = μρd,
где I0 - интенсивность генерируемого излучения,
I - интенсивность пропущенного излучения,
μ - коэффициент массового поглощения излучения, характерный для данного вещества,
ρ - плотность вещества,
d - дальность проникновения излучения сквозь вещество.

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

В Международной патентной заявке РСТ/ЕР 94/01320 рассматривается способ определения состава многофазной жидкости, добываемой из нефтяной скважины, путем пропускания через нее пучка фотонов и измерения уровня поглощения излучения при двух уровнях энергии излучения. В известном способе также принимается во внимание, что три доли: нефти, воды и газа в сумме составляют единицу, т.е. около 100% жидкости, протекающей по трубопроводу, так что можно составить систему из трех математических уравнений, посредством которых данные измерений обрабатываются для получения информации о составе жидкости.

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

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

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

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

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

В преимущественном варианте осуществления способ согласно изобретению содержит:
измерение уровня поглощения излучения жидкостью путем измерения потока фотонов в пропущенном излучении при трех уровнях энергии излучения, логарифмы отношений интенсивности генерируемого излучения Ion к интенсивности пропущенного излучения Intrans на трех уровнях энергии излучения, разделенное на дальность d проникновения излучения через материал обозначаются I1, I2, и I3 и алгоритм вычисления предполагает решение следующего матричного уравнения:

где ρ обозначает плотность компоненты жидкости в кг/м3,
μ обозначает коэффициент массового поглощения в м2/кг,
α обозначает долю объема компоненты в жидкости,
R обозначает возрастание плотности воды при возрастании солености на единицу,
ΔS обозначает изменение солености в кг/м3,
нижний индекс w относится к соленой воде,
нижний индекс о относится к сырой нефти,
нижний индекс g относится к газу,
нижний индекс Н2О относится к пресной воде; и
нижний индекс ini относится к первоначальному состоянию; и
блок обработки данных программируется для решения данного матричного уравнения.

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

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

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

Предположение о том, что соленость добываемой воды практически постоянна в течение промежутка времени порядка нескольких часов избавляет от необходимости решать матричное уравнение точно. Подходящим способом отыскания решения данного уравнения является нахождение статистически оптимального решения в смысле минимизации χ2 на основе методики минимизации функции. Эта методика описана в руководстве "Computational methods in optimization", написанном Е. Полаком и опубликованном издательством Академик пресс в Нью-Йорке в 1971 г.

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

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

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

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

На левой шкале фиг.1 и 2 откладывается фазовая доля а воды w, нефти о, газа g и соли s, а на правой шкале - соленость S в кг/м3. Горизонтальная шкала представляет собой шкалу времени, где время t выражается в секундах (с).

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

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

Кроме того, в течение шести месяцев соленость для некоторых скважин изменилась более чем на 10 кг/м3.

Твердотельная методика детектирования с высоким разрешением, описанная в Международной патентной заявке РСТ/ЕР 94/01320, позволяет разрешать уровни энергии 13,9, 17,8, 22, 26,3 и 59,5 кэВ пучка фотонов, излучаемых источником гамма-лучей Аm-241, и следовательно, эти уровни пригодны для определения солености.

Соли растворяются в добываемой воде, а раз так, они не характеризуются долей объема.

Линейный коэффициент поглощения в соленой воде представляет собой сумму линейных коэффициентов поглощения в соли и в пресной воде:

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

Плотность соленой воды рw возрастает вместе с соленостью приблизительно по линейному закону. Коэффициент возрастания различается для каждой соли и находится в пределах 0,6-0,9. Отсюда плотность соленой воды выражается формулой:

где R обозначает возрастание плотности воды при возрастании солености S на единицу.

Комбинация этих уравнений дает:

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

Измеряя поглощение излучения при трех уровнях энергии 1, 2 и 3 соответственно и принимая во внимание, что скважине для добычи сырой нефти доли нефти, газа и воды в сумме дают практически единицу, получаем следующее матричное уравнение:

или
Это система уравнений, имеющая решение. Измерив поглощение пропущенных гамма-лучей и решив это матричное уравнение, получаем фазовые доли и изменение солености ΔS. Можно показать, что матрица, обратная М, является слабо обусловленной, так что статистическая погрешность измерения может приводить к большим ошибкам в измеряемых фазовых долях и солености. Чтобы оценить выполнимость измерения состава, использующего три или более уровня энергии, весь процесс измерения был смоделирован при помощи компьютерной модели, описанной ниже.

В компьютерной модели использовалась основная ситуация, описанная в таблице (см. в конце описания).

Моделирование состояло из следующих шести этапов.

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

Этап 2
Выбор комбинации долей воды, нефти и газа, первоначальной солености Sini и изменения солености ΔS.
Этап 3
Ожидаемые значения потока фотонов вычисляются с использованием матричного уравнения. Необходимо знать потоки фотонов для пустой трубы при трех уровнях энергии. Они были выбраны реалистически, согласно полевым испытаниям:
Уровень энергии (кэВ) - Поток фотонов (с-1)
17,8 - 21000
26,3 - 2000
59,5 - 80000
Этап 4
Процесс подсчета подвергается статистическим флуктуациям. Измеренные потоки фотонов будут колебаться вокруг вычисленных потоков фотонов согласно распределению Пуассона. В компьютерной модели измеренные потоки фотонов генерируются методом случайной выборки из распределения Пуассона, математическое ожидание и дисперсия которого задается вычисленным потоком фотонов. Обычно моделировалось 3600 наборов из 3 измеренных потоков фотонов, каждый из которых моделирует время измерения 3 секунды.

Этап 5
Эти смоделированные потоки фотонов затем преобразовывались в доли воды, нефти и газа аw, ao, аg и солености S путем решения матричного уравнения согласно изобретению.

Графическое представление результатов этой процедуры изображено на фиг. 1. Входные фазовые доли в этом случае не изменяются со временем, причем доля воды аw=0,15, доля нефти аo=0,15 и доля газа ag=0,70. Входная первоначальная соленость Sini была выбрана равной нулю, а изменение солености ΔS было взято 100 кг/м3. Фиг.1 изображает большие флуктуации долей воды и нефти аw, ao и солености S. Доля газа ag проявляет относительно малую чувствительность к статистическим флуктуациям потоков фотонов. Был сделан вывод, что измеренные фазовые доли (конкретно, ватеркат ≡ аw/(аow)) и соленость становятся крайне чувствительными к статистическим колебаниям измеренных потоков фотонов. Отсюда следует, что в данном случае решение уравнений для каждого измерения не дает практического измерения состава.

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

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

Хи-квадрат , подлежащий минимизации, определяется следующей формулой:

в которой n - это число наборов данных, L - число доступных уровней энергии, причем L равно как минимум трем, и Imji и Icji - логарифмы измеренных и, соответственно, вычисленных потоков фотонов при уровне номер j для набора данных i. Icij вычисляется из фазовых долей и солености согласно матричному уравнению, σij - это стандартная дисперсия при измерении Imij. (В моделях σij, согласно статистике Пуассона, принимается равной квадратному корню из Imij). n наборов фазовых долей и единое постоянное значение для солености, которое минимизирует хи-квадрат, - это наилучшее решение для измерения в статистическом смысле.

Каждый набор данных имеет 2 независимых параметра. Таким образом, для минимизации χ2 всего имеется 2n+1 независимых параметров. Согласно матричному уравнению Icij нелинейно относительно своих независимых параметров, что делает процесс минимизации обременительным упражнением. Процедура минимизации значительно упрощается благодаря разделению линейной и нелинейной частей. Для этой частной задачи минимизации такое разделение естественным образом вытекает из определения χ2.

Отсюда возникает следующая двухэтапная процедура.

Этап 1
Icij является функцией согласно

что следует из матричного уравнения путем исключения ag. Выбирается пробная соленость ΔS. Для этого ΔS варьируют чтобы приблизиться к минимуму χ2i

. Поскольку Ic линейна относительно фазовых долей, минимизацию χ2i
можно осуществить с использованием разложения особого значения, которое представляет собой неитерационный метод и предполагает только действия с матрицами. В процессе минимизации фазовые доли заключены в интервале от нуля до единицы. Поэтому ищут набор фазовых долей в интервале от нуля до единицы, который минимизирует χ2i
для ΔSпроб.
Затем общий хи-квадрат χ2общ
вычисляется путем суммирования n отдельных χ2i
с последующим нормированием.

Этап 2
Общий хи-квадрат χ2общ

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

Использование вышеописанного алгоритма вычисления фазовой доли для того же набора моделированных потоков фотонов, что на фиг.1, приводит к результатам для фазовых долей и солености, изображенным на фиг.2. Фиг.2 демонстрирует, что использование статистического алгоритма минимизации χ2 дает статистически приемлемое определение фазовой доли, которое не чувствительно к изменениям солености.

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

Включение как можно большего числа уровней энергии в процесс измерения имеет очевидное преимущество, поскольку при этом минимально необходимое время для сбора информации существенно уменьшается. Включение уровня 22 кэВ источника гамма-лучей Аm-241 дополнительно к упомянутым ранее уровням 17,8, 26,3 и 59,5 кэВ уменьшает минимальное время сбора информации в 3-4 раза. Включение уровня 13,9 кэВ в качестве пятого уровня дает выигрыш еще в 7-8 раз.

Еще одно преимущество включения как можно большего числа уровней энергии состоит в том, что при этом возрастает статистическая точность измерения состава, поскольку при определении фазовой доли используется в целом больше счетчиков фотонов. При низкой жидкостной нагрузке поток фотонов, получаемый на уровне 17,8 кэВ, преобладает и введение уровней 13,9, 22 и/или 26,3 кэВ не привносит большой статистической точности. При высокой жидкостной нагрузке поток фотонов в пике 17,8 сравним с потоками для уровней 22 и 26,3 кэВ.

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

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

Аналогичная серия аргументов применима к опорному потоку фотонов для нефти. Исходя из потоков фотонов для пустой трубы опорные потоки для нефти и газа можно вычислить на основании плотностей и массовой доли углерода для нефти и газа (что касается газа, необходимо также знать долю других газов, например, N2, CO2, H2S).

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

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

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

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

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

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

Похожие патенты RU2184367C2

название год авторы номер документа
УСТРОЙСТВО ДЛЯ ИЗМЕРЕНИЯ СОСТАВА ФЛЮИДА 1994
  • Александр Мейнхарт Схерс
RU2122724C1
СПОСОБ И УСТРОЙСТВО ДЛЯ ОПРЕДЕЛЕНИЯ СОСТАВА МНОГОФАЗНОГО ПОТОКА СКВАЖИННОЙ ПРОДУКЦИИ 2006
  • Якимов Михаил Николаевич
  • Коркин Роман Владимирович
RU2334972C2
РЕГУЛИРОВАНИЕ ПОТОКА МНОГОФАЗНОЙ ТЕКУЧЕЙ СРЕДЫ, ПОСТУПАЮЩЕЙ ИЗ СКВАЖИНЫ 2005
  • Экен Адриан Николас
RU2386016C2
СПОСОБ И СИСТЕМА СОЗДАНИЯ СТВОЛА СКВАЖИНЫ В ПОЧВЕННОЙ ФОРМАЦИИ 1995
  • Робин Адрианус Хартманн
  • Элвира Хендрика Мюлдер
RU2131975C1
СПОСОБ ОПРЕДЕЛЕНИЯ НАПРАВЛЕНИЯ СКВАЖИНЫ (ВАРИАНТЫ) 1994
  • Николсон Джеймс Уилльям[Us]
RU2109943C1
СПОСОБ ИНТЕРПРЕТАЦИИ СИГНАЛОВ ЯМР ДЛЯ ПОЛУЧЕНИЯ РЕЗУЛЬТАТОВ ИЗМЕРЕНИЙ ПОТОКА МНОГОФАЗНОЙ ТЕКУЧЕЙ СРЕДЫ В СИСТЕМЕ ГАЗ/ЖИДКОСТЬ 2014
  • Аппель Матиас
  • Фримен Джон Джастин
RU2666124C1
СПОСОБ ПОЛУЧЕНИЯ ЭЛЕКТРИЧЕСКОЙ ЭНЕРГИИ ИЗ ПРИРОДНОГО ГАЗА С ИСПОЛЬЗОВАНИЕМ ТОПЛИВНОГО ЭЛЕМЕНТА НА ТВЕРДОМ ОКИСЛЕ 1998
  • Хейнс Майкл Ренни
RU2199172C2
СПОСОБ РЕГУЛИРОВАНИЯ СВОЙСТВ СМАЗКИ ИЛИ РАБОЧЕЙ ЖИДКОСТИ ДЛЯ ИСПОЛЬЗОВАНИЯ В СИСТЕМЕ С ИЗМЕНЯЮЩИМИСЯ УСЛОВИЯМИ ЭКСПЛУАТАЦИИ, СПОСОБ СМАЗКИ СИСТЕМЫ, СПОСОБ ОБЕСПЕЧЕНИЯ РАБОЧЕЙ ЖИДКОСТЬЮ СИСТЕМЫ, УСТРОЙСТВО ДЛЯ ОСУЩЕСТВЛЕНИЯ СПОСОБОВ И ПРИМЕНЕНИЕ КОМПОЗИЦИИ В КАЧЕСТВЕ ОБРАТИМОГО РАЗБАВИТЕЛЯ БАЗОВОЙ ЖИДКОСТИ 2003
  • Вейне Франк Дейвид
RU2319892C2
МУЛЬТИБЛОЧНЫЕ ГИДРОГЕНИЗОВАННЫЕ ПОЛИМЕРЫ ДЛЯ КЛЕЕВ 1994
  • Гленн Рой Хаймс
RU2160755C2
СПОСОБ ОПРЕДЕЛЕНИЯ ХАРАКТЕРИСТИКИ ПОРИСТОЙ ПОДЗЕМНОЙ ФОРМАЦИИ 1995
  • Ридван Аккурт
  • Пьер Назарет Тутуньян
  • Харолд Дж. Вайнегар
RU2134894C1

Иллюстрации к изобретению RU 2 184 367 C2

Реферат патента 2002 года СПОСОБ И ИЗМЕРИТЕЛЬНЫЙ ПРИБОР ДЛЯ ОПРЕДЕЛЕНИЯ СОСТАВА МНОГОФАЗНОЙ ЖИДКОСТИ

Изобретения относятся к области исследования объектов радиационными методами и могут быть использованы для определения состава многофазной жидкости. Состав многофазной жидкости с изменяющейся соленостью определяется путем пропускания пучка фотонов через жидкость и измерения уровня поглощения жидкостью излучения при как минимум трех уровнях энергии излучения и подачи полученных таким образом данных о поглощении излучения на обработчик данных, который выдает данные, касающиеся состава жидкости, включая содержание соли, если есть, на основании алгоритма вычисления фазовой доли, в котором принимается во внимание поглощение излучения солью. Техническим результатом изобретений является увеличение точности измерений, возможность учета соленой составляющей жидкости, упрощение калибровки измерителей. 2 с. и 6 з.п. ф-лы, 1 табл., 2 ил.

Формула изобретения RU 2 184 367 C2

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

при этом w, H2O и s относятся к соленой воде, пресной воде и соли соответственно;
ini - к первоначальной солености;
ρ - плотность, кг/м3;
μ - коэффициент массового поглощения, м2/кг;
R обозначает возрастание плотности воды при увеличении солености на единицу;
ΔS обозначает изменение солености, кг/м3,
выбирают линейный коэффициент поглощения воды, имеющей начальную соленость ρwμw|ini, и определяют изменение солености из решения матричного уравнения.
2. Способ по п. 1, в котором уровень поглощения излучения жидкостью измеряют путем измерения потока фотонов пропущенного излучения при трех уровнях энергии излучения: логарифмы отношений интенсивности генерируемого излучения Ion к интенсивности Intrans пропущенного излучения при трех уровнях энергии излучения (n= 1, 2, 3), поделенное на дальность d проникновения излучения сквозь вещество, обозначаются I1, I2, I3, и уравнение расчета содержит следующую схему матричного уравнения:

где ρ - обозначает плотность компоненты жидкости, кг/м3;
μ - обозначает коэффициент массового поглощения, м2/кг;
α обозначает долю объема компоненты в жидкости;
R обозначает возрастание плотности воды при возрастании солености на единицу;
ΔS - обозначает изменение солености, кг/м3;
нижние индексы: w относится к соленой воде; s относится к соли; o относится к сырой нефти; g относится к газу; Н2О относится к пресной воде; ini относится к воде, имеющей первоначальную соленость,
и блок обработки данных программируют на решение этого матричного уравнения.
3. Способ по п. 2, в котором уровень поглощения излучения измеряют неоднократно и блок обработки данных программируют так, что он учитывает, что содержание соли в водной фазе жидкости остается практически постоянным на протяжении промежутка времени, ограниченного, по крайней мере, несколькими часами и так, что он решает указанное матричное уравнение путем отыскания статистического оптимума. 4. Способ по п. 3, в котором блок обработки данных программируется так, что он ищет указанный статистический оптимум на основании статистического алгоритма минимизации хи-квадрат. 5. Способ по любому из предшествующих пунктов, в котором многофазная жидкость истекает из скважины для добычи углеводородной жидкости, которая выдает сырую нефть совместно с природным газом и/или соленой водой, и фазовые доли указанных составляющих флуктуируют и колеблются в пределах от 0 до 100%, при условии, что скважина выдает указанные составляющие в вялотекущем режиме. 6. Измерительный прибор для определения состава многофазной жидкости, включающей соленую воду, содержащий источник излучения для пропускания пучка фотонов сквозь жидкость, детектор излучения для измерения уровня поглощения излучения жидкостью при, как минимум, трех уровнях энергии излучения и блок обработки данных, который программируется с возможностью вычисления, в соответствии с алгоритмом вычисления, фазовой доли на основании данных о поглощении излучения и возможностью выдачи данных, касающихся состава жидкости, на основании указанных вычислений, при этом в блок обработки данных вводится алгоритм вычисления, отличающийся тем, что введенный алгоритм вычисления устраивается с возможностью принятия во внимание воздействие содержания соли в соленой воде на поглощение жидкостью пучка фотонов при вычислении фазовой доли, при этом при введенном алгоритме вычисления вычисляется матричное уравнение, относящееся к измеренным данным поглощения излучения, к фазовым долям многофазной жидкости и к изменению в солености соленой воды в отношении воды, имеющей первоначальную соленость, и линейный коэффициент поглощения ρwμw соленой воды определяется на основании следующего равенства:

при этом w, H2O и s относятся к соленой воде, пресной воде и соли соответственно;
ini - к первоначальной солености;
ρ - плотность, кг/м3;
μ - коэффициент массового поглощения, м2/кг;
R обозначает возрастание плотности воды при увеличении солености на единицу;
ΔS - обозначает изменение солености, кг/м3,
при этом алгоритм вычисления для определения изменения в солености вычисляется путем решения матричного уравнения, используется выбранный линейный коэффициент поглощения воды, имеющей первоначальную соленость, ρwμw|ini.
7. Измерительный прибор по п. 6, в котором источник излучения пучка фотонов содержит рентгеновскую трубку.
8. Измерительный прибор по п. 6, в котором источник излучения пучка фотонов содержит источник гамма-лучей.

Документы, цитированные в отчете о поиске Патент 2002 года RU2184367C2

ФОРМИРОВАТЕЛЬ СТРОБИМПУЛЬСОВ ДЛЯ ШИРОКОПОЛОСНОГО ОСЦИЛЛОГРАФА 0
SU236623A1
СУШИЛЬНАЯ БЫТОВАЯ СВЧ-ПЕЧЬ 1994
  • Кирсанов Ю.А.
  • Тишков В.С.
RU2088050C1
US 4618975 A, 21.10.1986
RU 94008179 A1, 10.02.1996.

RU 2 184 367 C2

Авторы

Схерс Александер Мейнхарт

Слейкерман Вальтер Фредерикус Йоханнес

Даты

2002-06-27Публикация

1997-04-25Подача