Наукастинг появления геоэффективных потоков протонов в околоземном космическом пространстве



Наукастинг появления геоэффективных потоков протонов в околоземном космическом пространстве
Наукастинг появления геоэффективных потоков протонов в околоземном космическом пространстве
Наукастинг появления геоэффективных потоков протонов в околоземном космическом пространстве
Наукастинг появления геоэффективных потоков протонов в околоземном космическом пространстве
Наукастинг появления геоэффективных потоков протонов в околоземном космическом пространстве
Наукастинг появления геоэффективных потоков протонов в околоземном космическом пространстве
Наукастинг появления геоэффективных потоков протонов в околоземном космическом пространстве
Наукастинг появления геоэффективных потоков протонов в околоземном космическом пространстве

 


Владельцы патента RU 2485557:

Буров Вячеслав Анатольевич (RU)
Очелков Юрий Павлович (RU)

Использование: для прогноза космической погоды, определяемой вспышечной активностью Солнца. Сущность: заключается в том, что осуществляют непрерывный сбор и анализ информации о потоке рентгеновского излучения Солнца в диапазоне длин волн 0.1-0.8 нм, при этом в процессе формирования прогноза в режиме реального времени осуществляют поиск всплесков рентгеновского излучения с интенсивностью, превосходящей значение 5*10-5 Вт/м2, и для каждого такого всплеска определяют параметр «протонности» всплеска по заданному соотношению и в зависимости от его значения и гелиодолготы всплеска определяют вероятность наблюдения в околоземном космическом пространстве солнечного протонного события с энергией протонов больше 30 МэВ с интенсивностью в максимуме события выше заданных уровней. Технический результат: обеспечение возможности более высокой заблаговременности прогноза, а также обеспечение возможности лучшего учета временной эволюции всплеска и возможность прогноза появления после всплесков протонных событий разного уровня интенсивностей. 8 ил.

 

Изобретение относится к области прогноза космической погоды, определяемой вспышечной активностью Солнца, и может быть использовано для прогноза геоэффективных последствий солнечных вспышек, в частности явлений нарушения коротковолновой радиосвязи, ухудшение определения местоположения по данным ГЛОНАСС/GPS навигации; повышение радиационной опасности для экипажей и пассажиров высотных самолетов с трассами полета в полярных областях, а также сбоям в работе бортовых космических приборов и возрастанию опасности радиационного поражения экипажей пилотируемых космических аппаратов.

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

Прогноз осуществляется по данным наблюдения рентгеновского излучения Солнца в диапазоне длин волн 0.1-0.8 нм на космическом аппарате (КА) GOES с минутным временным разрешением. Используются данные с ежеминутно обновляемого файла по адресу http://www.swpc.noaa.gov/ftpdir/lists/xray/Gp_xr_1m.txt.

Способ включает в себя решение двух задач: мониторинг всплесков рентгеновского излучения с интенсивностью больше 5*10-5 Вт/м2 (всплески класса больше М5 по общепринятой классификации всплесков) и определение их параметров; прогноз интенсивности в максимуме солнечного протонного события с энергией протонов выше 30 МэВ, сопровождающего всплеск в соответствии с установленным решающим правилом. Под прогнозом понимается определение вероятности наступления после всплеска события СПС с энергиями протонов больше 30 МэВ с интенсивностями в максимуме больше чем 0.5 част/см2сср, 10 част/см2сср, 100 част/см2сср и 1000 част/см2сср.

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

1. A.Belov "Properties of Solar X-ray flares and proton event forecasting". Advances in Space Research, Vol.43, Issue 4, 2009, p.467-473 [1].

2. M.Laurenza, E.W.Cliver, J.Hewitt, M.Storini, A.G.Ling, C.C.Balch, M.L.Kaiser "A technique for short-term warning of solar energetic particle events based on flare location, flare size, and evidence of particle escape". Space Weather, Vol.7, S04008, doi:10.1029/2007SW000379, 2009 [2].

Способ прогноза, изложенный в работе [1], в настоящее время реализован в Институте земного магнетизма и распространения радиоволн РАН (ИЗМИРАН). Существенное отличие предлагаемого способа прогноза - использование (дополнительно к таким параметрам, как величина потока рентгеновского излучения в максимуме всплеска (Jx,m) и гелиодолгота всплеска) нового параметра - Tr,1 -средней ширины временного профиля рентгеновского всплеска в логарифмических координатах по интенсивности на фазе роста потока рентгеновского излучения от значения в 0.1 Jx,m до Jx,m. Использование этого параметра существенно улучшает оправдываемость прогноза. Кроме того, предлагаемый метод позволяет обеспечить непрерывный круглосуточный прогноз в реальном времени, т.е. в режиме наукастинга.

Способ прогноза, изложенный в работе [2], предполагается использовать в США в Центре прогноза космической погоды (SWPC NOAA). Существенное отличие предлагаемого здесь способа - использование других параметров всплеска для прогноза и, соответственно, другого решающего правила. В [2] в качестве основного параметра, характеризующего «протонную опасность» (вероятность наблюдения геоэффективного солнечного протонного события после всплеска), предлагается использовать интегральный поток рентгеновского излучения на фазе развития всплеска от 1/3 Jx,m на фазе роста до 1/3 Jx,m на фазе спада. Нами предлагается использование других параметров - поток рентгеновского излучения в максимуме, который характеризует масштаб события, и параметра, характеризующий временную эволюцию всплеска на фазе роста - среднюю ширину временного профиля рентгеновского всплеска в логарифмических координатах по интенсивности на фазе роста потока рентгеновского излучения от значения в 0.1 Jx,m до Jx,m. Впервые предлагается такой параметр всплеска, характеризующий его «протонную опасность» как: Р=Tr,l+0.35Lg(Jx,m*103).

Преимущества предлагаемого способа следующие: более высокая заблаговременность прогноза (используются данные только для фазы роста всплеска до максимума), лучший учет временной эволюции всплеска (предлагаемый временной параметр описывает также форму всплеска), возможность прогноза появления после всплесков протонных событий разного уровня интенсивностей: больше чем 0.5 част/см2сср, 10 част/см2сср, 100 част/см2сср и 1000 част/см2сср (в [2] предлагается только способ прогноза появления после всплесков протонных событий с энергией протонов больше 10 МэВ с интенсивностью выше, чем 10 част/см2сср).

Реализация предлагаемого способа

Предлагаемый метод реализован, во-первых, установлением решающего правила прогноза, во-вторых разработкой алгоритма непрерывного поиска всплесков рентгеновского излучения с интенсивностью больше 5*10-5 Вт/м2 (всплески класса М5 и больше по общепринятой классификации всплесков) и определения их параметров, используемых в решающем правиле. Решение этих задач позволило разработать программу непрерывного краткосрочного прогноза (наукастинга) появления в околоземном космическом пространстве геоэффективных потоков протонов.

1. Определение решающего правила прогноза

Метод и проблемы

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

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

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

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

Другим важным вопросом, который требует изучения, является вопрос о значительном разбросе (дисперсии) интенсивностей СПС от вспышек с одинаковыми значениями интенсивности электромагнитного излучения. Этот вопрос имеет отношение к взаимосвязи между количеством ускоренных электронов, ответственных за генерацию излучений, и протонов. Особенно большая дисперсия имеет место для тепловых излучений, например, таких как мягкое рентгеновское излучение (наблюдаемая дисперсия достигает 4-х порядков величины).

Существование большой дисперсии значительно затрудняет изучение взаимосвязи интенсивности СПС и мягких рентгеновских всплесков. Причины такой дисперсии в настоящее время не ясны. Одним из важных факторов является тепловая природа мягких рентгеновских всплесков, в то время как генерация протонов - нетепловой процесс. Интенсивность мягкого рентгеновского излучения определяется температурой, до которой нагревается область излучения, и масштабами этой области. Связь интенсивности СПС с нетепловыми излучениями солнечных вспышек, такими как микроволновое излучение (механизм - магнитотормозное излучение нетепловых электронов, ускоренных во вспышках) и жесткое рентгеновское излучение (механизм - тормозное излучение нетепловых электронов) значительно лучше, то есть дисперсия интенсивности СПС значительно меньше (хотя также велика). Однако для практических целей прогноза интенсивности СПС использование мягкого рентгеновского излучения обладает рядом преимуществ. Во-первых, данные о мягком рентгеновском излучении в настоящее время доступны в реальном времени со спутника GOES с минутным временным разрешением (и круглосуточно), во-вторых, имеется достаточно большой объем однородных данных наблюдений мягкого рентгеновского излучения за несколько циклов солнечной активности. Жесткое рентгеновское излучение наблюдалось к настоящему времени только в отдельных научных экспериментах. Данные по нему в оперативном режиме недоступны. Данные по микроволновому излучению основываются на наземных наблюдениях сети станций. Сбор и получение данных в оперативном режиме в этом случае затруднен.

Выбор параметров для прогноза

В качестве основного параметра, используемого для прогноза СПС по излучению вспышек, естественно выбрать поток рентгеновского излучения в максимуме события, так как эта величина определяет масштаб солнечного события («синдром большой вспышки» по Kahler S.W. The role of the big flare syndrome in correlations of solar energetic proton fluxes and associated microwave burst parameters, J.Geophys.Res.,1982, 87, 3439-3448). Эта величина изменяется в диапазоне от 10-7 Вт/м2 до 10-3 Вт/м2.

Однако использование только этого параметра делает качество прогноза очень низким.

В ряде работ отмечалась зависимость вероятности наблюдения СПС после всплеска электромагниного излучения вспышек от временных параметров всплеска (см., например: Croom, D.L Solar Phys.1971, 19, 171-185; E.T.Sarris, S.D., S.D.Shawhan. Solar Physics.1973, v.28, p.519-532; Cane et al. Astrophys.J., 1986, 301, 448-459). Однако количественных результатов для такой зависимости получено не было. В данном методе впервые предлагается в качестве временного параметра использовать параметр Tr,1, характеризующий временное развитие всплеска на фазе роста - среднюю ширину логарифмического временного профиля рентгеновского всплеска на фазе роста для интервала времени (t1, tm), в течение которого поток рентгеновского излучения изменяется от значения 0.1 Jm до значения Jm, где Jm - поток рентгеновского излучения в максимуме всплеска:

Эта величина имеет ряд преимуществ перед временными параметрами, используемыми в других работах (обычно время роста и спада всплеска). Во-первых, она не является локальной, взятой в одной точке временного профиля, а характеризует весь профиль, в том числе его форму, а не только ширину, взятую в одной точке профиля; во-вторых, она описывает временное развитие всплеска на фазе роста. Использование только временного параметра фазы роста всплеска для прогноза СПС предлагается впервые. В этом существенное отличие от прогноза, предлагаемого в других работах, например в [2]. Такой прогноз обладает преимуществом, по сравнению с прогнозом по фазе спада, так как имеет более высокую заблаговременность. Прогноз по фазе спада в достаточно большом количестве случаев лишен смысла из-за нулевой и даже отрицательной заблаговременности. В [2] для увеличения заблаговременности при прогнозе по фазе спада предлагается экстраполировать временной ход на фазе спада по наблюдению фазы спада на начальной стадии. Для целого ряда событий прогноз следует проводить только по фазе роста потока рентгеновского излучения, так как на фазе спада прогноз уже не имеет смысла из-за близкой к нулю заблаговременности.

Исследование зависимости «протонности» всплесков (т.е. вероятности наблюдения СПС после всплеска) от основных параметров проводилось для всплесков, произошедших на гелиодолготах, для которых гелиодолготное ослабление мало и им можно пренебречь, то есть в области гелиодолгот 30-90°W. Были использованы данные по СПС с энергиями протонов больше 30 МэВ, так как потоки протонов с такими энергиями значительно меньше возмущены межпланетными ударными волнами (МУВ), и в большинстве случаев инжекцию протонов можно считать мгновенной. Для увеличения статистической обеспеченности в рассмотрение были включены события, вплоть до значений потоков в максимуме 5.0*10-1 част/см2 с ср. В предлагаемом методе прогноза рассматривается прогноз СПС от всплесков класса М5 и выше. Оценки, сделанные для всплесков от западной половины диска, показывают, что после всплесков класса М5 и выше наблюдается 83% СПС с интенсивностью в максимуме больше 100част/см2сср и 94% СПС с интенсивностью больше 101част/см2сср.

Диаграммы для нахождения решающего правила

Для исследования зависимости вероятности СПС от выбранных параметров нами были впервые построены диаграммы зависимости Tr,1 от lgJx для всех всплесков класса больше М5 (поток рентгеновского излучения в максимуме превышает 5*10-5 Вт/м2) для сектора гелиодолгот 30-90°W отдельно для событий 22 и 23 циклов солнечной активности и для событий 22 и 23 циклов совместно. Эти диаграммы представлены на фигурах 1-3.

Как следует из фигур 1-3, существует зависимость вероятности протонного всплеска от времени Tr,1. Поскольку существует также зависимость вероятности от интенсивности всплеска, то вероятность растет наиболее быстро при росте параметра, возрастающего как с ростом Tr,1, так и с ростом Jx,m. Предлагается использовать в качестве такого параметра величину P=Tr,1/60+0.35Lg(Jx,m*103), где Tr,1 в минутах, a Jx,m в Вт/м2, которую мы будем называть параметром «протонности» всплеска. Пространство (Tr,1;LgJx,m), как следует из фиг.3 ( а также из фигур 1-2 ), можно разбить на 4 области, которые на фигурах разделяются наклонными пунктирными линиями: область D1, ограниченная осями координат и прямой: Tr,1(мин)=-21.2*lg(105*Jx,m(Вт/м2))+27.6; область D2, ограниченная прямыми: Tr,1(мин)=-21.2*lg(105*Jx,m(Вт/м2))+27.6 и Tr,1(мин)=-21.2*lg(105*Jx,m(Вт/м2))+36.4; область D3, ограниченная прямыми:Tr,1(мин)=-21.2*lg(105*Jx,m(Вт/м2))+36.4 и Tr,1(мин)=-21.2*lg(105*Jx,m(Вт/м2)+48; область D4, выше прямой: Tr,1(мин)=21.2*lg(105*Jx,m(Вт/м2)+48. Значения параметра Р являются постоянными, если параметры всплеска удовлетворяют приведенным уравнениям прямых. В каждой из областей D вероятность наблюдения протонного события разная.

Диаграммы, построенные отдельно для событий 22 и 23 циклов, обладают такими же особенностями в распределении двумерных точек событий (Tr,1, LgJx,m).

Данные, приведенные на диаграмме фиг.3, позволяют рассчитать условные вероятности протонного события при условии наблюдения всплеска рентгеновского излучения с заданными параметрами W(Ai/Bk).

В нашем случае события B1, B2, В3, В4 - наблюдение рентгеновского всплеска с параметрами Tr,1 и lgJx,m такими, что точки с координатами Tr.1, lgJx,m принадлежат областям D1, D2, D3, D4, соответственно. Принадлежность параметров события к области D определяется значением параметра «протонности» Р. Если Р<-0.24, то точка соответствующая всплеску, попадает в область D1, если -0.24≤Р<-0.09 - то в область D2, если -0.09≤Р<0.10 - то в область D3, если 0.10≤Р - то в область D4.

События B1, B2, В3, В4 образуют полную группу несовместимых событий:

События А1, А2, А3, А4, А5 - наблюдение после рентгеновского всплеска протонного события с энергиями протонов больше 30 МэВ с потоком в максимуме в интервалах:

lgJp(част/см2сср)<-0.3, -0.3≤lgJp(част/см2сср)<1, 1≤lgJp(част/см2сср)<2, 2≤lgJp(част/см2сср)<3, lgJp(част/см2сср)≥3. События А1, А2, A3, А4, А5 также образуют полную группу несовместимых событий:

Вероятность - наблюдение события Ai при условии события Bk - W(Ai/Bk) в нашем случае означает условную вероятность наблюдения протонного события с потоком в максимуме в одном из пяти интервалов после всплеска со значениями параметров, принадлежащих заданным областям Dk. Эту вероятность можно обозначить также как:

w (ni≤lgJp<ni+1/(Tr,1;lgJx)⊂Dk), где n1=-0.3, n2=1, n3=2, n4=3, n5=∞.

Для краткости изложения введем следующие обозначения: w(i/k)- условная вероятность наблюдения СПС с уровнем потока протонов с энергией больше 30 МэВ балла i при условии k. Значение индекса i=-1 соответствует уровню потока протонов меньше фона (lgJp<-0.3), значение индекса i=-0.3 соответствует значениям потока протонов в пределах -0.3≤lgJp<1, значение индекса i=1 соответствует значениям потока протонов в пределах l≤IgJp<2, значение индекса i=2 соответствует значениям потока протонов в пределах 2≤lgJp<3, значение индекса i=3 соответствует значениям потока протонов, удовлетворяющим условию lgJp≥3. Индекс k обозначает условие принадлежности точки с параметрами всплеска Tr,1 и Jx,m к области D1 (к=1), к области D2 (к=2), к области D3 (к=3), к области D4 (к=4). Например, w(2/3) означает условную вероятность: w(2≤lgJp<3/(Tr,1;lgJx)⊂D3).

Для вычисления вероятностей событий СПС с интенсивностью выше заданного уровня следует проводить суммирование вероятностей w(i/k) по индексу i. Например, вероятность w(>1/k)=w(1/k)+w(2/k)+w(3/k).

При расчете вероятностей были использованы всплески с гелиодолготами 30-90°W, а также всплески за западным лимбом, происшедшие в 22 и 23 циклах СА: Таких всплесков класса ≥М5 произошло 179 (из них 78 событий класса ≥X1).

Из них:

1.Число «протонных» всплесков (с наблюдением последующего СПС с энергий протонов больше 30 МэВ и потоком в максимуме больше чем 0.5 част/см2сср) - 54;

2. Число всплесков с неопределенной «протонностью» (из-за повышенного фона от предшествующих СПС) - 20

3. Число условно «непротонных» всплесков (поток протонов от которых не превышал значение 0.5 част/см2сср) - 105.

При расчете вероятностей наблюдения событий СПС события с неопределенной «протонностью» также учитывались. Это делалось следующим образом:

Пусть - реальная условная вероятность (в случае, если «протоность» всех всплесков известна ), Np,r - число «протонных» всплесков в области (реальное), Np+ - число наблюдаемых «протонных» всплесков, Np- число «непротонных» всплесков, Np? - число всплесков, «протонность» которых неопределена.

- наблюдаемая условная вероятность «протонного» всплеска;

- условная вероятность наблюдения всплеска с неопределенной «протонностью».

Тогда получим уравнение:

Отсюда следует:

Решающее правило прогноза

Решающее правило прогноза представлено таблицами 1-3

Таблица 1.

Значения условных вероятностей, наблюдения СПС с потоком в максимуме выше заданного уровня в зависимости от значений параметров всплеска рентгеновского излучения для сектора гелиодолгот: 0°≤φ≤90°W

k 1 2 3 4
w(>-0.3/k) 0.09 0.29 0.82 0.91
w(>1/k) 0.00 0.11 0.50 0.82
w(>2/k) 0.00 0.00 0.18 0.55
w(>3/k) 0.00 0.00 0.14 0.18

Таблица 2.

Значения условных вероятностей, наблюдения СПС с потоком в максимуме выше заданного уровня в зависимости от значений параметров всплеска рентгеновского излучения для сектора гелиодолгот: -30°E≤φπ0°

k 1 2 3 4
w(>-0.3/k) 0.00 0.12 0.60 0.86
w(>1/k) 0.00 0.00 0.20 0.36
w(>2/k) 0.00 0.00 0.04 0.00
w(>3/k) 0.00 0.00 0.00 0.00

Таблица 3.

Значения условных вероятностей, наблюдения СПС с потоком в максимуме выше заданного уровня в зависимости от значений параметров всплеска рентгеновского излучения для сектора гелиодолгот: -90°Е≤φπ-30°Е

k 1 2 3 4
w(>-0.3/k) 0.00 0.00 0.18 0.54
w(>1/k) 0.00 0.00 0.14 0.17
w(>2/k) 0.00 0.00 0.00 0.00
w(>3/k) 0.00 0.00 0.00 0.00

Отметим здесь, что вероятности для сектора -30°Е≤φπ0° соответствуют вероятностям для сектора 0°≤φ≤90°W при увеличении интенсивности прогнозируемого потока протонов в 25 раз, а вероятности для сектора - 90°Е≤φπ-30°Е соответствуют вероятностям для сектора 0°≤φ≤90°W при увеличении интенсивности прогнозируемого потока протонов в 100 раз, что соответствует гелиодолготному ослаблению потоков протонов в СПС в 25 и 100 раз, соответственно (см.также М.Н.Беловский, Ю.П.Очелков. Солнечные данные, 1976. т.12, с. 73).

Эти результаты получены нами при исследовании диаграмм lgJp - LgJx,m для секторов гелиодолгот: 0°≤φ≤30°W, 30°≤φ≤90°W, -30°E≤φπ0°, -90°E≤φπ-30°E.

Проверка решающего правила

Нами была проведена проверка решающего правила по данным событий 23 цикла из сектора гелиодолгот 0°-30° W, которые не были использованы для нахождения решающего правила. Всего было рассмотрено 37 событий. Результаты проверки представлены на диаграмме Tr,1-lgJx,m, приведенной на фиг.4.

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

2. Алгоритм мониторинга всплесков рентгеновского излучения класса М5 и выше и определения их параметров.

Алгоритм мониторинга всплесков, потенциально опасных с точки зрения появления после них геоэффективных потоков протонов в ОКП, предполагает несколько режимов работы: режим контроля уровня рентгеновского излучения, режим фазы роста потока рентгеновского излучения, режим фазы спада, режим расчета параметров всплесков, необходимых для прогноза. На фиг.5 изображены режимы работы на разных фазах временной эволюции потока рентгеновского излучения. Общая блок-схема алгортма изображена на фиг.6. Блок схема режима фазы роста представлена на фиг.7. Блок схема режима прогноза представлена на фиг.8.

Режим контроля осуществляется при значении потока рентгеновского излучения меньше чем 4*10-6 Вт/м2. В этом режиме осуществляется сравнение значения потока с этой величиной и переход к режиму фазы роста в случае превышения этого значения. В режиме фазы роста учитывается возможность двух ситуаций: рост до максимального значения потока, не превышающего значения 5*10-5 Вт/м2, и рост до значения потока превышающего значение 5*10-5 Вт/м2. Поскольку заранее неизвестно, какая ситуация реализуется, то в режиме фазы роста всегда происходит запоминание данных по потоку рентгеновского излучения для каждой минуты фазы роста. Если реализуется первая ситуация, то массив данных стирается, если - вторая, то он используется для нахождения параметров прогноза. В первом случае реализуется режим фазы спада. В этом режиме последующее значение потока сравнивается с предыдущим. После достижения минимального значения, в зависимости от его значения реализуется либо режим фазы роста, либо режим контроля.

Если максимальное значение для фазы роста превышает значение 5*10-5 Вт/м2, то определяются параметры всплеска, используемые в решающем правиле прогноза. Вычисляется величина максимального значения потока во всплеске, значение Jmax/10 сравнивается с полученным ранее значением Jmin.

Затем вычисляется средняя ширина временного профиля в логарифмических координатах для потока рентгеновского излучения. В случае если предыдущее максимуму минимальное значение потока, после которого начинается рост до максимума Jmin, удовлетворяет неравенству: Jmin<Jmax/10, для вычислений используется формула:

в которой целочисленный параметр i0 ( номер отсчета), такой, что Jr(i0), значение потока, наиболее близкое к Jmax/10, a Jr(imax)=Jmax (индекс r означает, что значения потока берутся на фазе роста); Tr,1 в минутах.

В случае если предыдущее максимуму минимальное значение потока, после которого начинается рост до максимума Jmin, удовлетворяет неравенству: Jmin>Jmax/10 (это означает, что всплеск является сложным), для вычислений используется формула:

в которой целочисленный параметр i0, такой, что Jr(i0) для найденного k, значение наиболее близкое к Jmax/10k. При этом параметр Tr,1 вычисляется по формуле:

,

где к=0.1, 0.3, 0.5, 0.7, 0,9 f(0.1)=25, f(0.3)=5.5, f(0.5)=2.7, f(0.7)=1.65, f(0.9)=1

Значения функций f получены нами на основе изучения временной эволюции всплесков на фазе роста и равны отношению средней ширины логарифмического временного профиля всплеска на фазе роста от 10-1 Jm до Jm, к средней ширине логарифмического временного профиля на фазах роста от 10-kJm до Jm.

После вычисления Tr,1 по приведенным формулам вычисляется параметр «протонности» Р.

Если параметр «протонности» лежит в интервале Р<-0.24, то для расчета вероятностей используются данные первого столбца таблиц вероятностей 1, 2, 3. Если значение параметра лежит в интервале - -0.24<Р<-0.09, то - данные второго столбца, если значение параметра лежит в интервале:-0.09<Р< 0.10, то - данные третьего столбца, если Р>0.1,то- данные четвертого столбца.

В зависимости от гелиодолготы всплеска выбираются данные таблиц 1, 2, 3 для заданного уровня потока протонов в максимуме СПС: первая строка для lgJp≥-0.3, вторая строка для lgJp≥1, третья строка для lgJp≥2, четвертая строка для lgJp≥3 (Jp в част/см2сср).

Пример применения способа прогноза.

Рассмотрим работу предлагаемого способа прогноза на примере рентгеновского всплеска, произошедшего 13 декабря 2006 года на гелиодолготе W24. Расчеты параметров всплеска дают следующие значения:

Tr,1=12.762 мин. Jx,m=3.76*10-4 Вт/м2, Р=0.064 (что соответствует сектору D3 и высокому уровню протонной опасности). В соответствии с результатами расчета прогноз вероятности события СПС с интенсивностью в максимуме больше 101част/см2с - 82%, с интенсивностью в максимуме больше 102част/см2с - 50%. Наблюдалось протонное событие с потоком в максимуме 4*102част/см2сср.

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

3. Краткое описание чертежей

Фиг.1

Диаграмма зависимости Tr,1 от lgJx. для всплесков сектора 30-90° W класса ≥М5 для 23 цикла СА. Темные кружки - «протонные» всплески (Jp>0.5 част/см2сср), светлые - условно «непротонные» всплески (Jp<0.5 част/см2сср), большие светлые кружки - события с «неопределенной протонностью» из-за фона СПС от предыдущих вспышек. Фиг.2

Диаграмма зависимости Tr,1 от lgJx. для всплесков сектора 30°-90° W класса ≥М5 для 22 цикла СА. Темные кружки - «протонные» всплески (Jp>0.5 част/см2сср), светлые - условно «непротонные» всплески (Jp<0.5 част/см2сср), большие светлые кружки - события с «неопределенной протонностью» из-за фона СПС от предыдущих вспышек.

Фиг.3

Диаграмма зависимости Tr,1 от lgJx для всплесков из сектора гелиодолгот 30-90° W класса ≥М5 для 22 и 23 циклов СА. Светлые малые кружки - события, в которых поток протонов в максимуме Jp<0.5 част/см2сср, большие светлые кружки - события, в которых уровень потока протонов меньше фона от СПС, малые черные кружки - события, в которых поток протонов в максимуме 0.5≤Jp<10 част/см2сср, средние черные кружки - события, в которых поток протонов в максимуме в пределах 10≤Jp<100 част/см2сср, большие черные кружки - события с потоком протонов в максимуме в пределах 100≤Jp<1000 част/см2сср, самые большие черные кружки - события, в которых поток протонов в максимуме Jp≥1000 част/см2сср.

Фиг.4

Диаграмма зависимости Tr,1 от lgJx для всплесков из сектора гелиодолгот 0°-30° W класса ≥М5 для 23 цикла СА. Светлые малые кружки - события, в которых поток протонов в максимуме Jp<0.5 част/см2сср, большие светлые кружки - события, в которых уровень потока протонов меньше фона от СПС, малые черные кружки - события, в которых поток протонов в максимуме 0.5≤Jp<10 част/см2сср, средние черные кружки - события, в которых поток протонов в максимуме в пределах 10≤Jp<100 част/см2сср, большие черные кружки - события, с потоком протонов в максимуме в пределах 100≤Jp<1000 част/см2сср, самые большие черные кружки - события, в которых поток протонов в максимуме Jp≥1000 част/см2сср.

Фиг. 5

Режимы алгоритма мониторинга всплесков и прогноза появления потоков протонов на разных фазах временной эволюции потока рентгеновского излучения. Нижняя кривая - временной ход потока рентгеновского излучения для диапазона длин волн 0.05-0.3 нм, верхняя кривая - для диапазона 0.1-0.8 нм. Сплошные горизонтальные линии - уровни потока рентгеновского излучения для диапазона длин волн 0.1-0.8 нм 4*10-6 Вт/м2 и 5*10-5 Вт/м2, соответствующие режимам фазы роста и прогноза.

Фиг.6

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

Фиг.7

Блок-схема алгоритма наукастинга для режима фазы роста.

Фиг.8

Блок-схема алгоритма наукастинга для режима прогноза.

Способ прогнозирования появления геоэффективных потоков протонов в околоземном космическом пространстве, заключающийся в непрерывном сборе и анализе информации о потоке рентгеновского излучения Солнца в диапазоне длин волн 0,1-0,8 нм, отличающийся тем, что в процессе формирования прогноза в режиме реального времени осуществляют поиск всплесков рентгеновского излучения с интенсивностью превосходящей значение 5·10-5 Вт/м2 и для каждого такого всплеска определяют параметр «протонности» всплеска по соотношению:
P=Tr,1+0,35Lg(Jx,m·103),
где Jx,m - поток рентгеновского излучения в максимуме всплеска, Вт/м2;
Tr,1 - средняя ширина временного профиля рентгеновского всплеска в логарифмических координатах по интенсивности на фазе роста потока рентгеновского излучения от значения в 0,1·Jx,m до Jx.m, ч;
и в зависимости от его значения и гелиодолготы всплеска определяют вероятность наблюдения в околоземном космическом пространстве солнечного протонного события с энергией протонов больше 30 МэВ с интенсивностью в максимуме события выше заданных уровней: 0,5 част./см2сср, 10 част./см2сср, 100 част./см2сср, 1000 част./см2сср.



 

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

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

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

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

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

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

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

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

Изобретение относится к области метеорологии, а более конкретно к предсказанию состояния погоды в микроклиматической зоне. .

Изобретение относится к области метеорологии. .

Изобретение относится к области метеорологического приборостроения и может использоваться в составе средств обеспечения посадки боевых многоцелевых самолетов в дневное время.

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

Изобретение относится к области метеорологии и может быть использовано для определения зон возможного обледенения воздушных судов в режиме реального времени. Согласно заявленному способу проводится регистрация фактических значений вертикального профиля температуры приземного слоя атмосферы n раз при помощи наземного температурного профилемера, который устанавливают в заданном районе наблюдения, а по данным наземных наблюдений определяют приземное значение относительной влажности воздуха, приземное значение температуры точки росы и значение высоты нижней кромки облачности. Затем осуществляют математическую обработку метеорологических данных, используя для расчетов формулу Годске или метод, который предложен в NCEP. Способ может быть использован в первую очередь на аэродромах, где отсутствует регистрация фактических значений вертикальных профилей температуры и влажности воздуха. Технический результат - повышение достоверности определения обледенения воздушных судов. 3 ил.

Изобретение относится к области метеорологии и может быть использовано для прогнозирования погоды. Сущность: выбирают из множества элементов информации о погоде, которые относятся к областям и моментам времени и которые включают в себя, по меньшей мере, температурные данные, множество наборов информации о погоде, относящихся к множеству моментов времени в течение фиксированного периода, касающихся первой области, содержащей местоположение, в котором размещается устройство использования воздуха. Решают с помощью выбранного множества наборов информации о погоде в качестве входных данных дифференциальные уравнения, выражающие информацию о погоде на основе моделей анализа, используемых для проведения моделирований погоды. Формируют множество первых наборов информации о погоде в локальной области, относящихся к множеству вторых областей, которые расположены в пределах первой области и которые меньше, чем первая область. Выбирают второй набор информации о погоде в локальной области, касающийся второй области, содержащей местоположение устройства использования воздуха, среди сформированного множества первых наборов информации о погоде в локальной области. Формируют распределение накопленной повторяемости температуры или распределение вероятности превышения температуры в течение фиксированного периода с помощью данных о температуре, содержащихся во втором наборе информации о погоде в локальной области, для того, чтобы вычислять расчетную температуру устройства использования воздуха. Технический результат: прогнозирование погоды. 3 н. и 15 з.п. ф-лы, 14 ил.

Изобретение относится к области метеорологии и может быть использовано для прогнозирования погоды. Сущность: выбирают информацию о погоде, включающую в себя данные температуры, связанную с временами и областями. Выбирают информацию о погоде, связанную с областью, содержащей местоположение, в котором размещается устройство использования воздуха, связанную с множеством моментов времени в течение некоторого периода. Посредством решения с помощью информации о погоде в качестве входных данных дифференциальных уравнений, выражающих информацию о погоде, на основе моделей метеорологического анализа, используемых для проведения моделирований погоды, формируют информацию о погоде в первой малой области, связанной с областями, меньшими области, соответствующей информации о погоде. Выбирают информацию о погоде во второй малой области относительно области, содержащей местоположение устройства использования воздуха, из информации о погоде в первой малой области. Формируют накопленное температурное распределение или распределение вероятностей превышения температуры за некоторый период посредством использования данных температуры, содержащихся в информации о погоде во второй малой области, для вычисления расчетной температуры устройства использования воздуха. Технический результат: прогнозирование погоды. 3 н. и 12 з.п. ф-лы, 14 ил.

Способ исследования изменений климата Земли заключается в том, что измерительную систему, включающую два идентичных оптических телескопа, располагают на видимой поверхности Луны. Во время движения Луны вокруг Земли оптические телескопы последовательно производят измерения энергии отраженного и рассеянного Землей солнечного излучения по всем направлениям и суммарной энергии собственного теплового излучения поверхности и атмосферы Земли. Затем определяют значение альбедо Бонда Земли, величину отклонения энергетического баланса Земли от равновесного состояния с использованием заранее известных данных солнечной постоянной. И на основании полученных данных оценивают изменения энергетического состояния Земли и степени пропускания атмосферой теплового излучения земной поверхности в космическое пространство, а также состояние климата. Технический результат - повышение точности и надежности исследования климатических изменений, происходящих на Земле. 2 н. и 5 з.п. ф-лы, 3 ил.
Наверх