МЕТОД МОДЕЛЮВАННЯ РУХУ ТЕМПЕРАТУРНОГО ФРОНТУ ЗА НЕІЗОТЕРМІЧНОЇ ФІЛЬТРАЦІЇ

Динаміка об’єктів з розподіленими параметрами описується диференціальними рівняннями в частинних похідних параболічного типу, які з крайовими умовами є мате- матичними моделями багатьох нестаціонарних нелінійних процесів. Математичними моделями тепломасопереносу є системи рівнянь параболічного типу з такими ж гранич- ними умовами. Усі реальні процеси, як правило, є нелінійними. Вибір оптимального методу розв’я- зання тієї або іншої задачі теорії поля і технічного засобу для її реалізацій є складним питанням. У наш час найбільше поширення при математичному моделюванні складних об’єк- тів з розподіленими параметрами одержали методи дискретизації математичної моделі шляхом просторово-тимчасового квантування. Представлення математичної моделі об’єктів з розподіленими параметрами системами звичайних диференціальних або алгебраїчних рівнянь дозволяє моделювати їх на аналогових і цифрових обчислю- вальних машинах. Можна прийняти, що час роботи циркуляційної системи обмежений часом досягнення температурним фронтом експлуатаційної свердловини. Проведеними дослідженнями [1] встановлено, що теплоприток від гірського масиву, що оточує шар, у реальних пласто- вих умовах не виявляє істотного впливу на час роботи циркуляційної системи в постій- ному температурному режимі. Тому в розрахунках теплопритоком нехтуємо. У добуванні геотермальної енергії має місце напірна фільтрація, при якій величина μ має значення порядку 10-6 м-2. У зв’язку з цим система виходить на стаціонарний режим за час, малий у порівнянні з часом її роботи. У статті пропонується метод моделювання руху температурного фронту з вико- ристанням диференціальної моделі з переходом до кінцево-різницевої. Після обчислення першого наближення значення швидкості руху холодної води це значення уточнюється з використанням ітерацій за різними параметрами моделі.