Є ще 5 сторінок.

Дивитися все сторінки або завантажити PDF файл.

Формула / Реферат

Спосіб визначення фактичних характеристик парку відцентрових нагнітачів газотранспортної системи, що включає вимірювання тиску і температури газу на вході і виході за кожним нагнітачем, а також вимірювання числа обертів вала кожного нагнітача та густини компримованого газу, який відрізняється тим, що одночасно за кожним k-им нагнітачем експлуатованого парку з NK газоперекачувальних агрегатів виконують NL серій з NI періодичних вимірювань зазначених параметрів агрегатних автоматик, формують для кожного і-го вимірювання l-ої серії вектор вимірювань , що включає значення температури на вході і виході нагнітача , , тиску газу на вході і виході нагнітача ,  ,частоти обертання вала нагнітача  та густини газу :

,

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

,

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

,

де  - розрахункове значення величини j-го узагальненого дефекту на момент початку l-ої серії вимірювань k-го нагнітача,  - фільтроване значення швидкості збільшення величини j-го узагальненого дефекту на l-ій серії вимірювань  k-го нагнітача за час його наробітку між і-им i і-l вимірюваннями,  - час наробітку ВЦН, що вимірюють від початку серії, при цьому значення  і  визначають по закінченні кожної серії вимірювання і=NI, розв'язуюючи наступну систему рівнянь

,

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

для j=1,...,c-1, c+l,...ND, далі шляхом лінійної інтерполяції розрахункових характеристик вибраного виду узагальненого нормованого дефекту політропного коефіцієнта корисної дії  та зведеного ступеня стиску  для обчисленого значення , яке дорівнює , та масиву значень зведеної об'ємної продуктивності , де  і NT - кількість елементів масиву, визначають відповідні масиви значень  та , на основі яких розраховують апроксимуючі поліноми фактичних зведених характеристик парку ВЦН газотранспортної системи .

Текст

Спосіб визначення фактичних характеристик парку відцентрових нагнітачів газотранспортної системи, що включає вимірювання тиску і температури газу на вході і виході за кожним нагнітачем, а також вимірювання числа обертів вала кожного нагнітача та густини компримованого газу, який відрізняється тим, що одночасно за кожним k-им нагнітачем експлуатованого парку з NK газоперекачувальних агрегатів виконують NL серій з NI періодичних вимірювань зазначених параметрів агрегатних автоматик, формують для кожного і-го вимірювання l-ої серії вектор вимірювань Xk, i, l , що включає значення 3 83696 де Dk,0,l, j - розрахункове значення величини j-го узагальненого дефекту на момент початку l-ої сеф рії вимірювань k-го нагнітача, VD - фільтроваk,l, j не значення швидкості збільшення величини j-го узагальненого дефекту на l-ій серії вимірювань VD k,l, j k-го нагнітача за час його наробітку між іим i і-l вимірюваннями, tнар - час наробітку ВЦН, що вимірюють від початку серії, при цьому значення Dk,0,l, j і VD k,l, j визначають по закінченні кожної серії вимірювання і=NI, розв'язуюючи наступну систему рівнянь ì ¶Fk,l -1, j =0 ï ï ¶Dk,0,l, j , í ¶F k,l -1, j ï =0 ï ¶VD k,,l -1, j î де величина F ,l -1 j k , дорівнює NI ф виFk,l -1 j å (D k,i,l - 1 j - D r = )2 , а VD , , k,i,l -1, j k,l, j = 1 i значають зі співвідношення ф ф = b1 × VD VD k,l, j + b 0 × VD , де b1, b0 k,l, j k, l-1 j , коефіцієнти фільтра, потім, по закінченні NL серій періодичних вимірювань для кожного k-го нагнітача, вибирають вид найбільш ймовірного узагальненого нормованого дефекту с на основі критерію Винахід відноситься до вимірювань у газотранспортній області, а саме до знаходження необхідних даних для ефективного управління газотранспортною системою і призначена для визначення фактичних зведених характеристик відцентрових нагнітачів (ВЦН) газоперекачувальних агрегатів (ГПА) на компресорних станціях (КС) магістральних газопроводів. Ефективне керування газотранспортною системою можливо тільки на основі оперативної інформації про фактичні характеристики експлуатованого енергоустаткування, що включає в першу чергу фактичні характеристики ВЦН. В даний час фактичні характеристики визначаються спеціалізованими організаціями при проведенні теплотехнічних іспитів ГПА. Проведення цих іспитів регламентується ГОСТ 20440-75 "Установки газотурбинные. Методы испытаний" та відповідними інструкціями з визначення технічного стану ГПА. Більшість інструкцій, також як і ГОСТ 2044075, випущені більше 25 років тому і не забезпечують необхідні для ефективного керування значення точності та вірогідності зумовлених показників та характеристик. До того ж проведення теплотехнічних іспитів є дорогою процедурою. Їх проведення прив'язують до виконання окремих видів ремонтів. Тому необхідна оперативна інформація про 4 мінімуму середньоквадратичного відхилення узагальненого нормованого дефекту від його прогнозованого значення Dr ,i,l, j , який задовольняє наk ступну умову NL NI NL NI l = 1i = 1 l = 1i = 1 å å (Dk,i,l,c ( X k,i,l ) - Dr ,i,l,c )2 < å å (Dk,i,l, j ( Xk,i,l ) - Dr ,i,l, j )2 k k для j=1,...,c-1, c+l,...ND, далі шляхом лінійної інтерполяції розрахункових характеристик вибраного виду узагальненого нормованого дефекту політропного коефіцієнта корисної дії hр k,c (Qпр, Dk, c ) та зведеного ступеня стиску пол eр k, c (Q пр, Dk,c ) пр для обчисленого значення Dk,c , яке дорівнює Dk,0,NL,c , та масиву значень зведеної об'ємної продуктивності {Qпр }, t де t = 1 NT і NT - кількість елементів масиву, визна, чають відповідні масиви значень e пр та { } {hпол }, на основі яких розраховують апроксимуk, t k, t ючі поліноми фактичних зведених характеристик парку ВЦН газотранспортної системи ìhp ü, ìep (Q )ü, k 1, NK . = í пол k (Qпр )ý í пр k пр ý î þ î þ фактичні характеристики на поточний або заданий момент часу всіх експлуатованих відцентрових нагнітачів відсутня. Відомі методи визначення фактичних характеристик відцентрових нагнітачів, які ґрунтуються на ряді припущень про незмінність їхніх окремих характеристик або визначеному характері зсуву фактичних характеристик відносно паспортних, також як і методи, що ґрунтуються на експериментальне виявлених взаємозв'язках зміни окремих показників процесу політропного стиску газу в нагнітачі зі зміною його газодинамічних характеристик, не забезпечують необхідну точність. Відомий метод ідентифікації фактичних (реальних) характеристик ВЦН в умовах КС на основі використання співвідношень, які поєднують коефіцієнт напору y , коефіцієнт витрати j та політропний коефіцієнт корисної дії (ККД) hпол ВЦН у будь-якій точці газодинамічних характеристик y y* k1 = -q = f1(q) , hпол h *пол k 2 = y - q × y * = f2 ( q) , q = j / j * , де k1,k2 - комплекси, що є функціями відносного зсуву режимів за витратою q для різних типів ВЦН. Ці комплекси 5 83696 визначають по всій сукупності газодинамічних характеристик всіх експлуатованих типів ВЦН. Індекс "*" відноситься до параметрів, що відповідають фіксованій точці максимального політропного ККД на газодинамических характеристиках. Метод припускає, що об'ємна продуктивність Q * nn та коефіцієнт витрати j в точці характеристики, що відповідає максимальному політропному ККД ВЦН відомі. Значення об'ємної продуктивності Q n визначають виходячи з припущення відносної стабільності характеристики ВЦН у координатах "зведена відносна внутрішня потужність é N i ù ê ú ë r н û пр зведена об'ємна продуктивність Q nnn " é Ni ù ê ú = j1(Qn nn ) . Вимірюючи тиск і температур у ë rн û пр газу на вході і ви ході нагнітача Рн,Тн,Рк ,Тк , частоту обертання його вала n та знаючи газову постійну компримованого газу R, коефіцієнт стискальності газу на вході ВЦН Zn та номінальну частоту обертання n0 обчислюють коефіцієнт А 2 k -3 æ n ö , ç ÷ ç n ÷ k - 1 z n R( Tk - Tn ) è 0ø k 1 1 де , = × k - 1 h пол m T A = 10 -1 æ T ö æ æ P öö k - показник m T = lgç k ÷ × ç lgç k ÷ ÷ , ç T ÷ ç ç P ÷÷ k -1 è н ø è è н øø псевдоізоентропи, m Т - температурний показник політропи. Потім розв'язують систему двох рівнянь ìé N i ù = A × Q n nn ïê ú ïë r н û пр ï , з якої визначають знаí ïé N i ù = j 1(Q n nn ) ïê r н ú ïë û пр î чення Qпр. Значення об'ємної продуктивності наn гнітача розраховують як Q n = Q nnn . Далі, n0 використовуючи співвідношення, що зв'язують j і Q n , а також y і ступінь стиску e визначають j*, y*, h *пол . Задаючи ряд значень qj, за допомогою виразів, що визначають значення комплексів k1,k2 , обчислюють відповідні значення коефіцієнтів витрати j j , коефіцієнтів напору y j і політропного ККД hпол j . Потім, використовуючи вирази, що пов'язують безвимірні та звичайні координати газодинамічних характеристик, розраховують характеристики ВЦН у координатах e - Q n , hпол - Q n та N i - Q n [1]. r Даний метод визначення фактичних характеристик ВЦН в умовах КС, так само як і спосіб, що 6 заявляється, включає вимірювання тиску і температури газу на вході і виході за кожним нагнітачем, а також вимірювання числа обертів валу кожного нагнітача та густини компримованого газу. Однак цей метод недостатньо точний через використання цілого ряду припущень: про незмінність значень комплексів, що пов'язують коефіцієнти напору, витрати і політропного ККД при наявності різних видів проточної частини ВЦН; про відносну стабільність характеристики ВЦН у координатах "зведена відносна внутрішня потужність - зведена об'ємна продуктивність" та ряду інших. Наведені значення оцінок розбіжності розрахункової витратно-напірної характеристики відносно експериментальної в 4% є неприпустимими. Відомий також метод визначення газодинамічних характеристик нагнітачів на КС на основі розрахунку безвимірних характеристик нагнітача, політропного ККД h пол та політропного коефіцієнта напору y пол від коефіцієнта витрати f . Безвимірні характеристики розраховуються при іспитах нагнітачів на повітрі на декількох режимах за витратою при різних еквівалентних частотах обертання. Вхідною інформацією для розрахунку є дані вимірювань наступних параметрів: температура зовнішнього повітря, перепад тиску повітря на витратомірній шайбі, розрідження повітря перед нагнітачем та його тиск за нагнітачем, температура повітря перед і за нагнітачем, барометричний тиск, частота обертання ротора і перепад тиску на конфузорі, встановленому у всмоктувальній камері нагнітача. Для обчислення газодинамічних характеристик нагнітача за безвимірними характеристиками виконують наступну послідовність дій. Діапазон значень коефіцієнта витрати від мінімального значення (помпажу) до максимального (режим максимальної продуктивності) розбивають на 12 точок. Для кожного значення коефіцієнта витрати Ф j , j = 112 , визначають за безвимірними , характеристиками нагнітача відповідні значення hпол(Ф j) і Yпол( Ф j) . Потім обчислюють витрату газу на вході в нагнітач Q j = pu2D 2Ф j / 4 , де u2 2 кругова швидкість на периферії робочого колеса, D2 - діаметр робочого колеса та визначають різницю ентальпій газу на вході інj і виході ikj як u2 . Далі, за допомогою Di j - ikj - iнj = yпол j × 2 hпол j ітераційної процедури за значеннями тиску на виході нагнітача Рk, температурі газу на вході в наh гнітач Tн і політропного ККД полj із використанням рівнянь політропного процесу стиску в нагнітачі підбирають значення тиску на вході в нагнітач Pнj таким чином, щоб значення ентальпії ikj у відповідній точці (Pk,T k) забезпечувало знайдене значення Di j = ikj - iнj . Кінцевою операцією обчислення газодинамічних характеристик є розрахунок зведеного ступеня стиску та політропного ККД у функції 7 83696 від зведеної об'ємної продуктивності за результатами зазначених обчислень для всіх 12 точок коефіцієнта витрати Фj [2]. Даний метод визначення газодинамічних характеристик нагнітачів на КС за допомогою безвимірних характеристик нагнітача також як і спосіб визначення фактичних характеристик парку відцентрових нагнітачів газотранспортної системи, що заявляється, включає вимірювання тиску і температури газу на вході і виході за кожним нагнітачем, а також вимірювання числа обертів вала кожного нагнітача і густини компримованого газу. Однак цей метод має низьку точність і низьку вірогідність визначення фактичних характеристик через те, що безвимірні характеристики нагнітача визначаються одноразово та у процесі експлуатації їхня стабільність у часі, незалежність їх від розвитку де фектів проточної частини не підтверджені. Найближчим за технічною суттю аналогом, який обрано у якості прототипу, є метод визначення газодинамічних характеристик відцентрового нагнітача відповідно з яким знаходять два незалежних коефіцієнта, що враховують зсув фактичних характеристик відносно паспортних при постійному значенні однієї з термогазодинамічних величин відповідно до використовуваної однопараметричної термогазодинамічної моделі процесу стиснення газу. У якості незалежної змінної вибирають зведену об'ємну продуктивність ВЦН Qnp та за умови Qnp=idem визначають наступні коефіцієнти: - коефіцієнт, що враховує "зсув" характеристики зведеної внутрішньої відносної потужності П æ Ni ö ç ÷ çr ÷ è l ø nр æ Ni ö ; ç ÷ K Ni = çr ÷ H æN ö è l ø пр i÷ ç çr ÷ è l ø пр - коефіцієнт, що враховує "зсув" характеристиh ки політропного ККД hпол hпол = пол , де симhП пол вол "п" означає визначення параметра за відповідною паспортною характеристикою ВЦН при значенні Qnp=idem. Коефіцієнти зсуву за іншими характеристиками ВЦН знаходять на підставі цих двох коефіцієнтів. Значення параметрів æ N i ö ç ÷ çr ÷ è l ø пр та hпол визначають на основі термогазодинамічної моделі процесу стиснення газу, що використовує наступні вхідні параметри: значення температури газу на вході і ви ході ВЦН t1 ,t 2, тисків газу на вході і виході ВЦН Р 1, Р 2, частоти обертань його ротора n, компонентну сполуку газу та об'ємну продуктивність ВЦН. Термогазодинамічну модель процесу стиснення газу описують наступною системою співвідношень: - різниця знтальпій газу на виході і вході ВЦН (повний напір) Dh=f1{P1 ,P2,t1 ,t2); потенційна робота стиску газу (політропний напір) 8 w1,2=f 2{P1,P 2,t1,t2); w, політропний ККД hпол 12 ; Dh масова продуктивність Q=f3(Q,P1 ,t1), потужність ВЦН Ni=Q×Dh; зведена внутрішня відносна потужність éNi ù ê ú = f 4 (Ni , P1, f1, П) ; ë r l û пр паспортний політропний ККД hП = f5 (Q, П) ; пол паспортне значення зведеної внутрішньої відП æN ö носної потужності ВЦН ç i ÷ = f 6 (Q 3 , П ) . ç r ÷ è ø Фактичні характеристики політропного ККД та зведеної внутрішньої відносної потужності визначають як добуток відповідних коефіцієнтів "зсуву" на значення паспортних характеристик [3]. Даний метод визначення газодинамических характеристик відцентрового нагнітача також як і спосіб визначення фактичних характеристик парку відцентрових нагнітачів газотранспортної системи, що заявляється, включає вимірювання тиску і температури газу на вході і виході за кожним нагнітачем, а також вимірювання числа обертів вала кожного нагнітача і густини компримованого газу. Однак цей метод має низьку точність і низьку вірогідність визначення фактичних характеристик за наступними причинами: передбачається використання у якості вхідного параметра інформацію про об'ємну продуктивність ВЦН, яка в умовах компресорної станції відсутня, в той же час обчислення об'ємної продуктивності за допомогою звужуючого пристрою на вході ВЦН пов'язано з нерегламентованою похибкою, що може досягати до 20% від змінюваної величини; фактичні характеристики визначають за разовими вимірюваннями його параметрів та при фіксованому значенні зведеної об'ємної продуктивності, що не дозволяє зменшити вплив випадкових складових похибок каналів вимірювання параметрів ВЦН на результати обчислень; фактичні характеристики обчислюють за допомогою коефіцієнтів "зсуву", постійних для всього діапазону зміни зведеної об'ємної продуктивності ВЦН без урахування виду та величини дефекту проточної частини, який викликав відхилення фактичних характеристик від паспортних. В основу корисної моделі поставлена задача в способі визначення фактичних характеристик парку відцентрових нагнітачів газотранспортної системи шляхом послідовних серій періодичних вимірювань параметрів ВЦН із використанням статистичної обробки і фільтрації результатів обчислень, вибору виду найбільш ймовірного узагальненого нормованого дефекту на основі критерію мінімуму середньоквадратичного відхилення узагальненого нормованого дефекту від його прогнозованого значення, забезпечити підвищення точності та вірогідності визначення фактичних характеристик парку ВЦН газотранспортної системи. Поставлена задача вирішується за рахунок того, що у відомому способі визначення фактичних характеристик парку відцентрових нагнітачів газо 9 83696 транспортної системи, що включає вимірювання тиску і температури газу на вході і виході за кожним нагнітачем, а також вимірювання числа обертів вала кожного нагнітача і густини компримованого газу, відповідно до корисної моделі одночасно за кожним k-им нагнітачем експлуатованого парку з NK газоперекачувальних агрегатів роблять NL серій з NI періодичних вимірювань зазначених параметрів агрегатних автоматик, формують для кожного і-го вимірювання 1-ої серії вектор вимірювань Xk, i, l , що включає значення температури на вході і ви ході нагнітача Tвх k,i,l , Tвих k,i,l , тиску газ у на вході і ви ході нагнітача Рвх k,i,l , Рвих k,i,l частоти обертання валу нагнітача n k ,i,l та густини газу rk ,i,l : { чення зведеного ступеня стиску eпр( Xk,i,l ) і політропного коефіцієнта корисної дії hпол(Xk,i, l) відповідно до політропного методу Шульца, для кожного k-го нагнітача, і-го вимірювання l-ої серії, і потім на основі розрахункових характеристик політропного коефіцієнта корисної дії hр (Q , D ) полk , j пр k, j зведеного ступеня стиску р e (Q , D ) , пр k , j пр k, j j = 1,ND , де ND - кількість видів узагальненого нормованого дефекту, які задаються в функції від величини і виду узагальненого нормованого дефекту Dl , у якості якого приймають зважену суму конкретних видів дефектів проточної частини відцентрового нагнітача та що однозначно визначає відхилення розрахункових характеристик від відповідних паспортних для різних видів і величин дефектів проточної частини нагнітача складають систему з двох рівнянь ìe ( X ) = eр (D ,Q ) ï прк k, i, l прk, j k,i, l, j пр k,i,l, j , í ïhполк ( Xk, i,l ) = hр (Dk, i,l, j, Qпр k,i, l, j ) полk , j ï î з якої за всією множиною векторів вимірювань Xk,i,l визначають множину значень узагальнених { } нормованих дефектів {Dk,i,l,j}, Dr ,i,l, j = Dk ,0,l, j + tнар × VD фl, j k k, , де Dk,0,l, j - розрахункове значення величини jго узагальненого дефекту на момент початку l-ої серії вимірювань k-го нагнітача, VD ф - фільтроk,l, j ване значення швидкості збільшення величини j-го узагальненого дефекту на l-ій серії вимірювань VD k,l, j k-го нагнітача за час його наробітку між іим і-l вимірюваннями, tнар - час наробітку ВЦН, що вимірюють від початку серії, при цьому значення Dk,0,l, j і VD k,l, j визначають по закінченню кожної серії вимірювання і=NI, розв'язуюючи наступну систему рівнянь ì ¶Fk,l -1, j =0 ï ï ¶Dk,0,l, j , í ¶F k,l -1, j ï ï ¶VD k,,l -1, j = 0 î } Xk,i,l = Tвх k,i,l , Tвих k,i,l , Pвх k, i, l, Pвих k, i, l, n k,i,l , rk,i,l , розраховують коефіцієнти стискальності, коефіцієнти ізобаричної стискальності, поправки до теплоємності на основі модифікованого рівняння стану Бенедикта-Вебба-Рабіна, далі з урахуванням розрахованих значень коефіцієнтів стискальності і поправок до теплоємності обчислюють зна та 10 де k= 1 NK , , i= 1 NI , l= 1 NL , j= 1,ND та величини їхніх про, , ü гнозованих значень ìD r í k,i,l, j ý , які розраховують î þ виходячи з пропорційного збільшення значення узагальненого нормованого дефекту у функції від часу його. наробітку, вимірюваного від початку серії, за допомогою виразу де величина F ,l -1 j k , дорівнює NI ф виFk,l= å (D k,i,l - 1 j - D r )2 , а VD -1 j , , k,i,l -1, j k,l, j = 1 i значають зі співвідношення ф ф , де b1, b0 = b1 × VD VD k,l, j + b 0 × VD k,l, j k, l-1 j , коефіцієнти фільтра, потім, по закінченні NL серій періодичних вимірювань для кожного k-го нагнітача вибирають вид найбільш ймовірного узагальненого нормованого дефекту с на основі критерію мінімуму середньоквадратичного відхилення узагальненого нормованого дефекту від його прогнозованого значення Dr , який задовольняє наk,i,l, j ступній умові NL NI NL NI l =1i =1 l =1i =1 å å (Dk,i,l,c ( Xk,i,l ) - Dr ,i,l,c )2 < å å (Dk,i,l, j( X k,i,l ) - Dr ,i,l, j )2 k k для j=1,...,c-1, c+1,...ND, далі шляхом лінійної інтерполяції розрахункових характеристик обраного виду узагальненого нормованого дефекту політропного коефіцієнта корисної дії hр пол k,c (Qпр, Dk, c ) та зведеного ступеня стиску eр пр k, c (Q пр, Dk,c ) для обчисленого значення Dk, c , яке дорівнює Dk,0,NL,c та масиву значень зведеної об'ємної { } продуктивності Qпр , де t = 1 NT і NT - кількість , t елементів масиву, визначають відповідні масиви значень e пр та hпол на основі яких розра { k, t } { k, t } ховують апроксимуючі поліноми фактичних зведених характеристик парку ВЦН газотранспортної ü ì p ü системи ìhp í пол k (Qпр )ý, íeпр k (Qпр )ý, k = 1, NK . î þ î þ Технічний результат, якого можна досягти при використанні винаходу, виражений у тому, що за 11 83696 безпечується підвищення точності та вірогідності визначення фактичних характеристик парку ВЦН газотранспортної системи. Причинно-наслідковий зв'язок між сукупністю ознак винаходу і те хнічним результатом простежується в тому, що нові ознаки: використання послідовних серій періодичних вимірювань параметрів ВЦН із використанням статистичної обробки та фільтрації результатів обчислень, вибір виду найбільш ймовірного узагальненого нормованого дефекту на основі критерію мінімуму середньоквадратичного відхилення узагальненого нормованого дефекту від його прогнозованого значення, при взаємодії з відомими ознаками, а саме вимірюванням тиску і температури газу на вході і виході за кожним нагнітачем, а також вимірюванням числа обертів вала кожного нагнітача і густини компримованого газу забезпечують прояв нових технічних властивостей, таких як зменшений вплив випадкових складових похибок каналів вимірювання параметрів ВЦН на результати вимірювань, малу похибку обчислення об'ємної продуктивності за рахунок не використання звужуючого пристрою, врахування виду та величини узагальненого дефекту проточної частини ВЦН. Це дозволяє одержати очікуваний технічний результат, а саме: підвищення точності та вірогідності визначення фактичних характеристик парку ВЦН газотранспортної системи. При цьому підвищення точності та вірогідності визначення фактичних характеристик ВЦП газоперекачувального агрегату в запропонованому способі забезпечується за рахунок наступних ознак винаходу. 1. З метою зменшення впливу випадкових складових похибок каналів вимірювання параметрів ВЦН на результати обчислення його фактичних характеристик згідно з формулою винаходу "роблять NL серій з Nl періодичних вимірювань зазначених параметрів агрегатних автоматик, формують для кожного i-го вимірювання l-ої серії вектор вимірювань Xk,i,l". Вимога про незмінність значення приведеної об'ємної продуктивності ВЦН, так само як вимога про незмінність інших його параметрів під час проведення серії вимірювань не пред'являється. 2. Змінюється склад використовуваних параметрів відцентрового нагнітача за якими розраховуються витрати компримованого газу. Згідно з формулою винаходу кожний поточний і-ий вектор вимірювання Xk,i,l включає "значення температури на вході і виході нагнітача Tвх k,i,l, Твы х k,i,l, тиску газу на вході і виході нагнітача Рвх k,i,l, Рвы х k,i,l частоти обертання валу нагнітача nk,i,l та густини газу r k,i,l: { } Xk,i,l = Tвх k,i,l, Tвых k,i,l, Рвх k ,i, l, Рвых k,i,l, n k,i,l, r k,i,l ". Інформацію про об'ємну продуктивність відцентрового нагнітача не використовують. В умовах КС об'ємну продуктивність ВЦН обчислюють за допомогою звужуючого пристрою на його вході з похибкою до 20%. Така величина похибки не дозволяє використовувати цей параметр для визначення фактичних характеристик ВЦН із необхідною похибкою в декілька відсотків. 3. Використання "розрахункових характеристик політропного коефіцієнта корисної дії 12 h p (Q пр , D k, j ) полk , j та зведеного ступеня стиску p ( Q пр , D k, j ), j = 1, ND , де ND - кількість видів пр k , j узагальненого нормованого дефекту, які задаються в функції від величини і виду узагальненого нормованого дефекту Dk,j, у якості якого приймають зважену суму конкретних видів дефектів проточної частини відцентрового нагнітача, що однозначно визначає відхилення цих характеристик від відповідних паспортних для різних видів і величин дефектів проточної частини нагнітача" у відповідності з формулою винаходу, дозволяє відмовитися від припущення прямого "зсуву" фактичних характеристик ВЦН від відповідних паспортних. Введення в розрахунок припущень, які не ґрунтуються на експериментальних випробовуваннях та теоретичних дослідженнях з оцінкою похибки, внесеної в результати розрахунку, істотно знижує вірогідність результатів визначення фактичних характеристик ВЦН. 4. Вибір найбільш ймовірного узагальненого нормованого дефекту проточної частини ВЦН із загальної кількості дефектів введених у розгляд, що охоплює всі можливі комбінації конкретних видів дефектів (збільшення зазору в ущільненні покриваючого диска, підрізу лопаток робочого колеса нагнітача і т.д.) ґрунтується на фізичних процесах зміни останніх. Так за сукупністю обробки всіх вимірювань параметрів ВЦН "вибирають вид найбільш ймовірного узагальненого нормованого дефекту с на основі критерію мінімуму середньоквадратичного відхилення узагальненого нормованого дефекту від його прогнозованого значення", що "розраховують виходячи з пропорційного збільшення значення узагальненого нормованого дефекту у функції від часу його наробітку, вимірюваного від початку серії". Такий вибір виду найбільш ймовірного узагальненого нормованого дефекту є одним з основних факторів, що забезпечують максимальну вірогідність визначення фактичних характеристик ВЦН. 5. Величину найбільш ймовірного узагальненого дефекту ВЦН обчислюють за результатами обробки декількох послідовних серій періодичних вимірювань його параметрів з використанням статистичної обробки і фільтрації результатів обчислення в різних режимах експлуатації. Фактичні зведені характеристики ВЦН визначають у функції як від виду та величини узагальненого нормованого дефекту, так і від величини зведеної об'ємної продуктивності "шляхом лінійної інтерполяції розрахункових характеристик обраного виду узагальненого нормованого дефекту політропного коефіe та hp (Q , D ) полk, c пр k,c приведеного ступеня стиску ep ( Q , D )". прk,c пр k, c Визначення фактичних характеристик ВЦН за його розрахунковими є уніфікованим способом для всіх експлуатованих типів парку ВЦН газотранспортної системи, для яких вплив дефектів проточних частин на положення фактичних характеристик відносно паспортних може істотно відрізнятися. цієнта корисної дії 13 83696 14 На кресленнях наведені: Фіг.1 - система, що реалізує запропонований спосіб (приклад); Фіг.2 - схема послідовності дій запропонованого способу; Фіг.3 - модуль обробки серії векторів вимірювань Xk, i, l . Система, що реалізує пропонований спосіб, містить модуль узгодження паспортних характеристик ВЦН на основі використовуваної газодинамичної моделі стиску газу 1, модуль визначення розрахункових характеристик k-го ВЦН 2, датчики агрегатної автоматики k-го ВЦН 3, блок обробки вхідної інформації і формування вектора вимірювань X k, i,l 4, блок визначення поточних значень значень Qпр k,i, l, j , Dk,i,l, j 5, блок розрахунку прогнозованих льшення дефектів ìV Dфk,l, jü j = 1 ND та обчис, í ý î þ лення поточних прогнозованих значень узагальне { } значень Dr k,i,l, j 6, блок визначення найбільш ймовірного виду узагальненого нормованого дефекту k-го ВЦН D k,c 7, блок розрахунку фактичних зведених характеристик k-го ВЦН 8, блок розрахунку коефіцієнтів стискальності та поправок до теплоємності 9, блок візуалізації та архівування 10 (Фіг.1). Модуль узгодження паспортних характеристик ВЦН на основі використовуваної газодинамичної моделі стиску газу 1 призначений для уточнення положення однієї з трьох паспортних характеристик за двома іншими характеристиками на основі модифікованого рівняння стану Бендикта-ВеббаРобина та визначення показників процесу стиску газу за допомогою політропного методу Шульца. Модуль визначення розрахункових характеристик k-го ВЦН 2 призначений для визначення розрахункових зведеної характеристики ступеня стиску та політропного ККД функції від зведеної об'ємної продуктивності та величин і видів узагальнених нормованих дефектів. Датчики агрегатної автоматики k-го ВЦН 3 забезпечують періодичні вимірювання значення температури газу на вході і ви ході нагнітача Tвх k,i,l, Твих k,i,l, тиску газу на вході і виході нагнітача Рвх k,i,l, Рвих k,i,l частоти обертання валу нагнітача nk,i,l та густини газу r k,i,l. Блок обробки вхідної інформації і формування вектора вимірювань Xk,i,l 4 призначений для зчитування, аналізу вірогідності, відновлення недостовірної, фільтрації інформації датчиків агрегатної автоматики та для формування вектора вимірювань X k, i,l { } Xk, i, l = Tвх k, i, l, Tвих k, i, l, Рвх k, i,l , Р вих k, i, l, n k, i, l, r k, i, l . Блок визначення поточних значень Qпр k,i, l, j , Dk,i,l, j 5 призначений для розрахунку зведеного ступеня стиску e пр ( X k,i,l ) та політропного коефіцієнта корисної дії hпол( Xk, i, l ) за поточним вектором вимірювання Xk, i, l , знаходження поточних значень Qпр k,i, l, j та Dk,i, l, j на основі розрахункових характеристик за знайденим значенням e пр ( X k,i,l ) та hпол( Xk, i, l ) шляхом рішення системи рівнянь та формування множини зведеної {Q пр k, i,l, j }= ,i об'ємної продуктивності та 1 NI,= 1 ND, = 1 NK,= 1 NL, , j , k , l , множини значень узагальненого нормованого де { } фекту D k, i,l, j , i = 1 NI, j = 1, ND, k = 1, NK, l = 1, NL . , Блок розрахунку прогнозованих значень Dr k,i,l, j 6 призначений для визначення значень узагальнених нормованих дефектів D k= 1, ND та значень швидкості їхнього j ,0,k, j { } { } збільшення V D k, l, j j = 1, ND на момент початку l-ої серії вимірювань за результатами розрахунку значень цих дефектів на попередній серії вимірювань, обчислення фільтрованих значень швидкості збі { r } них нормованих дефектів D k , i, l . Блок визначення найбільш ймовірного виду узагальненого нормованого дефекту k-го ВЦН D k,c 7 призначений для вибору найбільш ймовірного узагальненого нормованого дефекту Dk,c, за обчисленими множинами значень та D k, i, l, j { } ì Dr ü í k,i,l, j ý на основі критерію мінімуму величини î þ середньоквадратичного відхилення узагальненого нормованого дефекту D k,i,l, j від його прогнозованого значення Dr k,i,l, j . Блок розрахунку фактичних зведених характеристик k-го ВЦН 8 призначений для їх визначення на основі вибору найбільш ймовірного виду узагальненого нормованого дефекту Dk,c та розрахункових характеристик ep (Qпр , Dkc ), hp (Qпр , Dkc ) . пол прk Блок розрахунку коефіцієнтів стискальності та поправок до теплоємності 9 призначений для розрахунку коефіцієнтів стискальності, коефіцієнтів ізобаричної стискальності, поправок до теплоємності на основі модифікованого рівняння стану Бенедикта-Вебба-Рабіна. Блок візуалізації та архівуванния 10 призначений для візуалізації та архівування поточної вхідної інформації та результатів розрахунку фактичних зведених характеристик ВЦН. У відповідності зі схемою послідовності дій (Фіг.2) запропонований спосіб реалізується таким чином. Дії по знаходженню фактичних характеристик парку ВЦН газотранспортної системи в умовах компресорних станцій розподіляються на дії, які виконуються одноразово та періодичні дії. За кожним типом ВЦН одноразово виконуються наступні дії: - узгоджуються паспортні характеристики ВЦН на основі використовуваної газодинамічної моделі стиску газу; - визначаються розрахункові характеристики ВЦН; 15 83696 Необхідність узгодження паспортних характеристик, реалізована в модулі 1, диктується тим, що усі використовувані характеристики ВЦН в способі, що заявляється: паспортні, розрахункові, фактичні повинні ґрунтува тися на єдиному рівнянні стану природного газу і на тому ж самому методі визначення показників процесу стиску газу в нагнітачі. Без узгодження паспортних характеристик ВЦН зростає методична помилка визначення фактичних зведених характеристик. Узгодження паспортних характеристик полягає в розрахунку однієї з трьох паспортних характеристик за двома іншими на основі модифікованого рівняння стану Бенедикта-Вебба-Рабіна та політропного методу Шульца визначення показників процесу стиску газу в нагнітачі: - витратно-напірної характеристики за характеристиками політропного ККД та зведеної відносної внутрішньої потужності; - характеристики політропного ККД за характеристиками витратно-напірної та зведеної відносної потужності; - характеристики зведеної відносної потужності за характеристиками витратно-напірної та політропного ККД; Розрахункові характеристики ВЦН визначають на основі натурних випробовувань та теоретичних досліджень зміни характеристик нагнітача для різних видів та величин дефектів проточної частини. За отриманими даними вказаних випробовувань та досліджень у модулі 2 визначають множини розрахункових характеристик політропного коефіцієнта корисної дії hp пол та зведеного ступеня стиску eр функції від зведеної об'ємної продукпр тивності Qпр та величини та виду узагальненого нормованого дефекту D kj, ìhp ü , , í полkj (D kj, Qпр)ý, k = 1 NK, j = 1 ND; î þ ìep (D , Q )ü, k = 1 NK, j = 1 ND , , , í пр kj kj пр ý î þ де NK - кількість газоперекачувальних агрегатів, ND - кількість видів узагальнених нормованих дефектів. У якості узагальненого нормованого дефекту приймають зважену суму конкретних видів дефектів проточної частини відцентрового нагнітача, що однозначно визначає відхилення цих характеристик від відповідних паспортних для різних видів і величин дефектів проточної частини нагнітача (зазор в ущільненні покриваючого диска, підріз лопаток робочого колеса та лопаткового дифузора та ін.). Також знаходять об'ємний номінальний показник політропи функції від об'ємної зведеної продуктивності ml ' (Q пр ). н Виконувані періодичні дії у свою чергу підрозділяються на обчислення, які передбачають обробку кожного вимірювання в серії з Nl вимірювань та на дії з результатами обчислень по закінченню кожної серії. По завершенню обробки NL-ої серії вимірювань визначаються фактичні зведені характеристики ВЦН (Фіг.2). 16 У модулі обробки серії векторів вимірювань Xk, i, l (Фіг.3) періодично за кожним і-им вимірюванням виконують наступні дії. За допомогою блоку обробки вхідної інформації і формування вектора вимірювань Xk,i, l 4 з датчиків агрегатної автоматики k-ого ВЦН періодично зчитують наступну інформацію: температура газу на вході та ви ході нагнітача Tвх k,i,l, Твих k,i,l, тиск газу на вході і виході нагнітача Рвх k,i,l, Рвих k,i,l частоти обертання вала нагнітача nk,i,l та густина газу r k,i,l. Цю ін формацію аналізують на достовірність, відновлюють недостовірна та фільтрують з відповідних датчиків агрегатної автоматики. Формують вектор вимірювань { } = X k,i, l Xk, i, l {Tвх k,i,l, Tвих k,i,l, Рвх k,i,l, Рвих k,i,l, n k,i,l, r k,i,l} . Далі в блоці розрахунку коефіцієнтів стискальності та поправок до теплоємності 9 розраховують коефіцієнти стискальності на вході та виході нагнітача Zвхk,і,l, Zвихk,i,l середній коефіцієнт стискальності Zсерk,i,l, середній коефіцієнт ізобаричної стискальності Vcер , середнє значення теплоємності газу на основі модифікованого рівняння стану Cp сер Бенедикта-Вебба-Рабіна. З урахуванням розрахованих значень коефіцієнтів стискальності та теплоємності обчислюють політропний напір Hn ( Xk, i,l ) та об'ємний і температурний показники політропи m V ( X k, i,l ) , m T ( Xk,i,l ) Hnпр 2 æ nн ö ç ÷ , [ Дж / кг]; = Hn ( X k, i, l ) × ç nk, i, l ÷ è ø Hn ( Xk,i,l ) = æ Zв х k ,i,l × R × Tв х k,i,l ç æ Pв их k,i,l çç mn ( Xk,i,l) ç ç Pв х k ,i,l è è m (X ) ö ö n k, i,l ÷ ÷ - 1÷, [ Дж / кг ]; ÷ ÷ ø ø æ Zвх k,i,l × Tвих k, i, l ö æ Pвих k,i, l ö ÷ / lgç ÷; mn ( Xk, i, l) = lgç ç Zвх k,i, l × Tвх k,i,l ÷ ç Pвх k,i,l ÷ è ø è ø æ Tвих k, i, l ö æ Pвих k, i, l ö ÷ / lgç ÷; mT ( Xk, i,l ) = lgç ç Tвх k,i,l ÷ ç Pвх k,i, l ÷ è ø è ø 1,204 × Rпов , [кгс × м / кг × К ] . rk, i, l Далі для кожного виду узагальненого нормованого дефекту D k,j виконують наступні дії. Розраховують значення об'ємної зведеної продуктивності Q S для обчислення об'ємного пр номінального показника політропи m V (Q пр ) як R= н середнє значення мінімального Qпр m in та максимального значень Qпр m ax зведеної об'ємної продуктивності. Шляхом лінійної інтерполяції розраховують об'ємний номінальний показник політропи m V (Q пр ) для Q пр = Q S . пр н У блоці визначення поточних значень Qпр k,I,l,j, Dk,I,l,j, 5 розраховують значення зведеного ступеня стиску e пр ( X k,i,l ) та політропного коефіцієнта 17 83696 корисної дії hпол( Xk, i, l ) відповідно до політропного методу Шульца для потокового вектора вимірювань X k, i,l за формулами: 1 æ m V н (Qпр k ,i,l ) × H n пр ö m V (Qпр k ,i ,l ) e пр (X k,i,l ) = ç + 1÷ н ; ç Zпр × R пр × Tвх пр ÷ è ø æ ö Zсер i × R ç ÷ hпол(X k,i,l ) = ç ÷ ç Cp - Zcее × R × Vcее і × mT - 1(X k,i,l ) ÷ è сер і ø -1 × 1 m T (X k,i,l ) . Далі знаходять Dk,I,l,j, Q пр k,і,l,j з рішення наступної системи рівнянь: ì e p ï пр k ( Xk, i,l ) = e пр k, j (D k, i,l, j , Qпр i, j ) . í p ïhпол k ( Xk, i,l ) = h пол k, j (Dk,i,l, j, Qпр i, j ) î За умови коли абсолютна величина різниці між Q S та Qпр i, j більше заданої константи Qс пр значенню Q S надають значення Qпр i, j та попр вторюють дії, починаючи з розрахунку об'ємного номінального показника політропи. В іншому випадку визначають множину значень зведеної об'ємної продуктивності { } = 1 NI, l 1 NL j 1 ND та мно= , Q= , , пр k, i, l, j , k 1 NK, i , жину значень узагальненого нормованого дефекту {D k,i,l, j}, k = 1, NK, i = 1, NI, l = 1, NLj = 1, ND . Вище вказані дії повторюють для всіх видів узагальненого нормованого дефекту {D k,= , i,l, j }k 1 NK, i = = 1, ND та для кож= 1 NI, l 1 NL j , , , ного вимірювання з серії i = 1 NI . , По закінченні (l-I)-ої серії вимірювань у блоці 6 для кожного j-го виду узагальненого нормованого дефекту обчислюють розрахункове значення його величини Dk,0,l,j та швидкості його збільшення VDk,l,j на момент початку наступної l-ої серії вимірювань, виходячи з рішення наступної системи рівнянь ì ¶Fk,l -I, j =0 ï ï ¶Dk,0,l, j , í ¶F ï k,l -I, j = 0 ï ¶VD k, l-I, j î де VD k,l, j = b1 × VD фk,l, j + b 0 × VD фk,l - l, j . З метою зменшення впливу похибок датчиків агрегатної автоматики на розрахунок прогнозованих значень узагальнених нормованих дефектів ìDr ü í k,i,l, j ý обчислення швидкості фільтрують за доî þ помогою фільтра першого порядку NI Fk,l = å (D k,i,l - 1 j - D r )2 , де b1, b0 - кое-1 j , , k,i,l -1, j = 1 i фіцієнти фільтра. Прогнозовані значення Dr розраховують k,i,l, j виходячи з пропорційного збільшення узагальненого нормованого дефекту у функції від часу його 18 наробітку tнар, вимірюваного від початку серії, за допомогою виразу Dr ,i,l, j = Dk,0,l, j + t нар × VD ф . k k,l, j По завершенню обчислення значень узагальнених нормованих дефектів та їхніх прогнозованих значень для всіх NL серій вимірювань у блоці визначення найбільш імовірного виду узагальненого нормованого дефекту k-го ВЦН Dk,c 7 вибирають найбільш імовірний узагальнений нормований дефект Dk,c , який однозначно характеризує положення фактичних характеристик відцентрового нагнітача відносно відповідних паспортних, за { } обчисленою множиною значень Dk, i,l, j , k = 1 NK , , i = 1 NI , l = 1 NL , j = 1,ND на основі критерію міні, , муму величини середньоквадратичного відхилення узагальненого нормованого дефекту від його прогнозованого значення Dr : k,i,l, j NL NI NL NI l=1i =1 l =1i =1 å å (Dk, i, l, c ( Xk,i,l ) - Dr ,i, l, c )2 < å å (Dk, i,l, j ( Xk, i, l ) - Dr ,i,l, j)2 k k для j=1,...,c-1,з+1,...ND... Фактичні зведені характеристики k-го ВЦН f eпр k (Qпр), hf k (Qпр) визначають шляхом ліпол нійної інтерполяції розрахункових характеристик обраного виду узагальненого нормованого дефекту с політропного коефіцієнта корисної дії hр k ( Qпр)= (Qпр, Dkc ) та зведеного ступеня пол p стиску e f (Q ) = h ( Q пр , D kc ) для значенпр k пр пр k , c ня Dkc, що дорівнює Dk,0,N L,c. Для визначення фактичних зведеної витратно-напірної e f (Qпр), пр k політропного ККД hр та зведеної внутполk (Qпр) рішньої відносної потужності f [Nu / r]прk (Q пр ) характеристик у вигляді апроксимуючих поліномів необхідного ступеня. Діапазон від мінімального значення зведеної об'ємної продуктивності Qпр min k до максимального Qпр max k розбивають на рівні відрізки та формують множину значень зведеної { } об'ємної продуктивності Q прt , де t = 1 NT та NT , кількість елементів масиву. Для сформованої множини значень Q прt за характеристиками { } f eпр k (Qпр) та hk k (Qпр) визначають відповідно пол масиви значень {eпр k, t } {hполk, t } t = 1,NT на ос нові яких розраховують апроксимуючі поліноми відповідно фактичних витратно-напірної та політропного ККД характеристик. Апроксимуючий поліном фактичної характеристики зведеної внутрішньої відносної потужності [Nu / r]f (Q пр ) прk розраховують за масивами значень e пр k, t та {hполk, t } t = { } 1 NT на основі використовуваної в , 19 83696 способі, що заявляється, газодинамічної моделі процесу стискання газу. Джерела інформації: 1. Галиуллин З.Т., Леонтьев Е.В. Интенсификация магистрального транспорта газа. -М.: Недра. 1991. - 272с. 2. Определение газодинамических характеристик нагнетателей на КС. А.М. Проскуряков, А.И. 20 Черников, В.И. Лысюк, М.Б. Письман. Газовая промышленность, № 5, 2000г, с.48-50. 3. Поршаков Б.П., Лопатин А.С., Назарьина A.M., Рябченко А.С. Повышение эффективности эксплуатации энергопривода компрессорных станций. - М.: Недра, 1992. - 207с.: ил. 21 83696 22 23 83696 24 25 Комп’ютерна в ерстка Т. Чепелев а 83696 Підписне 26 Тираж 26 прим. Міністерство осв іт и і науки України Держав ний департамент інтелектуальної в ласності, вул. Урицького, 45, м. Київ , МСП, 03680, Україна ДП “Український інститут промислов ої в ласності”, вул. Глазунова, 1, м. Київ – 42, 01601

Дивитися

Додаткова інформація

Назва патенту англійською

Method for determination of real characteristics of stock of centrifugal superchargers of gas transportation system

Автори англійською

Khymko Myroslav Petrovych, Pryschepo Oleksandr Oleksiiovych, Ped'ko Borys Ivanovych, Izmalkov Borys Ivanovych, Ilchenko Borys Samuilovych, Kolodiazhnyi Valerii Vasyliovych, Kucheruk Mykola Vasyliovych, Nalisnyi Mykola Borysovych, Shurunova Natalia Shyoivna

Назва патенту російською

Способ определения фактических характеристик парка центробежных нагнетателей газотранспортной системы

Автори російською

Химко Мирослав Петрович, Прищепо Александр Алексеевич, Педько Борис Иванович, Измалков Борис Иванович, Ильченко Борис Самуилович, Колодяжный Валерий Васильевич, Кучерук Николай Васильевич, Налисный Николай Борисович, Шурунова Наталья Шиоивна

МПК / Мітки

МПК: F04D 27/00, G01F 1/34, G01M 5/00, F04D 17/00

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

Код посилання

<a href="https://ua.patents.su/13-83696-sposib-viznachennya-faktichnikh-kharakteristik-parku-vidcentrovikh-nagnitachiv-gazotransportno-sistemi.html" target="_blank" rel="follow" title="База патентів України">Спосіб визначення фактичних характеристик парку відцентрових нагнітачів газотранспортної системи</a>

Подібні патенти