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

книги из ГПНТБ / Светов Б.С. Теория, методика и интерпретация материалов низкочастотной индуктивной электроразведки

.pdf
Скачиваний:
2
Добавлен:
24.10.2023
Размер:
11.24 Mб
Скачать

Рис. 1. Эллипс поляризации на прост­ ранственной (вверху) и фазовой (внизу) плоскостях

Характеризующие их коми-

v

лексные величины # ж е _ г ф и

Нуе-*ч изображаются на комп­ лексной (временной или фазо­ вой) плоскости векторами, вра­ щающимися с угловой частотой ©=ф/^. Если частота считается заданной, то напряженности поля &€х и У могут быть пол­ ностью охарактеризованы их комплексными амплитудами

Нх = Нхе-*ч'х

и

tfj

• Ну&-

изображаемыми

на

комплекс­

ной

фазовой

плоскости

(рис.

1, внизу)

неподвижными

векто­

рами,

модули

которых

равны

Нх

и

Ну,

а

аргументы

(углы

наклона к оси абсцисс) —вели­ чинами ф ж и —ц>у. Оси абс­ цисс и ординат на этой плос­ кости соответствуют начальным

я

фазовым углам ф = 0 и ф = — —

и обозначаются ORe и 01т . Следует отметить, что при при­ нятой нами комплексной записи изменения поля во времени в

виде Fe~iwt знаки у аргументов комплексных выражений поля про­ тивоположны знакам соответствующих фаз колебаний.

Комплексные

амплитуды компонент

V

V

можно рас­

поля Нх

и Ну

сматривать [4]

как ортогональные пространственные

координаты

комплексного вектора амплитуды напряженности (или

просто на­

пряженности) магнитного поля Н:

 

 

 

Н=х0Нх0Ну0Нхе-1**+у0Нуе-Ъ.

 

 

(II.6

При этом мгновенное значение вектора

определяется равенст­

вом

 

 

 

 

3Vf=Re[[x0Hx-\-yQH,

= Re

[не

(II.7)

20

В соответствии с равенством (П.7) комплексный вектор Н опре­ деляет мгновенное значение вектора напряженности поля &6<$ по лю­ бому направлению и в любой момент времени, т. е. полностью ха­ рактеризует переменное эллиптически поляризованное поле задан­ ной частоты оз. С этой точки зрения изучение эллиптически поляризованного поля эквивалентно изучению поля комплексного

вектора Н.

v

Следует отметить, что комплексный вектор Я не является векто­ ром в обычном смысле этого слова. Поскольку его «координаты»

v v

Нх и Ну имеют различные фазы (аргументы), их можно рассматри-

v

вать как проекции комплексного вектора Н на пространственные

оси координат лишь условно,

а сам вектор Н не допускает пред-

V

V

ставления ни в форме # = # е -

г ф ° , ни в форме Н = Нп (п — единич­

ный вектор с заданным направлением). Поэтому нельзя говорить ни о модуле и аргументе (фазе) комплексного вектора, ни о его величине и направлении. На комплексные векторы распростра­ няются основные векторные операции, применимые к обычным век­ торам. В частности, их можно подвергать любым линейным преоб­ разованиям, если окончательный результат находится в соответст­ вии с равенством ( I I . 7 ) . Нелинейные преобразования комплексных векторов требуют специального определения.

Выделив в комплексном векторе Н его действительную и мни­ мую части, можно представить его в виде, аналогичном выражению (П.6):

 

 

^ = R e ^ + / I m ^ = ^ R e + / 7 / I m )

(II.6')

где векторы # к е и # i m являются действительными векторами, зани­

мающими

некоторое

положение в пространстве, а 1 = У—1 указы­

вает направление оси ординат на фазовой плоскости.

 

Подставляя

выражение ( I I . 6 ' ) в равенство ( I I . 7 ) , получим:

 

i

^ =

R e

[(/7R e +/77I r a ) e - '' < p ]=F R e cosc P +77 I m sin T .

(II.

Это равенство является параметрическим уравнением эллипса поляризации, записанным в векторной форме в косоугольной си­ стеме координат, определяемой векторами #в.е и Н\т. Из равен­ ства (II.8) следует, что векторы # н е и Him совпадают с вектором ей/ф в моменты времени,фиксируемые значениями мгновенной фазы

21

<р'=0±2л;п и ф" = — — ± 2 л п соответственно. Их ортогональными координатами являются действительные (активные) и мнимые (ре­

активные) части комплексных амплитуд декартовых

компонент

поля &6х и &6у. В силу этого #Re и Н\ш

называют активным и реак­

тивным векторами напряженности поля.

 

 

 

 

Исключая

из параметрических

уравнений

эллипса

поляризации

( I I . 1 ) , (П.2)

параметр

<р, получим

уравнение

эллипса

в

декартовой

•ортогональной системе координат Оху:

 

 

 

 

 

 

&в\

3€у

1&€х&€у

c o s ч)ху

Л,

(П.9)

tf2Sin2^y

^ s i n 2 ^

 

НхНуЫп^ху

y

 

 

 

гдефжу^фу — ц>х.

Повернем систему координат Оху на такой угол уа, при котором уравнение (П.9) запишется в каноническом виде. Это произойдет тогда, когда cos<pxy станет равным нулю, а амплитуды напряженностей поля по осям Ох и Оу — равными амплитудам поля по на­ правлениям полуосей эллипса поляризации:

ij2

= 1 .

( Н Л О )

А/2 I

 

 

Здесь &6а и &6ъ —мгновенные значения напряженностейполя по на­ правлениям полуосей, а На и Нь — их амплитуды, равные соответ­ ственно большой и малой полуосям. Так как в равенстве (11.10) cosqXy = 0, то напряженности поля ач&6ъ сдвинуты по фазе на

угол -Ц- (квадратурны). Направления полуосей эллипса поляриза­ ции — единственные направления, по которым напряженности поля

•одновременно взаимно ортогональны и взаимно квадратурны. Уравнения (II.1) и (II . 2) задают напряженности поля в произ­

вольный момент времени, фиксируемый фазой ф, по двум ортого­ нальным фиксированным направлениям — Ох и Оу. С их помощью можно найти напряженность поля по любому направлению и в лю­ бой момент времени, т. е. полностью охарактеризовать эллиптиче­ ски поляризованное поле. Можно поступить и наоборот: зафикси­ ровать некоторые моменты времени (например, для значений мгно-

венной фазы ф = 0 ± 2 я п и ф" = — — ± 2 п п ) и задать в эти моменты

напряженности поля по любому направлению, определяемому теку­ щим углом у:

J £ R e =

tfReCOS(T-TRe);

(н.Г)

^ I m = t f I m c o s ( T - T l m ) .

(И-2 ')

Здесь индексы Re и I m указывают на выбранные моменты времени. Поскольку в эти моменты мгновенный вектор напряженности поля

22

совпадает с активным # R e

и реактивным # i m векторами поля, то-

Нле= |Яке|, # i m = \ Him\,

а углы уке и Yim определяют расположе­

ние этих векторов в пространственной системе координат Оху (см. рис. 1,а). С учетом равенства (II.8) уравнения (П.1') и (П . 2'), так же как и (П.1) и (П.2), полностью характеризуют переменное эл­

липтически поляризованное поле. Если зафиксировать

угол

у,

то

напряженности поля

<Шв.е и <Жш можно

рассматривать

как

вре-

 

 

V

на фазовой

менные координаты некоторого «текущего» вектора

плоскости:

 

 

 

 

 

ж,=зеъ+тХа=Нъ

C o s ( T - T R e ) + / t f l m • c o s ( T - T i m ) .

(н.з>

При изменении угла у от 0 до 2я вектор

У описывает на фазо­

вой плоскости эллипс, для которого выражения {ИЛ')

и ( I I . 2 ' )

яв­

ляются параметрическими уравнениями с параметром у. В отличие от ранее рассмотренного пространственного эллипса поляризации этот эллипс можно назвать фазовым эллипсом поляризации. В соот-

V

ветствии с равенствами (П.З') и (П.6') текущий вектор & € v пред­ ставляет собой комплексную амплитуду напряженности поля по на­ правлению, составляющему угол у с осью Ох, поэтому обозначение векторности величины на фазовой плоскости ORelm не отличается от ранее введенного обозначения комплексности « V » -

v

можно записать

выражение,.

Для комплексной амплитуды &6У

аналогичное (II.8) и представляющее

комплексную амплитуду

в пространственной плоскости Оху:

 

 

M^=HxQOS^+Hys\n-i.

(II.8')

Равенство (11.8') можно рассматривать как параметрическое уравнение фазового эллипса поляризации, записанное в векторной форме в косоугольной системе координат, определяемой располо-

V V

жением векторов Нх и Ну на фазовой плоскости.

Исключая из уравнений (П.1') и (П.2') параметр у, получим уравнение эллипса поляризации в декартовых прямоугольных коор­ динатах на фазовой плоскости:

где YReim = YRe—Yim- Будем поворачивать систему координат ORelm до тех пор, пока уравнение (П.9') можно будет записать в канониче-

ском виде. Это произойдет при условии У=~2~-

В этом случае уста­

новятся равенства

е%?в.е = Ша, (fflim=

<Шъ, Нъеа,

# i m = #b и

уравнение фазового

эллипса

поляризации

полностью

совпадает

с равенством (11.10). Из

совпадения

этих уравнений следует

идентичность пространственного и фазового эллипсов поляризации

23

•с точностью до ориентации в соответствующих

системах координат,

т. е. с точностью до взаимной замены Ya—*—фа

(фа фаза

большой

полуоси).

 

 

 

 

 

При

переходе

от параметрических уравнений

( I I . 1 ) ,

(П.2) и

( I I . 1 ' ) ,

( I I . 2 ' ) к уравнениям пространственного и

фазового эллип­

сов поляризации

в форме

(11.9) и ( I I . 9 ' ) утрачиваются

сведения

о начале отсчета фазовых

(11.9) или пространственных (П.9') углов.

В соответствии с этим индексы векторов ей?Ф и &6у, описывающих эти эллипсы, становятся неопределенными, а сами векторы уже не дают полной характеристики эллиптически поляризованного поля.

Вектор е$?(р устанавливает связь между направлениями вектора на­ пряженности поля в некоторый момент времени и его величиной, но не указывает значения мгновенной фазы ф, которому этот момент

соответствует. Точно так же вектор Ш<$ связывает амплитуду на­ пряженности поля по некоторому направлению с ее начальной фа­ зой, но не содержит сведений об угле у, которому соответствует это направление.

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

Шу в плоскости ORelm или между фазой ф в плоскости ORelm и век­ тором &6у в плоскости Оху. Для того чтобы по углу y в плоскости Оху найти значение амплитуды поля по этому направлению (модуль

Ну), достаточно спроектировать на это направление пространствен-

v

фазовой

ный эллипс поляризации. После этого по модулю Ну на

v

 

плоскости может быть найден сам вектор Ну, угол наклона

которого

к оси абсцисс даст соответствующее значение начальной фазы — ц>у.

v

И наоборот, по модулю вектора Ну можно найти в плоскости Оху такое направление, для которого проекция эллипса поляризации

равна | Ну |. Угол между этим направлением и осью Ох равен углу у. При этих преобразованиях следует помнить, что углы у и —cpv располагаются в одном и том же квадранте пространственной и фазовой систем координат, определяемых соответственно векторами

На, Нь и На, Нь. Нетрудно путем

аналогичных рассуждений уста­

новить взаимно однозначную связь

между вектором

и соответ­

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

24

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

v

жение комплексного вектора напряженности поля Н.

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

ORelm:

/

}

 

{

 

C 0 S t

f x y

 

[

Я* s i n 2

о '

tf2sinV

'

нх"уsin2

?ху

или Оху:

\

•*

ху

у

ху

 

 

 

 

 

 

 

 

 

(

 

1

 

1

C 0 S 7 R e l m

Вместо этих коэффициентов могут быть взяты другие эквивалент­ ные им тройки независимых величин, инвариантных на фазовой ПЛОСКОСТИ х, Ну К Ц)Ху) ИЛИ В пространстве (#Re, Him И YRelm). Наконец, в качестве таких троек могут быть приняты и наиболее

непосредственные

характеристики

эллиптически

поляризованного

поля, инвариантные

относительно

поворота системы ORelm — #„,.

Нь и уа

или системы Оху На, Нъ и фа . Среди этих величин особое

место занимают полуоси эллипса поляризации На

и Нь. Они не за­

висят от ориентации

как пространственной (Оху),

так и временной

(ORelm)

систем

координат и в этом смысле могут быть названы

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

Таким образом, форма эллипса поляризации поля устанавли­ вается с помощью двух величин, инвариантных относительно начала отсчета как пространственных, так и фазовых углов — На и Нь. Для определения формы и ориентации пространственного (фазового) эл­ липса поляризации можно использовать три величины, инвариант­ ные только на фазовой (пространственной) плоскости, —Нх , Ну>. Фжу (Нпе, Him, YReim) или совокупность этих величин с полными ин­ вариантами — На, Нь, уа а, Нь, фа). Наконец, полная характери­ стика эллиптически поляризованного поля может быть получена с помощью четырех величин, часть из которых инвариантна в про­ странстве, а часть — на фазовой плоскости (например, На, Нь,,

Уа, фа) • Следует отметить, что эллиптически поляризованное поле мо­

жно полностью охарактеризовать и величинами, не инвариантными

как в пространстве,

так

и по фазе.

К

таким величинам относятся,

например, активные

и

реактивные

составляющие ортогональных

компонент поля — Re

 

I m Нх, Re Нv,

I m Ну.

Представляет интерес установить некоторые соотношения между ортогональными компонентами поля в пространстве и на фазовой плоскости и непосредственными характеристиками эллиптически поляризованного поля а, Нь, фа и уа).

25-

Опуская несложные промежуточные выкладки [58], получим:

, = 4 ~ [1/"я|+Я 2 у +2Я^Я у 51П <?Х у +

Y~H2x+H2y-2HxHySin<?Xy];

 

 

 

(И.11)

+ К я | е + Я ? т -

2 Я К е

Я 1 т sinT R e I m ] ;

( I I . 1 Г )

=4" [VHx+H2y-\-2HxHysm

<?ху - 1 ^ / 4 + / / , - 2 Я , / / , sin

;

 

 

 

(Н.12).

/ / 6 = ^ - [ l A / / R e - T - / / P m + 2 / / R e / / I m S i n T R e I m

-

L +

Я 2

т -

2 Я К е Я 1 т

sinT R e I m

J;

 

 

H l s l n

+ Я У s l n 2<?y

 

^

= Я

>

в

2 ^

+ Я у СОз2, у

5

* 2 Т - -

^ e c o s 2 T R e + H L c o s 2 T ] m '

 

 

 

 

2 -r~-

cos

 

 

 

t g 2 T a =

И х

2

— ;

 

 

 

 

 

1

2JL

 

 

 

 

 

2 £ l i 2 - c o s T R e I m

 

 

t g 2 ? « =

Я к е

2

5

 

 

 

 

 

1

I m

 

 

Я^ = Я ^ cos2 Т а + я Ь т 2 Т а ; Hle=Hl cos2 сра й 2 з1п 2 сра;

Я2 = Я 2 з т 2 Т а + Я 2 с о 8 2 Т а ;

Я , 2 т = Я 2

sin 2

? а + Я 6 2

cos2 сра ;

tg

 

*T « +

ff"« ;

 

"

1

<gT R e l m =

# -

+

.

 

 

1

"

-

(Ц.12')

( 1 1 Л З )

( П Л З )

(U.14)

(Н.14')

(11.15)

(Ц.15')

(11.16)

(Ц.16')

(Ц.17)

(11.17')

Здесь параллельно выписаны соотношения на пространственной и фазовой плоскостях. Эти формулы не только позволяют произво­ дить взаимные пересчеты одних величин в другие, но и составляют основу для проектирования аналоговых (или дискретных) измери­ телей элементов электромагнитного поля. Так, например, измере­ ния большой и малой полуосей эллипса поляризации осуществля­ ются в аппаратуре ЭПП путем приема поля на систему из двух ортогональных датчиков и последующего вычисления в измеритель­ ной аппаратуре формул (11.11), (11.12). Формула (11.17) указывает на возможность определения фазового сдвига <$х у без проведения фазовых измерений. Аналогично формула (11.17') показывает пути и возможности измерения угла между активным и реактивным век­ торами поля без задания ориентировки системы координат в про­ странстве.

СВЯЗЬ ПЕРВИЧНОГО И ВТОРИЧНОГО ПОЛЕЙ С СУММАРНЫМ ЭЛЛИПТИЧЕСКИ ПОЛЯРИЗОВАННЫМ ПОЛЕМ

Первичным

полем некоторого возбудителя

называется

его

поле

в однородном

непроводящем пространстве.

Отклонение

от

этого

поля, обусловленное наличием проводящей среды, удобно

приписать

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

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

ризующих вторичное поле,

и служит единицей его измерения.

В связи с этим представляет

интерес получение некоторых общих

соотношений между вторичным и наблюдаемым полями. Предполо­ жим первичное Щй и вторичное & € 1 поля линейно поляризован­ ными и совместим пространственную плоскость Оху с плоскостью

векторов

и е^?1, а оси

абсцисс пространственной

Ох и фазовой

ORe систем координат направим по первичному полю.

 

Тогда

 

 

 

 

 

 

(II.18>

 

I f

C O S ^ + c p ' ) ,

(II.19)

причем вектор Я 1 составляет угол Y 1 С первичным полем (осью Ох).

27

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

первичное

и вторичное поля

0 , Я 1 , у1 и ф1 )

посредством

формул

 

 

Hl=(fff+(HX¥

 

cos2

т 1 + 2 Я ° Я !

cos т' cos ср1;

(11.20)

 

 

Я 2 > е = ( Я 0 ) 2 + ( / / 1 ) 2

cos2

с р ' + г я 0 / / 1

cos f

cos т';

(11.20')

 

 

 

,

 

Я 1

sin у1 cosт1

 

 

 

 

 

r>i\

 

 

 

t g

 

НО + Я1 cos tp1 cos 71

'

 

 

 

 

 

 

 

 

to-Y

 

Я ' cos yi sin f

 

'

 

 

 

(U9V\

 

 

 

l e l R e — яо +

Я ' c o s 71 cos <pl

 

 

 

V i ' - ^ W

 

 

 

 

Я у = Я >

sinT ' ;

 

 

 

 

 

(11.22)

 

 

 

Я 1 г

а = Я >

sin <pi;

 

 

 

 

 

(11.22')

 

 

 

 

 

сру =срь

 

 

 

 

 

 

(11.23)

 

 

Яа =4"

 

 

T , m = T ' ;

 

 

 

 

 

 

 

(H.23')

 

 

( Я ° ) 2 + ( Я ' ) 2

+ 2 Я ° Я ' cos (f1

-

 

 

 

 

 

 

+ " | / " ( Я ° ) 2 + ( Я 1 ) 2 + 2 Я ° Я 1

c o s ( f ' - f ср')];

 

 

(П.24)

 

 

Яй =4" [ К ( Я ° ) 2 + ( Я 1

) 2

+ 2 Я ° Я 1 cos (т1 -

91 )

-

 

 

 

 

- V ( Я ° ) 2 + ( Я 1

) 2 + 2 Я « Я 1

cos (т' +

Т1 ) ] ;

 

 

(И-25)

 

 

( Я ' ) 2 5 1 п 2 < р 1 + 2 Я 0 Я 'со 5 Т 1 5 1 п у 1

 

.

,

.

 

 

L e ^ a

(Я0)2+

(tfi)2 cos 2yi + 2Я0Я1 cos fi cosy1

'

^

'

 

 

t a 9 v . .

( Я 1 ) 2 5 1 п 2 7 1 + 2 Я 0 Я ' с о 5 у Ч 1 п Т 1

 

 

 

,

 

 

 

1 ё ^ Т а

( Я « ) 2 +

( Я 1 )2 cos 2^1 + 2 Я 0 Я 1 cos 71 cosyi

 

yu-ti)

Из этих выражений видно, что суммарное поле становится эл­

липтически поляризованным

только в том

случае,

если

одновре­

менно выполняются неравенства У1ФО и ф ^ О .

 

 

 

 

 

 

На практике больший интерес представляют равенства, обрат­

ные

(11.20)—(11.27), они определяют

параметры

вторичного

поля

Я 1 ,

у1 , ф1

через элементы

суммарного

эллиптически

поляризован­

ного поля (первичное поле предполагается известным). Нетрудно показать, что для полной характеристики линейно поляризованного вторичного поля достаточно располагать любой тройкой независи­ мых параметров, характеризующих эллипс поляризации суммарного поля на пространственной или фазовой плоскости [58].

28

Если в качестве таких троек взять, например, Нх, Ну, срж и HRe,

Him, Уяе, то выражения для Я 1 , ср1

и Y 1 могут быть записаны в виде

W=y(H°)2+H2x+H2y-2H°Hxcos<fx

 

;

(11.28)

Н1 = У(Н°У+Н1

 

+//L-2tf7/R e cosT R e ;

(11.28')

to ср' =

 

^ s i n ^

 

(Н 29)

 

 

Я„

Sin 7„

 

 

 

t g T ' =

R e

^ — - ;

 

(11.29')

tg т ' =

 

H

y

;

 

(H.30)

] / ( Я ° ) 2

+ я 2 - 2 Я 0 Я ^ с о з ^

 

 

tgcpi =

 

Я " "

r

- .

(11.30')

Аналогичные выражения

для Я 1 , ср1, у1

могут

быть

получены и

через другие тройки независимых параметров. Такие тройки можно составить также из величин, инвариантных в пространстве ( Я а , Нь, фа ) или на фазовой плоскости ( Я а , Нь, уа)- Однако соответствующие формулы оказываются достаточно сложными. Поскольку и прямые выражения для этого случая (11.24) — (11.27) тоже сложны и нена­ глядны, целесообразно их упростить, рассмотрев практически важ­ ный частный случай малых аномалий, когда выполняется условие

Н°Э>Н\

Ограничиваясь членами первого порядка малости, найдем:

На да Я ° + Я 1 cos ср1

cos т ^ Я 0 +

Re Нхх;

 

(11.31)

tfj^Z/'sliVsliiT^InitfJ;

 

 

 

 

 

(11.32)

 

/ / 1

 

 

^

Im Н\

 

 

(11.33)

n ^ T T j S i n c p ' c o s T 1

я о

;

 

 

 

H l

1

- 1

R e / / y

 

 

/ I I

ЪЛ\

T a ^ 7 7 F C O S ( P

s m T 1 = = я о

 

 

(

И - 3 4 )

Из полученных выражений

следует, что графики распределения

поля в пространстве

для большой полуоси эллипса поляризации и

ее начальной фазы будут сходны с графиками

компоненты вторич­

ного поля, параллельной первичному, а графики

малой

полуоси

и угла наклона — с графиками

ортогональной

к первичному полю

компоненты вторичного поля. Частотные характеристики

аномаль­

ных значений большой полуоси и угла наклона будут

аналогичны

частотной характеристике реальной части вторичного

поля, а ча­

стотные характеристики малой полуоси и начальной фазы

большой

полуоси — частотным

характеристикам

мнимой

части

вторичного

29

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