автореферат диссертации по разработке полезных ископаемых, 05.15.11, диссертация на тему:Разработка методики для выбора параметров крепи вертикальных выработок области многолетней мерзлоты

АКАДЕМИЯ НАУК СССР ОРДЕНА ЛЗГЛНА СИБИРСКОЕ ОТДЕШШЕ

ИНСТИТУТ ГОГхЮГО ДЕЛА СЕВЕРА

На правах рукописи

Иудин .Михаил Михайлович

УДК 622:28-551.34

РАЗРАБОТКА МЕТОДИКИ ДЛЯ ВЫБОРА ПАРАМЕТРОВ КРЕПИ ВЕРТИКАЛЬНЫХ ВЫРАБОТОК ОБЛАСТИ МНОГОЛЕТНЕЙ МЕРЗЛОТЫ

Специальность: 05.15.11 - Физические процессы горного

производства

Автореферат диссертации на соискание ученой степени кандидата технических наук

Якуток 1991

Работа выполнена в -Институте горного дела Севера Сибирского отделения АН СССР

Научный руководитель - доктор технических наук, профессор

Скуба Валентин Николаевич

Официальные оппоненты - доктор технических наук, профессор

Насонов Илья «Дмитриевич

доктор технических наук, профессор Батугин Сергей Андриановпч

Ведущее предприятие - Всесоюзный научно-исследовательский институт золота и редких металлов В1Ш-1

Защита состоится 1991 г. в часов на'

заседании Специализированного совета К.003.44.01 при Институте горного дела Севера СО АН СССР (677891, г.Якутск, ГСП, ул. 50 лет ВЛКСМ, 26

С диссертацией можно ознакомиться в библиотеке Института горного дела Севера СО АН СССР

Автореферат разослан " У " Я 1991 г.

Ученый секретарь Специализированного совета, доктор технических наук, профессор

В.Ю.Изаксон

ОБЩАЯ ХАРАКТЕРИСТИК РАБОТЫ

Актуальность. При проектировании й разработке месторождений полезных ископаемых ва:;гнсе место занимает строительство и эксплуатация вертикальных выработок, выбор парамзтроз кропи которых существенно зависит от горногеологлчсоких и геокриологических условий. Однако,как показывает практика, расчет кропи отзолов осуществляется по методикам и рекомендациям, которуэ не учитывают специфику и особенности деформирования мзрзлых перед. Опыт эксплуатации выработок выявил, что основное влияние на безопасную эксплуатацию ствола оказывают тепловые процессы в приконтурном слое мерзлых пород. Ни в одной чо известных работ влияние тепловых процессов на механизм совместного деформирования крепи и пород не учитывается в полной мере, что приводит в одних случаях к завышенным, в других - к заниженным результатам при определении параметров крепи стволов. Этим объясняются многие'случаи раэруше!шя крепи сырабсток, эксплуатирующихся в многолзтиемерэлых породах.

Работа выполнена в соответствии с планом научных исследований ИГДС СО АН СССР по теме "Разработать научные основы крепления и поддержания капитальных и подготовительных выработок з условиях многолетней мерзлоты" (К> гос .регистрации 8С070053) и по договорной теме "Разработать рекомендации по управлению тепловш/м и механическими процессами при строительстве и эксплуатации алмазных рудников" О.' гос .регистрации 0286.0045157).

Идея работы: Формирование нагрузки на крепь вертикачьного ствола, пройденной в многолетнемерзлнх коренных горных породах, происходит под воздействием механических процессов, протекающих v приконтурном слое массива, и связанных с изменением деформационных характеристик мерзлых пород в зоне прстаивання.

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

Задачи исследований;

- обобщить опыт.эксплуатации и провести анализ существующих

методов расчета крепи вертикальных выработок, пройденных в мно-голетнемерзлкх породах;

- экспериментально исследовать тепловые процессы в мерзлых породах вокруг вертикального ствола под воздействием теплового реиа..а выработки, установить динамику образования зоны оттаивания;

- исследовать влияние горногеологлческих и геокриологических условий на характер взаимодействия крепи вертикального ствола с массивом горных пород;

- разработать метод расчета нагрузки на крепь вертикального ствола, проДленного в мерзлых породах.

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

Научные положения. защищаемые автором:

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

- установлено, что основное прирацспие нагрузки на крепь вертикального ствола происходит при размере зоны противания в 0,1 ... 1,5 радиуса выработки и степени упругой неоднородности пород,- изменяилейся от 1,2 до 1,0.

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

Наушая новизна работы состоит в том, что

- научно обоснована методика расчета нагрузки на крепь вертикальной выработки, учитываяюцая влияние изменения де^рыацион-ных свойств мерзлых пород в зоне протаивания;

, - установлено, что уменьшение модуля упругости мерзлых пород в зоне протаивания является основным фактором в формировании нагрузки на крепь ствола при взаимодействии крепи с оттаи-ващиы массивом мерзлых пород;

- определены условия и количественно оценено влияние горно-

геологических и геокриологических факторов на интенсивность формирования нагрузки на крепь ствола.

Личный вклад автора состоит:

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

- в изучении особенностей взаимодействия крепи вертикального ствола с оттаивающим массивом мерзлых пород;

- в разработке методики расчета крепи верткглльного ствола, пройденной в кноголетнемерзлкх горных породах;

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

Практическая ценность работы состоит в еозг.;ол-лости использования полученных результатов для выбора рацион,гшшх параметров крепи вертикального ствола в'условиях Севера. Результаты исследований могут слузить теоретической основой для разработки нормативной документации и могут быть использованы проектными организациями. . .

Реализация работы. Результаты выполненных автором исследований вошли составной чаотью в работу "Методика определения размеров протаивания вокруг горных выработок ш&хт Северо-Востока СССР при полом5тельном тепловом режиме" (Якутск, 1982 г.). Рекомендации по рациональному регпму растепления ЛПО переданы для-использования институту "Якутнии промалмаз" ШО "Якуталмаз" с ожидаемым экономическим эффектом в 44 тыс.руб. в год на один вертикальный ствол.

Результаты диссертационной работы внедрены в учебном процессе по специальности "Подземная разработка месторовдений полезных ископаемых" в Якутском госуниверситете, а такяе использованы при. составлении учебного пособия.

Апробация работы: основные положения диссертации докладывались на I конференции молодых ученых и специалистов ИГД.С (г.Якутск, 1981 г.), на УП Всесоюзной конференции по механике горных пород (г.Днепропетровск, 1981 г.), УШ Всесоюзном семинаре по исследовании горного давления и. способов охраны капитальных и подготовительных выработок (г.Якутск, 1982 г.), на ХУЛ Всесоюзном семинаре по горной теплофизике (г.Якутск, 1986 г.).

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

Объе:.; рай о 1ч;. диссертационная работа состоит из введения, пяти глав, заключения и содержит 133 страницы машинописного текста, II рисунков, <0 таблиц, списка литература из 90 наименований. '

С<£таШИ2 РАБОТЫ

В перзой главе приведен аначиз опыта крепления вертикальных выработок в ыноголеткемерзлых породах и методов расчета параметров крепи стволов. Проанализировано влияние тепловых процессов, протекающих в прел;: и массиве пород па устойчивость вертикальных стволов. Выявлены основные причины разрушения крепи, обуслоален-ные совместна действие}.: неравномерных нагрузок, вызванных силой тяжести оттаявших мерзлых грунтов, и разрушающих термических напряжений в материале крепи. Отмечено различие в механическом ззашодействии мерзлых грунтов и мерзлых осадочных пород с крепью выработки. Обосновала необходимость дифференцированного -подхода по расчету крепи ствола а зависимости от типа вмещаадкх пород при проталвании мерзлого массива.

Во второй глазе изложены результаты экспериментальных исследований по динамике протаицания мерзлых лород вокруг вертикального ствола.

Исследования проводились на физической модели,'основными элементами которого являются металлически'! стенд, модель ствола, эквивалентны!: материал и контрольно-пзмерптелькач аппаратура. Линейный геометрический масштаб при моделировании . составлял 1450. Диаметр модели.ствола был равен 130 ш. Замораживание,• выстоГжа модели массива пород, поддержание граничных условий осуществлялась холодным воздухом, который■циркулировал в замкнутой системе из теплоизоляционной рубрики и воздуховодов. Из- ' аарение температуры массива пород выполнялось медь-константа-ковшл термопара.-,:::. При моделировании варьировались температура выстойкп подели массива пород и температура воздуха в стволе. •

Регулируемыми параметрами является температурный и времен- ' ной масштабы, определяемые через критерии К0, р0 .

Обработка экспериментальных данных метода:/,и статистики позволила получить уравнения регрессии с коэффициентами корреляции 0,56-0,28, а проверка достоверности корреляционной связи

подтверждает зависимость глубины протаиваяия от времени в л:у;е

$ (1п Ъ) • Бели получены регрессионные зависимости для наиболее характернцх э теплсфизическом отношении пород: известняк, ангидрит, песок, доломит. Совместная обработка этих результатов в критериально;; форме позволила получать следующую сормулу для расчета зоны протапвания вокруг 2ерг:::;а;ън0Х'0 ствола:

аКт [0.53 + ехр{-0.Ц.Ко)}(1.*-0.5-К8) . , ... (I)

Яо'0-57 1 - 0,922 - 8 + 0,(52 - 8 г Ьп^о-^»

где лКт - глубина противалпя, м; Ко- радиус ствола, и; Ко -критерий Коссовича; К{ - критерий Козкера, Р0 - критерий Фурье; 8 - безразмерна температурнп!: симплекс.

В третьей главе рассмотрен процесс формирования нагрузи! на крепь вертикального ствола при оттаивании мерзлых -пород в стадии их совместного деформирования.

В начале главы выполнен расчет напря>г.енно-де$ор:ирозакного состояния оттаивающего мерзлого массива пород вокруг вертикального ствола. В мерзлом массиве пород проводится вертикальная выработка круглого сечения. Образование полости приводит к перераспределению начального напряденного состояния. Постоянная крепь возводится с отставанием от забоя ствола во времени и в пространстве и имёет прочный контакт с породным обнажением. По истечении некоторого времени под воздействием теплого воздуха начинается процесс оттаивания мерзлых пород вокруг закрепленного ствола. Образуется зона талых пород, в которой механические свойства мзрзлцх пород изменились. Новое механическое состояние приконтурного слоя массива приводит к вторичному перераспределения НДС массива вокруг ствола. Это перераспределение приводит к совместному деформирования крепи ствола с массивом пород, что вызывает соответствующее нагрунение крепи выработки.

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

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

Осноыюе уравнение совместности перемещений контура крепи и породного обнажения имеет вид:

где и, - смешение породного контура до ввода постоянной крепи в эксплуатации, иг - смещение внешнего контура крепи к моменту установления равновесия з системе крепь-массив; и5 - смещение дородного контура к моменту установления равновесия в системе крепь-ыасскв.

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

где - коэффициент Пуассона, модуль упругости мерзлых по-

род; С|, - напряжение нетронутого массива пород.

Смещение внешнего контура крепи иг равно:

где ¡>0 . Е0 - козй^ициект Пуассона, модуль крутости материала крепи; г0- радиус ствола в свету; Р - нагрузка на крепь ствола; вызываемая совместным деформированием с протаявшим мерзлым массивом пород.

Смещение породного контура ц5 определяем из решения ЩС оттаявхего мерзлого массива пород п с помощью условия (2) с учётом (3), (4) найдем величину нагрузи: на крепь ствола:

и,» иг«и5

(2)

(3)

(4)

р3 с(°|-аг)

(5)

где

1

ггт

1 М^Ч). Е.

а-.

I

е,-(г* -1) * i: \irit V»

р,, Е, - коэс^фициент Пуассона, модуль упругости талкх: пород;

гт- безразмерный радиус протаиванья.

При слабой зависимости смещения породного контура и3 от величины давления Р нагрузка на крепь' ствола определяется:

р. Ев (|»1>2) («-4)-(гг-|)(п,-У<) \ . Егг* ' * (6)

где п,= £г/Е, , "-Ь,- < < .

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

Е.С.Ержанов показал, что деформирование большинства пород описывается ядром ползучести в виде степенной функции. Путем подстановки временных операторов £ , V в выражения (5), (6) м№.-ло получить зависиности нагрузки на крепь ствола, учитывающие проявления свойств ползучести горными породами.

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

На величину и характер деформационных процессов лрпконтур-ного слоя массива существенно.влияэт геокриологические и горногеологические факторы, х которым откосятся уровень напряженности

Э

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

Выполненные расчеты показывают, что основное прирадение нагрузки на крепь происходит при изменении степени упругой неоднородности от 1,1 до 1,8. С увеличением глубины протаивания влияние на интенсивность нагрузки снижается. Например, при л!?г»|ц величина йР составляет 80,1$, а при дКт = 20 м - уяе 76,6$. Другим существенным фактором является размер зоны протаиЕания. Анализ показывает, что наиболылой прирост нагрузки на крзпь- ствола наблюдается при ширине зоны протаивания, равной 0,1 ... 1,5 радиуса выработки. Известно, что коэффициент Пуассона практически не влияет на формирование нагрузки на крепь. При условии равенства коэффициента Пуассона мерзлых и талых пород расчеты под-творвдаот этот вывод. В случае расхождения этих коэффициентов влияние на изменение интенсивности нагрузки сказывается заметно, что необходимо учитывать при расчете нагрузки на крепь ствола.

Таким образом наиболее существенными факторами, влияющими на формирование нагрузки на крепь, являются уровень напряженности нетронутого массива, степень упругой' неоднородности и размер зоны протаивания. .

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

- расчёт теплового рекиыа ствола, включающий прогноз, температуры воздуха и массива пород;

- установление расчетной схемы, исходных механических параметров проектируемой конструкции крепи и массива горных пород;

- определение сил взаимодействия крепи и массива, внутренних усилий, напряжений и деформаций в крепи;

- проверка несущей способности крепи и подбор конструкции необходимой и достаточной несущей способности для конкретных горнотехнических условий проектирования.

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

Глубину протаивания вокруг ствола на протяненной части вы-

чиоляют по формуле

йЙт» дЯг0 Кг-К5-К«-К,-К1 , (7)

где Kt.Ks.Kw , К, . Кг - эмпирические коэффициенты, опре-

деленные по таблицам 1,2,3.

Нормативная нагрузка на крепь ствола равна:

рн* 2.22 к, -¡и, (а)

где К*, К„ , Кр - эмпирические коэффициенты (табл.4, 5); у -плотность горных пород; Н - глубина ствола.

Расчетная нагрузка я толщина крепи вычисляются по формуле:.! СНиП 11-94.80.

Таблица I

Коэ£- !т7ЯГ,„,рп,п!_Срок службы ствола, год. .. .

фициент! Пар«ие .р,-5-*-10' --37

а К? 2,2 3,0 3,7 4,2

V, 5, £

30 0,92 0,24 0,95 0,56

40 1,03 1,06 1,04 1,03

50 1,22 * 1,16 1,12 1,И

I 0,52 0,53 0,54 0,65

5 1,62 1,49 1,48 1,45

10 2,21 2,1'6 2,12 2,08

* 20 2,52 2,78 2,70 2,65

• Таблица 2

Злакнос ть, V* , % 1 ^ 4 6 8 10

Ку/ 1,38 1,0 0,85 0,79 0,77

Таблица 3

Коэффициент ! I Вт/мК

! 1,0 . 1.5 2,0 2,5 3,0

Ч 1,18 1,09 1,0 0,91 0,82

*2 0,28 0,63 0,82 0,93 1,0

....... ......Таблица 4.....

Коэффициент ! ! ДКг'Авт/йо

! ! 0,2 0.6 1,0 2.0 4,0 6,0

-•к« 0,375 0,714 0,862 1.0 1,071 1,09

Таблица 5

Коэффициент I 1 Коэффициент упругой неоднородности .

! ! 1Д 1,2 1.3 1,4 1,6 1,8

к„ 0,35 0,62 0,33 1,0 1,24 ,1,38..

Примеры расчетов по СНиП 11-94-80 и по методике покаэываш преимущество разработанной методики, которая учитывает влияние таких факторов как степень упругой неоднородности и размеры зоны протаивания на формирование нагрузки на крепь ствола при оттаивании приконтурного слоя массива ыерзлнх пород, что позволяет рационально подобрать параметры крепи вертикального ствола.

В пятой главе проанализировано терыонапряженное состояние многослойной крепи вертикального ствола в замороженных породах и даны рекомендации по рациональному режиму растепления ледопо-родного ограаяения.

Длительность периода оттаивания ЛПО оценивалась по экспериментальной зависимости, полученной на основе результатов физического моделирования, а таюзэ по формуле В.Т.Седова.

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

.та* . .

¿в 4 «пр (9)

где Кор - расчетное сопротивление сжатию материала слоя крепи.

Анализ термокапряхешюго состояния чугунно-бетонной конструкции крепи ствола было выполнено при температуре воздуха в выработке, 20, 40 , 60°С.

• Во всех случаях значителыгае тангенциальные напряжения возникает в слое из чугунных тюбингов, но условие прочности всегда соблюдается. 3 других слоях (бетонном и цементио-песча-ном) конструкции при температурах воздуха 40-60°С уровень тер-ыонапряяений значительно выше несущей способности этих материалов, что приведет к образованию зон трециноватости и разрушений в данных слоях конструкции крепи.

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

На первой .ступени воздух нагревается до 40°С и подается в выработку в течение 11-14 суток, обеспечивая прогрев конструкции крепи до пологаиелькых температур. На второй ступени воздух подогревается до 60°С с длительностью подачи 11-21 суток. Расчетная продслгяителыгасть оттаивания .НПО вокруг вертикального ствола составит 25-35 суток.

ЗАКЛЮЧЕНИЕ

В работе дано новое решение актуальной каучной задачи выбора равдональных параметров крепи вертикальных выработок, проводимых в коренных многолетнемератах породах.

1. Суцествувдие представлений о механизме формирования нагрузки на крепь вертикального ствола не учитывает1 особенностей взаимодействия элементов системы "крепь - зона протаивания мерзлых пород". .

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

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

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

5. Исследовано терыонапрягенное состояние многослойной крепи ствола при оттаивании ледопородного ограждения. Определена продолжительность оттаивания ЛПО при различных режимах нагрева вентиляционного воздуха. Исходя из прочности крепи ствола выбрал оптимальный ступенчато возрастающий режим оттаивания ЛПО. Ожидаемый экономический эффект от внедрения рекомендаций по оттаиванию ЛПО составит 44 тыс. рублей.

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

7. Разработана методика расчета параметров крепи вертикального ствола в многолетнемерзлых породах при положительном тепловом режиме.

Перечень работ, в которых опубликованы научные результаты, включаемые в диссертации:

1. Изаксон В.10.', Сенук Д.П., Захамалдш А.И., Иудин UM. Напряженное состояние системы "порода-крепь" вертикальных стволов зоны многолетней мерзлоты при пологдтельном тепловом реглкме // Тезисы докладов Всесоюзной конференции по механнко гори« пород. - М. : ИГД км.А.А.Скочинского, IS8I г.

2. CavoxnH A.B., Иудин LI..VI., Шурпш Б.В. Установка для физического моделирования теплового решила в мерзлых породах вокруг горных выработок // Проблемы горного дела Севера. -ЯкутскСО АН СССР, 1981 г.

3. Методика определения размеров протаизанкя вокруг горных выработок иахт Северо-Востока СССР при положительно:.", тепловом режиме (Иудин U.M. раздел 6. Расчет глубины протаивания вокруг вертикальных выработок).- Якутск: ЯЗ СО АН СССР, 1982 г.

' 4. Иудин М.М. Методика расчета крепи капитальных выработок области многолетней мерзлоты // Крепление, поддержание . и охрана горных выработок,- Новосибирск: КГЦ. СО Ail СССР, 1983 г.

5. Иудин М.М. Расчет нагрузок на крепь ствола при оттаивании мерзлых 'пород // Совершенствование, технологии подземной разработки россыпных и рудных месторождений севера.- Якутск: ЯФ СО АН СССР, 1987 г.

6. Иудин U.M. Термонапряженное состояние многослойных ¡конструкций крепи вертикального ствола при размораживании ЛПО // Повыззние эффективности освоения рудных месторождений Севера.- Якутск: ЯНЦ СО АН СССР, 1989 г.