CC BY
43
ЭЛЕКТРОФИЗИКА, ЭЛЕКТРОФИЗИЧЕСКИЕ
УСТАНОВКИ
УДК 536.75
ЗАВИСИМОСТЬ МЕРЫ КУЛЬБАКА ФЛУКТУАЦИЙ НАПРЯЖЕНИЯ НА ЭЛЕКТРОЛИТИЧЕСКИХ ЯЧЕЙКАХ ОТ МЕТЕОРОЛОГИЧЕСКИХ ФАКТОРОВ
А.Н. Морозов
МГТУ им. Н.Э. Баумана, Москва, Российская Федерация e-mail: amor59@mail.ru
Приведены результаты обработки долговременных записей флуктуаций напряжения на электролитических ячейках двух независимых экспериментальных установок в двух сериях экспериментов. Установлено наличие корреляции значений мер Кульбака флуктуаций напряжения на электролитических ячейках и следующих метеорологических факторов: температуры воздуха на улице; температуры точки росы; скорости ветра; относительной и абсолютной влажности воздуха; давления насыщенного пара. Показано, что значения мер Кульбака запаздывают относительно метеорологических факторов на периоды времени от 18 до 51 ч. Наибольшая корреляция наблюдается при воздействии вариаций температуры точки росы, абсолютной влажности и давления насыщенного пара. Установлено, что указанное влияние не может быть объяснено температурным воздействием на экспериментальные установки.
Ключевые слова: флуктуации напряжения, электролитическая ячейка, мера Кульбака, температура, точка росы, относительная влажность, абсолютная влажность, давление насыщенного пара, скорость ветра, коэффициент корреляции.
CORRELATION BETWEEN KULLBACK'S MEASURE OF ELECTROLYTIC CELL VOLTAGE FLUCTUATIONS AND METEOROLOGICAL FACTORS
A.N. Morozov
Bauman Moscow State Technical University, Moscow, Russian Federation e-mail: amor59@mail.ru
The article discusses the results of the persistent data processing of electrolytic cell voltage fluctuations conducted in two separate experimental setups during two sets of experiments. It presents the determined correlation between the electrolytic cell voltage fluctuations and the following meteorological factors: outside air temperature, dew-point temperature, wind speed, relative humidity, absolute humidity, and vapor pressure. It is found out that Kullback's measure parameters lag behind the meteorological factors up to 18-51 hours. The maximum correlation is observed from the exposure to dew-point temperature variations, absolute humidity, and vapor pressure. This correlation can not be explained by the exposure of the experimental setups to temperature variations.
Keywords: voltage fluctuation, electrolytic cell, Kullback's measure, temperature, dew-point, relative humidity, absolute humidity, vapor pressure, wind speed, correlation factor.
Результаты измерений меры Кульбака флуктуаций напряжения на электролитических ячейках, полученные при долговременных экспериментах, выполненных в трех экспериментальных сериях (с 10 декабря 2007 г. по 16 марта 2008 г., с 19 мая 2009 г. по 1 сентября 2009 г. и с 20 марта 2011 г. по 18 декабря 2012 г.) приведены в работах [1,2]. Установлено, что меры Кульбака, рассчитанные по данным измерений флуктуаций напряжения на двух независимых экспериментальных установках, коррелируют между собой.
Цель настоящей работы — установление зависимости мер Кульбака от внешних метеорологических факторов. При анализе будут использованы результаты двух серий долговременных экспериментов: ранее проведенных в период с 20 марта 2011 г. по 18 декабря 2012 г. (15 360 ч, или 640 сут.) и новых — с 25 сентября 2013 г. по 29 декабря 2014 г. (11 064 ч, или 461 сут.). Общая продолжительность экспериментов составила 26424 ч, или 1101 сут.
В первой серии экспериментов на электролитические ячейки подавалось постоянное напряжение 1,2... 6,8 В, во втором — электролитические ячейки не подключались к внешнему источнику напряжения. Объем электролита в тонких каналах электролитических ячеек в двух сериях экспериментов был примерно равен V = 2 • 10-14 м3. Конструкция установок и методика расчета мер Кульбака флуктуаций напряжения на электролитических ячейках описаны в работах [3-5].
По результатам измерений флуктуаций напряжения со скважностью 1 мин осуществлялись вычисления дисперсий (Д1, Д2) и мер Кульбака (Н1, Н2) для сигналов, зарегистрированных на двух установках. Одновременно проводилось измерение температуры на улице (Т1) и температуры внутри пассивного термостата (Т2), где размещались установки. Временные ряды далее усреднялись и прорежались на периоде, равном 1 ч.
Мера Кульбака для описания объема информации была введена в работах [6, 7], в настоящее время она широко используется в теории передачи информации и для описания неравновесных термодинамических систем [8, 9]. В работах [1, 2] она применена как параметр, характеризующий флуктуации напряжения на электролитических ячейках.
Записанные со скважностью 1 ч указанные выше ряды подвергались фильтрации в соответствии с описанной в работе [2] методикой, что обеспечивало вычитание сигналов с периодами менее 48 ч (2 сут.) и более 1200 ч (50 сут.). Затем методом скользящей регрессии из полученных значений вариаций мер Кульбака вычитались вариации дисперсии флуктуаций напряжения и температуры установок. Указанная процедура позволяла очистить полученные значения мер Кульбака от их возможной зависимости от температуры установок. Отметим
0,10 - V у
0,15-1-1-1-1-1-
-3000 -2000 -1000 0 1000 2000 t, ч
R(Hh Td)
0,15
-0,20 -
-0,25 - '
-0,30 -1-1-1-1-1-
-3000 -2000 -1000 0 1000 2000 ч
б
Рис. 1. Зависимость коэффициента корреляции значений мер Кульбака для двух независимых установок от времени (а) и зависимость коэффициентов корреляции Я(И1 (1) и Я(Й2 (2) значений мер Кульбака и вариаций температуры точки росы от времени (б)
также, что в целях подавления возможных мешающих факторов экспериментальные установки были экранированы от электромагнитных наводок, а электролитические ячейки — герметизировались и теплоизолировались от окружающей среды.
Зависимость усредненного для двух экспериментальных серий коэффициента корреляции Я(Н1, Н2) мер Кульбака (Н1 и Н2) от времени представлена на рис. 1, а. Согласно приведенной кривой, наблюдается небольшая корреляция Я(Н1 ,Н2) = 0,185 ± 0,053 мер Кульбака (Н1 и ) для двух независимых установок при практически нулевом сдвиге этих временных рядов. Дисперсия ошибки а# = 0,053 значения корреляции Я(Н1, Н2) вычислялась по формуле, полученной в работе [10]. Из приведенных результатов следует, что отношению сигнал/шум для корреляции Я(Й1 ,Н2) равно 3,5. В соответствии с распределением Стьюдента вероятность наличия корреляции значений мер Кульбака Н1 и Н2 между двумя независимыми установками составляет Р = 0,9995 [11].
Для выявления причины возникновения взаимной корреляции мер Кульбака и их зависимости от метеорологических факторов использовались данные по температуре воздуха Т, температуре точки росы ТА, скорости ветра V, относительной влажности воздуха Як и атмосферному давлению Р, соответствующие двум сериям экспериментов. Перечисленные данные взяты с сайта "Погода и климат" (www.pogodaiklimat.ru) для метеостанции, расположенной в Москве на территории ВДНХ (индекс WMO: 27612). Кроме того, использованы данные о потоке космических нейтронов К1 (http://cr0.izmiran.rssi.ru/mosc/main.htm) и солнечной активности Бг (ftp://ftp.ngdc.noaa.gov/STP/SOLAR_DATA/SOLAR_RADIO/FLUX/ Penticton_Observed/daily/DAILYPLT.OBS).
Зависимость коэффициента корреляции Я(Н,,,Тё), г = 1, 2, вариаций значений мер Кульбака (Н1 и Н2) для двух установок и вариаций температуры точки росы на улице ТА, усредненных для двух серий экспериментов, от времени приведена на рис. 1, б. Наблюдается отрицательная корреляция значений мер Кульбака Н,, г = 1, 2, и вариаций температуры точки росы ТА Я(Н1,Тд) = -0,256 ± 0,068; Я(Н2,ТА) = -0,241 ± 0,059 с запаздыванием значений мер Кульбака относительно температуры точки росы на 27 ч для первой экспериментальной установки и на 33 ч для второй. Вероятность наличия корреляции для двух установок составляет Р1 = 0,9998 и Р2 = 0,99996. Согласно приведенным кривым, для двух независимых установок наблюдаются близкие значения корреляционных функций Я(Н1,Т<£) и Я(Н2,Тд), что свидетельствует об одинаковом характере воздействия указанного фактора на флуктуации напряжения в малых объемах электролитических ячеек.
Зависимость коэффициента корреляции среднего значения мер Кульбака Н12 = + Н^ у/2 для двух установок и вариаций температуры точки росы ТА показана рис. 2, а: Я (н1>2, Т^ = —0,317±0,070 при запаздывании на 30ч; отношение сигнал/шум 4,5; вероятность указанной корреляции Р = 0,999994.
Для определения зависимости от времени корреляции средних значений мер Кульбака Н12 и вариаций температуры точки росы ТА по методике, изложенной в работе [2], вычислялось значение Бнт = Н12Т А, усредненное за период времени 720 ч (30 сут.). Для первой серии экспериментов приведена зависимость произведения Бнт при сдвиге значений меры Кульбака Н12 и температуры точки росы ТА на величину запаздывания по времени, равную 30 ч (рис. 2, б). Произведение этих величин во всем диапазоне периода проведения эксперимента имеет отрицательное значение.
Я №,2, т2)
Td,° С
20 10 0 -10 -20 -30
0 2000 4000 6000 8000 10000 12000 14000 ч
в
Рис. 2. Зависимость коэффициента корреляции средних значений мер Кульбака и вариаций температуры точки росы от времени (а), зависимость произведения средних значений мер Кульбака и вариаций температуры точки росы для первой серии экспериментов от времени (б), зависимость температуры точки росы от времени за период наблюдений с 20 марта 2011 г. по 18 декабря 2012 г. (в)
Усредненное за весь период наблюдения в первой серии экспериментов произведение средних значений мер Кульбака Н\,2 и вариаций температуры точки росы Т6 равно Бит = -2,36 • 10-7 К.
Явно выраженные минимумы кривой, приведенной на рис. 2, б, наблюдаются при значениях времени от начала временного ряда 1314, 5 556, 8 025, 9 942, 10 374, 12 153 и 14 100 ч, что соответствует 13 мая и 5 ноября 2011г., 16 февраля, 6 и 25 мая, 7 августа и 22 октября
2012 г. Указанным датам, за исключением 16 февраля 2011 г. и 7 августа 2012г., соответствуют моменты времени колебания температуры точки росы Т^ в диапазоне значений 0... +100С (рис. 2, в).
Зависимость произведения средних значений мер Кульбака и вариаций температуры точки росы от времени для второй серии экспериментов, аналогичная зависимости, приведенной на рис. 2, б, представлена на рис.3,а. Среднее значение произведения £>нт = -3,31 х х 10-7 К, что совпадает со значением, полученным для первой серии экспериментов. Это указывает на независимость воздействия на флуктуации напряжения от приложенного к электролитическим ячейкам внешнего напряжения, так как во второй серии экспериментов, в отличие от первой, оно отсутствовало.
Минимумы кривой, приведенной на рис. 3, а, наблюдаются для значений времени от начала временного ряда 5 247 и 9 425 ч, что соответствует 1 мая и 22 октября 2014 г. Так же, как и для первой серии экспериментов, в эти даты наблюдались колебания температуры точки росы Т^ в диапазоне значений 0... +100С (рис. 3, б).
Рис. 3. Зависимость произведения средних значений мер Кульбака и вариаций температуры точки росы от времени для второй серии экспериментов (а) и зависимость температуры точки росы от времени за период наблюдения с 25 сентября 2013 г. по 29 декабря 2014 г. (б)
Были рассчитаны корреляции средних значений мер Кульбака Н 12 от других параметров, указанных выше. Установлена несколько меньшая, чем для температуры точки росы, корреляция с температурой на улице по данным метеостанции, расположенной в Москве Т (Д(Нх,2,Т) = -0,236 ± 0,071 при запаздывании 27 ч, Р = 0,9995), и температурой Т, измеренной на улице рядом с установками (Я(Н\,2,Т1) = -0,208 ± 0,069 при запаздывании 18 ч, Р = 0,997).
Согласно кривой, приведенной на рис. 4, а, наблюдается небольшая корреляция вариаций средних значений мер Кульбака и вариаций скорости ветра V на улице Я(Д\,2, У) = 0,195 ±0,058 с запаздыванием средних значений мер Кульбака Н1}2 относительно скорости ветра V на улице на 51 ч. Вероятность наличия корреляции равна Р = 0,9993.
Наблюдается корреляция меры Кульбака Н1}2 с измеренной относительной влажностью на улице ЯН Я(Д\,2, ЯН) = -0,195 ± 0,064 при запаздывании 33 ч (Р = 0,998) и с рассчитанной в соответствии с формулами, приведенными в работах [12, 13], абсолютной влажностью Яо Я(Н1,2, Яо) = -0,352 ± 0,070 при запаздывании 30 ч (рис.4,б). В последнем случае вероятность корреляции составляет
Рис. 4. Зависимости коэффициента корреляции средних значений мер Кульбака и вариаций скорости ветра (а) и абсолютной влажности (б) от времени
Рис. 5. Зависимость коэффициентов корреляции значений мер Кульбака и температуры установок от внешних метеорологических факторов
Р > 0,999999. Аналогичный результат получается и в случае нахождения корреляции с давлением насыщенного водяного пара Ро Д(Н\,2,Ро) = -0,351 ± 0,071 при запаздывании 30ч. Указанные корреляции являются самыми сильными из всех обнаруженных при корреляционном анализе.
Корреляционный анализ показал отсутствие значимых корреляций вариаций средних значений мер Кульбака НН1,2 с вариациями температуры установок Т2, атмосферного давления Р, потока космических нейтронов КI и солнечной активности Бг.
Зависимость абсолютных значений коэффициентов корреляции мер Кульбака Д(Н\,2,Х) и температуры установок Я(Т2,Х) от
внешних метеорологических факторов X (Яо,Тб,,Т,Т1 ,ЯК) приведена на рис. 5. Явно прослеживается обратная зависимость коэффи-
циентов корреляции
R(H 1,2Д) и R(T2,Х)
, которая указывает на
увеличение влияния метеорологического фактора на вариации мер Кульбака Н12 при уменьшении указанного влияния на температуру установок Т2. С учетом этого можно сделать вывод о том, что механизм влияния внешних метеорологических процессов на изменение значения меры Кульбака не связан с непосредственным температурным воздействием на экспериментальные установки. Сведем полученные выше результаты в таблицу.
Согласно данным, приведенным в таблице, запаздывание значений мер Кульбака относительно скорости ветра несколько больше (51ч) по сравнению с запаздыванием относительно температуры на улице (27 ч), а корреляция имеет противоположный знак, что объясняется зависимостью процесса охлаждения атмосферного воздуха от скорости ветра, приносящего холодные атмосферные массы.
Меньшее запаздывание температуры на улице, измеренной в непосредственной близости от установок, расположенных на территории МГТУ им. Н.Э. Баумана (18 ч), по отношению к данным метеостанции, расположенной в Москве на территории ВДНХ (27 ч), связано
Полученные результаты
Коэффициент корреляции Отношение сигнал/шум Вероятность корреляции Опережение (запаздывание), ч
R(ffbff2) = 0,185 ± 0,053 3,5 0,9995 0
Д(ЯЬ TTd) = -0,256 ± 0,068 3,8 0,9998 -27
R(H2, Td) = -0,241 ± 0,059 4,1 0,99996 -33
R(Hlt2,Td) = -0,317 ± 0,070 4,5 0,999994 -30
R(Hi,2, TT) = -0,236 ± 0,071 3,3 0,999 -27
R(Hi,2, Ti) = -0,208 ± 0,069 3,0 0,997 -18
R(ff1l2,T) = 0,195 ± 0,058 3,4 0,9993 -51
R(ff1>2, Rh) = -0,198 ± 0,064 3,1 0,998 -33
R(Hlt2,Po) = -0,351 ± 0,071 4,9 >0,999999 -30
R(H 1>2, Ro) = -0,352 ± 0,070 5,0 >0,999999 -30
с локальным прогревом воздушных масс в черте плотной городской застройки.
Полученные результаты показывают наличие корреляции мер Кульбака флуктуаций напряжения на электролитической ячейке и метеорологических факторов. Максимальная корреляция мер Кульбака наблюдается с вариациями абсолютной влажности воздуха (давлением насыщенного пара). Эта корреляция несколько выше в сравнении с аналогичным значением для температуры точки росы, которая также немного больше корреляции с температурой воздуха по данным метеостанции, находящейся в Москве. Из этого следует, что наибольшее влияние на экспериментальные установки оказывают вариации абсолютной влажности атмосферного воздуха (давления насыщенного пара), но конкретный механизм этого влияния не установлен.
ЛИТЕРАТУРА
1. Морозов ^.^Предварительные результаты измерений меры Кульбака флуктуаций напряжения на электролитической ячейке // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2011. № 2. С. 16-24.
2. Морозов А.Н.Применение меры Кульбака для оценки долговременных изменений флуктуаций напряжения на электролитической ячейке // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2013. № 3. С. 52-61.
3. Морозов А.Н.Необратимые процессы и броуновское движение: Физико-технические проблемы. М.: Изд-во МГТУ им. Н.Э. Баумана, 1997. 332 с.
4. Коротаев С.М., Морозов А.Н., Сердюк В.О., Сорокин М.О. Проявление макроскопической нелокальности в некоторых естественных диссипативных процессах // Известия вузов. Физика. 2002. № 5. С. 3-14.
5. Experimental study of macroscopic nonlocality of large-scale natural dissipative processes / S.M. Korotaev, A.N. Morozov, V.O. Serdyuk, V.A. Machinin, J.V. Gorokhov // NeuroQuantology. 2005. Iss. 4. P. 275-294.
6. KullbackS., Leibler R.A. On information and sufficiency // Ann. Math. Statist. 1951. Vol. 22. P. 79-86.
7. Кульбак С. Теория информации и статистика; пер. с англ. М.: Наука, 1967. 408 с.
8. ЗариповР.Г. Новые меры и методы в теории информации. Казань: Изд-во Казан. гос. тех. ун-та, 2005. 364 с.
9. Климонтович ЮЛ.Турбулентное движение и структура хаоса: Новый подход к статистической теории открытых систем. М.: Наука, 1990. 320 с.
10. Макс Ж. Методы и техника обработки сигналов при физических измерениях; пер. с франц. Т. 1. М.: Мир, 1983. 312 с.
11. Янке Е., Эмде Ф., Леш Ф. Специальные функции. Формулы, графики, таблицы; пер. с нем. М.: Наука, 1977. 344 с.
12. Sonntag D. Advancements in the field of hygrometry // Meteorol. Z., N.F. 1994. Vol. 3. P. 51-66.
13. Murphy D.M., Koop T. Review of the vapor pressures of ice and supercooled water for atmospheric applications // Quart. J. Royal Met. Soc. 2005. Vol. 31. P. 1539-1565.
REFERENCES
[1] Morozov A.N. Preliminary results of recording the Kullback measure of voltage fluctuations on electrolytic Cell. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Estestv. Nauki [Herald of the Bauman Moscow State Tech. Univ., Nat. Sci.], 2011, no. 2, pp. 16-24 (in Russ.).
[2] Morozov A.N. Application of the Kullback measure for estimation of long-term variations in voltage fluctuations on the electrolytic cell. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Estestv. Nauki [Herald of the Bauman Moscow State Tech. Univ., Nat. Sci.], 2013, no. 3, pp. 52-61 (in Russ.).
[3] Morozov A.N. Neobratimye protsessy i brounovskoe dvizhenie: Fiziko-tekhnicheskie problem [Irreversible Processes and Brownian Motion: Physical and Technological Problems]. Moscow, MGTU im. N.E. Baumana Publ., 1997. 332 p.
[4] Korotaev S.M., Morozov A.N., Serdyuk V.O., Sorokin M.O. Experience of Macroscopic Non-Locality in Some Natural Dissipative Processes. Izv. Vyssh. Uchebn. Zaved., Fizika [Proc. Univ., Physics], 2002, no. 5, pp. 3-14 (in Russ.).
[5] Korotaev S.M., Morozov A.N., Serdyuk V.O., Machinin V.A., Gorokhov J.V., Machinin V.A. Experimental study of macroscopic nonlocality of large-scale natural dissipative processes. Neuro Quantology, 2005, iss. 4, pp. 275-294.
[6] Kullback S., Leibler R.A. On information and sufficiency. Ann. Math. Statist., 1951, vol. 22, pp. 79-86.
[7] Kullback S. Information Theory and Statistics. New York, John Wiley & Sons, Inc.,
1959, 409 p. [Russ. ed.: Kullback S. Teoriya informatsii i statistika. Moscow, Nauka Publ., 1967. 408 p.].
[8] Zaripov R.G. Novye mery i metody v teorii informatsii [New Measures and Techniques in Information Theory]. Kazan', Kazan. Gos. Tekh. Univ. Publ., 2005. 364 p.
[9] Klimontovich Yu.L. Turbulentnoe dvizhenie i struktura khaosa: Novyy podkhod k statisticheskoy teorii otkrytykh system [Turbulent Motion and Structure of Chaos: A New Approach to the Statistical Theory of Open Systems]. Moscow, Nauka Publ., 1990. 320 p.
[10] Max J. Metody i tekhnika obrabotki signalov pri fizicheskikh izmereniyakh [Fr. ed.: Mothodes et techniques de traitement du signal et applications aux mesures physiques. Paris, Masson, 238p. Eng. ed.: Methods and Techniques of Signal Processing in Physical Measurements]. Vol. 1. Moscow, Mir Publ., 1983. 312 p.
[11] Jahnke E., Emde F., Lesh F. Spetsial'nye funktsii [Ger. ed.: Jahnke E., Emde F., Losch F. Tafeln Höherer Funktionen. Stuttgart: B.G. Teubner Verlagsgesellschaft,
1960. Eng. ed.: Tables of higher functions]. Moscow, Nauka Publ., 1977. 344 p.
[12] Sonntag D. Advancements in the field of hygrometry. Meteorol. Z., N.F., 1994, vol. 3, pp. 51-66.
[13] Murphy D.M., Koop T. Review of the vapor pressures of ice and supercooled water for atmospheric applications. Quart. J. Royal Met. Soc., 2005, vol. 31, pp. 1539-1565.
Статья поступила в редакцию 27.01.2015
Морозов Андрей Николаевич — д-р физ.-мат. наук, профессор, заведующий кафедрой "Физика" МГТУ им. Н.Э. Баумана. Автор более 200 научных работ в области прецизионных измерений и физической кинетики.
МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.
Morozov A.N. — Dr. Sci. (Phys.-Math.), professor, head of "Physics" department of the Bauman Moscow State Technical University. Author of more than 200 publications in the field of high precision measuring systems and physical kinetics theory. Bauman Moscow State Technical University, 2-ya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.
Просьба ссылаться на эту статью следующим образом:
Морозов А.Н. Зависимость меры Кульбака флуктуаций напряжения на электролитических ячейках от метеорологических факторов // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2015. № 3. C. 47-57. Please cite this article in English as:
Morozov A.N. Correlation between Kullback's measure of electrolytic cell voltage fluctuations and meteorological factors. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Estestv. Nauki [Herald of the Bauman Moscow State Tech. Univ., Nat. Sci.], 2015, no. 3, pp. 47-57.
