Спосіб вимірювання швидкості направленого потоку рідини або газу

Реферат: UA 104765 C2 (12) UA 104765 C2 Заявлений спосіб вимірювання швидкості направленого потоку рідини або газу, заснований на використанні перенесення міток потоком, полягає у використанні первинного вимірювального перетворювача міток в потоці з рівномірно розподіленою вимірювальною базою довжиною L , розміщеною в потоці під відповідним кутом до напряму потоку. Подають вихідний сигнал Rt  первинного вимірювального перетворювача через вторинний вимірювальний перетворювач на спектроаналізатор. Обчислюють поточну за час T функцію спектральної щільності SR t  сигналу Rt  . Визначають частоти ƒi , на яких функція SR ƒ  має мінімум. Нумерують частоти ƒi в порядку зростання від ƒ0  0 до ƒi  1 2, ..., n і обчислюють середню за час T швидкість , потоку VT за формулою VT  2 n  Lƒi  ƒi-1cos  . n i1 UA 104765 C2 5 10 15 20 25 Винахід належить до вимірювальної техніки, призначений для використання в гідрометеорології і може бути використаний для вимірювання витрати рідини або газу в трубопроводах. Відомі кореляційні способи вимірювання швидкості потоку, засновані на оцінці часу перенесення міток (неоднорідності і (або) суспензії) потоком на відстань між двома фіксованими точками по максимуму кореляційної функції [1]. Як мітки використовуються просторові неоднорідності температури, швидкості, пульсації тиску, коефіцієнта теплообміну датчика з середовищем, густина середовища, електропровідність, концентрації різних розчинених речовин і завислих частинок. Недоліком кореляційних способів вимірювання швидкості потоку є труднощі технічної реалізації точкових датчиків міток з високою швидкодією. Відомі спектральні способи вимірювання швидкості потоку в точці по допплерівському зсуву частоти відображеного від міток акустичного сигналу [2]. Недоліком допплерівських способів вимірювання швидкості потоку є залежність результату вимірювання від швидкості звуку в середовищі, що обмежує точність вимірювання. Аналогів заявленого винаходу заявником не виявлено. У основу винаходу поставлено задачу створення способу вимірювання швидкості направленого потоку рідини або газу, сукупністю суттєвих ознак якого досягається технічний результат - підвищення точності, спрощення технічної реалізації способу вимірювання швидкості потоку і розширення областей застосування. Поставлена задача вирішується тим, що використовують первинний вимірювальний перетворювач міток в потоці з рівномірно розподіленою вимірювальною базою довжиною L (контактний або безконтактний), розміщують вимірювальну базу в потоці під відомим кутом α (0 ≤ α ≤ 90°) до напряму потоку, подають вихідний сигнал Rt  первинного вимірювального перетворювача (ПВП) через вторинний вимірювальний перетворювач (ВВП) на спектроаналізатор, обчислюють відомим способом поточну за час T функцію спектральної щільності сигналу SR ƒ  , визначають частоти ƒi , на яких функція SR ƒi  має мінімуми, нумерують частоти ƒi в порядку зростання від ƒ0  0 до i  1, n , обчислюють середню за час T швидкість потоку VT за формулою 30 35 40 45 50 2 n  Lƒi  ƒi-1cos  . (1) n i1 Суть способу пояснюється кресленнями, на яких зображено: - на фіг. 1 - розміщення рівномірно розподіленої бази первинного вимірювального перетворювача міток в направленому потоці; - на фіг. 2 - схема перенесення потоком просторових міток (неоднорідностей) уздовж вимірювальної бази; - на фіг. 3 - схема трансформації просторової неоднорідності міток x  в часовий процес VT  t  в первинному вимірювальному перетворювачі, що здійснює ковзаюче усереднювання; - на фіг. 4 - графік квадрата модуля передавальної функції розподіленого первинного вимірювального перетворювача при постійній швидкості потоку; - на фіг. 5 - графік функції спектральної щільності S  ƒ  , яку мав би сигнал на виході точкового безінерційного первинного вимірювального перетворювача, що сприймає неоднорідності міток, які переносяться потоком; - на фіг. 6 - графік квадрата модуля передавальної функції розподіленого первинного вимірювального перетворювача при нестаціонарній швидкості потоку; - на фіг. 7 - графік функції спектральної щільності SR ƒ  сигналу Rt  на виході розподіленого первинного (і вторинного) вимірювального перетворювача; - на фіг. 8 - трансформація частотного діапазону сигналу Rƒ  в робочому діапазоні вимірюваних стаціонарних швидкостей; - на фіг. 9 - трансформація функції спектральної щільності SR ƒ  сигналу Rt  в робочому діапазоні нестаціонарних швидкостей; - на фіг. 10 - узагальнена структурна схема типового пристрою для здійснення способу; - на фіг. 11 - структурна схема пристрою для реалізації способу при використанні міток температури і коефіцієнта теплообміну датчика з середовищем з розподіленим первинним вимірювачем температури, що підігрівається; 1 UA 104765 C2 5 10 15 20 25 - на фіг. 12 - структурна схема пристроїв для вимірювання профілю швидкості потоку з розподіленим термопрофілеміром з підігрівом (а) і без підігріву (б); - на фіг. 13 - структурні схеми пристроїв для реалізації способу при використанні міток пульсацій тиску з двома датчиками абсолютного тиску (а) і одним датчиком різницевого тиску (б); - на фіг. 14 - структурна схема пристрою для реалізації способу при використанні міток неоднорідності електропровідності з розподіленим чотириелектродним датчиком електропровідності; - на фіг. 15 - структурна схема пристрою для реалізації способу при використанні міток неоднорідності коефіцієнта ослаблення направленого світла на базі прозороміра; - на фіг. 16 - структурна схема пристрою для реалізації способу при використанні неоднорідності коефіцієнта ослаблення направленого звуку з джерелом і приймачем акустичного сигналу; - на фіг. 17 - типові графічні залежності необхідної смуги частот t вимірника і мінімального часу виміру T залежно від довжини L вимірювальної бази при реалізації способу для виміру швидкостей в конкретному діапазоні. Суть пропонованого способу вимірювання швидкості потоку полягає в наступному. У способі використовується властивість розподіленого первинного вимірювального перетворювача виконувати ковзаюче усереднювання в часі переносимих уздовж нього потоком просторових неоднорідностей вимірюваної величини (міток). Розподілена вимірювальна база (РВБ) первинного вимірювального перетворювача міток розміщується в потоці під кутом α до  напряму вектора швидкості V (0 ≤ α ≤ 90°) (фіг. 1). При цьому швидкість міток в потоці уздовж вимірювальної бази V складе V  V cos 1  . (2) Просторовий розподіл міток x  переноситься уздовж вимірювальної бази з швидкістю V (фіг. 2), усереднюється на базі L і перетвориться в часовий сигнал Rt  на виході первинного вимірювального перетворювача відповідно до рівняння згортки (фіг. 3) 1L 1 t 1   x, t  dx    t  dt   r t     t  dt , (3)  L0  t   де  t  - прямокутна функція, рівна 1 на інтервалі t  , t  , і нулю поза цим інтервалом Rt   r 30 35    LV 1 . (4) У операторній формі sin ƒ Rp  Hp p  expjƒ p . (5) ƒ Первинний вимірювальний перетворювач міток з рівномірно розподіленою базою є ланкою ковзаючого середнього на інтервалі часу  , квадрат модуля амплітудно-частотної характеристики (АЧХ) якого рівний H2ƒ  2 sin 2ƒ   2ƒ  2   sin 2ƒLV1 2ƒLV1 2 . (6) Як видно з графіка Hƒ  2 на фіг. 4, модуль АЧХ має нулі на частотах ƒ i  L з інтервалом 2 ƒi  ƒi  ƒi-1,i  1 2,..., причому , 40  2 ƒ i   2 ƒ  i LV 1 . (7) Звідси Vi  V cos   2ƒ iL cos   2Lƒi  ƒi-1cos  . (8) 2 і визначення ƒi , вхідний для первинного вимірювального Hƒ  перетворювача сигнал міток t  повинен мати спектр S  ƒ  , що покриває одну або декілька смуг ƒi  ƒi-1 (фіг. 5). Для цього перед первинним вимірювальним перетворювачем в потоці при необхідності встановлюють спеціальний турбулізатор. Для ідентифікації 45 2 UA 104765 C2 Незалежно від спектра S  ƒ  , спектр вихідного сигналу SR ƒ  матиме нулі на тих же частотах ƒi , що і Hƒ  2 . Оскільки ідентифікація ƒ зажадає деякого часу T , те значення швидкості перенесення міток   уздовж вимірювальної бази за цей час не буде постійним і 5 замість нулів H2ƒ 2 матимуть місце локальні мінімуми (фіг. 6), відповідні яким частоти ƒi ; слід визначати при аналізі функції спектральної щільності SR ƒ  (фіг. 7) вихідного сигналу первинного вимірювального перетворювача R(t), причому 2 SR 2ƒ  H2ƒ S 2ƒ . (9) У цьому загальному випадку ƒi , будуть випадковими величинами, і доцільно узяти їх середнє значення по виразу 10 1n  ƒi  ƒi-1  , n i1 де i  1 2,..., n . , Для середньої швидкості потоку, враховуючи вирази (1) і (8), отримаємо ƒi  2L cos  n  ƒi  ƒi-1  . (11) n i1 Діапазон вимірюваних швидкостей Vmax  Vmin визначає діапазон можливих частот в VT  15     одному вікні ƒi max  ƒi min  V max  V min L1  Vmax  Vmin cos1 L1 . Як видно з останнього виразу, широкосмугість вхідного для ПВП сигналу θ(t) можна   забезпечити за рахунок установки кута α, при   90  отримаємо ƒ1 max  ƒ1min   . Оскільки розрішення по частоті dƒ  T -1 , то для заданого числа рівнів N по вимірюваній швидкості можемо записати 25 30 35    20  N  ƒ1 max  ƒ1min T 1 (12) Для необхідного часу вимірювання отримаємо N (13) T ƒ1 max  ƒ1min З виразу (13) видно, що із збільшенням ƒ   отримаємо T  0 . З виразів (8) і (11) видно, що, в принципі, для визначення швидкості потоку досить визначити ширину будь-якого одного з n частотного вікна ƒi  ƒi  ƒi-1 . Типові просторові спектри основних параметрів водних середовищ в спадаючій частині добре апроксимуються статечними законами вигляду k m , де показник m приймає значення від 1 до 5, тобто є нескінченними [3]. При перенесенні міток потоком щодо первинного вимірювального перетворювача ці просторові спектри трансформуватимуться в тимчасові спектри, також нескінченні. Безконтактні первинні вимірювальні перетворювачі, що використовують випромінювання, можуть бути зроблені майже безінерційними. Отже, для підвищення швидкодії вимірювання швидкості потоку пропонованим способом немає принципових обмежень, окрім чутливості і швидкодії електроніки обрамлення первинного вимірювального перетворювача. Якщо первинний вимірювальний перетворювач має модуль G2ƒ  передавальної функції Gp  , то замість виразу (6) слід використовувати 2 sin ƒ  G2ƒ .(14) ƒ При монотонно спадаючої G2ƒ , що характерний для всіх інерційних перетворювачів, H12ƒ  2 40 положення мінімумів функцій H12ƒ  2 і H2ƒ 2 не зміниться і, отже, частоти ƒi , залишаться тим самим. Розглянемо варіанти реалізації способу при використанні різних сигналів міток в потоці. Узагальнена структурна схема типового пристрою для здійснення способу представлена на фіг. 10. 3 UA 104765 C2 5 10 15 До складу пристрою (фіг. 10) входять: огорожа 1 первинного вимірювального перетворювача (ПВП) 2 з розподіленою базою, вторинний вимірювальний перетворювач (ВВП) 3, частотний фільтр (ЧФ) 4, аналого-цифровий перетворювач (АЦП) 5 і обчислювач (О) 6. Огорожа 1 служить для фіксації напряму потоку і розміщення вимірювальної бази ПВП 2 в потоці під заданим кутом α до напряму вектора швидкості. Контактний ПВП 2 перетворить середнє значення сигналу міток на розподіленій вимірювальній базі в електричний сигнал. У безконтактному ПВП 2 проводиться модуляція мітками потоку інформативного параметра променя якого-небудь випромінювання, наприклад, оптичного або акустичного. Безконтактний ПВП 2 містить джерело і приймач випромінювання, які можуть бути суміщені при використанні відображеного сигналу. Вторинний вимірювальний перетворювач 3 перетворить інформативний параметр сигналу ПВП 2 в сигнал, зручний для фільтрації частотним смуговим або нижніх частот фільтром ЧФ 4 і, якщо це необхідно, подальшого аналого-цифрового перетворення АЦП 5. При використанні фільтра нижніх частот його частота зрізу повинна бути поблизу ƒi max справа, при використанні смугового фільтра - декілька ширше за смугу ƒi  ƒi-1 . Частота дискретизації ƒ 0 АЦП повинна задовольняти нерівності ƒ0  2ƒ i або ƒ0  2nƒ . Вибір max значення n ≥ 1 проводять з наявності спектральних складових сигналу міток в потоці в робочій смузі частот. Обчислювач 6 служить для визначення функції спектральної щільності SR 2ƒ і 20 25 30 35 40 45 50 55 частот ƒ i , i  1, n . У пристрої може бути використаний як 0 6 аналоговий спектроаналізатор для обчислення SR ƒ  , в цьому випадку АЦП 5 відсутній. При використанні цифрового спектроаналізатора як 0 6 на нього може бути покладена частина завдання смугової фільтрації сигналу Rt  . Всі вузли електроніки 3-6 спільно утворюють блок електроніки БЕ 7, який входить до складу всіх пристроїв, що здійснюють викладений спектральний спосіб вимірювання швидкості потоку. Робота типового пристрою полягає в ковзаючому усереднюванні сигналу міток на розподіленій вимірювальній базі, обчисленні функції спектральної щільності сигналу з виходу первинного вимірювального перетворювача, визначенні частот ƒi , на яких ця функція має мінімум, і обчисленні середньої швидкості потоку за формулою (11). Різні варіанти пристроїв для реалізації способу залежать від середовища потоку, природи сигналу міток в потоці, виду первинного вимірювального перетворювача і представлені на фіг. 11 - 16. Технічні характеристики пристроїв для реалізації способу залежать від діапазону вимірюваних швидкостей, параметрів перетворювачів і обчислювача, що входять в їх склад. Розглянемо очевидні для фахівців приклади пристроїв для здійснення запропонованого способу вимірювання швидкості потоку рідини або газу, що відрізняються природою використовуваного сигналу влучний в потоці і типом первинного вимірювального перетворювача цього сигналу в електричний. Відомо [4, 5], що потоки рідин і газів мають неоднорідності температури і швидкості, які можуть служити мітками. Для знімання сигналу міток неоднорідності температури використовують розподілений первинний вимірювальний перетворювач температури 2, який встановлюють на вимірювальній базі під відомим кутом α до вектора швидкості потоку (фіг. 10). Для знімання сигналу міток неоднорідності швидкості використовують встановлюваний розподілений первинний вимірювальний перетворювач температури, що також підігрівається робочим струмом або вбудованим нагрівачем 2 (фіг. 11), який, будучи чутливим до зміни коефіцієнта теплообміну, сприйматиме також мітки неоднорідності фізичних параметрів середовища: щільності, теплоємності, теплопровідності, кінематичній в'язкості. При вимірюванні профілю швидкості по перетину потоку як первинний вимірювальний перетворювач використовують розподілені термопрофілемір 2 з підігрівом або без, який встановлюють аналогічно розподіленому вимірникові температури (фіг. 12). Гідродинамічні неоднорідності в потоці рідини як мітки виявляються в пульсаціях гідростатичного тиску, для вимірювання середнього значення яких на розподіленій вимірювальній базі на її межах розміщують два датчики абсолютного гідростатичного тиску (фіг. 13 а) або два входи одного датчика різницевого тиску (фіг. 13 б). На відміну від відомих вимірників швидкості [5] з вимірюванням перепаду тиску до і після діафрагми, в даному вимірнику обидва датчики встановлюють без або після діафрагми. Якщо середовище потоку має електропровідність, то просторові неоднорідності температури і фізичних параметрів середовища в потоці приведуть до локальних 4 UA 104765 C2 5 10 15 неоднорідностей електропровідності, які як мітки вимірюють розподіленим первинним вимірювальним перетворювачем електропровідності, наприклад, чотириелектродним, потенційні електроди (ПЕ) якого встановлюють на межах вимірювальної бази, а струмові електроди (ТЕ) поза нею, як це показано на фіг. 14. Фактично це відомий чотириелектродний датчик электропроводимости в трубковому виконанні [6], використаний для реалізації способу вимірювання швидкості течії. Гідродинамічні неоднорідності в потоці сприяють створенню неоднорідності розчинених і зважених речовин, що є мітками для випромінювань, що проходять через потік, наприклад, оптичних і акустичних, інформативним параметром яких є коефіцієнт ослаблення світла в певних спектральних вікнах або ослаблення звуку на певних частотах. Використовувані в даний час вимірники показника ослаблення направленого світла (прозороміри) [6], мають розподілену вимірювальну базу між джерелом і приймачем направленого світла або між джерелом-приймачем і триппель-призмою. Такі прозороміри можна використовувати для здійснення спектрального способу вимірювання швидкості потоку, розмістивши їх джерело (И) і приймач (П), вимірювальну базу в огорожі потоку і модуль електроніки (ME) 8 перед БЕ 7 по схемі на фіг. 15. Ослаблення направленого світла в середовищі відбувається згідно із законом Бугера [7], відповідно до якого інтенсивність сигналу I t  на вході (виході) приймача пов'язана з інтенсивністю сигналу на виході джерела I0 t  за формулою 20 I t   I0 t e  x dx , (15) де x  - коефіцієнт ослаблення світла. З виразу (15) можемо записати L  x, t dx  1n I t  / I0 t   yt  . (16) 0  I t   Оскільки I t   I0 t  , то yt   1n 0  - вихідний сигнал прозороміру, причому  I t   25 1L 1   x, t  dx   r t   t  dt . (17)  L0   Переходячи до операторної форми запису виразу (17), отримаємо вираз y t   sin 2ƒ  ƒ   exp j   p . (18) 2ƒ  4  Отже, квадрат модуля амплітудно-частотної характеристики вимірювального перетворювача прозороміра з вихідним сигналом y(t) відповідає виразу (6), а функція спектральної щільності S y 2ƒ  сигналу y(t), аналогічно виразу (9), рівна yp  30 35 2 S y 2ƒ  H2ƒ S 2ƒ (19) і має мінімуми (нулі) на частотах, що фіксуються ƒi , i  1 2,..., n . Далі справедливі формули , (7), (8), (10), (11). При реалізації акустичного вимірника швидкості потоку за спектральним способом передавач і приймач сигналу розміщуються на межах вимірювальної бази. Для вузького акустичного променя інтенсивність звуку на вході приймача I t  пов'язана з інтенсивністю звуку на виході передавача I0 t  за формулою L I t   I0 t  k x, t dx , (20) 0 40 де k x  - просторовий розподіл коефіцієнта загасання. По аналогії з вищевикладеним отримаємо I t  1     t k t dt (21) I0 t     r або в операторній форміy t   yp  sin 2ƒ  ƒ   exp j  k p . (22) 2ƒ  4  5 UA 104765 C2 Для ФСП сигналу y(t) отримаємо S y 2ƒ  H2ƒ Sk 2ƒ , (23) 2 5 10 яка матиме мінімуми (нулі) на частотах, що фіксуються ƒi , i  1 2,..., n . Далі справедливі , формули (7), (8), (10), (11). У структурну схему акустичного пристрою вимірювання швидкості потоку за спектральним способом (фіг. 16) замість вторинного вимірювального перетворювача входить модуль електроніки (ME) 8 у складі приймачів-передавачів, що забезпечує вироблення сигналу yt  для БЕ 7. Здійснюють заявлений спосіб наступним чином. Приклад 1 Вимірювання швидкості течії природних вод, типовий діапазон швидкостей 0,02  5 м  с 1 . Задамося розрішенням за швидкістю 1% (N=100), приймемо кут нахилу розподіленої вимірювальної бази до напряму вектора швидкості потоку α = 45°. Використовуючи вираз (8) можемо визначити перший частотний діапазон ƒ1 max  ƒ1min , ширину частотного вікна  15 20 25 30 35 40 45 50  ƒ  ƒ1 max і по виразу (12) необхідний час вимірювання в першому частотному діапазоні (n=1) при використанні можливих вимірювальних баз різної довжини L від 0,01 до 10 м. Графіки цих функцій ƒв  L і Tв  L в логарифмічному масштабі показані на фіг. 17. З міркувань спадкоємності просторового усереднювання вимірювань з існуючими приладами (наприклад, гідрометричними вертушками) і зручності експлуатації приймаємо L = 0,1 м. З графіка (фіг. 17), і по формулах, визначаємо потрібне ƒв  71 4 Гц і Tв  1 4 с. Вимоги , , наявності смуги частот пропускання більше 71,4 Гц можуть бути задоволені використанням таких первинних вимірювальних перетворювачів: розподілений датчик температури, що підігрівається, розподілений термопрофілемір з підігрівом, датчик різницевого тиску, розподілений датчик електропровідності (для морської води), датчик прозороміру, вимірник коефіцієнтів ослаблення акустичного сигналу. Типовий прийнятий час усереднювання при вимірюванні течій (20 ÷ 100) с може бути забезпечений накопиченням відліків. Приклад 2 Вимірювання швидкості вітру в типовому діапазоні швидкостей 12  60 м  с 1 . , Аналогічно викладеному вище в прикладі 1 приймаємо N=100 і α=0°, розраховуємо функції ƒ   L  , T  L  і будуємо ці графічні залежності на фіг. 17. З міркувань практичної доцільності приймаємо довжину вимірювальної бази L = 0,3 м. Визначаємо ƒ   200 Гц і T  0,51 с Вимога наявності смуги частот пропускання вище 200 Гц може бути задоволене використанням як первинні вимірювальні перетворювачі датчиків різницевого тиску, прозороміру і вимірника коефіцієнта ослаблення акустичного сигналу. Основною гідністю заявленого способу, при будь-якому його технічному втіленні, є простота апаратурної реалізації в обмін на відносну складність обробки вимірювальних сигналів, прийнятну на сучасному рівні техніки. Список джерел інформації: 1. Лангс Ф. Корреляционная электроника: Основы и применение корреляционного анализа в современной технике связи, измерений и регулирования. - Л.: Госсудпромиздат, 1969. - 342 с. 2. Мясников В. И. Измерение расхода и объема воды. Мир измерений. - 2001. -№ 3, 4. - С. 49. 3. Бендат Дж. Основы теории случайных шумов и ее применения. М., «Наука», 1965. - 465 с. 4. Хинце И. О. Турбулентность. М. - Л., Физматгиз, 1063. - 680 с. 5. Макс Ж. Методы и техника обработки сигналов при физических измерениях. В 2-х томах. Пер. с франц. - М.: «Мир», 1983. - Т. 1. - 312 с. 6. Петухов Б. С. Опытное изучение процессов теплопередачи. ГЭИ., М. - Л. - 1952. - 344 с. 7. Дерюгин К. К., Степанюк И. А. Морская гидрометрия. - Л.: Гидрометеоиздат, 1974. - 285 с. 6 UA 104765 C2 ФОРМУЛА ВИНАХОДУ 5 10 1. Спосіб вимірювання швидкості направленого потоку рідини або газу, заснований на використанні перенесення міток потоком, полягає в тому, що використовують огорожу, в якій розташований первинний вимірювальний перетворювач міток в потоці з рівномірно розподіленою вимірювальною базою довжиною L , розміщеною в потоці під кутом 0    90  до напряму потоку, при цьому подають вихідний сигнал Rt  первинного вимірювального перетворювача через вторинний вимірювальний перетворювач на аналоговий спектроаналізатор, обчислюють поточну за час T функцію спектральної щільності SR t  сигналу Rt  , визначають частоти ƒi , на яких функція SR ƒ  має мінімум, нумерують частоти ƒi в порядку зростання від ƒ0  0 до ƒi  1 2, ..., n і обчислюють середню за час T швидкість , потоку VT за формулою 2 n  Lƒi  ƒi-1cos  . n i1 2. Спосіб за п. 1, який відрізняється тим, що як мітки використовують неоднорідності температури в потоці, а як первинний вимірювальний перетворювач - розподілений датчик температури. 3. Спосіб за п. 1, який відрізняється тим, що як мітки використовують неоднорідності температури і коефіцієнта теплообміну в потоці, а як первинний вимірювальний перетворювач розподілений датчик температури, що нагрівається. VT  15 20 7 UA 104765 C2 8 UA 104765 C2 9 UA 104765 C2 10 UA 104765 C2 11 UA 104765 C2 Комп’ютерна верстка В. Мацело Державна служба інтелектуальної власності України, вул. Урицького, 45, м. Київ, МСП, 03680, Україна ДП “Український інститут промислової власності”, вул. Глазунова, 1, м. Київ – 42, 01601 12