MISCELLANEOUS / Geophysics / Geophysics Methods Костицын В. И
..pdfВведение
Таким образом, разведочная (прикладная) геофизика и скважинная (промысловая) геофизика – это многогранные на-
учно-прикладные дисциплины со сложной структурой. Существуют разные подходы к их классификации:
1)по геофизическим полям (гравиразведка, магниторазведка, электроразведка, сейсморазведка, терморазведка и ядерная геофизика);
2)по технологиям и месту проведения работ (аэрокосмические, полевые, акваториальные, подземные методы и геофизические исследования скважин);
3)по прикладным направлениям и решаемым задачам (глубинная, региональная, разведочная, инженерная и экологическая геофизика);
4)по видам деятельности (теоретическая, инструментальная, экспериментальная, вычислительная и интерпретационная геофизика).
История развития геофизики. В период с середины
XIX в. до 1920-х годов возникла общая геофизика. Отдельные ее отрасли зародились и развивались в пределах других дисциплин: физики, механики, астрономии, геодезии. Возникли гравиметрия, учение о земном магнетизме, сейсмология, электрометрия. Активно геофизические методы исследования недр начали развиваться с 1920-х гг. Однако ее физико-математические основы заложены значительно раньше. Так же давно началось использование физических полей Земли для практических целей.
Ранее других методов возникла магниторазведка. Первые сведения о применении компаса для разведки магнитных руд в Швеции относятся к 1640 г. Теория гравиметрии берет свое начало с 1687 г, когда Ньютон сформулировал закон всемирного тяготения. В 1753 г. М.В.Ломоносов высказал мысль о связи значений силы тяжести на земной поверхности с внутренним строением Земли и разработал идею газового гравиметра. Его же работы в области сейсмологии, атмосферного электричества можно считать первыми, относящимися к геофизическим исследованиям Земли. Первыми работами по электроразведке можно считать наблюдения Р.Фокса (Англия) в 1830 г. естественной поляризации сульфидных залежей и Е.И.Рогозина (Россия), который в 1903 г. дал первое изложение основ этого метода. В
19
В.К. Хмелевской, В.И. Костицын
1913 г. К. Шлюмберже (Франция) разработал метод электроразведки постоянным током, а в 1918 г. К. Зунберг и Н. Лунберг (Швеция) предложили электроразведку переменным током. Со времени установления Кулоном закона взаимодействия магнитных масс (1785 г.) начинается развиваться теория земного маг-
нетизма.
Первыми систематическими разведочными работами в России и в мире были съемки Курской магнитной аномалии (КМА), начатые профессором МГУ Э.Е. Лейстом в 1894 г., а также проведенные магнитные съемки на Урале Д.И. Менделеевым и в районе Кривого Рога И.Т. Пассальским в конце XIX в. В 1919 г. были начаты магнитные съемки на КМА. Именно эти работы можно считать началом развития отечественной разведочной геофизики. Теоретические работы Э. Вихерта (Германия) и Б.Б. Голицына (Россия) в начале XX в. в области сейсмологии имели самое непосредственное отношение к созданию сейсморазведки.
За период развития советской и российской геофизики
внесли большой вклад в области гравиразведки: Б.А. Андреев, А.Д. Архангельский, Б.К. Балавадзе, В.М. Березкин, Ю.Д. Буланже, Е.Г. Булах, К.Е. Веселов, Г.М. Воскобойников,
Г.А. Гамбурцев, |
К.В. Гладкий, Ю.Н. Годин, |
Г.Я. Голиздра, |
|
В.И. Гольдшмидт, |
Н.П. Грушинский, |
|
М.С. Жданов, |
А.А. Заморев, В.К. Иванов, Н.И. Идельсон, |
И.Г. Клушин, |
||
С.С. Красовский, |
З.А. Крутиховская, |
М.М. Лаврентьев, |
|
О.К. Литвиненко, |
П.И. Лукавченко, |
|
Е.Н. Люстих, |
В.А. Магницкий, |
А.К. Маловичко, |
|
А.А. Михайлов, |
Е.А. Мудрецова, |
Л.Д. Немцов, |
А.А. Непомнящих, |
|
П.М. Никифоров, |
П.С. Новиков, |
В.М. Новоселицкий, |
|
Б.В. Нумеров, |
В.Л. Пантелеев, |
С.А. Поддубный, |
М.А. Садовский, М.У. Сагитов, Н.Б. Сажина, П.К. Соболевский, Л.В. Сорокин, С.И. Субботин, К.Ф. Тяпкин, Д.Г. Успенский, В.В. Федынский, Э.Э. Фотиади, А.В. Цирульский, О.А. Шванк, П.Ф. Шокин, В.И. Шрайбман, А.А. Юньков, М.И. Юркина, А.Я. Ярош; в области магниторазведки: А.Д. Архангельский, И.П. Бахурин, В.П. Боронин, Б.К. Вейнберг, Ф.М. Гольцман, Г.И. Гринкевич, А.П. Казанский, В.В. Колюбакин, А.А. Логачев, Д.С. Миков, Н.В. Розе, В.А. Сельский, Т.Н. Симоненко,
20
Введение
Б.М. Яновский; в области электроразведки: Л.М. Альпин, И.М. Блох, В.Р. Бурсиан, Л.Л. Ваньян, А.В. Вешев, В.Н. Дахнов,
А.И. Заборовский, |
Е.Н. Каленов, А.П. Краев, |
А.А. Огильви, |
|
И.К. Овчинников, |
А.А. Петровский, |
А.М. Пылаев, |
|
А.С. Семенов, |
А.Г. Тархов, |
А.Н. Тихонов, |
В.А. Фок, |
Ю.В. Якубовский; |
в области |
сейсморазведки: |
И.С. Берзон, |
Г.Н. Боганик, В.С. Воюцкий, Е.И. Гальперин, Г.А. Гамбурцев, Б.Б. Голицын, С.В. Гольдин, И.И. Гурвич, Е.В. Карус, О.К. Кондратьев, П.П. Лазарев, В.И. Мешбей, В.П. Номоконов,
Г.И. Петрашень, |
М.К. Полшков, |
Н.Н. Пузырев, |
Ю.В. Ризниченко, |
Л.А. Рябинкин, |
Е.Ф. Саваренский, |
А.К. Урупов, О.Ю. Шмидт; в области |
ядерной геофизики: |
В.И. Баранов, Ю.П. Булашевич, В.И. Вернадский, В.В.Ларионов,
Е.М.Филиппов; в области геофизических исследований скважин (промысловой геофизики): Л.М. Альпин, М.И. Бальзамов,
О.К. Владимиров, В.Н. Дахнов, Г.В. Горшков, А.Н. Граммаков, Д.В. Жабрев, С.С. Итенберг, С.Г. Комаров, А.М. Нечай, Н.А. Савостьянов, В.А. Сельский, А.С. Семенов, М.М. Соколов, В.А. Фок, В.В. Шаскольский, В.А. Шпак.
В настоящее время по уровню развития теории обработки и интерпретации данных геофизических методов и практическому их использованию отечественная геофизика занимает передовые позиции в мире. Ведущими учеными в этих направлениях являются по гравиразведке: Ю.Н. Авсюк, Ю.В. Антонов,
В.И. Аронов, |
П.И. Балк, А.С. Варламов, |
А.Г. Гайнанов, |
В.А. Гершанок, |
В.М. Гордин, А.С. Долгаль, |
Л.Я. Ерофеев, |
Г.И. Каратаев, А.И. Кобрунов, Ю.П. Конценебин, В.Н. Конешов, В.И. Костицын, М.И. Лапина, В.В. Ломтадзе, П.С. Мартышко, В.С. Миронов, В.Р. Мелихов, В.О. Михайлов, И.Н. Михайлов, В.А. Романюк, С.А. Серкеров, З.М. Слепак, В.И. Старостенко, В.Н. Страхов, В.Г. Филатов, М.С. Чадаев; по магниторазведке: Ю.И. Блох, Л.А. Гершанок, Ю.С. Глебовский, Ф.М. Гольцман, В.М. Гордин, А.М. Городницкий, Л.А. Золотая, Т.Б. Калинина, В.Е. Никитский, Д.К. Нургалиев, М.Г. Попов, Г.А. Трошков, В.В. Филатов; по электроразведке: М.Н. Бердичевский, В.А. Богословский, Д.С. Даев, В.П. Колесников, Б.К. Матвеев,
И.Н. Модин, |
П.Ю. Пушкарев, |
А.А. Редозубов, |
В.М. Сапожников, |
Б.С. Светов, В.В. Спичак, В.К. Хмелевской, |
|
|
|
21 |
В.К. Хмелевской, В.И. Костицын
В.А. Шевнин, М.Н. Юдин, А.Г. Яковлев; по сейсморазведке:
Ю.П. Ампилов, |
Л.Д. Бовт, |
В.И. Бондарев, |
|
М.Л. Владов, |
Г.Н. Гогоненков, |
Е.А. Козлов, О.Л. Кузнецов, |
Т.И. Облогина, |
||
В.Б. Писецкий, |
О.А. Потапов, |
М.Б. Рапопорт, |
В.И. Рыжков, |
И.А. Санфиров, Б.А. Спасский, Ю.В. Тимошин; по ядерной
геофизике: |
В.А. Мейер, |
Р.А. Резванов, |
А.Г. Талалай, |
||
В.И. Уткин; по |
геофизическим исследованиям скважин: |
||||
Н.С. Березовский, |
А.М. Блюменцев, |
|
Р.А. Валиуллин, |
||
Б.Ю. Вендельштейн, Г.С. Возжеников, |
И.Н. Гайворонский, |
||||
Ю.И. Горбачев, |
|
С.И. Дембицкий, |
|
В.М. Добрынин, |
Ф.Х. Еникеева, В.Ю. Зайченко, Г.М. Золоева, В.И. Иванников,
Л.Е. Кнеллер, В.Ф. Козяр, Ю.И. Кузнецов, |
В.В. Лаптев, |
|
Э.Е. Лукьянов, |
А.А. Молчанов, |
А.С. Некрасов, |
В.М. Сапожников, |
И.Г. Сковородников, |
В.Г. Фоменко, |
М.И.Эпов, Г.Г. Яценко; по комплексированию геофизических методов: В.В. Антонов, С.В. Аплонов, В.М. Бондаренко, В.В. Бродовой, А.А. Булычев, Г.С. Вахромеев, В.Н. Глазнев, В.В. Калинин, А.А. Кауфман, Н.Я. Кунин, Ф.М. Ляховицкий, А.А. Никитин, В.Е. Никитский, Г.П. Новицкий, В.В. Попов, Н.И. Селиверстов, В.С. Сурков, А.Г. Тархов, С.В. Шалаев, Н.В. Шаров, С.А. Шихов, В.К. Хмелевской.
В основу данного издания положены учебники и учебные пособия В.К. Хмелевского, изданные им в разные годы в Московском государственном университете (МГУ), Международном университете природы, общества и человека «Дубна» и лекции В.И. Костицына, читаемые студентам Пермского государственного университета (ПГУ) по специальности «Геофизика».
Авторы благодарят Наталью Марковну Ржевитину за долгий, кропотливый и нелегкий труд по подготовке компьютерного макета учебника.
Предложения и замечания можно направлять по адресу: 614990, г. Пермь, ул. Букирева, 15. Пермский государственный университет, кафедра геофизики. Костицыну Владимиру Ильичу.
Факс: (342) 237-16-11, (342) 239-68-32. E-mail: geophysic@psu.ru; kostitsyn@psu.ru.
22
Гравиразведка
ГЛАВА 1. ГРАВИРАЗВЕДКА
Гравиметрическая, или гравитационная разведка (со-
кращенно гравиразведка) – это геофизический метод исследования земной коры и разведки полезных ископаемых, основанный на изучении распределения аномалий силы тяжести на земной поверхности, акваториях, в воздухе. Существование поля силы тяжести обусловлено в основном ньютоновским притяжением Землей всех тел, обладающих массой. Так как Земля сферически неоднородна, да еще вращается, то поле силы тяжести на земной поверхности непостоянно. Изменения эти малы и требуют высокочувствительных приборов для их изучения. Основными измеряемыми параметрами гравитационного поля являются ускорение силы тяжести и градиенты (изменения ускорения по разным направлениям). Величины параметров поля силы тяжести зависят, с одной стороны, от причин, обусловленных притяжением и вращением Земли (нормальное поле), а с другой стороны – от неравномерности изменения плотности пород, слагающих земную кору (аномальное поле). Эти причины изменения силы тяжести на Земле послужили основой двух на-
правлений гравиметрии: геодезической гравиметрии и грави-
тационной разведки.
От других методов разведочной геофизики гравиразведка отличается сравнительно большой производительностью полевых наблюдений и возможностью изучать горизонтальную (латеральную) неоднородность Земли. Гравиразведка применяется для решения самых различных геологических задач с глубинностью исследований от нескольких метров (например, при разведке окрестностей горных выработок) до 200 км (например при изучении мантии).
1.1. Сила тяжести
Силой тяжести (F) называют равнодействующую двух сил – силы ньютоновского притяжения всей массой Земли (Fн) и центробежной силы, возникающей вследствие суточного вращения Земли (P). Отнесенные к единице массы, эти силы харак-
23
В.К. Хмелевской, В.И. Костицын
теризуются ускорениями силы тяжести g = F/m, ньютоновского притяжения f = Fн/m и центробежным p = P/m. Ускорение силы тяжести равно геометрической сумме ускорения притяжения и центробежного ускорения (рис. 1.1). Обычно в гравиметрии, когда говорят «сила тяжести», подразумевают именно ускорение силы тяжести.
Рис. 1.1 Ускорение силы тяжести и его составляющие
Единицей ускорения в системе СИ является м/с2. В гра-
виметрии традиционно используют единицу – Гал, равную 1 см/с2. В среднем на Земле значение силы тяжести составляет 981 Гал. В практике гравиразведки применяется величина в 1000 раз меньшая, получившая название миллигал (мГал), и в 1 000 000 раз меньшая – микрогал (мкГал).
Сила притяжения какой-либо массы (m) всей массой Земли (M) определяется законом всемирного тяготения Ньютона:
F G |
mM |
, |
(1.1) |
н |
r 2 |
|
где r – расстояние между центрами масс m и M, т.е. радиус Земли; G – гравитационная постоянная, равная 6,67·10-11 м3/кг·с2. Сила притяжения единичной массы (m = 1) равна f ≈ GM/r2 и направлена к центру Земли.
24
Гравиразведка
Центробежная сила (P) направлена по радиусу, перпендикулярному оси вращения (R), и определяется формулой
P mR 2 , |
(1.2) |
где ω – угловая скорость вращения Земли.
Величина P изменяется от нуля на полюсе (R = 0) до максимума на экваторе. Отношение P/F ≤ 1/288, поэтому сила тяжести почти целиком определяется силой притяжения, а ускорение
силы тяжести практически равно ускорению притяжения g ≈ f ≈ GM/r2.
Земля в первом приближении является эллипсоидом вращения, причем экваториальный радиус = 6378 км, а полярный c = 6357 км, - c = 21км. Разная величина радиуса Земли на полюсе и экваторе наряду с изменением центробежной силы приводит к увеличению g на полюсе (gп = 983 Гал) по сравнению с g на экваторе (gэ’= 978 Гал). По известным g и r были определены масса Земли М = 5,98 1024 кг и ее средняя плотность σЗ = 5,51·103
кг/м3 (5,51 г/см3).
1.2. Потенциал силы тяжести
Потенциал силы тяжести (W) был введен в теорию гра-
виметрии для облегчения решения теоретических задач. В точке А, расположенной на расстоянии rA от центра Земли, выражение для потенциала принимается равным WA = GM/rA, а в любой
точке |
B, |
расположенной |
на |
продолжении радиуса r, |
||
WB |
GM / rA |
r . Поэтому разность потенциалов будет рав- |
||||
на |
|
|
|
|
|
|
|
W |
WB |
WA GM / |
|
r |
. |
|
|
|
||||
|
rA rA |
rB |
||||
|
|
|
|
|
||
|
В пределе при малом |
r имеем |
|
|||
|
|
|
|
|
|
25 |
В.К. Хмелевской, В.И. Костицын
W GM r / r 2 g r ,
отсюда g = – dW/dr, т.е. сила тяжести есть производная потенциала силы тяжести по направлению к центру Земли.
С другой стороны, работа, которая может быть произведена при движении притягиваемой точки по отрезку r , равна A = g r . Поэтому W = – A, или работа силы тяжести по перемещению единичной массы на отрезке r равна разности зна-
чений потенциала на концах этого отрезка.
При перемещении точки в направлении, перпендикулярном силе тяжести, dW = 0. Это означает, что W = const. Поэтому гравитационное поле можно представить в виде набора бесконечного числа поверхностей, на которых потенциал остается постоянным, а сила тяжести направлена перпендикулярно этой поверхности. Такие поверхности называют эквипотенциальными, или уровенными. В частности, поверхность жидкости на Земле, например моря, совпадает с уровенной поверхностью. У Земли есть одна уникальная уровенная поверхность, которая совпадает с невозмущенной волнениями поверхностью океанов. Она называется геоидом.
Таким образом, геоид – это условная уровенная поверхность, совпадающая со средним уровнем океанов, открытых морей и продолженная под континентами так, чтобы она в каждой точке была горизонтальна, т.е. перпендикулярна линиям силы тяжести.
1.3. Производные потенциала силы тяжести
Производные потенциала силы тяжести по трем координатным осям gx = W/ x, gy = W/ y, gz = W/ z однозначно определяют его полный вектор.
Если ось z направить к центру Земли, то |
W |
|
W |
0 |
, а |
|
|
|
|||
x |
|
у |
|||
|
|
|
|
g Wz .
26
Гравиразведка
В гравиметрии кроме первых производных изучаются
вторые производные потенциала или их разности:
|
2W |
2W |
2W |
|
2W |
2W |
2W |
2W |
|
2W |
|
|
|
|||||||
|
|
, |
|
, |
|
, |
|
, |
|
, |
|
, |
|
|
|
. |
|
(1.3) |
||
|
x y |
x z |
y z |
x2 |
y 2 |
z 2 |
x2 |
y 2 |
||||||||||||
Физический смысл этих выражений легко получить, если |
||||||||||||||||||||
иметь в виду, что g = |
W/ z. Так, вторая производная |
2W |
|
g |
|
|||||||||||||||
x z |
|
x |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
указывает на скорость изменения силы тяжести по оси х, т.е.
является горизонтальным градиентом силы тяжести.
Аналогичный смысл имеют вторые производные
2W/ x z и 2W/ z2.
Вторые производные 2W/ x y, 2W/ x2- 2W/ y2 харак-
теризуют форму уровенной поверхности (геоида), изучаемую в геодезической гравиметрии. Практической единицей измерения градиента силы тяжести принимается 1 этвеш (Е) = 10-9/c2, что соответствует изменению силы тяжести в 0,1 мГал на 1 км.
1.4. Нормальное значение силы тяжести
Нормальным значением силы тяжести ( 0) называется сила тяжести, обусловленная суточным вращением и притяжением Земли, в предположении, что она состоит из однородных по плотности концентрических слоев.
Принимая Землю за сфероид, Клеро получил следующую приближенную формулу:
0 |
g |
э |
1 |
sin 2 , |
|
|
|
где gэ – сила тяжести на экваторе; – географическая широта пункта наблюдения; – коэффициент, зависящий от угловой скорости вращения и сжатия сфероида.
Однако Земля соответствует по форме в большей степени геоиду, тогда нормальные значения силы тяжести для его поверхности рассчитываются по формуле
27
В.К. Хмелевской, В.И. Костицын
0 gЗ 1 |
sin 2 |
1 sin 2 2 |
2 cos2 cos2 , |
(1.4) |
где – географическая долгота пункта наблюдения. Коэффициенты , 1 и 2 зависят от формы Земли, ее уг-
ловой скорости вращения, распределения масс. Наиболее широ-
ко используется нормальная формула Гельмерта
0 |
978030(1 0,005302sin 2 |
0,000007sin 2 2 ), |
(1.5) |
|
|
|
где – широта гравиметрического пункта. Имеются и другие формулы, полученные исследователями в разных странах и в разные годы.
При наличии координат пунктов нормальные значения могут быть найдены по таблицам Бурдюкова. В условиях, когда неизвестны точные географические координаты, а известна лишь некоторая средняя широта участка работ, нормальное значение для I -й точки может быть определено по формуле
i |
оп |
|
xi , |
|
|
|
|
|
|
||||
|
|
l |
|
|
||
где оп – нормальное значение на опорном пункте, |
|
|
– гори- |
|||
|
l |
|||||
|
|
|
|
|
|
|
зонтальный градиент изменения нормального поля, |
|
хi – рас- |
стояние по широте от опорной точки до текущей в метрах. Таким образом, геоид является поверхностью относимо-
сти, по отношению к которой рассчитываются аномалии.
1.5. Редукции силы тяжести
В наблюденные значения силы тяжести вводятся поправки (редукции). Введение поправок необходимо потому, что нормальные значения относятся к поверхности геоида, которая совпадает с уровнем океана, а измеренные значения относятся к действительной (реальной) земной поверхности. Для того чтобы все наблюдения силы тяжести были сопоставимы, их приводят к одной поверхности – уровню геоида, т.е. как бы опускают точку
28