Моделирование электромагнитных явлений и программное обеспечение при геодинамических исследованиях

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

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

Как известно, почти вся территория Кыргызского Тянь-Шаня расположена в 8-9 балльных зонах, с древнейших времен здесь по различным источникам отмечены многочисленные землетрясения, иногда катастрофические по своим последствиям (Джанузаков и др., 2003). Только за последние сто лет в пределах Кыргызского Тянь-Шаня и в приграничных районах произошли такие катастрофические землетрясения, как 9-10- балльное Беловодское 1885 г., 9-10-балльное Верненское 1887 г., 10- балльное Чиликское 1889г., 10-11 балльное Кеминское 1911г., 8-9 балльное Кемино - Чуйское 1938г., 9-10 балльное Чаткальское 1946г., 8-9 балльное Сарыкамышское 1970г., 8- балльное Исфара - Баткенское 1977г., 8-9- балльное Жаланаш - Тюпское 1978г., 8- балльное Дараут- Курганское 1978, 8- балльное Кочкор- Атинское 1991г., 8-9 балльное Суусамырское 1992г. и большое количество более слабых подземных толчков. Очаги землетрясений Тянь- Шаня залегают в земной коре на глубинах 5-20 км, максимальное их число залегает на глубинах 12-15км. Очаги с глубиной 20-25 км выявлены в Южно-Чуйском, Северо-Ферганском участках и Южно-Ферганской зоне, 25-30 км отмечены в Кемино-Чиликском и Фергано-Чаткальском сейсмоактивных участках. Глубже 30 км очаги встречаются редко.

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

Северо-Тянь-Шаньская зона включает в себя северный склон Кыргызского хребта, полосу адыров (предгорий) и узкую полосу равнинной части Чуйской впадины, начиная от западной границы республики и до Боомского отрезка Чуйской долины. Эта зона связана с важнейшими разломами, которые по своему характеру являются краевыми и по которым граничат три разнородных по своему строению и режиму структурных элемента. Район проведения геомагнитных и электромагнитных исследований в диапазоне сверхнизких частот по изучению вариаций удельного электрического сопротивления горных пород в координатах 4230' ? 4250' по широте и 7400 '? 7500' по долготе, располагается в пределах Северо-Тяньшанской сейсмогенной зоне.

Геомагнитные ?Z-предвестники были зарегистрированы на расстояния от 50 до 300 км, амплитуды предвестников перед землетрясениями составляли: 14 нТл (Алма-Атинский полигон - Жаланаш-Тюпское 1978 г.)], 130 и 135 нТл (Фрунзенский геофизический полигон - Кочкорское 1974 г.)]. Длительность геомагнитного предвестника ?Z составляла от 13 суток до более чем один год.

Перед Суусамырским землетрясением 19 августа 1992 г. Кр=17,0, Н=25 км ц_N = 42,07⁰, л_E =73,63⁰, силой в эпицентре 910 баллов, землетрясение произошло без форшокового проявления и только после события наступил «взрыв афтершоков». Все сейсмические события с К ? 11,6 с 1986 по 1992 г. до Суусамырского землетрясения были самостоятельными землетрясениями или были форшоками и афтершоками других событий. На отдельных пунктах наблюдений перед Суусамырским землетрясением 19 августа 1992 г. (К=17) происходили уменьшения значений ДТ с 2 нТл до 13 нТл, на других повышение значений ДТ до 3 - 6 нТл в течение 57 месяцев или длительные трендообразные снижения значений ДТ от 2 нТл до 4 нТл.

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

В результате многолетних исследований среднегодовых изменений составляющих H, Z геомагнитного поля и полного вектора Т геомагнитного поля проводилось по данным 12 пунктов. Результаты анализа среднегодовых значений компонент геомагнитного поля показывают: по Н-составляющей за период с 1989 по 2007 гг. по всем пунктам характеризуется понижением среднегодовых значений, особенно в Северо-Тяньшаньской сейсмогенной зоне; по Z-составляющей отмечается линейное повышение среднегодовых значений (рис. 13,14).

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

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

Рис.13. Графики среднегодовых значений Н, Z- составляющих геомагнитного поля.

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

Рис. 14. А - графики расчётов среднегодовых значений полного вектора Т геомагнитного поля. Б - графики расчётов среднегодовых значений полного вектора Т геомагнитного поля.

По всем выше перечисленным геомагнитным пунктам расположенным на территории Кыргызстана и Казахстана с 1995 по 2005 гг. отмечается линейное повышение среднегодовых значений полного вектора Т геомагнитного поля.

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

Анализ исследования геомагнитного поля показывает, что основную территорию Кыргызского Тянь-Шаня занимают тела, создающие положительные аномалии магнитного поля. Локальные высокоинтенсивные отрицательные аномалии отмечаются в виде остатков небольших размеров, которые по видимому связаны с высокотемпературными участками земной коры, их динамикой и сейсмичностью. Анализ аэромагнитной съемки проведенные в различные годы, позволяет выделить многочисленные участки геомагнитных аномалий, а также крупные участки геомагнитных аномалий, а также крупные участки геомагнитных аномалий контактово- метасоматического характера (рис.15).

Рис.15. Карты аномального геомагнитного поля Кыргызского Тянь-Шаня для высот: Н=60 км; Н=100 км, в координатах л = 70 - 80⁰ в.д. и ц = 39 - 40⁰ с.ш. Сплошные изолинии - положительные значения, пунктирные отрицательные поля Т, мЭ. (Карты составлены на основе измерений спутником Земли).

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

Серединным массивам соответствуют дифференцированные геомагнитные поля преимущественно положительного знака, складчатым областям - в основном, линейно вытянутые отрицательные. На фоне региональных отрицательных и положительных аномалий выделяются локальные отрицательные кольцевые аномалии интенсивные по напряженности. Аномальные изменения модуля полного вектора напряженности магнитного поля Земли, были зарегистрированы перед землетрясениями различной энергии. Амплитуды Т- предвестников составляли единицы, реже 20 нТл ( Алайское М= 6,8, 1978 г; Джигатальское М= 5,9, 1984г; Сууамырское, М =7,3 1992г;)(Мавлянов и др. 1978; Сковородкин, 1980; Ким, 1981; Муминов и др. 1984; Шакиров и др., 1986 ; Аббдулабеков, 1989; Курскеев, 1990; Максудов, 1996; Бакиров К.Б и др.,2006). Геомагнитные ДZ - предвестники были зарегистрированы на расстояния от 50 до 300 км, амплитуды предветсников перед землетрясениями составляли: 14нТл (Алма - Атинский полигон - Жаланаш - Тюпское 1978г) ( Ержанов и др ., 1979), 130 и 135 нТл ( Фрузенский геофизический полигон - Кочкорское 1974 г) ( Шакиров и др., 1978). Сильные вариации геомагнитного поля наблюдались за несколько дней до землетрясения в Заалайском хребте 2 ноября 1978 г., аномалии порядка несколько нТл фиксировались в зоне эпицентра Заалайского землетрясения.

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

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

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

По профилю Оттук-Конкино, между перевалом Кескен - Бел и п. г.т Боконбаево, отмечаются четыре разлома вертикального характера, интенсивность магнитного поля Т варьирует от 180 до 250 нТл (рис.16 ).

Рис. 16. Геомагнитный профиль Оттук-Конкино.

А: 1 - разностное поле ДТ относительно геомагнитной обсерватории «Таш-Мойнок»; 2 - МПЗ для высоты «0» без учета влияния неоднородностей земной коры по модульным и компонентным измерениям спутником «МАГСАТ».

Б: 1 - кремнистые породы; 2 - крупно- и среднезернистые граниты; 3 - граниты и гранодиориты; 4 - граниты, гнейсовые порфировидные гранодиориты тоналиты.

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



Рис. 17. Геомагнитный профиль Балыкчи-Кен-Суу.

А: 1 - разностное поле ДТ относительно геомагнитной обсерватории «Таш-Мойнок»; 2 - МПЗ для высоты «0» без учета влияния неоднородностей земной коры по модульным и компонентным измерениям спутником «МАГСАТ». Б: 1 - песчаники, известняки, сланцы; 2 - кремнистые породы; 3 - граниты и гранодиориты.

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

По наблюдениям повторных наблюдений установлено, что между 1987 и 1990 гг., на данном объекте наблюдались понижения поля в пределах от 10 до 12 нТл в период подготовки Чолпон - Атинского землетрясения с К= 12,9.

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

Проведенными экспериментальными циклами повторных измерений геомагнитного поля на 22 пунктах в пределах Бишкекского и Иссык-Кульского геодинамических полигонов, установлено, что степень магнитной активности в разных пунктах была различной (рис. 18).

А. ПК на северном берегу Иссык-Куля

Б. ПК на южном берегу Иссык-Куля

В. Чуйская впадина и ее горное обрамление

Рис. 18. Графики разностных значений по пунктам повторных наблюдений геомагнитного поля относительно базовой геомагнитной обсерватории Таш-Мойнок с 1986 по 1992 г.

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

Таблица 1 - Землетрясения с К ? 11,6 произошедшие с 1986 по 1999 г. в координатах ц_N = 4244⁰, л_E =7376.⁰

Дата, г д м	t0 ч м сек	, N	, E	Kp	h,км	Msr	Примечание
1988 13 03	13 48 26.9	42,10	75,47	12,6	7	4,8	Кочкорское 6 б
1989 05 03	13 48 41.8	42,57	74,73	11,7	5	4,3
1989 08 03	12 43 40.0	42,57	74,73	12,0		4,3	Таш-Башатское

Таблица 2 - Землетрясения с К ? 11,6 произошедшие с 1986 по 1990 г. в координатах ц_N = 4244⁰, л_E =73,3082⁰.

№	Дата, г д м	t0 ч м сек	, N	, E	Kp	h, км	Msr	Примечание
1	1986 25 01	16 00 41.2	43,10	77,38	12,6	15	4,8
2	1986 14 02	01 52 01.0	44,00	78,20	12,7	20	4,8	Шалакское
3	1986 10 05	12 47 41.0	43,20	78,20	12,4	10	4,7
4	1986 18 06	04 30 07.0	43,00	79,80	11,5		4,3
5	1987 05 01	22 52 46.2	42,20	81,20	13,7	18	5,7
6	1987 26 05	13 44 20.8	42,87	77,93	11,8	10	4,3
7	1987 20 09	03 54 04.6	42,90	77,60	11,9	10	4,7
8	1988 13 03	13 48 26.9	42,10	75,47	12,6	7	4,8	Кочкорское, 6 б
9	198817 06	13 30 43.6	42,93	77,40	12,9	21	5,0	Чолпон-Атинское, 6-7б
10	1989 05 03	13 48 41.8	42,57	74,73	11,7	5	4,3
11	1989 08 03	12 43 40.0	42,57	74,73	12,0		4,3
12	1990 02 02	14 04 24.4	42,28	76,48	11,8	5	4,3
13	1990 03 05	10 02 19.4	42,92	76,83	12,2	1015	4,5
14	1990 12 11	12 28 51.4	42,93	77,93	15,3	15	6,3	Бойсоорунское, 8 б
15	1990 28 12	00 26 14.4	42,98	77,95	11,7	15-20	4,3

Моделирование, исследование изменения магнитного поля в связи со взрывом на участке Тору -Айгыр (Ким и др.,1981), анализ результатов сравнения измерений величин поля Т на ряде точек до и после взрывов разной мощности исследователи получили изменения поля интенсивностью от 5 до 10 нТл.

Непосредственно при участии автора данной работы были проведены экспериментальные исследования с применением геофизического аппаратурного комплекса (магнитометры типа М-33, ММП-203) по изучению изменения геомагнитного поля до и после заполнения Токтогульского водохранилища (Шакировым, Погребным, Бакировым К.Б и др., 1986).Моделирование и изучение модуля полного вектора Т показало, что локальные изменения магнитного поля под действием нагрузки - разгрузки составили ~ 6- 10 н Тл.

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

Исследования проводились в разные годы, с ежегодным дополнением новых данных, по восьми крупным разрезам в пределах Бишкекского и Иссык-Кульского геодинамических полигонов. Данные Аламединского разреза показали, что в неоднородных массивах могут освобождаться различные по величине энергии (рис.19, 20, 21, 22).

Иссык - Атинский разрез по скоростным и плотностным свойствам более или менее квазиоднороден. Сосновский разрез имеет плотность 2,7 г/см³. Скорость продольной волны изменяется от 5,2 до 6,1 км/с, а поперечной - от 3,1...3,5 км/с. Результаты исследований позволяют предположить, что глубинные горизонты интрузивных тел поддаются детальному расчленению и корреляции.

Терскейский магматический комплекс пород характеризуется плотностью 2,65...2,7 г/см³ .. Кунгейская зона. В пределах всей длины интрузивного комплекса плотность пород изменяется от 2,55...2,57.. до 2,69...2,67..2,66 г/см³ , Vp - 4,8...5,1...6,2...4,4 км/с, Vs - 3,1...2,9 км/с.

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



Рис. 19. Схематическая карта скоростных свойств массивов горных пород Северного Тянь-Шаня (Бакиров К.Б., 2005):1- геологический возраст интрузивных пород; 2 - значение скорости поперечной волны км/с;3 - значение скорости продольной волны км/с.

Рис. 20. Схематическая карта магнитных и плотностных свойств массивов горных пород Северного Тянь-Шаня (Бакиров К.Б., 2005). 1 - направление склонения вектора остаточной намагниченности; 2 - значение плотности 1.03 г/см3; 3 - геологический возраст интрузивных пород.



Рис. 21. Схематическая карта прочностных свойств массивов горных пород Северного Тянь-Шаня (Бакиров К.Б., 2005).1 - песчаники; 2 - эффузивы; 3 - известняки; 4, 6, 7 - граниты; 8 - диориты; 9 - значение предела прочности пород (МПа); 10 - геологический возраст интрузивных пород.

Рис. 22. Схематическая карта отрицательных магнитных аномалий

Северного Тянь-Шаня (Бакиров К.Б.). 1 - отрицательные аномалии МПЗ (нТл), 2 - геологический возраст интрузивных пород.

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

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

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

Анализ многолетних исследований геомагнитных материалов в Северном Тянь - Шане позволил выявить кольцевые магнитные аномалии на 24 участках. В пределах Кыргызского хребта выделен отрицательный аномальный участок, интенсивность аномалии не превышает 200 нТл. На северном крыле профиля обнажаются андезитовые порфириты, туфы, андезито - базальтовые порфириты, которые отмечаются положительным полем интенсивностью не более 500 нТл. Иссык - Атинское меридиональное отрицательное поле глубинного характера интенсивностью 200 нТл совпадает с выходами гнейсовидных адамеллитов гранодиоритов. На северном фланге этой аномалии, выделяются две отрицательные аномалии интенсивностью 360- 600 нТл. И в районе Кемино- Актюзского рудного поля - пять локальных кольцевых аномалий. Интенсивность поля варьирует от 400 до 700 нТл.

Между реками Каракунгей и Турасу оконтурены в виде цепочки три кольцевые аномалии, интенсивность которых варьирует от 200 до 400 нТл. В сейсмическом отношении район характеризуется проявлением землетрясений с К= 10 - 13. Между перевалами Тон и Барскоон в приосевой части Терскейского хребта цепочкой расположены с запада на восток пять отрицательных кольцевых аномалий. Вокруг кольцевых аномалий, происходили землетрясения с К = 15,8 12,9 и 10.Южнее п. Теплоключенка и Нововознесеновка на северном склоне Терскей Ала -Тоо оконтурена обширная змееобразная по форме зона отрицательных аномалий. Вокруг отрицательных кольцевых аномалий происходили землетрясения с К = 16, а вокруг аномалии землетрясения с К= 11, 12 и 10. Интенсивность - от 400 до 500 нТл. В восточной части этой аномалии между реками Талды - Суу и Кен - Суу 24 км марта 1978 г. произошло сильное Жаланаш - Тюпское землетрясение (л- 78,90, ц -42,90) К= 15,6, М =6,6, н = 20 км. Ближе к середине этой аномалии - 12 ноября 1990г. (через 12 лет) произошло Бойсоорунское землетрясение (ц- 42,96, л - 78,05) К = 15,3, М= 6,3, Н= 3 км.

Зона отрицательной магнитной аномалии в своих краевых частях характеризуется различным проявлением сейсмичности. Северная и западная ее части более сейсмически активны, чем южная. Только с 1989 по 2004 г. на юге в районе пос. Чон - Орукты произошло более 40 землетрясений. Район г. Чолпон - Ата также находится в зоне отрицательной магнитной аномалии. Здесь она связана с выходами высокотемпературных источников. К северу от г. Чолпон - Ата на обоих бортах р. Чон - Кемин расположены две небольшие кольцевые аномалии, где произошло более 15 землетрясений К? 10.

На основании анализа по данным аэромагнитных и наземных съемок в Северном Тянь-Шане установлено, что интрузии развиты в виде отдельных изолированных участков, неправильных в плане тел .В Суусамырском, Терксейском и Кунгейском хребтах они образуют мощные, цепочкообразные, широтно- вытянутые поля протяженностью более 600км. К ним относятся Аксу - Коккийская аномальная зона, интенсивность аномалий варьирует от 100 до 400 нТл, отмечены Джетыогуз-Теримторская. Бурхан- Арабельская . Каракаманская. Капкаташская. Каракуджурская. Куганидская, Калмакашуйская, Акташская Западно - Сонкельская, Суусамыр - Каракольская и Арамсуйские магнитные аномалии.

В контактовых частях гранитоидной интрузии Северного Тянь - Шаня установлено 15 крупных магнитных аномалий контактового типа.

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

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

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

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

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

Анализ экспериментальных измерений аудиомагнитотеллурического (АМТ) мониторинга показывают, что наиболее высокие кажущиеся сопротивления отмечены на участках: Кызыл-Туу, расположенного на среднечетвертичных пролювиально - аллювиальных отложениях 300 - 700 Ом м, Кызыл - Бирдик 80-300 Ом м и Уч - Эмчек 40 - 90 Ом м. На участке Таш- Башат, приуроченного к средневернедевонским интрузиями, кривые кажущегося сопротивления р_к на частотах выше 100 Гц отличаются. Значения кажущегося сопротивления (р_к) в данной точке участка меняется в пределах 200- 600 Ом м, что свидетельствует о наличии неоднородностей в приповерхностной части разреза и на глубине.

Результаты мониторинга в тензочувствительной точке участка Таш -Башат с 01.06 по 06.06 2005 г показали, что для азимута MN = О^о суточные изменения р_к были менее заметные, чем для азимута МN = 90^о, но при этом наблюдается прямая корреляция графиков р_к и вертикальной приливной деформации dH.

Результаты измерений методом АМТЗ показывают, что наиболее высокие кажущиеся сопротивления отмечены на участках: Кызыл-Туу расположенного на среднечетвертичных пролювиально - аллювиальных отложениях 300 - 700 Ом м, Кызыл - Бирдик 80-300 Ом м и Уч - Эмчек 40 - 90Ом м.

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

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

Наибольшие приливные деформации характерны для вертикальной составляющей dH. На участке Таш Башат, разрез представлен гранодиоритами на глубине около 10 м и перекрывающими осадочными отложениями. По данным АМТЗ установлено, что значения кажущихся сопротивлений, на этом участке меняются от 200 до 600 Ом м. По результатам мониторинга, выполненного 04.06.2005 г. На участке Таш Башат на частоте 7 Гц, установлена прямая корреляция графиков приливной деформации dH и кажущихся сопротивлений (рис.24).

Величина вариаций составляет 2-3 % при изменениях dH до 30 см. Различный характер зависимости между графиками с_к и dH в Бишкекском регионах может быть объяснен различной степенью водонасыщения горных пород

. При мониторинге мы наблюдаем прямую корреляцию графиков с_к и dH. Отсутствие приливных вариаций с_к для однородных водонасыщенных осадочных пород иллюстрируется данными мониторинга на участке Уч Эмчек (Бишкекский регион).

На этом участке разрез сложен мощными (сотни метров) осадочными отложениями с пониженными значениями удельного сопротивления (по результатам АМТЗ до 200 Ом м).

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

Рис. 24. Сопоставление значений кажущегося сопротивления на различных частотах и вертикальной составляющей приливной деформации (dH) за период с 01.06.- 06.06.2005 г. на участке Таш-Башат: А - для азимута МN = 0⁰ и Б - для азимута МN = 90⁰ ,расчеты произведены по программе SM-27.

Кривые с_к и _к для направлений MN 0 град. и 90 град (рис. 25) близки между собою, что свидетельствует об однородности геологической среды.

Рис.25. Кривые кажущегося сопротивления и фазы импеданса двух точек на участках: а - Кызыл-Туу, б - Уч-Эмчек. Кривые ₁ и ₁ соответствуют азимуту 0, а кривые ₂ и ₂ - 90

На участке Таш Башат с высокой тензочувствительностью при изменениях dH на 30 см наблюдаются заметные вариации с_к2 (для азимута MN 0 град.) на трех частотах, в то время как на участке Уч Эмчек с низкой тензочувствительностью при изменениях dH на 35 см значения с_к2 не меняются в течение дня. Данные мониторинга иллюстрируют высокую точность измерений кажущегося сопротивления (не хуже 1 %) и возможность контроля небольших (2-3 %) изменений кажущегося сопротивления. Это позволяет сделать вывод о возможности применения метода АМТЗ, программным обеспечением, с использованием современной аппаратуры для мониторинга напряженно-деформированного состояния горных пород и использования в дальнейшем для поисков предвестников землетрясений.

Измерения естественного электромагнитного поля в аудиочастотном диапазоне позволили установить наличие корреляции между вариациями кажущегося сопротивления и приливными деформациями земной коры. Вариации с_к с амплитудой до 3-8% связаны с вертикальными приливными деформациями dH с амплитудой 15-30 см.

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

По результатам мониторинга на участке с анизотропными породами установлен различный характер корреляции с_к и dH, для случаев различной поляризации: прямая корреляция для Н-поляризации, и обратная - для Е-поляризации. Более заметные изменения с_к характерны для Е поляризованного поля.

Высокая точность измерений кажущегося сопротивления (не хуже 1%) с использованием современной аппаратуры АМТЗ позволяет регистрировать небольшие (2-3%) вариации кажущегося сопротивления. Поэтому метод АМТЗ может использоваться для мониторинга напряженно-деформированного состояния горных пород и поиска предвестников землетрясений.

На Иссык-Кульском полигоне полевые работы методом АМТ были проведены на четырех точках в пределах северной части - пункты Григорьевка и Семеновка, и одном в пределах юго-западной части Иссык-Кульского полигона - пункт Оттук. Положение точек зондирования было выбрано с учетом имеющейся геологической информации и помеховой обстановки. Участок Орнёк расположен в зоне развития осадочных отложений, участки Григорьевка и Семеновка приурочены к области развития интрузивных пород кислого состава палеозойского возраста. Участок Оттук расположен на мощной толще осадочных отложений кайнозойского возраста. Участки Орнёк, Григорьевка и Семёновка расположены в зоне влияния новейших краевых разломов - Культорского, Талдысуйского и Южно-Аксуйского. Участок расположен на значительном удалении от известных разломов.

На каждой точке осуществлялась регистрация сигналов горизонтальных электрических и магнитных антенн в диапазоне частот от 5 до 400 Гц. Обработка данных осуществлялась с использованием программ SM-27 и SM+. При обработке осуществлялось вращение тензора импеданса. Значения кажущихся сопротивлений составляют десятки - сотни Ом м. По расчетам суммарная мощность рыхлых проводящих осадочных отложений составляет около 1600 м. Особенности кривых зондирования для точки на участке Григорьевка указывают на то, что данная точка, расположенная в Григорьевском ущелье, пригодна для проведения электромагнитного мониторинга для поисков предвестников землетрясений. Уровень кажущихся сопротивлений составляет сотни - тысячи Ом•м. Наблюдается значительный разбег кривых кажущегося сопротивления и фазы импеданса для разных поляризаций, что свидетельствует о высокой степени неоднородности среды. Точка на участке Семёновка расположена в Семеновском ущелье, в 4 км восточнее от участка Григорьевка. Данная точка расположена в зоне распространения осадочных пород кайнозоя. Это проявляется в снижении уровня кажущегося сопротивления по сравнению с участком Григорьевка до десятков-сотен Ом•м.

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

Рис. 26. Кривые кажущегося сопротивления и фазы импеданса для точек Иссык-Кульского полигона.



Рис. 27. Геоэлектрический разрез, полученный по результатам 2D инверсии данных АМТ зондирования для точек Оттук-Орнёк-Григорьевка (расчеты проведены по программе SM-27).

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

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

На участке Джарды-Суу измерения выполнены на двух точках, для них характерны пониженные значения кажущегося сопротивления для обоих значений (сотни Ом•м). Пункт расположен в зоне развития осадочных отложений. Глубина до фундамента по оценочным расчетам составляет около 1,5 км. Для п. 1, расположенной на 4,5 км южнее п. 2, характерны повышенные значения удельного сопротивления (первые тысячи Ом•м) (рис.28).

Следует отметить, что данная точка расположена в зоне влияния Шамси-Тюндюкского разлома. В целом, данная точка перспективна для проведения мониторинга электромагнитных предвестников при поисках землетрясений. Пункт 1 участка Белогорка расположена в зоне влияния Шамси-Тюндюкского разлома. Для неё характерны повышенные значения удельного сопротивления (тысячи-десятки тысяч Ом•м) и высокая степень неоднородности среды. Данная точка также перспективна для проведения электромагнитного мониторинга.

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

Рис.28. Кривые кажущегося сопротивления и фазы импеданса для точек Бишкекского полигона.

В ходе проведения экспериментальных исследований на Бишкекском и Иссык-Кульском полигонах, основываясь на полученных данных и их анализе, были выбраны тензочувствительные пункты для проведения мониторинга электромагнитных предвестников землетрясений. Для Иссык-Кульского полигона наиболее перспективным является участок Григорьевка, для Бишкекского полигона - все изученные участки. Проведение АМТ-зондирований для оценки степени неоднородности точек наблюдений и мониторинга для оценки возможностей регистрации приливных вариаций кажущегося сопротивления позволили выбрать тензочувствительные пункты для проведения мониторинга сейсмической активности в пределах Бишкекского прогностического полигона. В результате проведения работ впервые на Бишкекском полигоне были зарегистрированы приливные вариации кажущегося сопротивления.

При изучении аномалий ЭМИ производится вычисление автоспектров и динамических спектров

Анализ полученных материалов показывает, что в диапазоне частот 200-300 Гц наблюдались локальные изменения уровня поля, которые потенциально могут быть связаны с излучением локального источника вблизи точки наблюдения (зоны разлома и др.). В диапазоне частот 8-40 Гц хорошо проявлены сигналы Шумановских резонансов естественного электромагнитного поля в волноводе земля - ионосфера. Динамические спектры являются более информативными по сравнению с автоспектрами. Они дают возможность не только анализировать спектральный состав зарегистрированного за определенный интервал времени сигнала, но и получить представление о появлении или исчезновении кратковременных сигналов в пределах этого интервала. Анализ особенностей этих кратковременных сигналов может быть наиболее информативным при выделении аномалий ЭМИ перед землетрясениями

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

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

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

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

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

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

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

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

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

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

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

2.Выделены региональные области вариаций магнитного поля интенсивностью от 1,0 до 10 нТл, имеющие значительные размеры, которые вероятнее всего, обусловлены региональными процессами, протекающими под земной корой или с течением электрических зарядов, образующихся в результате вязкого течения подкорового вещества. Выявлены более мелкие по размерам локальные изменения магнитного поля интенсивностью от 3 до 17 нТл. Такие изменения магнитного поля, возможно, обусловлены деформационными процессами в земной коре, сопровождающими подготовку и проявление сейсмических событий.

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

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

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

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

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

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

9.Разработаны электромагнитные методы на основе моделирования процессов подготовки землетрясений, опробованы различные подходы и намечены пути совершенствования методики изучения аномалий ЭМИ. Впервые установлены закономерности поведения электромагнитного поля в различных частотных интервалах сверхнизкочастотного диапазона, определены частотные интервалы.

10.Впервые на территории Бишкекского и Иссык-Кульского геодинамического полигона создана периодически действующая экспериментальная сеть АМТ-мониторинга сейсмической активности, с программно-управляемым функционированием с последующим использованием для постоянно действующего мониторинга. На Бишкекском полигоне были зарегистрированы приливные вариации кажущегося сопротивления и практически показаны возможности мониторинга изменений напряженно-деформированного состояния геологической среды.

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

ОСНОВНОЕ СОДЕРЖАНИЕ ДИССЕРТАЦИИ ОТРАЖЕНО В ПУБЛИКАЦИЯХ

В монографии

1. Бакиров, К.Б. Геомагнитные явления и сейсмотектонические процессы Кыргызского Тянь-Шаня [Текст]: монография /Бакиров К.Б./ ? Бишкек: Илим, 2005. - С.144

Статьи в научных журналах

2. Бакиров, К.Б. Результаты наземных геомагнитных исследований. Геологическое строение и сейсмичность Тотогульского и Фрунзенского полигонов [Текст]/Бакиров К.Б., Шакиров Э.Ш //. Матер. междунар. аэроксм. эксперимент. Тянь-Шань-Интеркосмос-88. - Фрунзе: Илим, 1988 - С. 46-47.

3. Бакиров, К.Б. Результаты исследований по региональной сети магнитных пунктов Северной Киргизии [Текст]: тез. докл. /Бакиров К.Б. // Мат-лы конф. молод. учен. - Фрунзе: Илим,I989. - С. 17-18.

4. Бакиров, К.Б. Магнитометрические исследования вблизи Токтогульского водохранилища [Текст]/Бакиров К.Б., Шакиров Э.Ш.// Сейсмичность Тянь-Шаня. - Фрунзе: Илим 1989. - С. 247-264.

5. Бакиров, К.Б. Возбужденная сейсмичность чаши Токтогульского водохранилища в связи с изменением объёма воды [Текст]/Шакиров Э.Ш., Бакиров К.Б., Кожоева Г.М.// Сейсмичность Тянь-Шаня. - Фрунзе: Илим, 1989. - С. 238-246.

6. Бакиров, К.Б. Магнитометрические исследования в сейсмогенных зонах Фрунзенского прогностического полигона [Текст]: тез. докл. /Бакиров К.Б.//Материалы Всесоюзной конф. молод. учен. - Фрунзе: Илим, 1990. - С. 18-I9.

7. Бакиров, К.Б. Глубинное строение Таласо-Ферганского разлома в районе Токтогульского водохранилища по геомагнитным данным [Текст]: тез. докл. /Бакиров К.Б., Шакиров Э.Ш. // Всесоюзной школы-семинар 5-14 сентября Геолого-геофизические исследования в сейсмоопасных зонах СССР. - Фрунзе: Илим, 1989. - С. 82-83.

8. Бакиров, К.Б. Магнитометрические исследования на территории Киргизии [Текст] / Шакиров Э.Ш., Бакиров К.Б., Погребной В.Н., Ким Л.Е., Янкевич Т.Н. // Всесоюзный съезд по геомагнетизму. - Москва-Суздаль, 1991.

9. Бакиров, К.Б. Изучение локальных геомагнитных вариаций в сейсмоактивных районах [Текст]/Бакиров К.Б. // Геодинамика и металлогения полезных ископаемых. /Материалы междунар. науч.-технич. конф. ? Бишкек, 1999. ? С. 21 ? 25

10. Бакиров, К.Б. Условия образования отрицательных магнитных полей, локальный разогрев и связь с эпицентрами сильных землетрясений [Текст]/Бакиров К.Б., Шакиров А.Э. // Экологический вестник Кыргызстана. - Бишкек, 2000. ? № 4. ? С. 13-14.

11. Бакиров, К.Б. Результаты исследований наклономерно-деформометрических процессов в приделах Бишкекского прогностического полигона. Экология и безопасность жизнедеятельности в горных условиях. [Текст]/Бакиров К.Б.//Материалы науч.-практ. конф. ? Бишкек: Илим, 2000. ? С. 32 ? 36.

12. Бакиров, К.Б. Исследования геомагнитных и гидрогеохимических полей в целях прогноза сильных землетрясений [Текст] /Бакиров К.Б., В.В. Гребенникова, А.А. Усенова. // Наука и новые технологии. - № 3. ? Бишкек, 2000. ? С. 22 - 25.

13. Бакиров, К.Б. Отрицательные магнитные аномалии. [Текст]/Бакиров К.Б. //Недра гор Кыргызстана - народу. Материалы науч.-тех. конф. .? Бишкек, 2003. ? С. 66 ? 69.

14. Бакиров, К.Б. Исследование локальных изменений геомагнитного поля /Бакиров К.Б., Берикова Г.К. // Недра гор Кыргызстана - народу. Материалы науч.-тех. конф. .? Бишкек, 2003. ? С. 72-76.