2.3. Методи локалізації та трасування вимірюваного гравітаційного поля
Нам треба виконати трасування поля вздовж розламів під різними кутами під орієнтиром на північ.
Нам потрібно знайти місце, де розлами пересікаються.
Запишемо формулу для пошуку місць пересікання розламу :
Do j = 2.46; Do i = 2.46
(2.3.1)
Якщо ми візьмемо середнє поле по 9 точкам, то в результаті отримаємо карту поля майже геологічного однорідного масиву.
Трасування поля під кутом 45º на північний схід :
(2.3.2)
Трасування залишкового поля :
(2.3.3)
Трасування поля під кутом 45º на північний захід :
(2.3.4)
Трасування поля під кутом 22,5º до північний захід :
(2.3.5)
Трасування поля під кутом 22,5º до північний схід:
(2.3.6)
Залишкова аномалія аналог третьої похідної :
(2.3.7)
(2.3.8)
(2.3.9)
Залишкова аномалия є аналогом другої похідної.
Перераховується поле на висоту за допомогою інтеграла Пуассона перетвореного в кінцеву суму.
(2.3.10)
2.4. Розв’язок задачі за допомогою мов програмування
2.4.1. Робочі елементи програми.
Зовнішній файл в якому зберігається загальна карта поля
25 format((2f8.0,f15.11))
open(11,file='d:\work\temp\testgo11.dat',status='unknown',form='formatted')
read(11,*)ZecZT
rewind 11
print '(1x,i3,1X,2f9.0,f9.4)',((k,(ZecZT(i,j,k),i=1,3),j=1,1),k=1,511,10)
close(11)
Виписування частини файлу для досліджуваної карти поля
do i5=1, 47
do j5=1,47
gtolerg(i5,j5)=ZecZT(3,5*i5,5*j5) !-ZecsigT(3,i5,j5)*0.5 +25.0
enddo ;end do
побудова залишкових аномалій трасованих карт гравітаційного поля під різними кутами до напрямку на північ
do i=1,47
do j=1,47
gtole(i,j)=0.20
end do;end do
do i=3,45
do j=3,45
gtole(i,j)=-(gtolerg(i-2,j+2)+gtolerg(i-1,j+1)+gtolerg(i,j)+ & !NW45
gtolerg(i+1,j-1)+gtolerg(i+2,j-2))/5+gtolerg(i,j)
end do;end do
do i=3,45
do j=3,45
gtole(i,j)=-(gtolerg(i-2,j-2)+gtolerg(i-1,j-1)+gtolerg(i,j)+ & !NE45
gtolerg(i+1,j+1)+gtolerg(i+2,j+2))/5+gtolerg(i,j)
end do;end do
do i=3,45
do j=3,45
gtole(i,j)=-(gtolerg(i-1,j+2)+gtolerg(i,j) & !NW 22.5
+gtolerg(i+1,j-2))/3+gtolerg(i,j)
end do;end do
do i=3,45
do j=3,45
gtole(i,j)=-(gtolerg(i-1,j-2)+gtolerg(i,j) & !NE 22.5
+gtolerg(i+1,j+2))/3+gtolerg(i,j)
end do;end do
do i=3,45
do j=3,45
gtole(i,j)=-(gtolerg(i,j+2)+gtolerg(i,j+1)+gtolerg(i,j)+ & !Nord 0
gtolerg(i,j-1)+gtolerg(i,j-2))/5+gtolerg(i,j)
end do;end do
open(603,file='d:\work\temp\testgo603.dat',status='unknown',form='formatted')
write(603,33) gtole
rewind 603
read(603,33) gtole !Гравика (Х0;У0)=(+10;+10)
rewind 603
!print '(1x,10f7.2)',gtole