Аппаратурно-методические аспекты
детерминистического
подходак прогнозу землетрясений
Д.Г. Таймазов
Институт геологии ДНЦ РАН, старший научный сотрудник
Несмотря на то, что число известных предвестников землетрясений перевалило уже за 300 [1], не выработано еще надежных алгоритмов практического прогноза землетрясений, под которым подразумевается предопределение с приемлемой точностью силы, места и времени будущего толчка, т. е. простое увеличение числа контролируемых предвестников не приводит к адекватному повышению надежности прогноза. Основная причина такой ситуации, которую по упомянутому признаку можно назвать кризисной, заключается в том, что большинство известных предвестников землетрясений, включая сейсмические, геофизические, геохимические и гидрогеологические, несут по сути трансформированную геологической средой вторичную информацию о предвестниковых деформациях земной коры. Она осложнена неоднородностью самой среды, выражающейся в широком разбросе значений, случайном распределении в пространстве и изменчивости во времени коэффициентов трансформации (тензочувствительности) среды, и
фоновыми асейсмическими физико-химическими процессами. Вследствие этого сигналы-предвестники как по величине, так и по знаку неадекватны энергетическим и пространственно-временным характеристикам готовящегося землетрясения, что предопределяет их статистический характер, а также сугубо вероятностный характер основанного на них прогноза.Из вышеизложенного вытекает целесообразность перехода на наблюдения за ограниченным числом детерминированных предвестников, в первую очередь деформационных, находящихся с процессами подготовки землетрясений в неопосредованной причинно-следственной связи. Разумеется, при этом наблюдения должны быть достаточно представительными, т.е. площадными и длиннобазисными.
Главное преимущество деформационных предвестников перед вторичными обусловлено тем, что, согласно теории, между компонентами деформационного поля на дневной поверхности, измеряемыми в 3-х и более азимутах, и компонентами тензора деформаций в объеме охваченного измерениями массива существует взаимно-однозначная связь. Это означает, что любые деформационные процессы, происходящие в потенциальном очаге и связанные с подготовкой землетрясения, приводят к количественно интерпретируемым изменениям деформационного поля в приповерхностном слое земной коры, где располагаются измерительные приборы [2].
Реализация прямых площадных разномасштабных деформационных наблюдений для прогноза землетрясений, необходимость которых обоснована выше, упирается в проблему создания ДС, пригодных по метрологическим и технико-экономическим характеристикам для решения этой задачи, и разработку методики обработки результатов для надежного выделения слабых предвестниковых сигналов на фоне случайных помех на достаточном удалении от источников параметрических аномалий [3].
В предлагаемой работе дано теоретическое и технико-экономическое обоснование одного из возможных направлений решения проблемы деформационного мониторинга земной коры для прогноза землетрясений. Оно основано на ряде предложенных автором аппаратурно-методических разработок, защищенных авторскими свидетельствами. Основные из них вкратце описаны ниже.
1.
Газо-жидкостный гравиметр. В нем роль рабочего тела играет ртутный столб, занимающий нижние части двух широких резервуаров, разнесенных по высоте, и соединяющую их трубку. Верхний резервуар вакуумирован, а нижний заполнен газом. В первом варианте [4] во ртуть в верхнем резервуаре частично погружен цилиндрический груз, подвешенный к струне частотного преобразователя. Во втором [5] индикация уровня ртути осуществляется емкостным преобразователем перемещений.Основные преимущества перед известными аналогами
: наличие эффективной термокомпенсации [6] в сочетании с термостатированием, независимость показаний от малых наклонов, частотный выходной сигнал, высокая степень демпфирования, отсутствие дрейфа нуля. Ожидаемая погрешность ± (1 ё 0,1) мкГал, диапазон измерений 5000 мГал.Область возможного применения: Регистрация приливных и неприливных изменений силы тяжести в стационарных условиях, регистрация вертикальных движений земной коры (совместно с градиентометром), морская гравиметрия, гравиразведка.
2.
Вертикальный гравитационный градиентометр [7]. В нем пробное тело выполнено в виде жестко связанных между собой поплавка и сосуда со ртутью, метацентры которых разнесены по вертикали. Поплавок погружен во ртуть в верхнем сосуде и связан со струнным преобразователем, измеряющим разбаланс между действующими на пробное тело в противоположных направлениях гидростатическими силами, обусловленный вертикальным градиентом силы тяжести.Основные преимущества перед известными аналогами: независимость показаний от малых наклонов, высокая степень демпфирования, большая производительность съемочных работ, частотный выходной сигнал, наличие эффективной термокомпенсации. Ожидаемая погрешность ± 0,2 Е в диапазоне 2Ч 10 4 Е.
Область возможного применения: Регистрация временных вариаций вертикальных градиентов силы тяжести, гравиразведка, регистрация вертикальных движений земной коры (в сочетании с гравиметром), морская гравиметрия, аэрогравиразведка.
3.
Гравитационный вариометр [8]. Для измерения вторых производных гравитационного потенциала в нем используются три разнесенных по координатным осям двухкоординатных наклономера, установленных в одном корпусе. Влияние наклона корпуса устраняется составлением разностей показаний разнесенных по каждой координате датчиков, а влияние неидентичности установки датчиков в корпусе устраняется повторными измерениями с поворотом корпуса на 180° и составлением алгебраических сумм показаний прибора при двух взаимно противоположных ориентациях. Вариометр впервые позволит измерять компоненты Wxx и Wyy раздельно и по уравнению Лапласа вычислять вертикальный градиент Wzz. Ожидаемая погрешность ± (1ё 0,1) Е.Основные преимущества перед известными аналогами: Возможность измерения всех 6 независимых компонентов вторых производных гравипотенциала, большая производительность, частотный выходной сигнал.
Область возможного применения: Регистрация временных вариаций вторых производных гравипотенциала, наземная гравиразведка.
4.
Двухкоординатный наклономер-сейсмограф [9]. В этом варианте наклономера пробное тело выполнено в виде сферического резервуара со ртутью. Для уменьшения требуемого количества ртути и расширения функциональных возможностей в него помещена без механического контакта со стенками жесткая полая сфера, которая связана с чувствительной системой, настроенной на другой динамической диапазон, например сейсмический. Таким образом, устройство выполняет функции наклономера-сейсмографа с широким динамическим диапазоном. Ожидаемая погрешность измерения наклонов ± 5" · 10 -5 в диапазоне до 5" .Основные преимущества перед известными аналогами: неподвижность рабочего тела, отсутствие перекрестного влияния датчиков, высокая степень демпфирования, частотный выходной сигнал, наличие двух чувствительных систем с частично перекрывающимися динамическими диапазонами, большой динамический диапазон (220 дцБел), дистанционная система арретирования, возможность установки в скважине, автоматическое переключение режимов.
Область возможного применения: регистрация приливных и тектонических наклонов земной коры для прогноза землетрясений, регистрация малых горизонтальных ускорений гравитационной природы в прецизионных физических экспериментах, регистрация сейсмоакселерограмм.
5.
Двухкоординатный скважинный деформограф [10]. Он представляет собой цилиндрический корпус из нержавеющей стали, разделенный перегородкой на воспринимающий и измерительный отсеки. В воспринимающий отсек коаксиально ему помещен инварный цилиндр. Зазор между цилиндрами разделен на 3 разноориентированные секции и заполнен рабочей жидкостью (силиконовым маслом). Каждая секция сообщается с отдельным сильфоном в измерительном отсеке, снабженным емкостным преобразователем перемещений его свободного конца. Сопоставляя показания датчиков деформаций трех разноориентированных секций, можно определить величины и направление приращений главных напряжений в исследуемом объекте. Расчетная чувствительность прибора ~ 2Ч 10-12, диапазон измерений ~ ± 10-3 (без перестройки шкалы).Основные преимущества перед известными аналогами: возможность измерений как объемных, так и линейных деформаций в двух ортогональных азимутах, наличие эффективной термокомпенсации, малые габариты, цифровой выходной сигнал.
Область возможного применения: регистрация объемных и линейных деформаций земной коры для прогноза землетрясений, контроль за состоянием крупных инженерных сооружений.
6.
Емкостный преобразователь перемещений для геофизических приборов [11]. Он относится к типу преобразователей с изменяемой площадью перекрытия статорных и роторных обкладок, которые выполнены в виде множества параллельно соединенных поперечных полосок с промежутками, причем для наиболее полного использования преимуществ дифференциального метода измерений предусмотрена возможность автоматического подбора числа полосок каждого ряда, включаемых в плечи дифференциальной схемы исходя из требования максимальной близости их емкостей. По соотношению чисел этих полосок осуществляется грубый отсчет значения фазы перекрытия статорных полосок роторными, а по измеряемой дифференциальной схемой разности емкостей - тонкий отсчет. Расчетная погрешность ± 5Ч 10-5 мкм в диапазоне ± 100 мм, что позволяет отказаться от всех других преобразователей, используемых в геофизических приборах. При измерении угловых разворотов, например в теодолитах, соответствующая погрешность преобразователя при диаметре лимба 10 см составит ± 1" · 10-5 против ± 1" в известных кодовых теодолитах. Конструктивно теодолит с вышеописанным емкостным преобразователем несравненно проще существующих кодовых теодолитов.Из других областей использования преобразователя следует назвать гравиметры, наклономеры, крутильные весы, гравиградиентометры, гидронивелиры, экстензометры, пьезометры, установки для эталонирования гравиметров наклоном, астрометрические приборы, прецизионные станки и т.п..
7.
Проволочный экстензометр [12]. В нем между постаментами натянуты три одинаковых параллельных проволоки, свободные концы которых перекинуты через блоки. Один из блоков установлен на первом постаменте, а два других - на втором постаменте. Смежные концы, соответственно, первой - второй и второй - третьей проволок попарно скреплены между собой, а свободный конец третьей проволоки проходит через регистрирующее устройство, например емкостный датчик перемещений. Специальным подбором весов грузов, подвешенных к концам проволок, обеспечивается полная компенсация их крипа, а подбором объемов грузов – компенсация барических изменений длины проволоки L за счет всестороннего сжатия – расширения. Для учета изменений L за счет термического расширения проволоки параллельно, на тех же постаментах, устанавливается второй экстензометр с проволокой из материала с другим коэффициентом термического расширения, например из графита.Таким образом, в предлагаемой системе ПЭ все основные инструментальные погрешности поддаются или прецизионной компенсации, или аналитическому учету. При
L = 100 м относительная погрешность измерений составит ± 5Ч 10-13 в диапазоне ± 10-3.8.
Гидростатический нивелир [13]. Он представляет собой замкнутую герметизированную систему, состоящую из двух цилиндрических гидростатических сосудов (ГС) и двух эластичных трубок, соединяющих соответственно нижние и верхние части ГС. Нижняя трубка и нижние части ГС заполнены тяжелой жидкостью, а верхняя трубка и верхние части ГС - более легкой непроводящей жидкостью, не смешивающейся с тяжелой. Путем подбора параметров, в предлагаемом ГН практически устранены влияния градиентов температуры по трассе нивелирования, колебаний интегральной температуры измерительной системы и перепадов атмосферного давления между станциями (за счет герметичности системы и отсутствия в ней свободных объемов), локальных изменений температуры ГС и их наклонов (за счет специального расположения соединительных трубок вблизи ГС и подвешивания статорных обкладок к корпусу), гидродинамических помех (за счет автоматической установки уровней) и погрешностей индикации уровней (за счет использования предложенного нами емкостного преобразователя перемещений).9.
Глубинный геодезический репер [14]. Он представляет собой полую инварную трубу нулевой плавучести, закрепленную в забое гидроизолированной водонаполненной скважины и проходящую вблизи ее осевой линии без механического контакта со стенками. Благодаря этому любые поперечные гравиинерциальные силы, действующие на репер, уравновешиваются противонаправленными гидростатическими (архимедовыми) силами, т.е. на репер не действует никаких изгибающих усилий. Верхний конец репера, выходящий в траншею с деформографами и гидронивелирами или на дневную поверхность (для геодезических наблюдений), в реальном масштабе повторяет движения породы в месте закрепления в скважине его нижнего конца. Кроме того, репер может быть использован как эталон длины вертикального штангового экстензометра, для чего достаточно снабдить его прецизионным измерителем вертикальных движений относительно обсадной трубы скважины, например емкостным. Таким образом, использование предлагаемого глубинного репера устраняет метеовоздействия на базы и постаменты измерительных приборов, избавляя от необходимости размещения их в глубоких подземных выработках, и тем самым устраняет также искажения деформационного поля за счет эффекта полости, т.е. существенно улучшают как метрологические, так и экономические характеристики ДС.В целом ДС с описанным выше гравиинерциальным комплексом будет представлять собой проложенные в направлениях С-Ю и В-З две субгоризонтальные пересекающиеся траншеи с перекрытием длиной порядка 50-100 м и сечением около 2,5
ґ 1 м 2, в которых установлены проволочные экстензометры и гидростатические нивелиры (рис. 1). Для повышения информативности дополнительно может быть проложена и оборудована диагональная траншея. На пересечениях траншей пробурены три субвертикальные скважины глубиной 100-200 м (до скальных пород), в которых установлены двухкоординатные скважинные деформографы [10], наклономеры-сейсмографы [9] и глубинные геодезические репера [14], верхние концы которых выходят в траншею (к экстензометрам и гидронивелирам) или на дневную поверхность (для геодезических измерений) и в реальном масштабе повторяют движения пород в местах закрепления в скважинах их нижних концов. Одновременно репер, представляющий собой инварную трубу, служит эталоном длины вертикального штангового экстензометра, для чего он снабжен емкостным преобразователем вертикальных перемещений относительно обсадной трубы скважины.На одном из пересечений траншей располагается лабораторное помещение станции, где устанавливаются гравиметр, гравиградиентометр, вариометр и аппаратура для предварительной обработки и передачи информации в центральный пункт системы прогнозных наблюдений.
По предварительным расчетам предлагаемая ДС в 20
ё 50 раз дешевле известного аналога штольневого типа, что делает реальным создание в сейсмоопасных районах сети ДС с интервалами порядка 20-30 км . Создание такой сети существенно расширит круг контролируемых параметров земной коры, в число которых войдут [15]: 1) локальные объемные деформации пород; 2) локальные линейные деформации по двум горизонтальным координатам; 3) градиенты линейных и объемных деформаций (по разнесенным в плане скважинным деформографам);
4) локальные наклоны по двум ортогональным азимутам;) сейсмоакселерограммы по двум горизонтальным координатам в динамическом диапазоне 120 децибел; 6) длиннобазовые (
~ 100 м) линейные деформации по трем ортогональным координатам; 7) градиенты среднемасштабных линейных деформаций (по разнесенным в плане ДС); 8) среднемасштабные наклоны земной коры в двух ортогональных азимутах (гидронивелиры); 9) крупномасштабные горизонтальные движения земной коры (по геодезическим измерениям с использованием глубинных реперов); 10) вертикальные движения земной коры (по гравиметрическим наблюдениям); 11)неприливные изменения силы тяжести; 12) амплитуды и фазы приливных гармоник во временных изменениях деформаций, наклонов и силы тяжести. Становится возможным регулярное разномасштабное картирование всех перечисленных параметров, что означает качественно новый уровень контроля за динамикой земной коры.Для наиболее полного выделения из результатов площадных наблюдений информации, детерминированной процессами подготовки землетрясений, мы считаем целесообразным такой подход к обработке данных, при котором используются только наиболее общие признаки, присущие разным моделям подготовки землетрясений. Такими признаками являются следующие: 1) энергия потенциального очага землетрясений заключена в упругой составляющей накопленных в нем сдвиговых деформаций ; 2) величина
и скорость ее убывания с удалением от эпицентра определяются магнитудой готовящегося землетрясения и глубиной залегания очага; 3)
в потенциальном очаге землетрясений наращиваются скачкообразно и изолинии амплитуд этих скачков оконтуривают потенциальный очаг [16]. Следовательно, эти изолинии содержат статистически значимые кольцевые составляющие, параметры которых определяются магнитудой и глубиной очага [17]. Таким образом, обратная задача (задача прогноза) сводится к поиску математическими методами кольцевых структур в изолиниях распределения скачков сдвиговых деформаций, которые согласно [16] следуют за подвижками на разломе и геометрические центры которых совпадают с эпицентрами потенциальных очагов. Это может быть осуществлено, например, по методике, предложенной в работе [18]. По сечениям изолиний в кольцевых структурах можно оценить магнитуды будущих толчков и глубины очагов уже на докритической стадии их подготовки, соответствующей этапу среднесрочного прогноза. Определенность координат очаговой области потенциального главного толчка, где располагаются и очаги форшоков, делает возможным распознавание последних в реальном времени, что существенно повысит надежность краткосрочного прогноза сильных землетрясений по форшокам, в том числе и одиночным.
Поскольку организация непрерывных площадных деформационных наблюдений – трудная научно-техническая задача, то в предлагаемом методе средне- и краткосрочного прогноза землетрясений могут быть использованы и данные площадных наблюдений эндогенной составляющей микросейсмической эмиссии и/или скоростных характеристик среды, которые также характеризуют ее напряженно-деформированное состояние.
Итак, первый шаг, предлагаемый нами в направлении повышения прогностической информативности деформационного мониторинга земной коры, состоит в повышении чувствительности и помехозащищенности используемых для этого приборов (гравиметров, гравиградиентометров, вариометров, наклономеров, гидронивелиров, деформографов). Это достигается за счет гидростатического усиления измеряемых гравиинерциальных сил (гравиметр, вертикальный градиентометр, наклономер, микроакселерометр), прецизионной компенсации влияния внешних факторов (температурных, барических и т.п.) адекватным и противонаправленным воздействием тех же факторов (гравиметр, градиентометр, скважинный деформограф, гидронивелир, проволочный экстензометр), использования слабонагруженных струнных преобразователей (гравиметр, градиентометр, наклономер, акселерометр) и прецизионных емкостных преобразователей (гравиметр, деформограф, экстензометр, гидронивелир, теодолит), уменьшения влияния наклонов (гравиметр, градиентометр, вариометр, гидронивелир), подавления динамических помех жидкостным демпфированием (гравиметр, градиентометр, наклономер, акселерометр, вариометр).
Второй шаг заключается в использовании глубинных реперов, практически устраняющих влияние метеопомех на базы и постаменты измерительных приборов (теодолитов, светодальномеров, экстензометров, гидронивелиров) и существенно уменьшающих стоимость ДС.
Улучшение метрологических характеристик приборов и подавление влияния метеопомех на их базы и постаменты создают предпосылки для третьего шага - удаления ДС от источников параметрических аномалий, например от разломных зон во внутренние участки геоблоков, где деформационные поля имеют более регулярный характер и адекватно отражают динамику накопления упругих деформаций в их объеме (в этих деформациях и заключена вся тектоническая энергия, поскольку она может аккумулироваться только в консолидированной среде). Уменьшение стоимости ДС открывает возможность площадных наблюдений деформационного поля с интервалами между станциями порядка 20-30 км.
И, наконец, четвертый, завершающий шаг, предлагаемый нами как один из возможных, состоит в выделении математическими методами из данных площадных наблюдений компонентов деформационного поля статистически значимых кольцевых составляющих, которые по предложенной в [16] модели подготовки очага должны отражать скачкообразный характер накопления упругих деформаций в очаговой зоне. По параметрам кольцевых структур возможно определение магнитуды и координат очаговой зоны землетрясения на докритической стадии его подготовки, соответствующей этапу среднесрочного прогноза, а по форшоковой и микросейсмической активности этой зоны - определение времени главного толчка.
В соответствии с вышеизложенным и предложенной моделью подготовки очага стратегия прогноза землетрясений должна включать следующие основные этапы.
1. Выявление активных разломов и оценка их тектонических потенциалов (детальное сейсморайонирование).
2. Выявление на активных разломах участков с наименьшими скоростями контрастных движений, соответствующих в принятой модели зонам максимальной концентрации тектонической энергии взаимодействия блоков (долгосрочный прогноз).
3. Организация в окрестностях выявленных участков площадных деформационных наблюдений с максимально возможным охватом внутренних (удаленных от разлома) участков взаимодействующих геоблоков.
4.
Математическая обработка площадных наблюдений для выделения в них кольцевых структур и определения по ним магнитуды и координат очаговой зоны готовящегося землетрясения (среднесрочный прогноз).5. Определение времени главного толчка по форшоковой и микросейсмической активности очаговой зоны с использованием также других предвестников, как деформационных, так и вторичных (краткосрочный прогноз).
СПИСОК ЛИТЕРАТУРЫ
1. Основы прогноза землетрясений. Душанбе-Москва: Дониш, 1984. 216 с.
2. Латынина Л.А., Кармалеева Р.М. Деформографические измерения. М.: Наука,
1978. 154 с.
3. Соболев Г.А. Основы прогноза землетрясений. М.: Наука, 1993. 313 с.
4. Таймазов Д.Г. А.С. СССР № 1251701. Гравиметр Таймазова. 1986.
5. Таймазов Д.Г. А.С. СССР № 1241885. Гравиметр Д.Г. Таймазова. 1986.
6. Таймазов Д.Г. А.С. СССР № 1153692.Термокомпенсатор для газовых
гравиметров. 1985.
7. Таймазов Д.Г. А.С. СССР № 1238565. Вертикальный гравитационный
градиентометр.1986.
8. Таймазов Д.Г. А.С. СССР № 1373202. Способ измерения вторых производных
гравитационного потенциала. 1987.
9. Таймазов Д.Г. А.С. СССР № 1242713. Наклономер Д.Г. Таймазова. 1986.
10. Таймазов Д.Г. А.С. СССР № 1640544. Двухкоординатный деформограф
Таймазова. 1989.
11. Таймазов Д.Г. Емкостный преобразователь перемещений для геофизических
приборов. Тез. докл. межд. симпоз. “Мониторинг и прогнозирование
чрезвычайных ситуаций”. Махачкала, 1997. С. 127-128.
12. Таймазов Д. Г. Проволочный экстензометр для сейсмопрогностических
наблюдений. Там же, с. 107-108.
13. Таймазов Д. Г. Прецизионный гидростатический нивелир для
сейсмопрогностических наблюдений. Там же, с. 112-114.
14. Таймазов Д. Г. О путях повышения эффективности сейсмопрогностических
наблюдений в Дагестане.- В кн.: Наука и социальный прогресс Дагестана.
Махачкала, 1997. С. 50-55.
15. Таймазов Д. Г. Состояние и перспективы сейсмопрогностических
исследований в Дагестане. – в кн.: Будущее Дагестана. Махачкала,
1994. С. 178-181.
16. Таймазов Д. Г. Наведенная сейсмичность и тектоника плит. –Там же,
С. 278-282.
17. Таймазов Д. Г. О возможности прогноза землетрясений по скачкообразным
изменениям деформационного поля. Тез. докл. межд. симпоз. “Мониторинг и
прогнозирование чрезвычайных ситуаций”. Махачкала, 1997. С. 25-26.
18. Овчаренко А. В. ДАН. Т.342. №4, 1995. С.537-539.
© Д.Г. Таймазов