CC BY
57
Определение устойчивости территории зон тектонических нарушений методом
спутникового нивелирования
И.Г. Гайрабеков, Ю.И. Пимшин
Грозненский государственный нефтяной институт, г. Грозный
Ростовский государственный строительный университет, г. Ростов-на-Дону
В последние годы радикальным образом пересматриваются взгляды по вопросу влияния геодинамических процессов на эксплуатационную надежность зданий и сооружений. Фактор геодинамической (сейсмотектонической) активности геологической среды, отражающийся на земной поверхности в форме аномальных смещений разной направленности не учитывался при анализе аварийных ситуаций. Равнинные (платформенные) области ранее относились к территориям малой активности, то есть современные вертикальные движения земной поверхности (СВДЗП) рассматриваемых территорий принимались в пределах 1-5 мм в год. На основе экспериментальных данных учеными [1-4] получены результаты, свидетельствующие о том, что имеют место современные суперинтенсивные деформации (СД) земной поверхности со скоростями до 50-70мм. в год, приуроченные к зонам тектонических нарушений. Были выявлены пространственная локализованность СД - процессов (от 0,1км до 1-2км), их знакопеременный и (или) пульсационный характер проявлений (период СД - импульсов составил от 2 раз в год до 1 раза в 2-3 года). Причем, максимальная интенсивность геодинамических аномалий выявлено в зонах асейсмичных разломов в равнинных областях. Выявлен новый класс геодинамических движений в разломных зонах и с периодом 30-60 сек., и 40-60 мин [1]. Экспериментально подтверждено, что всем рассматриваемым геодинамическим движениям, наряду с трендовой составляющей свойственны знакопеременная направленность и пульсационный характер. Так, в зонах тектонических нарушений периодически происходят порывы трубопроводов, аварии на скважинах- более 70% аварийных ситуаций на трубопроводах приурочены к зонам тектонических нарушений (ЗТН) [2,3] Как известно, одним из явных признаков наличия ЗТН являются водотоки. Территории населенных пунктов и промышленных объектов, как правило, тяготеют к таким водотокам, и соответственно являются территория благоприятными для проявления аномальных локализованных деформаций земной поверхности, негативно влияющих на эксплуатационную надежность зданий и сооружений. Если амплитуда рассматриваемых аномальных деформаций земной поверхности (ДЗП) в ЗТН превысит допустимые или предельные деформации сооружений, то это приведет к соответствующим негативным последствиям, в том числе возможно и его разрушению. Следует отметить, что опасность представляют не только сверхдопустимые деформации, но и колебательные деформации меньшей величины, например, с периодом проявления от одной до шестидесяти минут. В таком случае негативные последствия будут связаны с проявлением усталостных эффектов. Результаты геодезических исследований ДЗП в ЗТН легли в основу нормативного документа «Инженерно-геодезические изыскания для строительства» (СП 11-104-97)
В п.10.64 СП 11-104-97 излагается необходимость проведения геодезических наблюдений за ДЗП в районах развития современных разрывных тектонических смещений (РТС). При этом решаются задачи как выявления РТС, количественной оценки и прогнозирования их развития, так и слежения за РТС в период строительства и эксплуатации технически особо сложных и уникальных (I и II уровней ответственности согласно ГОСТ 27751 -88) предприятий и сооружений. Рекомендуется геодезические наблюдения проводить на территории построенных объектов при негативных воздействиях ряда факторов (предположительно тектонической природы) на их устойчивость и надежность.
Причем наблюдения за РТС рекомендуется выполнять как в горных районах, так и
в равнинно-платформенных областях. В горных районах геодезические сети для исследований развития РТС создаются путем локальных плановых и высотных построений по линиям, пересекающим в крест каждое разрывное нарушение. Протяженность таких линий может составлять от сотен метров до нескольких километров. В равнинно-платформенных районах рекомендуется создавать нивелирные построения в виде сплошной сети полигонов с периметром не менее 20 километров.
1. Для выявления и изучения короткопериодических колебаний земной поверхности в ЗТН целесообразнее всего использовать спутниковые методы определения деформаций. Следует отметить, что выявлять и исследовать суперинтенсивные колебания в ЗТН (с периодом от нескольких минут до нескольких часов) только спутниковыми методами и возможно. Если точность определения координат Х и У не вызывает затруднений, то точность результатов спутникового нивелирования остается невысокой для решения ряда научных и производственных задач, в том числе и для наблюдения за деформациями земной поверхности.
В статье [5] приведена формула вычисления геодезической высоты пункта по результатам спутниковых измерений, минуя итерационный процесс вычисления геодезической широты:
Г —2 2 —11/2
Н=|_Х2 +У2 + 22(1-а1Я2г-3)- (1+Ъх2V“3) ] -
п-1/2 , (1)
-а
1-
е2 22 Я ~2
22 Я ~2 +(1-а1Я 2г “3 )2 (1+Ъ122 г -3 )-2
3/2
где а1 = ае2 (1- е2) Ъ1 = Ъе'2;
а и Ъ - полуоси отсчетного эллипсоида;
2 г2
е и е - первый и второй эксцентриситета;
Я=4ХГ+У2 ;
г=^(:1-е2)(X2 + У2)+22 .
Выполним оценку точности вычисления геодезической высоты по результатам спутниковых измерений.
Для оценки точности вычисления геодезической высоты по результатам спутниковых измерений формула (1) представлена в следующем виде:
Н =Р-И, (2)
где N - радиус кривизны первого вертикала, выраженный через прямоугольные координаты;
Р=у1 X2 +У2+22(1-а1Я2г~3)"2(1+Ъх22г-)2 .
(3)
Полный дифференциал выражения (1) имеет вид:
дР дИ дР дИ дР дИ
с1Н=—йХ-—йХ+—йУ-—йУ+—й2-—й2. (4)
дХ дХ дУ дУ д2 д2
Для упрощения дальнейших преобразований введены обозначения
t1 =1-а1Я2 г_3; =1+Ъ122 г_3,
следовательно,
Р=^Х2 + У2 + 2 \%2 .
После нахождения частных производных окончательно полный дифференциал (4) имеет вид:
йН —
Х 3Ъ1(1-е2 )2Х24 2a1t22Х22 _ 3а1 (1-е2)t22Х22Я2 Р „ , „ „ „ ,
ае2 Х22
Pt12г5
Pt13г3
Р^г5
С ( - <3 )3/2 Я
а«3 Х22
4/3 3/2 4 2 4/3 3/2 3 2 4/3 3/2 5 3 4/3 3/2
«4 (4 «3) Я 2 4 («4 «3) г 2«4 («4 «3) г «2«4 («4 «3) г
2aa1t1tз Х + 3aa1t1t3 (1-е2 )ХЯ2 3аЪ^\ «3 Х2
3/2 3
3/2 5
3/2 5
йХ+
+
У_ 3Ъ1(1-е2 )t2У24 2a1t22У22 3а1(1-е2 )<22У22Я2 Р „ , „ „ „ ,
ае 2У22
Р^г5
Р^г3
Р«3г 5
а^У22
4/3 3/2 4 2 4/3 3/2 3 2 4/3 3/2 5 3 4/3 3/2
«4 («4 «3) Я 2 4 («4 «3) г 214 («4 t3) г «2«4 («4 t3) г
2aa1t1t3У + 3aa1t1t3 (1-е2 )УЯ2 3аЪ1«12«3У7
(5)
3/2 3
3/2 5
3/2 5
йУ+
+
2-73
2t2 3Ъ1«22 3а1Я 2
ае2 2
[_ Р«12 Р«12 г5 Р«3г 3 ^ ( - tз )3/2 Я2 «г (t4 - tз ) Я
й2, где t3 =
а«3 2
3/2
3aa1t1t3 2Я
ФГ ( - tз )3'2 г5 ФГ ( - tз ) г5
Воспользовавшись тем обстоятельством, что коэффициенты 1Х и «2 отличаются от
3aЪ1t113 23
3/2
е2 22
единицы на малую величину, член формулы (5) виде:
2«2
—2“^ можно представить в следующем Р tl
2х\ 2 (1-а1Я2 г-3 )2 2 2а12Я2 2Ъ123
Р12«12 Р(1+Ъ12 V3 )2 Р Рг3 Рг3
С учетом (6) формулу (5) можно записать как
(Х V г У Л
= ---+ Л
I Р
йН--
йХ+
—+в
йУ+
(6)
(7)
где
Л = 3Ъ1 (1-е2 )«2Х2 4 + 2а1«2 Х22 3а1 (1-е2 )«22Х2 2 Я2 +
Л — -Т—7--I--~~~—7-------+
ае2 Х22
Р^г5
Р«3г 3
Р«3г 5
«г (t4 - tз )3/2 я 4
2aa1t1t3 Х + 3aa1t1t3 (1-е2 )ХЯ2 заЪ«2^ Х2
at3 Х22
(8)
4/3 3/2 4 2 4/3 3/2 3 2 4/3 3/2 5 3 4/3 3/2 5
¿4 («4 «з ) Я «2 4 («4 «3 ) г 214 («4 «3 ) г «2 «4 («4 «3 ) г
в = 3Ъ1(1-е2 )«2У24 + 2а1«22У22 3а1(1-е2 )«22У2 2Я2 + ае2У22
Р^г5
Р«3г3
Р«3г5
Г ( - «3 ) Я4
2аа1«1«3У + 3аа1«1«3 (1-е2 )УЯ2 3аЪ1«12«3У7
4/3 3/2 4 2 4/3 3/2 3 2 4/3 3/2 5 3 4/3 3/2
«4 («4 «3) Я «2 «4 («4 «3) г «2 «4 («4 «3) г «2 «4 («4 «3) г
3/2 3
3/2 5
3/2 5
2гу3
С=-
2а12Я2 2Ьх2 3Ь1ї22і 3>ахЯ2
ае 2
Рг3
а?3 2
Рг3 Р^2 г5 Р1\г5 ^/3 (4 - ?3 )3/2 я2
3аа1і1і32Я2 3аЬ1?12?323
+
«4/3/7 « У/2 Я2 «2,4/3/« « \3/2 г5 + «2,4/3/« « \3/2 г5
«4 («4 «3) Я «2 «4 («4 «3) г «2 «4 («4 «3) г
Переходя от дифференциалов к конечным приращениям, а от них к средним квадратическим ошибкам, из (7) следует:
\3/2 5 '
т2н =
т2 +
2+с р
\
т
2-
(11)
(У V
—+в
V Р у V р у
Анализ точности определения координат пунктов, используя современные спутниковые измерения, показывает, что в большинстве случаев тХ — тУ — т2 — тк . В
таком случае формула оценки точности вычисления геодезической высоты пункта может быть представлена в виде:
тн = тк,
X2 + У2 + 22 2АХ 257 2С2
+
Р
+
Р
+
Р
+А2 + В2 + С2.
(12)
Анализ формулы (12) показал, что с ошибкой, не превышающей 4%, она может
Х2 +У2 + 22 ,
быть представлена в виде: тН — тк, так как------ ----«1, а остальные члены малы.
Например, в средней полосе В = 45°, Ь= 37, Н=200 м; Х= 3 608 020 м ; У= 2 718 839 м; 2= 4 487 489 м;
X2 + У2 + 22
2 АХ 2ВУ 2С2
= 1,023; -------= -0,0048; ------= 0,0089; ------= 0,0397;
Р2 Р Р Р
А2 = 0,000018; К=1,0237; В2 = 0,000108; С2 = 0,000791.
На экваторе В =0°, 1=37°, Н=200 м; Х=5 093 965м; 7=3 838 578м; 2=0;
X2 +У2+22 2АХ 2 ВУ 2С2 а2 п2 ^2
-= 1; -----= 0;-------= 0; -----= 0; А = 0; В = 0; С = 0.
Р
Р
Р
В Заполярье В =72°, 1=37°, Н=200 м; Х= 1 578 909 м; 7=1 189 793м; 2=6 043 875м;
Х2 +У2 + 22
2 АХ
1,0012; --------= 0,0027;
2ВУ
Р
-= 0,0036; К=1,0593;
2С2
Р
-= 0,1163;
р2 ■ ■ р
А2 = 0,0060; В2 = 0,00009; С 2 = 0,0037.
Несмотря на сложную и нелинейную зависимость геодезической высоты от координат пункта, определенных по результатам спутниковых измерений (1), средняя квадратическая ошибка вычисления геодезической высоты практически не зависит от положения пункта при равноточных декартовых координатах.
Литература
1. Панжин, Сашурин
2. Касьянова, Кузьмин Ю.О. современная аномальная геодинамика недр и ее влияние на объекты нефтегазового комплекса. - М.: «Геоинформмарк», 1996-55с.],
3. Сидоров В. А., Кузьмин Ю.О. Современные движения земной коры осадочных бассейнов. - М.: Наука, 1986.-183с.
4. Клюшин Е.Б., Кравчук И.М. Спутниковое нивелирование. Приложение к журналу Известия ВУЗов. «Геодезия и аэрофотосъемка». Сборник статей по итогам международной научно-технической конференции, посвященной 230-летию основания МИИГАиК. Выпуск 2 (в двух частях). Ч. II. Изд-во МИИГАиК. М.: - 2009, с. 40-42.
