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

8

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

НАЦІОНАЛЬНА АКАДЕМІЯ НАУК УКРАЇНИ

ІНСТИТУТ ГЕОФІЗИКИ ім. С.І. СУББОТІНА

ДУБОВЕНКО ЮРІЙ ІВАНОВИЧ

УДК 550.831:550.838

ВИЗНАЧЕННЯ КОНТАКТНОЇ ГРАНИЦІ ЗА ЗНАЧЕННЯМИ ПОХІДНИХ ЛОГАРИФМІЧНОГО ПОТЕНЦІАЛУ НА ІСТОТНО ОБМЕЖЕНИХ МНОЖИНАХ (ЛІНЕАРИЗОВАНА ПОСТАНОВКА)

04.00.22 Геофізика

АВТОРЕФЕРАТ

дисертації на здобуття наукового ступеня

кандидата фізико-математичних наук

Київ 2004 р.

Дисертацією є рукопис.

Робота виконана у відділі глибинних процесів Землі і гравіметрії Інституту геофізики ім. С.І. Субботіна Національної академії наук України.

Науковий керівник: доктор фізико-математичних наук, професор

ЧОРНИЙ Арнольд Володимирович

Інститут геофізики, головний науковий співробітник

Офіційні опоненти:

доктор фізико-математичних наук,

КОРЧАГІН Ігнат Миколайович,

Інститут геофізики, провідний науковий співробітник

кандидат фізико-математичних наук,

ПЕТРОВСЬКИЙ Олександр Павлович,

науково-технічна фірма „Біпекс Лтд”, директор, старший науковий співробітник

Провідна установа: Національна гірнича академія, кафедра геофізики, МОН України, м. Дніпропетровськ

З дисертацією можна ознайомитись у бібліотеці Інституту геофізики ім. С.І. Субботіна НАН України.

Вчений секретар

спеціалізованої вченої ради,

доктор геологічних наук М.І. Орлюк

ЗАГАЛЬНА ХАРАКТЕРИСТИКА РОБОТИ

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

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

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

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

Гранична задача про відновлення потенціалу за модулем його градієнта та її використання в геодезичній гравіметрії і в теорії інтерпретації гравітаційних і магнітних аномалій (1996-2000 рр., № держ. регістр. О 196 U 004773). логарифмічний потенціал контактний задача

Класи єдиності розв'язку обернених задач теорії потенціалу (2001-2005 рр., № держ. регicтр. 0101 U 000286).

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

Основні завдання досліджень:

1. Визначення інтегральних співвідношень для “лінеаризованої” оберненої задачі структурної гравіметрії та класу єдиності її розв'язків за умов збурених похибками вхідних даних у рамках теорії наближеного вирішення умовно коректних задач.

2. Аналіз характерних властивостей “лінеаризованого” інтегрального рівняння для контакту, яке дає розв'язок оберненої задачі, узгоджений з даними про похибки вхідних даних.

3. Конструктивне з'ясування необхідних і достатніх умов однозначного розв'язання задачі в аналітичному вигляді за умов малих збурень.

4. Реконструкція пропозицій Малкіна, Сенька і Андреєва щодо уточнення наближень Нумерова на базі сучасних уявлень про розв'язання некоректних задач.

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

6. Розробка алгоритмів, позбавлених залежності розв'язку від деякої середньої глибини до поверхні контакту та програмного забезпечення для пп. 5-6.

7. Перевірка якості і надійність отриманих алгоритмів на тестових даних та вивчення впливу похибок методу і вхідних даних на точність відновлення контакту.

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

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

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

Особистий внесок автора. Самостійно отримано результати, що описані у всіх розділах, за винятком п.п. 1.2.3, 2.1, 2.5, 3.3. Окремі зображення ітераційних процесів для обчислення наближень контакту (п.п. 2.4, 3.1), дослідження їх властивостей (п. 1.3.1) та умов існування розв'язку (п. 2.3.1) задачі запозичено з попередніх публікацій А.В. Чорного. Результати, отримані в п.п. 4.2, 4.3, розвинуто на основі спільних з ним публікацій. Програмна реалізація та тестові обчислення й аналіз повністю здійснені автором.

Апробація результатів дисертації. Положення досліджень, що несуть основне смислове навантаження та результати чисельного моделювання, доповідались на 28-й сесії Міжнародного семінару ім. Д.Г. Успенського “Вопросы теории и практики геологической интерпретации гравитационнх, магнитнх и лектрических полей” (м. Київ, 29 січня - 2 лютого 2001 р.), ІІ-й Міжнародній конференції для молодих вчених “Геофізичний моніторинг небезпечних геологічних процесів та екологічного стану геологічного середовища” (Київ, 8-10 жовтня 2001 р.).

Публікації. Основні результати за темою дисертаційного дослідження опубліковано в 4 статтях в фахових журналах “Геофізичний журнал”, “Доповіді НАН України” та 2-х збірках тез.

Структура і обсяг роботи. Дисертація складається із вступу, чотирьох розділів, висновку, переліку використаної літератури, бібліографічного коментаря та двох додатків. Загальний обсяг роботи 205 сторінок, в тому числі 167 найменувань літературних джерел, 30 малюнків та 3 таблиці, обсяг додатків - 35 сторінок. Наукові інтереси автора формувались протягом 1996-2003 р.р. в процесі плідної співпраці з науковим керівником А.В. Чорним, якого автор шанує як вчителя і щиро вдячний йому за діяльну участь й цінні поради на усіх етапах роботи, які допомогли автору вийти за межі власних пристрастей й інтересів. Автор воліє висловити глибоку вдячність академіку НАН України В.І. Старостенку за увагу і всебічну підтримку, складає щиру подяку докт. фіз.-мат. наук В.М. Шуману та канд. фіз.-мат. наук А.М. Тупчієнку за консультації на заключному етапі та канд. фіз.-мат. наук О.В. Легостаєвій, канд. геол.-мін. наук Т.П. Єгоровій і наук. співр. К.П. Барановій за щире ставлення та сприяння у підготовці цєї праці. Окремо варто подякувати сім'ї та любій дружині Владиславі за розуміння і терпіння.

ЗМІСТ РОБОТИ

1. МАТЕМАТИЧНА МОДЕЛЬ ЛІНЕАРИЗОВАНОЇ ЗАДАЧІ ВИЗНАЧЕННЯ КОНТАКТУ

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

,

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

Для плоскої контактної задачі на основі введеної моделі середовища справедливе рівняння, що пов'язує параметри середовища й поля

. (1)

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

(2)

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

, , (3)

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

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

Висновок 1. У класі функцій нелінійне інтегральне рівняння (2) 1-го роду можна замінити з точністю до квадрату ухилення лінійним інтегральним рівнянням Фредгольма 1-го роду

, . (4)

Обравши у (4) замість визначеного у фіксованій точці значення будь-яке інше, скажімо, значення , при , отримаємо нове лінійне рівняння

, , . (5)

Рівняння (4), (5) містять повільно спадаючі при ядра Пуассона. Для забезпечення задовільної точності розв'язків потрібно задавати функцію на “майже нескінченній” множині, а виміряне поле отримане на істотно коротких профілях. Для подолання цієї вади є альтернатива:

здійснити високоточну екстраполяцію значень з короткого на довгий інтервал, після чого обчислити наближений розв'язок рівняння;

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

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

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

,

, (6)

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

Теорема 1. Якщо розв'язок нелінійного інтегрального рівняння (2) для контактної поверхні належить до , тo функція , що визначається як сума розв'язків лінійних інтегральних рівнянь 1-го роду

, , (7)

при , відрізняється від розв'язку не більше, ніж на квадрат відхилення , .

Але можна визначити контакт на короткому відрізку, виходячи лише з одного рівняння.

Висновок 2. Якщо розв'язок нелінійного інтегрального рівняння (2) належить до класу , то його лінійне наближення визначається з точністю до квадрату відхилення , , як розв'язок лінійного інтегрального рівняння Фредгольма 1-го роду:

,

. (8)

Ще один альтернативний шлях для визначення лінійного наближення дає вираз

, (9)

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

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

, , (10)

відрізняється від розв'язку не більше, ніж на квадрат відхилення .

Властивості лінійних операторів системи (10) характеризує така теорема.

Теорема 3. Якщо функція контакту , то лінійний інтегральний оператор задачі (7) обмежений, неперервний і компактний у просторах , .

Аналогічний результат справедливий і для операторів (4), (5). Множину допустимих розв'язків лінеаризованої контактної задачі характеризують такі обмеження:

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

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

3) апріорі обрано гіпотетичні оцінки і , , ;

4) контактна границя розділяє два однорідні середовища з відомим перепадом густин ;

5) спостережені в m точках значення поля містять випадкову похибку заданої інтенсивності ; .

2. МЕТОДИКА РОЗВ'ЯЗУВАННЯ ЛІНЕАРИЗОВАНОЇ ЗАДАЧІ ЗА ПОЛЕМ, ЗАДАНИМ НА НЕОБМЕЖЕНІЙ МНОЖИНІ

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

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

Вдосконалено класичні способи розв'язання контактної задачі (2). Спосіб Нумерова має низьку точність і забезпечує нульове наближення , для якого отримано оптимізовані поправки.

Теорема Андреєва. Якщо розв'язок нелінійного інтегрального рівняння для контактної границі (2) належить до класу , то його лінеаризоване наближення виражається у вигляді

, (11)

і від істинного відхиляється не більш, ніж на при .

Глибокий аналіз формули (11) породив наступні за нульовим наближення контакту, що мають вигляд

,

,

,

(12)

.

Теорема Малкіна. Корекція кожного з лінеаризованих наближень (12) розв'язку рівняння (2) для контактної границі на величину поправки , де

, (13)

підвищує точність наступних наближень на порядок порівняно з початковими, тобто до величини . Практично достатньо точне наближення поправки забезпечує величина

,

де деяке “середнє” значення поля на осі .

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

, (14)

генерує наближення , які ухилятимуться від точного розв'язку не більш, ніж на величину .

Як узагальнення класичних способів Сенька та Маловичка уточнення контакту запропоновано новий ітераційний процес, названий процесом Нумерова-Маловичка:

, (15)

, .

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

Розв'язок лінійного інтегрального рівняння (5) відшукуємо у вигляді границі послідовно сті , генерованої ітераційним процесом

, , (16)

який будемо називати процесом Лаврентьєва-Андреєва. Для спрощення дослідження проблем розв'язності згаданого процесу ітерацій вдамось до попереднього еквівалентного перетворення виразу (16).

Лема 1. Якщо послідовність функцій , генерована ітераційним процесом (16), то еквівалентна їй послідовність генерується за допомогою послідовних наближень

. (17)

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

де .

Ця теорема стверджує існування розв'язку , а швидкість збіжності спадає з наростанням глибини h і номера ітерації. Для узагальненого рівняння (4) існує узагальнений процес ітерацій Лаврентьєва-Андреєва,

, , (18)

для якого справедливі результати, описані лемою 1 та теоремою 5, але замість h підставлено значення в деякій фіксованій точці .

Для змістовної практичної реалізації в рамках першого постулату альтернативи (значення поля сили тяжіння , на майже “нескінченних” інтервалах) запропоновано ефективну схему екстраполяції значень поля за межі профілю, отриману через розклад в ряд Тейлора шуканої функції, яка входить до складу інтегралу Пуассона:

, (19)

де . Для відновлення контакту за полем, заданим на майже необмежених множинах, потрібно розв'язати два лінійні інтегральні рівняння Фредгольма 1-го роду - (19) та (16) або (18).

3. ВИЗНАЧЕННЯ КОНТАКТУ ЗА ПОЛЕМ, ЗАДАНИМ НА ІСТОТНО ОБМЕЖЕНИХ МНОЖИНАХ

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

, , (20)

, ,

який надалі називатимемо процесом Лаврентьєва-Чорного. Так само, як і для процесу (16), вдаємось до еквівалентного зображення ітерацій.

Лема 2. Якщо послідовність функцій породжується ітераційним процесом (20), то еквівалентна їй послідовність генерується за допомогою послідовних наближень

. (21)

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

,

зі швидкістю геометричної прогресії, яка оцінюється нерівністю

, де .

Прямо можна довести лише збіжність наближень , оскільки нерівність справедлива для будь-яких й , а виконується лише за обмеження , тому збіжність наближень проілюстровано лише чисельно. Співставлення ітераційних процесів Лаврентьєва-Андреєва (16) і Лаврентьєва-Чорного (20) виявляє, що другий процес збігається швидше за перший майже вдвічі, але їх швидкості майже однакові за порядком. Практично процес (20) має перевагу над “андреєвським”: завдяки швидкому затуханню по горизонталі ядер перетворень (20) для відновлення з певною точністю наближень потрібно знати значення поля на істотно коротших профілях, ніж при використанні інтегралу Пуассона.

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



, ,

(22)

для якого справедливі такі ж твердження, що описані лемою 2 й теоремою 6 для ітерацій (20). Для останнього рівняння модуль неперервності має оцінку

е мажоранта виразу

за умови , а аналітично продовжене значення гармонічної функції з рівня на рівень .

Висновок 3. Кожен з послідовності компактних операторів

...

на обмеженій множині функцій майже проектує на компакт .

Звести до єдиного знаменника “швидкісні” характеристики ітерацій (16), (18) та (20), (22) можна, якщо прискорити збіжність двох перших. Ця процедура поліпшує чисельну стійкість алгоритмів, до яких її застосовують. Нехай деяка послідовність наближень є функцією знаменника геометричної прогресії, , а q є, власне, знаменником прогресії, наприклад для послідовності ітерацій (20). Для прискорення ітерацій пропонується така схема . Для ітерацій (17) отримаємо “прискорені” наближення у вигляді

(23)

, .

Цю схему легко поширити на інші ітераційні процеси. Збіжність наближень (23) описує теорема.

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

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

Теорема 8. Якщо за малого поля, породжені “паралельними” контактами з класу , незначно ухиляються одне від одного в смислі нерівності , то самі границі незначно відрізняються одна від одної в смислі , де .

4. ОСОБЛИВОСТІ ЧИСЕЛЬНОГО РОЗВ'ЯЗАННЯ КОНТАКТНОЇ ЗАДАЧІ НА ОБМЕЖЕНИХ МНОЖИНАХ

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

Характеристики прямих операторів задачі

(24)

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

Скінченновимірну апроксимацію системи (24) здійснено методом проекцій. Це дозволяє зменшити розмірність системи, оскільки вхідна матриця симетрична і додатно визначена. В підсумку для операторів рівнянь (16) і (20) відповідно маємо:

,

. (25)

Похибки апроксимації операторів операторами оцінюються, відповідно, як

, .

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

,

. (26)

В ході регуляризації за Тихоновим ця система редукується до вирішення рівняння Ейлера

(28)

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

, (29)

Значення отриманого поля, збурені 5%-ним шумом, служать вхідними даними для вирішення оберненої задачі. Під час тестування вивчено чисельну збіжність другого рівняння системи (7) (рис. 1) та поведінку інтегральних ядер (6). Останні нерівнозначні за точністю і швидкодією Працездатність алгоритмів (16), (18) проілюстровано на модельних прикладах з параметрами: довжина профілю ; крок , коливання контакту ; перепад густин ; апріорна глибина контакту ; амплітудні величини ; ; .

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

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

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

Ітерації закінчують по досягненні заданої середньоквадратичної нев'язки чи по виконанні заданої кількості кроків ітерації. Похибки отриманих розв'язків не перевищують заданих неточностей поля. Здійснено порівняння розроблених в роботі і відомих класичних алгоритмів Чорної, Нумерова, Маловичка. Більша ефективність алгоритмів Лаврентьєва-Андреєва, Лаврентьєва-Чорного порівняно з класичними схемами зумовлена кращим ступенем обумовленості систем алгебраїчних рівнянь. Чисельні експерименти виявили, що запропонований алгоритм (16) переважає за точністю розв'язок (2) Чорної майже удвічі, а (20) майже утричі. Реалізація ітерацій на основі рівнянь (8) і (9) дозволяє досягти таких же результатів за вдвічі меншого обсягу обчислень. Окремі здобутки в цьому плані ілюструє рис. 2. Дослідження мають перш за все теоретичний інтерес, але алгоритми можна застосувати у складі автоматизованих систем для моделювання і експрес-аналізу гравіметричних даних на маловивчених територіях з неускладненою будовою.

ОСНОВНІ РЕЗУЛЬТАТИ РОБОТИ

У дисертації наведено теоретичне узагальнення і нове вирішення оберненої задачі потенціалу для контактної границі в лінеаризованій постановці. Головні наукові і практичні досягнення роботи:

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

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

3. Для розв'язання рівнянь кожного типу задач запропоновано збіжні ітераційні процеси Лаврентьєва-Андреєва (16) і Лаврентьєва-Чорного (20) та їх узагальнені аналоги (18) і (22). Практичною перевагою ітерацій (20), (22) є швидша збіжність та можливість задання значень контакту на короткому профілі. Збіжність усіх ітерацій обґрунтовано за новим способом, що базується на еквівалентному зображенні підінтегральних виразів в ітераційних процесах через узагальнені функції.

4. В ході дослідження вищезгаданих ітерацій поліпшено ітераційні способи уточнення нумерівського наближення контакту, запропоновані свого часу Андреєвим, Малкіним і Сеньком та запропоновано новий спосіб Нумерова-Маловичка і обґрунтовано стійкість уточнень. Ці способи можна використати для аналітичних апроксимацій поля в автоматизованих системах інтерпретації.

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

6. Для підвищення надійності обчислень запропоновано ефективний спосіб прискорення збіжності ітерацій на підставі математичних характеристик їх збіжності.

7. Для чисельного вирішення задачі розроблено регуляризуючі алгоритми на базі концепції регуляризації Тихонова. Модернізовано знаходження оптимального параметра регуляризації , мінімізовано похибки округлень. Створено пакет програм і перевірено на тестових задачах. Порівняння чисельних розв'язків довело перевагу алгоритму (20) над алгоритмом (16). Досягнута точність розв'язку при похибках прямого оператора і спостережень задовільна для практичних застосувань.

ОСНОВНІ ПОЛОЖЕННЯ ДИСЕРТАЦІЇ ОПУБЛІКОВАНО В НАСТУПНИХ РОБОТАХ:

Дослідження оберненої задачі потенціалу для контактної поверхні // Геофіз. журн. - 2002. - 24, № 3. - С. 77-92 (співавтор А.В. Чорний).

До теорії структурної задачі гравіметрії в комплексній площині // Доп. НАН України. - 2002, № 4. - С. 145-149 (співавтор А.В. Чорний).

Відновлення контактної границі в шаруватому середовищі // Геофіз. журн. - 2002. - 24, № 6. - С. 36-41.

Уточнення деяких способів наближеного визначення контактної межі // Доп. НАН України. - 2002, № 12. - С. 99-103 (співавтор А.В. Чорний).

Линейная граничная задача восстановления контакта по значениям поля, заданного на существенно ограниченных множествах // Вопросы теории и практики геологической интерпретации гравитационных, магнитных и электрических полей: Материалы 28-й сессии Международного семинара им. Д.Г. Успенского, Киев, 29 января - 2 февраля 2001 г. - Москва: ОИФЗ РАН, 2001 г., - С. 35-36.

Про структурну задачу гравіметрії // Геофізичний моніторинг небезпечних геологічних процесів та екологічного стану геологічного середовища: Матеріали ІІ-ої Міжнародної конференції для молодих вчених, Київ, 8-10 жовтня 2001 р. - Київ: Видавництво КУ, 2001 р., - С. 60-61 (співавтор А.В. Чорний).

АНОТАЦІЇ

Дубовенко Ю.И. Определение контактной границы по значениям производных логарифмического потенциала на существенно ограниченных множествах. - Рукопись.

Диссертация на соискание ученой степени кандидата физико-математических наук по специальности 04.00.22 геофизика. - Институт геофизики им. С.И. Субботина НАН Украины, Киев, 2004.

Диссертация посвящена разработке теории и алгоритмов решения задачи восстановления контактной границы при условии, что поле представлено значениями производных логарифмического потенциала, заданными как на неограниченных так и на существенно ограниченных профилях. Исследованы основные математические модели обратной задачи логарифмического потенциала для контактной поверхности. Отмечено, что параметрическая зависимость решения задачи от плотности активного слоя и от средней глубины или асимптоты контакта есть следствие рассматриваемых моделей слоистой среды, а не принципиальной характеристикой задачи для неограниченных тел в отличие от задачи для тел ограниченных. В случае поля первого типа предложены эффективные схемы экстраполяции значений поля за пределы профиля и ускорения сходимости итерационных процессов Для поля, заданного на ограниченных множествах, обоснован ряд предложенных способов восстановления контактной границы как граничных функций последовательностей решений линейных интегральных уравнений первого рода с компактными операторами, имеющими быстро спадающие интегральные ядра типа Шварца. Установлена корректность решения полученных уравнений на специально введенном классе Нумерова и разработан регуляризирующий алгоритм с повышеной численной устойчивостью. Осуществлена серия численных экспериментов для изучения влияния погрешностей на дискретные аналоги полученных алгоритмов.

Ключевые слова: логарифмический потeнциал, контактная граница, обратная задача гравиметрии, интегральное уравнение, регуляризирующий алгоритм, итерация, ускорение сходимости, приближение.

Дубовенко Ю.І. Визначення контактної границі за значеннями похідних логарифмічного потенціалу на істотно обмежених множинах. - Рукопис.

Дисертація на здобуття наукового ступеня кандидата фізико-математичних наук за спеціальністю 04.00.22 геофізика. - Інститут геофізики ім. С.І. Субботіна НАН України, Київ, 2004.

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

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

Dubovenko Yu.І. Соntact interfaces definition by the inversion of the logarithm potential derivatives given on the essentially restricted sets. - Manuscript.

Thesis for a candidate's degree by speciality 04.00.22 geophysics. - Institute of Geophysics named after S.I. Subbotin, NASU, Kyiv, 2004.

The theses for a candidate's degree are devoted to development of a theory and algorithms concerned with the problem for contact interface recovery provided that the gravity field is represented by the logarithm potential derivatives values given as on the infinited profiles as on the essentially restricted ones. In case of first type field the effective schemes for field values extrapolation beyond the profile limits and acceleration of iteration convergency have been proposed. As for the field, received on the restricted sets, for the contact border renovation the series of methods proposed treated as limit functions of the succession of the solutions of 1^st kind linear integral equations with compact operators having the fast decreasing integral cores of Schwarz's type are substantiated. A correctness of received equations solutions is established on the specially introduced Numerov class and regularising algorithm with enforced numerical stability is worked out.

Key words: logarithm potential, contact interface, gravity inversion, integral equation, regularising algorithm, iteration, convergency acceleration, approximation.

Размещено на Allbest.ru

...