Астрономический вестник, 2020, T. 54, № 6, стр. 520-528

Электрофизическая модель грунта Луны

О. В. Юшкова a*, И. Н. Кибардина a, Т. Н. Дымова a

a Фрязинский филиал Институт радиотехники и электроники им. В.А. Котельникова РАН
Фрязино, Россия

* E-mail: o.v.y@mail.ru

Поступила в редакцию 28.02.2020
После доработки 11.07.2020
Принята к публикации 13.07.2020

Полный текст (PDF)

Аннотация

На основе анализа имеющихся в литературе данных о лабораторных измерениях диэлектрических характеристик образцов лунного грунта, доставленных на Землю в ходе миссий “Луна” и “Apollo”, разработана электрофизическая модель верхнего слоя грунта Луны. Модель учитывает температурные условия в районе проведения измерений, изменение температуры и распределение плотности грунта в верхнем слое реголита. Модель предназначена для численного моделирования экспериментов, направленных на исследование свойств лунного грунта с помощью электромагнитных волн.

Ключевые слова: Луна, реголит, коренные породы, диэлектрическая проницаемость, электромагнитные волны, электрофизическая модель

ВВЕДЕНИЕ

Исследования Луны с помощью космических аппаратов начались 2 января 1959 с запуска автоматической межпланетной станции Луна-1. С тех пор спутник Земли исследовали космические аппараты России, США, Японии, Индии, Европы, Китая. Российская лунная программа предусматривает запуск к Луне четырех автоматических станций. Среди экспериментов, направленных на изучение строения и состава верхнего слоя грунта запланировано выполнить георадарные измерения с лунохода, моно и бистатическую радиолокацию поверхности, каротажное зондирование с посадочного модуля. При создании научной аппаратуры для проведения таких измерений, для прогноза их результатов и отработки технологии решения обратных задач необходимо уметь численно моделировать натурные эксперименты. Ключевым элементом при моделировании радиофизических измерений является электрофизическая модель лунного грунта. Как было показано в ряде работ (Heiken и др., 1991), диэлектрическая и магнитная проницаемости грунта Луны зависят от состава грунта, температуры и частоты электромагнитного поля. Основной магнитный минерал лунных пород – чистое железо, общее содержание примесей которого примерно одинаково и составляет 0.5–0.6% веса для всех мест посадок. Кроме ферромагнетиков в состав лунного грунта входят вещества, относящиеся к парамагнетикам и антиферромагнетам (Флоренский и др., 1975). Антиферромагнетики становятся парамагнетиками при очень низких температурах. Парамагнетики приобретают магнитные свойства под действием внешнего магнитного поля и относятся к слабомагнитным веществам, их магнитная проницаемость незначительно отличается от единицы, поэтому ниже при построении электрофизической модели мы ограничимся рассмотрением только диэлектрической проницаемости лунного грунта. Модель основана на теории диэлектриков, учитывает результаты лабораторных измерений электрофизических характеристик образцов коренной породы и реголита, доставленных на Землю с помощью миссий Луна и Apollo, и может быть использована при моделировании в диапазоне частот 100 МГц–10 ГГц при температуре поверхности от 80 до 380 К.

Грунт Луны представляет собой сложную систему, образованную из скальных и осадочных пород с геологически сложившейся структурой. Дистанционные и контактные наблюдения показали, что грунт Луны имеет выраженную слоистую структуру (Черкасов и др., 1975): коренные породы, в основном базальтовые в морских районах и анортозиты в континентальных районах, покрыты слоем реголита. Реголит представляет собой рыхлый однородный грунт, сформированный из гомогенизированных обломков коренных пород, разрушенных под воздействием температурных колебаний и метеоритных бомбардировок. В справочнике (Heiken и др., 1991) приведено описание образцов коренных пород и реголита, доставленных на Землю. Лабораторные исследования образцов (Heiken и др., 1991) подтвердили гипотезу, высказанную В.Д. Кротиковым и В.С. Троицким (Кротиков, Троицкий, 1963) о том, что реголит по электрофизическим свойствам похож на диэлектрики. Электрические характеристики диэлектрика зависят от минералогического состава, температуры и частоты электромагнитного поля. Поэтому далее мы рассмотрим существующую информацию о динамике суточной температуры на поверхности и внутри верхнего слоя грунта, распределении плотности грунта по глубине и результаты лабораторных исследований диэлектрических характеристик образцов. Все эти факторы будут включены в модель.

ТЕМПЕРАТУРА НА ПОВЕРХНОСТИ ЛУНЫ

На Луне нет атмосферы, поэтому перепады температуры на ее поверхности велики. Температура поверхности измерялась непосредственно во время миссий Apollo-15, -17. С 2007 г. начался этап дистанционного картографирования температурных условий на поверхности Луны. На спутнике Chang’e-1 был установлен пассивный микроволновый радиометр (MRM) с четырьмя частотными каналами 3, 7.8, 19.35 и 37 ГГц (Zheng и др., 2012). Измерения, выполненные в течение восьми полных лунных дней, были использованы в качестве основы для создания глобальной карты яркостной температуры Луны. В 2009 г. на борту Lunar Reconnaissance Orbiter (LRO) был размещен девяти канальный радиометр Diviner Lunar Radiometer Experiment (DLRE) (Paige и др., 2010). Семь тепловых каналов работали в диапазоне длин волн от 7.6 до 400 мкм для изучения яркостной температуры поверхности Луны. Узкий угол диаграммы направленности радиометра и низкая орбита аппарата (50 км) обеспечили высокое пространственное разрешение измерений (приблизительно 200 м) и позволили построить детальную карту суточного изменения радиояркостной температуры лунной поверхности (Williams и др., 2017). Приповерхностный слой грунта, доступный для изучения MRM, составляет несколько сантиметров, а DLRE принимает излучение практически с поверхности. С учетом натурных измерений была построена теоретическая модель для описания суточного изменения температуры слоя грунта в зависимости от координат и зенитного угла Солнца (Vasavada и др., 1999) и (Huang, 2008). На рис. 1 приведен суточный ход температуры поверхности (Vasavada и др., 1999). Анализ результатов работ, посвященных температуре грунта Луны, позволил сделать следующие выводы:

Рис. 1.

Суточный ход температуры на поверхности в месте посадки Apollo-15. Сплошная линия – модель, точки – измерения DLRE (Vasavada и др., 1999).

1. Значения физической температуры, пересчитанные из радиояркостной температуры, совпадают с теоретическими оценками и результатами контактных измерений, проведенными во время посадки Apollo-15. Точками на рис. 1 обозначены суточные измерения температуры на посадочной площадке Apollo-15, сплошной линией приведены результаты ее моделирования с учетом свойств реголита, теплового потока и контактных измерений физической температуры.

2. Наибольшая разница между дневными и ночными температурами наблюдается на экваторе и достигает 300 К. Самая низкая температура, зарегистрированная прибором DLRE, составляет 25–35 К ночью в зоне ударных кратеров Южного полюса.

3. Широтные изменения ночной температуры – незначительны.

Согласно работе (Кротиков, Троицкий, 1963), реголит обладает высокой теплоизоляцией. В приближении моделирование теплового режима в грунте сводится к решению нелинейного одномерного уравнения теплопроводности в конечной области с учетом граничных условий. На основе теории теплового баланса был разработан алгоритм расчета глубинного профиля температуры в лунном грунте (Маров и др., 2010). В этой работе конечная область оценена как 10–15 см, а при наших расчетах мы полагаем, что суточные колебания температуры затрагивают слой около 1 м (согласно статьям (Шевченко, 1990; Минчин, Улубеков, 1972)). Этот параметр модели, так же как теплопроводность и теплоемкость реголита и основных пород, может быть изменен в зависимости от накопленной информации о свойствах грунта и положения конкретного места проведения измерений. На рис. 2 приведены результаты расчета профиля температуры в поверхностном слое реголита Луны для различных поверхностных температур. Горизонтальная ось графика указывает глубину в метрах, а вертикальная ось – температуру в градусах Кельвина. Линия 1 соответствует температуре поверхности 380 К, линия 2 – 300 К, 3 – 230 К, 4 – 150 К, 5 – 80 К. Температура грунта ниже глубины 1 м постоянна и равна 230 К.

Рис. 2.

Глубинный профиль температуры реголита Луны. Температура на поверхности 1 – 380 К, 2 – 300 К, 3 – 230 К, 4 – 150 К, 5 – 80 К.

ГРУНТ ЛУНЫ

Плотность грунта. Информация о плотности коренных пород и реголита стала доступна после проведения лабораторных анализов доставленных на Землю образцов лунного грунта. Подробный обзор физико-механических характеристик грунта Луны приведен в (Слюта, 2014). Из мест посадок космических аппаратов были доставлены буровые керны грунта с нетронутой структурой. Исследование кернов позволило вывести формулу для расчета плотности реголита как функцию глубины z (Olhoeft, Strangway, 1975):

(1)
$\rho \left( z \right) = {{\rho }_{0}} + n~{\text{ln}}\left( {z + 1} \right).$

В этой формуле $\rho \left( z \right)$ измеряется в г/см3, z – в сантиметрах, значение $n$ определено эмпирически и равно 0.121. Значение ${{\rho }_{0}}$ соответствует плотности реголита на поверхности. Формула (1) верна для расчета плотности реголита в слое толщиной до 100 м.

Толщина реголита. Толщина слоя реголита зависит от его возраста и расположения, она может быть до нескольких десятков метров, при этом слой мощнее в гористой местности (Heiken и др., 1991). Средняя толщина слоя реголита составляет около 4–5 м для морей и 10–15 м для горной местности. Согласно измерениям георадара LPR (Lunar Penetrating Radar), который был установлен на луноходе Yutu, толщина реголита в месте посадки КА Chang’E-3 составляла 3–5 м (Xiao и др., 2015). Для определенности при расчетах мы полагали, что толщина реголита равна 4 м, хотя не ограничивая общности, толщину реголита в модели можно менять в зависимости от расположения района исследования и по мере накопления о нем более точной информации.

Диэлектрические характеристики образцов лунного грунта. В этом разделе предложено математическое описание диэлектрических характеристик образцов лунного реголита в зависимости от температуры реголита и частоты радиоволн. В работе (Кибардина, Юшкова, 2017) проведен анализ имеющихся в литературе данных о диэлектрической проницаемости лунных пород, в частности о $\varepsilon {\kern 1pt} '$ (действительная часть диэлектрической проницаемости) и тангенсе угла потерь . Не смотря на то, что на Землю было доставлено более двух тысяч образцов лунного грунта, в работе (Heiken и др., 1991) приведены только семь, у которых $\varepsilon {\kern 1pt} '$ и ${\text{tg}}\delta $ были измерены в нескольких частотных диапазонах: 0.1–10, 450 МГц, 9.375 ГГц, при температуре, примерно равной 290 К, и похожих условиях хранения. Реголит представлен четырьмя образцами: 10 084 (fines, Apollo-11), 12 070 (fines, Apollo-12), 14163 (fines, Apollo-14), 66 041 (fines, Apollo-16), коренные породы – образцами 14 310 (crystal rock, Apollo-14), 60 025 (anorthosites, Apollo-16) и 72 441 (basalt, Apollo-17). Плотность и процентное содержание двуокиси титана и оксида железа этих образцов даны в табл. 1.

Таблица 1.  

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

  N %${\text{Ti}}{{{\text{O}}}_{{\text{2}}}}$ [5 ] %${\text{FeO}}$ [5 ] $\rho $, г/см3
Реголит 10 084 7.56 15.94 1.94
12 070 2.81 16.40 1.74
14 163 1.77 10.45 1.2
66 041 0.63 5.80 1.5
Коренные породы 14 310 1.30 7.70 7
60 025 0.02 0.67 2.1
72 441 1.53 8.68 1.65

Измерения диэлектрических характеристик проводились несколькими научными группами, поэтому один образец грунта делился на фрагменты, которые имели похожую, но не одинаковую плотность. Нам надо было построить график тангенса диэлектрических потерь образцов для частот диапазона от 10 Гц до 10 ГГц. Так как в работе (Кротиков, Троицкий, 1963) высказано предположение, согласно которому для сухих грунтов величина ${{{\text{tg}}\delta } \mathord{\left/ {\vphantom {{{\text{tg}}\delta } \rho }} \right. \kern-0em} \rho }$ является инвариантной относительно плотности, а в работе (Ржевский, 1976) эта гипотеза была экспериментально подтверждена для мелкозернистых лунных пород для описания потерь в грунте, нами была выбрана именно эта, не совсем стандартная для радиофизики, зависимость. На рис. 3 приведены результаты измерений ${{{\text{tg}}\delta } \mathord{\left/ {\vphantom {{{\text{tg}}\delta } \rho }} \right. \kern-0em} \rho }$ образцов реголита: треугольные маркеры соответствуют образцу 10 084, круглыми маркерами даны результаты измерений для образца 12 070, звездочки соответствуют образцу 14 163; квадратные маркеры – образцу 66 041. На рис. 4 представлены графики частотной зависимости ${{{\text{tg}}\delta } \mathord{\left/ {\vphantom {{{\text{tg}}\delta } \rho }} \right. \kern-0em} \rho }$ для коренных грунтов. Квадратными маркерами отмечены экспериментальные измерения, соответствующие кристаллическому образцу 14 310; треугольные – анортозитовому образцу 60 025, круглые – образцу 72 441. На этом рисунке по горизонтальной оси отложены значения $\lg f$, где частота $f$ дана в Гц, по вертикальной оси – . Для рассмотренных образцов диапазон изменения значений ${{{\text{tg}}\delta } \mathord{\left/ {\vphantom {{{\text{tg}}\delta } \rho }} \right. \kern-0em} \rho }$ осадочных пород значительно шире диапазона изменения ${{{\text{tg}}\delta } \mathord{\left/ {\vphantom {{{\text{tg}}\delta } \rho }} \right. \kern-0em} \rho }$ коренных пород. Качественное поведение графиков для всех образцов одинаково. В диапазоне частот ${{10}^{2}}{\kern 1pt} - {\kern 1pt} {{10}^{6}}$ Гц функция ${{{\text{tg}}\delta } \mathord{\left/ {\vphantom {{{\text{tg}}\delta } \rho }} \right. \kern-0em} \rho }$ убывает, что, связано с электропроводимостью грунта. На частотах 6–8 МГц значения функции минимальны; последующий ее рост на частотах выше 8 МГц, по-видимому, связан с релаксационными эффектами (Тареев, 1982). В табл. 2 дано соответствие экспериментальных значений на графике и ссылок на публикацию.

Рис. 3.

Зависимость нормированного тангенса потерь реголита от частоты. Маркерами отмечены экспериментальные измерения: треугольные соответствуют образцу 10 084, круглые – 12 070, звездочки – 14 163, квадратные – 66041.

Рис. 4.

Зависимость нормированного тангенса потерь грунта коренных пород от частоты. Маркерами отмечены экспериментальные измерения: квадратные соответствуют образцу – 14 310, треугольные – 60 025, круглые – 72 441.

Образцы 14163 и 14310 доставлены из одного района Луны, как и образцы 66041 и 60025. Во второй паре образец 60025 является анортозитом, характерной породой для места приземления космического аппарата Apollo-16, а образец 14310 – кристаллической породой, по плотности и иным свойствам он отличался от других пород, доставленных Apollo-14 из того же места посадки. Поэтому именно пара 66041 и 60025 была выбрана нами для моделирования Луны в виде: геологически-однородного слоя реголита, лежащего на однородной коренной породе. На рис. 5 приведены результаты лабораторных измерений в частотном диапазоне 100 Гц–9.375 ГГц при температуре 290 К $\varepsilon {\kern 1pt} '$ и ${{{\text{tg}}\delta } \mathord{\left/ {\vphantom {{{\text{tg}}\delta } \rho }} \right. \kern-0em} \rho }$ только для образца 66041. На приведенном графике треугольные маркеры соответствуют экспериментальным измерениям $\varepsilon {\kern 1pt} '$ из (Katsube, Collett, 1973), круглые – ${{{\text{tg}}\delta } \mathord{\left/ {\vphantom {{{\text{tg}}\delta } \rho }} \right. \kern-0em} \rho }$ из (Katsube, Collett, 1973), квадратные – $\varepsilon {\kern 1pt} '$ и ${{{\text{tg}}\delta } \mathord{\left/ {\vphantom {{{\text{tg}}\delta } \rho }} \right. \kern-0em} \rho }$ из (Gold и др., 1973). Далее приводится анализ графика с точки зрения теории диэлектриков.

Рис. 5.

Частотные зависимости $\varepsilon {\kern 1pt} '$ и ${{{\text{tg}}\delta } \mathord{\left/ {\vphantom {{{\text{tg}}\delta } \rho }} \right. \kern-0em} \rho }$ образца лунного грунта 66041 (Apollo-16). Треугольные маркеры – экспериментальные измерения $\varepsilon {\kern 1pt} '$ из (Katsube, Collett, 1973), круглые – ${{{\text{tg}}\delta } \mathord{\left/ {\vphantom {{{\text{tg}}\delta } \rho }} \right. \kern-0em} \rho }$ из (Katsube, Collett, 1973), квадратные – $\varepsilon {\kern 1pt} '$ и ${{{\text{tg}}\delta } \mathord{\left/ {\vphantom {{{\text{tg}}\delta } \rho }} \right. \kern-0em} \rho }$ из (Gold и др., 1973). Черные линии – результаты расчета по формуле (2), серая – по формуле (5).

В теории диэлектриков для описания комплексной диэлектрической проницаемости диэлектриков с ионной структурой, аналогов реголита на Земле, обладающих и проводимостью и поляризацией (Сканави, 1949), используется соотношение

(2)
$\varepsilon \left( \omega \right) = {{\varepsilon }_{\infty }} + \frac{{\left( {{{\varepsilon }_{0}} - {{\varepsilon }_{\infty }}} \right)}}{{1 + {{\omega }^{2}}{{\tau }^{2}}}} - i\left[ {\frac{{{\omega \tau }\left( {{{\varepsilon }_{0}} - {{\varepsilon }_{\infty }}} \right)}}{{1 + {{\omega }^{2}}{{\tau }^{2}}}} + \frac{\gamma }{\omega }} \right].$

Здесь ${{\varepsilon }_{0}}$ – значение диэлектрической проницаемости вещества при предельно низких частотах, ${{\varepsilon }_{\infty }}$ – предельное значение диэлектрической проницаемости вещества при высоких частотах, $\omega = 2\pi f$, $f$ – частота приложенного электромагнитного поля в Герцах, $\tau $ – макроскопическое время релаксации в с, ${{i}^{2}} = - 1$, $\gamma $ – проводимость в 1/c. При ионной проводимости число ионов находится в экспоненциальной зависимости от температуры и задается формулой (Фрелих, 1960):

(3)
$\gamma = {{\gamma }_{0}}{\text{exp}}\left( {{{ - B} \mathord{\left/ {\vphantom {{ - B} T}} \right. \kern-0em} T}} \right).$

В этой формуле параметры ${{\gamma }_{0}}$ и $B$ определяются минералогическим составом диэлектрика, $T$ – температура в градусах Кельвина. Время релаксации, входящее в соотношение (2), тоже экспоненциально связано с температурой выражением (William, Thomas, 1998):

(4)
$\tau = {{\tau }_{0}}{\text{exp}}\left( {{F \mathord{\left/ {\vphantom {F T}} \right. \kern-0em} T}} \right),$
где значения ${{\tau }_{0}}$ и $F$ также определяются минералогическим составом диэлектрика.

В работе (Юшкова, Кибардина, 2017) подробно описан метод определения параметров, входящих в формулы (2)(4). Для образца 66041 они равны: ${{\varepsilon }_{0}} = 2.7$, ${{\varepsilon }_{\infty }} = 1.4$, ${{\gamma }_{0}} = 1971$ 1/c, $B = {\text{550}}$, ${{\tau }_{0}} = 1.04 \times {{10}^{{ - 16}}}$ c, $~F = {\text{2900}}$, для образца 60 025 – ${{\varepsilon }_{0}} = 5$, ${{\varepsilon }_{\infty }} = 4.6$, и ${{\gamma }_{0}} = {\text{4088}}$ 1/c, $~B = {\text{550}}$, ${{\tau }_{0}} = 6.8 \times {\text{1}}{{{\text{0}}}^{{ - {\text{15}}}}}$ c, $F = {\text{2900}}$.

На рис. 5 черными линиями приведены результаты расчета по формуле (2), с учетом параметров (3), (4) для образца 66041.

В работе (Кибардина, Юшкова, 2017) показано, что для всех семи образцов лунного грунта зависимость ${{{\text{tg}}\delta } \mathord{\left/ {\vphantom {{{\text{tg}}\delta } \rho }} \right. \kern-0em} \rho }$ от частоты может быть аппроксимирована функцией:

(5)
${{{\text{tg}}\delta } \mathord{\left/ {\vphantom {{{\text{tg}}\delta } \rho }} \right. \kern-0em} \rho } = D{\text{1}}{{{\text{0}}}^{{a{{{\left( {{\text{lg}}f - {\text{lg}}{{f}_{0}}} \right)}}^{2}}}}},$
где $f$ – частота в Гц, ${{f}_{0}}$ в Гц – частота, при которой для данного образца принимает наименьшее значение, равное $D = {{10}^{\varphi }}$, $\rho $ – плотность образца. Константы $\varphi $ и $a$ связаны с минералогическими особенностями образца, в частности, с содержанием в грунте оксидов металлов. Для образца 66 041: $a = 0.{\text{074}}$, $\varphi = - 2.9$, ${\text{lg}}{{f}_{0}} = 6.69$; для образца 60 025: $a = 0.{\text{0684}}$, $\varphi = - 3.{\text{271}}$, ${\text{lg}}{{f}_{0}} = 6.7$. Результаты расчета по формуле (5) для tgδ/ρ представлены на рис. 3 и 4 сплошной черной линией, а рис. 5 – серой линией.

Как видно из рис. 5, формула (2) дает грубое приближение в окрестности минимума функции ${{{\text{tg}}\delta } \mathord{\left/ {\vphantom {{{\text{tg}}\delta } \rho }} \right. \kern-0em} \rho }$: значение занижено на три порядка, но учитывает связь диэлектрических параметров грунта с его температурой, в отличие от формулы (5). Ниже в модели были скомбинированы функции (2) и (5): для низких и высоких частот использовалось приближение (2), для средних – (5). Для этого при расчете диэлектрической проницаемости реголита для других температур, в формулу (5) введена температурная зависимость коэффициентов ${{f}_{0}}\left( T \right)$, $a\left( T \right)$, $D\left( T \right)$ так, чтобы частота ${{f}_{0}}$, на которой наблюдается локальный минимум функции (5), менялась с температурой так же как для функции, заданной формулой (2). Для функции (2) такое значение находится из условия и равно ${{f}_{0}} = {{{{\omega }_{{{\text{min}}}}}} \mathord{\left/ {\vphantom {{{{\omega }_{{{\text{min}}}}}} {2\pi }}} \right. \kern-0em} {2\pi }}$, где ${{\omega }_{{{\text{min}}}}} = \sqrt {\frac{{A - 2 - \sqrt {{{A}^{2}} - 8A} }}{{2{{\tau }^{2}}\left( {A + 1} \right)}}} $, $A = \frac{{{{\varepsilon }_{0}} - {{\varepsilon }_{\infty }}}}{{{\gamma \tau }}}$ (Юшкова, Кибардина, 2017). Значения $\tau $ и $\gamma $ зависят от температур согласно выражениям (3) и (4). По аналогии с приближением (2), значение функции $D\left( T \right)$ в (5) принято пропорциональным величине $\sqrt {{\gamma \tau }} $, т.е.

$D\left( T \right) = D{\text{exp}}\left( {{{\left( {F - B} \right)} \mathord{\left/ {\vphantom {{\left( {F - B} \right)} {2T}}} \right. \kern-0em} {2T}}} \right).$

Параметр $a\left( T \right)$ задан из условия равенства функций (2) и (5) значению в точке ${{\omega }_{{{\text{max}}}}}$:

$\begin{gathered} a\left( T \right) = {\text{lg}}{{\left( {\frac{{{{\varepsilon }_{0}} - {{\varepsilon }_{\infty }}}}{{2\varepsilon {\kern 1pt} '\rho D\left( T \right)}}} \right)} \mathord{\left/ {\vphantom {{\left( {\frac{{{{\varepsilon }_{0}} - {{\varepsilon }_{\infty }}}}{{2\varepsilon {\kern 1pt} '\rho D\left( T \right)}}} \right)} {{\text{l}}{{{\text{g}}}^{2}}\left( {\frac{{{{\omega }_{{{\text{min}}}}}}}{{{{\omega }_{{{\text{max}}}}}}}} \right)}}} \right. \kern-0em} {{\text{l}}{{{\text{g}}}^{2}}\left( {\frac{{{{\omega }_{{{\text{min}}}}}}}{{{{\omega }_{{{\text{max}}}}}}}} \right)}}, \\ {\text{где}}\,\,{{\omega }_{{{\text{max}}}}} = \sqrt {\frac{{A - 2 + \sqrt {{{A}^{2}} - 8A} }}{{2{{\tau }^{2}}\left( {A + 1} \right)}}} . \\ \end{gathered} $

На рис. 6 приведены результаты такого комбинированного расчета $\varepsilon {\kern 1pt} '$ и ${{{\text{tg}}\delta } \mathord{\left/ {\vphantom {{{\text{tg}}\delta } \rho }} \right. \kern-0em} \rho }$ в зависимости от частоты для образца 66041 при температурах на поверхности 150, 200, 350 К. При уменьшении температуры, область дисперсии смещается к низким частотам, при которых поглощение увеличивается, при этом поглощение коротких волн уменьшается.

Рис. 6.

Диэлектрические характеристики $\varepsilon {\kern 1pt} '$ и ${{{\text{tg}}\delta } \mathord{\left/ {\vphantom {{{\text{tg}}\delta } \rho }} \right. \kern-0em} \rho }$ реголита образца 66 041 в зависимости от частоты: линии 1 соответствует температура 350 К, линии 2 – 200 К, 3 – 150 К.

ЭЛЕКТРОФИЗИЧЕСКАЯ МОДЕЛЬ ГРУНТА

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

1. Слой однородного реголита толщиной 4 м, лежит на однородной коренной породе.

2. Плотность реголита меняется по глубине по формуле (1), плотность на поверхности принята равной ${{\rho }_{0}} = 1.5$ г/см3, плотность коренного грунта равна 2.1 г/см3.

3. Температура на поверхности меняется от времени суток согласно рис. 1 (Vasavada и др., 1999).

4. Температура верхнего слоя зависит от температуры на поверхности и меняется по глубине в соответствии рис. 2. Толщина слоя, подверженного суточным изменениям, считается равной 1 м.

На рис. 7–9 показаны результаты расчетов $\varepsilon {\kern 1pt} '$ и tgδ верхнего слоя грунта Луны с учетом условий 1–4. Заметим, что плотность грунта меняется с глубиной, поэтому графики на этих рисунках построены для ${\text{tg}}\delta $, а не для tgδ/ρ как на графиках выше. На рис. 7 показано распределение по глубине действительной части диэлектрической проницаемости $\varepsilon {\kern 1pt} '$ и тангенса угла потерь ${\text{tg}}\delta $ верхнего слоя грунта Луны для частот 1 МГц (линия 1) и 1 ГГц (линия 2) при температуре на поверхности (а) 80 К и (б) 380 К. Из расчетов следует, что изменение плотности реголита по глубине не существенно влияет на глубинное распределение диэлектрических характеристик, в отличие от изменения температурного режима. В геологически однородном верхнем слое в течение лунного дня действительная часть диэлектрической проницаемости меняется от 1.4 до 2.7, а ${\text{tg}}\delta $ – от 10–6 до почти 1. При ночной температуре внешний слой реголита вымораживается, значения $\varepsilon {\kern 1pt} '$ и ${\text{tg}}\delta $ в этом слое снижаются. Для волн разной частоты поведение диэлектрической проницаемости $\varepsilon {\kern 1pt} '$ и ${\text{tg}}\delta $ с глубиной – различны.

Рис. 7.

Распределение по глубине $\varepsilon {\kern 1pt} '$ и ${\text{tg}}\delta $ для образца 66041 при температуре на поверхности (а) 80 К, (б) 380 К. Линии 1 – 1 МГц , 2 – 1 ГГц.

Рис. 8.

Глубинные профили $\varepsilon {\kern 1pt} '\left( z \right)$ верхнего слоя реголита: (а) частота 1 МГц, (б) частота 1 ГГц. Температура на поверхности: 1 – 380 К, 2 — 300 К, 3 – 230 К, 4 – 150 К, 5 – 80 К.

Рис. 9.

Глубинные профили (z) верхнего слоя грунта: (а) частота 1 МГц; (б) частота 1 ГГц. Температура на поверхности 1 – 380 К, 2 — 300 К, 3 – 230 К, 4 – 150 К, 5 – 80 К.

На рис. 8а приведены более детальные глубинные профили действительной части диэлектрической проницаемости верхнего слоя реголита толщиной 1 м для разного времени суток на частоте 1 МГц, на рис. 8б – 1 ГГц. На рисунке 9а приведены глубинные профили мнимой части диэлектрической проницаемости верхнего слоя реголита толщиной 1 м для разного времени суток на частоте 1 МГц, на рис. 9б – 1 ГГц. По горизонтальной оси графиков отмечена глубина в метрах, по вертикальной – безразмерные величины $\varepsilon {\kern 1pt} '$ (а) и (б). Линия 1 соответствует температуре на поверхности 380 К, линия 2 – 300 К, 3 – 230 К, 4 – 150 К, 5 – 80 К.

ЗАКЛЮЧЕНИЕ

В работе изложена методика построения электрофизической модели верхнего слоя грунта Луны на основе имеющихся данных о его минералогическом составе и представлений о свойствах диэлектрической проницаемости диэлектриков в диапазоне частот ${{10}^{2}}{\kern 1pt} - {\kern 1pt} {{10}^{{10}}}$ Гц. Параметры, используемые в формулах при расчете диэлектрических характеристик, выведены на основании анализа лабораторных измерений диэлектрической проницаемости, проведенных с образцами грунта, доставленными на Землю различными авторами (Heiken и др., 1991). Методика иллюстрируется графиками, построенными для образцов реголита 66041 (Apollo-16) и анортозита 60025 (Apollo-6). Модель учитывает зависимость диэлектрической проницаемости и тангенса угла потерь от плотности и температуры грунта на поверхности, а также от частоты зондирующего электромагнитного поля.

Из анализа расчетов следует, что диэлектрическая проницаемость верхнего слоя грунта Луны зависит от температуры на поверхности и от частоты электромагнитного поля, изменением плотности грунта по глубине при расчете диэлектрических свойств верхнего покрова можно пренебречь. Основной вывод, полученный в результате анализа расчетов и положений теории диэлектриков, состоит в том, что геологически однородный слой реголита, а иногда и находящиеся под ним коренные породы, в течение лунного дня меняют свои диэлектрические свойства за счет изменения температуры. Действительная часть диэлектрической проницаемости меняется от 1.4 (значение ${{\varepsilon }_{0}}$) до 2.7 (${{\varepsilon }_{\infty }}$), а ${\text{tg}}\delta $ – от 10–6 до почти 1 (${{\left( {{{\varepsilon }_{0}} - {{\varepsilon }_{\infty }}} \right)} \mathord{\left/ {\vphantom {{\left( {{{\varepsilon }_{0}} - {{\varepsilon }_{\infty }}} \right)} 2}} \right. \kern-0em} 2}$). Изменения происходят на восходе и закате Солнца, скорость изменения диэлектрической проницаемости и толщина слоя, в котором эти изменения проявляются, зависят от частоты взаимодействующего электромагнитного поля. При этом на частотах волн, длина которых превышает толщину теплоизолирующего слоя в 2–3 раза, суточные изменения диэлектрических характеристик верхнего покрова не заметны, а для более коротких волн изменения действительной части диэлектрической проницаемости и ${\text{tg}}\delta $ с изменением температуры сравнительно равномерны, что позволяет изучать теплоизоляционные свойства грунта при проведении измерений в одном месте в течение суток. Для зондирования Луны радиоволнами мегагерцового диапазона следует учитывать частотную изменчивость диэлектрических параметров не только верхнего слоя толщиной 1–2 м, но и грунта, лежащего ниже слоя реголита, а при использовании радиоволн гигагерцового диапазона помнить, что при низких температурах $\varepsilon '$ и ${\text{tg}}\delta $ уменьшаются настолько, что реголит можно считать для них “прозрачным”.

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

Работа выполнена в рамках государственного задания Института радиотехники и электроники им. В.А. Котельникова РАН.

Список литературы

  1. Кибардина И.Н., Юшкова О.В. Определение плотности и процентного содержания оксидов металлов грунта Луны по радиолокационным данным // Космич. исслед. 2017. Т. 55. № 3. С. 201–206.

  2. Кротиков В.Д., Троицкий В.С. Радиоизлучение и природа Луны // Успехи физ. наук. 1963. Т. LXXXI. № 4. С. 589–639.

  3. Маров М.Я., Колесниченко А.В., Мануйлов К.К., Осипов В.П. Эксперимент “ТЕРМОФОБ”: Прямые исследования теплофизических свойств грунта Фобоса// Астрон. вестн. 2010. Т. 44. № 5. С. 393–402.

  4. Минчин С.Н., Улубеков А.Т. Земля-космос-Луна. М.: Машиностроение, 1972. с. 244.

  5. Ржевский В.В., Шварев В.В., Силин А.А. Исследование физических свойств грунта “Луны-20” и его земных аналогов // Космич. исслед. 1976. Т. XIV. В. 2. С. 187–292.

  6. Сканави Г.И. Физика диэлектриков. Область слабых полей. М.-Л.: ГИТТЛ, 1949. С. 1–500.

  7. Слюта Е.Н. Физико-механические свойства лунного грунта // Астрон. вестн. 2014. Т. 48. № 5. С. 358–382.

  8. Тареев Б.М. Физика диэлектрических материалов. М.: Энергия, 1982. с. 302.

  9. Черкасов И.И., Шварев В.В. Грунт Луны. М.: Наука, 1975. С. 1–144.

  10. Флоренский К.П., Базилевский А.Т., Николаева О.В. Лунный грунт: свойства и аналоги // Академия наук СССР: Институт геохимии и аналитической химии им. В. И. Вернадского, 1975. С. 1–50.

  11. Фрелих Г. Теория диэлектриков. М.: Издательство иностранной литературы, 1960. С. 1–251.

  12. Шевченко В.В. Луна // Физическая энциклопедия / Советская энциклопедия, 1990. Т. 2. С. 613–615.

  13. Юшкова О.В., Кибардина И.Н. Диэлектрические свойства верхнего покрова Луны // Астрон. вестн. 2017. № 2. С. 136–141.

  14. Bassett H.L., Shackelford R.G. Dielectric properties of Apollo 14 lunar samples at micrometer and millimeter wavelengths// Proc. Lunar Sci. Conf. 3rd. 1972. P. 3157–3160.

  15. Chung D.H., Westphal W.B., Olhoeft G.R. Dielectric properties of Apollo 14 lunar samples // Proc. Lunar Sci. Conf. 3rd. 1972. P. 3161–3172.

  16. Frisillo A.L., Olhoeft G.R., Strangway D.W. Effects of pressure, temperature and density on the electrical properties of lunar samples 72441,12, 15301,38 and a terrestrial basalt// Earth Planet. and Sci. Lett. 1975. P. 345–356.

  17. Gold T., Bilson E., Yerbury M. Grain size analysis, optical reflectivity measurements, and determination of high-frequency electrical properties for Apollo 14 lunar samples // Geochim. et Gosmochim. Acta. 1972. 3. Suppl. P. 3187–3193.

  18. Gold T., Bilson E., Yerbury M. Grain size analysis and determination of high frequency electrical properties of Apollo 15 and 16 samples // Proc. Lunar Sci. Canf. 4th. 1973. P. 3093–3099.

  19. Gold T., Bilson E., Baron R. L. Electrical properties of Apollo 17 rock and soil samples and a summary of the electrical properties of lunar material at 450 MHz frequency // Proc. Lunar Sci. Conf. 7th 1976. P. 2593–2603.

  20. Heiken G., Vaniman D., French B.M. Lunar Sourcebook: A User’s Guide to the Moon. Cambridge: Cambridge Univ. Press, 1991. P. 1–736.

  21. Huang S. Surface temperatures at the nearside of the Moon as a record of the radiation budget of Earth’s climate system // Adv. Space Res. 2008. V. 41. P. 1853–1860.

  22. Katsube T.J., Collett L.S. Electrical properties of Apollo 11 and Apollo 12 lunar samples// Proc. Lunar Sci. Conf. 2nd. 1971. P. 2367–2379.

  23. Katsube T.J., Collett L.S. Electrical characteristics of rocks and their application to planetary and terrestrial EM-sounding // Proc. Lunar Sci. Conf. 4th. 1973. P. 3111–3131.

  24. Olhoeft G.R., Strangway D. W. Dielectric properties of the first 100 meters of the Moon // Earth and Planet. Sci. Lett. 1975. V. 24. P. 394–404.

  25. Paige D.A., Foote M.C., Greenhagen B.T., Schofield J.T., Calcutt S., Vasavada A.R., Preston ·D.J., Taylor F.W., Allen C.C., Snook K.J., Jakosky B.M., Murray B.C., Soderblom L.A., Jau B., Loring S., ·Bulharowski J., Bowles N.E., Thomas I.R., Sullivan M.T., Avis C., De Jong E.M., Hartford W., McCleese D.J. The Lunar Reconnaissance Orbiter Diviner Lunar Radiometer Experiment // Space Sci. Rev. 2010. V. 150. P. 125–160.

  26. Vasavada A.R., Paige D.A., Wood S.E. Near-surface temperatures on Mercury and the Moon and the stability of polar ice deposits // Icarus. 1999. V. 141. P. 179–193.

  27. William G., Thomas D.K. Phenomenological and molecular theories of dielectric and electrical relaxation of materials // Novocontrol Application Note Dielectric 3, Novocontrol GmbH, 1998. www.novocontrol.de/pdf_s/APND3.PDF

  28. Williams J.-P., Paige D., Greenhagen B. T., Sefton-Nash E. The global surface temperatures of the Moon as measured by the diviner lunar radiometer experiment // Icarus. 2017. V. 283. P. 300–325.

  29. Xiao L., Zhu P., Fang G. et al. A young multilayered terrane of the northern Mare Imbrium revealed by Chang’E-3 mission // Science. 2015. V. 347. № 6227. P. 1226–1229.

  30. Zheng Y.C., Tsang K.L., Chan K.T., Zou Y.L., Zhang F., Ouyang Z.Y. First microwave map of the Moon with Change-1 data: The role of local time in global imaging // Icarus. 2012. P. 194–210. https://doi.org/10.1016/j.Icarus.2012.02.017

Дополнительные материалы отсутствуют.