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

Guruch. 16. Anormal magnit maydon bilan chiziqli yaqinlashishni yechishning iterativ jarayonining xarakteristikalari ( 4 7) _I Janubiy Xarayelax platosi: lgF 1 , IgF ₂ - logarifmlar sifat ko'rsatkichlari;
/ l - "simli tarmoq" ning nuqtalari soni;
I n - takrorlash raqami.
butun mintaqa P bo'yicha dalalarning irqiy farqi ; e ning quyidagi qiymatlari o'rnatiladi: e ₁ = 10 nT, e ₂ = 3 nT, e ₃ = 1 nT. To'g'ridan-to'g'ri keyin echimlar umumiy soniga mutanosib hisoblash vaqti N iteratsiyasini amalga oshirishda "ramka tarmog'i" yordamida - kottejlar, m x n x N x S / ( m x n x N - Xn) ~ 2,2 marta kamayadi va Intel-80486 protsessorli kompyuter uchun 2544 noma'lumli 2544 tenglamali tizimlar uchun taxminan 15 daqiqa. Koeffitsient S = N / N, N ₁ - erishish uchun zarur bo'lgan iteratsiyalar soni - SLAE Seidel an'anaviy usulini hal qilishda F1 va F2 qiymatlarini yaqinlashtirish .
Qayta tiklash maydoni A g * yoki A T * kosmosda oldindan belgilangan nuqtalarda ma'lum geometrik va fizik bug - metrga ega bo'lgan taxminiy strukturaning to'g'ridan-to'g'ri tortishish muammosini hal qilish uchun kamayadi . Transformantları hisoblash bo'lsa, tegishli operatorlar L ishlatiladi , operator o'rnini F so'z (3.10). Masalan, tortishish maydonining gorizontal gradientini hisoblashda V _xz : L = 3 ( Z - x) ( Z - z) / R ⁵ .
Magnit maydonni manbaga o'xshash yaqinlashtirish algoritmi tortishish maydonini yaqinlashtirish uchun yuqorida tavsiflangan algoritmdan tubdan farq qilmaydi. Ma'lumki, amplituda anomaliyalari bir necha ming nT dan oshmasa, A T parametri (geomagnit maydon kuchining umumiy vektorining birlik o'sishi) o'lchanadi - emy va zamonaviy kvant proton magnitometrlari asosan garmonik funktsiyadir [76].
Field A T taxminiy maydon A T * grid ekvivalent modeli m x n tashkil yupqa majmuasi, yarim cheksiz o'qi OZ (qutblar), kazh ostida joylashgan - Doi mos yozuvlar nuqtasi maydoni chuqurlikda h = 1,5 A x. Ushbu yaqinlashuvchi kesishma maydoni - (x, y, z) e P
koordinatalarida ruktsii formula bilan ifodalanadi:
mn
^D T * ⁽ x , y, z ) = XX ^J j ^T j (3.16)
где Т= a l-x- + m — + n-z- ^ R3R3R3= = 1 j = 1
- magnit maydonni aniqlaydigan garmonik funktsiyalar
Yagona yarim cheksiz vertikal yupqa novda (qutb) ning D T at
J _z = 0,01 A / m; a - kesma maydoniga qarab doimiy koeffitsient
tayoq; R = - x ) ² + ( n - y ) ² + ( Z - z ) ² - sr ning yuqori uchi orasidagi masofa ­
zhnya ( t, n, Z) va maydonni hisoblash nuqtalari; 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 ; J - vertikal rodlar magnitlanishi J _z (keyingi o'rinlarda - magnitlanish) bilan chiziqli bog'liq bo'lgan parametrlar .
Turbo Pascal 7.0 tilida yozilgan IBM - Combine - Mykh kompyuterlarida foydalanish uchun mo'ljallangan APGRV dasturidagi ReA - tortishish maydonining manbali yaqinlashuv algoritmi . Hisoblash jarayoni 2 bosqichni o'z ichiga oladi:
• ekvivalent manbalar (sferalar) to'plami bo'yicha Dd anomal tortishish maydonini yaqinlashtirish;
• maydonni qayta tiklash E d * yoki transformantlar V * har qanday berilgan sirtda S ' . SLAE yechimi, foydalanuvchining xohishiga ko'ra, to'rt xil usulda amalga oshirilishi mumkin.
iterativ usullar bilan:
• Zaydel usuli;
• dam olish usuli;
• oddiy takrorlash usuli bilan;
• Seidel usulining moslashuvchan modifikatsiyasi.
Rekonstruktsiya sohasida buzilishiga yoki chetiga bir transformantının bartaraf etish maqsadida - Ey maydoni qismlari P oldindan tashqari o'rganilgan gravity daladan - etsya doimiy fon sohasida maydoni atrofida ( "marginal zonasida" kengligi o'rtacha darajasini belgilaydi 4D x ). Fon qiymati fonconst.txt matn faylida saqlanadi.
Sferalar zichligi qiymatlarini tanlash sifati Chebishev metrikasida (3.13) nazorat qilinadi. Samarali zichligi Matrix ning sohalarda IME bilan o'zaro fayl saqlanadi - uni sig.bin. Sohalarini fazoviy tarqatish » zichlik beri, bir Geologik mazmunli yuk bermaydi noma'lumlar bo'yicha hech qanday cheklovsiz olingan tortishish qidiruvining teskari masalasining e - ekvivalent echimlaridan biridir. Infor - kerakli ma'lumotni ham sig.bin fayl bilan ta'minlaydi, har 5 yineleme, ustida yoziladi - proizoshed har qanday texnik buzilishlarida hrannost samarali axborot, - operatsiya davomida Som. Sig.bin fayli ma'lumotlarini dasturlarni qayta o'qish mumkin - agar kerak bo'lsa, oldindan tayyorlangan yechimni aniqlang .
Hisoblash jarayonining borishi displeyda ko'rsatiladi va bir vaqtning o'zida protokol.txt faylida qayd etiladi. Har bir iteratsiyada quyidagi ko'rsatkichlar paydo bo'ladi - echimlar-operatsiya va qo'shimcha parametrlar (12-jadval):
12-jadval