ISSN 222Б-3780
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ — ПРИКЛАДНЫЕ АСПЕКТЫ
УДК 519.6
А. С. Бешанова МАТЕМАТИЧНЕ МОДЕЛЮВАННЯ
ПРОСТОРОВОГО РОЗПОДІЛУ КОРИСНИХ КОПАЛИН З ВИКОРИСТАННЯМ ФОРМУЛ ШЕПАРДА І ЛИТВИНА
Тези посвячені розвитку методів математичного моделювання просторового розподілу корисних копалин з використанням даних про розподіл у похилих свердловинах за допомогою інтерполяційних формул Шепарда та Литвина.
Ключові слова: інтерполяційні функції Шепарда та Литвина.
1. Вступ
Серед методів побудови та дослідження математичних моделей розподілу корисних копалин найпоширенішими до недавнього часу були методи, що базувалися лише на мінімальній кількості інформації, яка знаходилась в кернах свердловин. Більш складними і фінансово більш затратними є методи сейсмічної томографії та методи, що використовують дані із штучних супутників Землі. Тому актуальна розробка та дослідження нових методів побудови таких математичних моделей на підставі даних з кернів свердловин і застосування сучасних методів теорії наближення функцій багатьох змінних.
2. Постановка проблеми
Аналіз методу побудови тривимірної моделі розподілу корисних копалин на підставі даних про розподіл корисних копалин у кожній точці заданої системи свердловин та методів інтерлінації функцій двох змінних. Детально досліджено формули інтер-лінації функцій трьох змінних з використанням узагальнень глобальних інтерполяційних формул Шепарда та Литвина. Наведено властивості побудованих математичних моделей та перспективи їх використання для розвідки корисних копалин.
3. Основна частина
3.1. Аналіз досліджень і публікацій. Значним є вклад в математичне моделювання об’єктів, процесів, явищ І. В. Сергієнка, В. С. Дейнеки, В. В. Скопецького, А. А. Самарського, Г І. Марчука, В. Л. Рвачова, О. М. Литвина та інших [1]. Основний вклад над цією проблемою висвітлюється у монографії Литвина О. М. математичні моделі розподілу корисних копалин, при відповідному їх узагальненні, можуть бути використані також для створення математичних моделей вмісту кори Землі за допомогою даних сейсмічного зондування [2].
3.2. Виклад основного матеріалу. Застосування наведених просторових математичних моделей у разі розвідки й розробки родовищ корисних копалин надає нові можливості для оптимізації вибору горизонтів, напрямків видобутку корисних копалин, оцінки їх запасів, дослідження меж зон вилуджування кам’яної солі тощо.
Зауважимо, що нерегулярно розміщені у просторі експериментальні дані використовуються у багатьох областях для відновлення неперервних поверхонь. У цих нерегулярно розміщених точках можуть бути задані значення різної фізичної, хімічної тощо природи: в метеорології їх отримують зі станцій спостереження погоди; в географії — встановлюючи координати розміщення об’єктів; при плануванні міст і регіонів — це центри зон збору інформації; в біології — точки спостережень тощо. При цьому припускається, що кожній заданій точці ставиться у відповідність одне число. Однією з загальних інтерполяційних формул, узагальнення яких можна використовувати для інтерлінації функцій трьох змінних у випадку нерегулярно розподілених прямих свердловин, є двовимірна глобальна інтерполяційна формула Шепарда.
У випадку, коли дані точки формують регулярну сітку, існує кілька розв’язків. До найважливіших розв’язків для прямокутної сітки належать: подання відповідної поверхні у вигляді гіперболічного параболоїда для кожних чотирьох заданих точок за допомогою білінійної інтерполяції у вказаних чотирьох точках; подання відповідної поверхні у вигляді полінома в околі п2 (п = 2,3,4,5,...) точок з використанням ньютонівських розділених різниць; подання відповідної поверхні за допомогою бікубічних сплайн-інтерполяційних формул [3]. Альтернативний метод полягає в побудові полі-номіальних або тригонометричних функцій м
Е (X У ) = Х ВА (х > у)
і=і
47------^
technology audit and production reserves — № 5/2(7), 2012, © A. Beshanova
математическое моделирование — прикладные аспекты
ISSN 222Б-37В0
з достатньою кількістю коефіцієнтів Ві, і = 1, М, щоб задовольнити всі умови
g (X;, Уі )= Zi , І = 1, N .
Як наслідок, для знаходження невідомих коефіцієнтів Ві , і = 1, М, отримуємо систему лінійних алгебраїчних рівнянь:
м ___
^ В^ (х(, у і) = Zi, і = 1, N.
І=1
Якщо базисні функції ^ (х,у), і = 1,М, є алгебраїчними або тригонометричними поліномами від двох змінних, то при великих N > М розв’язування таких систем вимагає великої кількості арифметичних операцій. У загальному випадку доцільно застосовувати не локальні формули наближення, а глобальні, які можуть передавати деякі особливості наближуваної функції на досить великі відстані від заданих точок. Такі наближення (а не локальні) можна використовувати для прогнозування розподілу корисних копалин за межами області, у якій розміщені свердловини.
Для випадку двох змінних х, у глобальну формулу Шепарда для сітки вузлів (X^,У), & = 1,М, і = 1, N, можна подати так:
3.3. Практичне значення одержаних результатів —
використання і реалізація ідей цієї роботи може підвищити науковий рівень геологорозвідувальних робіт та ефективність розробки уже розвіданих родовищ, оскільки лише тривимірні математичні моделі структури кори Землі найбільш повно можуть відображати структуру кори Землі. Наприклад, вугілля має пріоритет у паливоенергетичному комплексі України: вугілля — головний енергоносій, частка якого становить більше ніж 90 % запасів вітчизняних енергоресурсів (запаси кам’яного вугілля складають понад 43 млн. тис. т, а бурого — понад 2,5 млн. тис. т), згідно з Анциферовим А. В. Процес вугледобування належить до найбільш складних при розробці корисних копалин [1].
Література
1. Литвин О. М. Математичне моделювання розподілу корисних копалин методами інтерлінації та інтерфлетації функцій [Текст] / О. М. Литвин, Н. І. Штепа, О. О. Литвин. — Київ : Наукова думка, 2011. — 228 с.
2. Литвин О. М. Інтерлінація функцій та деякі її застосування [Текст] / О. М. Литвин. — Харків : Основа, 2002. — 544 с.
3. Литвин О. М. Інтерлінація функцій [Текст] / О. М. Литвин. — Харків : Основа, 1993. — 235 с.
Sm, n х( f; x, y)=
І £f (Xk Y )(V( x - Xk )2 + (y-Ye )2
k=0e=0
-,J( x - Xk )2 + (y - Y )2 * о V k=1, M, e=1,
м=0 N=0/1------- ---------2
ii(v( x - X )2 + (-Yj)
i=0 j=0 ^
f (,Yj), якщо З i є{1,..., M}, j є{1,...,N}: (x - Xi )2 + (y - Yj )2 = 0.
N,
Наведемо також глобальну інтерполяційну формулу Литвина, яку будемо використовувати для інтерлінації довільної функцій f (х, у, z), f єС (3), f(Х^,У^,2)=у£(2), & = 1,М. Побудуємо інтерліна-ційний оператор
м
Ом,х(І;х,у,z)=YJlkмх\(х,у), ^ 1, м = 2,3,...,
k=l
ек (х; у)
і м ,к,х (х; у) = ---—,
£ е (х; у)
І=1
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПРОСТРАНСТВЕННОГО РАСПРЕДЕЛЕНИЯ ПОЛЕЗНЫХ ИСКОПАЕМЫХ С ИСПОЛЬЗОВАНИЕМ ФОРМУЛ ШЕПАРДА И ЛИТВИНА
А. С. Бешанова
м
ek іx; y )=eM,k, хіx; y)= П
м di іx; y) i=1М k
d X м
I di іx;у)
i=1,i*k “а П di іx;y)Xk
i =1,їф k
di (x, y) =^(Xi - x)2 +(Yi - yf di,k =yj(Xi - Xk )2 + (Yi - Yk )2 ,
які для випадку уk (г) = уk = const, k = І,M, є інтерполяційними операторами на нерегулярній сітці вузлів, запропонованими Литвином в 1990 р. [1].
Тезисы посвящены развитию методов математического моделирования пространственного распределения полезных ископаемых с использованием данных о распределение в наклонных скважинах при помощи интерполяционных формул Шепарда и Литвина.
Ключевые слова: интерполяционные формулы Шепарда и Литвина.
Антонина Сергеевна Бешанова, магистрант кафедры математики в обучении Бердянского государственного педагогического университета, тел.: (067) 6098344, e-mail: [email protected].
MATHEMATICAL MODELING OF THE SPATIAL DISTRIBUTION OF MINERAL DEPOSITS USING DATA ABOUT THE FORMULS SHEPARD AND LYTVYN
A. Beshanova
The article is devoted to the development of the methods of mathematical simulation of spatial distribution of mineral resources using data distribution in inclined wells through interpolation formulas Shepard and the Lytvyn.
Keywords: globalm іnterpolyatsіynih funk^y Shepard that Lytvyn.
Antonina Beshanova, magistrant of department of mathematical
in studies Berdyanskiy state pedagogical university,
tel.: (067) 6098344, e-mail: [email protected].
С
48
технологический аудит и резервы производства — № 5/2(7), 2012