Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Курганов Р.А. Прогнозирование наклонного рассеивания радиоволн метеорными ионизациями

.pdf
Скачиваний:
7
Добавлен:
23.10.2023
Размер:
8.32 Mб
Скачать

регистрируемая электронная плотность переуплотнен­ ных следов:

 

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

î

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

Соседние файлы в папке книги из ГПНТБ