Розвиток теорії та моделювання термодинамічних процесів у корі та мантії Землі

Міждисциплінарна фундаментальна науково - дослiдна робота

Шифр теми  ОБ-35Ф

 1. ПІДСТАВИ ДЛЯ ВИКОНАННЯ РОБІТ

Наказ Міністерства Освіти і Науки України № 1043 від 17.11. 2008

Наказ ректора Львівського національного університету імені Івана Франка № Н-697 від 29.12. 2008

Виконавець:

кафедра теоретичної фізики, кафедра фізики Землі Львівського національного університету імені Івана Франка

Наукові керівники:

Ткачук Володимир Михайлович, професор кафедри теоретичної фiзики Львівського національного університету імені Івана Франка, доктор фiзико-математичних наук, професор

Фурман Віталій Васильович, завідувач кафедри фізики Землі Львівського національного університету імені Івана Франка, ступінь кандидат фізико-математичних наук, доцент

Вiдповiдальний виконавець:

Хом’як Микола Миколайович, канд. фiз.-мат. наук, доцент кафедри фізики Землі Львівського національного університету імені Івана Франка

Термiн виконання:

1.01.2009 — 31.12.2011

Сучасні дані свідчать, що структура верхньої мантії континентів має визначальне значення для еволюції Землі. Дотепер всі глобальні конвективні моделі розраховувалися, виходячи з припущення, що мантійні неоднорідності густини зумовлені винятково варіаціями поля температур. Однак сейсмічні і теплові дані вказують, що температура верхньої мантії під давніми кратонами різко знижена, принаймні, на 400-600°С у порівнянні з більш молодими континентальними й океанічними структурами. Отже, повинна існувати сила, що дестабілізує ці утворення, у результаті чого континентальна частина верхньої мантії матиме тенденцію до занурення, що значно обмежить їхній час існування. Проте міра цієї компенсації дотепер не визначена, також не з'ясовано, до яких глибин можуть поширюватися теплові і хімічні аномалії під континентами. Геодинамічні процеси на Землі відбуваються в умовах складного напруженого стану, що є інтегральним показником великої кількості джерел напружень різних ієрархічних рівнів. Дослідження конвективних плинів з використанням теоретичного моделювання глибинної геодинаміки повинні визначати структури двошарової конвекції та межі турбулентного режиму в горизонтальному шарі, що підігрівається знизу. Це дасть змогу оцінити енергетичні, часові й просторові параметри теплових плюмів і гарячих точок Землі.
Проблеми математичного моделювання в рамках існуючих підходів зумовлені довготривалістю і масштабністю фізичних процесів та факторами, що впливають на їхній перебіг у реологічно складних середовищах глибинних структур кори та мантії Землі. Різні підходи до пояснення рухів та деформацій земної кори і мантії на основі механіки суцільного середовища та механіки руйнування широко застосовуються в комп’ютерному моделювання нелінійних ефектів взаємодії, з урахуванням глобальних і регіональних структурних особливостей. Публікації останніх років свідчать про необхідність розробки і застосування досить гнучкого програмного забезпечення, реалізації потужних і ефективних методів розв’язування задач математичної фізики. На сьогодні в Україні розробки в цьому напрямі є недостатніми, і опираються переважно на іноземні, вузько-орієнтовані комерційні програми, що стримує їхнє впровадження як інструменту наукових досліджень геофізичних проблем.

Хоча мантія недоступна для безпосереднього вивчення, про неї вже накопичена значна інформація, що базується на застосуванні непрямих методів дослідження. Раніше передбачалося, що структура й еволюція мантійних потоків не залежить явно від структури і властивостей континентальної літосфери. Інтенсивність підігріву в мантії Землі істотно перевищує природні можливості теплообміну, що відповідає за виникнення конвективних потоків. За відсутності внутрішнього підігріву припускають процес остигання планети ззовні, а конвективні потоки могли б виникнути через нерівномірність остигання поверхні Землі. Але теплообмін залежить від градієнта температур і під час остигання він повинен знижуватися.
Отже, якщо з погляду фізики завданням є створення фізичної моделі досліджуваного об’єкта, то з погляду проблеми геофізичних досліджень опрацювання фізичної моделі – це перший етап, за яким настають етапи формування геофізичного зображення об’єкта, а потім вирішення завдання геологічної класифікації. Всі горизонти кори від чохла до межі Мохо за певних умов поводяться тектонічно активно, зазнаючи деформації. Тектонічне розшарування знижує міцність і монолітність літосферних плит у вертикальному розрізі. Актуальним є вивчення термодинамічних умов деформування порід у різних внутрішньоплитних обставинах, створення узагальнених моделей деформування для структур різного масштабу в земній корі, реконструкція і вивчення зон динамічного впливу, удосконалення моделей формування розшарованої нижньої кори й верхньої мантії із різними. Досягнення достовірного результату можливо за умови переходу до узагальнених реологічних моделей, що визначають залежність термодинамічних процесів від глибинних структур Землі та дають змогу оцінити їхню роль у формуванні фізичних полів під час геологічної еволюції планети.
Теоретичні напрацювання, зокрема авторів проекту, в розв’язуванні фазових рівнянь теорії розсіяння та обернених задач, у побудові ефективних алгоритмів на основі методу скінченних елементів для дослідження неоднорідних структур планується використати для побудови, обґрунтування і дослідження математичних моделей в задачах фізики Землі. Даний проект є продовженням попередньо виконаних тем “Розвиток теорії та методів дослідження фізичних полів Землі” (№ держ. реєстації 0100U001448), “Розвиток методів комп’ютерного моделювання геодинамічних характеристик структурних оболонок Землі” (№ держ. реєстації 0103U001920) і “Геодинамічні, термодинамічні і механічні характеристики глибинних структур та моделювання фізичних полів Землі” (№ держ. реєстрації 0106U001324). 

Сучасні фізичні та математичні підходи до створення моделей геодинамічних процесів та їхній взаємозв’язок із термодинамічними характеристиками глибинних структур Землі є надзвичайно важливими. Незважаючи на успіхи комп’ютерного моделювання мантійної конвекції, теорія тектоніки плит ще не завершена. Перехід у фундаментальних дослідженнях Землі до узагальнених моделей середовища і процесів дасть змогу визначати залежність геофізичних і геодинамічних процесів від термодинамічних характеристик глибинних структур Землі та їхню роль у геологічній еволюції планети. Створення термодинамічних моделей на основі встановлення рівнянь руху потоків надзвичайно в’язких фрагментів речовини мантії з неоднорідним розподілом густини та нелінійною взаємодією з іншими фрагментами мантії та кори Землі є важливим фактором досліджень структурних оболонок Землі на основі аналізу фізичних полів Землі.
Побудова самоузгодженої теплової та конвективної моделі мантії Землі є незалежним методом визначення розподілу температури і теплового потоку у верхній мантії Землі. Дослідження структури розподілу густини мантії має вирішальне значення для розуміння еволюції Землі і є рушійною силою мантійної конвекції. Однак наявні дані є недостатніми для повного розуміння природи термодинамічних та конвективних процесів, що відбуваються в мантії.
Використанням математичного моделювання і сучасних комп’ютерних технологій у вивченні геологічних і геофізичних процесів, реологічних моделей, знаходження термодинамічних характеристик та їхнє узгодження з геологічною будовою через удосконалення методики обробки і аналізу геофізичних даних дасть змогу значно ефективніше реалізовувати комплексну інтерпретацію, відтворювати термодинамічні умови та можливий прояв режимів розвитку палео-геотермодинамічних систем, зокрема, на прикладі Карпатського регіону та прилеглих територій.
Результати розробки будуть спрямовані на отримання функціональних залежностей для розподілів температури і теплового потоку в надрах Землі через розв’язання оберненої задачі сейсміки на основі рівняння Вільямсона–Адамса для градієнта густини реологічно неоднорідних середовищ кори та мантії Землі.
Результатом виконання проекту також має стати розвиток теоретичних основ математичного моделювання, їхня реалізація у вигляді програмного та інформаційного забезпечення, а також апробація підходу на реальних геологічних задачах. Це дасть змогу значно досконаліше за меншими видатками оцінювати прояв глибинних термодинамічних ситуацій – актуальної проблеми геологічних досліджень.
Очікувані результати вирізнятиме наукова концепція фізико-математичного моделювання термодинамічних характеристик реологічно неоднорідних середовищ глибинних структур Землі, що, безумовно, відповідає світовому рівню нового напрямку на стику сучасної фізики та геології.

Мета роботи: аналіз перспектив створення нових фізичних моделей, які описують структури Землі та їхню еволюцію, а також методичних підходів щодо класифікації структур і вимірювання параметрів фізичних полів Землі;
теоретичні дослідження структурних властивостей реологічно неоднорідних середовищ кори та мантії з урахуванням диференціації густини в мантії, пов'язаної з варіаціями термодинамічних характеристик, що є рушійною силою мантійної конвекції;
розробка нових принципів вирішення однієї з актуальних проблем побудови математичних моделей структурних оболонок Землі та теоретичних основ ефективних алгоритмів для моделювання термодинамічних характеристик геофізичних процесів.
Для досягнення цього необхідно вирішити такі завдання:
- розробити термодинамічні моделі реологічно неоднорідних середовищ глибинних структур Землі на основі встановлення відповідних рівнянь руху потоків надзвичайно в’язких фрагментів речовини мантії та кори Землі;
-створити для опису стану глибинних структур та фізичних полів Землі моделі теорії функціонала густини і термодинамічного аналізу суцільного середовища;
- у розрахунках узагальнених термодинамічнх характеристик врахувати у геодинамічних моделях неоднорідність реологічних властивостей глибинних оболонок Землі ;
- удосконалити методики комп’ютерного опрацювання інформації для точнішого вивчення структури, густини і розподілу визначальних характеристик для реалізації досконалішої комплексної інтерпретації геофізичних даних з метою реконструкції фізичних умов та характеру змін прояву геодинамічних характеристик літосфери та вірогідних режимів розвитку на прикладі палео-геодинамічних систем геологічних структур Карпат.

ЕТАПИ РОБIТ ПО ТЕМI

Шифр етапів

Назва етапiв виконання

Термiн виконання

Результати виконання робiт

Початок

Кінець

1

2

3

4

5

1.

Дослідження прояву фізичних полів Землі у формуванні теплового режиму та фізичної природи явищ у реологічно неоднорідних середовищах Землі

01.01.2009

30.06.2009

Технічне завдання. Модель теплової конвекції з урахуванням структурних особливостей середовищ Землі.

2.

Знаходження залежностей розподілів температури і теплового потоку в надрах Землі.

01.07.2009

31.12.2009

Математична модель деяких задач теплової конвекції.

3.

Термодинамічний аналіз реологічно неоднорідних середовищ на межі кори та мантії Землі.

01.01.2010

30.06.2010

Термодинамічні моделі властивостей реологічно неоднорідних середовищ кори та мантії Землі.

4.

Розробка термодинамічних моделей структур Землі із неоднорідним розподілом густини.

01.07.2010

31.12.2010

Методика комп’ютерного моделювання структур Землі з неоднорідним розподілом густини.

5.

Удосконалення методик комп’ютерного моделювання термодинамічних характеристик геологічних структур.

01.01.2011

30.06.2011

Розв’язки типових задач моделювання термодина-мічних характеристик гео-логічних структур.

6.

Розробка методик щодо врахування термодинамічних характеристик реологічно неоднорідних середовищ на межі кори для розрахунку процесів масотеплопереносу.

01.07.2011

31.12.2011

Методика комп’ютерного моделювання неоднорідних середовищ на межі кори та мантії Землі.

Результатом розробки має стати розвиток теоретичних основ математичного моделювання, їхня реалізація у вигляді програмного та інформаційного забезпечення, а також апробація підходу на реальних геологічних задачах. Це дасть змогу значно досконаліше за меншими видатками оцінювати прояв глибинних термо-динамічних ситуацій – актуальної проблеми геологічних досліджень. Отримані результати будуть спрямовані на отримання функціональних залежностей для розподілів температури і теплового потоку в надрах Землі через розв’язання оберненої задачі сейсміки на основі рівняння Вільямсона–Адамса для градієнта густини реологічно неоднорідних середовищ кори та мантії Землі.
Використання результатів різнопланове:
- під час підготовки фахівців високої кваліфікації з геології, екології, гідрогеології та геофізики – в курсах лекцій «Основи геофізики», «Фізика Землі», «Геофізичні методи досліджень», «Геодинаміка та тектонофізика», «Термодинаміка глибинних структур Землі», «Моделювання природних структур геологічного середовища», «Геофізичний моніторинг геологічного середовища», «Геотектоніка»;
- планується написання 1 підручник, 1 навчальний посібник для студентів геологічного та геофізичного профілю ВНЗ;
- для створення конкурентноспроможних методик експресної оцінки деформаційного, геодинамічного, термодинамічного стану глибинних та поверхневих структур літосфери і мантії в зонах вірогідних родовищ корисних копалин на різних етапах їхнього вивчення й освоєння.
Впровадження результатів передбачається у ВНЗ Міносвіти і науки України, Інституту геології i геохімії горючих копалин, Інституту геофізики НАН України імені С.І. Субботіна, Карпатського відділення інституту геофізики НАН України імені С. І. Субботіна, геологорозвідувальних експедиціях системи Державної геологічної служби України, ДГП “Захiдукргеологiя”, Західно-Українській геофізичній розвідувальній експедиції, Львівському відділенні Українського державного геологорозвiдувального інституту, Державному управлінні екобезпеки Львівської області.
Форма впровадження – через використання методичних рекомендацій за результатами досліджень зацікавленим організаціям під час планування, здійснення і корегування пошукових робіт та здійснення спільних науково-дослідних робіт у геологорозвідувальних установах та підприємствах Міністерства екології i природних ресурсів України.