ГОСТ 25645.126-85


Группа Т27

ГОСУДАРСТВЕННЫЙ СТАНДАРТ СОЮЗА ССР



ПОЛЕ ГЕОМАГНИТНОЕ


Модель поля внутриземных источников


Geomagnetic field.

Magnetic field model of internal originals

ОКСТУ 0080

Дата введения 1987-01-01

ИНФОРМАЦИОННЫЕ ДАННЫЕ


1. УТВЕРЖДЕН И ВВЕДЕН В ДЕЙСТВИЕ Постановлением Государственного комитета СССР по стандартам от 14.11.85 N 3609


ИСПОЛНИТЕЛИ


И.И.Алексеев, канд. физ.-мат. наук; А.В.Баюков, канд. техн. наук; Е.С.Беленькая, канд. физ.-мат. наук; Н.П.Бенькова, д-р физ.-мат. наук; Ю.А.Винтенко, канд. техн. наук; В.П.Головков, д-р физ.-мат. наук; Е.В.Горчаков, д-р физ.-мат. наук; М.С.Григорян; И.П.Иваненко, д-р физ.-мат. наук; В.В.Калегаев; Г.И.Коломийцева, канд. физ.-мат. наук; А.П.Кропоткин, д-р физ.-мат. наук; Е.Н.Лесновский, канд. техн. наук; В.М.Ломакин, канд. техн. наук; Ю.Г.Лютов; В.В.Мигулин, член-кор. АН СССР; Л.И.Мирошниченко, канд. физ.-мат. наук; В.Н.Никитинский; И.Я.Ремизов, канд. техн. наук; В.И.Степакин, канд. техн. наук; Л.Н.Степанова; И.Б.Теплов, д-р физ.-мат. наук; М.В.Терновская, канд. физ.-мат. наук; В.В.Хаустов, канд. техн. наук


2. СОГЛАСОВАНО с Государственной службой стандартных справочных данных (протокол от 16.06.85 N 18)


3. Срок первой проверки - 1989 г., периодичность проверки - 5 лет


4. ВВЕДЕН ВПЕРВЫЕ


5. ССЫЛОЧНЫЕ НОРМАТИВНО-ТЕХНИЧЕСКИЕ ДОКУМЕНТЫ


Обозначение НТД, на который дана ссылка


Номер пункта


ГОСТ 25645.127-85


1.1




6. ПЕРЕИЗДАНИЕ (декабрь 1989 г.) с Изменением N 1, утвержденным в сентябре 1989 г. (ИУС 12-89)


7. Проверен в 1989 г.



Настоящий стандарт устанавливает модель геомагнитного поля внутриземных источников на расстоянии от 100 до 40000 км от поверхности Земли.


Стандарт предназначен для использования в расчетах при определении условий функционирования технических устройств в космическом пространстве.




1. ОСНОВНЫЕ ПОЛОЖЕНИЯ


1.1. Вектор индукции магнитного поля в магнитосфере Земли вычисляют по формуле


, нТл, (1)

где - вектор индукции геомагнитного поля внутриземных источников;


- вектор индукции магнитного поля магнитосферных токов по ГОСТ 25645.127-85.


1.2. Магнитное поле внутриземных источников состоит из поля электрических токов в земном ядре (далее - главное поле), составляющего ~98% всего поля, и поля земной коры, являющегося полем магнетизма горных пород и составляющего ~2% всего поля.


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


1.3. Модель главного поля представлена рядами сферических гармоник в зависимости от географических координат. При длине ряда 10-13 гармоник погрешность вычисления геомагнитного поля на поверхности Земли составляет 2%.


В первом приближении геомагнитное поле является полем диполя, расположенного в центре Земли, и представляется первым членом сферического гармонического ряда.


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




2. МОДЕЛЬ ГЛАВНОГО ГЕОМАГНИТНОГО ПОЛЯ ВНУТРИЗЕМНЫХ ИСТОЧНИКОВ


2.1. Потенциал индукции геомагнитного поля внутриземных источников в точке пространства со сферическими координатами , , вычисляют по формуле


, нТл·км, (2)


где полюс сферической системы координат совпадает с географическим полюсом Земли;


- геоцентрическое расстояние, км;


- долгота от Гринвичского меридиана, ...°;


- дополнение до широты, , ...°;


- широта в сферических координатах, ...°;


- средний радиус Земли, км;


(3)

где - нормировочный множитель,


2 для 1 и 1;


, - сферические гармонические коэффициенты, нТл;


- степень сферических гармоник;


- порядок сферических гармоник;


10 - максимальная степень сферических гармоник.


2.2. Все экспериментальные данные и положения ИСЗ в пространстве представляют в географических (геодезических) координатах , , , основанных на аппроксимации поверхности Земли эллипсоидом вращения. В ряде задач в первом приближении эллиптичностью Земли пренебрегают, не делая разницы между сферическими и геодезическими координатами. Однако при более точных расчетах необходимо учитывать сжатие Земли. Для учета сжатия Земли и вычисляют по формулам:


, (4)

, (5)

где - географическая (геодезическая) широта точки в пространстве, ...°;


- высота точки над уровнем моря, км;


- большая полуось земного эллипсоида вращения, км;


- малая полуось земного эллипсоида вращения, км.


Долготы в сферических и геодезических координатах тождественны.


Примечание. Значения и приведены в рекомендуемом приложении 1.



2.3. Составляющие вектора индукции геомагнитного поля внутриземных источников : , и вычисляют по формулам:


, нТл, (6)

, нТл, (7)

, нТл. (8)


Составляющие , и используют для расчета вектора индукции по формуле (1).


(Измененная редакция, Изм. N 1).


2.4. Пространственно-временное распределение вектора индукции геомагнитного поля обычно описывают геомагнитными элементами:


прямоугольными составляющими , , , , нТл;


угловыми элементами и , ...°;


модулем вектора индукции , нТл.


Определения геомагнитных элементов приведены в приложении 2.


2.4.1. В точке пространства с координатами , , прямоугольные составляющие вектора индукции в геодезической системе координат рассчитывают по формулам:


; (9)

, (10)

; (11)

. (12)


(Измененная редакция, Изм. N 1).


2.4.2. Угловые элементы и модуль вектора индукции вычисляют по формулам:


; (13)

; (14)

. (15)


2.4.3. Значения элемента поля для точки пространства при 0 получают линейной интерполяцией.


2.4.4. Значения сферических гармонических коэффициентов , для 1985 г. приведены в приложении 1, а результаты расчета поля на тот же год - в приложении 3.


Расчет поля на другие годы осуществляют с помощью векового хода. Пример программы для расчета геомагнитных элементов приведен в приложении 4.


(Измененная редакция, Изм. N 1).


2.5. Вековой ход геомагнитного поля определяют потенциалом , который вычисляют по формуле


, нТл км/год, (16)

где , - сферические гармонические коэффициенты, нТл/год.


2.5.1. При расчетах векового хода не учитывают сжатие Земли и пренебрегают различием между сферическими и задаваемыми географическими координатами (полагают , ).


2.5.2. Вековой ход элементов геомагнитного поля рассчитывают по формулам:


, нТл/год, (17)

, нТл/год, (18)

, нТл/год, (19)

, нТл/год, (20)

, ...’/год, (21)

, ...’/год, (22)

, нТл/год, (23)

где элементы поля , , , , , , вычисляют по формулам (9-15);


при 0 определяют линейной интерполяцией.


2.5.3. Сферические гармонические коэффициенты , определяют по экспериментальным данным для различных временных интервалов. Значения , для 1985-1990 гг. приведены в приложении 5. Пример расчета векового хода приведен в приложении 3.


(Измененная редакция, Изм. N 1).


2.6. Главное поле на любой заданный год вычисляют способами, приведенными в пп.2.6.1 и 2.6.2.


2.6.1. Геомагнитные элементы на заданный год вычисляют по формуле


, (24)

где - любой из элементов поля ( , , , , , и ) на год ;


- элемент поля, рассчитанный по , (известным на год ) по формулам (6-15) с учетом формул (4-5);


- вековой ход элемента поля, рассчитанный по , по формулам (17-23).


2.6.2. вычисляют по формулам (6-15), в которых , заменяют на


; (25)

*. (26)

________________

* Формула соответствует оригиналу. - Примечание изготовителя базы данных.


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


2.6.4. Пример расчета на 1989 г. по п.2.6.2 дан в программе, приведенной в приложении 4.


2.6.5. Примеры расчета поля на 1988 г. обоими способами приведены в приложении 3.


2.6.4, 2.6.5. (Измененная редакция, Изм. N 1).




3. ПАРАМЕТРЫ ГЕОМАГНИТНОГО ДИПОЛЯ


3.1. Дипольное геомагнитное поле соответствует полю, представленному первым членом сферических гармоник. Составляющие дипольного члена рассчитывают по формулам:


;


;


.


3.2. Координаты полюсов дипольного поля (геомагнитных полюсов) и его магнитный момент рассчитывают по формулам:


, (27)

, (28)

, Тл·м , (29)

где - географическая широта геомагнитного полюса, ...°;


- географическая долгота геомагнитного полюса, ...°.


3.3. Параметры геомагнитного диполя для 1985 г. приведены в приложении 6. Пример расчета дипольного поля приведен в приложении 3.


(Измененная редакция, Изм. N 1).




ПРИЛОЖЕНИЕ 1

Рекомендуемое


Сферические гармонические коэффициенты , , нТл, для 1985 г.



1


0


-29877


0


1


1


-1903


5497


2


0


-2073


0


2


1


3045


-2191


2


2


1691


-309


3


0


1300


0


3


1


-2208


-312


3


2


1244


284


3


3


835


-296


4


0


937


0


4


1


780


233


4


2


363


-250


4


3


-426


68


4


4


169


-298


5


0


-215


0


5


1


356


47


5


2


253


148


5


3


-94


-155


5


4


-161


-75


5


5


-48


95


6


0


52


0


6


1


65


-16


6


2


50


90


6


3


-186


69


6


4


4


-50


6


5


17


-4


6


6


-102


20


7


0


75


0


7


1


-61


-82


7


2


2


-26


7


3


24


-1


7


4


-6


23


7


5


4


17


7


6


9


-21


7


7


0


-6

8

0

21


0


8

1

6

7


8

2


0


-21


8


3


-11


5

8

4

-9

-25


8


5

2


11


8


6


4


12


8


7


4


-16


8


8


-6


-10

9


0


5


0


9


1


10


-21


9


2


1


16


9


3


-12


9

9


4


9


-5


9


5


-3


-6


9


6


-1


9

9


7


7


10

9


8


2


-6

9


9


-5


2


10


0


-4


0

10


1


-4

1


10


2


2


0


10


3

-5


3

10


4


-2

6


10


5


5

-4


10


6


3

0


10


7


1


-1


10


8


2


4


10


9


3

0


10


10


0


-6