Изобретение относится к способу и измерительному прибору для определения состава многофазной жидкости путем пропускания через нее пучка фотонов.
Конкретно изобретение относится к способу и измерительному прибору для определения состава многофазной жидкой смеси, выдаваемой одной или несколькими скважинами для добычи сырой нефти, где сырой нефти обычно сопутствует некоторое количество природного газа и/или воды. Такое измерение предоставляет информацию о составе жидкости и схему потока в трубопроводе от каждой скважины, что дает понимание, необходимое для оптимального управления системой добычи и нефтеносными пластами, откуда добывают сырую нефть.
Существует несколько способов определения состава подобной многофазной жидкости путем пропускания через нее пучка фотонов, где пучок фотонов обычно порождается источником гамма- и/или рентгеновских лучей.
Все эти способы исходят из того, что поглощение пучка фотонов в любом веществе выражается формулой: 1n I0/I = μρd,
где I0 - интенсивность генерируемого излучения,
I - интенсивность пропущенного излучения,
μ - коэффициент массового поглощения излучения, характерный для данного вещества,
ρ - плотность вещества,
d - дальность проникновения излучения сквозь вещество.
В этих способах также учитывается, что при конкретном уровне энергии излучения каждое вещество имеет соответствующий коэффициент μ массового поглощения излучения, который характерен для данного вещества и, кроме того, что при изменении уровня энергии излучения коэффициент μ массового поглощения излучения также изменяется по закону, характерному для данного вещества.
В Международной патентной заявке РСТ/ЕР 94/01320 рассматривается способ определения состава многофазной жидкости, добываемой из нефтяной скважины, путем пропускания через нее пучка фотонов и измерения уровня поглощения излучения при двух уровнях энергии излучения. В известном способе также принимается во внимание, что три доли: нефти, воды и газа в сумме составляют единицу, т.е. около 100% жидкости, протекающей по трубопроводу, так что можно составить систему из трех математических уравнений, посредством которых данные измерений обрабатываются для получения информации о составе жидкости.
В описании европейского патента 236623 идет речь о том, что имеет определенный смысл измерять уровень поглощения излучения пучка фотонов при более чем двух уровнях энергии для получения более расширенного набора данных о поглощении излучения, из которого можно вывести состав жидкости. Этот известный способ также дает возможность учитывать поглощение излучения, обусловленное частицами песка, погруженными в жидкую смесь, и содержанием серы в сырой нефти за счет использования дополнительных систем излучения, которые излучают пучки фотонов с другими уровнями энергии.
Хотя эти известные способы измерения предоставляют полезную информацию, было обнаружено, что известные измерительные приборы требуют обширной калибровки, которую требуется регулярно обновлять, но и это не избавляет процесс измерения от некоторой неточности.
В основу настоящего изобретения положена задача увеличения точности определения состава многофазной жидкости и предоставление более простых и менее частых процедур калибровки измерителей многофазного состава.
Способ согласно изобретению содержит определение состава многофазной жидкости путем пропускания через нее пучка фотонов и измерения уровня поглощения излучения жидкостью при как минимум трех уровнях энергии излучения и передачу данных измерения поглощения излучения на блок обработки данных, который программируется так, что он осуществляет вычисления в соответствии с алгоритмом вычисления фазовой доли на основании упомянутых данных о поглощении излучения и на основании упомянутых вычислений выдает данные, касающиеся состава жидкости, причем жидкость содержит соленую воду и алгоритм расчета фазовой доли включает в себя этапы определения солености воды.
Способ согласно изобретению базируются на понимании того факта, что содержание соли, если таковая имеется, в воде, выдаваемой, например, скважиной для добычи сырой нефти, может оказывать существенное воздействие на поглощение жидкостью пучка фотонов.
В преимущественном варианте осуществления способ согласно изобретению содержит:
измерение уровня поглощения излучения жидкостью путем измерения потока фотонов в пропущенном излучении при трех уровнях энергии излучения, логарифмы отношений интенсивности генерируемого излучения Ion к интенсивности пропущенного излучения In trans на трех уровнях энергии излучения, разделенное на дальность 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/(аo+аw)) и соленость становятся крайне чувствительными к статистическим колебаниям измеренных потоков фотонов. Отсюда следует, что в данном случае решение уравнений для каждого измерения не дает практического измерения состава.
Чтобы усовершенствовать расчет фазовых долей, был разработан статистический алгоритм минимизации χ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 варьируют чтобы приблизиться к минимуму χ
Затем общий хи-квадрат χ
Этап 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 согласно изобретению может вводиться в обработчик данных как программа вычисления и обработчик данных может производить отображение информации, демонстрируя доли нефти, газа, воды и соленость добываемой многофазной жидкости или набор этих характеристик или их производных.
Наконец, следует принять во внимание, что способ и измерительный прибор, соответствующие изобретению, полезны для определения состава любой соленой многофазной жидкости, в котором требуется учитывать поглощение излучения соленой составляющей, с учетом того, что для жидкостей, отличных от тех, что вытекают из углеводородных скважин, алгоритм матричного вычисления и алгоритм минимизации хи-квадрат подлежат соответствующей минимизации.
Изобретения относятся к области исследования объектов радиационными методами и могут быть использованы для определения состава многофазной жидкости. Состав многофазной жидкости с изменяющейся соленостью определяется путем пропускания пучка фотонов через жидкость и измерения уровня поглощения жидкостью излучения при как минимум трех уровнях энергии излучения и подачи полученных таким образом данных о поглощении излучения на обработчик данных, который выдает данные, касающиеся состава жидкости, включая содержание соли, если есть, на основании алгоритма вычисления фазовой доли, в котором принимается во внимание поглощение излучения солью. Техническим результатом изобретений является увеличение точности измерений, возможность учета соленой составляющей жидкости, упрощение калибровки измерителей. 2 с. и 6 з.п. ф-лы, 1 табл., 2 ил.
при этом w, H2O и s относятся к соленой воде, пресной воде и соли соответственно;
ini - к первоначальной солености;
ρ - плотность, кг/м3;
μ - коэффициент массового поглощения, м2/кг;
R обозначает возрастание плотности воды при увеличении солености на единицу;
ΔS обозначает изменение солености, кг/м3,
выбирают линейный коэффициент поглощения воды, имеющей начальную соленость ρwμw|ini, и определяют изменение солености из решения матричного уравнения.
где ρ - обозначает плотность компоненты жидкости, кг/м3;
μ - обозначает коэффициент массового поглощения, м2/кг;
α обозначает долю объема компоненты в жидкости;
R обозначает возрастание плотности воды при возрастании солености на единицу;
ΔS - обозначает изменение солености, кг/м3;
нижние индексы: w относится к соленой воде; s относится к соли; o относится к сырой нефти; g относится к газу; Н2О относится к пресной воде; ini относится к воде, имеющей первоначальную соленость,
и блок обработки данных программируют на решение этого матричного уравнения.
при этом w, H2O и s относятся к соленой воде, пресной воде и соли соответственно;
ini - к первоначальной солености;
ρ - плотность, кг/м3;
μ - коэффициент массового поглощения, м2/кг;
R обозначает возрастание плотности воды при увеличении солености на единицу;
ΔS - обозначает изменение солености, кг/м3,
при этом алгоритм вычисления для определения изменения в солености вычисляется путем решения матричного уравнения, используется выбранный линейный коэффициент поглощения воды, имеющей первоначальную соленость, ρwμw|ini.
7. Измерительный прибор по п. 6, в котором источник излучения пучка фотонов содержит рентгеновскую трубку.
ФОРМИРОВАТЕЛЬ СТРОБИМПУЛЬСОВ ДЛЯ ШИРОКОПОЛОСНОГО ОСЦИЛЛОГРАФА | 0 |
|
SU236623A1 |
СУШИЛЬНАЯ БЫТОВАЯ СВЧ-ПЕЧЬ | 1994 |
|
RU2088050C1 |
US 4618975 A, 21.10.1986 | |||
RU 94008179 A1, 10.02.1996. |
Авторы
Даты
2002-06-27—Публикация
1997-04-25—Подача