книги из ГПНТБ / Курганов Р.А. Прогнозирование наклонного рассеивания радиоволн метеорными ионизациями
.pdfрегистрируемая электронная плотность переуплотнен ных следов:
|
U\ -гх-г2 |
(г, — r2) (1 — cos2 ß-sins Ф) |
|
Т 2 |
|||
а |
|
|
|||||
|
|
|
|
|
|
|
|
|
4 « / |
\ е |
/ |
32и4 |
G r G 2 - c o s 2 |
fx |
|
|
|
|
|
||||
|
|
|
|
|
|
|
(21) |
Зная ao m o , можно |
рассчитать |
плотность |
дифферен |
||||
циального |
распределения |
длительностей |
р(Т) |
для |
|||
элементарного метеорного потока |
d<o, регистрируемого |
||||||
в объеме |
dV = dty-dQ. |
|
|
|
|
||
Для этого весь диапазон регистрируемых длитель |
|||||||
ностей разбивается |
на |
интервалы: 0—0,1—0,2—0,4 |
0,8-1,6-3,2-6,4-12,8-25,6-51,2-102,4 - 204,8 сек. Наименьшая регистрируемая в объеме dV максималь ная линейная электронная плотность следа аопгѵ со ответствует случаю отражения от точки вблизи
максимума ионизации |
следа, |
находящегося на высоте |
||||||
homv сигнала, |
равного |
пороговому уровню. Для отра |
||||||
жения от следов |
с такой |
электронной |
плотностью |
|||||
время превышения |
порогового уровня |
ТА |
= 0 |
и диапа |
||||
зон высот, при отражении в пределах |
которого |
|||||||
отраженный сигнал выше порогового уровня, |
||||||||
Для o.mv > aomv |
сигнал с |
длительностью |
превыше |
|||||
ния порогового уровня ТА^0 |
|
будет |
создаваться при |
|||||
отражении от следов с зеркальной |
точкой, находя |
|||||||
щейся на удалении до + AhlroЫі2Т0 |
от высоты мак |
|||||||
симальной ионизации |
следа hmv, |
лежащей |
в диапазоне |
|||||
высот от /іІТ0{а.тѵ) |
до |
hm){amv), |
причем |
А ш — А2 Л ) = |
||||
= ДЛП) по определению. |
|
|
|
|
|
|||
Значения |
hiT0, |
Іі2Т0 |
определяются |
методом после |
||||
довательных |
приближений, для чего |
вычисляется ве |
личина ТА — времени превышения сигналом порогово го уровня при его отражении от точки, находящейся на некотором расстоянии ДАЛ от высоты hmv макси мальной ионизации следа с amv.
50
Если |
принять Герлофсонское |
распределение |
элек |
|||||
тронной |
плотности |
вдоль |
следа, то электронная |
плот |
||||
ность в точке отражения |
на высоте |
hmv |
+ àhA |
будет |
||||
|
|
_ |
_9_ |
|
А |
. |
Л |
|
|
|
|
|
|
|
|
||
При а,, |
, д / , |
< а,0, |
т. е. |
критической |
электронной |
|||
тѵ |
А |
|
|
|
объема |
|
|
|
плотности для данного элемента |
|
|
||||||
|
at 0 = |
7,5-1013-exp |
"тѵ+ |
AhA |
|
|
||
|
ЗА2 |
sec2 |
Ф |
|
|
|||
|
|
|
|
|
|
след в точке отражения является недоуплотненным. Следовательно, амплитуда отраженного сигнала падает под действием амбиполярной диффузии по экспонен циальному закону и время превышения порогового уровня *
|
|
|
|
Wrh |
+ ДЛ |
|
|
|
|
|
|
п,тѵ |
+ à"A |
|
7V = х -In |
|
|
|
|
|
|
|
|
• exp |
X2 |
-sec2 Ф |
|
|
|
|
|
|
||
При |
а. ть |
А > ac0 след |
в точке |
|
отражения — пере |
|
уплотненный |
и время превышения |
порогового уровня |
||||
при |
условии |
разрушения |
следа только |
за счет амби |
||
полярной диффузии |
|
|
|
|
^ = 1 , 1 3 . 1 0 - 1 4 . х г а Л / п г ) + Д Л л .
|
X2 «sec2 Ф |
Здесь х„ = |
постоянная распада отражения за |
16я2-£>
счет амбиполярной диффузии. Коэффициент амбиполяр ной диффузии по данным [37] с погрешностью менее 10% аппроксимируется выражением
0,078-(Лт г , + ДЛд) — 6,7
D = e x P |
Шз |
* |
4* |
51 |
При одновременном разрушении следа за счет эффекта прилипаний с постоянной распада
|
•сл = 40-ехр |
|
9 6 - ( / * , „ „ + ДЛ А ) |
|
|
||||
|
|
н |
|
|
|
||||
|
|
|
|
|
|
|
|
||
время превышения |
порогового |
уровня |
уменьшается |
||||||
до величины Т'„ , определяемой |
из |
уравнения |
|||||||
|
|
|
Тп |
|
• ехр |
|
|
|
|
по формуле |
|
|
|
|
|
|
|
|
|
|
|
Тп |
|
2,3 - х„ + Т |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
дающей |
решение |
этого |
|
трансцендентного |
уравнения |
||||
|
|
|
|
|
|
|
|
|
т |
с погрешностью менее 20% для значений 0< |
< 100 , |
||||||||
соответствующих |
всему |
диапазону |
|
|
in |
||||
регистрируемых |
|||||||||
на трассах линейных электронных плотностей. |
|||||||||
Выше, на основании |
|
анализа |
экспериментального |
||||||
материала, установлено, |
|
что значительная |
часть ме |
||||||
теорных радиоотражений длительностью |
больше 3 сек |
||||||||
является |
федингующей. Процентный вклад |
числа фе |
|||||||
дингующих отражений является |
в |
первом |
приближе |
||||||
нии функцией длительности отражения |
Тп. |
|
|||||||
|
/ = 1 -exp[(6.10- 3 -X.sec<£ + |
|
|
||||||
, |
+ 0,63). Г , |
-2,2 . }0~2 . X • sec Ф + 1,73 |
|
(22) |
|||||
|
|
|
|
]. |
|
а относительное время превышения порогового уров ня федингующим сигналом является функцией его амплитуды и формы фединга
/ = 1 |
• f <hmv + àhA ) |
•ехр |
|
X2.sec2<I> |
|
amv |
|
|
|||
a |
\ 0,217 |
„ о с |
Л Л Л |
|
|
"omv \ |
0,65—jj- |
1 |
(23) |
||
|
|
|
|
Тогда искомое время превышения сигналом порогово го уровня будет равно средней длительности превы-
52
ці-ения порогового уровня |
отражением, часть |
которо |
|||||||
го федингует |
|
|
|
|
|
|
|
|
|
|
T^T'n[(\~f) |
+ t.f\. |
|
|
|
|
(24) |
||
Очевидно, что hm |
и Л 2 Г о ~ э т о |
высоты hmv |
+ |
AhlTÙ |
|||||
и hmv — Дй2 Г 0 , |
при которых |
отраженный |
сигнал |
имеет |
|||||
длительность |
Г л = 0. |
|
Все вычисленные |
в процессе |
|||||
приближения |
значения |
ТА |
0 являются |
полезными, |
|||||
так как соответствующие |
им высоты |
hmv |
+ àhA |
есть |
|||||
значения h1TA, |
h2TA, |
а |
разность |
hXTA |
— h2TA, |
|
равная |
ДЛЛ — шагу, принятому для h в процессе приближения, является толщиной метеорного слоя для отражений от метеоров элементарного,потока длительностью от
Тд |
до ТА+1. |
Тогда часть элементарного потока |
час |
|||
тиц с максимальной электронной плотностью amv |
соз |
|||||
даваемых |
следов, отражающая сигналы длительности |
|||||
от ГА до 7Ѵ,+ 1 , будет |
|
|
|
|||
|
cfyV=Q(aZ 0 ).a|-1 • cos^'Z- sinß-a-*. |
|
|
|||
|
|
|
{hmv + |
ДА « ) ' Д А л |
|
|
|
.p(v) |
m v |
А—à--da-dvde-du. |
(25) |
||
|
|
|
cos2 e-si-n2^ |
v |
' |
|
Так как отражающие точки распределены |
по высоте |
|||||
равномерно, а при достаточно малом ДАД можно |
счи |
|||||
тать |
ТА |
линейно-зависящей от h, то dN будет равно |
||||
мерно распределено по Т в интервале от TА |
до |
Г л + 1 |
||||
и должно |
быть занесено в графы таблицы |
дифферен |
циального распределения длительности в форме вкла
дов dNk = |
— d j V ' g f e |
* |
для k1 < k < k2, |
где kx соответ- |
|
TA + l ~ T |
A |
табличному значению Th > |
|
ствует наименьшему |
||||
> т і п [ Г Л , |
TA+Ï\, |
k2 |
— соответствует |
наименьшему |
значению |
|
|
|
|
для |
kl<k<.k2, |
okx = 7\ — min [ГА , Гд-! ], |
53
Сумма вкладов в „£" графу таблицы отражений
для всех ДЛл , « т |
> % т ѵ , г», ѳ , |
ß, ф равна |
УѴА — числу |
от |
||
ражений на трассе с длительностью от |
до |
7"Ä. |
||||
Общее |
число |
отражений |
на трассе, |
т. е. |
часовое |
|
|
12 |
|
|
|
|
|
число |
N—Y*Nk, |
плотность |
дифференциального |
рас- |
||
пределения длительностей |
|
|
|
|
||
|
|
р(Т )= |
N k |
|
|
|
k' # . ( 7 * - 7 * _ і )
исредняя длительность отражений
12 |
• |
LÀ
î
i n
Изложенная методика расчета может быть приме нена для прогноза на трассах длиной больше 1000 км величин часового числа отражений выше заданного уровня, количества отражений выше заданной длитель ности, а следовательно, величины средней длитель ности и коэффициента заполнения, разброса длины пробега отраженных сигналов, а также угловой на правленности принятых сигналов. Время расчета этих параметров для одного часа на ЭВМ М-20 порядка 10—12 минут.
1.2.3. Метод прогноза объемной плотности отражающих точек при метеорном распространении радиоволн
Наиболее общим методом прогноза основных пара метров, характеризующих метеорное распространение радиоволн, применимым для трасс любой длины, а также позволяющим исследовать тонкую структуру
. распределения в пространстве плотности отражающих точек, является метод расчета объемной плотности отражающих точек. Объемная плотность отражающих точек есть число метеорных следов, создающих в единице объема, расположенного вокруг данной точ ки метеорной зоны ионосферы, отражения, превышаю щие пороговый уровень в точке приема на время, больше некоторого заданного.
54
Примем в качестве основной декартовую систему координат, связанную с опорной плоскостью — плос костью, проходящей через точки приема и передачи перпендикулярно земному радиусу, проведенному в точку центра трассы. Ось х направим по линии трас сы в сторону приемника, ось _у — перпендикулярно оси X в опорной плоскости, ось z перпендикулярно опорной плоскости. Начало координат пометим в центре трассы.
Объемная плотность отражающих точек элемента объема dVс координатами центра х, у, z может быть определена как частное от деления суммарного ме теорного потока, входящего за единицу времени через поверхность Е, ограничивающую элементарный объем, и дающего отражения внутри него на величину эле ментарного объема
A ( J C | у, |
z ) = d " ^ У. ') , |
|
||
v |
^ |
' |
dV |
|
где |
|
|
|
|
d N = = z И*"" * ' d w ' d a ' |
' P(v)dv- |
С2 6 ) |
||
Здесь t • s-dm• do'ti(s) • p(v)• |
dv — элементарный |
ме |
||
теорный поток |
метеоров |
данного направления |
и |
скорости, входящий в объем dV через бесконечно малый элемент его поверхности do и создающий в объеме отражения, превышающие пороговый уровень в точке приема на время больше заданного, dm — вели чина телесного угла, в пределах которого падает элементарный поток, s—вектор радианта элементарного
потока, |
t — нормаль |
к |
элементу |
поверхности do, |
||
n(s) — плотность метеорного |
потока данного |
направ |
||||
ления |
и скорости, создающего в элементе объема |
|||||
следы, |
отражающие из объема в точку приема |
радио |
||||
сигналы, превышающие |
пороговый |
уровень: |
|
|||
|
|
|
со |
|
|
|
|
n(s) |
= b (s)• J |
m~sdm, |
|
||
где tn0 |
— минимальная регистрируемая масса метеоров |
данного направления и скорости, входящих в элемент объема через do и b (s) — плотность радиантов для
55
направления |
s. Величина минимальной |
регистрируе |
|
мой массы метеоров данного направления |
и |
скорости |
|
определяется |
условием равенства мощности |
или ам |
плитуды радиосигнала, отраженного от следа данного метеора, некоторой пороговой мощности или ампли туде в точке приема. Следовательно, решение задачи нахождения объемной плотности метеорных радиоот ражений сводится к определению минимально ре гистрируемой массы т0, задаче энергетической и
определению метеорного потока известной плотности, имеющего отражающие точки в данном объеме, зада че чисто геометрической.
Для каждой точки метеорной зоны, характеризуемой координатами je, у, z, существует только одна зеркаль ная плоскость—плоскость, касательная к одному из се мейства эллипсоидов, имеющих фокусами точку пере дачи и приема. Пересечение зеркальной плоскости с плоскостью местного горизонта дает зеркальную по луплоскость, лежащую над местным горизонтом. Лю бой метеор, траектория которого лежит в этой зер кальной полуплоскости, создает отраженный сигнал в точке приема.
Амплитуда отраженного сигнала зависит от линей ной электронной плотности в точке отражения а (1), ä последняя пропорциональна массе метеорной части цы (2).
Произведя в выражениях (1), (2) замену
, |
-г -у ( |
{ |
Л-ѵ2А |
} |
\ 3 |
, _ |
|
« = p3,.cos |
Z |
|
— |
|
|||
_ JL |
|
|
|
_ |
|
|
|
Р = е " . |
Рт = е |
" . |
|
получим
а = 5,13-КГ1 4 |
cos 2Z-p |
(pn-j-.pj.v9. |
|
|
(27) |
Приравняв это выражение к величине минимальной регистрируемой в элементе объема электронной плот-
56
ности (20), (21) и решив полученные уравнения отно сительно рт, найдем
c o s 2 Z |
4((хЯ) |
0,5 |
• (28) |
|
р.ѵ9 |
|
АЛ |
|
|
9-10-'5 |
|
|
||
|
2-l.g |
|
|
Очевидно, что метеоры данного направления ß, траек
тории которых |
лежат в зеркальной плоскости, соот |
ветствующей |
выбранным координатам х, у, z и х, у, |
h, создадут в точке приема сигнал, превышающий по |
|
роговый уровень UQ, если давление в точке их макси |
|
мального испарения р , „ > р т о , где ^ — с о о т в е т с т в у е т |
летящему со скоростью ѵ метеору минимальной мас сы т0, создающему сигнал, равный пороговому.
Значение |
тп определяется как |
|
|
|
||||||
|
|
|
/ |
Л . Л 4 \ 3 _ _ _ _ з |
g |
|
|
|
||
|
|
|
А.у^А V |
|
|
|
|
(29) |
||
|
т |
0 |
= V |
'21-g ) ' COS |
|
|
|
|
||
|
|
|
|
|
|
|
||||
Тогда плотность метеорного потока |
данного |
направ |
||||||||
ления |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
оо |
|
|
|
|
|
|
|
|
n{s) = b- |
^m~s-dm |
|
|
|
|||
|
b(s) |
|
АА |
2 |
|
- і З ( і - і ) |
|
|
||
n{s) |
|
|
|
\_2-l-g |
V |
PmOJ |
|
|
(30) |
|
|
|
|
cosZ |
|
|
|||||
Для определения метеорного |
|
потока |
dN (х, у, z)> |
|||||||
связанного |
с элементом |
объема |
dV, |
окружающим |
||||||
точку (х, у, |
z), |
воспользуемся методикой, |
предложен |
|||||||
ной Н. Каррара [91]. |
примем „объем части эле |
|||||||||
За элемент |
|
объема dV |
||||||||
ментарного |
конуса |
с апертурой |
|
2üfa, |
соосного |
с бис |
||||
сектрисой угла рассеяния 2Ф, соответствующего |
центру |
элементарного объема, ограниченный двумя сфери ческими поверхностями радиусов р и р + dp, имеющими центр в точке пересечения биссектрисы угла рассеи вания с осью трассы (рис. 10)". Величина телесного угла, в пределах которого падает метеорный поток направления s, проходящий через бесконечно малый элемент поверхности аз и дающий отражения в пре делах элементарного объема, du) = e ' ( ß ' , P)d§'. Здесь
57
|
Рис. |
10. |
|
ß'— угол между нормалью |
t к поверхности элемента |
||
do и вектором |
радианта s, |
лежащим в |
плоскости іг, |
касательной к |
зеркальному |
эллипсоиду, |
проходящему |
через точку Р центра элемента dV. s' — угол в пло скости, перпендикулярной плоскости тс и проходящей через направление на радиант s, в пределах которого
расположены метеорные |
следы, дающие отражения |
||
в пределах объема |
dV. |
|
|
Тогда |
выражение |
(26) |
приобретает вид |
dM(x, у, |
z) = |
р2_ |
^d3§s-ts'($',P)-d$^n(s)-p(v)-dv,(3l) |
|
где s - 7=cosß' . При пренебрежении величинами вто рого порядка малости в [91] получено следующее вы ражение для е'
|
e' = |
2ufa.cosß' 1 |
- / • О - |
£ ) ] • |
|
< 3 2 > |
|||
где а и Ь — полуоси |
зеркального |
эллипсоида, прохо |
|||||||
дящего |
через |
точку |
Р, |
/ — направляющий |
косинус |
||||
вектора |
радианта s — U + mj+ nk. |
|
|
|
|||||
Интегрированием |
по |
ß' |
в пределах ф\ <$г |
<.ф'2, со |
|||||
ответствующих метеорному |
потоку, входящему в эле |
||||||||
мент объема через |
элемент |
поверхности |
do- опреде |
||||||
ляется связанный с da метеорный |
поток, |
создающий |
|||||||
в пределах элементарного |
объема |
отражения, превы |
|||||||
шающие |
пороговый уровень |
в точке приема. |
|
Используя выражения (31), (32), решим задачу об определении объемной плотности отражающих точек следующим образом.
58
Свяжем с центром элемента объема dV вспомога тельную декартовую систему координат, имеющую
[n-Z] |
|
орты и, v. Орт и =- cos Zr |
, где п — нормаль к по |
верхности зеркального эллипсоида, проходящего через точку Р,
п = —( — і + -I- / + — k),
x, y, z — координаты точки P, Z — единичный вектор зенита в точке Р (рис. 11)
х-і +y-j + f/?3-cos |
+ z |
Рис. 11.
Z m l n — минимальный зенитный угол для метеорных следов, лежащих в данной зеркальной плоскости. Очевидно, что и совпадает с линией пересечения зер кальной плоскости, проходящей через точку Р, с плоскостью горизонта.
59