Tog‘larda gravimetrik va magnitik suzish ma’lumotlarini qayta qilish va tavsib etish uchun kompyuter texnologiyalari

Guruch. 10. Daralar ustidagi yuqori aniqlikdagi er magnit tadqiqoti natijalari. Udmurtiya. (V.V.Brodoviyga ko'ra, 1984).
Eslatma: datchikning er yuzasidan balandligi 2 m; relyefni tashkil etuvchi jinslarning magnitlanishi (12-60) x10 ^-6 birlik. SI.
Yuqori tuzatishlarni hisoblash er va aeromagnit tadqiqotlar ma'lumotlarini ofisda qayta ishlash jarayonida amalga oshiriladigan majburiy operatsiyalarning bir qismi emas. Magnit-qidiruvda relefning ta'siriga tuzatishlarni hisoblash uchun pallet usullari va dasturlari ishlab chiqilgan, ammo amalda juda kam qo'llaniladi [62, 76, 88, 94, 95, 97 va boshqalar].
2.2.Magnit relyefning ta'siri uchun tuzatishlarni hisoblash algoritmi
Hisoblash algoritmi muallif tomonidan ishlab chiqilgan topopopravok 5 T _p [47, 49] vertikal to'rtburchaklar ko'p dan to'g'ridan-to'g'ri magnit ma'lumotlar muammoni hal asoslangan - o'zboshimchalik darajasida z da yuzasi va maqsad o'rtasida kiritilgan jinslarning taxminan qatlamini DAVLAT parallelepipeds. ₀ = const gorizontal tekislik. L x L parallelepipedlar asoslarining kattaligi foydalanilgan DTM ning D x qadamiga mos keladi . Hisob-kitoblar magniyda qabul qilingan to'rtburchak koordinatalar tizimida amalga oshiriladi - torazvedke: geografik shimolga yo'naltirilgan OX o'qi ; OU o'qi - sharqqa; OZ o'qi - vertikal pastga.
DTM yordam balandligi bir Matrix bo'lgan H olgan, m satr va n ustunlar. DTM tugunlari parallelepipedlarning yuqori asoslari markazlariga to'g'ri keladi, ularning ustida bT _{p ni} aniqlash nuqtalari _joylashgan . FIELD hisoblangan - to'g'ridan-to'g'ri zarba relyefi D 2R tomoni bo'lgan kvadrat bilan cheklangan ( R - "mozaika" radiusi). Qiymat hisoblangan nuqta dT _tumani , kvadratning markazida joylashgan.
Birlik parallelepiped D T _p ning anomal ta'siri Yerning T 0 normal magnit maydonidan ancha zaifdir , shuning uchun uni anomal magnit maydon komponentlarining T ₀ vektor yo'nalishi bo'yicha proyeksiyasi sifatida ifodalash mumkin :
(2.1)^A T ^(x , y, z ) = S ⁽ IW _XZ + mW _yz + nW _zz )
bu yerda l, m, n - maydonning yo'nalish kosinuslari T 0 : l = cosI sin D, m = cosI cos D, n = sinI; I, D - vektorning moyilligi va moyilligi T 0 ; W, W _z W _zz - parallelepipedning tortishish potensiali W ning ikkinchi hosilalari ; c - doimiy koeffitsient [19].
Nuqta va markaz hisoblash topopopravki asosi orasidagi masofaga qarab - Nia parallelepiped qiymatlari W _xz , W, W _zz aniq yoki taxminiy - - qo'shilgan formulalardan hisoblanadi; ikkinchi holatda anomal vertikal ta'siri bilan belgilanadi - Nogo yupqa novda, unda hisoblash tezligi 10 dan ortiq marta oshdi va magnit maydon amplitudasi xatosi ± 1 nT dan oshmaydi.
Hisoblash aniqligini oshirish uchun dT _P baland bo'yli dan interpolyatsiya tufayli nuqtasi hisoblash topopopravki yaqinida avtomatik ravishda oshiriladi semiz DTM tarmog'i - teglar yordam H . Bunday holda, "palitra" ning markaziy parallelepipedi L / 3 asoslari tomonlari bilan 9 ta kichik parallelepipedga bo'linadi ; va qolgan tanasi 24, joylashgan - kvadrat hajmi 5L x 5L ichida schiesya - har bir kichik parallelepiped ustida 4, tayanch tomoni L / 2 hajmi . 1.2-bo'limda keltirilgan topopopravki 8 g _p markaz mintaqasining algoritmini aniqlashdan farqli o'laroq (subdomain D 1 ), bu holda yn - balandliklarning interpolyatsiyasi garmonik funktsiyadan foydalangan holda amalga oshiriladi [4].
TRRM dasturi aero- va Nazem - yadroviy magnit tadqiqotlar uchun 8 T _p topopopravok hisoblash imkonini beradi . Petromagnit parametrlari magmatik va metamorfixlar - FIR jinslari yuqori lateral o'zgaruvchanlikka ega; zamonaviy erotik yilda - sek Tor bu qoyalar deyarli nonmagnetic cho'kindi bilan aloqada bo'lishi mumkin kesib - DAVLAT konlarini. Shuning uchun, umumiy holatda relyefning magnit xususiyatlarini tavsiflashda, jami magnitlanish arra tegirmonlari J fazoviy koordinatalarning funktsiyasi - kuchlanish DTM: J = J (x, y) . J vektorining yo'nalishini vertikal deb hisoblash maqsadga muvofiqdir: J = J _z tabiiy to'shakdagi jinslarning magnit xususiyatlarini aniqlashda odatiy holdir. Ma'lumki, izotrop jinslarning induktiv magnitlanishi doimo T 0 maydoni bo'ylab yo'naltiriladi ; tabiiy qoldiq magnitlanish vektorining yo'nalishi
ness jinslari Kaynozoyning yaqin yo'nalishi vektori T ₀ bo - Li qadimgi jinslar oxir-oqibatda dastlabki qoldiq - xlorid magnitlanishini yo'qotishga moyildirlar [76].
Sof matematik nuqtai nazardan, modulning kattaligi bilan emas, balki J vektorining yo'nalishi bilan ham ishlash qiyin emas, uni har bir birlik hujayra ichida har xil bo'lishini o'rnatish. Juda qiyin geologik-muhit petromagnetic kutish xususiyatlarini batafsil bayoni sifatida geofizik tadqiqotlarning amaliyotda - ronefroz. Aks holda, [97] da hazil bilan ta'kidlanganidek: Institutning ba'zi miqdori - jinslarning jismoniy xususiyatlarini shakllantirish "ma'nosiz gravirazvedochnye (magnitorazvedochnye) operatsiyasi bo'lib, uning maqsadi zich taqsimlanishini o'rganishdir - STI (magnitlanish). er yuzida."
Topopopravki dT _p er usti yoki aeromagnit (AMC) sharoitlari uchun hisoblanishi mumkin:
• berilgan parvoz balandligi matritsasi uchun H _qavat ;
• tadqiqot parvozlarining doimiy haqiqiy balandligida;
• to'liq relyef oqimi bilan AMS sharoitlarini simulyatsiya qilishda.
Nuqta balandligi H _t gorizontal koordinatalarini (x, y) , hisoblash bo'lgan d t _p amalga oshiriladi , quyidagi munosabatlar bilan belgilanadi:
• yerni tekshirish uchun: H _t (x, y) = H (x, y) + 1 m ;
• To'liq relyefni o'rash uchun AMC: H _t (x, y) = H (x, y) + H (x, y) da oldindan bo'yalgan - ritelno tayyorlangan matritsa H _qavat ;
• doimiy haqiqiy balandlikda AMS uchun: H = H = barcha poyga qismi uchun const - juft 5 T ;
R'
• relyef atrofida to'liq oqim bilan AMS sharoitlarini simulyatsiya qilishda: H _t (x, y) = H (x, y) + H _qavat , bu erda H _qavat = butun hisoblash bo'limi b T uchun const.
Yuqoridagi ifodalarda H (x, y) - koordinatalari (x, y) bo'lgan nuqtadagi relefning mutlaq balandligi .
TPPM dasturlash tilida Turbo Paskal 6.0 yozilgan dastur, oldindan - to'plami harakat IBM-mos keluvchi kompyuterlar foydalanish uchun tayinlangan - tizimli yo'l operatsion tizimi MS-DOS (yoki musobaqa rejimi MS-DOS bilan). Dastur DIAL - vy operatsiyasiga ega, u yuqori tezlik va ishlov beriladigan ma'lumotlar miqdori bo'yicha amaliy - aniq chegaralarning yo'qligi bilan ajralib turadi .
Tuzatishlarni hisoblashning namunaviy misolini ko'rib chiqing dT _p : DTM ajratilgan Rel - ephah m = 100 qatorni, n = 100 ustunni o'z ichiga oladi ; tarmoq 200x200 m; er balandligi belgisi 100 - 1085 m oralig'ida o'zgarib turadi, o'rtacha 418,4 m va o'rtacha - standart og'ish 3 m ± 266. (11-rasm). TPPM dasturi bilan har xil turdagi magnit yozuvlar uchun kT _p topografik o'zgarishlar hisoblab chiqilgan - Yer, 250 m balandlikda er ustida to'liq oqim bilan AMC va 2000 m doimiy bosim balandligidagi AMC tortishish parametrlari haqiqiy amalga oshirishga mos keladi. Norilsk hududida magnitorazvedochnyh ish turlari. Qabul qilingan o'zgarmas lateral vertikal magnitlanish arra tegirmonlarini hisoblashda J = 3,5 A / m marta - yuqoriga mozaika 4,4 x 4,4 km edi; natijalar rasmda ko'rsatilgan. 12 va tab. 3.
3-jadval