Файл: Богданов, В. И. Вычисление гравитационных аномалий от трехмерных тел (графические способы).pdf

лос и их частей, ограничивающих измененные размеры тела по простиранию, на число площадок двухмерной палетки, цена деления которых уменьшена в два раза (гравитационный эффект от полубесконечных призм). Изложенное поясняет рис. 12. Двух­ мерные эффекты от элементарных призм составляют 2.0 мгл каж­ дый. Трехмерный эффект от призмы длиной 25=6.0 км составляет

2-0.905 = 1.81 мгл, а от призмы длиной 25=10 к м — 2-0.45= =0.90 мгл.

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

(*. ■')

Здесь Рк — количество площадок равного действия но двухмерной палетке, выраженное в миллигалах; dk — масштабный коэффи­ циент, снимаемый с нижней части палетки; q. — коэффициент, определяющий число элементарных телесных углов, под которыми из центра палетки видны размеры тела по простиранию на гори­ зонте Z=4.5 Az. Процесс вычислений повторяется для каждой точки профиля.

Трехмерная палетка Г. А. Гамбурцева может быть использо­ вана также при вычислениях аномалий силы тяжести от тел сложной формы по их горизонтальным сечениям. В этом случае строятся наборы плоских гиперболических палеток для каждого горизонта Az отдельно. Реальное тело аппроксимируется ма­ териальными листами на тех же глубинах, и для каждого листа в отдельности подсчитывается гравитационный эффект аналогично тому, как это описано В. Б. Наугольниковым [33]. Вместо набора палеток можно воспользоваться какой-нибудь одной дляфикси­ рованного горизонта, а при вычислениях на других горизонтах изменять цены делений, определяемые коэффициентом 1/тп (табл. 2).

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

31

На рис. 13 приведен пример распета гравитационного эффекта от горизонтального кругового прямого цилиндра, симметричного относительно плоскости XOZ. Как видно из рисунка, точность вычислений на теоретической модели вполне удовлетворительная.

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

раций.

Рис. 13. Расчет аномалии силы тяжести над горизонтальным прямым круговым цилинд­ ром (<з=1.0 г/см3, 2Ь=20 км) по трехмерной палетке Г. А. Гамбурцева.

j— теоретическая кривая силы тяжести; 2 — зна­ чения силы тяжести, вычисленные по палетке.

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

Некоторые другие способы

Уравнение (45) допускает различные способы деления полу­ пространства на равные телесные углы. Рассмотрим некоторые из них, в частности те, которые позволяют построить наборы плоских палеток.

Разделим нижнее полупространство на равные телесные углы системой коаксиальных круговых конических поверхностей с осью симметрии OZ и с общей вершиной в начале координат. В этом случае на полусфере радиуса R вырезается несколько сферических поясов, высоты которых равны между собой и составляют h=Az. Проведем также через равные интервалы А2 параллельные гори­ зонтальные плоскости. В пересечении конических поверхностей с- горизонтальными плоскостями образуются усеченные коаксиаль­

32

ные конические кольца, создающие в начале координат одинаковый гравитационный эффект.

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

и палетка для уровня Z=4.5 Az приведены на рис. 14,а. По­ строение палеток весьма просто. Для этого необходимо выбрать интервал Az и радиус сферы R. При делении плоскопараллельной пластины телесными углами иа десять частей, каждая из которых создает в начале координат одинаковый гравитационный эффект

ния горизонтальных линий (следов пересечения плоскопараллель­ ных пластин с плоскостью чертежа) с окружностью радиуса R соединяются с началом координат. При вращении полученных треугольников вокруг оси OZ полупространство разделяется кони­ ческими поверхностями на равные телесные углы (46).

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

секториальные площадки. Цена деления одной секториальиой площадки плоской палетки определяется формулой (50), однако в ней необходимо заменить т = 4к/Д шна ;?1=2тг/Дш. На рис. 14, а

т = 10, а 71 = 16, поэтому е=0.25 мгл.

Процесс работы с набором таких палеток описан в литературе [8, 33]. Геологическое тело аппроксимируется горизонтальными материальными листами на глубинах интервалов Дz, а грави­ тационный эффект, создаваемый каждым листом на дневной по­ верхности, определяется как произведение числа секториальных площадок соответствующей палетки, попадающих в контур ма­ териального листа, на цену деления, избыточную плотность и множитель, устанавливающий отличие масштабов изображения тела и палетки. Суммарный эффект от всех материальных листов в данной точке поверхности определяет гравитационное дейст­ вие тела.

Вместо набора палеток можно воспользоваться одной, для ка­ кого-либо фиксированного уровня Z, например для Z=4.5 Дг, а при вычислениях на других Zk или перестраивать масштаб изо­ бражения, что не всегда удобно, или же изменять цену деления в пределах каждого концентрического кольца. Изменение цены деления легко может быть установлено аналитически или графи­ чески из рис. 14 через изменение телесных углов, под которыми на различных горизонтах видны отдельные концентрические кольца на горизонте Z=4.5 Az. В табл. 3 приведены соответ­ ствующие значения цены деления для такой палетки. Расчеты произведены по формуле

где i = l, 2, 3. . . — порядковый номер плоского угла, определяю­ щего соответствующий телесный угол Дсо., а следовательно, и номер концентрического кольца плоской палетки; h. — высота сферического пояса, вырезаемого на сфере радиуса R тем же телесным углом.

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

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

34

Т а б л и ц а 3

Изменение цены деления концентрических колец r 4—r <+j плоской палетки для уровня 2= 4 .5 Дг (рис. 14, а) при расчетах ею на горизонтах Z k, мгл

h

Ti-T-i+l