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

Topografik tuzatishlar va balandlik belgilari o'rtasidagi munosabatlarning xarakteristikasi

Parametrlar	Rad etilgan ballar,%	Magnit maydon o'lchovlarining balandligi, m
		0	yuz	250	500	750	1000
Koeffitsient	0	0,787	0,747	0,828	0,963	0,972	0,945
korrelyatsiyalar	10	0,863	0,802	0,875	0,973	0,976	0,958
F - mezon	0	8.44	8.56	12.32	30.89	35.67	18.18
	10	15.02	18.47	29.59	61.50	95.37	71.16

Eslatmalar:
bitta) ... Ikki dumli testga ko'ra, korrelyatsiya koeffitsienti statistik ahamiyatga ega bo'lib, 0,321 dan boshlanadi (p = 0,1%; n = 100).
2) ... 1% ahamiyatlilik darajasida va V 1 = V 2 = 60 da F - mezonning nazariy qiymati 1,667 ni tashkil qiladi.
10-jadval
Magnit relyef va balandlik belgilarining ta'siri uchun tuzatishlar o'rtasidagi regressiya munosabatlarini hisoblashning to'g'riligini baholash.

Magnit maydon o'lchov balandligi, m	Rad etilgan nuqtalarning boshqa soni M uchun regressiya tenglamasini tuzishda ildiz o‘rtacha kvadrat xatosi , nT
	M = 0	M = 5	M = 10	M = 20	M = 30
0	97.1	83.2	73.2	57.4	46.5
	75.9	63.8	56.9	44.5	36.7
yuz	86.8	70.6	63.4	52.1	40.9
	62.4	47.7	42.5	34.4	26.2
250	47.8	37.3	34.2	28.7	23.3
	34.1	24.1	22.0	16.9	14.1
500	14.4	13.7	10.7	8.9	7.9
	12.7	10.4	9.0	7.0	5.3
750	8.5	6.9	5.8	4.5	3.7
	7.8	6.1	5.2	4.5	3.9
1000	9.1	7.8	6.4	4.6	3.3
	6.5	5.2	4.6	3.7	2.8

Eslatma: 2-ustunlarda - 6-jadval - hisoblagich taxminiy xatoliklarga bog'liqlikni tavsiflaydi - balandliklar orasidagi ko'prik va 1-darajali relyefli polinom uchun o'zgartirilgan: ST = b 1 h + b _on ; maxraj 2-darajali ko'phad: ST = b 2 h ² + b _t h + b 0 .
10-jadvaldagi natijalar afzalliklari parabolik ko'rsatadi - barcha uchun chiziqli vakillik oldin S T _p va h o'rtasidagi munosabatlarning cal yaqinlashishi - lennogo dala o'lchov oralig'i balandliklar. Qabul qilingan ro'yxatga olishni tasdiqlash uchun - ressionnyh bog'liqliklari va kirish ma'lumotlari ANOVA apparati ishlatiladi.
Jadvalda keltirilgan materiallarni tahlil qilishda. 9, 10 quyidagilar aniq:
• Topopopravkami S T _p va siz o'rtasida juda yaqin bog'liqlik mavjud - yuzlab relyef belgilari h , barcha balandliklarda aslida magnit maydonlarining o'lchovlari (0-1000 m) o'tkazildi. Korrelyatorning chiziqli koeffitsientlarining qiymatlari - lyatsii cho'qqisi 500-750 m balandlikda otish, keyin esa damping tufayli - maydonning balandligi bilan, bu munosabatlar zaiflasha boshlaydi.
• 8 T _p = b ₂ h ² + b ₁ h + b ₀ kvadratik regressiya tenglamasi yuqori tuzatishlar va relyefning balandliklari o'rtasidagi munosabatlarning miqdoriy tavsifi uchun juda mos keladi, buning olingan qiymatlari tasdiqlaydi. RMS va F- mezoni. Tajribada ishtirok etuvchi maydonning 90% nuqtalari uchun 500 m balandlikdan boshlab, kvadratik regressiyaga bog'liqlikni qurishda xatolik AMS aniqligi bilan taqqoslanadigan bo'ladi.
Keyinchalik u dastlabki oldin empirik topilgan - shakllantirish 8 T _p = Lev Nikolaevich 8 T _p hisobga jadal mıknatıslanmış jinslarning [30] ustida magnit maydon qiymatlarining logarifmik normal taqsimlanishini va juft qurish aniqligini oshiradi regressiya funktsiyasi.
Ushbu ma'lumotni P.S. tomonidan berilgan ma'lumotlar bilan solishtirish qiziq. Revyakin, V.V. Brodov, E.A. Revyakina , yuqori aniqlikdagi yerga asoslangan magnit tadqiqotlarda zaif magnit ( J = 0,01-0,02 A / m) relyefini hisobga olish muammosi bo'yicha . Izohlab kuni - D oltin-sulfid konlari tabiiy belgi tasodif kuzatiladi - spektral Vazifalar qirlar qadriyatlar va xato magnit maydonini ( D T) _va past chastotalarda. Past xatolik sifatida o'rnatiladi 8 = ± 0,8 nT to'rt - bu parametrlar o'rtasidagi kvadratik regressiya munosabatlari individual pro - Lyax uchun (maydon amplitudasi ( A T) _va 6-8 nT gacha). Ushbu turdagi qaramlik EC tomonidan taklif qilinadi - relyefning ta'sirini hisobga olish uchun foydalaning [93].
Bo'limda 2.3 axborot, tahlil xatolar dalolat - yordam modeli magnit maydon matic hisoblash, ob'ektlarini har xil - ronefroz omillar lol qoldirardi mavjud, balandligi orasidagi chiziqli bo'lmagan regresyon, munosabatlar qurish xato oshmasligi mumkin relyef va topo - o'zgartirilgan. Bu uchastkaning yuqori qismining magnit xususiyatlarining yuqori lateral o'zgaruvchanligi bilan subvertik magnitlangan jinslardan tashkil topgan relyefning ta'sirini hisobga olishning statistik usulini qo'llashning asosiy shartidir.
Bu holda D T anomal magnit maydonining miqdoriy tavsifi quyidagicha ifodalanishi mumkin:
D T = D T _L + 8 T _R + D T _f (2.3)
bu yerda D T _{l -} maydonning mahalliy komponenti; 8 T _p -komponent maydoni mage sabab - netic er; D T _f - maydonning fon komponenti.
Keling, quyidagi taxminlarni kiritaylik:
• Ploscha - di hisob-kitoblarda ishlatiladigan o'lchamlarga mos keladigan etarlicha kichik qismlar oralig'ida T _p komponenti , toymasin oyna doimiy;
• butlovchi D T _l butlovchi ko'ra ko'proq chastota xarakterga ega 8 T _p ;
• kuzatilgan sohadagi hissa komponenti 8 T _p komponent D T _L hissasidan sezilarli darajada oshadi .
Agar ushbu shartlar bajarilsa, magnit anomaliyalarning aniq fazoviy korrelyatsiyasi va kun yuzasi rel'efining xususiyatlari qayd etiladi. 8 T _p (yuqori tuzatishlar) maydonining komponentini aniqlash uchun quyida tavsiflangan moslashuvchan filtrlash algoritmi taklif etiladi.
Algoritm quyidagi operatsiyalarni o'z ichiga oladi:

1. 
   Shakllantirish namuna anomal magnit maydon qadrlaydi DT = { D m} va cho'qqilarni
   yordam ning H = {soat . _j } tomoni L boʻlgan surma oyna ichida , shu jumladan n x n nuqta ( n — foydalanuvchi tomonidan koʻrsatilgan toq son). Konvertatsiyani amalga oshirish
   

D T = Lev Nikolaevich D T .

2. 
   Kvadrat regressiyaning b ₂ , b ₁ , b ₀ koeffitsientlarini aniqlash : D T = b 2 h 2 + b 1 h + b ₀ normal tenglamalar tizimini Gauss usuli bilan yechish:
   

nn nn nn ~
b 0 n ² + b l II h j + b ₂ II h = CD T
i = 1 j = 1 i = 1 j = 1 i = 1 j = 1
nn nn nn nn
(2.4)b 0 II h jj + b 1 II h + b 2 II h j = 11 D T j h j
i = 1 j = 1 i = 1 j = 1 i = 1 j = 1 i = 1 j = 1
nn nn nn nn
b 0 II h + b 1 II h j + b 2 II h 4 = CP Th
l = 1 j = 1 i = 1 j = 1 i = 1 j = 1 i = 1 j = 1
Ushbu tizimni shakllantirishda aniqlanmagan maydon qiymati (D T j = 999) o'rgatmaydi - Quvur liniyasi.

3. 
   Magnit maydonning umumiy o'zgaruvchanligini hisoblash:
   

n n ти(Д ~ 1 i=1 j =1
Д T )2Cp )
(2.5)
3 

bu erda D T = - ¹ - CD T- oynadagi maydonning o'rtacha qiymati;

^p n i = 1 j = 1
shuningdek, regressiya tenglamasidan og'ishlar:
1 nn / ~ \ 2
(2.6)D = ----C (D T - bh ² - bh - b _n )
to'xtatish > - l * - ¹ U 2 ij 1 ij 0 /
n - 4 i = 1 j = 1
va Fisher mezoni F = D / D oyna ichida.

4. 
   Magnit maydonining bir qiymatini tahlil qilingan parametrlar sonidan chiqarib tashlash D T , buning uchun regressiyadan og'ish D _ocm maksimal qiymatni oladi. Ushbu qiymat D T ga 999 nol qiymat kodi beriladi. Run - NII da 2-o'tish nuqtasida quyidagi shartlardan biri amalga oshiriladi:
   

T d dam > s 0 (2.7)
k ₀ (2.8)
bu erda S 0 - regressiyaga bog'liqlikni qurishning talab qilinadigan aniqligi; k ₀ - rad etilgan nuqtalarning ruxsat etilgan maksimal soni.

5. 
   Ko'rsatkichlar (i, j) bilan toymasin oynaning markaziy nuqtasi uchun anomal magnit maydonni hisoblash (yuqori tuzatish ) :
   

(2.9)
^d e ^f a = bhb. + ^b ! ^h ij + ^b 0 .

6. 
   Sürgülü oynani fazoning keyingi nuqtasiga o'tkazish va 1 nuqtaga o'tish.
   
   

Shunday qilib, hisob-kitoblar natijasida nafaqat yuqori tuzatish 8 T _p , balki uni hisoblashning sifat ko'rsatkichlari ham aniqlanadi: regressiyaga bog'liqlikni qurishning ildiz o'rtacha kvadrat xatosi d / D _ost va F-. mezon. Ushbu algoritm yuqori shovqin immuniteti bilan tavsiflanadi.
TPPM dasturi yordamida namunaviy misolni ko'rib chiqaylik, tarmoqdagi 6400 ballni 200x200 m masofada vertikal - magnitlangan jinslar J = 3 A / m bo'lgan geologik qismning yuqori qismidan hisoblab chiqilgan . AMS uchun sharoit 250 m balandlikdagi relyef atrofida to'liq oqim bilan taqlid qilingan."Palitra" ning maydoni 16 kv.km (15.A-rasm). Dalada u qo'shimcha ravishda qo'shilgan tasodifiy shovqin teshiklar bo'ylab taqsimlangan - nol o'rtacha va rms bilan normal qonun Mos kelmaslik - neniem nT ± 50 (15.B-rasm).
Yuqoridagi algoritmdan foydalanib, kunduz yuzasining topografiyasi tufayli yuzaga kelgan anomal effekt taxminan aniqlandi (15B-rasm). Morfologiya natijasi - maydon topografiyasini tuzatish va kvadratik regressiya funktsiyasi asosida hisoblangan ln R T = b ₂ h ² + b ₁ h + b ₀ , va ularning statistik xarakteristikalari etarli darajada - lekin bir-biriga yaqin (11-jadval.). Amplituda ildiz-o'rtacha-kvadrat farq bt _p , shovqin huzurida tahliliy (TPPM dasturiy ta'minot) va statistik usullar bilan aniqlanadi, NT, masalan ± 20,9 hisoblanadi Bu xatolar uchun teng bo'ladi D ^T , yo'q - amaliy muammolarni hal muqarrar mavjud. Dastlabki magnit maydonda shovqin bo'lmasa, turli xil algoritmlar yordamida olingan yuqori tuzatishlarning RMSD 8 T _p , faqat ± 9,4 nT ni tashkil qiladi.
Шкала раскраски поля, нТлI I I II I I II IО 200400600

Taxminan km 2 km 4 km