Добавил:

ivanov666 Опубликованный материал нарушает ваши авторские права? Сообщите нам.

Вуз:

Башкирский Государственный Аграрный Университет

Предмет:

[НЕСОРТИРОВАННОЕ]

Файл:

книги из ГПНТБ / Сохранов Н.Н. Машинные методы обработки и интерпретации результатов геофизических исследований скважин

.pdf

Скачиваний:

Добавлен:

25.10.2023

Размер:

14.73 Mб

Скачать

☆

<<< < Предыдущая 1 2 3 4 5 6 7 8 9 1011 / 2411 12 13 14 15 16 17 18 19 20 21 22 23 24 > Следующая >>>

5. При выполнении условия (28) определяют глубину z'u метки, ограничивающей интервал, сумма слов в котором удовлетворяет требованиям синхронизации преобразования, и записывают сумму

s
2 NM,	глубину z'u и число Вм	кратности суммы слов номинальному
значению Лги. В случае невыполнения условия				(28)—(29)	к про
анализированной сумме добавляют следующее значение NN					и повто
ряют	анализ.
6. Определение глубины zx		первого слова записи каротажных
данных для первого цикла анализа меток (і'			=	1) по формуле
	2Х =	2 „ + А ^ „1			(30)
и проверка соответствия ее глубине забоя (zx			z3 a 6	и zx — AzM	z3 a C ).

7.Запись данных в таблицу ТМ при совпадении глубин с точ ностью до 2—3 м.

8.Вычисление глубины zM первой метки и запись данных в таб лицу ТМ.

9. При невыполнении		обоих условий		соответствия					глубин zx
и z3 a 6 по пункту 7 останов		ЭВМ и запись			штрафного				очка в таб
лицу ТК.				(і =		іы)
10. После окончания анализа меток							из	таблицы ТМ
вызывают первую пару значений 2 NM			и		и	определяют число a't
избыточных (al > 0)	или недостающих		слов.
11. Перезапись цифровых данных М\,				соответствующих						интер
валу с порядковым номером і'', с пропуском а\							слов	(а[ >		0) или
добавлением а\ слов	(а. <	0). Исключение			слов и добавление их

производят через равные участки интервала. Перезапись производят

поочередно		во всех интервалах		с а\ =±= 0				до выполнения условия
Ï — і'ы, где ï — число записей				2	-^м в	таблице ТМ.
Рассмотрим алгоритм исправления ошибок и увязки данных по
глубине с учетом смещения точек записи							(рис. 39).
1. Распаковка слов массива Mj цифровых данных, представля
ющих кривую геофизической величины с шифром /'.
2.	Нахождение		экстремумов		кривой		и	запись	их координат
(^ЭКС. І1	^экс. I
		) в таблицу Т (аэ к с ).							средних с одно
3. Сглаживание кривой по способу скользящих
временным		нахождением экстремумов				сглаженной			кривой а'экс. ѵ-
Экстремумы		аэкс. с	записывают	в	таблицу		Т (аэкс_ г)- Сглаживание
кривой производят вычислением средних значений									а\ для m точек
кривой по		формуле			і Уп

			a ''= i!r2a "						( 3 1 )

смещая каждый раз интервал вычисления на шаг дискретизации А. Значение m целесообразно брать равным трем. Сглаженная кривая 100существенно меньше искажена помехами измерения, и поэтому

Определение адреса h; (32) записи кривой

при большом влиянии помех может


быть	использована			для	интерпрета
ции	вместо		исходных		кривых.
4.	Поочередное			сравнение орди
нат сглаженной и исходной						кривых,
взятых		соответственно			в	точках ï
и 1 =		1 + - ^ .
5.	При		большом		расхождении
ординат		кривой результат				измере

ния геофизической величины в точ ке і кривой ah соответствующей максимуму расхождения Да,- = | at —

— а\>\ — max, заменяют средним значением измерений щ в соседних

точках,	если	выполняются			следу
ющие условия:
а) Д а т а х ï >		га[, при		совпадении
глубин	экстремумов		обоих		кривых
(Д2Э	ЭКС. (		Р И	•1 = 0);
б) Д а т а х . і ' > е 2 а і П				небольшом
расхождении		глубин		экстремумов
кривых	а'ѵ и а г ( Д г э к с . г			=	\|z3 K C ., —

ïm — 1 \

—^экс.Г I ^ — — ) •

6.Указанные операции повто ряют до тех пор, пока не будут обработаны все точки кривой / (см.

рис. 39; ï = N„ — m +і,			где N„ —
полное	число	слов в кривой).
7.	После	исправления	ошибок

в записи цифровых данных опреде

ляют		начальный адрес		Aj разме
щения		этих	данных в МОЗУ		или
на	магнитном		барабане	(магнитной
ленте).
	8.	Преобразование		данных	по
формуле
		а,-= па£ А(аі ),
где	aL — геофизическая			величина,

выраженная в физических единицах;

к (а£) — поправочный			коэффициент
на нелинейность		измерений; п —
масштаб записи.
9.	Запись	отредактированных
данных	на временное		хранение.

Вход

Нахождение О •

ïI і+т

Нахождение Q3KCjJ

а--=па1К(ас)

Запись [Qj)j по Aj

Рис. 39. Блок-схема алгоритма испра вления ошибок и увязки данных по глубине.

101

Из рассмотрения описанных алгоритмов можно рекомендовать

следующую		схему	редактирования цифровых	каротажных данных
на ЭВМ.
1.	Ввод	данных	в ЭВМ.
2.	Обработка контрольных записей нуль-			и стандарт-сигналов

с целью уточнения масштаба записи (цены разряда цифровых дан

ных) и положения нуля.
3.	Увязка меток глубины и числа слов (циклов) между ними.
4.	Распаковка слов, исправление ошибок в записи отдельных чисел
(сбоев	преобразования, записи и ввода) и увязка данных по глубине
с учетом разной глубины точек		записи геофизических величин.
По	окончании редактирования	числовых массивов M j геофизиче

ских данных на печать выдаются результаты оценки качества их с указанием характера обнаруженных погрешностей (сбоев) и ча стота их встречи. Можно рекомендовать следующий образец таблицы

для	печати		(табл.			5).					Т а б л и ц а 5
											Т а б л и ц а 5
		Оценка		качества отредактированных геофизических данных
		Характер погрешности (штрафа)							Число	Процент	Оценка (балл)
		Характер погрешности (штрафа)							ошибок	ошибок	Оценка (балл)
Ошибка		установки			масштаба				0	—	3
									1		2
									2		1
Расхождение			числа		слов в интервале			меток,	0		3
									1		2
										Sa 20	1
Сбои	сннхронизащтп					при записи слова		(цик	0	3=0,2	3
ла)										3=0,2	2
											1
Ошибка		измерения				(записи)	геофизической		0		3
величины										3=0,2	2
										> 1	1
Несогласование				измерений с			контрольными		0		3
данными									1		2
										> 2 0	1
Общая		оценка		качества					—	—	Средний
											балл

14. ТРАНСФОРМАЦИЯ ГЕОФИЗИЧЕСКИХ Д А Н Н Ы Х

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

102

их за влияние вмещающих пород, скважины и аппаратурных помех. Это влияние обычно выражается в искажении формы кривых и изменении существенных значений, отсчитываемых против пласта.

Поэтому

трансформация

заключается

в следующем:

а)

выделение всех прослоев по экстремумам или точкам перегиба;

б) отсчет существенных значений против

прослоев;

в)

исправление

существенных

значений;

г) присвоение исправленных значений всем точкам в интервале

прослоев и запись трансформированных

данных по точкам в форме,

принятой для исходных геофизических дан

ных:

l j 11' a

12'

а/1

211

u Ci;

., a

12m2i

j'lmn

22i

«/Vil

•

(33)

где

/ ll

— Я/12'

a/i m ,

— исправлен-

ное существенное значение для первого

прослоя;

а/ Ѵ і

= ауѵ „

аjvmv

ТО ЖѲ,

для ѵ-го прослоя;

— номер

параметра;

ту — число

точек

в ѵ-м прослое.

Возможна

более

плотная

(экономичная)

форма

записи

трансформированной

кривой:

Й Д> тѵ> а /2і т2І

;

ßjy тпѵі • • •

(34)

Для выделения прослоев необходимо при

менять способы отбивки

границ,

автомати

чески

обеспечивающие

фильтрацию

помех

измерения

регистрации

(например,

из-

зубренность кривой). Таким способом яв

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

гра

Рис. 40.

Примеры

отсчета

ниц.

средних

значений о

кривых

прослоев

зависимости

от

бокового

(а) и

радиоактив

Против

ного (б) каротажа.

вида каротажа и мощности прослоя

(пла

ста)

отсчитывают

средние

значения

а с р ,

экстремальные

(макси

мальные

а т а х ,

минимальные

а т і

п ) .

Средние значения

отсчиты

вают в прослое, если мощность его больше граничного значения НТ.

Величину		Нг устанавливают для каждого вида каротажа.				Среднее
показание		обычно	вычисляют	не во всем	интервале	прослоя
Z K — z n	( Z K '	z n — глубины соответственно кровли			и подошвы), а в той
части	его	z'K—z'n,	где влияние	вмещающих	пород незначительно

(с учетом зоны экрана). Примеры отсчета средних значений для кри вых Б К и НГК показаны на рис. 40. Экстремальные значения от считывают при мощности прослоя меньше граничного значения.

Существенные значения исправляют на влияние следующих фак торов: мощности прослоя, скважины (КС, ГК и НГК), удельного сопротивления прослоя и параметров зоны проникновения (ПС), помехи интегрирующей ячейки (НГК, ГК и др.) и т. п. Влияние мощности интегрирующей ячейки учитывают только в тонких

103

прослоях (# - < Нт). Исправление существенных значений целесо образно производить умножением их на поправочные коэффициенты Kj, определяемые по кривым зависимости Kj — ф (Сх, С2), где Сг и С2 — величины, влияющие на результаты геофизических измерений.

По исправленным существенным значениям (а'ср, a'max и a'min) определяют приведенные' значения геофизических величин по фор муле

	—	I а'іч —a / l		I
	°/v = п. -a-		'	( 3 5 )
	— "vi
где	a;-v — исправленное существенное	значение		против ѵ-го пласта;
ajx	и ay, — исправленные существенные значения против опорных

пластов, например против глин (ау 1 ) и плотных (или чистых) пород (а/г)-

Приведенные существенные значения присваивают всем точкам в интервале пласта и записывают эти значения по форме (33) или (34).

15. У В Я З К А Д А Н Н Ы Х КАРОТАЖА ПО Г Л У Б И Н Е

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

1. Увязка каротажных кривых взаимной корреляцией. Способ

взаимной корреляции кривых	[46] заключается		в вычислении		функ
ции взаимной корреляции	k+n
	k+n
R (ô)0 , -	—	(	3	6	)
4 ' 0 l	f h+n	k+n			4 '

'i=/t i=k

для двух последовательностей чисел: аоі (основной кривой) и а{і (согласуемой кривой). Одна из них сдвигается относительно другой каждый раз на величину шага дискретизации А по глубине до полу чения максимума функции взаимной корреляции. Число сдвигов ô = ±1 одной последовательности чисел относительно другой до положения максимума функции взаимной корреляции соответствует величине смещения Az(- = Al кривых по глубине. Эта величина учи тывается при интерпретации. Интервал корреляции (количество чисел ІѴК) обычно выбирается в зависимости от характера коррели руемых кривых. Для Туймазинского месторождения этот интервал достаточно взять равным 20—25 м (7ѴК = 100—125).

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

бован на фактическом материале. Результаты опробования в скв. 718

104

ПІкапово (см. рис. 33) при А = 0,2 м приведены в табл. 6. Как видно, величины смещения кривых разнородных геофизических величин определяются методом взаимной корреляции со значительными ошибками (до 2 м). Из этого следует, что способ взаимной корреля ции пригоден для увязки по глубине кривых только однородных геофизических величин (ПС, КС, записанных градиент-зондами с АО < Я і др.) и не годится для увязки по глубине в общем случае.

					Т а б л и ц а 6
Результаты определения смещения кривых различных видов каротажа
	относительно кривой КС А0,8М0,Ш		методом взаимной корреляции
		(Д = 0,2	ы)
		Смещение в величинах Д			Ошибка в величинах
			определен!тое на ЭВМ
	Вид измерения	фактиче	без транс-	с транс	без транс	с транс
		ское	без транс-	с транс	формации	форма
			формацші	форма		цией
			формацші	цией
				цией
К С ,	з о н д A4M0.5N	6	5	4	1	2
К С ,	з о н д В7,5А0,75М	6	10	3	4	3
КС,	з о н д АБК - 3	—4	- 1	- 5	- 3	—1
К С, м н к р о з о и д		2	8	3	6	1
ПС		0	10	2	10	2
Н Г К		0	7	0	7	0
КС,	з о н д А0,4М0,Ш	2	2		0

2. Увязка по глубине трансформированных кривых. Результаты определения величины смещения кривых при помощи

взаимной корреляции улучшаются, если кривые предварительно трансформировать в графики прямоугольной формы (см. рис. 33). Погрешность определения смещения (увязки кривых по глубине) в этом случае значительно меньше (см. табл. 6). Это показывает, что увязку геофизических данных по глубине на ЭВМ можно осуще ствлять с небольшой погрешностью (не более 2А) путем корреляции трансформированных кривых по формуле (36) или сопоставлением последовательности глубин z0 ѵ и zjv границ пластов, выделенных соответственно по основной и привязываемым к ней кривым. Для не которых каротажных кривых можно сопоставлять глубины харак терных точек (экстремумов, перегибов и др.).

3. Корреляция границ пластов.		глубин z o v и zjv
Алгоритм корреляции последовательностей					гра
ниц пластов заключается в следующем:
а) определяют положения границ пластов		(характерных		точек)
на коррелируемых кривых и	запоминают их	глубины	z0 ѵ	и	Zjv;
б) для каждой границы z0 ѵ ,	выделенной ка основной		кривой, на

ходят идентичную ей границу z;-v- на другой кривой, расположенной

105

к ней ближе других в		интервале z o V		± Az. Величина Az задается
в зависимости от предполагаемого смещения кривой;
в) вычисляют	абсолютные		разности	глубин \| z o v —	Z/v'\| для	всех
идентичных пар	границ	и их	сумму
	к
•s(ô) = 2l2ov - (z/ v + A-o)\|,				ô = o , i , . . . ,	I;	(37)
здесь к — число пар		границ.

Сумму вычисляют после каждого сдвига границ z;-ѵ< на величину шага квантования до получения минимума суммы S (ô) при б = ± Z. Величина сдвига I последовательности глубин границ zjv относи тельно последовательности глубин границ z0 ѵ основной кривой определяет величину смещения кривых по глубине.

4. Увязка каротажных данных по опорной кривой.

Увязка цифровых данных по опорным кривым целесообразна при многоканальных измерениях в скважине. В этом случае для увязки кривых по глубине необходимо в каждом многоканальном цикле измерений записывать одну и ту же (опорную) кривую. Взаим ная корреляция опорных кривых, записанных в разных циклах измерений, позволяет увязать между собой данные многоканальных измерений [18]. В песчано-глпнистом разрезе опорной кривой может служить кривая ПС. Применение этого простого и надежного спо соба увязки разновременных измерений по глубине в настоящее время встречает большие затруднения в связи с тем, что при проектиро вании комплексных приборов в большей части случаев запись опор ной кривой не предусматривается.

16. ОЦЕНКА КАЧЕСТВА Д А Н Н Ы Х Б К З И ДРУГИХ ЗОНДОВ КАРОТАЖА СОПРОТИВЛЕНИЙ

Точность результатов интерпретации материалов бокового каро тажного зондирования (БКЗ) в значительной степени определяется качеством исходных данных. Перед обработкой материалов Б К З необходимо проверить их качество: установить погрешность опреде ления удельного сопротивления бурового раствора и в случае боль шого ее значения уточнить это сопротивление, выявить кривые со противления с низким качеством измерений (неудовлетворительные замеры) и, если невозможно их исправить, то исключить из обработки.

При определении на ЭВМ удельного сопротивления пластов по данным Б К З и комплекса измерений разнотипными зондами проверка качества исходных данных имеет еще большее значение, чем при ручной обработке. Способы контроля качества в обоих случаях могут быть одними и теми же.

Основными источниками погрешности измерений отдельными зондами являются: неточность установки масштаба регистрации со противления (определение цены разряда кода чисел при цифровой записи), утечки тока, индуктивные помехи и дрейф нуля измеритель

на

ных и регистрирующих систем. На результаты измерений зондами Б К З против пластов высокого сопротивления основное влияние ока зывают ошибки в установке масштаба регистрации сопротивления. Эти ошибки встречаются нередко.

Утечки тока, индуктивные помехи и дрейф нуля, вызывающие смещение точки отсчета показаний, оказывают влияние преимуще ственно на результаты измерения кажущегося удельного сопротивле ния иизкоомиых пластов (за исключением индукционного каротажа). Величину этих помех устанавливают обычно по размеру в колонне и сопоставлением показаний разных зондов против опорных пластов низкого сопротивления. Часто в разрезе нет опорных пластов низкого сопротивления, и замер КС в колонне остается единственным сред ством для оценки влияния указанных помех.

Для исключения	ошибок, вызванных смещением точки отсчета,
из каждого элемента	массива чисел, представляющих кривую КС,
зарегистрированную	зондом с порядковым номером /, вычитают

величину a.jK, соответствующую замеру КС этим зондом в колонне (обычные зонды). В индукционном каротаже смещение точки отсчета исправляют тем, что из кажущейся проводимости вычитают показа ния, зарегистрированные против пластов очень большого удельного сопротивления.

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

1.Наилучшим способом контроля правильности удельного со противления бурового раствора и выявления погрешности, связан ной с неточной установкой масштаба регистрации, является сопо ставление значений кажущихся удельных сопротивлений, получен ных различными зондами [17] против достаточно мощных (Н > 5 м ) высокоомных однородных опорных пластов без зоны проникно вения.

2.Критериями качества измерений различными зондами каро тажа сопротивлений и точности определения удельного сопротивле

ния бурового раствора р с	является наличие и характер	пересечения
интервалов [А/, В,] значений модулей условного р у /р е		сопротивле
ния, границы которых А,	и Ву- получены для различных зондов /

по палеткам для пластов неограниченной мощности без проникнове ния при заданном относительном отклонении ô существенных зна чений КС от истинных (см. главу I X , раздел 23).

При выборе высокоомных опорных пластов признаком двуслой ной кривой зондирования может служить пересечение интервалов [Aj, В j] для всех зондов или для всех зондов, кроме одного. При этом исключают из рассмотрения зонды с самым большим ради усом исследования, показания которых могут отличаться от зна чения, соответствующего пласту неограниченной мощности, из-за неоднородности среды и влияния соседних пластов. При оценке качества и уточнении удельного сопротивления бурового раствора

107

предпочтение должно отдаваться опорным пластам с пересечением всех интервалов.

3. Погрешность определения удельного сопротивления бурового раствора устанавливают по величине интервала высокоомных опор ных пластов для ряда задаваемых значений (приближений) сопроти вления раствора, полученных из исходного умножением [или делением его на некоторый коэффициент Сх (например, Сх = 1,1). В случае большой погрешности исходных значений сопротивления вели чина его корректируется по величине интервала пересечения. Дис персия Ос значения рс , определенного по опорным пластам, не должна превышать некоторого предельного значения, например 15%.

4. В качестве критерия качества измерений зондами каротажа сопротивлений можно использовать наличие пересечения интервалов при различной величине относительного отклонения существенных значений КС от истинных. Пересечение интервалов [Aj, В;-] для всех зондов при е < 10% указывает на хорошее качество измерений по пласту; пересечение при е = 2 0 % — на удовлетворительное каче ство измерений по пласту; отсутствие пересечения — на плохое качество исходных данных по одному из зондов. Непересечение может быть вызвано как погрешностями измерения, так и погреш ностью отсчета существенных значений кажущегося удельного со противления. Качество исходных данных по каждому пласту оце нивают баллом Б , соответственно равным 2, 1 и 0 при хорошем, удовлетворительном и плохом качестве. Качество измерений в целом характеризуют средним баллом Б с р , подсчитанным по всем опорным пластам.

На основе изложенных положений составлен алгоритм контроля качества данных каротажа сопротивлений [17], блок-схема которого показана на рис. 41. По этому алгоритму составлена программа контроля качества данных каротажа сопротивлений на машине Минск-22.

Предполагается, что контроль качества данных каротажа сопро тивлений будет проводиться после того, как отмечены границы пла стов, отсчитаны существенные значения сопротивлений [46] и по результатам измерений малыми градиент-зондами определено удель ное сопротивление р с 0 ) (верхний индекс рс — порядок приближения) бурового раствора. В соответствии с этим исходными данными для контроля качества данных каротажа сопротивлений будут являться

номер пласта	ѵ, глубина его	подошвы	zn ѵ , мощность Нѵ, диаметр
dcv скважины	против пласта,	удельное	сопротивление рс % заполня

ющего ее бурового раствора, исправленные на влияние мощности пла

ста существенные значения		кажущихся удельных	сопротивлений
Pjv (/ — номер	зонда).
Алгоритм	разделяется на	три блока — проверки	правильности

определения удельного сопротивления раствора по отдельным пла стам (I), оценки качества измерений (II) и уточнения удельного сопротивления раствора и определения поправочных коэффициентов для показаний отдельных зондов (III) .

108

Б л о к I . Начиная с первого пласта (ѵ = 1 ) , последовательпо пе ребирают все пласты, и из них выбирают опорные. Опорными счи

таются пласты, удовлетворяющие						следующим условиям: Нѵ^>5			м;
dv =S l,05d„		(d„ — номинальный				диаметр	скважины);	кажущееся
удельное сопротивление, полученное зондом со вторым по величине
радиусом исследования				(/ =	2; обычно это		градиент-зонд с АО		=
= 1,05	м),	больше	15	р с " ;	число	?г3 градиент-зондов с Lj sc 0,85Я
(к ним добавляется и зонд индукционного							каротажа,	если рк ѵ	<5
<: 100	Ом-м) равно трем или больше. Для каждого из опорных пла
стов выполняются			следующие операции.
1. Для градиент-зондов с Lj						=s 0,85ff	и зонда индукционного
каротажа определяют границы А						и Вj интервалов значений модулей

условных сопротивлений р у у / р с , подсчитанных по двуслойной палетке

при

Б =

20%

удельном

сопротивлении

бурового

раствора

р с

= р о

0 ) , и

находят интервал

их

пересечения

[max i l , min В]

ве

личину рш

(min В

— max А)ш.

Указанные

операции проводят

по универсальной программе определения

удельного

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

пластов

[16].

2. Повторяют указанные в предыдущем пункте операции в пред

положении,

что

удельное

сопротивление

бурового

раствора

равно

p c n =

CjPo0' (эта

операция

условно характеризуется

значением кри

терия

0).

При

этом

получается

рш

(min В — max

А)ш

(первое приближение,

s =

1).

Сравнивают полученные значения pw

и р(1).

Если выполняется

условие

р ш > Р ( в \

(38)

то процесс повторяется со значениями удельного сопротивления буро

вого	раствора	рс 2 > =	СіРс0 ';	рс 3 ) =	СІРс"	и т. д. до тех пор,			пока
				p(s+1>>p(s\					(39)
Здесь s — порядковый номер операции и соответственно								показа
тель	степени	при коэффициенте			Сѵ
Значения p c s ) и ри),			при которых условие (39)				перестает	выпол
няться, запоминают.						т. е. ра)	< р(0>,
4.	Если условие (38) не			выполняется,		т. е. ра)	< р(0>,	то	про

водят указанные в первом пункте операции в предположении, что удельное сопротивление бурового раствора

Рс" = Рс°7С1; Рс2 , -Рс°ѴС! и т. д.

до тех пор, пока соблюдается соотношение (39). Эта операция условно

характеризуется		значением	критерия	R	ф 0.
Значения	Pc" ц и рЩ* (и. — номер			опорного пласта), при которых
неравенство	(39)	перестает	выполняться,		запоминают и передают
в блок I I для дальнейшего			анализа.

109

<<< < Предыдущая 1 2 3 4 5 6 7 8 9 1011 / 2411 12 13 14 15 16 17 18 19 20 21 22 23 24 > Следующая >>>

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