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

APP dasturidan foydalangan holda magnit maydon balandligiga qayta hisoblashning aniqligini miqdoriy baholash bo'yicha hisoblash eksperimenti natijalari

Magnit maydon	Statistik parametrlar, nT
	Eng kam	Maksimal	O'rtacha	RMS
D T ^250	137.8	263.3	-0,9	± 75,0
D T ^500	-83,8	173.6	-0,1	± 49,4
d T ^250	-151.3	266.2	-1.2	± 76,2
D T ^250va D T ^250maydonlari orasidagi farq	-16.7	18.4	0.3	± 2,5

Shaklda keltirilgan ma'lumotlardan dalolat beradi. 18 va tab. 13, qayta aniqligi - magnit maydon hisob D T turli darajalarda juda yuqori.
Foydalanilgan kompyuterga asoslangan geofizik maydonlarni o'zgartirishning an'anaviy usullariga nisbatan manba jihatidan yaqinlashuvning afzalliklari - Nuh matematikasining raqamli integratsiya va differentsiallash usullari:
| i | = n | j | = n
^V ( X ^, y ^, z ) = XYo CU j (3.17)
i = 0 j = 0
Bu erda U .. - gorizontal tepada joylashgan kvadrat tarmoq tugunlaridagi asl maydon qiymati - NOSTA va koeffitsientlar C .. transformatsiya yadrosi K ga bog'liq :

Om 5000m 10000m
Guruch. 18. Turli darajadagi (balandliklarda) D T magnit maydonini qayta hisoblashning to'g'riligini baholash .
A - 250 m balandlikda qayd etilgan dastlabki magnit maydon; B - A maydonini 500 m darajaga qayta hisoblash natijasida olingan magnit maydon; B - magnit maydon B maydonini 250 m darajasida konvertatsiya qilish natijasida olingan magnit maydon Eslatma: A va B maydonlari o'rtasidagi farq kutilayotgan 0,3 nT va srednekvad bilan tavsiflanadi - raqobatbardosh og'ish ± 2,5 nT deb hisoblanadi.
.. x i + D x / 2 y i + A x / 2
(3.18)C j = - jj K ( x -% , y -n , z ) d & n ^2 P x i -A x / 2 y i -A x / 2
namunaviy misolni ko'rsatadi. Prizmaning tortishish maydonida D d (19-rasm) tasodifiy shovqin 8D d o'rnatildi , [-0,25 mGal, 0,25 mGal] oraliqda bir xil taqsimlangan; keyin z = -500 m tekislikda V maydonning birinchi vertikal hosilasi ikki xil usulda hisoblangan (20-rasm).
Hisoblash maydonidagi prizma maydoni 0,01 dan 2,96 mGal gacha, o'rtacha qiymati 0,30 mGal va RMSD ± 0,54 mGal bo'lgan. Kuzatish yuzasida D d prizmasining model maydoni va A \ g aralashuvi bilan murakkablashgan bir xil maydon o'rtasidagi standart og'ish ± 0,520 mGal edi; model maydoni D d va qayta tiklangan maydon D d * o'rtasida - ± 0,146 mGal. Binobarin, D d ^ D d * konvertatsiyasida foydalaniladigan funksiyalarning garmonik tabiati tufayli 8D d shovqinining amplitudasi deyarli 4 barobar kamayadi .
Rasmda ko'rsatilgan misol uchun taxminiy tuzilmaning parametrlarini tanlash protokoli (protokol.txt faylini chop etish). 20, Zaydel usulining moslashtirilgan modifikatsiyasi quyida keltirilgan:
** YAKINLASHTIRISH QURILISHNI TANLASH PROTOKOLI ** Maydon qiymatlarining fayl nomi: po1.e.dtb
Balandlik fayl nomi: gtob.dgb
SLAE yechim usuli: 4

N1	F2 -	= 0,14	F1	0,09	FM =	0,55	KZ -	0	E =	0,27
N 2	F2 -	- 0,28	F1	0,27	FM =	= 0,95	KZ -	92	E =	0,27
N 3	F2	- 0,07	F1	0,05	FM =	- 0,52	KZ -	0	E =	0,09
N 4	F2	- 0,10	F1	0,09	FM =	- 0,42	KZ -	88	E =	0,09
N 5	F2	- 0,10	F1	0,09	FM =	- 0,71	KZ -	94	E =	0,09
N 6	F2	- 0,03	F1	0,02	FM =	- 0,29	KZ -	0	E =	0,03
N 7	F2	- 0,04	F1	0,03	FM =	- 0,20	KZ -	82	E =	0,03
N 8	F2	- 0,03	F1	0,03	FM =	- 0,31	KZ -	93	E =	0,03
N 9	F2	- 0,02	F1	0,01	FM =	- 0,15	KZ -	0	E =	0,00
N10	F2	- 0,01	F1	0,01	FM =	0,06	KZ -	0	E =	0,00
N 11	F2	- 0,01	F1	0,00	FM =	- 0,08	KZ -	0	E =	0,00

Dasturning ishlash vaqti 14 soniya edi
Eslatma: SLAE o'lchami 1681 ta noma'lumli 1681 tenglama edi; Glu - loviya og'irligi 150 m (1,5D x ); e ₀ = 0,03 mGal parametri o'rnatildi.
Xuddi shu masalani Zaydel usullari bilan yechish, gevşeme (relaksatsiya parametri p = 0,5 bilan) va oddiy iteratsiya vaqtni sarflagan holda 11, 8,7 iteratsiyalarda muammoni hal qilish sifatining yaqin yakuniy ko'rsatkichlarini olish imkonini berdi. mos ravishda 26 sek, 20 sek, 19 sek.
Ko'p hisoblash tajribalar natijalari dastur dasturi APP APGRV va Trans hisoblash aniqligini oshiradi, deb ko'rsatadi - formant 5-10 marta yoki undan ko'proq. Bu, bu holda, quyidagilar bilan izohlanadi:
• dastlabki maydon nuqtalari va qayta hisoblash balandligida joylashgan gorizontal tekislik orasidagi masofalardagi haqiqiy farqlar hisobga olinadi;
• balandlikdagi potentsial maydonning zaiflashuvining haqiqiy qonuni qo'llaniladi;
• asl - birinchi maydonning garmonik xarakterini buzadigan shovqinlarni bekor qilish amalga oshirildi ;
• sxemasini hisoblashda reydlar tashqari dala tabiati haqida hech qanday ma'lumot yo'q - STI P uning ish (sahna o'rnatilgan ekvivalent MANBA - taxalluslar yaqinlashtirish sohasida allaqachon maydon ekstrapolyatsiya: U ^ 0 masofada sodir bo'ladi . R ^ ~).
3 

Taxminan km 0,6 milya
3 


Guruch. 19. Namunaviy misol: egri sirtdagi prizmaning tortishish maydoni Ag.
A - kuzatish sirtining relyefi; B - prizmaning Ag maydoni.
1 - prizma konturi (o = 0,5 g / kub sm); 2 - prizma asosining markazi.
Рис. 20. Результаты расчета V призмы на уровне 500 м по гравитационному полю, осложненному случайной помехой 8Ag.
Olingan maydon V z : A - sonli farqlash usuli bilan (TRAP dasturi); B - analitik yaqinlashish usuli bilan (APGRV dasturi).
1 - prizma asosining markazi; 2 - atrof-muhitning ekvivalent mesh modeli.
Eslatmalar: TRAP dasturi tomonidan hisoblangan maydon va aniq yechim o'rtasidagi RMS nomuvofiqligi ± 1,46E; APGRV dasturi tomonidan hisoblangan maydon va aniq yechim o'rtasidagi ildiz-o'rtacha kvadrat farqi ± 0,25 E.

Ом1ОООм2000м

Guruch. 21. V tortishish maydonidan V tortishish potensialini hisoblashda shovqin anomaliyalarini bostirish .
A - tasodifiy shovqin bilan murakkablashgan sharning V _z maydoni ; B - egri kuzatuv yuzasida sharning V maydoni; B - z = - 1200 m tekislikdagi sharning V maydoni.
Shuni qo'shimcha qilish kerakki, yuqoridagi omillarning birinchi ta'siri bo'lmagan taqdirda ham (ya'ni, maydonni to'g'ridan-to'g'ri Terr - Is yuzasiga aylantirishda) potentsial maydonlarni o'zgartirishning taxminiy usulining aniqligi - "dan sezilarli darajada yuqori" klassik" usullar (3.17) (3.18) ...
Shovqin anomaliyalarini kuchli bostirish, shuningdek, A g va A T (V yoki A g) antiderivativlarini hisoblashni ta'minlaydi . Simulyatsiya natijalari guvohi - mavjud (21-rasm) gravitatsiyaviy potentsial V hisoblash asoratlarning intensivligi - xlorid tasodifiy shovqin 8A g (anomaliya amplitudasining 50% gacha yetib borishi) maydon A g shar (G _z . = 0,3 g / cc). ), shovqin darajasining keskin pasayishi va foydali signalning kontrastini tanlashni ta'minlaydi. Shu bilan birga, interferensiyaning buzuvchi ta'siri 100 - 1077 m balandlikdagi farq bilan relyefda qayd etilgan A g maydonini maksimal balandlikdan yuqorida joylashgan gorizontal tekislikka qayta hisoblashda ham saqlanib qoladi , bu potentsialning ma'lum bir siljishida ifodalanadi. anomaliya V sharning markazidan.
Simulyatsiya natijalari metrologik uc uchun - tokoobraznoy taxminan magnitlanish vektor J = J _x + J _y + J _z real magnitovozmuschayuschih ob'ektlar yo'nalishi haqida hech qanday qo'shimcha ma'lumot ko'rsatadi . Xususan, anomaliya A T obliquely mıknatıslanmış obyekti dastur hamjihatlikni dan - Agar tasdiq tuzilishi, vertikal magnit iborat J = J , J = J = 0 Y juda yuqori aniqlik qayta vertikal yarim cheksiz asolarini kafolat beradi - magnit maydon A hisob T (shakl 22.). Ushbu xususiyat setochning amaliy ahamiyati - DAVLAT modellari muhiti aniq.
Keling, amaliy misolni ko'rib chiqaylik: rasmda. 23 magnit maydon isodynamic liniyalari xaritasi ( A T) _bir Sharqiy Sajana haqida olingan oltin yo'nalishlaridan biri - ozod atrofida to'liq oqim bilan 75 m balandlikda 25 000,: Havodagi geofizik tadqiqot ko'lamli 1 integratsiya beriladi . Relyef qismi tog', balandlik diapazoni 440 m dan 1120 m gacha tekislikka havola Z = -1200 m magnit maydoni ( A T) _a , kuzatilgan maydon bilan solishtirganda, oddiy morfologiya va zaiflashgan armatura bilan tavsiflanadi - kommunikatsiyalar mahalliy anomaliyalar shakllari. kunlik yengillik... Balandligi z = - 1200 m yaqin
3 

Рис. 22. Магнитное поле ДТ наклонно намагниченной (Jx = Jy = Jx = 2 А/м) призмы на горизонтальной плоскости z = -1200 м, полученное: путем решения прямой задачи гравиразведки (А); путем пересчета с помощью программы АРР (Б).

1 - magnit prizma; 2 - atrof-muhitning ekvivalent mesh modeli.
Eslatma: A va B maydonlari o'rtasidagi o'rtacha kvadratik tafovutlar ± 1,13 nT.
Рис.23. Аппроксимационное преобразование аномального магнитного поля (ДТ)а. Восточный Саян.

0 км
3 

A - kuzatilgan maydon (DT) izo-dinamikasining xaritasi _a ; B - magnit maydon (DT) _va z = -1200 m tekislikda izo-dinamikaga xarita ; B - z = -1200 m tekislikdagi psevdogravitatsion maydonning Dd izoanomal xaritasi ;

G - relyefning izogips xaritasi;
Guruch. 24. Oltin rudasi mintaqasining potentsial konini turli usullarda o'zgartirish natijalarini taqqoslash. APGRV va TRAP dasturlari.
1000 m balandlikdagi V ₂ dala , hisoblangan: A - yaqinlashtirish usuli bilan, B - A. K. Malovichko usuli bilan; 1000 m balandlikdagi 1-vertikal hosila V _zz , hisoblangan: B - yaqinlashtirish usuli bilan, D - sonli farqlash usuli bilan.
Maydonning izoanomaliyasi: 1 - musbat, 2 - nol, 3 - salbiy; 4 - oltinning konlari va rudalari. hudud relyefining maksimal balandligigacha.
Bu, shuningdek, aylantirish psevdogravitatsionnoe magnit amalga oshirish tavsiya etiladi - texnik oddiygina Asosiy: oshiriladi dolzarb bir tekisligi, kuni o'tmish Fi - lytically magnitlanish ma'lum qadriyatlardan J hamjihatlikni inshootlar - yozildi. Xaritalar psevdogravitatsionnogo maydoni A g magnit maydon xaritalar ko'ra aniqroq beradi, o'rganish territo geologik tuzilishi katta rasmga qarang - Ree, chunki ularning qurilishida ijobiy va salbiy haybatli ta'sirini bartaraf - qo'shni anomaliya domenga. Gravitatsiya anomaliyalarini talqin qilish uchun ishlab chiqilgan barcha usullar A g anomaliyalarini talqin qilishda qo'llaniladi [99].
Yana bir bor ta'kidlash kerakki, yaqinlashtiruvchi to'r modellarini qurishda noharmonik xarakterdagi interferentsiya quyidagi sabablarga ko'ra bostiriladi: on - o'lchash xatosi; geologik muhitning er yuzasiga yaqin bir xilligi; U ning raqamli modellari maydonini shakllantirish jarayonida yuzaga keladigan xatolar . Tajriba shuni ko'rsatadiki, katta kabi individual punktlarida bunday xatolar kattaligi - qoida, geofizik tadqiqotlarning aniq bir necha marta katta. Xususan - NOSTA yuqoridagi amaliy misol uchun maksimal farq mage - Netic maydon A T va A T * kuzatish yuzasida S -387 nT, o'rtacha 1,6 va standart og'ish nT ± 3 nT haqida. Binobarin, qadriyatlarini foydalanish U * dan - xaritalash mahalliy anomaliyalar aralashmaslik tartiblash va talqin sohalarda asl raqamli modellari foydalanish afzallik beriladi U .
Sharqiy Sayanning markaziy qismida o'tkazilgan 1: 25 000 miqyosdagi gravimetrik suratga olish ma'lumotlarining o'zgarishlari rasmda ko'rsatilgan. 24. bu ko'rsatkichni ko'rib u veso yordamida an'anaviy tarzda ijro 1000 m va APGRV dasturi balandlikda uchun og'irlik natijalari dala recalculations yaqqol farq - toymasin oynada yig'ish qiymati qadriyatlarni umidsiz. Shuningdek, Utilized ikkinchi holda - (transformable dala qirlar o'rtacha mutlaq darajasida darajasida joylashgan tekisligiga o'rnatilgan, deb taxmin Rudgeofizika - tanishtirish keng tarqalgan ulan dasturi Virjiniya ishlab z = -600 m). Maksimal poygalar - yurish maydonlari periferiya qismida belgilanadi va modulda farq qilishi mumkin, Meas - 0,5 - 5,0 mGal umumiy maydonning 19,5% zonasida.
Rasmda ko'rsatilgan xaritalarni taqqoslash. 24.V, janob ishonchli guvoh - yaqinlashtirish usuli ko'proq qarama-qarshi ano beradi, deb mavjud - V zz tortishish maydonining Malia vertikal lotin soni farqlash ko'ra - TRAP dasturida amalga vanie. Ko'rinib turibdiki, bu yanada ishonchli bo'lishi mumkin - lekin oltin konlari va hodisalarining joylashishini nazorat qiluvchi sulfidlangan tektonik zaiflashgan zonalar uchun mas'ul chiziqli chiziqlarni ajratish.
Transformatorlardagi farqlarning miqdoriy xarakteristikalari 14-jadvalda aks ettirilgan.
Ushbu materiallarda moslashtirilgan Meto-ning afzalliklarini ko'rsatish mumkin - ha, chiziqli o'lchamni hal qilish, bu holda 3878 noma'lum 3978 tenglamani tashkil etdi, £ 0 = 0,05 mgal bilan: 14-jadvalda keltirilgan protokolni tanlash.
Xuddi shu muammoni Seidel va gevşeme usullari bilan hal qilganda, F2 = 0,04 mGal qiymatiga mos ravishda 15 va 13 iteratsiyadan so'ng erishildi. Yechim vaqti 202 soniya 169 soniya edi. Tanlovning eng sekin konvergentsiya jarayoni oddiy takrorlash usuli uchun olinadi: ish dasturidan 390 soniya o'tgach va F2 = 0,10 mgalning 30 ta takrorini bajaring .
qaror dasturi APGRV SLAE katta hajmi tezligi haqida fikr olish uchun, quyidagi amaliy vazifalarni hal natijalari - chi (25-shakl.) Vologochansky maydoni (Norilsk viloyati) uchun. O'lchovlar boshlang'ich mat - matritsalar 141 qator x 201 ustun edi, ya'ni. 28341 noma'lum adaptiv timsoli Seidel usuli bilan neniya - tenglama 28 341 SLAE yechim o'tkazildi . Shunday qilib ispol'uet - Qo'ng'iroq IBM-mos Pentium-I kompyuter soat operatsiya protsesso-
133 MHz va 16MB RAM (jadval 15.).
14-jadval
Turli usullar bilan olingan gravitatsion maydon o'zgarishlarining statistik xarakteristikalari

Gravitatsiya maydonining dastlabki ma'lumotlari va o'zgarishlari	Eng kam	Maksimal	O'rtacha	RMS
Dastlabki ma'lumotlar: 1. Kuzatilgan maydon	-10,45 mGal	13,05 mGal	0,009 mGal	± 5,62 mGal
2. Relyefning balandlik belgilari	440 m	1045 m	637 m	± 102 m
Maydonni 1000 m balandlikda qayta hisoblash: 3. Taxminan hisoblash usuli	-8,02 mGal	13,04 mGal	0,08 mGal	± 4,85 mGal
4. Raqamli integrasiya	-8,54 mGal	8,41 mGal	0,09 mGal	± 5,11 mGal
5. 3 va 4-maydonlar orasidagi farq	-1,78 mGal	5,01 mGal	0,07 mGal	± 0,53 mGal
1000 m da 1-vertikal hosila: 6. Taxminlash usuli	-5.93 E	35.85 E	0,139 E	± 2,77 E
7. Sonli farqlash	-4.40 E	7.36 E	0,002 E	± 1,86 E
8. 4 va 5 maydonlar orasidagi farq	-2,72 E	28.49 E	0,137 E	± 1,30 E

** YAKINLASHTIRISH QURILISHNI TANLASH PROTOKOLI **
Maydon qiymati fayl nomi: gravn.grd Balandlik
qiymati fayl nomi: vis.grd
SLAE yechim usuli: 4

N1	F2 = 1,10	F1 = 0,78	FM = 9.08	KZ = 0	E = 0,45
N 2	F2 = 0,98	F1 = 0,71	FM = 8.93	KZ = 43	E = 0,45
N 3	F2 = 0,83	F1 = 0,61	FM = 7.60	KZ = 58	E = 0,45
N 4	F2 = 0,68	F1 = 0,54	FM = 6.14	KZ = 70	E = 0,45
N 5	F2 = 0,57	F1 = 0,49	FM = 4.41	KZ = 81	E = 0,45
N 6	F2 = 0,50	F1 = 0,46	FM = 3.04	KZ = 90	E = 0,45
N 7	F2 = 0,46	F1 = 0,45	FM = 2.03	KZ = 95	E = 0,45
N 8	F2 = 0,28	F1 = 0,18	FM = 2.19	KZ = 0	E = 0,15
N 9	F2 = 0,21	F1 = 0,17	FM = 1,55	KZ = 56	E = 0,15
N 10	F2 = 0,18	F1 = 0,16	FM = 1.36	KZ = 80	E = 0,15
N 11	F2 = 0,16	F1 = 0,15	FM = 1.15	KZ = 92	E = 0,15
N 12	F2 = 0,13	F1 = 0,09	FM = 1.34	KZ = 0	E = 0,05
N 13	F2 = 0,09	F1 = 0,07	FM = 0,90	KZ = 38	E = 0,05
N 14	F2 = 0,07	F1 = 0,06	FM = 0,60	KZ = 64	E = 0,05
N 15	F2 = 0,06	F1 = 0,05	FM = 0,50	KZ = 80	E = 0,05
N 16	F2 = 0,05	F1 = 0,05	FM = 0,38	KZ = 90	E = 0,05
N 17	F2 = 0,05	F1 = 0,05	FM = 0,30	KZ = 95	E = 0,05
N 18	F2 = 0,07	F1 = 0,05	FM = 0,56	KZ = 0	E = 0,00
N 19	F2 = 0,05	F1 = 0,04	FM = 0,41	KZ = 0	E = 0,00
N 20	F2 = 0,04	F1 = 0,03	FM = 0,27	KZ = 0	E = 0,00

Dasturning ishlash vaqti 110 soniya edi.

БРис. 25. Трансформация потенциального поля (программа APGRV). Норильский район.

A - boshlang'ich maydoni Ag; B - V ning 1-vertikal hosilasi ; B - tortishish potensiali V.
Maydonning izoanomallari: 1 - ijobiy; 2 - nol; 3 - salbiy.
Характеристика итерационного процесса подборафизических параметров аппроксимирующих поле источников
Izohlar: V _zz va V transformatorlari gorizontal tekislikda z = -250 m hisoblanadi.

Takrorlash raqami	Muammoni hal etish sifati ko'rsatkichlari
	F2, konv. birliklar	F1, konv. birliklar	FM, konv. birliklar
bitta	0,87	0,73	5.5
5	0,55	0,52	3.77
10	0,16	0,12	1.20
15	0,06	0,05	0,38
o'n sakkiz	0,03	0,02	0,23

15-jadval

18 iteratsiyadan so'ng u kuzatilgan va modelning yuqori darajadagi bir-biriga mos kelishiga erishdi - Oyoq maydonlari; masalani yechish vaqti 6 soat 32 minut 22 soniya edi.

3.3. Noto'g'ri tarmoq tugunlarida berilgan tortishish va magnit maydonlarining manbaga o'xshash yaqinlashuvi

Muammoning bir umumiy bayonot uchun sourcewise taxminiy o'z ichiga oladi - murojaat manbai salohiyati maydon qadriyatlar U non-gorizontal yuzalar = z ustida z ( x, y ) tartibsizlik, panjara tugunlarida. Amalda, yagona saytlar tarmog'ida maydonning pol - chit qiymatini geofizik o'lchash sohasi juda qiyin, shuning uchun U matritsasi har doim yangi sintezlangan tasvirlarning aniqligini aniqlaydigan ma'lum bir interpolyatsiya algoritmidan foydalangan holda shakllantiriladi. maydonining U. Aynan shu vakolatxonalardir - Lenia keyin izohlash va xaritalash uchun ishlatiladi. Manbaga ko'ra ann - roksimatsiya oldingi bo'limda ko'rsatilganidek , U maydonini nafaqat turli xil o'zgarishlarni amalga oshirishga, balki po- qiymatlarini tiklashni ham amalga oshirishga imkon beradi.

muntazam tarmoq tugunlarida potentsial maydon U * . Bu interpolyatsiya usuli uchun asosiy talablarning bajarilishini ta'minlaydi - "maydonning tasodifiy komponentlari ta'siriga nisbatan barqarorlik va tarmoqqa kirish ma'lumotlarini aniqlash - natijalarning hozirgi aniqligini o'zgartirish" [5].
Quyida geopotenni manba bo'yicha yaqinlashtirishning asl algoritmi - ijtimoiy maydonlar, ularga interpolyatsiya qilish va kesishda - heterojen tarmoqlarga o'tkazish va ko'p sonli o'lchash nuqtalarini (100 000 yoki undan ortiq) kuzatish imkonini beradi . Yaratilgan algoritm V.I. tomonidan taklif qilinganni ishlab chiqadi. Aronovning inter - polatsiyaviy garmonik funktsiyalariga yondashuvi [4, 6], lekin ilgari hosil bo'lgan bo'shliqdan farq qiladi - toymasin oynada o'rtacha hisoblash va moslashtirilgan protsedura yordamida dala komponentlarini bosqichma-bosqich chiqarish, taxminiy dizaynning geometrik parametrlarini aniqlash. Taxminlash muammosi SLAE sabab bo'lgan horo - sho ketma-ketligini hal qilish orqali xarakterli jismoniy - Shomil ko'pligi anomalieobrazuyuschih jismlarini aniqlashga qisqartiradi . Dastlabki potentsial maydon muammosini hal qilish jarayonida U , M mintaqasida tasodifiy joylashgan nuqtalarda sobit P , taxminan - Xia maydoni U * kuzatuv yuzasi ostida joylashgan manbalar.
Ekvivalent manba sifatida bu holda sferalar (Yagona-AF - nye og'irligi) apriori berilgan fazoviy tartibga ega bo'lgan q.s. - U va U * maydonlarining dimaya yaqinlik darajasiga ularning fi - fe- parametrlarini p (zichlik) aniqlash orqali erishish mumkin. yoki magnitlanish) chiziqli - xlorid teskari masalani yechish orqali. Muayyan murakkablik pro - fazoviy joylashuvning aniq ta'rifidir (koordinatalar {E, Z, n } markazlar) ammo ikkita qarama-qarshi talabni hisobga olish kerak:
| Chiziqli tarmoq muhiti modeli (3.12) uchun teskari muammoni barqaror hal qilish faqat manba chuqurligi va tugunlar orasidagi masofa 2 martadan ko'p bo'lmagan farq qiladigan holatda mumkin (bu nisbat butun P maydoni uchun saqlanishi kerak , matritsaning aks holda ahvol soni original sohasida va jismoniy bug'laridan o'rtasidagi munosabatlarni belgilashda buyon - keskin metr manbalari ortadi);
| ekvivalent manbalarning joylashuvi chuqurliklari H maydon anomaliyalarining avtokorrelyatsiya radiuslari bilan qattiq bog'liq U , ya'ni. samarali mintaqaviy interpolyatsiya qilish uchun - xlorid olgan qismlariga nisbatan katta qiymatlari zarur N , va albatta - Kalnoy - kichik.
Bu talablarni qondirish uchun S 1 , S 2 , S 3 , ... turli darajalarda joylashgan dekompensatsiyalangan - zitsiyu muammosi - manbalarni o'tkazish , ularni ketma-ket S yuzasiga olib chiqish va yuqori chastotali barcha komponentlarni tanlash tavsiya etiladi. iterativ jarayon davomida maydon. Anomal ta'sir U * hamjihatlikni con - qurilish tomonidan taqdim etiladi:
k N
u * (x, y, z ) = £ i P j ^ph y (3.19)
i = 1 j = 1
bu erda k - sharlarni joylashtirish darajalari (chuqurliklari) soni; F - mos ravishda o = 1 g / sm ³ yoki J = 0,01 A / m da birlik sharining magnit yoki tortishish maydoni U * ni aniqlaydigan harmonik funktsiyalar ; p - massa yoki IDMAGN - ness maydonlari bilan chiziqli bog'liq bo'lgan parametrlar; N = N (i) - i- darajadagi sharlar soni : i = 1,2,3, ..., k .
Masalan, Ag tortishish maydoni uchun:
F = ( Z - z ) / [« - x ) ² + ( P - Y ) ² + ( Z - z ) ² ] 32 ⁽³ - ²⁰⁾
Muammoni parchalash uchun, shuningdek , M 1 , M 2 , ..., M k anomal ta'sir U * bilan bog'langan shartli ravishda aniqlanishi mumkin bo'lgan asl maydon nuqtalarining kichik to'plamini tanlash kerak - lennymi darajalari ekvivalent manba joylari N 1 , N 2 , ... N k . Keyin yaqinlashish muammosi SLAE ketma-ketligini hal qilish uchun qisqartirilishi mumkin: