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

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

0046

7941

Авроров Сергей Александрович

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

Специальность: 05.13.17 - Теоретические основы информатики

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

1 6 ЛЕН 2010

Новосибирск - 2010

004617941

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

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

Хайретдинов Марат Саматович

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

Попков Владимир Константинович Кандидат физико-математических наук Мурзин Федор Александрович

Ведущая организация - Московский физико-технический институт

Защита состоится 23 декабря 2010 г. в 16 часов на заседании Диссертационного совета Д 219.005.02 при ГОУ ВПО "Сибирский государственный университет телекоммуникаций и информатики", по адресу. 630102, г. Новосибирск, ул. Кирова, д. 86, ком. 625.

С диссертацией можно ознакомиться в библиотеке ГОУ ВПО "СибГУТИ".

Автореферат разослан ноября 2010 г Учёный секретарь

диссертационного совета Д 219.005.02 кандидат технических наук доцент

т-

И. И. Резван

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

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

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

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

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

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

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

Решению указанных задач способствуют фундаментальные исследования в этой области отечественных ученых: Алексеев A.C., Адушкин В.В., Бугаев A.C., Глинский Б.М., Гущин В.В., Гордиенко В.А., Еманов А.Ф., Ковалевский В.В., Михайленко Б.Г., Овчинников A.M., Омельченко O.K., Хайретдинов М.С., Филина А.Г. и ряд др., а также зарубежные ученые: Карлтон П.Н., Фу К.С., Херрин Ю., Чен С.Н. и др.

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

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

В соответствии с целью определены следующие задачи:

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

2. разработка алгоритмов обнаружения и измерения параметров сейсмических волн на фоне шумов в поточном режиме. Исследование их помехоустойчивости и точности оценивания параметров;

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

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

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

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

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

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

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

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

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

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

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

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

® мест падения с воздуха боевых средств поражения;

• эпицентров скрытых подземных ядерных взрывов по авторшоковой эмиссии (инспекция на месте в интересах Договора о всеобъемлющем запрещении ядерных испытаний) и др.

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

Реализация и внедрение результатов, работы. Результаты диссертации составляют основу программного инструментария, применяющегося в отделе геофизической информатики Института вычислительной математики и математической геофизики (ИВМиИГ СО РАН) при проведении лабораторных исследований и полевых экспедиционных работ.

Диссертационные исследования выполнялись в рамках проектов Российского фонда фундаментальных исследований (РФФИ) №07-07-00214а, 09-07-12075-офи-м и междисциплинарного проекта СО РАН №16_3.

Исследования алгоритмов обнаружения, измерения параметров и локации тестовых полигонных взрывов выполнялись в рамках хоздоговора "Полесовщик" между НГТУ и ОАО "Завод Ленинец" (г. Санкт-Петербург). Полученные результаты исследований используются в составе технологии автоматизированной локации полигонных взрывов падающих боеприпасов и карьерных взрывов.

Достоверность полученных результатов подтверждается проведёнными натурными экспериментами и математическим моделированием, а также экспертизами отчетов по хоздоговору и грантам РФФИ.

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

-Международной научной конференции «БГОГОСОЫ» (Новосибирск, 2008 г., Иркутск, 2010 г.).

-Всероссийской научной конференции молодых ученых «Наука. Технологии. Инновации» (г. Новосибирск, 2006, 2008 гг.).

-Конференции молодых учёных ИВМиМГ СО РАН (г. Новосибирск, 2008, 2009 гг.).

-Международных научных конгрессах «ГЕО-Сибирь» (г. Новосибирск, 2007, 2008 гг.).

-Международной конференции «Мониторинг ядерных испытаний и их последствий» (п. Боровое, Казахстан, 2008 г.).

-VI Всероссийской научно-практической конференции студентов, аспирантов и молодых ученых (Томск, 2008 г.).

-Международном семинаре "Акустика неоднородных сред-Х" (Новосибирск, 2009 г.).

-Научной сессии НГТУ (г. Новосибирск, 2008, 2009 гг.).

Публикации. По теме диссертации опубликовано 8 научных статей, а также тезисы докладов. Опубликована 1 статья в рецензируемом сборнике «Научный вестник НГТУ» № 2, 2010 г. Результаты исследований отражены в отчётах по грантам РФФИ, СО РАН и хоздоговору.

Основные положения диссертации, выносимые на защиту.

Комплекс алгоритмов и программ, включающих в себя:

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

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

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

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

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

5. Экспериментальные оценки точности локации сейсмического источника на основе измерения параметров сейсмических и акустических волн, порождаемых малоудалёнными полигонными взрывами, а также Р- и Б-волн от удалённых промышленных взрывов.

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

Структура и объём диссертации. Работа состоит из введения, четырех глав, заключения и списка литературы из 65 наименований. Основное содержание представлено на 144 страницах, содержит 43 рисунка и 16 таблиц.

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

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

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

При решении задачи геоакустической локации выделяют ряд основных этапов:

• выделение сейсмических и акустических волн в шумах и измерение их времен вступлений;

• идентификации (распознавание) типа источника;

» определение параметров очага;

• отображение координат источника на цифровой карте.

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

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

500 т в тротиловом эквиваленте. При этом перед началом каждого сеанса регистрации необходимо расстановка профиля регистрирующей аппаратуры в направлении прихода ожидаемой волны. Такая возможность обеспечивается за счёт использования мобильных регистрирующих систем ( МСГ). МСГ располагается на относительно небольшой площади (порядка 3 на 3 км), что позволяет организовывать систему сбора данных на основе обычных проводных каналов и при этом получить точность определения параметров сейсмических событий на уровне, сопоставимым с обычными региональными сейсмическими сетями при значительно более высокой оперативности получения результатов.

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

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

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

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

Т=т}(Х,в) + е, (1)

где / =((,,(,,- вектор времен пробега сейсмических сигналов,

п(х,в) - Л'-мерный вектор вычисляемых времен пробега (теоретический го до-

7*

граф) или функция регрессии, е = .....¿Гд,) - вектор невязок,

в = (х, у, г, V, ()Г - /и-мерный вектор оцениваемых параметров, X = (5) ,х2,...,хн) - матрица координат датчиков (или точек излучения), Ы- число датчиков (или точек излучения).

В качестве оцениваемых параметров выступают пространственные координаты источника-*,^, скоростная характеристика среды v и время в источнике I. В ряде случаев скорость в среде является известной.

Решение уравнения (1) сводится к решению обратной задачи. При этом точность решения зависит, в первую очередь, от ошибок оценивания вектора времен вступлений волн t=(i,,i2.....tn). скорости волн v, ошибок измерения

е =(el,...,cN)T. выбора геометрии расстановки датчиков на дневной поверхности

Земли. В частности, по отношению к полярной системе координат дисперсия ошибки определения азимута Az на источник и расстояния R «источник-приемник» с помощью триады сейсмостанций (случай N=3) определяется соотношениями:

о!я=а!Р,(1#), cri=aÍF/i.v) (2)

где ст/ - ошибки оценивания времён <, ц/ - вектор параметров, характеризующий геометрию расстановки сейсмической группы.

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

-тл-, (3)

2Af(2E/NJ

где г,., - длительность волнового импульса, Af- ширина его спектра, E/No -отношение энергии импульса к спектральной плотности внешнего шума.

Как видно из (2) от точности измерений времён вступлений волн зависит точность определения параметров источника. Трудность решения задачи усугубляется влиянием шумов, имеющих техногенную природу.

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

В частности, если сигнал на выходе первого датчика принять равным u, = F3 -coscot, тогда, в случае п датчиков максимальный выигрыш по чувствительности для мультипликативной обработки составляет:

TK~F;-¥Jk.e,d)t (4)

и

П 1 2

для аддитивной обработки:

и1 =«, +м2 + ... +ип ~ п-Рэ (5)

где с! - расстояние между соседними датчиками, Е, - характеристика чувствительности по направлению одного датчика, 6 - угол между направлением прихода волны и перпендикуляром к линии расположения сейсмоприемников, к = 2п/Х - волновое число.

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

где и: - сигнал на ¡-м датчике сейсмической антенны, п - количество датчиков, Л^ - шаг перебора по ожидаемым углам прихода волны, параметр] определяет диапазон перебора углов от 0 до М.

На множестве сигналов (сейсмограмм) (О, где ¡=1, ...п, в плоскости О, О оптимальному направлению соответствует положение годографа - линии вступления волн. Угол наклона годографа к оси времен определяет скорость распространения данного типа волны. Процедура поиска истинного годографа осуществляется путём последовательного перебора. Учитывая различие в траекториях распространения основных типов сейсмических и акустических волн по положению годографа возможно осуществлять избирательную регистрацию волн того или иного типа.

При помощи моделирования произведены оценки качества работы алгоритмов (6), (7) по показателю «сигнал/шум» в зависимости от дисперсии шума по отдельным каналам.

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

с параметрами а-1, п=I, а=20, /о=42 Гц. Данный набор параметров функции Берлаге приближённо соответствует частотно-временным свойствам сейсмических волн от полигонных взрывов, рассматриваемых в данной работе. Это соответствует одной из прикладных задач геоакустической локации.

В результате модельного эксперимента алгоритм мультипликативной обработки (6) показал более высокое значение сигнал/шум на выходе, по сравнению с традиционно применяемым аддитивным (7). Кроме того, в соответствии с (3) мультипликативный алгоритм обеспечивает меньшую ошибку оценки времени вступления импульса за счёт более широкого спектра результирующего импульса (импульс после мультипликативной обработки (6) более "острый" по сравнению с результатом аддитивной обработки).

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

(7)

(6)

у(1) = (а1пе"' )г*т(2х и)

(8)

(7) вдоль предполагаемых годографов волн, соответствующих предполагаемым скоростям основных типов волн. При этом времена вступлений продольных волн Р и поперечных Б соответствуют максимумам обеих функций. На рис. 16 в увеличенном масштабе приведен фрагмент записи, содержащий взрывную сейсмограмму. На рис. 1в приведен результат ее выделения с помощью аддитивного метода, на рис. 1г - тот же результат, но полученный с помощью мультипликативной обработки.

• ............. /........г Волна Р -v___ ........i................г ......................i..... 1 1 1 f Волна S , }

20 43 ее so ico 120 1«) leo tso

г ,.......1 1 . | ;; ■

■'i во во «Ю 120 140 isú 1SC

<МЬ i j

OOi Волна Р Волна S

(б)

(в)

(г)

Рис. 1. Результат обработки реальной многочасовой сейсмической записи: (а) исходная сейсмограмма; (б) - участок сейсмограммы в увеличенном масштабе, содержащий искомую волну; (в)- результат выделения волны методом аддитивной обработки; (г)-то же с помощью мультипликативной обработки

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

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

ции распределения энергии анализируемого сигнала по уровням разложения с учётом порядка вейвлета. Применительно к* описанной модели сейсмической волны при использовании вейвлета Добеши 12 порядка более 90% энергии сигнала сосредотачивается на 2-м уровне разложения.

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

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

Алгоритм STA/LTA производит обработку уже отфильтрованного сигнала при помощи двух скользящих во времени окон («скользящее среднее»): короткое усредняющее окно (Short Time Average window — STA) и длинное усредняющее окно (Long Time Average — LTA). STA вычисляет оценку значения «мгновенной» амплитуды полезного сигнала (например, Р-волны промышленного взрыва), LTA оценивает среднее значение шума на длительном участке. Далее производится вычисление соотношения STA/LTA двух значений, вычисленных для каждого нового отсчёта сигнала, поступающего на вход алгоритма.

Рассмотрим цифровую однокомпонентную сейсмограмму y(t), где t=ih (¿=1,2,...) - дискретное время, h - интервал дискретизации. Короткое и длинное временные окна синхронно скользят по волновой форме и в каждой позиции вычисляются средние амплитуды:

r.=STA,/LTAn

где NS и NL- число отсчетов в коротком и длинном временных окнах, соответственно (NS«NL), гi - отношение STA/LTA. Момент времени Т0 = jh, когда

rj>C>l (9)

(С - порог обнаружения) объявляется началом полезного сигнала.

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

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

держащая 50 модельных импульсов с последующим наложением шума. На первом этапе производилась многоканальная (аддитивная и мультипликативная) обработка, далее применялся алгоритм БТА/ЬТА с параметрами 8ТА=0.1 с, ЬТА=1 с. Значение порога обнаружения импульсов для аддитивного метода составляет 1.1, для мультипликативного - 2, для одноканального метода с использованием вейвлет-фильтрации -1.1. Для каждого отдельно взятого модельного

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

■ Дздотв-ьй кежщ (9 кэ^гсез) !

ВЭ ГДтьтити^туеньй мзгсд (9 кэкегсе) | в 1\Л?ттве№гетч1ътътаал'И(11С* ГГТ1

16.00 8.00 4.00 200 1.00 0.50 0.25 В>од-ое "си>агУщ»/Г

Рис. 2. Диаграмма погрешностей вычисления времён вступлений волн а (Дг) для многоканальной обработки (аддитивного и мультипликативного методов) и одноканального метода вейвлет-фильтрации при различных входных значениях сигнал/шум; данные модельного эксперимента.

Результат вейвлет-фильтрации модельного сигнала приведён на рис. 3. Результат вычисления функций (6) и (7) применительно к модельному сигналу приведён на рис.4.

4

-<1—•-^-1||-1м-^

\ г

Рис. 3. Результаты одноканальиой обработки модельного сигнала: (а) исходная запись (соотношение сигнал/шум составляет 1/4), (б) - результат вейв-лет-фильтрации, (в) - результат сглаживания алгоритмом БТА/ЬТА с параметрами 0.1 с/1 с.

133 7 13*1 135 1 135 8 136»

Рис. 4. Результаты многоканальной обработки модельного сигнала по 9 каналам: (а) исходная запись (соотношение сигнал/шум составляет 1/4), (б) - результат аддитивной обработки, (в) - результат сглаживания алгоритмом БТА/ЬТА с параметрами 0.1 с/1 с; (г-д) - результат мультипликативной обработки с последующим сглаживанием БТА/ЬТА с параметрами 0.1 с/1 с соответственно.

На рис. 4в,д видно, что мультипликативная обработка выдаёт более контрастный момент вступления искомого сигнала по сравнению с аддитивной, что является экспериментальным подтверждением выражения (3): дисперсия оценки времени вступления импульсного сигнала обратно пропорциональна ширине спектра импульса.

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

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

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

В результате изучения особенностей сейсмических записей взрывов была выбрана система из четырёх признаков [4-5]:

• отношение амплитуд сейсмических и акустических волн (рис.5а);

• отношение средних частот этих волн РаЛ^ (рис. 56);

• задержка между временами вступления сейсмических и акустических волн Та-ТБ (рис. 5г);

• длительности сейсмических и акустических волн.

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

Рис. 56 отражает результаты анализа средних частот преобладающих спектров для обоих типов волн. Видно, что основные частоты сейсмических волн лежат в диапазоне 15-18 Гц, а акустических волн - в диапазоне 25-40 Гц. Это явление связано с тем, что высокие частоты затухают в земле быстрее, чем в воздухе.

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

2000,000 4000,000 5 ООО. ООО

Время обнаружения, с

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

(б)

5

О -

О. ООО

¿ОСО.РОО 4000,000 5 ООО. ООО Бремя обнаружения, с

Длительности сейсмических и акустических волн

Экспериментальная зависимость расстояния Я от разности Та*Т$

ю 2.5 £

I .54

г

I "

О Семей. ^ А руст.

|(В)

2000,000 4000,060 6000,000 Время обнаружения, с

Та-Тз.с

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

На Рис. 5г показана экспериментальная зависимость расстояния «источник-приёмник» II от разности времён вступлений сейсмических и акустических волн.

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

По аналогии с рассмотренными был осуществлен выбор информативных признаков по отношению к удаленным промышленным (карьерным) взрывам [6-9]:

® отношение амплитуд продольных Р и поперечных Б волн;

© отношение средних частот этих волн- Рр/Тэ;

® задержка между временами вступления Р и Б волн -Та-Тв (рис. 5г);

• длительности Р и 5 волн.

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

параметру - соотношения средних частот спектрограмм волн Р и Б. Преобладающие значения спектров Р-волн промышленных взрывов сосредоточены в районе 6-9 Гц, Б-волн- в районе 2-5 Гц. В среднем преобладающая частота Р-волн в 1.5-2 раза выше, чем частота Б-волн. Указанные различия обусловлены нелинейными процессами распространения волн в земле, характеризуемой неидеальной упругостью.

Метод разделяющих замкнутых разделяющих поверхностей. Рассматриваемый в рамках данной работы класс сейсмических событий (промышленные и полигонные взрывы, падающие на землю тела и др.) характеризуются многообразием волновых форм - откликов среды [10-11]. Такое многообразие характеризуется как параметрами самого события - его спектрально-временной функцией, мощностью источника, а также расстоянием «источник - приёмник», типами волн - сейсмические и акустические, сейсмические Р- и Б-волны, поверхностные волны и др.

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

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

Здесь М -число эталонов, аппроксимирующих данную область; Щ — пороговое значение дляу'-го эталона; т\'\ сг)!) —соответственно математическое ожидание и среднеквадратическое отклонение /-го признакау-го эталона.

Решение о принадлежности к данной области принимается при выполнении одного из условий:

Объем эталона в (10) определяется выбором порогового значения Щ, в

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

» ( - <>) Л2

£ о-=1.2.....ао.

|-| ^ ст у

(10)

(/ = 2.....М).

(П)

то для вычисления соответствующего заданным А' и а, необходимо знать распределение статистики, определяемой левой частью уравнения (12). Данное распределение описывается выражением вида

(l/NXN + l)Г¿_t,

где Т,; , = ^ Л ,_„(аг) есть Т7 — распределение с N-1 степенями свободы;

N - 1 - п

1\ Л1п(а) — известное ^-распределение с п и Ы-1-п степенями свободы, которое для заданного а может быть найдено из таблицы. С учетом сказанного пороговое значение определяется соотношением

Г,* *•«>>. <»>

N N-1-я

Построение эталонов происходит следующим образом. Первая точка х, является центром координат первого эталона — гиперсферы радиуса Л;. Если вторая точка на попадает в первый эталон, то она образует второй эталон в виде гиперсферы того же радиуса. В противном случае корректируются координаты первого эталона. Третья точка определяется на принадлежность одному из двух сформированных эталонов и корректирует параметры х,и т, того из них, в который попадает. В противном случае образуется третий эталон и т. д. Если число точек в каком либо эталоне достигает величины , то в дальнейшем он определяется как гиперэллипсоид с параметрами х, и т, и порогом Л*. Принадлежность к этому эталону очередной точки х[м + ^устанавливается из условия

+ (14)

Если условие (14) выполнено, то корректируются параметры х, и т, данного эталона и т. д.

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

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

Расстояние «источник-приёмник» для удалённых (более сотни километров) сейсмических источников рассчитывается по формуле

8 = У,У,/(УР-У,)АТ, (15)

где Ур, - скорости сейсмических Р и 8 волн соответственно, ЛТ - разность времён вступлений Р и Б волн. Применительно к полигонным источникам: 5 ~(У$УЛ )/(У$ - V. )ЛТ, где Кх, УА - скорости распространения сейсмических и акустических волн соответственно, ЛТ- разность времён вступлений сейсмических и акустических волн.

•л

л.

Отношение частот волн

Длительность первой волны, с

ЩР - Полигонные взрывы (диапазон расстояний: 570-750 м) ф - Карьерные короткозамедленные взрывы (диапазон расстояний: 130-320 км)

Рис. 6. Иллюстрация результата покрытия эталонами областей принятия решения для двух типов событий: полигонных и карьерных промышленных взрывов

Азимут на источник определяется по направлению максимума диаграммы направленности регистрирующей линейки сейсмоприемников, выступающей в качестве приемной антенны. Диаграмма направленности вычисляется путём расчёта коэффициента направленного действия (КНД) сейсмической антенны.

Вычисление коэффициента направленного действия производится по формуле:

К.Н.Д. = Есум/Етм, где (16)

= £ и (О+Л С. + до + ■ • • + Л (/. + Л' • (к -1)))2, Етх = \/к- • £/„2

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

Ei = / R,, (17)

где <7r - стандартное отклонение невязок между вычисленными и истинными координатами событий, Л, - истинное расстояние "источник-приёмник" для каждого г'-го события (вычисленное по GPS).

Для полигонных взрывов с помощью многоканальной обработки и обработки в одноканальном режиме (с использованием вейвлет-фильтрации) погрешность (17) составила не более 3.5%.

Погрешность локации промышленных взрывов (удалённых на несколько сотен километров) составила не более 5%.

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

Программный комплекс спроектирован в двух вариантах, ориентированных на различные применения:

1. Система обработки данных в режиме реального времени - «SeismoDetec-tor». Представляет собой консольное приложение, осуществляющее обмен данными с внешним программным окружением посредством стандартных интерфейсов и протоколов.

2. Интерактивная среда обработки сейсмических данных с графическим интерфейсом пользователя - «SeismoDetector GUI». Включает в себя все разработанные программы обработки данных, которые являются универсальным инструментом оперативной цифровой обработки сейсмических сигналов с визуальным отображением результатов.

Система «SeismoDetector» разработана с учетом возможностей современных вычислительных средств, а также исходя из соображений универсальности взаимодействия с программным окружением POSIX-совместимых ОС:

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

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

® Обеспечена возможность переносимости программной системы на различные ОС и различные аппаратные архитектуры.

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

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

« вейвлет-фильтрации сейсмических данных;

• вычисления огибающей сигнала на основе метода STA/LTA;

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

• определения координат источника (в полярной системе координат) на основе вычисления азимута на источник (коэффициент направленного действия сейсмической антенны) и расстояния «источник-приёмник».

Реализованы следующие сервисные функции:

• подсистема протоколирования процесса обработки сейсмических данных на всех этапах решения задачи локации;

• вывод результатов обнаружения, измерения параметров волн и локации сейсмических источников в БД (MySQL, Postgres SQL, и др.);

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

• визуализации координат источника на цифровой карте местности на основе вывода результатов локации на стандартном языке KML, используемом в целом ряде профессиональных ГИС-систем, таких как: ESRI ArcGIS Explorer, GRASS, MAPinfo, OpenStreetMap, Google Earth и др.

Интерактивная среда обработки сейсмических данных с графическим интерфейсом пользователя - «SeismoDetector GUI» представляет собой графическую оболочку к разработанному комплексу программ и является универсальной средой обработки сейсмических данных. Программа SeismoDetector GUI включает в себя алгоритмы: многоканальной обработки, вейвлет-фильтрации, построения амплитудной огибающей, скользящего среднего, построения функции STA/LTA, обнаружения, распознавания сейсмических событий, вычисления азимутального направления источника, корреляционного анализа данных, интерполяции данных, спектрального анализа, генерирования модельных сейсмических сигналов различного типа и формы.

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

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

• использование человеком-оператором в ручном или полуавтоматическом режиме;

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

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

ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ

Основные результаты, полученные в работе:

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

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

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

3 Выполнен анализ эффективности применения алгоритма вейвлет-фильтрации для выделения волновых форм для последующего высокоточного определения параметров сейсмических волн. Доказана его высокая помехоустойчивость при выделении волновых форм для соотношений амплитуд волн и уровней шумов до 1/8.

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

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

6. Разработано и исследовано программное обеспечение для оперативной обработки широкого класса сейсмических сигналов и шумов в Р051Х-совместимых ОС и графических средств визуализации сейсмических источников. Оценена эффективность выполненной разработки по отношению к задачам обнаружения и распознавания полигонных и карьерных взрывов в режиме реального времени.

Основные результаты диссертации опубликованы в работах [1-11].

ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

1. Хайретдинов М.С, Авроров С.А. Нелинейный метод выделения и измерения параметров сейсмических волн. ГЕО-Сибирь-2006. //Сборник материалов международного научного конгресса «ГЕО-Сибирь-2006». Т.З. Мониторинг окружающей среды, геоэкология, дистанционные методы зондирования Земли и фотограмметрия. Новосибирск: Изд-во СГГА, 2006., Часть 2- С. 8-13.

2. Авроров С.А. Разработка программной системы для поточного обнаружения и распознавания сейсмических сигналов. НАУКА. ТЕХНОЛОГИИ. ИННОВАЦИИ //Материалы всероссийской научной конференции молодых ученых в 7-ми частях. Новосибирск: Изд-во НГТУ, 2006. Часть. 2 - С. 64-66.

3. Авроров С.А. Поточное обнаружение и распознавание сейсмических волн в шумах. ГЕО-Сибирь-2007. //Сборник материалов международного научного конгресса «ГЕО-Сибирь-2007». Т.З. Дистанционные методы зондирования Земли и фотограмметрия, мониторинг окружающей среды, геоэкология. Новосибирск: Изд-во СГГА, 2007,- С. 269-273.

4. Авроров С.А., Хайретдинов М.С. Распределённая обработка данных в иерархических системах и сетях. Новосибирский государственный технический

университет. Научный вестник НГТУ /[редкол.: А.С. Востриков (гл. ред.) и др.].-Новосибирск: Изд-во НГТУ.-№ 2.-2010, С.3-12.

5. Авроров С.А. Исследование характеристик обнаружения и распознавания импульсных сейсмических источников // Труды конференции молодых учёных,- Новосибирск: ИВМиМГ СО РАН, 2009- С. 3-10.

6. Хайретдинов М.С, Авроров С.А. Автоматизированная процедура поточного обнаружения и идентификации сейсмических событий. // Вестник НЯЦ РК, 2008. -С.70-75.

7. Авроров С.А. Программная система сейсмического мониторинга окружающей среды. НАУКА. ТЕХНОЛОГИИ. ИННОВАЦИИ //Материалы всероссийской научной конференции молодых ученых в 7-ми частях. Новосибирск: Изд-во НГТУ, 2008. Часть. 1 - С. 159-162.

8. Авроров С.А. Программная система сейсмического мониторинга окружающей среды. Молодежь и современные информационные технологии. Сборник трудов VI Всероссийской научно-практической конференции студентов, аспирантов и молодых ученых. Томск, 26-28 февраля 2008 г, С.453-454.

9. S.Avrorov, M.Khairetdinov. Stream-Handling Seismic Waves Detection and Recognition System. Proceedings of 2008 IEEE region 8 International conference on computational technologies in electrical and electronics engineering, Novosibirsk, Russia, July 21-25, 2008, pp. 313-317.

10. Khairetdinov, M. Avrorov, S. Computational technology in seismic monitoring system. Proceedings of 2010 IEEE region 8 International conference on computational technologies in electrical and electronics engineering, Irkutsk, Russia, July 1115, 2010, pp. 86-90.

11. Khairetdinov M.S., Avrorov S.A., Livenets A.A. A computational technology in seismic monitoring. // Bulletin of the Novosibirsk Computing Center. Series: Mathematical modeling in geophysics. Is. 13 (2010). NCC Publisher, Novosibirsk, 2010 (in press).

Авроров Сергей Александрович

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

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

Подписано в печать 19.11.2010 г. Формат 60x84/16. Усл. печ. л. 1.5. Уч.-изд. л. 1.5. Тираж 100 экз. Заказ №125

Отпечатано в ООО «Омега Принт» 630090, Новосибирск, пр. Лаврентьева, 6

Введение.

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

1.1 Общие сведения.

1.2 Проблема повышения точности локации сейсмических источников

1.3 Технология калибровки сейсмических трасс с использованием мобильных сейсмических групп.

1.4 Проблема инспекции на месте в рамках Договора о всеобъемлющем запрещении ядерных испытаний.

1.5 Задача определения параметров скважинного источника.

1.6 Задача определения координат падающих ступеней ракет-носителей

1.7 Обзор источников по теме исследования.

1.7 Выводы к главе 1.

2. Задача обнаружения сейсмических сигналов.

2.1 Постановка задачи.

2.1.1 Многоканальная обработка сейсмических данных.

2.1.2 Пространственная фильтрация сейсмических волн.

2.1.3 Математическое моделирование и сравнительный анализ алгоритмов многоканальной обработки.

2.2 Алгоритмы фильтрации импульсных сейсмических сигналов.

2.2.1 Особенности анализа сейсмических волновых форм.

2.2.2 Дискретное вейвлет-преобразование.

2.2.3 Выбор параметров вейвлет-фильтра.

2.2.4 Математическое моделирование процедуры вейвлет-фильтрации в шумах.

2.3 Обнаружение сейсмических событий.

2.3.1 Оптимизация выбора порога обнаружения (Алгоритм «БТА/ЬТА»).

2.3.2 Модельный эксперимент по оценке точности вычисления времён вступления сейсмических сигналов.

2.4 Выводы к главе 2.

3 Распознавание сейсмических событий и локация сейсмоисточника.

3.1 Метод разделяющих замкнутых поверхностей.

3.2 Обнаружение и распознавание полигонных взрывов.

3.3 Обнаружение и распознавание промышленных взрывов на угольных разрезах.

3.4 Вычисление координат источника.

3.4.1 Геоакустическая локация малоудалённых полигонных взрывов.

3.4.2 Решение задачи локации с помощью многоканальной обработки .97 3.4.4 Сейсмическая локация с использованием вейвлет-фильтрации данных.

3.5 Оценивание точности локации промышленных взрывов.

3.6 Помехи обнаружения и распознавания.

3.7 Оценка времени выполнения алгоритмов.

3.8 Выводы к главе 3.

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

4.1 Требования к реализации и общая характеристика программного комплекса.

4.2 Программа обработки сейсмических данных в реальном времени -SeismoDetector.

4.2.1 Общая структура программы SeismoDetector.

4.2.2 Описание функционирования.

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

4.2.4 Подсистема протоколирования.

4.2.5 Особенности реализации алгоритмов обработки сейсмических данных.

4.2.6 Сборка и запуск программы SeismoDetector.

4.2.7 Настройка параметров обработки сейсмических данных.

4.3 Графическая среда обработки сигналов SeismoDetector GUI.

4.4 Модуль сопряжения с картографической системой.

4.5 Выводы к главе 4.

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

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

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

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

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

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

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

В соответствии с целью определены следующие задачи:

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

• разработка алгоритмов обнаружения и измерения параметров сейсмических волн на фоне шумов в поточном режиме. Исследование их помехоустойчивости и точности оценивания параметров;

• разработка алгоритма распознавания сейсмических событий на основе выбора информативных параметров сейсмических и акустических волн;

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

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

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

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

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

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

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

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

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

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

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

• мест падения с воздуха боевых средств поражения;

• эпицентров скрытых подземных ядерных взрывов по авторшоковой эмиссии (инспекция на месте в интересах Договора о всеобъемлющем запрещении ядерных испытаний) и др.

• короткозамедленных мощных (до 100-300 тонн) карьерных взрывов, регистрируемых на расстояниях в сотни километров. Они являются 8 ценными событиями с точки зрения калибровки сейсмических трасс на различных азимутальных направлениях и расстояниях, что, в свою очередь, позволяет повысить точность локации сейсмических событий сетью глобального сейсмического мониторинга.

Реализация и внедрение результатов работы. Результаты диссертации составляют основу программного инструментария, применяющегося в отделе геофизической информатики Института вычислительной математики и математической геофизики (ИВМиИГ СО РАН) при проведении лабораторных исследований и полевых экспедиционных работ.

Диссертационные исследования выполнялись в рамках проектов Российского фонда фундаментальных исследований (РФФИ) №07-07-00214а, 09-07-12075-офи-м и междисциплинарного проекта СО РАН №163.

Исследования алгоритмов обнаружения, измерения параметров и локации тестовых полигонных взрывов выполнялись в рамках хоздоговора "Полесовщик" между НГТУ и ОАО "Завод Ленинец" (г. Санкт-Петербург). Полученные результаты исследований используются в составе технологии автоматизированной локации полигонных взрывов падающих боеприпасов и карьерных взрывов.

Достоверность полученных результатов подтверждается проведёнными натурными экспериментами и математическим моделированием, а также экспертизами отчетов по хоздоговору и грантам РФФИ.

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

-Международной научной конференции «SIBIR.COM» (Новосибирск, 2008 г., Иркутск, 2010 г.).

-Всероссийской научной конференции молодых ученых «Наука. Технологии. Инновации» (г. Новосибирск, 2006, 2008 гг.).

-Конференции молодых учёных ИВМиМГ СО РАН (г. Новосибирск, 2008, 2009 гг.).

-Международных научных конгрессах «ГЕО-Сибирь» (г. Новосибирск, 2007, 2008 гг.).

-Международной конференции «Мониторинг ядерных испытаний и их последствий» (п. Боровое, Казахстан, 2008 г.).

-VI Всероссийской научно-практической конференции студентов, аспирантов и молодых ученых (Томск, 2008 г.).

-Международном семинаре "Акустика неоднородных сред-Х" (Новосибирск, 2009 г.).

-Научной сессии НГТУ (г. Новосибирск, 2008, 2009 гг.).

Публикации. По теме диссертации опубликовано 8 научных статей, а также тезисы доклада. Опубликовано 3 статьи в рецензируемых сборниках, рекомендуемых ВАК.

Структура и объем диссертации. Работа состоит из введения, четырех глав, заключения, списка литературы из 65 наименований и приложения. Основное содержание представлено на 143 страницах, содержит 43 рисунка и 16 таблиц.

Основные результаты, полученные в работе:

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

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

3. Выполнен анализ эффективности применения алгоритма вейвлет-фильтрации для выделения волновых форм для последующего высокоточного определения параметров сейсмических волн. Доказана его высокая помехоустойчивость при выделении волновых форм для соотношений амплитуд волн и уровней шумов до 1/8.

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

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

6. Разработано и исследовано программное обеспечение для оперативной обработки широкого класса сейсмических сигналов и шумов в РСШХ-совместимых ОС и графических средств визуализации сейсмических источников. Оценена эффективность выполненной разработки по отношению к задачам обнаружения и распознавания полигонных и карьерных взрывов в режиме реального времени.

Заключение

1. Science plan for a new global seismographic network // Incorporated Research Institutions for Seismology, 1984.-130 p.

2. А. Шульц, С. Майер, Д. Свенсон, М. Фланаган и др. Интеграция различных калибровочных результатов для улучшения сейсмической локации. Ливерморская национальная лаборатория им.Лоуренса. Контракт № -7405-ENG-48.

3. Активная сейсмология с мощными вибрационными источниками /отв. ред. Г.М. Цибульчик.- Новосибирск: ИВМиМГ СО РАН, Филиал "Гео" издательства СО РАН, 2004. 387 с.

4. Создание технологий калибровки сейсмических станций и сейсмических трасс с использованием мощных сейсмических вибраторов. Отчет по проекту МНТЦ №1067, часть 2. ИВМиМГ СО РАН, Новосибирск, 2001.

5. Григорюк А.П. Распределённая система сбора данных реального времени. // Сб. докл. междунар. научно-технич. конф. «Информационныесистемы и технологии». Новосибирск, Изд. НГТУ, 2000. - Т. 3. - С. 456460.

6. Шорохов М.Н. Вибросейсмическая измерительно-регистрирующая станция. // Сб. докл. междунар. научно-технич. конф. «Информационные системы и технологии». Новосибирск, Изд. НГТУ, 2000. - Т. 1. — С. 144148.

7. Тухватулин Н.Т., Беляшова Н.Н. Деятельность национального ядерного центра в поддержку договора о всеобъемлющем запрещении ядерных испытаний. Вестник НЯЦ РК, 2001, вып.2, С.5-7.

8. Договор о всеобъемлющем запрещении ядерных испытаний (ДВЗЯИ). Вена: Подготовительная комиссия Организации Договора о всеобъемлющем запрещении ядерных испытаний, 1996 г. - 165 с.

9. Пассивная сейсмометрия при проведении инспекции на месте: этап поиска и этап локализации события // Вестник НЯЦ РК, выпуск 2, 2004. -С.112-117.

10. М. S. Khairetdinov, О. К. Omelchenko and others. Numerical algorithms and results of experiments to determine the parameters of the borehole bottom and medium. Bulletin of the Novosibirsk computing center, issue 10 (2005), pp. 127134.

11. Авроров С.А., Хайретдинов M.C. Автоматизированная процедура поточного обнаружения и идентификации сейсмических событий. //Вестник НЯЦ РК, 2008-2. С.70-75.

12. Дробжева Я.В., Краснов В.М., Маслов А.Н. Акустическое поле на земле при взрыве ракеты-носителя //Вестник НЯЦ РК, 2006-2. С.79-85.

13. Краснов В.М., Дробжева Я.В., Пак Г.Д., Салихов Н.М., Козловский В.А., Мухамеджанов Э.К., Маханов Д.И., Лазуркина В.Б. Инфразвуковое излучение при запусках ракет-носителей и его влияние на здоровье населения //Вестник НЯЦ РК, 2008-2. С.60-68.

14. М.Н. Власов, С.В. Кричевский. Экологическая опасность космической деятельности: Аналит. обзор // М: Наука, Центр экологической политики России. 1999.- 238 с.

15. Я.Н. Шойхет, В.Б. Колядо, И.Б. Колядо, С.В. Богданов, С.Н. Дикарев, Г.Я. Евлашевский. Заболеваемость населения территорий, прилегающих к районам падения отделяющихся частей ракет-носителей // Барнаул: Аз Бука. -2005.-200 с.

16. Drobzheva, Ya., The acoustic field in the atmosphere and ionosphere caused by a point explosion on the ground / Ya Drobzheva, V.M. Krasnov // J. Atmos. and Solar-Terr. Phys., 2003. V. 65. Issue 3. P. 369-377.

17. W.Leith, L.J.Kluchko, C.P.Knowles, D.A.Linger, L. Gabriel, Беляшова H.H., Тухватулин Ш.Т., Демин B.H., Коновалов В.Е. Калибровочные испытания ДВЗЯИ на Семипалатинском испытательном полигоне. 1997-2000 гг. //Вестник НЯЦ РК, 2000-2. С.31-35.

18. Хайретдинов М. С., Омельченко О. К., Родинов Ю. И. Автоматизированная технология локации сейсмического источника // Труды международной конференции Математические методы в геофизике. Ч. II. -Новосибирск: Изд. ИВМиМГ СО РАН, 2003. С. 529-535.

19. Freiberger, W. (1963). An approximate method in signal detection, J. Appl. Math. 20, 373-378.

20. Stewart, S. (1977). Real-time detection and location of local seismic events in central California, Bull Seism. Soc. Am. 67, 433-452.

21. Vanderkulk, W., F. Rosen, and S. Lorenz (1965). Large aperture seismic array signal processing study, IBM Final Report, ARPA Contract Number SD-296.

22. Кедров O.K., Поликарпова JI.A., Стеблов Г.М. Алгоритм обнаружения слабых короткопериодных сейсмических сигналов на основе частотно-временного анализа трёхкомпонентных записей в режиме реального времени //Физика земли, 1998. №2. - С.30-45.

23. Хайретдинов М.С., Клименко С.М. Программная система автоматизированной локации и визуализации сейсмических источников.

24. Труды III Межд. конф. "Мониторинг ядерных испытаний и их последствий", Боровое, 2004,-Вестник Ш1Ц РК, С.70-76.

25. Yuri V. Fedorenko, Eystein S. Husebye. First breaks automatic phase pickings of P- and S-onsets in seismic records. Geophysical Research Letters, Vol. 26, No. 21. (November 1999), pp. 3249-3252.

26. Годзиковская A.A. Задача распознавания карьерных взрывов и местных землетрясений // Вопросы инженерной сейсмологии. Вып. 28. М.: Наука. 1987.С.232-235.

27. Кедров O.K., ан ВА., Лаушкин ВА., Люкэ Е.И., Овчинников В.М., Поликарпова JI.A. Методы контроля подземных ядерных взрывов по сейсмическим данным на эпицентральных расстояниях свыше 500 км // Изд. АН СССР. Сер.Физика Земли.М. 1990. № 12. С. 31-46.

28. Savino J.M., Sykes L.R., Lieberman R.S., Molnar P. Exitation of seismic surface waves with periods of 15 to 70 sek for earthquakes and underground explosions. JGR. 1971. V. 76. № 32. P. 8003-8019.

29. Пасечник И.П., Коган С. Д., Султанов Д.Д., Цибульский В.И. Результаты сейсмических наблюдений при подземных ядерных и тротиловых взрывах // М.: Тру-ды ИФЗ АН СССР. 1960. № 15(182). С. 3-52.

30. Кейлис-Борок В.И. Различие спектра поверхностных волн при землетрясениях и подземных взрывах // М.: Труды ИФЗ АН СССР. 1960. № 15(182). С. 88-101.

31. Bennet TJ., Murphy J.R. Analisis of seismic discrimination capabilities using regional data from western United States events. BSSA. 1986. V. 76. № 4. P. 1069-1086.

32. Taylor SR., Sherman N.W., Denny M.D. Spectral discrimination between NTS explosions and western United States earthquakes at regional distance. BSSA. 1988. V. 78. № 4. P. 1563-1579.

33. Liu Н.Н Fu K.S. A syntactic approach to seismic pattern recognition. IEEE Trans. Pattern Analysis Mach. Intell- 1982, PAMI-4, pp. 136-140.

34. Зиневич A.H., Филина А.Г, Еманов А.Ф. Статистический подход к распознаванию промышленных взрывов и землетрясений // Методы изучения, строение и мониторинг литосферы: Материалы к междунар. конф. Новосибирск: Изд-во СО РАН, НИЦ ОИГГМ, 1998. С.205-213.

35. Асминг В.Э. Программный комплекс для автоматизированной обработки сейсмических записей "EL". В сб. Приборы и методика геофизического эксперимента. Мурманск, изд. ООО "МИП-999", 1997, с. 125-132.

36. М.Б.Штарк, Ю.К.Постоенко и д.р. САМАС системы автоматизации в экспериментальной биологии и медицине. Новосибирск, Изд-во «Наука» Сибирское отделение. 1979.- 17 с.

37. M.S. Khairetdinov, O.K. Omelchenko, G.F. Sedukhina, G.M. Voskoboynikova. Numerical algorithms and results of experiments to determine the parameters of the borehole bottom and medium. Bull. Nov. Сотр. Center, Math. Model. In Geoph., 10 (2005), 33-41.

38. A. Kijko, J. Ungvarai. The accuracy of teleseismic events location by a tripartite seismic network in Hungary // Acta Geodaet., Geophys. et Montanist. Acad. Sci. Hung. Tomus 15 (2-4), pp. 239-246.

39. Уидроу Б., Стирнз С. Адаптивная обработка сигналов: Пер. с англ. -М.: Радио и связь, 1989. 440 е.: ил.

40. Тихонов В.И. Оптимальный прием сигналов,- М.: Радио и связь, 1983.-644 с.

41. Гурвич И.И., Боганик Г.Н. Сейсмическая разведка: Учебник для вузов. 3-е изд., перераб. - М.:Недра, 1980, 551 с.

42. Хайретдинов М.С., Григорюк А.П., Дворецкая Л.Г. Частотно-временные характеристики вибросейсмических полей // Методы изучения, строение и мониторинг литосферы: материалы к междунар. конф. Новосибирск: Изд-во СО РАН, НИЦ ОИГГМ, 1998. С. 111-116.

43. Добеши И. Десять лекций по вейвлетам. М.: РХД, 2001. — 464 с.

44. P.J. Onincx. Automatic phase detection in seismic data using the discrete wavelet transform, CWI Report PNA-R9811, October 1998.

45. Donoho D.L. De-noising by soft thresholding // IEEE Trans. Inform. Theory. 1995. 41, N3. P. 613-627.

46. Donoho D.L., Johnstone I.M. Ideal spatial adaptation by wavelet shrinkage//Biometrika. 1994. 81, N3. P.425-455.

47. Андерсон Т. Введение в многомерный статистический анализ. М., Изд-во физ.-мат. лит., 1963.

48. Khairetdinov M.S., Avrorov S.A., Livenets A.A. A computational technology in seismic monitoring. // Bulletin of the Novosibirsk Computing Center. Series: Mathematical modeling in geophysics. Is. 13 (2010). NCC Publisher, Novosibirsk, 2010 (in press).

49. Кишкина С.Б. Параметры сейсмического эффекта массовых короткозамедленных взрывов. Труды III Межд. конф. "Мониторинг ядерных испытаний и их последствий", Боровое, 2004.- Вестник НЯЦ РК, 2004, вып.2, С.171-178.

50. Шоу Уолтон. Создаиие сетевых приложений в среде Linux.: Пер. с англ. М.: Издательский дом "Вильяме", 2001. - 464 с.

51. В.А. Галатепко. Программирование в стандарте POSIX. http://www.intuit.ru/department/se/pposix/. (15 июня 2004 г.).

52. ISO С working group, http://www.open-std.org/jtcl/sc22/wgl4/. October 2010.

53. М. Galassi et al, GNU Scientific Library Reference Manual (3rd Ed.) http://www.gnu.org/software/gsl/. (Октябрь 2010 г.).

54. А. Б. Сергиенко "Цифровая обработка сигналов" (второе издание). СПб, Питер, 2006. С. 302-308.

55. Matteo Frigo and Steven G. Johnson, "The Design and Implementation of FFTW3," Proceedings of the IEEE 93 (2), 216-231 (2005). Invited paper, Special Issue on Program Generation, Optimization, and Platform Adaptation.

56. RFC 3164 The BSD Syslog Protocol. http://www.faqs.org/rfcs/rfc3164.html. (Октябрь 2010 г.).

57. Бланшет Ж., Саммерфилд М. Qt4: программирование GUI на С++. Пер. с англ. 2-е изд., доп.- М.:КУДИЦ-ПРЕСС, 2008.-736 с.64. rsyslog project, http://www.rsyslog.com. (Октябрь 2010 г.).

58. Введение в документацию КМЬ. Электронное справочное руководство, http://code.google.com/intl/ru/apis/kml/documentation/ (Октябрь 2010 г.).

59. ОТКРЫТОЕ ПКЦИОНЕРНОЕ ОБЩЕСТВО НПУЧНВ-ТЕХНИЧЕСКИЙ ЦЕНТР "ЗОВОД ПЕНИНЕЦ"

60. Утверждаю Генеральный директора о. лого На №www.zavod-teninetz.ru e-mailrinfot^onegroup.ru тел. (812)327-90-99 факс (812) 324-61-00

61. Россия, 196084, Санкт-Петербург ул. Коли Томчака, д. 9

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

63. Результаты внедрялись при выполнении совместных работ по хоздоговору «Полесовщик-сейсмо» между ОАО НТЦ «Завод Ленинец» и ГОУ ВПО «Новосибирский государственный технический университет».

64. Председатель комиссии: Члены комиссии:• (A.C. Майков)

65. Учреждение Российской академии наук

66. Ковалевский В.В. зав. лабораторией геофизической информатики ИВМиМГ СО РАН, доктор технических наук.

67. Результаты диссертационной работы были использованы при выполнении работ по грантам РФФИ №№ 07-07-00214а, 09-07-12075-офи-м, 10-07-10000 к, 1007-00387, междисциплинарного проекта СО РАН №163.

68. Просп. Академика Лаврентьева, 6, Новосибирск, 630090 Тел.: (383)330-83-53; Факс (383)330-87-83 E-mail: mikh@sscc.ru ОКПО 03533843, ОГРН 1025403656420, И /т2- /Z.^ZojO № 15301/ ■f—Z.'/SS1. Ha №от1. АКТ

69. Председатель комиссии: Члены комиссии: