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

Размещено на http://www.allbest.ru/

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

математичний фізика газ трубопровід

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

Характерною особливістю таких задач є те, що в більшості випадків невідомі необхідні дані для завдання початкових і граничних умов в аналітичному вигляді. Так при моделюванні руху газу в підземних сховищах дані про температуру та розподіл тиску відомі в нерегулярних точках. Це приводить до формулювання некоректних за Тихоновим задач математичної фізики в умовах значної невизначеності. Для побудови розв'язку задач такого типу необхідно використовувати регуляризуючі алгоритми, що опираються на апріорну інформацію, яка відома в неформалізованій формі.

Крім цього в більшості випадків нестаціонарні задачі математичної фізики з постійними коефіцієнтами часто розв'язуються з використанням інтегрального перетворення Лапласа за часом. Однак тоді виникає проблема обернення, яка за своєю складністю не поступається складності вихідної задачі. Існуючі таблиці відповідності між оригіналами та зображеннями Лапласа, а також методи контурного інтегрування не завжди дають можливість знайти аналітичний розв'язок і тому необхідно будувати адаптивні числові методи обернення, які враховують апріорну інформацію про поведінку цих розв'язків.

Описання процесу переносу домішкової речовини в об'єктах природного середовища (пористих матеріалах) зводиться, як правило, до розв'язування системи взаємозв'язаних лінійних параболічних рівнянь другого порядку для концентрацій з конвективною складовою і джерелами, що пропорційні до шуканих функцій-концентрацій. Після приведення до безрозмірної форми ці системи містять малий параметр біля старшої похідної. У зв'язку з цим також виникають труднощі при застосуванні інтегральних перетворень та теорії збурень. При цьому наявність малого параметра робить малопридатними традиційні числові методи.

Зауважимо, що аналогічні проблеми виникають при розв'язуванні нелінійних нестаціонарних задач фільтрації газу в пористих середовищах і задач газової динаміки при русі газу в трубопроводах.

Одним із перспективних підходів розв'язування сформульованих проблем є застосування операційних та спектральних методів. В обчислювальних експериментах ці методи мають ряд переваг, які пов'язані з параметричним (символьним) розв'язуванням задач на різних їх етапах. Такі розв'язки дають можливість дослідити вплив параметрів фізичного процесу на його перебіг і вивчити стійкість розв'язку відносно малих збурень вхідних даних. В тих випадках, коли символьний розв'язок задачі не вдається отримати, цінність цих методів частково зберігається, оскільки можна побудувати розв'язок за зображеннями (в термінах відповідних операторів) в базисі так званих конструктивних елементів, а значення зображень тоді обчислюються з наперед заданою точністю.

Під конструктивними елементами розуміють функції, які з точністю до нормуючого множника задаються в деякому функціональному базисі скінченним набором цілих чи раціональних чисел. До такого класу елементів належать, наприклад, класичні ортогональні многочлени, кожен з яких має цілі коефіцієнти. Використання базису конструктивних елементів, разом з можливістю знаходження значень зображень, в свою чергу дозволяє організувати обчислювальний процес розв'язку задачі з гарантованою точністю. При цьому під обчисленням з гарантованою точністю розуміється така організація обчислювального процесу, при якій на всіх проміжних етапах виключаються похибки машинної (апаратної) реалізації арифметичних операцій такого порядку, що забезпечує точність кінцевого результату. За Д.К.Фадєєвим такий числовий метод розв'язку задачі називається регулярним.

При числовому знаходженні розв'язку задачі з гарантованою точністю обчислень головним чином залишаються тільки похибки методу. В тих випадках, коли не вдається організувати обчислення з гарантованою точністю, виникають додаткові ускладнення, пов'язані з нагромадженням машинних похибок. Якщо при цьому числові методи забезпечують достатню точність відновлення розв'язку задачі, то можна вважати, що зберігаються переваги символьного розв'язку шляхом застосування операційних методів в обчислювальному експерименті. Дійсно, в цьому випадку розв'язок задачі в зображеннях може вважатися точним аналітичним розв'язком при фіксованих вхідних даних.

Таким чином, у даний час актуальною є загальна проблема розробки ефективних, простих в реалізації й достатньо універсальних аналітико-числових моделей фізичних процесів та адаптивних методів розв'язування задач математичної фізики і обробки експериментальних даних, які орієнтовані на використання апріорної інформації.

В даній роботі розроблено аналітико-числові моделі руху газу в трубопроводах та масопереносу в природних пористих середовищах стосовно практичних проблем потокорозподілу та оптимізації процесів згідно до критеріїв споживання газу. Побудовано адаптивні алгоритми та методи для розв'язування типових нестаціонарних задач гетеродифузії у середовищах з мікроструктурою, фільтрації газу в пористих природних об'єктах, газової динаміки стосовно транспортування газу в трубопроводах. При цьому для конкретних прикладних задач розв'язано проблеми обробки експериментальних даних (апроксимації та фільтрації сигналів, стиску інформації, ідентифікації параметрів об'єктів), на основі спектральних представлень функцій в базисах ортогональних многочленів. Запропоновані аналітико-числові моделі згаданих процесів та схеми знаходження розв'язків дозволяють контролювати точність обчислень, виключати наростання машинної похибки, враховувати апріорну інформацію про поведінку шуканих полів в певних областях та використовувати отримані результати для керування процесів масопереносу.

Зв'язок роботи з науковими програмами, планами, темами. Дослідження за темою дисертації проводились відповідно до планових науково-дослідних робіт Інституту прикладних проблем механіки i математики iм. Я.С.Пiдстри-гача НАН України та його Центру математичного моделювання, в яких автор виступав як відповідальний виконавець та виконавець.

«Розробити математичне і програмне забезпечення типової регіональної системи обробки аерокосмічної інформації в інтересах народного господарства, виробити рекомендації відносно її створення і дальнішого введення в експлуатацію» (1984-1987, № держреєстрації 01.84. 00855510).

«Розробити математичні моделі, методи розв'язування прямих і обернених задач дистанційного зондування стосовно тіл простої конфігурації»(1988-1991, № держреєстрації 01.8.80054489).

«Розробка і дослідження математичних моделей механоелектродифу-зії в локально неоднорідних багатокомпонентних середовищах»( 1991-1995, Державний фонд фундаментальних досліджень, № держреєстрації 1/898).

«Розробити математичні моделі процесів переносу в багатофазних тілах, розвинути методи розв'язання крайових задач при неповних, достатніх або надлишкових даних» (1992-1995, № держреєстрації 0193U041575).

«Математичне забезпечення обробки даних дистанційної спектрометрії земних утворів для вирішення природоресурсних та екологічних задач» (1993, проект № 9.3/129 державної науково-технічної програми України).

«Розробка математичних моделей і програм розрахунку міграції забруднень у геологічних структурах з суттєво неоднорідними відкладами (Чорнобильська програма)» (1996-2000, № держреєстрації 0193U041577).

«Розробка та дослідження нелінійних математичних моделей механотермодифузійних процесів у багатофазних системах при локальних змінах стану компонент» (1997-2000, № держреєстрації 0199U000133).

«Розробити програмно-технічний комплекс для розв'язування задач економного і надійного газопостачання (НТП Міннауки)» (1997-2000,

№ держреєстрації 0199U000624).

«Математичні моделі та методи термодинамічного опису локально неоднорідних гетерогенних багатокомпонентних систем та об'єктів природного середовища» (1999-2003, № держреєстрації 0199U000626).

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

Для досягнення цієї мети у дисертації були поставлені та вирішені такі основні завдання:

сформулювати задачі математичної фізики стосовно руху газу в трубопроводах та вивчити вплив гідродинамічних параметрів газу та геометричних параметрів трубопроводів на цей процес;

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

сформулювати відповідні крайові задачі та дослідити процес перерозподілу домішкових речовин в приповерхневих шарах ґрунту з врахуванням їх конвективного руху;

побудувати регуляризуючі адаптивні алгоритми та методи розв'язування задач теорії масопереносу з використанням спектральних представлень в базисах ортогональних многочленiв;

з метою адаптації моделей, що вивчаються, до опису конкретних науково-технічних задач розробити алгоритми побудови початкових і граничних умов в параметричній формі;

проаналізувати основні проблеми числової реалізації розв'язку сформульованих задач математичної фізики.

Об'єкт дослідження. Процеси поширення газу в трубопроводах, фільтрації газу в пористих середовищах (підземних сховищах газу) та поширення домішкових речовин в приповерхневих шарах ґрунту.

Предмет дослідження. Математичні моделі процесів поширення газу в трубопроводах, фільтрації газу в пористих середовищах, поширення домішкових речовин в приповерхневих шарах ґрунту та аналітико_числові методи розв'язування задач математичної фізики, цифрової і статистичної обробки інформації з метою формування початкових і граничних умов для задач потокорозподілу та оптимізації процесів масопереносу.

Методи дослідження. Процеси масопереносу описуються, як правило, нелінійними диференціальними рівняннями (системами диференціальних рівнянь) в частинних похідних. Методи розв'язування таких рівнянь вивчені недостатньо. Тому для доведення моделей до числа використано аналітичні, числово-аналітичні (ітераційні) та числові методи розв'язування задач математичної фізики. В якості аналітичних методів застосовано, в основному, спектральні методи розв'язування задач в базисах класичних ортогональних многочленів з використанням теорії апроксимації функцій та асимптотичних розкладів. Шуканий розв'язок подано у вигляді ряду з невідомими коефіцієнтами, для яких на основі вхідної інформації та математичних залежностей між параметрами фізичного процесу, що моделюється, побудовано алгоритми обчислення. При цьому використано методи теорії функцій комплексної змінної, асимптотичних розкладів, функціонального аналізу, тощо.

Наукова новизна результатів роботи полягає в наступному.

Розв'язано важливу науково-технічну проблему в галузі математичного моделювання та обчислювальних методів - розроблення адаптивних аналітико-числових моделей масопереносу, які відповідають практичним проблемам оптимізації масопереносу за критеріями раціонального споживання, а також алгоритмів розв'язування відповідних крайових задач, побудованих на основі спектральних представлень функцій в базисах ортогональних многочленів, які орієнтовані на використання апріорної інформації про апроксимуючу функцію або її зображення Лапласа.

В процесі розв'язання вказаної проблеми отримано такі нові результати.

· В умовах значної невизначеності вхідних даних побудовано адаптивні моделі руху газу в трубопроводах та розроблено числово-аналітичні методи розв'язування задач газової динаміки для визначення розподілу тиску і масової швидкості в газопроводах з метою використання цих моделей для розрахунку параметрів режимів роботи газотранспортних систем.

· Сформульовано уточнені математичні моделі фільтрації газу в підземних сховищах, які враховують неканонічність геометрії пласту, наявні в ньому структурні порушення, задання вхідної інформації в нееквідистантних точках всередині пласту, а також досліджено вплив геометричних та гідродинамічних параметрів на розподіл тиску в них. Отримані результати дають можливість оцінювати колекторські властивості пластів підземних сховищах газу та керувати процесом відбирання-закачування газу.

· Досліджено шляхи та швидкість дифузії домішкової речовини в приповерхневих шарах ґрунту на основі отриманих аналітичних розв'язків відповідних задач математичної фізики при врахуванні конвективної складової та проведено оцінку коефіцієнтів процесу дифузії.

· Для визначення параметрів адаптації побудованих математичних моделей до реальних фізичних процесів запропоновано способи розв'язування обернених коефіцієнтних задач масопереносу та газової динаміки з використанням інтегрального перетворення Лапласа, асимптотичних розкладів та обробки вимірюваних даних.

· З метою параметричної побудови початково-крайових умов спектральними методами розроблено способи обчислення коефіцієнтів ортогональних розкладів для довільних їх порядків та проведено аналіз похибки відновлення вхідної функції. Виведено рекурентні формули для обчислення коефіцієнтів ортогональних розкладів, які можна використовувати як в прямому (за зростанням номерів) так і в оберненому порядку (за спаданням номерів).

· Розроблено єдиний підхід до обернення інтегрального перетворення Лапласа, побудови регуляризуючих алгоритмів для числово-аналітичного розв'язування крайових задач теорії масопереносу, а також інтегральних рівнянь типу згортки (Вольтеррівської та півбезмежної).

Запропоновані в роботі методи розв'язування прямих і обернених задач математичної фізики мають самостійне значення і є перспективними при застосуванні їх до дослідження фізичних процесів, які описуються лінійними або нелінійними системами диференціальних рівнянь в частинних похідних зі змінними коефіцієнтами в умовах невизначеності та інтегральних рівнянь типу згортки.

Практичне значення одержаних результатів.

В роботі розроблено математичні моделі процесів масопереносу та методи їх дослідження і побудови крайових умов, що дало можливість побудувати оптимізуючі алгоритми розрахунку режимних параметрів газотранспортних потоків, які приводять до зменшення енергетичних затрат в процесі транспортування газу, а також сформулювати рекомендації щодо керування основними гідродинамічними параметрами, що впливають на потокорозподіл в газотранспорних мережах. Результати теоретичних досліджень використані при розв'язанні наступних практичних задач:

Визначення розподілу тиску в газопроводах при врахуванні зміни гідродинамічних параметрів газу з метою побудови алгоритмів розрахунку газотранспортних мереж та оптимізації газових потоків (м. Київ, «Укртрансгаз»).

Визначення акумулюючої здатності газоносних пластів підземних сховищ газу з метою розрахунку режимів роботи системи пласт підземного сховища газу-газозбірний пункт (Дашавське, Червонопартизанське, Богородчанське та Опарське підземні сховища газу - м. Київ, «Укртрансгаз»).

Для розрахунку міграції забруднень у геологічних структурах з суттєво неоднорідними відкладами отримано аналітичні розв'язки задач конвективної гетеродифузії в приповерхневих шарах ґрунту на основі яких знайдені потоки агресивних домішок у ґрунтові води (Чорнобильська програма).

Розроблено пакет прикладних програм для визначення вмісту важких металів та нітратів в рослинних покровах Землі на основі обробки їх спектральних портретів (м. Київ, Центр аерокосмічних досліджень Землі Інституту геології НАН України, Національне космічне агентство України).

Побудовано математичне забезпечення та пакет прикладних програм для ідентифікації об'єктів в ЕПР-дозиметрії (м. Київ, Інститут фізики напівпровідників НАН України).

Побудовано математичне забезпечення та пакет прикладних програм для фільтрації мультиплікативних шумів в растровій інформації для ідентифікації об'єктів і розв'язування лідарних рівнянь з метою виявлення наявності та оцінки концентрації твердотільних фракцій викидів цементних заводів (м. Москва, НВО «Астрофізика»).

До дисертаційної роботи додаються акти про використання результатів прикладних досліджень для розробки програмних комплексів розрахунку режимних параметрів оптимальної роботи газотранспортної системи України та керування потоками газу з метою раціонального його використання, а також розрахунку режимів роботи системи пласт підземного сховища газу-газозбірний пункт з метою визначення акумулюючої здатності підземних сховищ газу і можливістю їх роботи в пікових режимах.

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

Апробація роботи. Основні результати досліджень доповідалися на Всесоюзній науково-технічній конференції "Автоматизація досліджень проектування й випробувань складних технічних систем" (Калуга, 1989), Регіональних науково-технічних конференціях "Моделювання й автоматизація проектування складних технічних систем" (Калуга, 1990) і "Автоматизація досліджень, проектування та випробування складних технічних систем i проблеми математичного моделювання" (Калуга, 1991), Першій Всеукраїнській конференції "Обробка сигналів i зображень та розпізнавання образів" (Київ, 1992), Всеукраїнській науковій конференції "Нові підходи до розв'язування диференціальних рівнянь" (Дрогобич, 1994), IV Міжнародній конференції з механіки неоднорідних структур (Тернопіль, 1995), Другій національній науковій конференції "Інформатика: теорія, технологія, техніка" (Одеса, 1995), Всеукраїнській науковій конференції "Розробка та застосування математичних методів в наукових дослідженнях", присвяченій 70-рiччю від дня народження професора П.К. Казiмiрського (Львів, 1995), Міжнародному колоквіумі EUROMECH Colloquium 356 "Transform methods in solid mechanics" (Munchen, 1996), Міжнародній науковій конференції «Обчислювальна математика і математичні проблеми механіки» (Дрогобич, 2001), Другій Європейській конференції з обчислювальної механіки «Solids, Structures and Coupled Problems in Engineering» (Краків, Польща, 2001), Українсько-Польському симпозіумі «Інфороматично-математичне моделювання складних систем» (Львів-Брюховичі, 2002), 1-й Міжнародній науково-технічній конференції "Развитие компьютерных комплексов моделирования, оптимизации режимов работы систем газообеспечения и их роль в диспетчерском управлении технологическими процессами в газовой отрасли" (DISCOM -2002), (Москва, 2002), 4-й Українсько-Польській конференції «Механіка середовища, методи комп'ютерних наук та моделювання» (Львів, 2004), Міжнародній конференції „International workshop on free boundary flows and related problems of analysis», (Київ, 2005).

У повному обсязі дисертаційна робота доповідалась і обговорювалась на семінарах Центру математичного моделювання Інституту прикладних проблем механіки i математики iм. Я.С.Пiдстригача НАН України, Інституту прикладних проблем механіки i математики iм. Я.С.Пiдстригача НАН України, Інституту гідродинаміки НАН України, Міжнародного математичного центру ім. С. Банаха, кафедри прикладної математики Львівського національного університету імені Івана Франка, кафедр математичної фізики та інформаційних технологій Київського національного університету імені Тараса Шевченка, Інституту кібернетики ім. В.М.Глушко-ва НАН України, Національного університету „Львівська політехніка».

Публікації та особистий внесок здобувача. Основні результати досліджень, відображених у дисертації, опубліковані у 55 наукових роботах, в тому числі у 33 статтях в наукових фахових журналах і збірниках наукових праць; 12 наукових статтях та препринтах, у 10 матеріалах міжнародних семінарів, конференцій, конгресів і симпозіумів. Всього за темою роботи опубліковано 83 праці.

Результати дисертаційної роботи автором отримано самостійно.

В спільних публікаціях автору належить: розробка методів розв'язування задач математичної фізики _ [2,9,32]; формулювання задач математичної фізики _ [17,22,46]; отримання ключових співвідношень _ [5,19,22]; проведення порівняльного аналізу методів розв'язування задач математичної фізики _ [14,28,31]; отримання числових результатів _ [20,24,38]; побудова методики визначення впливу параметрів газу на процес його руху в трубопроводах і пористих середовищах _ [23,27,38]; розробка аналітико_числових моделей процесів масопереносу _ [25,26,29]; побудова спектральних методів розв'язування диференціальних рівнянь _ [16,44]; побудова спектральних методів розв'язування інтегральних рівнянь _ [18]; проведення аналізу впливу процесу дискретизації на точність розв'язування диференціальних рівнянь в частинних похідних _ [15]; ідея використання швидкого перетворення Фур'є при сумуванні ортогональних рядів в базисі многочленів Якобі та виведення основних співвідношень _ [12]; аналіз отриманих результатів та формулювання висновків _ [21,24,28]; порівняльний аналіз моделей _ [30,37,45]; побудова спектральних методів розв'язування обернених задач математичної фізики та задач обробки інформації _ [39,43,49,53]; обернення інтегралу Лапласа числово_аналітичним методом _ [33,34,54,55]; ітераційна схема побудови початково_граничних умов, залеж-них від розв'язку задачі _ [47,50]; побудова алгоритмів програмних комплек-сів при моделюванні процесів масопереносу _ [51,52].

Структура роботи. Дисертаційна робота складається з вступу, восьми розділів, висновків та списку використаних джерел. Загальний обсяг дисерта-ції становить 290 сторінок машинописного тексту, з них 264 сторінки основ-ного тексту, 17 сторінок літератури (211 найменувань), 53 рисунків і 13 таблиць.

У вступі подано загальну характеристику роботи, розкрито суть і стан вивчення наукової проблеми, обґрунтовано актуальність теми дисертації, сформульовано її мету, відзначено новизну отриманих результатів, висвітлено практичну цінність отриманих результатів, наведено дані про апробацію результатів і публікації, що відображають основний зміст роботи, визначено особистий внесок здобувача у публікаціях, підготовлених до друку за участю співавторів.

У вступі подано загальну характеристику роботи, розкрито суть і стан вивчення наукової проблеми, обґрунтовано актуальність теми дисертації, сформульовано її мету, відзначено новизну отриманих результатів, висвітлено практичну цінність отриманих результатів, наведено дані про апробацію результатів і публікації, що відображають основний зміст роботи, визначено особистий внесок здобувача у публікаціях, підготовлених до друку за участю співавторів.

В першому розділі наведено огляд літератури за темою дисертації, окреслено її місце у вирішенні науково-прикладних проблем спектрального аналізу, операційного числення та методів розв'язування задач математичної фізики, зокрема задач газової динаміки в газотранспортних мережах, дифузії та обробки інформації.

Ряд фізичних процесів описуються нелінійними диференціальними рівняннями або системами нелінійних диференціальних рівнянь (звичайних або в частинних похідних). Такі рівняння, як правило, отримують у рамках дискретних чи континуальних моделей макроскопічної фізики, зокрема механіки суцільного середовища. Достатньо повний аналіз опису різних фізичних процесів приведено в роботі І.М. Федоткіна. Так рух газу в трубопроводі, що знаходиться в ґрунті, при відповідних початково-крайових умовах в нестаціонарному неізотермічному режимі описується взаємозв'яза-ною системою нелінійних диференціальних рівнянь в частинних похідних

Неусталена фільтрація реального газу в неоднорідному пористому середовищі (пласті підземного сховища газу) зводиться до задачі інтегрування рівняння.

Фізико-математичні моделі руху газу в трубопроводах та пористих середовищах сформульовані в роботах І.А. Чарного, А.В. Александрова, Є.І. Яковлєва, С.А. Бобровского, М.Г. Сухарєва та багатьох інших. Стосовно методів розв'язування рівнянь газової динаміки доцільно відзначити роботи М.А. Жідкової, Є.І. Яковлєва, А.В. Алєксандрова, Ф.Г. Темпеля, М.Г. Сухарєва та інших.

Зауважимо, що в деяких випадках пористі середовища, зокрема підземні сховища газу, мають неканонічну форму і початкові та граничні умови відомі в дискретних нееквідистантних точках з невеликою точністю (як правило, три-чотири значні цифри). Для прикладу на рис. 1 подано структурну карту Дашавського підземного сховища газу, де крапками позначено замірні та робочі свердловини, а жирними лініями - структурні та тектонічні зсуви пласту.

Аналітичний розв'язок конкретних задач математичної фізики, поставлених на цій основі, можна отримати тільки в окремих часткових випадках. Тому при побудові числових і числово-аналітичних методів їх розв'язування необхідно використовувати певні наближення, які повинні бути узгоджені з вибраними модельними положеннями.

Конвективна гетеродифузія домішкових речовин у середовищах з неоднорідною мікроструктурою, зокрема, пористих тілах, описується системою рівнянь

Основні положення побудови та дослідження теоретичних моделей механіки суцільного середовища розглянуто в працях А.Ерінгена, О.А.Ільюшина, Л.І.Сєдова, К.Трусдела та інших. Зокрема, моделі термомеханіки розвинуто в роботах Я.С.Підстригача, Г.С.Кіта та їхніх учнів. Значний внесок у теорію взаємозв'язаних фізико-механічних полів внесли роботи Я.Й.Бурака, Г.І.Баренблата, В.Г.Карнаухова, І.Т.Селєзова, А.Ф.Улітка та інших. Континуальна модель дифузії домішок двома шляхами з врахуванням переходів частинок розглядалась в роботах Е.Айфантіса, Я.Й.Бурака, Є.Я.Чаплі та інших.

Загальні аспекти математичних методів наближеного розв'язування таких задач розглянуто в роботах В.В.Іванова, Л.В.Канторовича, В.І.Крилова, М.М.Лаврентьєва, Г.І.Марчука, В.І.Митропольського, А.А.Самарского, О.С.Лимарченка, Ф.Г.Темпеля, І.А.Чарного, А.В.Ликова та інших.

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

Так, зокрема, описання процесу дифузії у середовищах із складною внутрішньою структурою приводить до системи диференціальних рівнянь з малим параметром та конвективною складовою, побудова розв'язку яких та їх числове дослідження існуючими методами пов'язане із значними труднощами. Крім цього, початково-крайові задачі масопереносу, особливо задачі газової динаміки і фільтрації газу в пористих середовищах, як правило, забезпечені початковими і крайовими даними в нерегулярних точках з невисокою точністю. Це приводить до потреби розробки таких схем їх розв'язування, які узгоджуються з існуючими експериментальними даними. При розв'язуванні температурних задач теорії пружності на основі перетворення Лапласа приходимо до зображень, обернення яких неможливе в аналітичному вигляді. Основою для розв'язування відповідних задач математичної фізики є спектральне подання функцій, що входять в математичну модель, та початкових і крайових умов у відповідних ортогональних базисах.

Суть спектральних методів розв'язування задач полягає в представленні відомих функцій та шуканих розв'язків ортогональними рядами в деяких базисах та побудові алгоритмів для обчислення коефіцієнтів цих рядів (узагальнених спектрів).

Спектральні методи розв'язування задач розвинуті в роботах С.Качмажа, Г.Штейгауза, Г.Сеге, П.К.Суєтіна й інших. Застосування спектральних методів в базисах тригонометричних функцій та базисах Хаара і Уолша наведені в роботах Кулі-Тюкі, Б.Голд, Ч.Рейдер, В.В.Солодовнікова, М.Д.Єгупова, В.К.Задираки та їх учнів.

Порівняльний аналіз спектральних методів в базисах Фур'є, Хаара та Уолша показав, що часто задовольняються не всі критерії, які ставляться до розв'язків сформульованих задач. Тому в літературі розвиваються спектральні методи в інших ортогональних базисах - Якобі, Чебишева-Лагерра, Ерміта тощо.

Теоретичні основи операційних методів закладено в працях В.А.Діткіна, А.П.Прудникова, Я.Мікусінського та їхніх учнів. При їх застосуванні часто виникає проблема обернення інтегралу Лапласа. Дослідженню цієї проблеми присвячені роботи В.А.Діткіна, А.П.Пруднікова, В.М.Амербаєва, В.І.Крилова, Н.С.Скоблі, Дж.Уіддера, С.А.Баймультдінової, В.М.Рябова, Дж.Влоха, Р.Беллмана, Р.Піссенса й інших.

Використання сучасних підходів до моделювання складних систем, сучасних технічних засобів збору, обробки та передачі даних, систем спряження їх із виконавськими пристроями та механізмами дає можливість в оперативному режимі приймати і реалізовувати обґрунтовані управлінські рішення. Для досягнення цього необхідно поставити і розв'язати комплекс оптимізаційних задач. Вказані проблеми розглядаються в роботах І.В.Бейка, Б.Н.Бублика, П.Н.Зінька, В.С.Дейнеки, І.В.Сергієнка., В.В.Скопецького, А.Г.Євдокимова, А.Д.Тевятева, В.В.Дубровського, М.Н.Кулика, Б.Л.Кучика, Ю.І.Максимова, Р.М.Клемфуса.

Роботи Б.В.Гери присвячені побудові математичних моделей оптимізаційного відтворення фізичних полів природних об'єктів за неповних даних та визначення теплофізичних характеристик і параметрів утворень земної поверхні. При цьому враховуються основні фізичні властивості полів земної поверхні та їх залежність від неоднорідностей і зміни характеристик середовища.

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

У другому розділі розроблено способи обчислення узагальнених спектрів як на основі аналітичного задання функцій, так і при відомих її дискретних значеннях. При моделюванні процесів масопереносу результати розділу використовуються для побудови параметричних зображень початково-граничних умов у вигляді ортогональних рядів, а в цифровій обробці інформації - для розв'язування задач апроксимації сигналів.

Формула (6) дає можливість обчислити коефіцієнти ряду (5) у випадку, коли функція задана аналітично. При відомих дискретних значеннях функції на деякій множині точок для обчислення елементів узагальнених спектрів використовуються квадратурні формули. Зокрема, якщо використати квадратурну формулу Гауса-Якобі, в якій множина є коренями многочлена порядку, то для обчислення має місце формула

Формула (7), як і інші квадратурні формули, дає можливість обчислити обмежену кількість перших коефіцієнтів. В загальному випадку, як правило, їх недостатньо для відновлення функції Тоді доцільним є використання асимптотичних розкладів для при великих. В цьому розділі для певних класів функцій на базі апріорної інформації про гладкість функцій й розриви першого роду, поведінку в околах граничних точок, поведінки зображення Лапласа на безмежності та в околах крайніх правих особливих точок отримано нові асимптотичні формули для узагальнених спектрів при великих, які з достатньою для практики точністю дають можливість обчислювати значення.

Аналіз існуючих розв'язків задач математичної фізики, зокрема задач газової динаміки та механіки деформівного твердого тіла, показує, що в околах кінців області визначення вони мають степеневу, показникову, синусоїдальну та логарифмічну поведінки за часом. У внутрішніх точках області визначення розв'язки можуть мати розриви першого роду. Поряд з цим, при дослідженні фізичних процесів в нестаціонарному випадку, що описуються диференціальними рівняннями або системами диференціальних рівнянь в частинних похідних зі сталими коефіцієнтами, використовується інтегральне перетворення Лапласа. Тому доцільно будувати асимптотичні розклади для узагальнених спектрів виходячи із знайдених зображень Лапласа шуканих розв'язків.

Якщо функція розкладається в ортогональний ряд за многочленами Чебишева-Лагерра

Для многочленів Чебишева-Лагерра має місце інтегральне представлення

Оскільки функція має мінімум в точці то, як слідує з теорії асимптотичних розкладів інтегралів цього типу, їх асимптотика визначається поведінкою при. За теоремою Ефроса перетворення Лапласа функції має вигляд

Отримано асимптотичні формули для узагальнених спектрів на основі відомого зображення Лапласа апроксимуючої функції. При використанні базису Якобі знайдено подання узагальнених спектрів Фур'є-Якобі у вигляді інтегралу типу Мелліна

Показано, що асимптотика узагальнених спектрів залежить від поведінки зображення Лапласа в околі безмежності та поведінки в околах крайніх правих особливих точок зображення. Формули (12) та (13) є вихідними для виведення асимптотик коефіцієнтів Фур'є-Якобі при.

Таким чином в розділі подані як відомі квадратурні формули, які дозволяють обчислювати узагальнені спектри функцій при невеликих порядках і є оптимальними в класі, так і побудовані нові асимптотики узагальнених спектрів, які дають можливість обчислювати їх для великих порядків. Отримані тут результати дозволяють знаходити узагальнені спектри для довільних порядків, тим самим будувати в задачах масепереносу початково-граничні умови в параметричній формі. Основні результати опубліковані в роботах [10,13,42].

В третьому розділі досліджено похибку апроксимації функції ортогональними рядами. Отримані результати використовуються для побудови регуляризуючих алгоритмів при розв'язуванні прикладних задач.

Недоліком мажорантних оцінок залишків ортогональних рядів є те, що вони можуть бути завищеними і неефективними у практичному використанні. Крім цього, вони отримані в припущенні абсолютної збіжності відповідних рядів, що звужує область їх застосування (в прикладних задачах абсолютна збіжність не завжди гарантується ). Тоді важливого значення набуває виведення асимптотичних формул для залишків ортогональних рядів, які поблизу граничних значень відповідних параметрів дають результати близькі до реальних. Тому для похибок типу виведені асимптотичні формули при великих значеннях порядку залишкового члена.

Для многочленів Чебишева-Лагерра має місце асимптотична формула типу Хільба.

На основі координат точок розриву та стрибків функцій, що задаються ортогональними рядами Фур'є-Лагерра, побудовано допоміжні функції, які їх враховують. Це дає можливість із загального ряду виділити його повільно збіжну частину та просумувати її. Ряд, що залишається, збігається значно швидше.

Запропоновано спосіб використання швидкого перетворення Фур'є для сумування рядів Фур'є-Якобі. Використання швидкого перетворення Фур'є дозволяє значно зменшити час проведення обчислень та кількість арифметичних операцій, а цим самим і зменшити накопичення машинної похибки при проведенні обчислень. Проведено обчислювальний експеримент на модельних задачах.

Виведені в цьому розділі мажорантні та асимптотичні оцінки залишкових членів ортогональних рядів разом з результатами попереднього розділу дали можливість побудувати регуляризуючі алгоритми апроксимації функцій ортогональними рядами. Отримано ряд нових результатів щодо їх сумування для неперервних функцій та функцій, що мають розриви першого роду. Запропоновано спосіб використання швидкого перетворення Фур'є при сумуванні ортогональних рядів Фур'є-Якобі. Результати розділу використовуються для апроксимації початкового розподілу тиску в пористих середовищах та розв'язування задач математичної фізики, що стосуються моделювання масопереносу. [12,19,20,39,43,48].

Четвертий розділ присвячений побудові адаптивної схеми відновлення оригіналу Лапласа за відомим його зображенням на базі ортогональних рядів. Інтегральне перетворення Лапласа за часом використовується при розв'язуванні диференціальних рівнянь в частинних похідних, якими описуються процеси масопереносу. Однак не завжди можна перейти від зображення до оригіналу аналітичним способом. Побудована схема дає можливість відновити оригінали шуканих розв'язків задач математичної фізики.

Суть побудованої схеми полягає в наступному. Числовому оберненню інтегралу Лапласа передують дослідження зображення для формування апріорної інформації про оригінал. Для цього функція розкладається в асимптотичні ряди в околах нескінченно віддаленої точки та крайніх правих особливих точок. На основі цих розкладів визначається поведінка оригіналу в нулі і на нескінченості (до отриманих розкладів зображень застосовуються теореми тауберового типу) та розриви першого роду (координати точок розривів та стрибки оригіналу знаходяться аналітичним чи числовим методом).

_ розрахункові коефіцієнти;- коефіцієнти неперервної частини оригіналу, обчислені за асимптотичними формулами;- коефіцієнти, обчислені за класичною схемою при точних значних цифрах в значеннях зображення;- коефіцієнти функції;

_ повна похибка методу, що складається з залишку ортогонального ряду неперервної частини оригіналу, залишку ряду, що виникає за рахунок використання асимптотичних формул для та похибки заокруглення (машинної похибки). Тут - ортогональні многочлени,- нормуючий множник, а - рубіж поганої обумовленості при значних цифрах зображення.

При побудові асимптотичних розкладів для, як правило, достатньо просто отримати перші їхні члени. Однак при використанні асимптотики для обчислення при великих виникає необхідність у знаходженні подальших членів асимптотичного розкладу. Тому в даному розділі розгляда-ється спосіб уточнення асимптотичних розкладів узагальнених спектрів методом невизначених коефіцієнтів.

На основі спектрального представлення оригіналу в базисах класичних ортогональних многочленів побудовано алгоритм його відновлення. Обчислення оригіналу зводиться до обчислення коефіцієнтів відповідних ортогональних рядів на основі відомого зображення Лапласа та відповідної апріорної інформації. Показано ефективність побудованого алгоритму обернення перетворення Лапласа для розв'язування прикладних задач. Відзначено, що використання асимптотичних формул для відновлення оригіналу та виділення повільно збіжної його частини дають можливість значно уточнити результати обчислень та зменшити час рахунку [7,8, 17,21,22,33,34,35,54,55].

П'ятий розділ дисертації присвячений побудові спектральних методів розв'язування інтегральних та диференціальних рівнянь (як звичайних, так і в частинних похідних) в базисах многочленів Якобі та Чебишева-Лагерра. При дослідженні процесів масопереносу результати цього розділу застосовують-ся для знаходження початкового розподілу тиску вздовж трубопроводу та задач обробки інформації (фільтрація мультиплікативних шумів).

Запропоновано модифікацію інтегрального перетворення Чебишева-Лагерра функції.

Тут - інтеграли від многочленів Якобі кратності.

Оскільки зображення відомі, а при зображення, то можна вважати, що, якщо, де достатньо велике натуральне число. Покладаючи в (20), , знаходимо зображення Якобі Значення шуканого розв'язку можна уточнити, якщо замість при використати їх значення, обчислене за побудованими асимптотичними формулами. Використання рекурентних формул дає можливість контролювати похибку обчислення.

На основі ортогональних рядів за многочленами Чебишева-Лагерра побудовано способи обернення одно- та двовимірних вольтеррівських згорток, тобто розв'язані одно- та двовимірні інтегральні рівняння

В цьому розділі запропоновано спектральний метод розв'язування звичайних диференціальних рівнянь та одно- і двовимірних інтегральних рівнянь типу згортки. Основою таких методів є побудова рекурентних алгоритмів для обчислення узагальнених спектрів шуканих розв'язків. Рекурентні формули для обчислення узагальнених спектрів можна застосовувати як в прямому, так і в оберненому порядку [6,9,11,36]. Результати розділу застосовуються для фільтрації замірних даних при побудові початково-граничних умов та адаптивних параметрів моделей процесів масопереносу.

В шостому розділі розв'язана задача конвективної гетеродифузії. Приймається, що переміщення домішок у пористому насиченому водним розчином середовищі здійснюється шляхом конвективної дифузії у поровому розчині та дифузії у частково зв'язному стані і супроводжується процесом взаємопереходу домішкових частинок між цими шляхами. Процеси такого роду описуються системою нелінійних некласичних взаємозв'язаних диференціальних рівнянь (3), яка в лінеаризованому варіанті в одновимірному випадку та безрозмірній формі має вигляд

Процес вертикального переміщення домішок в шарі товщини (), який описується системою (24), досліджено при , початкових і граничних умовах.

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

Якщо використати спектральний метод розв'язування на основі рядів в базисі класичних ортогональних многочленів, то розв'язки подаються рядом

Тут функції визначаються через коефіцієнти для, а. Алгоритм побудови розв'язку полягає в наступному: задаються достатньо великі значення та; для на основі результатів, отриманих в другому та третьому розділах, будуються асимптотичні розклади для при і відповідні регуляризуючі алгоритми; для та коефіцієнти визначаються з останньої системи рівнянь; за знайденими відновлюються шукані розв'язки.

Слід зауважити, що асимптотичні розклади для можна будувати як на основі відомих результатів з теорії диференціальних рівнянь, так і на основі зображень Лапласа, використовуючи теореми тауберового типу. Результати обчислень наведено в табличному та графічному виглядах.

В результаті проведеного аналізу випливають наступні висновки. При застосуванні інтегрального перетворення Лапласа основна трудність полягає в проблемі обернення. Однак якщо відомий розв'язок задачі в зображеннях, то для відновлення оригіналу можна застосувати побудовану в розд. 4 адаптивну схему. Поведінка оригіналу при великих та малих значеннях часів використовується для побудови асимптоти узагальнених спектрів, контролю інших способів розв'язування та вибору порядку стиску-розтягу шкали аргументів при застосуванні теорії збурення.

Якщо вхідні дані відомі в дискретних точках і в математичну модель входять малі параметри, то в цьому випадку доцільно використовувати спектральні методи. Похибка вхідних даних істотно впливає на швидкість збіжності рядів Фур'є. Тому в окремих випадках доцільно застосовувати розклади за класичними ортогональними многочленами. Тоді для досягнення однакової точності використовуються часткові суми набагато менших порядків. При цьому природнім чином може бути використана апріорна інформація.

При заданні вхідних даних в аналітичному вигляді доцільним є поєднання операційного методу розв'язування зі спектральним в базисі класичних ортогональних многочленів. Ефективність такого підходу полягає в тому, що навіть якщо задача не може бути розв'язана операційним способом в аналітичному вигляді, то існують способи знаходження асимптот при великих та малих значеннях аргументів, а це, в свою чергу, дозволяє побудувати асимптоти узагальнених спектрів для великих та малих їх порядків. В результаті цього можна побудувати адаптивний за точністю алгоритм розв'язування задач.

Побудова розв'язків на основі методів теорії збурення особливо ефективна для дослідження стаціонарних процесів та розв'язування відповідних характеристичних рівнянь задачі. Основні результати цих досліджень викладені в роботах [2,3,14,18,21,31,32,40,41,44,45,53]

Сьомий розділ дисертації присвячений побудові розв'язків нелінійних диференціальних рівнянь в частинних похідних, якими описуються нестаціонарні процеси руху газу в трубопроводах. Отримані результати використані для оптимізації процесу транспортування газу та гідравлічного розрахунку газотранспортних мереж. При цьому граничні та початкові умови відомі в дискретному вигляді на нерегулярній сітці і з невисокою точністю. Фізичні параметри моделі, як правило, відомі з точністю до порядку і в конкретних практичних випадках для їх визначення необхідно розв'язувати відповідні обернені задачі.

За ізотермічних умов рух газу в газопроводах достатньо повно описується наступною нелінійною системою диференціальних рівнянь

Задачі, сформульовані для системи рівнянь (25), розв'язані як в лінеаризованому виді, так і ітераційними методами. За початкове (нульове) наближення вибирався розв'язок лінеаризованої задачі. Оскільки з практики відомі границі зміни швидкості руху газу, тобто, то лінеаризація проведена заміною кривої січною, що проходить через точки і, та дотичною, паралельною до неї на даному проміжку. При розв'язуванні задач як невідомі функції вважаємо густину масового потоку і тиск.

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

На основі отриманих результатів розв'язано задачі оптимізації процесу транспортування газу в лінійних структурах та розрахунок газотранспортних мереж.

Таким чином у цьому розділі проаналізовано методи розв'язування задач такого типу стосовно їх точності, швидкості збіжності та часу виконання. Останній фактор має велике значення, оскільки при побудові газотранспортних мереж нестаціонарні задачі як прямі, так і обернені необхідно розв'язувати для кожного лінійного елементу газопроводу. Запропоновано спосіб лінеаризації нелінійних систем диференціальних рівнянь в часткових похідних [1,5,15,16,23_26,28,30,37,38,47,49, 50,52].

Восьмий розділ дисертації присвячений дослідженню фільтрації газу в пористих середовищах (зокрема підземних сховищах газу, рівняння (3)) та впливу на цей процес геометричних параметрів середовища та гідродинамік-них параметрів газу. Граничні умови на тиск задаються або нульовими, або будуються на основі обробки заміряних дискретних даних, а його початковий розподіл апроксимується ортогональним рядом, який будується за даними вимірювань. Параметри , , приймаються сталими. Для розв'язування відповідної задачі математичної фізики побудовано ітераційну схему, доведено її збіжність. За початкове наближення приймається розв'язок лінеаризованого варіанту рівняння (3), яке має вигляд

В останньому рівнянні постійні і залежать від параметрів , що входять в рівняння (3). При моделюванні процесів переносу об'єктів природного середовища існуючі способи експериментального збору інформації вимірювання дають можливість отримувати дані про розподіл шуканих величин (температури, вологості, тиску газу, тощо) у дискретному вигляді з невисокою точністю в нееквідистантних точках. Тому як перший етап дослідження відповідних математичних моделей процесів переносу є подання вхідної інформації в параметричній формі.

Для побудови початкового розподілу тиску газу в пласті підземного сховища використовується спектральне представлення функції двох змінних рядом

За заданими значеннями тиску в робочих, контрольних та заміряних свердловинах відновлюється початковий розподіл. До лінеаризованого рівняння застосовується інтегральне перетворення Лапласа-Карсона за часом. Розв'язок знайдено у виді.

На основі проведених досліджень побудовано алгоритм гідравлічної ув'язки системи пласт підземного сховища газу-газозбірний пункт, який дає можливість розв'язувати наступні основні задачі при розрахунку гідродинамічних параметрів процесів закачування та відбирання газу із підземних сховищ.

Задача 1. При відомому розподілі пластового тиску і заданій величині об'ємного відбору визначити тиск на вході в газозбірний пункт.

Задача 2. При відомому розподілі пластового тиску і заданій величині вхідного тиску на газозбірному пункті визначити об'ємний відбір газу зі сховища.

Задача 3. Аналогічно до задачі 2 при відомому розподілі пластового тиску і заданій мінімальній величині вхідного тиску на газозбірному пункті визначити максимальний об'ємний відбір газу зі сховища.

На основі отриманих результатів побудовано способи розрахунку гідродинамічних характеристик роботи підземних сховищ газу, узгоджених з параметрами роботи газотранспортної мережі [4,17,22,27,29,46,51].