Поправка за влияние окружающего рельефа

На результаты измерений с гравиметрами оказывает воздействие избыток или недостаток масс, расположенных вокруг точки наблюдения. Необходимость введения поправки за окружающий рельеф определяется в каждом конкретном случае степенью расчленённости рельефа и необходимой точностью работ. Согласно инструкции по гравиразведке, поправки за влияние рельефа местности вводятся на тех пунктах, где они превосходят 0,5 величины проектной среднеквадратичной погрешности определения аномалии силы тяжести. Как положительные формы рельефа, находящиеся вокруг точки наблюдения, так и отрицательные уменьшают наблюдённое значение силы тяжести. Поэтому поправка за влияние рельефа всегда положительна (рис. 2.25).

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

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

Затем, рассчитав таблицы для различных углов I и используя превышения Ah, определяют поправку за рельеф в центральной зоне.

 

 

g

 

Рис. 2.25. Влияние положительных и отрицательных форм рельефа

на поле силы тяжести

Для определения поправок за рельеф в ближней, средней и дальней зонах массы рельефа аппроксимируют либо набором кольцевых секторов (рис. 2.26а), либо набором вертикальных параллелепипедов (рис. 2.26 б).

а

б

Рис. 2.26. Способы аппроксимации рельефа при вычислении поправки.

При реализации первого способа порядок работы следующий:

1. Вычерчивают на кальке палетку в виде кольцевых зон, разделенных радиальными лучами на сектора (в зависимости от зоны 4, 6, 8 или 16 секторов).

 

Рис. 2.27 Палетка для расчета поправки за рельеф

Рис. 2.28. Палетка П.И. Лукавченко для зоны 1000-1500 метров

2. Например, если пределы ближней зоны от 250 м до 1000 м, вычерчивают две
окружности от одного центра в масштабе топокарты радиусами 250 и 1000 метров и делят их
на 8 зон (рис. 2.27).

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

4. По специальным номограммам П.И. Лукавченко (рис. 2.28) для соответствующей
зоны, зная вычисленные разности высот, определяют гравитационное влияние каждого
сектора.

		All, M	Ag, * 10 ^2мгп
	1000-	-33*
		-32
h,M A	*10'мпл"	- 31
800-		_	" ^3° Ah, M
	-25		- 29
	-24		- 28
	-23	900-	- 27
			- 26
	-22
			- 25
	-21
			- 24
700-	-20		- 23
	-19		- 22
		800-
	-is		- 21
	-17		- 20
			- 19
	-16
			- 18
600-	-15	-
			- 17
		700-
	— 14	_	- 16
	-13	-	- 15
	-12		- 14
500-	— 11		- 13
		600-	- 12
	-10
		-	-11
			- 10
	-8	-
		500-
400-	-7		- s
	-6	-	- 7
	-5	400-	- 6
			- 5
300-	-4
			- A
	-3	300-	- 3
200-	-2		- 2
	-1	200-	_ l
100-		100-
	—0		- 0

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

* = —^. (2.70)

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

Эта формула представляется в виде суммы конечных разностей по кубатурной формуле Симпсона:

(

I» i j (xf+yf)

где S - шаг ЦММ;

xi=iS; yi=y^s - координаты текущей точки ЦММ; h0 - высота точки, для которой выполняется расчет; Vy - коэффициент формулы Симпсона.

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

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

В полевой период составляется предварительная карта аномалий силы тяжести в редукции Буге. При этом для съемок масштаба 1:50000 и мельче поправку за рельеф вводить не обязательно, а для съемок масштаба 1:25000 и крупнее - обязательно на участках сложного рельефа.

Поправки за рельеф местности вводятся в результаты наблюдений на тех пунктах, на которых они превышают 0,5 величины е_а.

Для съемок масштаба 1:50000 и мельче радиус учитываемой области при вычислении поправок за рельеф принимается равным 50 км. Для горных районов, где влияние масс в зонах с радиусом более 50 км превышает 0,5 мгл., радиус учета влияния рельефа устанавливается 200 км.

Таким образом, в результате окончательной камеральной обработки вычисляются в каждой точке аномальные значения силы тяжести по общей формуле:

В.=Ви~-Го+ (°>³⁰⁸⁶ " °>⁰⁴¹⁹ • °)^Н + *,> (2-⁷³)

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

Гравиметрическая карта в редукции Буге Нормальная формула Гельмерта (1909 - 1911 г.) Плотность редуцирования 2,67 г/см^3

Рис. 2.29. Пример оформления гравиметрической карты

Инструкцией предусмотрено обязательное построение карт аномалий Буге с плотностью промежуточного слоя, равной 2,3 г/см³ и 2,67 г/см³ для мелкомасштабных съемок, включая съемки масштаба 1:50000. Карты с плотностью промежуточного слоя, 2,3 г/см³ составляется без введения поправки за рельеф. Для геологической интерпретации строятся карты и графики аномалий Буге с истинной постоянной или кажущейся плотностями промежуточного слоя. В этих случаях поверхность относимости не обязательно должна совпадать с поверхностью геоида.

Основным критерием оценки точности аномалий служит СКП аномальных значений:

_р _ + _p2,,,2222 (21 A\

где £_ряд; s_on.; s_H - пояснены ранее, е^ - СКП определения нормального значения силы тяжести; е_р - СКП поправки за рельеф.