Бесплатный автореферат и диссертация по геологии на тему

Разработка методики межсважинной сейсмической томографии с целью выявления залежей углеводородов

ВАК РФ 04.00.12, Геофизические методы поисков и разведки месторождений полезных ископаемых

Автореферат диссертации по теме "Разработка методики межсважинной сейсмической томографии с целью выявления залежей углеводородов"

ГОСУДАРСТВЕННЫЙ КОМИТЕТ ДО НАРОДНОЕ ошзошла_________ -

_________МОСКОВСКИЙ ОРДЕНА ОКТЯБРЬСКОЙ РЕВОЛЮЦИИ И

ОРДЕНА ТРУДОВОГО КРАСНОГО ЗНАМЕНИ ИНСТИТУТ НЕФТй й ГАЗА им. И.М.ГУБКИНА

На правах рукописи УДК - 550.831

БУЗЕХ АБДУ1РАХМАН

РАЗРАБОТКА МЕТОДИКИ МЕдСКШИННОЙ СЕЙСМИЧЕСКОЙ ТОМОГРАФИИ С ШЬЙ ВЫЯВЛЕНИЯ ЗАЛШй УГЛЕВОДОРОДОЗ

Специальность 04.00.12 - Геофизические мзтоды поисков к разведки месторождений полезных ископаемых

АВТОРЕФЕРАТ

диссертации на соискание ученой иедеаи кандидата технических наук

Москза - 1991

Работа выполнена в Московской институте нефти и газа имени П.¡¿.Губкина на кафедре полевой геофизики.

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

профессор М.Б.Рапопорт

Научнай консультант: кандидат технических наук,

доцент А.М.Курсов

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

О.А.Погапов, зав. отд. ВНИИГеофизйКИ

кандидат геолого-минерадоги-• ческих наук Е.А.Ефшова, с.н.с. МГУ

Ведущее предприятие: Институт физики Земли АН СССР

Защита диссертации состоится "2 5"" ¿/.¿^Н^ 1991 г. в/Г часов в ауд.^З-Зка заседании специализированного совета Д.053.27.08 по защите диссертаций ва соискание ученой степени доктора наук при .¡¿аш1В«ю14-инсти2уте"яефти~ и -газа имени И.¿¿.Губкина по адресу: 11791?, г. Москва, ГСП-1, Ленинский проспект, 65, т. 93СК)3-81.

С диссертацией можно ознакомиться в библиотеке ¿¡ИНГ имени И.М.Губкина.

Автореферат разослан

"?/и ЛоЛ 1991 г.

Ученик секретарь специализированного

совета, доцент < Н.Н.Кравко

ОБЩАЯ ХАРАК'ШИСТЙКА РАБОТЫ

Актуальность темы. В настоящее время, несмотря на широкое пряиенение геофизических поисково-разведочных райог на ._____и таг^зачастую не обеспечивается редейте задачи поисков и разведки неструктурных и небольших по размеру задезей углеводородов (УВ). особенно в районах с неблагоприятными, сейсмогеологичесыши условиями, при сложном строении верхней части разреза (ВЧР). Традиционнее методы полевой геофизики в многих случаях не позволяю! получить достаточно информации для из'чеяан агиенений фнзетеских свойств горных пород, связанные с залежами УВ и их изменением в процессе разработки. Методы ГИС обеспечивай; достаточно точные оценки таких свойств лишь в стволе скважины. •

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

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

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

1. Математическое моделирование волновых полей при меЕСква&ин-них наблюдениях.

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

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

Оценка точности и надежности способов восстановления изображений среды по методу межскважинной сейсмической томографии (ШСТ).

5. Разработка рекомендаций по интерпретации результатов ШЙТ. Научная новизна.

1. Показана возможность изучения в слоистой среде методом межскважинной томографии (ШСТ) неоднародностей, которые могут быть вызваны залежами У£.

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

моделированием,

3. Показана возможность определения параметров среды, кон-lypa неоднородности (залеаи УВ) и контроля за его изменением на основе разработанной автором методики Ш1СТ,

Практическая ценность работа заключается в том, что сформулированные критерии выделения и оконтуривания локальных неоднородностей методом МСТ обеспечивают повышение эффективности всего процесса поисков и разведки залекей УВ, способствуют открытию новых и доразведке ухе открытых месторождений при ограниченной-скоростной дифференциации пород, на-прииер, в теригенном разрезе. Сейсмическая томография может сдуаить инструментом контроля процесса интенсивной добычи УВ. Разработанный комплекс программ реиения прямой задачи на ЭВМ позволяет вычислять данныз метода ыежскважиивого сейсмического просвечивания (МОЕ) для различных моделей реальных слоистых сред, а пакет программ обработки полевых данных метода МСТ позволяет получить более полную и достоверней информацию с скоростях и поглоцаюздх свойствах слоисто-неоднородных сред. Выполненная программная реализация разных алгоритмов восстановления параметров волнового поля в ММСТ позволяет легко адаптировать пакет програмн обработки к ыеняюциыся реальным условиям.

Приведены рекомендации по интерпретации и применения ЫМСТ в различных сайсмогеологическиг условиях.

Защищаемые положения. I. Метод иеаскзажинной сейсмической томографии позволяет изучать залежи УВ в слоистой среде в случае ограниченной (д<? 2($) контрастности с/.оростей на границах слоев, что обычно наблюдается з терригенном разрезе.

2. Наиболее точное восстановление изображения в МНСТ обеспечивает алгоритм.одновременного алгебраического восстановле-ния*($ЖТ)'при условии, что его начальное приближение получается как результат воссгагэвления по методу обратных протекций, не требующего априорных знаний о модели среды.

Апробация работы и публикации. Основные результаты исследований докладывались на:

- Всесоюзной конференции "Роль молодежи в решении конкретных научно-технических проблем нефтегазового комплекса страны", МйНГ им.Губкина, 19В9;

- Региональном соЁевднии "Проблемы техногенного изменения геологической среды, и охраны недр в горнодобывающих регионах", Пермь, 1990;

- П Всесоюзной вколе-сезшнаре "Разработка месторождений нефти и газам4кд)й!енЕое состояние, проблемы, перспективы11, г.Звенигород, 199®.

Результаты исследований включены в научный отчет лаборатории- нефтегазовой сейсморазведки ШШГ за 1990 г.

«и

По теме диссертации опубликовано 3 работы.

Объем работы* Диссертация состоит из введения« четырех глав и'заключения и содержит О страниц машинописного текста, 3/ рисунков. Список литературы' включает^/ наименований.

Работа выполнена на кафедре полевой геофизики ШНГ ии.й.М^Губхшт.

Автор выражает глубокую благодарность за научное руководство сроф.Ы.Б.Рапспоргу, а-также: доцД.МЛукову, доц.В.&.Ба-кирову, доц.Е.Б.Варозу, доц.Г.А«Кар&петову е другак сотрудни-

?

каи кафедры за консультации, постоянную лоддергку и внимание к работе, ______________________________________*

СОДЕРЖАНИЕ РАБОТЫ

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

-Применительно к поискам нефти и газа.лзрвае наследования по мЗЕСКвааяннои «оиографаи проводились в 1971 году Институтом нефти Франции [Бойс и др. 1971, 1972s] с целью изу-

-ч -ч

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

Многие зарубеаные исследователи Дайне я др. 1979, йзан-сон и др* 1985, 1986, Брегман и др. 1990 и др.) применяли нотой кезсквааинной томографии в основной для изучения кристаллических пород и реаения задач рудной, угольной и инженерно-геологической геофизики» Есть примеры применения метода для контроля за изменением физических свойств горных пород при • закачке пара с целью повышения нефтеотдачи [Накридес а др. 1988, Брегман и др. 1990 ] . Однако в зтих работах констатировалась только возможность применения кетода для случаев, когда расстояние ыезду сквааиаами было 51 м и 200 м, но ¿о

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

Исследования по методу мексквааинного сейсмопросвечива-ния и сейсмотомографии встречаются в работах и многих советских ученых [Ефимова Е.А. 1982, 1989, Рудерман E.H. Файзуллин И.С. 1982, 1965, Силаев В.А. 1986, 1988, Караев.'й. Д. 1987, Бельфер'и.К. 1986, 1989 и"др.]"' .

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

Анализ работ по методу мехсквахинного просвечивания, осо-' бенно с• использованием томографии, показывает, что метод представляет большой интерес длн изучения месторождений УВ в ues-сквахинном пространстве.

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

Вторая глава содеркат постановку задачи в МУСТ для'решений геологических и технических проблем, которые сложно решать другими методами, когда условия в ВЧ? неблагоприятны для возбуждения и приема сейсмических колебаний.

Решение прямой задачи - расчет волновых полей - является важным этапом при исследования новых методов на стадии моделирования. В работе обоонован_выбор-моделей-реальных сред. Использовались расчетные числовые модели, описывакщие свойства к параметры геологической среды (геометрия пластов, форма залека, скорости и коэффициенты поглощения упругшс волн). Параметры скоростей и поглощения в залезах подобраны на основе реальных данных для нефтегазовых месторождений, когда ско^ рости распространения сейсмических волн отличаются от. законтурных областей не белее, чем на 20 а позызение коэффициентов поглощения происходит почти на порядок. Приведенные геологические модели выбраны на основе данных, приводимых в работе Гальперин и Мирзоян, (д{?б .

При моделировании расстояние мэвду сквазинами выбрано порядка 800 ы. Это расстояние является реальным для разрабатываемых нефтявтк иесторовдений во многих районах, Мощность слоев модели выбиралась около 40-70 метров, исходя из реальной разреоащей способности для сейсмических волн низкой часты (50-50 Гц) при скорости распространен:« волн 2,5-3,0 ха/с* Мощность залеяи разнялась 20-40 м$ что превышало 1/2 первой зоны Френеля, представляющей разрешающую способность метода.

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

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

случаев при ограниченной контрасте окоросгей (до ± 20#) пря-

1 *

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

Восстановление изображений среды (глава 5) вклэчает и оценку ее физических характеристик по их проекциям. Математическое описание такого восстановления выражается как

У « АХ

где У - известные выходные данные (измерения);

I - неизвестные входные данные (параметры среды);

А - известная матрица проецирования (т.е. приращение лучей в клетках).

Задача восстановления'с самого начала решается в дисвфетг ной форме. Введем в интересующей нас плоскости прямоугольную сетку таким образом, чтобы она покрывала всю исследуемую часть плоскости, подлежащей восстановлению. Если клетки имеют номера от I до N , а лучи, по которым производится интегри-розаняе, от I до Ж, то можно переписать систему линейных уравнений в дискретноы виде

N ■ . '

II

'V4

Если У/ представляет время пробега, то ij представляет обратную величину скорости, а если У/ - амплитуда, то Ij -определяет коэффициент поглощения._______________________ ______

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

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

Реальный путь рзшеная задачи восстановления возможен на основе методов разложения в ряд и решения системы уравнений итерационными способами. Разработаны алгоритмы решения такой задачи итерационные методами. Наиболее эффективными для решения представляются два алгоритма: а) Алгебраические алгоритмы реконструкции {АЦгбыхсс Recon^iucti'on Techm-pusi* п ЯI ) и алгоритм реконструкции с одновременными итерациями {Simufra-

neu-i Tteiattire Reconstruct/с* Technique - SIR i Jc

Алгебраический алгоритм восстановления AR. I предлохеа Городом и др. (Х9?0) для приближенного решения больной системы уравнений. Для достижения этой цели производится итерационная процедура, в которой используется некоторое начальное

♦ч /

приближение параметров, подлежащих восстановлению (Х°).

Такое начальное аначение можно задазать на основе апри-' орных данных о модели. На 9--Й итерации текущее приближение 1 ^ уточняется, что дает новое приближение Х*^ , при этом

рассматривается только один луч, .например <■ -й, а изменению

а.

подвергаются только те компоненты вектора X , где проходит этот луч. Величина невязки между измеренным значением У

Д V V-

и расчетной приближенной зеличиной проекции \ О-е] А/ ,

Ц.

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

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

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

Алгоритм восстановления 5разработан Гельбер-

том и др.(1972) для решения больших разряженных систем урав-■ * -

нений и основан на одновременной итерации по всей сетке для

минимизации нормы. Для повышения качества изображения целе-

собразно уменьшать дисперсно; N _ о

N

представляет собой среднее значение обратной скорости или. поглощения в дискретизированной плоскостей восстановления.

Значение X не корректируется до тех пор, пока поправка не определена для всех клеток исследуемой плоскости.

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

__________терпретации восстановленных изображении.

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

1. для каздого луча ( -с") и "для каздой клетка ) вычисляются предварительные значения У^с^ ~

где У - наблюдаемое .значение;

- длина -с -го луча от источника до приемника.

2. Вычисление окончательного значения % ^ в 3 -й шегке как средневзвешенного значения

где О-с'з - длина отрезка < -го луча в J -а клетке;

К - число лучей, проходящих через 1 -в клетку. При программной реализации алгоритма А ЯЛ" искомому значении двумерной функции X в } -й щетке присваивается некоторое начальное значение, гелатеаьно близкое к реальким значениям. В диссертаций предлоаено использовать е качестве начального приблиа'ення значение Х^ •' з получзютз це иеходу обратно а прочкция.

ч- ' '

. Вычлслвегся оцевза У; геометрического интеграла вдоль

луча < ( Я - номер текущей итерации):

где йц - длина -го луча в ц -й клетке.

Поскольку !у не являются точными значениями, то У^ Для того, чтобы У^ «= У/^^ , к X] при-

бавляется поправочный коэффициент ДХ^

Д у _ V У

1

j — 31 :

j=-t

Таким образом, в каждой клетке j получаем новое значение параметров

'?+/ ? ' >

Ху Г Ху -V AX;

я

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

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

При реализации S7RT , как ив AR I , итерации на-чияаются с каких-то начальных- значений X j , которые улучшаются в течение итерационного.процесса.

Одна итерация STRT п0 прямым лучам состсшг из следующих шагов:

I. Для каждого <■ -го луча рассчитывается время пробега воле от источника до приемника

* Р- г

Х- - I а О- ху j- < J

Я - номер итерации. 2ь Определяется разница между р_ассчи1аннмм-и-наблпдавныц"вре~" менем для <"-го"луча

П. ^

Рассчитываются поправки для каждой J -ой клетке по формуле

А*6 - и V !г/ .. I2

После обработки всех лучей для каждого значения X в каждой Л -й клетке добавляется поправка А Ху ч Я

X; - X; + Д*,-

В результате одной итерации получаем козые аенчениг. яо-зомнх параметров. Если они не удовлетворяют задаяыка критериям, итерация повторяется. Б качестве критерия сходимости реакция используем меру с (погрешность)

Д - 2

К**-*>]

.) = <

где Ы - число клеток в плоскости;

А

Х^ - значение функции при Я итепациу • - значение функции при Ч--Н итераций« В качестве критерия сходимости моано применять сумму невязка между наблодаеиым и расчетный значениями параметра М

а - число лучей»

С целью выбора оптимальных условий применения метода МСТ было исследовано влияние исходных данных, полученных как с учетом,так и без ^чета преломления на границах. Показано, что влияние преломления аналогично влиянию случайного иума (ошибок регистрации). Эффект от преломления в пластах с 20&

х

контрастом скорости на границах слоев эквивалентен, 10$ уровню

•ч

иума в наблюдаемых данных. Это позволяет использоьать алгоритмы томографии, основанные на прямолинейности лучей (без учета преломления при решении обратной задачи). Потеря точности восстановления отмечается в тех местах, где локальный контраст по скоростям или по поглощению мал по сравнению с соседними участками. Другим важным фактором, влияющим на точность восстановления, является степень покрытия изучаемой среды лучами: отсутствие лучей в некоторых клетках приводит к неоднозначности решения для них. Точность решения по методу ИСТ зависит такав от .отношения высоты аномалий в среде к их ширине. Чем это отношение больше, тем'решение по ИСТ хуке. Б нефтегазовой геологии ато отношение обычно небольшое, т.к. пласты и ловушки располагается в горизонтальном направлении. В целом, качественное решение' по системе мегскважянныг наблюдений позволяет определить геометрию аномалии' (месторождения) на стадии поио-. ков или на стадии разработки.

Значения физических- параметров определяются с точностью до 20/ь в зонах аномалий и с еще большей точностью в простых структурных средах.

На основе полученных результатов выработаны следующие рекомендации по применению методов ЙСТ для выявления залезей УВ;

I» Метод Ш£СТ целесообразно применять при поисках, разведке

и разработке залежей УВ для их изучения в цеаскваетнноа пространстве в тех районах,где^сейсыичеокие^наблюдения на поверхности ату задачу не позволяют решить, а контрастность скоростей на границах сейсмических слоев не слетком велика (до 20&). Последнее услозие обычно выполняется в терригенном разрезе.

2. Система наблюдений должна выбираться с учетом априорных данных об изучаемых залогах (сейсмических наблюдений на поверхности, ГИС). .

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

Число источников и.приемников в сквааинах доляао быть не меньше числа клеток по вертикали.

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

6. Для получения окончательного топографического изображения рекомендуется метод одновременного алгебраического восстановления (SIR "П.

7. При интерпретации результатов восстановления следует учитывать .возможность появления тени объекта (в зависимости от системы наблюдений).

8. На краевых участках изобрааения, перзеев&еши. небольшим числом лучей, наблюдается неоднозначность восс?п."сг-сзяя, г,о~ згоау задай участки исключаютсг при интерпретации результатов HMGl. , .

9. Применение процедур сглаживания мевду итерациям повывает качество восстановлении при налички ыума г. иохедянх дашшх

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

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

ЗАШЧЕНИЕ

Широкое внедрение метода межскважинной томографии в СССР сдерживается из-за отсутствия удобных скважинных источников сейсмических волн. Силаев В.А. использовал для возбуждения волн взрывы торпед, приводящие'"к разрушению скважины, что неприемлемо. Но в настоящее время в ряде организаций разрабатываются скважинные источники импульсного и вибрационного типа, что устранит препятствия к внедрению ШСТ. Применение ШСТ'в комплексе с ЕСП и ГИС с использованием материалов наземной сейсморазведки позволит детально изучать околоскваашшое пространство, наращивая столь необходимые запасы нефти и газа, а повторные исследования в ходе эксплуатация месторождений позволят контролировать активные воздействия на пласт, необходимые для повышения нефтеотдачи.

Б настоящее время обработка данных на ЭВМ производится в крупных центрах, на больших расстояниях от месторождений нефзи и газа, что непригодно при использовании скважинной и межскважинной сейсморазведки в технологии разработки месторождений. Но. в ШНЕ уже разработаны полевые вычислительные комплексы на базе персональных ЭВМ, что позволяет оперативно обрабатывать данные на месте.

Метод ММСТ при использовании новых скважинных источников, многоканальной приемной аппаратуры и обрабатывающих комплексов

удовлетворяет требованиям современной технологии,__что_обеспечи-------

________вает- его-широкое" внедрение.

Основные результаты диссертации опубликованы в следувацас работах;

1. Бланк A.M., Бузех А., Хананоь И.£.-Сейсмические исследования в процессе бурения скважин. Тезисы доклада на Всесоюзной конференции "Роль молодежи в решении конкретных научно-технических проблем нефтегазового комплекса страны" (6-8 ияня 1989, поп.Красный Курган, стр.2).

2. Бузех А., дуков A.M., Рапопорт Ы.Б. Возможность контроля состояния геологической среды методом сейсмического меисква-аивного прозвучивашш. Тезисы доклада на Региональном совещании "Проблемы техногенного изменения геологической среды и охраны недр в горнодобывающих рсрвонах" (23-25 октября'1990, стрЛ29--13Х).

3. Рапопорт Ц.Б«, Бузех А., Дон Б.Н., Рыжков В.Н. Прогноз зон промышленной продуктивности я контроль га разработкой по данным сейсиорззведки.'_Сб.трудов П Всесоюзной школы-сеианара "Разработка месторождений нефти а газа. Современное состояние, проблема и перспективы" (11-16 марта I99X, г.Звевэдюрод).