автореферат диссертации по разработке полезных ископаемых, 05.15.11, диссертация на тему:Разработка способа повышения эффективности дегазации шахт скважинами, пробуренными с поверхности (на примере Шахтерско-Торезского района Донбасса)

ВВЕДЕНИЕ

I.АНАЛИЗ СОСТОЯНИЯ ВОПРОСА И ЗАДАЧИ

ИССЛЕДОВАНИЯ

I.I.Параметры и режимы работы дегазационных скважин, пробуренных с поверхности

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

Задачи исследования

2.ОПИСАНИЕ МОДЕЛИ ПРОЦЕССА ДЕГАЗАЦИИ ПОДРАБОТАННОГО ГОРНОГО МАССИВА СКВАЖИНАМИ, ПРОБУРЕННЫМИ С

ПОВЕРХНОСТИ

2.1.Постановка задачи моделирования

2.2.Описание принципов организации пакета программ для расчета движения газа в подработанном горном массиве методом конечных элементов.

Выводы.

3.ЭКСПЕРИМЕНТАЛЬНЫЕ ИССЛЕДОВАНИЯ ГАЗОДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК ДЕГАЗИРУЕМОГО ГОРНОГО МАССИВА И ИДЕНТИФИКАЦИЯ МОДЕЛИ

3.1.Характеристика объекта исследований и методика экспериментов

3.2.Результаты экспериментов и идентификация модели.

Выводы.

ИССЛЕДОВАНИЕ И ВЫБОР ОПТИМАЛЬНЫХ ПАРАМЕТРОВ

СКВАЖИН,ПРОБУРЕННЫХ С ПОВЕРХНОСТИ

4.1.Расчет движения газа в подработанном горном массиве без учета дегазационных скважин

2.Расчет движения газа к скважине, пробуренной с поверхности.

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

Выв оды.

Основными направлениями экономического и социального развития СССР на I98I-I985 годы и на период до 1990 года"предусматривается увеличение добычи угля к 1985 году до 770-800 млн.тонн в год. Основной прирост добычи должен быть достигнут за счет повышения производительности труда. Выполнение этого задания Коммунистической партии Советского Союза требует широкого применения новой техники и прогрессивной технологии разработки угольных месторождений.

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

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

- улучшение поддержания подготовительных выработок (затраты на поддержание уменьшаются в 2-4 раза);

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

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

Переход шахт на столбовые системы разработки является главным направлением совершенствования технологии подземной добычи угля. В 1982г. только шахты Минуглепрома УССР применяли их в 853 лавах. При этих системах обеспечивается возможность увеличения нагрузки на очистные забои, оснащенные современными выемочными комплексами, до 5000 т/сут. Однако использование этих возможноетей в газовых шахтах ограничивается обильным выделением метана, сопровождающим горные работы. Так при газообильности шахт 30 и3/т нагрузка на очистной забой ограничена 400-600 т/сут. В Донбассе около 70 шахт имеют газообидьность более 30 м3/т, из них 10 шахт - более 100 цв/т. Количество их непрерывно возрастает. Если в 1960г. сверхкатегорные и опасные по внезапным выбросам угня и газа шахты составляли 45% от общего их количества, то в 1980г. - 60%, а к 2005 году по прогнозу ожидается увеличение количества таких шахт до 80%. Поэтому достижение запланированной добычи угля и повышение производительности труда невозможно без эффективной дегазации угольных месторождений. Широко применяемый для этого способ каптирования метана скважинами, пробуренными из горных выработок, не обеспечивает необходимого снижения метано-обильности. В условиях погашения вентиляционных выработок скважины работают в течение нескольких дней, пока они находятся впереди очистного забоя. После подработки устьев скважин содержание метана в отсасываемом газе быстро падает и скважины отключают. В связи с этим при столбовых системах разработки одновременно работают не более двух скважин и эффективность их не превышает 30%. Такой эффект недостаточен для предотвращения опасных скоплений метана на сопряжениях лав с вентиляционными выработками. Ликвидация их путем усиления проветривания в большинстве случаев невозможна, так как достигнут предел скорости воздуха в рабочих пространствах лав. Там же, где скорость ниже предельной, увеличение расхода воздуха требует больших затрат на реконструкцию вентиляционных систем. Поэтому наиболее реальным способом обеспечения нормального содержания метана в атмосфере горных выработок высокопроизводительных выемочных участков является эффективная дегазация сближенных пластов и выработанных пространств. Попытки обеспечить ее за счет изменения параметров и схем расположения скважин, включая бурение их из специальных выработок, не дали положительных результатов. Этим, по-видимому, объясняется все более широкое применение дегазации подработанных пластов и выработанных пространств скважинами, пробуренными с поверхности, В 1982г. в пяти производственных объединениях были пробурены 105 таких скважин общей протяженностью 38,7 км. Затраты на их бурение и оборудование составили около двух миллионов рублей. В то же время эффективность таких скважин изучена совершенно недостаточно. До настоящего времени нет научно обоснованного метода проектирования дегазации шахт с помощью вертикальных скважин. "Руководством по дегазации угольных шахт" рекомендованы конструкция скважин и параметры способа, основанные на обобщении опыта применения его в Карагандинском бассейне. При этом не учитываются такие важные факторы, влияющие на метанодобываемость скважин, как газоносность дегазируемого массива, расположение в нем угольных пластов, разрежение в скважинах. Пренебрежение ими приводит к существенным ошибкам и нерациональному расходованию средств.

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

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

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

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

Основной объем исследований выполнен в условиях Шахтерско-Торезского района Донбасса. Экспериментальная часть работы выполнена автором совместно с работниками МакНИИ в рамках работы № 0II5805000 "Разработать и внедрить метод комплексной дегазации высокопроизводительных выемочных участков", входящей в программу ПШТ СМ СССР по проблеме 0.05.07 и выполнявшейся МакНИИ совместно с ИПКОН АН. СССР по договору о социалистическом содружестве. Цель работы - совершенствование метода расчета и определение оптимальных параметров и режимов дегазации шахт скважинами, пробуренными с поверхности.

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

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

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

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

Донбасса.

Научная новизна работы заключается в следующем:

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

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

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

Достоверность научных положений и полученных результатов подтверждается:

- удовлетворительной сходимостью (в пределах 10-15 %) расчетных характеристик дегазационной системы с экспериментальными, замеренными в условиях ш."Стожковская" ПО "Шахтерскантрацит";

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

Практическая ценность работы заключается:

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

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

I.Оптимальные параметры скважин, определенные на основании расчетов, применены шахтами "Винницкая" и "Стожковская" ПО "Шахтерскантрацит" с положительным техническим эффектом: на шахте "Винницкая" эффективность дегазации увеличилась с 37 до 56$, на шахте "Стожковская" - с 50 до 71$. В обоих случаях содержание метана в извлекаемой газовой смеси возросло с 25-30$ до 40-50$.

2.Разработанный метод расчета параметров дегазационных скважин, пробуренных с поверхности, вошел составной частью в научно-исследовательскую работу МакНИИ № II5805000079 "Разработать и внедрить метод комплексной дегазации высокопроизводительных выемочных участков", выполняемую по программе ГКНТ при СМ СССР по проблеме 0,05.07. Апробация работы. Основные положения диссертации докладывались:

- на Международной конференции "Применение. ЭВМ и математических методов в горном деле", АПК0М-80 (1980 г.);

- на научно-технических конференциях молодых ученых ИПКОН АН СССР (1981, 1982 гг.).

Публикации. Основные положения диссертации отражены в шести публикациях.

Объем работы. Диссертация состоит из введения, четырех глав и заключения, изложенных на 96 страницах, включая 21 рисунок, 7 таблиц, список литературы из 65 названий и два приложения.

Выводы

I.Разработанный пакет программ применен для расчета реальной дегазационной системы в условиях шахты "Стожковская"? Сравнение результатов расчета с натурными данными показало совпадение; с ошибкой,не превышающей погрешностей определения исходных данных, что доказывает применимость пакета для расчетов оптимальных параметров дегазационных систем. 2.Результаты численного моделирования процесса дегазации

ЗАКЛЮЧЕНИЕ

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

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

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

1.Разработана программная модель движения газа в подработанном горном массиве, реализованная в виде пакета прикладных программ для ЕС ЭВМ, позволяющая исследовать процесс дегазации шахт скважинами, пробуренными с поверхности, и определять оптимальные параметры и режимы работы дегазационных систем на основе реше

••ния задачи фильтрации газа методом конечных элементов.

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

- расстояние между забоем скважины и разрабатываемым пластом - 20 -f 30 его мощностей;

- расстояние между скважинами - 140 - 150 м;

- разрежение в устье скважины - не менее 150 гПа.

3.Выбор параметров дегазации на основе расчетов по предложенной программе, по сравнению с нормативными дает следующие преимущества:

- эффективность дегазации увеличивается, в среднем, на 20%;

- объем буровых работ сокращается в два раза;

- увеличивается содержание метана в извлекаемой газовой смеси и количество метана, пригодного для утилизации.

4.Расчетный экономический эффект от оптимизации параметров на основе математического моделирования составляет, в условиях Шахтероко-Торезского района Донбасса, 66 тыс.руб. в год на одну лаву.

5.Разработан в соавторстве способ дегазации подработанных сближенных пластов вертикольно-горизонтальными скважинами, защищенный авторским свидетельством.

Основные положения и результаты работы изложены в следующих публикациях.

1.Программная модель дегазационной системы выемочного участка угольной шахты. - Техника безопасности, охрана труда и горноспасательное дело: Реферативный сборник /ЦНИЭИуголь.-М.,1981г., Jfc'3f стр.16.

2.Разработка и исследование на ЭВМ алгоритма квазиоптимального управления процессом дегазации спутников угольного пласта.-Применение ЭВМ и математических методов в горном деле: Труды 17-го Международного симпозиума.- М., Недра, 1982г.,т.3, 307 стр. Соавторы: В.Н.Павлыш, О.Й.Федяев.

3.Постановка задачи математического моделирования движения газа в подработанном горном массиве с учетом дегазационных скважин, пробуренных с поверхности.- Прогноз и предотвращение газопроявлений при подземной разработке полезных ископаемых: Сб. науч.тр./ИПКОН АН СССР.- М.,1982, C.II3-II7.

4. Пакет прикладных программ для расчета движения газа в пористой среде методом конечных элементов.- Физико-технические и технологические проблемы разработки и обогащения твердых полезных ископаемых: Сб.науч.тр./ИПКОН АН СССР.- М., 1982, с.54-59.

5.Расчет движения газа в подработанном горном массиве методом конечных элементов,- Физико-технические и технологические проблемы разработки твердых полезных ископаемых:Сб.науч.тр./ /ИПКОН АН СССР.- М., 1983, с.45-49.

6. А.б. I0II856 (СССР) Способ дегазации подработанного угольного пласта.- Опубл. в Б.И., 1983, № 14.

Соавторы: А.Т.Айруни, М.А.Иофис, А.В.Шестопалов.

1. Лидии Г.Д. Газовыделение в угольных шахтах и меры борьбы с ним. -М.Дглетехиздат, 1952.-223 с.

2. Лидин Г.Д. Газообильность каменноугольных шахт СССР.Т.I,Факторы, определяющие газообильность шахт.-М.-Л.,Изд-во АН СССР, 1949.-342 с.

3. Айруни А.Т. Теория и практика борьбы с рудничными газами на больших глубинах.-М.,Недра, 1983.-335 с.

4. Айруни А.Т.,Слепцов Е.И. Дегазация угольных пластов за рубежом. -М.,ЦНИЭИуголь, 1973.80с.б.Вильчицкий В.В.,Морев A.M. Вентиляция и меры предотвращения взрывов газа и пыли в угольной промышленности США.-М.,ЦНИЭИуголь, 1976.-56 с.

5. Левин М.М. Исследование и выбор параметров дегазациии подработанных пластов и выработанных пространств в условиях шахт Карагандинского бассейна: Автореф.Дис. . канд.техн.наук.-Караганда,1969.-18 с.

6. Ералин"И.Е. Исследование и определение параметров дегазации сближенных пластов ивыработанных пространств Карагандинского бассейна на примере шахт, отрабатывающих пласты "Вышесредний" и "Слоистый": Автореф.Дис. . канд.техн.наук.- Караганда,1970. -24 с.

7. Таушкин Г.Т. Исследование ивыбор эффективных способов дегазации при разработке пластов верхней свиты в условиях Чурубай-Нурин-ского района: Автореф.Дис. . канд.техн.наук.-Караганда,1976. -21 с.

8. Белоусов В.Ф. Расчет заложения вертикальных дегазационных скважин, пробуренных с поверхности.- Безопасность труда впромышленности, 1981, №3,с.38-40.

9. Ю.Гайбович Ф.М. Эффективность дегазации подрабатываемых спутников и выработанных пространств скважинами с поверхности В кн.: Управление газовыделением средствами вентиляции и дегазации в угольных шахтах. ВостШШ, 1980,с.17-24.

10. П.Оспанов И.Ж. Эффективность комплексной дегазации при отработке пласта К^ "Нижняя Мариана" в Карагандинском бассейне В кн.: Рудничная аэрология: Сб.науч.тр./ ИГД им.А.А.Скочин-ского, вып.206.-М.,1982, с.25-32.

11. Лапин Г.К. Условия целесообразности применения дегазации сближенных пластов и выработанного пространства скважинами с поверхности. В кн.: Вопросы аэрологии угольных шахт: Сб.науч. тр./ ИГД им.А.А.Скочинского, вып. I98.-M., 1979, с.58-63.

12. Касимов О.И., Скворцов В.Г.,Рыбалко А.А. О целесообразности дегазации выемочных участков скважинами, пробуренными с поверхности.- Уголь Украины, 1978,№7, с.48-49.

13. Лейбензон Л.С. Движение природных жидкостей и газов в пористой среде. М.-Л., Гостехиздат, 1947.- 244 с.

14. Хб.Лейбензон Л.С. Собрание трудов, т.2. Подземная гидрогазодинамика.- М.,Изд-во АН СССР, 1953. 544 с.

15. Musket М., Botset М. Flow of gases through porous materials-Physics, 1934, v.I, No.I.

16. Коллинз Р. Течение жидкостей через пористые материалы М., i Мир, 1964.- 340 с.

17. Х9Ларный И.А. Подземная гидрогазодинамика.- М.,Гостоптехиздат, 1961.- 200 с.

18. Bear J. Dynamics of Fluids in Porous Media. American Elsevier

19. Publishing Company Inc. New York, 1972*

20. Годунов С.К.,Рябенький B.C. Внедрение в теорию разностных v схем. -М., Физматгиз, 1962.- 440 с.

21. Price H.O. A Computer Model Study of Methane Migration of Coal Beds.- The Canadian Mining and Metallurgical Bulletin, 1975 (Sept.), p.103-112.

22. Borrie D., de Vries G. A Finite Element Bibliography. Plenum. Hew Xork, 1976.2б.0ден Д.Ж. Конечные элементы в нелинейной механике сплошных сред.-М., Мир, 1976.- 426 с.

23. Стренг Г., Фикс Дж. Теория метода конечных элементов.-М., I, 1977.- 460 с.г28.Зенкевич 0, Метод конечных элементов в технике .-М.,: Мир , 1975 -360 с.

24. Irons В. Techniques of Finite Elements, Wiley & Sons, BY, 1979.

25. Gray W.G., Finder G.F., Brebbia O.A. (eds.) Finite Elements in Water Resources I, Pentech Press, London, 1977*

26. Brebbia C.A., Gray W.G., Pinder G.F. (eds.) Finite Elements in Water Resources II, Pentech. Press, London, 1978.

27. Noorishad J., Ayatollahi M.S., Witherspoon P.A. A Finite Element Method for Coupled Stress and Fluid Flow Analysis in Fractured Rock Masses. IJRM & MS, August 1982, vol.19, Bo. p. 185.

28. Mohsen M.F.N., Farquhar G.J., Kouwen N. Modelling Methane Migration in Soil. Appl. Math. Modelling, 1978, vol2 (December), p.294-301.

29. З^.Ляшко И.И.,Сергиенко И.В.,Мистецкий Г.Е.,Скопецкий В.В. Вопросы автоматизации решения задач фильтрации на ЭВМ.- К., Наукова думка, 1977.-288 с.

30. Schaeffer H.G. MSC/NASTRAN primer. Static and normal mode analysis. Mont Vernon, New Hampsire, Schaeffer Analysis. (1979)» 365 p.

31. Зб.Шайдуров В.В. Модульный комплекс программ MOK-I. Краткое описание. СО АН СССР. Препринт №18. Красноярск, 1980. » 1-53 с.

32. Zienkiewicz O.G. and Phillips D.V. An automatic mesh generation scheme for plane and curved surfaces by "isoparametric" co-ordinates. Int. J. num. Meth. Engng., 3, 5X9-528 (1971)•

33. Babuska I. and Rheinbolt W.C. Adaptive Approaches and Reliability Estimations in Finite Element Analysis. Computer Methods in Applied Mechanics and Engineering, vols.17/18, 519-540 (1979).

34. Carey G.F. Adaptive Refinement and Uonlinear Fluids Problems. Computer Methods in Applied Mechanics and Engineering, vols. 17/18, 541-560 (1979)•

35. Haber R., Shephard M.S. et al. A General Two-Dimensional Graphical Finite Element Preprocessor Utilizing Discrete Trans-finite Mappings. Int. Journ. for Numerical Meth. in Engineering, vol.17 (1981), p.1015-1044.

36. Tracy F.T. Graphical pre- and post-processor for two-dimensional finite element method programs. SIGGRAPH'77, 11(2), 8-12, (1977).

37. Wu S.C., Abel J.F. and Greehberg. Ail interactiv computer graphics approach, to surface representation. Commun. ACM, 20 (10), 703-712 (1977).4.3.Irons B.M. A frontal solution program for element analysis. -Int. J. num. Meth. Engng., v.2, 5-32 (1970).

38. Abbas S.F. Some novel applications of the frontal concept.-Int. J. num. Meth. Engng., vol.15, 519-536 (1980).45. ?i na H.L.G. An algorithm for frontwidth reduction. Int. Journ. num. Meth. Engng.,v.17, 1539-1546 (1981).

39. Razzaque A. Automatic reduction of frontwidth for finite element analysis. Int. J. num. Meth. Engng., v.15, I3I5-I324 (1980).

40. Outhill E. and McKee T. Reducing the bandwidth of sparse symmetric matrices. Proc. 24th National Conf. Assoc. Comput. Mech. ACM Publication P-69 (1969).

41. Cuthill E. Several strategies for reducing the bandwidth of matrices. In Sparse Matrices and their Applications (Eds. D.J.Rose and R.A.Will), Plenum Press, 1972.

42. Kamel H.A. and Shanta P.J. A solid mesh generation and result display package. J. Pressure Vessel Tech., ASME, 96(3), 207312, 1974.

43. Морев A.M.,Евсеев й.И. Дегазация сближенных пластов.-М., : Недра, 1975.- 168 с.

44. Айруни А.Т.,М.А.Иофис, А.И.Шмелев и др. Использование защитных пластов на угольных шахтах. Обзор ДНИЭйуголь.-М., 1981.- 33 с.

45. Brace W.F. Permeability of crystalline and argillaceous rocks— Int. J. Rock* Mech. and Mining Sciences, v«17, No.5» (October), 1980, p.241-251.

46. Кузнецов С.В.,Кригман P.H. Природная проницаемость угольных пластов и методы ее определения.- М., Наука, 1978.- 122 с.

47. Тарасов Б.Г.,Деменева А.И. К методу расчетного определения пористости горных пород по их сдвижениям.- В кн.: Вопросы рудничной аэрологии. Кемерово, Кемеровское книжное изд-во, с.127-139. (КузПИ).

48. Сегерлинд Л. Применение метода конечных злементов.-М.,Мир, 1979»- 388 с.

49. Таблицы физических величин. Справочник.-М.,Атомиздат,1979. -546 с.

50. Дополнения к руководству по проектированию вентиляции угольных шахт".-Недра,1981.-77 с.

51. Бухны Д.И.,Лапин Г.К. Оценка эффективности дегазации вертикальными скважинами, пробуренными с поверхности.-Техника безопасности, охрана труда и горноспасательное дело. Сборник/ ЦНИЭйуголь.1983, $2,стр.34.

52. Лапко В.В., Федяев 0.И.,Касимов 0.И. Исследование параметров ' дегазации подработанных угольных пластов на математическоймодели.-Известия ВУЗов, Горный Журнал, 1980, №9,стр.63-69.

53. Касимов С.0.Расчет движения газа в подработанном горном массиве методом конечных элементов.-Физико-техниееские и технологические проблемы разработки твердых полезных ископаемых: Сб.науч.тр./ИПКОН АН СССР.-М.,1983, стр.45-49.

54. A.C. I0II865 (СССР) Способ дегазации подработанного горного массива /А. Т.Айруни,М.А.Иофис,А.В. Шестопалов, С. О.Касимов. -Опубл. в Б.И., 1983,Ж4.