Спосіб вивчення складнопобудованих середовищ сейсморозвідкою

Спосіб вивчення складнопобудованих (анізотропних) середовищ, що включає збудження сейсмічних хвиль і їх реєстрацію багатоканальними розстановками, обробку даних та одержання висновку про будову геологічного середовища, який відрізняється тим, що до виконання наземних спостережень на площі робіт відпрацьовують одну або декілька глибоких свердловин поляризаційним методом вертикального сейсмічного профілювання (ПМ ВСП), причому один пункт збудження розташовують поблизу устя свердловин, а решту уздовж радіальних профілів з азимутальним кутом в межах 45°-60° та відстанню пунктів збудження C2 2 (19) 1 3 74914 4 розміщенні джерел на відстанях від свердловини, кінематичні характеристики спостережених хвиль розрахованих спеціальним чином. поля градієнтного середовища в характеристики Далі реєструють падаючу і відповідну їй середовища із середньою швидкістю, для цього відбиту хвилі свердловинними сейсмоприймачами знаходять коректуючі часові поправки як різницю в межах заданого інтервалу профілю спостереміж кінематичними характеристиками, визначенижень при вибраних кутах падіння та виходу хвиль, ми по ГГМ, і обчисленими за середньо зміщують джерела і сейсмоприймачі в сусідні швидкісними характеристиками. Визначають поля позиції для послідовного огляду біля свердловинкоректуючих поправок як функцію часу і ного простору, в процесі якого багатократно віддалення траси від спільного пункту збудження і змінюють кути падіння падаючої хвилі і очікувані спільної глибинної точки, по системі сейсмічних кути виходу відбитої хвилі, та формують сейсмічне профілів для кожного відбиваючого горизонту зображення середовища, що досліджується. одержують усереднені значення поля часу Недоліком способу є неврахування анізотропії реєстрації (Тог) і значення поля ефективної швидкостей, що призводить до помилок у швидкості відбитих хвиль, на якому ізолініями визначенні положення відбиваючих границь. Крім рівних значень відображають форму відбиваючих того необхідність багатократно змінювати кути горизонтів в часовому масштабі (То/2, Vеф). Перепадіння падаючої хвилі суттєво здорожує польові творюють ці поля в структурну карту глибин по роботи. вертикалі і виключають спотворення, пов'язані з Відомі також результати моделювання хвильнахилом площини променів в умовах нестрого'!' ових полів та відповідної інверсії даних в разі орієнтації профілів вхрест і по простяганню порід. анізотропії порід в товщі, що перекриває більш Недоліком способу є врахування анізотропії глибокі горизонти з прямокутним виступом [2]. Посейсмічних швидкостей в процесі інверсії даних, казано, що при врахуванні анізотропії зображення що призводить до спотворення дійсного положенвиступу на динамічному модельному розрізі може ня сейсмічних горизонтів на часових розрізах, на декілька сот метрів зміщуватись відносно його внаслідок чого структурні побудови можуть виявифактичного положення на моделі. Проте результись хибними. Просте перетворення кінематики тати моделювання не ув'язані з реальною методихвильового поля градієнтного середовища в серекою сейсмічних досліджень і, як наслідок, не довище із середньою швидкістю не забезпечує адаптовані для обробки польових матеріалів. необхідної точності структурних побудов. Найбільш близьким до запропонованого з точВ основу винаходу поставлено завдання ки зору вивчення складнопобудованих середовищ підвищити точність та достовірність є спосіб обробки сейсморозвідувальданих даних сейсморозвідки в умовах складнопобудованих [3], що включає реєстрацію сейсмічних хвиль багасередовищ шляхом врахування на токанальними установками по профілях, вибір деталізаційному етапі анізотропії швидкостей та інтервалів профілів для обробки, сумісну обробку залучення повнохвильового моделювання. з даними спільної глибинної точки (СГТ), введення Поставлене завдання досягається тим, що у коректуючих поправок у вихідні дані по профілях, запропонованому способі вивчення складнопобувиділення відбитих хвиль і визначення часу їх дованих середовищ, який включає збудження реєстрації, кута підходу, амплітуди та ефективної сейсмічних хвиль та їх реєстрацію багатоканальшвидкості розповсюдження. В подальшому визнаними розстановками, обробку даних та судження чають надійність виділених хвиль та про будову геологічного середовища, згідно виназакономірність зміни ефективної швидкості шляходу, до постановки наземних спостережень на хом апроксимації значень і часу реєстрації, а таплощі робіт відпрацьовують одну або декілька кож інтервали з рівним значенням градієнтів зміни глибоких свердловин поляризаційним методом ефективної швидкості в часі, значення часу одновертикального сейсмічного профілювання (ПМ кратно відбитих хвиль „Т0"п в точках перетину ВСП), причому один пункт збудження розташовуінтервалів з різними градієнтами ефективної ють поблизу гирла свердловин, а решту - уздовж швидкості. Уздовж профілів у вибраних інтервалах радіальних профілів з азимутальним кутом в мез витриманими градієнтами чи такими, що жах 45-60° та відстанню пунктів збудження від закономірно змінюються, визначають значення гирла свердловини, визначеною за умови, що макінтервальної швидкості, природу зареєстрованих симальний кут падіння сейсмічних променів на хвиль і положення кратноутворюючих границь. границю поділу товщ з різними фізичними харакПотім в межах вибраних інтервалів корелюють теристиками не перевищує критичного, по отризначення часу однократно відбитих хвиль „То", що маних сейсмозаписах обчислюють локальні співпадають з графіком ефективної швидкості на швидкості та кути підходу хвиль до кожного пункту кореляційній схемі, із значеннями часу відбитих збудження, оцінюють параметри ТТІ анізотропії хвиль „Тог" на розрізі СГТ шляхом апроксимації (трансверсально-ізотропного середовища з нахидискретних значень, розподілених в заданому вікні леною віссю симетрії) Томсена і та кути науздовж профіля, отримують часові розрізи хилу осі анізотропії уздовж свердловини, потім спільного пункту збудження і спільної глибинної будують відповідні графіки та індикатриси, після точки, визначають поля ефективної швидкості узчого проектують наземні 2,5D- або 3D-сейсмічні довж профіля і виявляють закономірності їх зміни і дослідження та проводять повнохвильове модезалежність від геологічної будови. Одержані дані лювання з вирішенням прямої та оберненої задач комплексно інтерпретують і будують узгоджену із сейсміки на основі скалярного та вертикального стандартним розрізом СГТ геолого-геофізичну хвильового рівняння, оцінюють очікувані ефекти модель (ГГМ) середовища, перетворюють впливу анізотропії на точність та достовірність 5 74914 6 структурних побудов, в подальшому виконують відклади 15, сейсмічні горизонти 16, бокова стінка наземні сейсмічні дослідження з орієнтацією смуг солі 17, падаючі сейсмічні промені 18. спостережень уздовж довгої осі еліпса перетину Як відомо, в умовах трансверсальноеліпсоїду апроксимації індикатриси анізотропії ізотропного середовища побудова границь з викогоризонтальною площиною та апертурою, узгодристанням середньошвидкісного закону розповженою з результатами свердловинних досліджень, сюдження сейсмічних хвиль забезпечує необхідну отримані матеріали обробляють із залученням точність структурних побудов. міграційного енергетичного оператора Кірхгофа, а В разі більш складного типу симетрії зокрема остаточне положення сейсмічних границь на зобТТІ-анізотропії (трансверсально-ізотропне середораженнях встановлюють шляхом ітераційного навище з нахиленою віссю симетрії) структурні побуближення результатів міграції сейсмозаписів надови можуть суттєво відрізнятися від реальних, що земних спостережень і обробки матеріалів ПМ викличе серйозні проблеми при проведенні пошуВСП до даних глибокого буріння. ково-розвідувальних робіт. Особлива увага повинВ порівнянні з прототипом спосіб на приділятись врахуванню анізотропії швидкохарактеризується наявністю наступних відмінних стей. Ігнорування цієї обставини, судячи з даних ознак: моделювання [2], може призвести до значних - відпрацювання перед постановкою наземних зміщень структурних побудов по латералі. Викодосліджень глибокої свердловини поляризаційним нане моделювання в межах центральної частини методом вертикального сейсмічного Дніпровсько-донецької западини, де потужність профілювання в азимутальній модифікації з досоляного діапіра перевищує 1.5км, показало, що сить значним віддаленням пунктів збудження від зміщення структурних побудов по підсольових гирла свердловини дозволяє вивчити параметри горизонтах може сягати сотні метрів (фіг. 3). Зваанізотропії складнопобудованого середовища і жаючи на те, що точність закладання пошуковозастосувати їх як в процесі повнохвильового морозвідувальних свердловин складає перші сотні делювання так і на етапі обробки та інтерпретації метрів, неврахування впливу анізотропії призведе польових сейсмозаписів; до помилок в розміщенні точок глибокого буріння. Таким чином, основним завданням при - оцінка параметрів ТТІ анізотропії Томсена і реалізації способу є вірне визначення анізотропії і куту нахилу осі анізотропії уздовж свердлошвидкостей гірських порід. Вона може бути визнавини та побудова відповідних графіків та чена шляхом аналізу приведених вступів прямої індикатрис забезпечує обґрунтоване повнохвильохвилі з ближнього та виносних пунктів збудження, ве моделювання та проектування наземних 2,5Dабо більш точно на основі застосування алгоритму або 3D-сейсмічних досліджень; Есмерсоу [4], в якому для отримання швидкостей - проведення повнохвильового моделювання з поздовжніх та поперечних хвиль як функції глибивирішенням прямої та оберненої задач сейсміки ни та кута падіння використовуються на основі скалярного та векторного хвильового багатокомпонентні сейсмічні спостереження. рівняння дозволяє оцінити очікувані ефекти вплиВідомо, що при вертикальному сейсмічному ву анізотропії на точність та достовірність струкпрофілюванні (В СП) більша частина енергії, турних побудов; падаючої поблизу перших вступів, - виконання наземних сейсмічних досліджень з розповсюджується вертикально уздовж прийомної орієнтацією полос спостережень уздовж довгої вісі установки. Тому локальні позірні швидкості еліпсу перетину еліпсоїду апроксимації уздовж установки, обчислені по інтервальному індикатриси анізотропії горизонтальною площиною часу пробігу, дозволяють безпосередньо отримута апертурою, узгодженою з результатами свердвати швидкості в середовищі. ловинних досліджень, дає змогу побудувати При виносних пунктах збудження основна панайбільш достовірні розрізи (зображення); даюча хвиля падає на прийомну установку під - залучення при обробці міграційного енергеневідомим кутом , а позірна швидкість уздовж устичного оператора Кірхгофа забезпечує високу тановки не відповідає швидкості в середовищі. В точність сейсмічних зображень; цьому випадку для визначення останьої необхідно - застосування ітераційних процедур на етапі мати додаткові відомості про середовище. Надетальної обробки та інтерпретації даних приклад, швидкість в середовищі (на заданій дозволяє максимально наблизити результати до глибині) можна обчислити за допомогою позірної даних глибокого буріння і таким чином підвищити швидкості, якщо відомий локальний кут падіння детальність вивчення складнопобудованих хвилі. Така інформація може бути отримана за (анізотропних) середовищ. допомогою трикомпонентних сейсмоприймачів та Суть способу пояснюється рисунком (фіг. 1), здобута з поляризації руху частинок середовища. де прийняті наступні позначення: свердловина 1, Таким чином, аналіз поляризації трикомпонентних профілі 2, виносні пункти збудження коливань 3, даних в часовому вікні поздовжньої хвилі дозволяє ближній пункт збудження коливань 4, еліпси перевизначити значення локального кута падіння хвилі. тину 5 та 6 еліпсоїду апроксимації індикатриси В подальшому значення позірних швидкостей поанізотропії горизонтальною та вертикальною пломножують на косинуси кутів падіння і отримують щинами відповідно, індикатриси 7 та 8 анізотропії швидкості в середовищі. відповідно в горизонтальній та вертикальній плоОписана вище процедура може бути ефективщинах, кут 9 між профілями збудження коливань, но застосована лише до ізольованих вступів бурова вишка 10, точки реєстрації коливань 11 у хвиль. В часових вікнах, де поздовжні і поперечні свердловині, денна поверхня 12, надсольові SV-хвилі інтерферують, час пробігу і поляризацію відклади 13, сольові породи 14, підсольові 7 74914 8 руху частинок середовища при проходженні хвилі По отриманих сейсмозаписах визначають не можна визначить незалежно одна від одної. позірні швидкості та кути падіння поздовжніх і поПри виносних пунктах збудження навіть при перечних SV-хвиль з кожного пункту збудження і помірних кутах падіння поздовжні та поперечні SVотримують локальні швидкості в середовищі для хвилі переходять одна в другу на всіх границях кожного профілю. розділу. Це вимагає застосування коротких вікон Виконують сумісний аналіз локальних швидкоаналізу, що призводить до помилок при визначенні стей і будують відповідні індикатриси в величин. Для зниження можливих похибок горизонтальній та вертикальній площинах. необхідне застосування інтерактивних методів По отриманих матеріалах знаходять параметдослідження. Крім того, обмінні SV-хвилі можуть ри анізотропії розрізу. В якості таких можуть бути бути використані при визначенні швидкості і струквикористані напівкількісні способи оцінки основних тури. параметрів ТТІ-анізотропії (трансверсальноОдночасне обчислення швидкості, поляризації ізотропного середовища з нахиленою віссю і форми хвильових полів дозволяє належним чисиметрії), які найбільшою мірою впливають на реном обробляти різні вступи з хвильовими формазультати міграційних перетворень. До них ми, що перекриваються. відносяться в першу чергу параметри Томсена і В результаті отримують модель, яка предта кут нахилу осі анізотропії (р [5]. ставлена чотирма параметрами, що не залежать Наступним кроком є проектування наземних від частоти (швидкості і кути падіння Р- та SV2,5D- або 3D-досліджень, при цьому полоси хвиль), а також відповідними хвильовими форма(профілі) спостережень орієнтують уздовж довгої ми. В подальшому зазначені параметри моделі вісі еліпсу перетину еліпсоїду апроксимації уточнюють шляхом мінімізації найменшої індикатриси анізотропії горизонтальною площиквадратичної похибки між спостереженими даними ною, що дозволить отримати найбільш достовірні і значеннями, отриманими за допомогою моделі. розрізи (зображення). На фіг. 1а приведена схема проведення польОдночасно проводять повнохвильове модеових досліджень, що передбачає відпрацювання лювання в пакеті Tesseral із застосуванням раніше глибокої свердловини 1 з одного ближнього 4 та визначених параметрів анізотропії середовища. декількох виносних 3 пунктів збудження, розташоДля цього розраховують хвильові поля для кожнованих на трьох профілях 2 з азимутальним кутом 7 го пункту збудження з використанням реальної між ними в 60°. В разі вимог підвищеної схеми відпрацювання наземних профілів та детальності кількість профілів збудження може апріорної моделі середовища. бути збільшена до чотирьох, тобто азимутальний Здійснюють обробку модельних сейсмозаписів кут між профілями становитиме 45°. з використанням відомостей про анізотропію На фіг. 1б зображена типова модель швидкостей та виконують міграцію часових сейсмічного розрізу в умовах соляного діапіризму розрізів згідно середньошвидкісного закону, що та схема ходу променів 18 з пунктів збудження 3 дає змогу попередньо оцінити очікувані ефекти та 4. стосовно впливу анізотропії на точність сейсмічних Виконання досліджень по такій схемі, що побудов. передбачає наявність докритичних зондуючих Проводять наземні сейсмічні дослідження променів у свердловині на глибині кожного (2,5D або 3D) та виконують обробку даних із залусейсмічного горизонту, дозволяє побудувати ченням міграційного енергетичного оператора індикатриси анізотропії швидкості уздовж свердКірхгофа та врахуванням параметрів анізотропії. ловини, які можуть бути апроксимовані деяким Порівнюють результати обробки з даними ВСП та еліпсоїдом. Для зручності аналізу виконують перебуріння і в разі потреби корегують (переглядають) тини 5 та 6 еліпсоїду відповідно горизонтальною вихідні параметри. та вертикальною площинами. Параметри Остаточне положення сейсмічних границь на анізотропії середовища прямим чином пов'язані з розрізах (зображеннях) встановлюють шляхом анізотропією швидкостей, що дозволяє врахувати ітераційного наближення даних наземних спостеїх як в процесі повнохвильового моделювання так і режень і ПМ ВСП до даних глибокого буріння. при обробці польових сейсмічних записів. Приклад реалізації способу Спосіб здійснюється слідуючою послідовністю Для реалізації способу в повному обсязі операцій. необхідно спроектувати і провести польові В межах площі досліджень будують сейсмічні роботи та виконати відповідну обробку сейсмічний (або геологічний) розріз, що проходить даних. через глибоку свердловину, і розраховують (або По тих сейсмічних матеріалах, що є на моделюють) місцезнаходження виносних пунктів сьогодні, спосіб можна випробувати лише у скорозбудження за умови, що кожній сейсмічній границі ченому варіанті. Таке випробування виконане на повинен відповідати пункт, який забезпечує докриодній із площ в межах центральної приосьової тичний (біля критичний) зондуючий промінь. частини Дніпровсько-донецької западини. На денній поверхні розбивають три (або чотиВ процесі проведення пошуковори) профілі з азимутальним кутом в 60° (або 45°) розвідувальних робіт в околиці Тарасівського сота розміщують на них виносні пункти збудження. ляного штоку глибока свердловина, запроектована Відпрацьовують свердловину з трикомпонентна розкриття підсольових нижньопермських ною реєстрацією сигналів з ближнього та виносних відкладів, пройшла по солі і не виконала свого пунктів збудження. геологічного завдання. Такі випадки не поодинокі. 9 74914 10 Після завершення буріння у свердловині були При середньому нахилі осі анізотропії 25 проведені сейсмічні дослідження із п'яти пунктів значення параметра Томсена складає 0,25. збудження (ПЗ) уздовж профілю, який проходив Важливо зазначити, що кути нахилу довгих через свердловину і був орієнтований ортогональосей 23 анізотропії в надсольовій товщі та в солі но до межі соляного штоку (в плані). близькі між собою, в той же час вони близькі до За даними свердловинних досліджень будувакута нахилу осадових верств, що дозволяє при лись як субгоризонтальні так і субвертикальні формуванні швидкісної моделі змінювати кут награниці, які мали б більш детально охарактеризухилу осі анізотропії швидкостей у відповідності з вати геометричні контури соляного тіла. Проте, нахилом надсольових осадових верств. буріння нахиленої свердловини, незважаючи на Менш важливий параметр Томсена може значне відхилення стовбуру від гирла свердловибути заданий виходячи з існуючих кореляційних ни (по горизонталі більш ніж 800 м), не виправдазв'язків для різних типів порід [6], на рівні 0,08. ло очікуваних сподівань. Результати побудови модельних сейсмічних Це спонукало до переінтерпретації матеріалів зображень на основі використання пакету програм як свердловинних, так і наземних сейсмічних спомоделювання Tesseral приведені на фіг. 3, де стережень та відповідного моделювання. прийняті такі позначення: сейсмічні границі в Були детально проаналізовані хвильові поля надсольової товщі 25, сольова частина розрізу 26, за даними ПМ ВСП та зроблені висновки щодо сейсмічні границі в підсольовій частині розрізу, можливості застосування алгоритму Есмерсоу для бокова стінка солі 28, спектри амплітуд 29 та вивчення швидкісної характеристики розрізу, який швидкостей 30. В процесі побудови модельних на даний час розглядається як найбільш досконамігрованих розрізів були використані синтетичні лий та гнучкий стосовно складнопобудованих сейсмограми, з формовані за допомогою Еоб'єктів. Найбільш доцільним для вивчення швидоператора з урахуванням анізотропії (фіг. 3а) та Ккостей у верхній (надсольовій) частині розрізу, а оператора без її урахування (фіг. 3б). також безпосередньо у солі є використання даних Розгляд мігрованого розрізу на фіг. 3а свідчить по пунктах збудження 1 (L=65 м), та 3 (L=630 м). про те, що підсольові границі 27 зміщені в бік солі Результати визначень локальної швидкості 700-800 м по відношенню до бокової стінки 28 на поздовжніх хвиль приведені на фіг. 2, де прийняті модельному розрізі. наступні позначення: 19, 20 - швидкості, визначені При врахуванні анізотропії сейсмічне зобраза матеріалами ПЗ 1 та ПЗ 3 відповідно, 21 - кут ження (фіг. 3б) повністю відповідає моделі. нахилу вектора поляризації хвиль за матеріалами На фіг. 4 приведені результати сейсмічних поПЗ 3, 22 - апроксимація індикатриси швидкостей будов уздовж сейсмічного профілю в межах поздовжніх хвиль за допомогою еліпса, 23, 24 північної околиці Тарасівського штоку. Прийняті поздовжня та поперечні вісі еліпса відповідно. такі позначення: пісковики в нижньопермських Як видно з фіг. 2а, в інтервалі свердловини від відкладах 31, тектонічне порушення 32. 1,0 км до 1,5 км кут підходу 21 хвилі до свердлоФрагмент сейсмічного розрізу, отриманий за вини з пункту збудження 3 досить стабільний і по традиційною технологією, показаний на фіг. 4а. Як абсолютній величині дещо зменшується. На цьому уже зазначалось вище, глибока свердловина, проінтервалі, як і в подальшому, локальна швидкість бурена в районі пікету 28.50, не вийшла із солі і не 19, визначена для ближнього пункту, завжди нижрозкрила, як очікувалось, підсольові ча ніж та 20, що визначена для ПЗ 3. Це, врахонижньопермські горизонти 27. вуючи положення зазначеного пункту збудження, Варіант сейсмічного зображення з урахувансвідчить про те, що напрям короткої осі 24 еліпса ням анізотропії швидкостей в солі та в надсольовій 22 (фіг. 2б), яким апроксимована індикатриса частині розрізу приведений на фіг. 4б, де стрілкою швидкості поздовжніх хвиль, близький до нормалі показані підсольові горизонти. стосовно залягання верств. Порівняння фрагментів сейсмічного розрізу Дані добре узгоджуються з теоретичними уяв(фіг. 4а) та сейсмічного зображення (фіг.4б) леннями. Співвідношення швидкостей при кутах свідчить про їх суттєві відмінності, зокрема, на падіння хвилі від 0° до 20° дозволяє оцінити в останньому підсольові горизонти зміщені більш ніж надсольовій частині розрізу параметранізотропії на 1 км в бік від соляного штоку. Останнє не вхоТомсена 0,25, при куті нахилу осі 30 (фіг. дить в суперечність з результатами попереднього 2б). буріння нахиленої свердловини. В межах солі ділянкою із відносно стабільним Неврахування анізотропії може призвести такутом підходу хвилі до свердловини може вважакож до розриву суцільності зображень, що може тися інтервал 2,0 - 3,3 км, в межах якого значення бути ідентифіковано з тектонічними порушеннями і кута поляризації 21 змінюються від 15° до 0° т. ін. (фіг.2а). При зростанні кута зростає різниця між Таким чином, проведене скорочене випробушвидкостями, визначеними з ПЗ 1 та ПЗ 3. В вання запропонованого способу на моделях та інтервалі 2,8-3,2 км, де кут падіння хвилі при ПЗ 3 реальних сейсмічних матеріалах свідчить про наближаються до нуля, значення швидкостей доцільність його застосування на етапі детального практично збігаються. вивчення складнопобудованих (анізотропних) Приведені дані свідчать про наявність об'єктів. Повномасштабне впровадження способу анізотропії в солі. На фіг. 2в показана дозволить суттєво підвищити ефективність пошуапроксимація дійсної індикатриси швидкостей в ково-розвідувальних робіт на нафту і газ, особлимежах соляного штоку за допомогою еліпса 22. во в умовах солянокупольної тектоніки. 11 74914 12 Перелік використаних джерел інформації: profiles" Geophysics, vol.55 №1 (January 1990), p. 1. Патент України №10797А, G 01V1/00, опубл. 39-50. 25.12.1996 р., Бюл. № 4. 5. Isaac, J.H. and Lawyer, L.c., 1997, Image mis2. J.Helen Isaac and Don C. Lawton image mispositioning due to dipping ТІ media: A physical seispositioning due to dipping Timedia: Aphysical seismic mic modeling study: Geophysics, 64, 1230-1238. modeling study, Geophisics Vol.64, № 0.4 (Juli6. Ryan-Grigor S. Empyreal relationships beAugust 1999); p.1230 - 1238, 17 Figs. tween transverse isotropy parameters and Vp/Vs: 3. Патент РФ № 2148838, G 01V1/30, опубл. Implications for AVO. - Geophysics, V.62. -P.135910.05.2000 р., прототип. 1364. 4. Cengis Esmersoy "Inversion of P and SV waves from multicomponent offset vertical seismic 13 74914 14 15 Комп’ютерна верстка М. Клюкін 74914 Підписне 16 Тираж 26 прим. Міністерство освіти і науки України Державний департамент інтелектуальної власності, вул. Урицького, 45, м. Київ, МСП, 03680, Україна ДП “Український інститут промислової власності”, вул. Глазунова, 1, м. Київ – 42, 01601