Результат интеллектуальной деятельности: СПОСОБ ОПРЕДЕЛЕНИЯ УРОВНЯ РАЗДЕЛА ФАЗ В КАНАЛАХ

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

Известно большое количество уровнемеров. Работа уровнемеров основана на различных физических принципах. Известен гидростатический уровнемер (В.П. Преображенский. Теплотехнические измерения и приборы, третье издание, переработанное. - Москва «Энергия», 1978 г., с.530). Принцип работы гидростатического уровнемера заключается в измерении давления столба жидкости и определении на основе измеренного давления столба жидкости уровня раздела фаз.

Основной недостаток гидростатического уровнемера при использовании его для измерения уровня в каналах с двухфазным слоем состоит в том, что плотность теплоносителя, на основе которого определяется высота столба жидкой фазы отличается от плотности воды на линии насыщения, которая используется для расчета уровнемера (В.А.Демченко. О точности измерения уровня воды в парогенераторах АЭС. Теплоэнергетика, №2, 1999, с.56). Ошибка определения уровня может быть значительной (25-35%).

Наиболее близким по технической сущности и достигаемому результату к предлагаемому изобретению является способ определения уровня газовой фазы в канале заключающийся в том, что устанавливают датчик, выполненный, например, в виде электропроводной проволоки, в канале по направлению силы тяжести, нагревают датчик путем пропускания тока через датчик, измеряют электрическое сопротивление датчика R. (Детектирование нижнего уровня масла с помощью термисторов (Low oil level sensing with thermistors. Loisch J. SAE Techn. Pap.Ser., 1989, №891758, 67-70). Контрольно-измерительная техника, №4, 1991 г. Государственный комитет по науке и технике. Академия наук СССР. Всесоюзный институт научной и технической информации (ВИНИТИ).

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

Предлагается.

1. Способ определения уровня раздела фаз в каналах заключающийся в том, что устанавливают датчик, выполненный, например, в виде электропроводной проволоки, в канале по направлению силы тяжести, нагревают датчик путем пропускания тока через датчик, измеряют электрическое сопротивление датчика R, отличающийся тем, что измеряют ток I, проходящий через датчик, определяют приращение температуры датчика на участках датчика, контактирующих с паровой и жидкой фазами Δt_п=I²R/(πdLα_п), Δt_ж=I²R/(πdLα_ж).

Определяют удельное электрическое сопротивление датчика, контактирующего с жидкой ρ_ж и паровой ρ_п фазами ρ_п=ρ₀(1+βΔt_п), ρ_ж=ρ₀(1+βΔt_ж), определяют толщину парового h_п и жидкостного слоя h_ж.

h_п=(R·S-ρ_жL)/(ρ_п-ρ_ж), h_ж=L-h_п, где ρ_ж и ρ_п - удельное электрическое сопротивление датчика, находящегося в жидкой и паровой фазе соответственно; R - электрическое сопротивление датчика; I - ток через датчик; L - длина датчика; S - поперечное сечение датчика, ρ термический коэффициент сопротивления, d - диаметр датчика, ρ₀ - удельное электрическое сопротивление материала датчика при t=20°С, α_п, α_ж - коэффициенты теплоотдачи при взаимодействии датчика с паровой и жидкой фазами соответственно определяются в предварительных опытах.

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

Δt_п=I1²R1/(πdLα_п), Δt_ж=I1²R1/(πdLα_ж)

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

ρ_п=ρ₀(1+βΔt_п), ρ_ж=ρ₀(1+βΔt_ж),

определяют толщину слоя, соответствующего суммарной толщине парового и двухфазного слоя h_п1=h_п+h_дв, h_п1=(R1·S-ρ_жL)/(ρ_п1-ρж)

определяют толщину жидкостного слоя,

h_ж=L-(R1·S-ρ_жL)/(ρ_п1-ρ_ж), определяют электрическое сопротивление основного датчика, соответствующее длине h_п1 - суммарной толщине парового и двухфазного слоя, R_п1=R-ρ_жh_ж/S,

определяют толщину парового и двухфазного слоя

h_п=(R_п1·S-ρ_двh_п1))/(ρ_п-ρ_дв), h_дв=L-h_ж-h_п,

где R1 электрическое сопротивление дополнительного датчика; L - длина датчика; I1 - ток через дополнительный датчик; S - поперечное сечение дополнительного датчика; ρ_дв - удельное электрическое сопротивление датчика, находящегося в двухфазном слое, ρ_п в паровом слое, ρ_ж - в жидкостном слое; R - электрическое сопротивление основного датчика.

Технический результат, на достижение которого направлено изобретение, заключается в повышении точности определения уровня раздела паровой и жидкой фаз, что обеспечивается тем, что измеряют ток I, проходящий через основной датчик, определяют приращение температуры датчика на участках датчика, контактирующих с паровой и жидкой фазами Δt_п=I²R/(πdLα_п), Δt_ж=I²R/(πdLα_ж). Определяют удельное электрическое сопротивление датчика, находящегося в жидкой ρ_ж и паровой ρ_п фазах ρ_п=ρ₀(1+βΔt_п), ρ_ж=ρ₀(1+βΔt_ж), определяют толщину парового h_п и жидкостного слоя h_ж

h_п=(RS-ρ_жL)/(ρ_п-ρ_ж), h_ж=L-h_п, где ρ_ж и ρ_п - удельное электрическое сопротивление датчика, находящегося в жидкой ρ_ж и паровой фазе соответственно; R - электрическое сопротивление основного датчика; I - ток через датчик; L - длина датчика; S - поперечное сечение датчика, β - термический коэффициент сопротивления, d - диаметр датчика, ρ₀ - удельное электрическое сопротивление материала датчика при t=20°C, α_п, α_ж - коэффициенты теплоотдачи на поверхности датчика при взаимодействии с паровой и жидкой фазами соответственно - определяются в предварительных опытах.

В случае наличия в канале двухфазного слоя технический результат, на достижение которого направлено изобретение, заключается в повышении точности определения уровня раздела паровой, жидкой фаз и двухфазного слоя, что обеспечивается тем, что устанавливают дополнительный датчик, длина которого равна длине основного датчика, в виде электропроводной проволоки в сечение канала, в котором отсутствует двухфазный слой, нагревают дополнительный датчик путем пропускания тока, измеряют ток I1, измеряют электрическое сопротивление дополнительного датчика R1, определяют приращение температуры датчика на участках дополнительного датчика, контактирующих с паровой и жидкой фазами Δt_п=I1²R1/(πdLα_п), Δt_ж=I1²R1/(πdLα_ж), определяют удельное электрическое сопротивление дополнительного датчика, контактирующего с паровой и жидкой фазами ρ_п=ρ₀(1+βΔt_п), ρ_ж=ρ₀(1+βΔt_ж), определяют толщину слоя, соответствующего суммарной толщине парового и двухфазного слоя h_п1=h_п+h_дв,

h_п1=(R1·S-ρ_жL)/(ρ_п1-ρ_ж), определяют толщину жидкостного слоя, h_ж=L-(R1·S-ρ_жL))/(ρ_п1-ρ_ж), определяют электрическое сопротивление основного датчика, соответствующее длине h_п1 - суммарной толщине парового и двухфазного слоя, R_п1=R-ρ_жh_ж/S, определяют толщину парового и двухфазного слоя, h_п=(R_п1·S-ρ_двh_п1))/(ρ_п-ρ_дв),

h_дв=L-h_ж-h_п, где R1 - электрическое сопротивление дополнительного датчика; L - длина датчика; I1 - ток через дополнительный датчик; S - поперечное сечение дополнительного датчика; ρ_дв - удельное электрическое сопротивление датчика, находящегося в двухфазном слое, ρ_п в паровом слое, ρ_ж - в жидкостном слое; R - электрическое сопротивление основного датчика.

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

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

Способ определения уровня раздела фаз в каналах осуществляется следующим образом.

1. Устанавливают датчик, выполненный, например, в виде электропроводной проволоки, в канале по направлению силы тяжести

2. Нагревают датчик путем пропускания тока через датчик

3. Измеряют ток, проходящий через датчик

4. Измеряют электрическое сопротивление датчика R

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

Δt_п=q/α_п, Δt_п=q/α_ж

α_п, α_ж - коэффициенты теплоотдачи при нахождении датчика в паровой и жидкой фазах соответственно, кВТ/м2°С, q плотность теплового потока на поверхности датчика кВТ/м², q=I²R/П, где П обогреваемый периметр П=πdL, Δt_п=I²R/πdLα_п, Δt_ж=I²R/πdLα_ж

6. Определяют удельное электрическое сопротивление датчика, находящегося в жидкой ρ_ж и паровой ρ_п фазах - используется зависимость удельного электрического сопротивления материала датчика от температуры, ρ_т=ρ₀(1+βΔt), ρ₀ - удельное электрическое сопротивление материала датчика при t=20°С, β - температурный коэффициент сопротивления. Для нержавеющей проволоки ρ_т=0.71(1+0.00091 Δt)

7. Определяют толщину парового (газового) h_п и жидкостного слоя h_ж:h_п=(RS-ρ_жL)/(ρ_п-ρ_ж), h_ж=L-h_п, где ρ_ж и ρ_п - удельное электрическое сопротивление датчика, находящегося в жидкой ρ_ж и паровой фазе соответственно; R - электрическое сопротивление датчика; L - длина датчика; S - поперечное сечение датчика.

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

R_п=ρ_пh_п/S, R_ж=ρ_ж h_ж/S, тепловые потоки q_п=I²R_п/πdh_п, q_ж=I²R_п/πdh_п, Δt_п=I²ρ_п/πdSα_п, Δt_ж=I²ρ_ж/πdS α_ж и новые значения ρ_п, ρ_ж, h_п, h_ж. Расчет заканчивается при достижении определенной разности Δh_п=h_пn+1- h_пn≤δ, где h_пn+1, h_пn - n и n+1 значения толщины парового слоя h_n, δ = точность расчета.

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

1. Устанавливают датчик, выполненный, например, в виде электропроводной проволоки, в канале по направлению силы тяжести

2. Нагревают датчик путем пропускания тока через датчик

3. Измеряют ток I, проходящий через датчик

4. Измеряют электрическое сопротивление датчика R

5. Устанавливают дополнительный датчик в виде электропроводной проволоки в сечении канала, где отсутствует двухфазный слой

6. Нагревают датчик путем пропускания тока через датчик

7. Измеряют ток I1, проходящий через дополнительный датчик

8. Измеряют электрическое сопротивление дополнительного датчика R1

9. Определяют удельное электрическое сопротивление датчика, находящегося в жидкой ρ_ж и паровой ρ_п фазах соответственно,

Для определения ρ_ж, ρ_п определяют приращение температур на участках дополнительного датчика, контактирующих с паровой и жидкой фазами Δt_п=q/α_п, Δt_п=q/α_ж α_п, α_ж - коэффициенты теплоотдачи в паровой и жидкой фазах соответственно, кВТ/м2°С, q - плотность теплового потока на поверхности датчика кВТ/м², q=I1²R1/П, где П - обогреваемый периметр П=πdL.

10. Определяют удельное электрическое сопротивление датчика, находящегося в жидкой ρ_ж и паровой ρ_п фазах - используется зависимость удельного электрического сопротивления материала датчика от температуры, ρ_т=ρ₀(1+βΔt), ρ₀ - удельное электрическое сопротивление материала датчика при t=20°С, β - температурный коэффициент сопротивления. Для нержавеющей проволоки ρ_т=0.71(1+0.00091 Δt)

11. Определяют толщину слоя, соответствующего суммарной толщине парового и двухфазного слоя h_п1=h_п+h_дв

h_п1=(R1·S-ρ_жL)/(ρ_п1-ρ_ж), определяют толщину жидкостного слоя

h_ж=L-(R1·S-ρ_жL))/(ρ_п1-ρ_ж), определяют электрическое сопротивление основного датчика, соответствующее длине h_п1 - суммарной толщине парового и двухфазного слоя, R_п1=R-ρ_ж h_ж/S, определяют толщину парового и двухфазного слоя, h_п=(R_п1·S-ρ_двh_п1))/(ρ_п-ρ_дв), h_дв=L-h_ж-h_п, где R1 электрическое сопротивление дополнительного датчика; L - длина датчика; I1 - ток через дополнительный датчик; S - поперечное сечение дополнительного датчика; ρ_дв - удельное электрическое сопротивление датчика, находящегося в двухфазном слое, ρ_п в паровом слое, ρ_ж - в жидкостном слое; R - электрическое сопротивление основного датчика. Полученные значения h_п, h_ж, h_дв являются первым приближением. Далее значения h_п, h_ж, h_дв уточняются.

На фиг.1 показана схема размещения датчика в парогенераторе, где отсутствует двухфазный слой. В парогенерирующем устройстве 1 имеем паровой слой толщиной h_п и жидкостный слой толщиной h_ж. Через электроизолятор 2 введен датчик 3, выполненный в виде электропроводной проволоки. Длина датчика 3 равна L=h_п+h_ж. На фиг.2 показана схема размещения датчика в парогенераторе, где присутствует двухфазный слой. В парогенерирующем устройстве 1 имеем паровой слой толщиной hп, жидкостный слой толщиной h_ж и двухфазный слой толщиной h_дв. Через электроизоляторы 2 введен основной датчик 3 и дополнительный датчик 4. Датчики выполнены в виде электропроводной проволоки. Длина датчика 4 равна L=h_п1+h_ж. Паровой слой h_п1=h_дв+h_п. Пример определения уровня раздела фаз в парогенераторе АЭС при наличии двухфазного слоя. В качестве базовой длины принята длина 4 м. В качестве датчика используем нержавеющий провод диаметром 1 мм, длина 4 м. Ток, пропускаемый через основной и дополнительный датчик, I=I1=10А. Измеренное сопротивление основного датчика - R=4,05 Ом, дополнительного R1=4,15 Ом, S=πd²/4=0,78510^-6м², ρ₀=0,71 Ом мм²/м. Определим удельное электрическое сопротивление датчика, находящегося в жидкой ρ_ж и паровой фазе. Для этого определим среднюю плотность теплового потока на поверхности дополнительного датчика q=W/П=I1² R/πd L=100·4,15/πd L=100 4,15/π·1·4·10^-3=33,04 кВт/м², W, П - электрическая мощность, подведенная к проводнику и обогреваемый периметр проводника соответственно, d - диаметр проводника. Определяем приращение температуры Δt_п, Δt_ж на частях провода, находящихся в контакте с паровой и жидкой фазой., Δt_п=q/α_п, Δt_ж=q/α_ж, α_п, α_ж - коэффициенты теплоотдачи провода с паровой и жидкой фазах соответственно, кВТ/м²°С. Определяются либо в предварительных опытах, либо из литературных источников. Примем следующие значения α_п, α_ж.

α_п=100 Вт/м²°С, α_ж=5000 Вт/м²°С, Δt_п=q/α_п=33040/100=330,4°С, Δt_ж=q/α_ж=33040/5000=6,608°С, ρ_п=0,71(1+0,00091·330.4)=0,9237, ρ_п1=ρ_п, ρ_ж=0,71(1+0,00091·6,608)=0,7143, h_п1=(4,15·0,785-0,714·4)/(0,916-0,714)=1,98886, h_ж=4-1,98886=2,011

Определяем электрическое сопротивление основного датчика, соответствующее длине h_п1 - суммарной толщине парового и двухфазного слоя

R_п1=R-ρ_ж h_ж/S=4,05-0,714 2,011/0,785=2,2208

h_п=(R_п1·S-ρ_дв h_п1))/(ρ_п1-ρ_дв)=(2,2208 0,785-0,715 1,98886)/(0,916-0,715)=1,5996, h_дв=4-2,011-1,5996=0,3894,

где ρ_дв=0,71(1+0,00091·1,612)=0,715, Δt_дв=q/α_дв=1,612. α_дв - коэффициент теплоотдачи в двухфазном слое - принят равным 20000 Вт/м²°С.

Второе приближение - определяем электрическое сопротивление датчика (провода), находящегося в паровой и жидкой фазах, R_п=ρ_п h_п/S=0,9237 1,98886/0,785=2,3403, R_ж=1,809. Определяем тепловые потоки и, соответственно, приращение температур на части датчика, находящегося в паровой и жидкой фазах q_п=100 2,3403/3,14·1·1,9886=37149,5 Вт/м², Δt_п=371,49°С, Δt_ж=5,73°С. ρ_п=0,947, ρ_ж=0,7137

h_п1=(4,15·0,785-0,7137·4)/(0,947-0,7137)=1,727

h_ж=4-1,7274=2,273, R_п1=R-ρ_ж h_ж/S=4,05-0,7137 2,273/0,785=1,933,

h_п=(Rп1·S-ρ_двh_п1))/(ρ_п1-ρ_дв)=(1,933·0,785-0,711 1,727)7(0,947-0,7137)=1,2396, h_дв=4-2,273-1,273=0,454, где ρ_дв=0,71(1+0,00091·1,447)=0,711, Δt_дв=q/α_дв=1,447. Разница между первым и вторым приближением толщины β=(0,454-0,3894)-100/0.454=14,2%.

Третье приближение двухфазного слоя.

Определяем электрическое сопротивление датчика (провода), находящегося в паровой и жидкой фазах

R_п=ρ_п h_п1/S=0,947 1,727/0,785=2,0834, R_ж=0,7137 2,273/0,785=2,066, q_п=100 2,0834/(3,14 1 1,727)=38419, Δt_п=384,19°С, Δt_ж=5,789°С.

ρ_п=0,955, ρ_ж=0,7136. Определяем тепловые потоки и приращение температур на части датчика, находящегося в паровой и жидкой фазах.

q_п=100 2,3403/3,14·1·1,9886=37149,5 Вт/м², Δt_п=371,49°С, Δt_ж=5,73°С.

ρ_п=0,947, ρ_ж=0,7137. h_п1=(4,15·0,785-0,7136·4)/(0,955-0,7136)=1,6706

h_ж=4-1,7274=2,3294Rп1=R-ρж hж/S=4,05-0,7136 2,3294/0,785=1932,

R_дв=0,711 0,454/0,785=0,4112, q_дв=100 0,4112/(3,14 1 0,454)=28845. h_п=(R_п1·S-ρ_дв h_п1))/(ρ_п1-ρ_дв)=(1,9324·0,785-0,711 1,6706)/(0,955-0,7109)=1,352, h_дв=4-2,2304-1,35=0,4196,

где ρ_дв=0,71(1+0,00091·1,447)=0,711, Δt_дв=q/α_дв=1,442. Разница между вторым и третьим приближением толщины двухфазного слоя составляет β=(0,454-0,4196)-100/0,454=7,5%.

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