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

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

Київський національний університет імені Тараса Шевченка

УДК 550.552.53.553

04.00.22 - Геофізика

Автореферат дисертації

на здобуття наукового ступеня доктора геологічних наук

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

Вижва Сергій Андрійович

Київ - 2004

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

Робота виконана у Київському національному університеті імені Тараса Шевченка.

Науковий консультант: доктор фізико-математичних наук, професор Продайвода Георгій Трохимович, Київський національний університет імені Тараса Шевченка, професор кафедри геофізики.

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

- доктор геолого-мінералогічних наук, професор Лизун Степан Олексійович, ДП "Науканафтогаз", генеральний директор, м. Київ;

- доктор геолого-мінералогічних наук Орищенко Іван Васильович, Інститут геологічних наук НАН України, провідний науковий співробітник, м. Київ;

- доктор геолого-мінералогічних наук, професор Кузьменко Едуард Дмитрович, Івано-Франківський національний технічний університет нафти і газу, головний науковий співробітник, м. Івано-Франківськ.

Провідна установа: Інститут геофізики імені С.І. Субботіна НАН України, м. Київ.

Захист відбудеться "17" червня 2004 р. о 10 год. на засіданні спеціалізованої вченої ради Д026.001.32 Київського національного університету імені Тараса Шевченка за адресою: 03022, м. Київ, вул. Васильківська, 90.

З дисертацією можна ознайомитись у науковій бібліотеці Київського національного університету імені Тараса Шевченка за адресою: 01033, м. Київ, вул. Володимирська, 58.

Автореферат розісланий "14" травня 2004 р.

Вчений секретар спеціалізованої вченої ради Д026.001.32 кандидат геолого-мінералогічних наук, доцент А.В. Сухорада.

Загальна характеристика роботи

Актуальність теми. Наприкінці ХХ - на початку ХХІ ст. спостерігається зростання кількості природних і природно-техногенних катастроф та обсягів спричинених ними загальнолюдських та матеріальних збитків. В основі природних катастроф лежать закономірні геологічні та геофізичні процеси, що відбуваються в оболонках Землі (атмосфері, гідросфері, земній корі, мантії й ядрі) та зазнають впливу космічних факторів і техногенної діяльності. Останні чинники самі можуть спричинити або прискорювати чи посилювати природні катастрофічні явища. Загальновизнаних методів оцінки потенційної небезпеки виникнення природних катастроф на сьогодні не існує. Більш ніж 50-річна (а по землетрусах - більше 100-річна) історія дослідження небезпечних природних і техногенно-природних геологічних процесів (НПТПП) не принесла помітних успіхів щодо їх прогнозу та попередження. Разом з тим, незважаючи на прорахунки та відсутність помітних результатів, проблема моніторингу НПТПП з метою запобігання їх негативних наслідків є однією із найбільш актуальних геофізичних задач. Для успішного вирішення цієї проблеми необхідно шукати нові методологічні та теоретичні підходи.

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

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

Зв'язок роботи з науковими програмами і планами НДР. Дисертаційну роботу виконано в Київському національному університеті імені Тараса Шевченка; вона нерозривно пов'язана з виконанням планів науково-дослідних робіт Київського національного університету імені Тараса Шевченка та Міністерства освіти і науки України, Мінекоресурсів України, Держкоменергозбереження України. Більшість результатів отримані при виконанні науково-дослідної теми "Розробка нових геофізичних технологій пошуків родовищ корисних копалин, геофізичного моніторингу небезпечних геологічних процесів та екологічного стану середовища" (0101U002769) в рамках програми "Надра" Фонду фундаментальних досліджень Київського національного університету імені Тараса Шевченка. Основні теоретичні, методичні та алгоритмічні розробки автора викладено у звітах за темами (№№ держреєстрації 0195U004128, 0196U003806, 0197U003155, 0101U002769).

Удосконалення, впровадження та практична реалізація виконаних автором розробок здійснювались у межах госпдоговірних робіт з ДГП"Укргеофізика", УкрДГРІ, ЗАТ "Геологічний інститут Київського університету".

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

Основні завдання досліджень. Для досягнення поставленої мети необхідно було вирішити наступні задачі.

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

Розробити новий методологічний підхід до вирішення задач геофізичного моніторингу за НПТПП. Провести оцінку сучасного стану геофізичного моніторингу навколишнього природного середовища в Україні та розробити підходи щодо комплексної геодинамічної інтерпретації даних геофізичного моніторингу НПТПП.

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

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

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

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

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

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

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

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

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

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

В дисертаційній роботі використовуються та узагальнюються методи механіки стохастично неоднорідних середовищ, зокрема метод умовних моментних функцій та його узагальнення; методи статистичного моделювання випадкових процесів та полів (метод Монте-Карло) на базі спектральних розкладів; ітераційні методи чисельних розв`язків систем лінійних та нелінійних рівнянь, ультразвуковий інваріантно-поляризаційний метод акустичних досліджень зразків гірських порід; свердловинні радіоізотопні методи досліджень густини та вологості ґрунтів.

Наукова новизна отриманих результатів.

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

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

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

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

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

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

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

Вперше розроблено технологію сейсмоакустичної НПТПП.

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

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

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

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

Особистий внесок здобувача. Основні теоретичні та методичні результати, які винесено на захист, отримані здобувачем особисто. Особистий внесок здобувача в основні роботи, виконані у співавторстві, визначається наступним чином. У роботі [1] автору належить ідея, безпосередня участь в постановці задачі, одержанні, обробці та аналізі результатів польових геофізичних досліджень. У роботах [2, 8, 9] дисертанту належить безпосередня участь у постановці задачі, проведенні лабораторних ультразвукових досліджень, обробці й аналізі результатів. У роботах [3, 13, 17-19, 24, 26] дисертанту належать ідея, постановка задачі, аналіз та узагальнення отриманих даних, підготовка висновків. У роботах [4-7, 11] здобувач брав участь у постановці задачі та одержанні результатів. У цих роботах всі чисельні розрахунки та побудови виконані разом з Продайводою Г.Т. У роботах [10, 14, 16, 21, 22] здобувачеві належить ідея та безпосередня участь у постановці задачі, розробці алгоритмів та аналізі результатів. У роботах [12, 20, 25] дисертант брав участь у постановці задачі, зборі й аналізі матеріалів, підготовці висновків.

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

Науково-технічній конференції "Будівництво в сейсмічних районах України" (Ялта, 1999);

наукових конференціях професорсько-викладацького складу геологічного факультету Київського національного університету імені Тараса Шевченка (Київ, 2000-2001);

XXV Генеральній Асамблеї Європейського Геофізичного Товариства (Ніцца, 2000);

X Сесії Російського акустичного товариства (Москва, 2000);

Міжнародних наукових конференціях "Геофізичний моніторинг небезпечних геологічних процесів та екологічного стану середовища" (Київ, 2000-2002);

ІІ та ІІІ міжнародних конференціях по геоінформатиці (Київ, 2003-2004);

Науково-практичній нараді "Інженерний захист берегів Чорного та Азовського морів України" (Ялта, 2003);

Міжнародній геофізичній конференції "Геофизика ХХІ столетия. Взгляд в будущее" (Москва, 2003);

Міжнародній науковій конференції "Моніторинг небезпечних геологічних процесів та екологічного стану середовища" (Київ, 2003).

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

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

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

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

Основний зміст роботи

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

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

За даними, наведеними в звітах The Center for Research on the Epidemiology of Disasters (CRED), Міністерства України з питань надзвичайних ситуацій та у справах захисту населення від наслідків Чорнобильської катастрофи (МНС), працях провідних вчених відзначається, що кількість природних катастроф у світі за останні 40 років постійно зростає, збільшується кількість загиблих і постраждалих y природних катастрофах осіб (прирости відповідно склали, в середньому, 4,3% і 8,6%), ще швидшими темпами ростуть економічні втрати - 10,4% на рік.

Основними чинниками, які впливають на ріст природних катастроф є:

глобальні процеси космічного та літосферного походження;

зростання чисельності населення планети;

урбанізація;

зростання техногенного впливу на навколишнє природне середовище.

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

Отже, при розробці стратегії прогнозування стану довкілля України необхідно враховувати дві основні складові розвитку НППТП:

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

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

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

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

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

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

У підрозділах 2.1 та 2.2 наведені дані про застосування геофізичних методів до вирішення задач моніторингу небезпечних геологічних процесів. Зокрема найбільш розробленим методично, теоретично та апаратурно є сейсмічний моніторинг. Значний внесок у розвиток методологічної, методичної й апаратурної бази сейсмічних спостережень на території України внесли В.А. Дядюра, О.В. Кендзера, Г.І. Кузьменко, В.Є. Кульчицкий, С.В. Мостовий, Б.Г. Пустовітенко, О.А. Харитонов та інші. Незважаючи на певні досягнення, землетруси є природним стихійним лихом, що на сьогодні не піддається прогнозуванню. Окремо необхідно виділити напрямок, пов'язаний з контролем за геологічним середовищем на території розташування небезпечних промислових об'єктів (атомних електростанцій (АЕС), сховищ радіоактивних (РАВ) та високотоксичних відходів, крупних гідротехнічних споруд, рудних родовищ корисних копалин, підприємств нафтогазового комплексу та нафтогазотранспортної системи).

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

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

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

У підрозділі 2.4 розглядається традиційний підхід до проблеми комплексної геофізичної інтерпретації, в розвиток якого зробили значний внесок Г.Я. Голіздра, Ф.М. Гольцман, Г.І. Каратаєв, С.С. Красовський, Е.Д. Кузьменко, В.І. Старостенко, В.М. Страхов. Він базується на розв'язку прямих і обернених геофізичних задач у рамках прийнятої моделі середовища і поля, а також на кореляційних взаємозв'язках між фізичними параметрами геологічного середовища. Основні труднощі застосування класичного підходу до вирішення задач геофізичного моніторингу обумовлені різною просторовю самоспряженістю джерел геологічного середовища, що визначають геофізичні поля спостережень.

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

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

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

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

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

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

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

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

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

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

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

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

модуль оцінки стійкості об'єкта моніторингу;

модуль оцінки надійності прогностичних оцінок параметрів та можливих матеріальних втрат;

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

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

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

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

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

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

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

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

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

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

У підрозділах 3.4-3.9 наведені дані моделювання акустичних і пружних параметрів для моделей з орієнтованими системами тріщин, різною орієнтацією породоутворюючих мінералів.

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

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

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

Таблиця 3.1

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

Дослідження впливу кристалографічної орієнтації зерен рогової обманки на пружні властивості амфіболіту здійснювалися на прикладі моделі "амфіболіт". Встановлено, що орієнтовані вздовж координатних осей кристалографічні осі рогової обманки (для 40% зерен, ax1, bx2, cx3) формату =3 призводять до акустичної анізотропії 6,4%. Моделі тріщинуватого амфіболіту (SAT-1, SAT-2, SAT-3) з анізотропною матрицею та однаковою концентрацією мікротріщин (с1=0,002) відрізняються величиною коефіцієнта акустичної анізотропії, що коливається від 5,35% до 22,2% (див. табл. 3.1). Пружна симетрія текстур моделей SAT-1, SAT-2, SAT-3 - планальна ромбічна. Стереопроекції ізоліній фазових швидкостей повздовжних хвиль мають ромбічну симетрію, але спільний вплив кристалографічної орієнтації й орієнтованих систем мікротріщин, в залежності від їх концентрації, призводить як до зміни просторового положення локалізованих екстремумів, так і до зміни вказівних поверхонь фазових швидкостей, кількості акустичних осей і повздовжніх нормалей.

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

У розділі 4 описується розроблений метод математичного моделювання впливу тиску і температури на швидкості поширення пружних хвиль у геологічному середовищі.

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

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

Результати інверсії показали, що розраховані й експериментальні криві фазових швидкостей повздовжньої поляризації перидотиту при різних тисках і температурах співпадають. Концентрація мікротріщин спочатку різко зростає у інтервалі температур 25-200С і досягає найбільшої величини при температурі 400С. При подальшому зростанні температури (в інтервалі 400-700С) концентрація мікротріщин зменшується. Отже, уникнути повного припинення процесу зміни форми мікротріщин неможливо. Можна лише досягти деякого стаціонарного стану, що відповідає діючим при даному тискові і температурі локальним напругам і деформованому стану гірської породи та величині діючих термопружних напруг. Це означає, що пористість і проникливість, а також поровий тиск у гірському масиві змінюється протягом часу і залежить від геодинамічних процесів.

У підрозділі 4.3 розглядаються отримані методом умовних моментів чисельні результати математичного моделювання ефективних термопружних властивостей лерцолітової моделі океанічної верхньої мантії в ромбічному наближенні з урахуванням форми орієнтації еліпсоїдальних зерен олівіну й енстатиту вздовж координатних осей. Причому, кристалографічні осі олівіну орієнтовані таким чином: ax1, cx2, bx3, а енстатиту - cx1, bx2, ax3.

Пружна симетрія текстури моделі лерцоліту виявилася аксіальною ромбічною. Вона задовільно узгоджується з експериментальними даними для дуніту Хороман і дуніту Твін-Сістерс, пружна симетрія яких також аксіальна ромбічна.

Азимутальна анізотропія швидкості квазіповздовжніх хвиль океанічної верхньої мантії за даними моделювання подана на рис. 4.3. На цей графік винесені результати сейсмічних спостережень на ділянках "Квартет" і "Флора" (Тихий океан, на захід від Лос-Анджелеса) з роботи Р.Райта, що обчислені за формулою Бакса.

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

Лерцолітова модель:

1 - глибина 10 км; 2 - глибина 25 км; 3 - глибина 75 км; 4 - глибина 200 км. Польові спостереження: 5 - ділянка "Квартет"; 6 - ділянка "Флора".

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

Термопружні напруги розтягання на глибинах 40-120 км досягають величини 3,5-3,68 МПаК-1. Внаслідок релаксації термопружних напруг, навіть при відсутності зовнішніх зсувних напруг, у зернах олівіну і енстатиту можуть розвиватися пластичні деформації. Розраховані величини локальних термопружних напруг на односпрямовано орієнтованих еліпсоїдальних зернах енстатиту. В напрямку осей Х1 та Х2 діють девіаторні напруги розтягнення, які швидко зростають з глибиною, а в напрямку осі Х3 - девіаторні напруги стиснення. Величина останніх збільшується з глибиною швидше. Отже, для формування текстури лерцоліту в термодинамічних умовах океанічної верхньої мантії існують досить сприятливі умови. Якщо взяти до уваги, що понижені в'язкі властивості астеносфери обумовлені частковим плавленням (вагова частка розплаву 1%), то витримати такі напруги протягом тривалого часу лерцоліт не зможе. Зокрема, при в'язкості =2020 Пуаз, сила волочіння для переміщення плит оцінюється в 10,1 МПа. Діючі девіаторні напруги на еліпсоїдальних зернах енстатиту на глибині 200 км досягають величини 5,5 ГПа.

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

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

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

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

Пружна постійна третього порядку *2 для моделей кварцового пісковику при зміні концентрації мікротріщин формату =10-4 змінюється дещо інакше. Для сухих мікротріщин вона досягає максимуму при концентрації с1=0.1% (*2=200 ГПа), а потім різко зменшується. Однак, для водонасичених вона змінює знак і досягає мінімуму при концентрації с1=0.4% (*2=-11700 ГПа).

Проведені чисельні розрахунки акустопружних ефектів модулів плагіограніту із хаотично орієнтованими "дископодібними" мікротріщинами (із щільністю мікротріщин =0,05) в умовах різного просторового напруженого стану на глибині 15 км континентальної земної кори. Розглянуті три моделі:

а) модель плагіограніту G-1, для якої горизонтально стискаючі напруги були прийняті однаковими ., а вертикальна складова тензора напруг ., де PST - літостатичний тиск;

б) модель G-2, для якої горизонтально стискаючі напруги були прийняті неоднаковими, тобто . і ., a вертикальна складова тензора напруг .;

в) модель G-3, для якої горизонтальна напруга розтягання . була більша, ніж горизонтальна напруга стискання ., а ..

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

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

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

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

У розділі 6 проводяться дослідження можливостей застосування методів статистичного моделювання (методів Монте-Карло) випадкових процесів та полів до задач геофізики.

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

На основі отриманих середньоквадратичних оцінок наближення розроблено алгоритм для чисельного моделювання стаціонарних гауссівських періодичних випадкових процесів (підрозділ 6.1).

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

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

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

Визначено відповідні експериментальним даним значення параметра с=5 та обчислено спектральні коефіцієнти цієї функції:

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

При необхідності запропонована методика дозволяє збільшити кількість точок чисельного моделювання випадкового поля на потрібний порядок. Детальність такої побудови дозволяє прийняти рішення про необхідність проведення додаткових досліджень на території проммайдинчика РАЕС. При порівнянні даних моделювання з результатами експериментальних спостережень густини отримані відхилення знаходяться в межах точності польових методів.

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

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

У розділі 7 приводяться результати застосування розроблених алгоритмів статистичного моделювання випадкових процесів дво- та тривимірних полів для моніторингу та прогнозування карстово-суфозійних процесів на території Рівненської АЕС.

З точки зору інженерної геології район розміщення Рівненської АЕС характеризується несприятливими умовами і відповідно зі схемою карстологічного районування знаходиться в Центральному районі Поліської карстової області.

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

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

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

На території РАЕС протягом останніх років відмічається підвищення рівня підземних вод. Передусім, це відбувається за рахунок втрат з інженерних комунікацій. При значних відмінностях, порівняно з підземним водами, в складі та температурі, порушенні природного фільтраційного режиму швидкість розчинення крейди може змінюватися в багато разів. При зменшенні густини скелету до 1,3103 кгм-3 ймовірність процесів карстоутворення різко збільшується. На рис. 7.1 наведені результати прогнозу густини крейдяної товщі, виконаного на основі прогнозу зміни усереднених значень за 29 свердловинами. Прогноз по кожній із свердловин проводився окремо на період до 2006 року на основі наявних рядів спостережень, потім будувалася спільна карта за допомогою пакету Surfer 7.0. На прогнозній карті з'являються зони в північно-західній та східній частинах майданчика, які поки що при свердловинних вимірах не виявлені. Крім цього, підтверджуються наявні аномалії, зменшення розмірів та амплітуди яких в подальшому не передбачається.

Висновок

Основні результати дисертаційної роботи зводяться до наступного.

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

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

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

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

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

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

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

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

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

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

...