shimbirev_b_p_teoriya_figury_zemli
.pdfВыведем формулу для аномалии силы тяжести для случая, когда прямая имеет бесконечное простирание. Представим (Х1У.26) в виде
АСо
Т]2
Положив т]2 = |
г]х = —со, найдем |
Шо
|
> |
|
*кт(0,0, с; |
|
г |
Рис. 81 |
Рис. 82 |
Для вычисления производной Тхг в формуле (XIV. 14) положим
н = 0, у = 2 = 0, г = |
+ |
+ |
= |
Тогда
'X2
Т
(XIV. 27.
Р(х,0,0)
Подынтегральное выражение можно представить как
йц • ^ |
(1 —/?'2) ан' |
|
||
н*> |
|
а |
|
|
где по-прежнему В = ]А]2 + с2 |
, с2 |
= |
ж2 + Си- |
|
В таком случае |
|
|
|
|
г |
|
I |
(ЗД'-Д'3). |
|
КЬ |
' Зс4 V |
|
||
Следовательно,, |
|
|
|
|
1)2 |
|
|
3%, |
|
Г |
|
|
|
|
|
|
|
Уг)2+-2+$02 |
+ |
ЗТ)1 |
|
|
Т)3 |
|
.400
преобразуя полученное выражение, получим |
|
ГЬ |
1 |
<1ц |
|
|
|
|
|
1 + |
ч1 |
|
|
|
|
|
|
|
|
|
|
|
|
V 1- |
Ч |
|
ц2 |
|
|
Полагая теперь Т]2 = |
|
-ОО |
получим |
формулу для |
прямой |
||
бесконечного |
простирания |
гр |
_ 4/а^рХ |
|
|
|
|
|
|
|
|
|
(Х1У.28) |
||
|
|
|
|
(**-К.а)2" |
|
||
|
|
|
|
|
|
||
Формулы |
(XIV. 27) и (XIV. 28) могут быть использованы и для |
однород- |
|||||
ного горизонтального |
кругового |
цилиндра |
бесконечного простирания, если |
||||
линейную плотность |
К положить |
равной пЯ2 Дб, где |
Я — радиус |
сечения |
|||
цилиндра плоскостью, перпендикулярной к его оси; Дб — аномальная плотность, длина цилиндра принята равной единице.
На рис. 81 представлены кривые Д^ и Тхг для случая горизонтального
бесконечного цилиндра кругового сечения, построенные на основании формул
(XIV.27) и |
(Х1У.28). |
|
|
|
|
|
2. В е р т и к а л ь н а я в е щ е с т в е н н а я |
п р я м а я . |
В этом слу- |
||||
чае (рис. 82) следует |
положить |
|
|
|
|
|
|
~ Л |
= 0, у = г = |
0, г = У®я-К3. |
<1т = Х<11,. |
|
|
В соответствии с |
(Х1У.11) |
|
|
|
|
|
|
Ь2 |
|
|
|
|
|
|
|
с я |
Ф |
|
|
|
|
|
|
|
|
|
|
Если |
= оо, то |
|
А |
|
|
|
|
|
дёГ: |
|
|
(Х1У.29) |
|
|
|
|
|
|
||
3. П р я м о у г о л ь н ы й |
п а р а л л е л е п и п е д |
и |
п р я м о - |
|||
у г о л ь н а я п р и з м а б е с к о н е ч н о г о |
п р о с т и р а н и я . На |
|||||
рис. 83 изображен прямоугольный параллелепипед, грани которого параллельны координатным плоскостям. Аномальную плотность параллелепипеда
примем равной Дб = б2 — бх, где |
б2 |
|
— плотность самого параллелепипеда, |
|
а бх — плотность окружающих горных |
|
пород. Вычислим аномалию силы тя- |
||
жести в точке О — начале координат. Полагая в (XIV. 11) |
||||
найдем |
|
Ы |
||
$2 >12 |
||||
Д, = /Дб | | |
| |
С Я лц аъ |
||
(Б*-И2-Ка) '/2 |
||||
|< 111 |
|
|
||
|
|
|
||
26 Заказ 1379 |
4 0 1 |
После выполнения интегрирования по ? и т)
?2ТЪ Ъ |
|
|
|
|
^ - - / Д б | | |
+ |
+ |
+ |
(XI V > I |
Выполнив интегрирования по |
получим |
|
|
|
Д^ = - / Д б | | I [Б1П(Л + |
Г) + Л |
|
+ |
(XIV : I |
с.тнб. |
|
|
|
|
В случае бесконечной призмы с прямоугольным сечением, грани коте параллельны плоскостям ху ж уг (рис. 84), в формуле (XIV.30) пределы к »
Рис. 83
следует считать от — д о Обозначив через I достаточно большое чиедч.- будем иметь
|1п(т1+/т1® + 1а + С2 )Й^ = 1п |
У |
|
|
- 1 + У |
"Г" п |
|
ИВЫкли, ' |
Считая дробь (§2 + малой величиной, разложим числитель и знамгнатель полученного выражения в ряд по степеням этой дроби, ограничивая: л малыми первого порядка
. |
У |
|
|
^ |
^ |
т |
. |
и |
П |
, .ч |
1П |
. - . Л |
. |
— |
|
1~~|ЧЧ2 |
= |
( |
й^+ё2" |
/ ' |
|
|
+ ] / н |
|
|
2 |
1* |
|
|
|
|
|
или, пренебрегая единицей^ |
поскольку |
первый |
член |
в скобке |
очень велкг_ |
|||||
|
|
|
1 + |
/ 1 - ' |
|
^ |
1п • |
|
|
|
|
|
|
1п |
— |
|
— ^ |
|
|
||
402.
Таким образом, формула (Х1У.30) для случая прямоугольной призмы бесконечного простирания будет
Д^ = —/Аб |
1п- |
422 |
= —/Дб | $ 1п 4/2 |
•1п |
4/2 |
Е, 1=со |
|
ё2 -к2 |
г=оо ; |
|
+ ^ |
или
А^ = /Аб | 1п
|2 + С!2
Применяя формулу интегрирования по частям, найдем
11п |
+ |
^ = 11п (52 + |
С2) - 2 ^ |
= § 1п (52 + О |
- 2| + |
агс1В |
|
Следовательно, |
|
|
|
|
|
||
|
|
|
-11 |
|
6. |
|
|
= |
/б [121п (I2 + С22) + |
2С2 агс*ё |
- & 1п (Й + ?22) - |
2^2 агс1§ |
- |
||
- |
ё21п (51 + |
II) ~ Чг агс1? |
+ & 1п {Ц + Й) + |
агс1? |
(Х1У.32) |
||
или в сокращенной |
записи |
|
|
|
|
||
|
|
д^ = /Лб |
|
|
|
(Х1У.ЗЗ) |
|
Положим, что сечение плоскостью хг рассматриваемой призмы представляет прямоугольник МNР^ (см. рис. 84). Аномалия силы тяжести Ад, вызванная этой призмой, может быть представлена как разность двух сумм:
Ад± + Ад2 и А§3 + |
|
где Д^, |
А§2, |
Ад3, |
|
— аномалии силы тяжести |
|||||
призм, |
сечения |
которых |
соответственно |
будут: 0%,2Р1,2, О ^ М ^ , О |
и |
||||||
О ^N11, 1- Ддя вычисления величин |
|
Д^2, Д^3 и |
положим в формуле |
||||||||
(XIV.32) пределы |
и |
равными нулю, после чего получим |
|
|
|||||||
I» |
|
|
I» |
|
|
|
|
|
|
|
|
| 1п |
+ 0 |
<11- |
| 1п (52 + Й) |
= & 1п |
+ |
С1) + |
25, агс1§ |
- 1 2 1 п Ц |
= |
||
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
• |
агс1§ |
- |
1п |
1 + |
-р- |
+ 2Ь, агс*г |
. |
|
Аномалия силы тяжести в случае прямоугольной призмы бесконечного простирания, две грани которой лежат в плоскостях ху и уг, в точке О будет вычисляться по формуле
Д<^ = /Дб >1п |
Ьг + 2^3 |
апЛЯ' |
С2 |
|
|
|
27* |
403." |
Введем обозначение |
|
|
|
|
|
|
|
|
|
||||
тогда До^ предстанет как функция |
|
|
|
|
|
|
|
||||||
|
|
|
А0д = 2/ Дб^2 [ 4 - Р 1п (1 + |
± |
) + |
агс1В р ] . |
|
|
|
||||
Можно это записать короче |
|
|
|
|
|
|
|
||||||
где |
|
|
|
|
|
Д0^ = Д б ^ , |
|
|
|
|
(XIЛ'.З. |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Р = 2/ |
|
р 1п (1 + ± ) + |
агс18 р ] • |
|
|
|
|||
Последнее равенство позволяет искомое значение |
аномалии для |
призмы |
|||||||||||
с сечением МИР(1 выразить таким образом: |
|
|
|
|
|
|
|||||||
|
|
|
|
Аё = Дб |
(Рм - Рк) + & (Рр - |
Ро). |
|
|
|
|
|||
Венинг-Мейнес составил удобные для пользования таблицы, в которых |
|||||||||||||
по аргументу р |
даются значения функции Р. Выдержки из этих таблиц приво- |
||||||||||||
дятся |
в табл. 44. |
|
|
|
|
|
|
|
|
Т а б л и ц а |
|||
|
|
|
|
|
|
|
|
|
|
|
|
||
V |
0 |
0,1 |
0,2 |
0,3 |
0,4 |
0,5 |
|
0,6 |
0,7 |
|
0,8 |
0,9 |
|
0 |
0 , 0 0 |
4 , 4 1 |
6 , 9 8 |
8 , 8 8 |
1 0 , 3 6 |
1 1 , 5 5 |
1 2 , 5 3 |
1 3 , 3 4 |
|
1 4 , 0 2 |
1 4 , 6 0 |
||
1 |
1 5 , 1 0 |
1 5 , 5 3 |
1 5 , 9 1 |
1 6 , 2 4 |
1 6 , 5 3 |
1 6 , 7 9 |
1 7 , 0 2 |
1 7 , 2 3 |
|
1 7 , 4 2 |
1 7 . 5 9 |
||
2 |
1 7 , 7 5 |
1 7 , 8 9 |
1 8 , 0 2 |
1 8 , 1 4 |
1 8 , 2 5 |
1 8 , 3 5 |
1 8 , 4 5 |
1 8 . 5 4 |
|
1 8 , 6 2 |
1 8 , 7 0 |
||
3 |
1 8 , 7 7 |
1 8 , 8 4 |
1 8 , 9 0 |
1 8 , 9 6 |
1 9 , 0 4 |
1 9 , 0 7 |
|
1 9 , 1 2 |
1 9 , 1 7 |
|
1 9 , 2 2 |
1 9 , 2 6 |
|
4 |
1 9 , 3 0 |
1 9 , 3 4 |
1 9 , 3 8 |
1 9 , 4 2 |
1 9 , 4 5 |
1 9 , 4 8 |
|
1 9 , 5 2 |
1 9 . 5 5 |
|
1 9 , 5 8 |
1 9 . 6 0 |
|
5 |
1 9 , 6 3 |
|
|
|
|
|
|
|
|
|
|
|
|
П р и м е ч а н и е . |
Чтобы получить Д§ в миллигалах, |
в формуле (Х1У.34) |
необходимо |
выраж..-: |
|||||||||
в километрах. |
|
|
|
|
|
|
|
|
|
|
|
|
|
Форма кривых для Ад существенно зависит от соотношения между раз- |
|||||||||||||
мерами граней призмы. На рис |
85 приведены кривые аномалии силы тяжести |
||||||||||||
|
|
|
|
|
|
над призмой с |
глубиной |
залегания верх- |
|||||
|
|
|
|
|
|
ней грани |
10 м и |
размерами граней 2'.1 |
|||||
|
|
|
|
|
|
и 200 м. Аномальная плотность Аб при- |
|||||||
|
|
|
|
|
|
нята равной единице. В первом случае |
|||||||
|
|
|
|
|
|
(кривая 1) грань шириной 200 м лежит |
|||||||
|
|
|
|
|
|
горизонтально, |
во |
втором (кривая 2) — |
|||||
|
|
|
|
|
|
вертикально. |
|
|
|
|
|
||
|
|
|
|
|
|
4. В е р т и к а л ь н а я |
с т у п е н : |
||||||
|
|
|
|
|
|
(сброс). Вертикальная ступень (рис. 86, (. |
|||||||
|
|
|
|
|
|
может |
рассматриваться как |
бесконечная |
|||||
|
|
|
|
|
|
призма* у которой одна из вертикаль- |
|||||||
|
|
|
|
|
|
ных граней отодвинута в бесконечности |
|||||||
|
|
|
|
|
|
(рис. 86, б), поскольку и в том, и в дру- |
|||||||
|
|
|
Р и с . |
8 5 |
|
гом случае аномалия силы тяжести будет |
|||||||
Рассчитаем |
аномалию силы |
иметь одну и ту же величину, |
|
||||||||||
тяжести |
для |
вертикальной |
ступени |
вдол: |
|||||||||
оси х (рис. |
87). Принимая у |
2 |
= 0, Аб = |
б |
|
|
|
|
|
|
|||
404.
и пределы интегрирования по |
от 0 до |
|
по т): от — д о + 0 0 |
и по от |
||||
1,-у до |
на основании формулы (Х1У.11) получим |
|
|
|||||
|
|
4 |
А6 |
I |
Щ |
и |
|
|
|
|
«=>оО1- о»о Ехг И |
"йг: |
|
||||
|
|
|
|
1-со+оо |
|
|
|
|
После выполнения интегрирования и подстановки пределов |
|
|||||||
Д^ = / Дб |
1п |
+ 2 |
|
+ агс1§ |
) - |
^ |
+ а г с 1 ^ ) ] } • |
(Х1У.35) |
Р(х,0,0)
" в, ц§§§
шшшш
+ 0шши;
1
Рис. 86
Формулы для вторых производных потенциала имеют вид
Г ^ / Д б Ы ^
7,д = 7,2г = 2/Дб(агс^^ -
На рис. 86 схематически показан ход кривых А§, Тхг, Тгг над вертикальной ступенью.
Положим в формуле (XIV. 35) х — х = 0, х
Тогда
при х = — оо |
|
при ж = 0 |
Д# = /яДб(^2 — |
при а; = + оо |
Д^ = 2/я Дб — ^1)- |
Оказывается, что в бесконечно удаленных точках, а также в точке, расположенной над краем сброса, значение аномалии силы тяжести зависит только
от величины |
— |
и избыточной плотности Дб пород, образующих сброс, |
|||
и не зависит в отдельности ни от |
ни от |
Для двух ступеней с одинако- |
|||
выми значениями разности |
— ^ |
и Дб, но с различными значениями ^ и |
|||
4 0 5
т.е. расположенных на различной глубине от земной поверхности (рис. 88),
вуказанных выше точках соответствующие этим ступеням кривые Ад будут пересекаться. Более детальное исследование формулы для аномалии силы тяжести показывает, однако, что во всех остальных точках кривые совпадать не будут. Именно, более близкой к поверхности ступени будет соответствовать кривая с более крутым наклоном, а ступени, расположенной на большей глубине, кривая с более пологим наклоном.
Для двух точек с координатами +х и —х, расположенных по обе стороны от края ступени, если х 1,х% можно написать приближенное равенство
Ад ( + х ) - Ад (~х) ^ 2я/ Аб - |
(Х1У.36) |
т. е. разность аномалий приблизительно пропорциональна разности глубин
§ 90. Р Е Ш Е Н И Е О Б Р А Т Н О Й З А Д А Ч И Г Р А В И М Е Т Р И Ч Е С К О Й Р А З В Е Д К И М Е Т О Д А М И ,
О С Н О В А Н Н Ы М И Н А Р Е Ш Е Н И Я Х П Р Я М О Й З А Д А Ч И
Полученные в предыдущем параграфе формулы, позволяющие решить прямую задачу, используются в гравиметрической разведке для решения обратной задачи при помощи специально разработанных методов. Среди методов могут быть названы: метод «характерных точек», метод подбора и др., применение которых возможно только в случае, если известны некоторые данные об искомом аномальном теле, как, например, его форма, аномальная плотность, глубина залегания. Эти данные можно получить либо по материалам геологических изысканий, проводившихся в исследуемом районе, либо в результате применения других методов геофизической разведки.
1. М е т о д « х а р а к т е р н ы х т о ч е к » . В случае, когда на основании каких-либо материалов сделано более или менее обоснованное предположение о форме аномального тела, могут быть использованы формулы, дающие значения Ад, ТХ2 и других производных, в виде тех или иных функций координат
точки наблюдения.
Так как наиболее часто значения этих величин рассматривались только для точек, расположенных на оси х, то можно положить, что
406.
Написанные равенства можно рассматривать как уравнения некоторых кривых. Для ряда тел, имеющих простую геометрическую форму, такие уравнения известны. На этих кривых легко отметить некоторые «характерные точки», например: 1) максимума и минимума, 2) перегиба, 3) точки, в которых значение ординаты у равно половине, четверти и т. д. максимального значения орди-
наты кривой, 4) точки, в которых ордината равна нулю, и т. д. Для всех этих точек можно найти соответствующие значения абсцисс х и, зная аналитические выражения для (х), /2 (х), . . ., возможно получить соотношения, связывающие абсциссу х для какой-нибудь из этих точек с элементами, характеризу-
ющими размеры, массу тела и его положение. Такие соотношения и позволяют решить обратную задачу для тех или других частных случаев.
Рассмотрим несколько таких случаев.
Пусть по каким-либо соображениям принято, что аномальное тело по своей
форме близко к шару. В этом случае решение обратной |
задачи может быть |
||||||
произведено следующим образом. |
|
|
|
|
|
|
|
Напишем уравнение кривой |
= |
|
(х), |
график которой дан на рис. 79. |
|||
Поскольку в данном случае г = |
]/х2 |
+ |
|
то в соответствии с формулой |
|||
(Х1У.17) это уравнение будет иметь вид |
|
|
|
|
|||
Л , _ |
/ М а С о |
_ |
№ а |
( л |
|
||
* |
( * в + с . Т ' |
|
?02 |
I |
?о2 |
|
|
Очевидно, что |
§ будет являться |
максимальным значением аномалии^ |
|||||
которого функция Ад достигает при х = |
0. |
Введя обозначение |
|||||
|
= |
|
|
|
|
&0 |
(Х1У.37) |
получим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Дйпах _ |
/ 4 |
+ |
|
|
' • # |
( X I V . 3 8 ) |
Для удобства вычислений выгодно принять это отношение равным некоторому целому числу п. После этого на графике функции Л^ = (х) следует определить абсциссу «характерной точки», аномалия которой Ад удовлетворяет
условию
А?
#
Тогда уравнение (XIV. 38) перепишем в виде
Получили соотношение, связывающее абсциссу х «характерной точки»
с глубиной залегания центра шара |
решим его относительно |
2 |
|
И |
|
Бо |
|
407
Положив п = 2, найдем |
|
|
|
|
|
|
Б0 = * |
, |
/ , = 1,305а:./,, |
|
(Х1У.39) |
где под XI/г понимается абсцисса, |
при которой аномалия |
Д^ равна |
половине |
||
максимального значения Адтах. |
|
(XIV. 39), по формуле |
(XIV. 37) |
находим |
|
Определив |
из соотношения |
||||
аномальную массу шара |
|
|
|
|
|
|
М а = |
С' . |
|
(Х1У.40) |
|
Если еще известна аномальная плотность
Дб = ба —бь
где б3 — плотность шара; бх — плотность окружающей породы, то можно найти радиус шара
Я - " ] / |
Ма |
У4я Аб •
Зная К и найдем глубину к залегания поверхности аномального тела
|
Решим теперь обратную задачу для шара, используя график кривой Тхг = |
|||
= |
/2 (х), изображенный на рис. |
79. |
Уравнение этой кривой в |
соответствии |
с |
(XIV. 20) имеет вид |
|
|
|
|
Тхг = - |
т |
а 1 ^ т 7 . |
(Х1У.41) |
Исследование уравнения (Х1У.41) показывает, что максимальное значение функции Тхг будет при а; = — (1/2) а минимальное — при х = + (1/2) По-
следовательно,
т. е. глубина центра шара равна расстоянию между точками максимума и минимума кривой Тхг. Подставляя в (XIV.41) значение х = —(1/2) Но> найдем соотношение между (Тхг)тах и аномальной массой Ма
|
|
3 |
I2 |
|
|
|
|
|
|
Ьо |
|
|
|
|
|
(Тхг)тах = |
Х-— - |
Ео = , |
" |
=3/Ма |
/ |
Ч1/- = |
0,858-^, |
|
|
( 4 + Й Г |
• |
|
|
|
|
откуда
Ма = 0.0175Ц (Тхг)тах.
Если аномальное тело имеет форму горизонтального кругового цилиндра бесконечного простирания, то уравнение кривой Ад = (х) в соответствии
с (XIV. 27) напишется
= |
(Х1У.42! |
408.
Очевидно, что максимум функции будет иметь место при х = 0, т. е.
|
^ |
= |
( Х 1 У . 4 3 ) |
Положим, что А§тах/А§ = |
п, тогда (Х1У.42) представим как |
|
|
|
|
х* |
(XIV.44) |
|
|
1 + ^ - , |
|
|
|
Т2" |
|
|
|
ьо |
|
где под х понимается |
абсцисса |
«характерной точки», аномалия которой А§ |
|
в п раз меньше А@тах. |
Эта абсцисса находится по графику функции |
Д^ (см. |
|
рис. 81). Соотношение (Х1У.44), связывающее абсциссу х с глубиной залегания
центра цилиндра позволяет определить
При п = 2, |
= XI/2, где хчг — обозначает абсциссу «характерной точки», |
|||
при которой аномалия Д^ = 1/2 Л&пах. Определив, таким образом, |
нахо- |
|||
дим X из формулы |
(XIV.43) |
|
|
|
|
Я = - Ц ^ = 0 , 0 0 7 4 9 ^ т а х . |
|
||
Решим эту же задачу, используя |
график кривой Тхг = /2 (х) (см. рис. 81). |
|||
Уравнение кривой дается формулой (XIV.28); максимум имеет место при |
||||
а минимум при |
|
|
V з |
|
|
|
Но |
|
|
|
V • |
= |
|
|
|
^т 1 п |
|
Уз |
|
Обозначив |
|
|
|
|
|
|
|
|
|
найдем, что |
^т 1п |
-^тах ~ |
|
|
|
|
|
|
|
Подставляя в (XIV.28) вместо х |
;гтах, получим |
|
||
|
,Т ч |
|
9А |
|
откуда |
*, = 0,0Щ2 |
(Г,2)тах. |
|
|
|
|
|||
Рассмотрим случай вертикальной ступени, особенно часто встречающейся в практике.
Итак, предположим, что аномалия А§ обусловлена действием одной поверхности раздела пород с плотностью б, (верхняя толща) и 62 (нижняя толща).
Допустим, что в точке с абсциссой -\-х (см. рис. 86) глубина указанной
поверхности равна |
а |
в точке с абсциссой —х глубина равна |
Тогда, предполагая, |
что изменение аномалий вызывается лишь рельефом |
|
разделяющей поверхности, можно использовать приближенное равенство (Х1У.36). Допустим, что известны Дб = б2 — б1( глубина контактной поверх-
409.
