APGRV dasturi tomonidan ko'rsatilgan SLAE ni yechishning iterativ jarayonini tavsiflovchi parametrlar
Parametrni belgilash
Parametr xarakteristikasi
Yo'q.
iteratsiya raqami
F2
maydonlarning o'rtacha kvadratiga mos kelmasligi
F1
o'rtacha mutlaq marja
FM
maydonlarning maksimal (modul) farqlanishi
KZ
"bog'langan" nuqtalar foizi
E
.. joriy qiymat £ 0
Eslatmalar:
biri). KZ parametri faqat SLAE ni Seidel usulining moslashuvchan modifikatsiyasi bilan echishda o'rnatiladi.
2). F1 va F2 qiymatlarini aniqlashda "qat'iy" nuqtalarda (c j = 1) A g va Ag * maydonlari orasidagi farqlar joriy qiymat £ 0k ga mos keladi deb taxmin qilinadi.
Dastlabki maydonni ko'rsatish yuzasida S , yoki manbalar ustida joylashgan ixtiyoriy sirt S ' da (gorizontal tekislikda, shu jumladan z = const ) siz A g * maydonini tiklashingiz yoki uning quyidagi transformatorlarini A x qadam bilan hisoblashingiz mumkin. :
• V xz - OX o'qi bo'ylab gorizontal maydon gradienti ; • V yz - OU o'qi bo'ylab gorizontal maydon gradienti ;
• V ^^ xz + Vyz "to'liq gorizontal maydon gradient moduli;
• V zz - maydonning 1-vertikal hosilasi;
• V zzz - maydonning 2-vertikal hosilasi;
• V A = (V yy - V xx ) - tekislik yuzasining egriligi;
• V - tortishish potensiali.
Barcha o'zgartirishlar zichlikni aniqlashdan iborat - NOSTA ekvivalent manbai eng uzoq operatsiyani bajaradigan yagona kompyuter bilan original maydonda amalga oshiriladi . A g maydonidagi ma'lumotlar va mahalliy darajada raqamli balandlik modeli - STI DTM matn (ASCII kodida) yoki ikkilik fayl ko'rinishida ko'rinadi - mate GIS SURFER (grd - fayl). Bu matritsalarning har biridagi elementlar soni 32765 dan, qatordagi sonlar soni 5000 dan oshmasligi kerak.
Magnit maydonning manba bo'yicha taxminiy algoritmi A T amalga oshirilgan dastur - IBM - mos keluvchi kompyuterlarda foydalanish uchun mo'ljallangan APP dasturi. Dastur Turbo Pascal 6.0 da yozilgan. Hisoblash jarayoni ham ikki - hetapnym, SLAE hal adaptiv modifikatsiyasi D ning Seidel usuli - sheniya SLAE (rod magnitlanish J ) ikkilik faylda bir xil js.bin - STKOM kompyuter diskida saqlanadi.
Js.bin fayl pro sodir har qanday texnik muvaffaqiyatsizlikka axborot samarali saqlanishini ta'minlaydi, har 5 iterasyonlar yoziladi - yon berish ishlari. Js.bin fayl ma'lumotlar ixtiyoriga bilan dasturi tomonidan qayta-o'qib bo'lishi mumkin - hal oydinlik go.
APGRV dasturida bo'lgani kabi, insult hisoblari dis - Lida ko'rsatiladi va bir vaqtning o'zida protokol.txt log fayliga yoziladi.
Asl maydon tayinlash yuzasida S , yoki gorizontal tekislik ustida joylashgan manba z = const magnit maydonini tiklashi mumkin A T yoki calcd -
o'g'ri qadam D x bilan kvadrat tarmoq tugunlari uning transformantlar quyidagi :
• vertikal komponent Z ;
• gorizontal komponent H ; dT dzbirinchi vertikal lotin
d 2 T • ikkinchi vertikal hosila - y ' dz 2 ( dT К dx )to'liq gorizontal gradient moduli
• psevdo-gravitatsion maydon D d .
D T maydoni va raqamli relef modeli DTM haqidagi ma'lumotlar VIRG-Rud-Geophysics tomonidan ishlab chiqilgan LOBAS DBMS formatidagi matnli (ASCII kodlarida) fayllar ko'rinishida taqdim etilgan. DBMS LOBAS yordamida bu format osongina matn yoki ikkilik GIS SURFER fayllariga (grd - fayllar) aylantiriladi. Ushbu matritsalarning har biridagi elementlar soni 32765 dan, qatordagi sonlar soni 3000 dan oshmasligi kerak.
Keyinchalik bir qator modellashtirish va amaliy misollar taqdim etiladi - magnit va tortishish maydonlarining konvertatsiya aniqligini taxmin qilish; egri sirtda qayd etilgan anomaliyalarni gorizontal tekislikka olib kelish maqsadga muvofiqligini tasvirlab bering; filtratsiyani ko'rsatish - nye imkoniyatlarni manbaga yaqinlashtirish va hokazo.
Model misolini ko'rib chiqaylik (17-rasm): Magnit maydon qiymatlari D T vertikal magnitlangan kub J z = 10 A / m o'lchamdagi 400x400x400 m o'lchamdagi hrebtoobda o'rnatilgan - 100 - 1077 m balandlikdagi farqli relyefning turli shakllari Yuqori yuzi. kub dengiz sathida ( z = 0 belgisida ); normal mage - Netic maydon T 0 Noril'sk viloyati hisob-kitoblarda ishlatiladigan parametrlari : moyillik I = 82 °, og'ish D = 20 °.
Ko'rinib turibdiki, rasm magnit maydonni "ko'rmoqda" bu juda oddiy mo - Dehli sirt topografiyasi o'lchovlari ta'sirida sezilarli darajada buzilgan. Bunday holda, noto'g'ri ekstremal D T paydo bo'ladi , dala manbasining kun yuzasiga proektsiyasidan tashqarida joylashgan.
Mahalliy ravishda ajratilgan relyef bo'yicha dalaning fazoviy taqsimlanishining buzilganligi - STI qazib olish jarayonida kerakli ob'ektlarning yo'qolishiga olib kelishi mumkin - aniqlangan anomaliyalar epitsentrida burg'ulash ra - bot va natijada - o'rganilayotgan hududning istiqbollari to'g'risida noto'g'ri xulosalar chiqarish. Markaziy maydon D T gorizontal tekislikda z = -1200 m atrof-muhitning ekvivalent modelidan foydalanib, 39x39 = 1521 vertikal yarim cheksiz novda, qiymatlarni tiklashni ta'minlaydi mage - netic maydon, deyarli bir xil natijalar to'g'ridan-to'g'ri muammoni hal qilish (o'rta - kvadratik 17B-rasm va 17.D-rasmda ko'rsatilgan amplitudali maydonlardagi nomuvofiqlik ± 2 nT dan kam).
Teng ta katta bilan turli darajadagi aylantirish sohalarda to'g'riligiga miqdoriy baholash - değerlikli Grid modellari, shuningdek, boshqa yo'l bilan erishish mumkin, men foydalanish emas - to'g'ridan-to'g'ri muammo hal [3] etkazish natijalari. Bunday baholash bir misolni ko'rib chiqaylik, biz xulosa - sistent quyidagi operatsiyalar keyingi ijro:
z = -250 m tekislikda berilgan magnit maydonining D T 250 ga yaqinlashishi vertikal rodlar maydoni ( F2 ~ ± 0,5 nT).
Magnit maydonni taxminiy tekislikdagi dizayndan tiklash - STI z = -500 m, bundan keyin D T 500 deb belgilanadi .
