RU2105145C1 - Способ измерения расхода фаз газожидкостного потока - Google Patents
Способ измерения расхода фаз газожидкостного потока Download PDFInfo
- Publication number
- RU2105145C1 RU2105145C1 RU96114284A RU96114284A RU2105145C1 RU 2105145 C1 RU2105145 C1 RU 2105145C1 RU 96114284 A RU96114284 A RU 96114284A RU 96114284 A RU96114284 A RU 96114284A RU 2105145 C1 RU2105145 C1 RU 2105145C1
- Authority
- RU
- Russia
- Prior art keywords
- gas
- flow rate
- liquid
- oil
- signal
- Prior art date
Links
Images
Landscapes
- Measuring Volume Flow (AREA)
Abstract
Назначение: изобретение относится к нефтегазодобывающей промышленности для раздельного измерения расхода двухфазных потоков эксплуатационных нефтяных, газовых и газоконденсатных скважин. Сущность изобретения: способ основан на вычислении среднеквадратических значений сигнала пульсаций давления в двух информативных диапазонах частот. Дополнительно определяют функциональные коэффициенты расходного газосодержания и по произведению соответствующего среднеквадратического значения пульсации давления на соответствующий функциональный коэффициент определяют расход фаз газожидкостной смеси. Применение способа позволяет повысить точность раздельного измерения дебита нефтегазовых скважин по жидкости и газу в условиях изменения расходного газосодержания потока в широких пределах. 4 ил.
Description
Предлагаемое изобретение относится к области нефтегазодобывающей промышленности и может быть использовано при измерении дебита двухфазных потоков эксплуатационных газовых, газоконденсатных и нефтяных скважин.
Известен способ определения дебита скважин, по которому осуществляют измерение пульсаций давления газожидкостного потока и определение среднеквадратического значения сигнала пульсаций давления, по интегральному значению которого судят о величине массового расхода нефти в потоке [1].
Недостаток этого способа заключается в том, что он обеспечивает определение только расхода нефти и не дает информации о расходе газа в потоке.
Наиболее близким к предлагаемому является способ раздельного определения дебита скважины по жидкости и газу путем измерения пульсаций давления потока продукции с последующей их фильтрацией в диапазонах нижних и верхних частот [2] . При этом о расходе смеси жидкости и газа и о расходе газа судят по интегральной величине среднеквадратического значения сигнала пульсаций давления в указанных диапазонах частот.
Недостатком известного способа является то, что он не учитывает многообразия режимов работы нефтяных, газовых и газоконденсатных скважин, связанных с изменением газового фактора в широких пределах, вследствие чего возникают недопустимые погрешности измерения расхода фаз.
Задачей изобретения является повышение точности раздельного измерения расхода жидкости и газа в потоке смеси при изменении расходного газосодержания в широких пределах.
Поставленная задача достигается тем, что в способе измерения расхода фаз газожидкостной смеси путем измерения пульсаций давления потока, фильтрации сигналов пульсаций давления в двух информативных диапазонах частот и определения среднеквадратических значений сигнала пульсаций давления в первом и втором частотных диапазонах, согласно изобретению, дополнительно вычисляют отношение среднеквадратических значений сигнала пульсаций давления второго и первого частотных диапазонов и отношение среднеквадратических значений сигнала пульсаций давления первого и второго частотных диапазонов, по которым определяют функциональные коэффициенты расходного газосодержания, соответственно для жидкости и газа, и расход жидкости находят как произведение среднеквадратического значения сигнала пульсаций давления в первом частотном диапазоне и функционального коэффициента расходного газосодержания для жидкости, а расход газа - как произведение среднеквадратического значения сигнала пульсаций давления во втором частотном диапазоне и функционального коэффициента расходного газосодержания для газа.
Сущность предлагаемого способа заключается в следующем.
Обобщенная информационная модель расхода газожидкостной смеси при фиксированном сужающем устройстве представлена уравнением:
где Q - объемный расход газожидкостной смеси;
U - среднеквадратическое значение сигнала пульсаций давления;
ρ - плотность газожидкостной среды;
Ψ(β) - безразмерная функция, учитывающая влияние коэффициента расходного газосодержания β ;
K - постоянный коэффициент.
где Q - объемный расход газожидкостной смеси;
U - среднеквадратическое значение сигнала пульсаций давления;
ρ - плотность газожидкостной среды;
Ψ(β) - безразмерная функция, учитывающая влияние коэффициента расходного газосодержания β ;
K - постоянный коэффициент.
Так как плотность среды определяется давлением потока продукции и расходным газосодержанием, информационная модель сводится к следующему виду:
Q = K*•U1/3•Ψ1(P)•Ψ2(β), (2)
где K* - постоянный коэффициент;
Ψ1(P) - функция, учитывающая влияние давления потока продукции на устье скважины;
Ψ2(β) - функция, учитывающая влияние коэффициента расходного газосодержания.
Q = K*•U1/3•Ψ1(P)•Ψ2(β), (2)
где K* - постоянный коэффициент;
Ψ1(P) - функция, учитывающая влияние давления потока продукции на устье скважины;
Ψ2(β) - функция, учитывающая влияние коэффициента расходного газосодержания.
С учетом известных закономерностей поведения спектра выходного сигнала пьезокерамического датчика, обусловленного особенностями движения газожидкостной смеси через сужающее устройство, информационные модели расхода фаз можно описать такими уравнениями:
Qж, Qг - расходы, соответственно, жидкости и газа;
U1, U2 - среднеквадратические значения сигнала пульсаций давления, вычисленные в первом и втором частотных диапазонах наибольшего влияния жидкости и газа, соответственно;
α, γ - показатели степени;
f1(P), f2(P) - функции, учитывающие влияние давления потока продукции на устье скважины, соответственно, в первом и втором частотных диапазонах;
f3(β), f4(β) - функции, учитывающие влияние коэффициента расходного газосодержания, соответственно, в первом и втором частотных диапазонах.
Qж, Qг - расходы, соответственно, жидкости и газа;
U1, U2 - среднеквадратические значения сигнала пульсаций давления, вычисленные в первом и втором частотных диапазонах наибольшего влияния жидкости и газа, соответственно;
α, γ - показатели степени;
f1(P), f2(P) - функции, учитывающие влияние давления потока продукции на устье скважины, соответственно, в первом и втором частотных диапазонах;
f3(β), f4(β) - функции, учитывающие влияние коэффициента расходного газосодержания, соответственно, в первом и втором частотных диапазонах.
Влияние расходного газосодержания в уравнениях расхода (3, 4) можно оценить с помощью функций f3(U21), f4(U12), аргументы которых составлены из соотношений:
Тогда информационные модели расхода фаз, учитывающие изменение расходного газосодержания потока в широких пределах, имеют вид:
где f1(P)=1+a1P+b1P2,
f2(P)=1+a2P+b2P2,
f3(U21), f4(U12) - функциональные коэффициенты расходного газосодержания, соответственно, для жидкости и газа;
a1, b1, c1, d1
a2, b2, c2, d2 - постоянные параметры полиномиальных моделей;
K1, K2 - постоянные коэффициенты.
Тогда информационные модели расхода фаз, учитывающие изменение расходного газосодержания потока в широких пределах, имеют вид:
где f1(P)=1+a1P+b1P2,
f2(P)=1+a2P+b2P2,
f3(U21), f4(U12) - функциональные коэффициенты расходного газосодержания, соответственно, для жидкости и газа;
a1, b1, c1, d1
a2, b2, c2, d2 - постоянные параметры полиномиальных моделей;
K1, K2 - постоянные коэффициенты.
Блок-схема измерительной установки, реализующей заявляемый способ, приведена на фиг. 1. Измерительная установка состоит из измерительного модуля 1, в состав которого входят пьезокерамический датчик 2 пульсаций давления потока и согласующий усилитель-корректор 3, и вторичного прибора 4. В состав вторичного прибора входят управляемый масштабирующий усилитель 5, фильтры нижних и верхних частот 6 и 7, коммутатор 8, аналого-цифровой преобразователь 9, микропроцессорный контроллер 10, цифровые индикаторы расхода жидкости 11 и газа 12, клавиатура 13. Для коррекции вычислений при изменении давления потока продукции на устье скважины установлен датчик давления 14.
Измерительный модуль 1 устанавливается на трубопроводе 15 на заданном расстоянии от специального сужающего устройства 16, предназначенного для более интенсивной турбулизации потока и формирования его структуры.
Пьезокерамический элемент 2 осуществляет преобразование пульсаций давления потока в пропорциональный электрический сигнал (далее по тексту -сигнал), который поступает на согласующий усилитель - корректор 3, служащий для предварительного усиления сигнала, согласования высокоомного сопротивления пьезокерамического датчика с входным сопротивлением вторичного измерительного прибора 4 и коррекции амплитудно - частотной характеристики. Сигнал с выхода согласующего усилителя - корректора 3 поступает на вход масштабирующего усилителя 5. Оптимальный коэффициент усиления этого усилителя задается автоматически от микропроцессорного контроллера 10 либо оператором с помощью клавиатуры 13. С выхода масштабирующего усилителя 5 сигнал поступает на фильтры нижних и верхних частот, соответственно, 6 и 7. Фильтр нижних частот 6 выделяет первую информационную полосу частот (первый частотный диапазон), интенсивность сигнала в которой, в основном, связана с расходом жидкой фазы. С выхода фильтра 6 сигнал поступает на первый вход коммутатора 8, на второй вход которого поступает сигнал с фильтра верхних частот 7, выделяющего вторую информационную полосу частот (второй частотный диапазон), интенсивность сигнала в которой, в основном, связана с расходом газовой фазы. На третий вход коммутатора 8 поступает сигнал от датчика давления 14.
С выхода коммутатора сигнал поступает на аналого-цифровой преобразователь (АЦП) 9, где преобразуется в двоичный код. Выходной сигнал АЦП поступает на вход микропроцессорного контроллера 10, который вычисляет значения расходов жидкости и газа в соответствии с формулами (5, 6).
Результаты измерений индицируются на цифровых индикаторах 11 (жидкая фаза) и 12 (газовая фаза).
Алгоритм работы микропроцессорного контроллера 10 приведен на фиг.2. Он содержит следующие основные операторы:
17 - пуск;
18 - подпрограмма самотестирования;
19 - подпрограмма инициализации ресурсов системы (задание по умолчанию значений коэффициента усиления K масштабирующего усилия 5 и количества циклов измерений N);
20 - опрос клавиш клавиатуры изменения K и N;
21 - ввод новых значений K или N;
22 - опрос клавиши "измерение";
23 - проверка окончания времени одного цикла измерений;
24 - переключение коммутатора 8 на первый вход;
25 - чтение микропроцессорным контроллером 10 из АЦП 9 кода U1i, соответствующего значению сигнала на выходе фильтра 6;
26 - переключение коммутатора 8 на второй вход;
27 - чтение микропроцессорным контроллером 10 из АЦП 9 кода U2i, соответствующего значению сигнала на выходе фильтра 7;
28 - переключение коммутатора 8 на третий вход;
29 - чтение микропроцессорным контроллером 10 из АЦП 9 кода P, соответствующего значению давления от датчика 14;
30 - расчет оптимального значения коэффициента усиления K;
31 - корректировка коэффициента усиления К масштабирующего усилителя 5;
32 - расчет
33 - проверка количества циклов измерений (i=N ?);
34 - расчет значений U и U , соответственно, для первого и второго частотных диапазонов по формулам:
35 - расчет отношений U21 и U12;
36 - расчет f1(p), f2(p), f3(U21), f4(U12);
37 - вычисление Qж и Qг в соответствии с формулами 5 и 6;
38 - вывод значений Qж и Qг;
39 - конец.
17 - пуск;
18 - подпрограмма самотестирования;
19 - подпрограмма инициализации ресурсов системы (задание по умолчанию значений коэффициента усиления K масштабирующего усилия 5 и количества циклов измерений N);
20 - опрос клавиш клавиатуры изменения K и N;
21 - ввод новых значений K или N;
22 - опрос клавиши "измерение";
23 - проверка окончания времени одного цикла измерений;
24 - переключение коммутатора 8 на первый вход;
25 - чтение микропроцессорным контроллером 10 из АЦП 9 кода U1i, соответствующего значению сигнала на выходе фильтра 6;
26 - переключение коммутатора 8 на второй вход;
27 - чтение микропроцессорным контроллером 10 из АЦП 9 кода U2i, соответствующего значению сигнала на выходе фильтра 7;
28 - переключение коммутатора 8 на третий вход;
29 - чтение микропроцессорным контроллером 10 из АЦП 9 кода P, соответствующего значению давления от датчика 14;
30 - расчет оптимального значения коэффициента усиления K;
31 - корректировка коэффициента усиления К масштабирующего усилителя 5;
32 - расчет
33 - проверка количества циклов измерений (i=N ?);
34 - расчет значений U
35 - расчет отношений U21 и U12;
36 - расчет f1(p), f2(p), f3(U21), f4(U12);
37 - вычисление Qж и Qг в соответствии с формулами 5 и 6;
38 - вывод значений Qж и Qг;
39 - конец.
Предлагаемый способ позволяет производить раздельное измерение расхода жидкости и газа в потоке продукции скважин с требуемой точностью в условиях широкого изменения расходного газосодержания потока.
Повышение точности раздельного измерения дебита нефтегазовой скважины при применении заявляемого способа измерения показано на примере, представленном в виде графиков на фиг. 3, 4. Эти графики отражают результаты измерений жидкости и газа в газожидкостных потоках, моделируемых на водонефтегазовом метрологическом полигоне при гидродинамических условиях, соответствующих режимам работы реальных нефтегазовых скважин. По оси ординат отложены фактические расходы газа (a) и жидкости (б), измеренные после разделения газожидкостной смеси на жидкость и газ в сепараторе.
По оси абсцисс отложены расходы газа (а) и жидкости (б), измеренные способом без учета изменения расходного газосодержания потока (фиг.3) и с учетом изменения расходного газосодержания потока предлагаемым способом (фиг.4).
Из графиков фиг.3 видно, что приведенная погрешность измерения расходов достигает 20% по газу и 25% по жидкости. При измерении расхода предлагаемым способом (фиг. 4) точность измерений существенно возрастает. Приведенная погрешность измерения расхода газа (фиг. 4, а) не превышает 5%, а при измерении расхода жидкости (фиг. 4, б) погрешность измерения не превышает 8%.
1. Патент США N 3834227, кл. 73-155, 1974.
2. Патент РФ N 1060791, МПК E 21 B 47/00, 1991.
Claims (1)
- Способ измерения расхода фаз газожидкостной смеси путем измерения пульсаций давления потока, фильтрации сигналов пульсаций давления в двух информативных диапазонах частот и определения среднеквадратических значений сигнала пульсаций давления в первом и втором частотных диапазонах, отличающийся тем, что дополнительно вычисляют отношение среднеквадратических значений сигнала пульсаций давления второго и первого частотных диапазонов и отношение среднеквадратических значений сигнала пульсаций давления первого и второго частотных диапазонов, по которым определяют функциональные коэффициенты расходного газосодержания соответственно для жидкости и газа, и расход жидкости находят как произведение среднеквадратического значения сигнала пульсаций давления в первом частотном диапазоне и функционального коэффициента расходного газосодержания для жидкости, а расход газа как произведение среднеквадратического значения сигнала пульсаций давления во втором частотном диапазоне и функционального коэффициента расходного газосодержания для газа.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU96114284A RU2105145C1 (ru) | 1996-07-17 | 1996-07-17 | Способ измерения расхода фаз газожидкостного потока |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU96114284A RU2105145C1 (ru) | 1996-07-17 | 1996-07-17 | Способ измерения расхода фаз газожидкостного потока |
Publications (2)
Publication Number | Publication Date |
---|---|
RU96114284A RU96114284A (ru) | 1998-01-27 |
RU2105145C1 true RU2105145C1 (ru) | 1998-02-20 |
Family
ID=20183291
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU96114284A RU2105145C1 (ru) | 1996-07-17 | 1996-07-17 | Способ измерения расхода фаз газожидкостного потока |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2105145C1 (ru) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2490449C2 (ru) * | 2011-09-08 | 2013-08-20 | Общество с ограниченной ответственностью "Газпром добыча Надым" (ООО "Газпром добыча Надым") | Способ гидрогазодинамических исследований скважин |
RU2527525C1 (ru) * | 2013-04-02 | 2014-09-10 | Открытое акционерное общество "Северо-Кавказский научно-исследовательский проектный институт природных газов" (ОАО "СевКавНИПИгаз") | Способ газодинамического исследования скважины |
CN110863812A (zh) * | 2018-08-14 | 2020-03-06 | 中国科学院沈阳自动化研究所 | 一种基于油套分离比的油井产气量计量方法 |
-
1996
- 1996-07-17 RU RU96114284A patent/RU2105145C1/ru not_active IP Right Cessation
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2490449C2 (ru) * | 2011-09-08 | 2013-08-20 | Общество с ограниченной ответственностью "Газпром добыча Надым" (ООО "Газпром добыча Надым") | Способ гидрогазодинамических исследований скважин |
RU2527525C1 (ru) * | 2013-04-02 | 2014-09-10 | Открытое акционерное общество "Северо-Кавказский научно-исследовательский проектный институт природных газов" (ОАО "СевКавНИПИгаз") | Способ газодинамического исследования скважины |
CN110863812A (zh) * | 2018-08-14 | 2020-03-06 | 中国科学院沈阳自动化研究所 | 一种基于油套分离比的油井产气量计量方法 |
CN110863812B (zh) * | 2018-08-14 | 2021-07-09 | 中国科学院沈阳自动化研究所 | 一种基于油套分离比的油井产气量计量方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Tu et al. | Fully developed periodic turbulent pipe flow. Part 1. Main experimental results and comparison with predictions | |
US5051922A (en) | Method and apparatus for the measurement of gas/liquid flow | |
CA2539640C (en) | Detection and measurement of two-phase flow | |
US5597961A (en) | Two and three phase flow metering with a water cut monitor and an orifice plate | |
JP4904339B2 (ja) | 気体流材料における液体流の割合を決定するための計量器電子装置及び方法 | |
JP5250035B2 (ja) | 流量計システム及び3相の流れの流れ特性を測定する方法 | |
EP1666850A1 (en) | Three-phase flow regulating means for oil, gas and water, three-phase flow measuring apparatus for oil, gas and water and measuring method thereof | |
JP2006500557A5 (ru) | ||
CA2262444A1 (en) | Method and apparatus for measuring pressure in a coriolis mass flowmeter | |
CN102216739A (zh) | 用于测量振动计中流体参数的方法和设备 | |
CA2536341C (en) | Flow meter filter system and method | |
EP1520154A2 (en) | Monitoring of two-phase fluid flow using a vortex flowmeter | |
WO2000057138A1 (en) | Method and apparatus for measuring the mass flow rate of a fluid | |
RU2105145C1 (ru) | Способ измерения расхода фаз газожидкостного потока | |
US11940309B2 (en) | Method for determining a fluid flow parameter within a vibrating tube | |
AU672762B2 (en) | Acoustic displacement flow meter | |
RU2103502C1 (ru) | Устройство для контроля дебита газовых, газоконденсатных и нефтяных скважин | |
RU2103503C1 (ru) | Устройство для контроля дебитов компонентов продукции скважин | |
RU2101698C1 (ru) | Устройство для измерения концентрации взвешенных веществ в жидкости | |
RU1789859C (ru) | Способ определени массового расхода среды в трубопроводе | |
RU1795287C (ru) | Способ определени массового расхода газа | |
RU2439502C2 (ru) | Система измерителя потока и способ для измерения параметров трехфазного потока | |
JP3498785B2 (ja) | 渦流量計 | |
RU2279641C2 (ru) | Способ измерения весового расхода газожидкостной смеси и устройство для его осуществления | |
RU2262082C2 (ru) | Установка для измерения расхода газожидкостной текучей среды |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20110718 |