Analysis of Electrical Safety Conditions Taking into Account Soil Conductivity Determined on the Basis of Fuzzy Logic
Manusov V.Z. 1, Zaytseva N.M. 2
^Novosibirsk State Technical University, Novosibirsk, Russian Federation, 2Innovative Eurasian University, Pavlodar, Republic of Kazakhstan
Abstract. The goal of this work is to prove a possibility of determining soil parameters that influence its conductivity being the basis of grounding, step voltage and touch voltage calculation. This in its turn increases the safety level of electric equipment operation. The article is devoted to development of new, no conventional models of soil conductivity using the theory of fuzzy sets and fuzzy logic. The description of the solution includes the following sections: fuzzy models of specific electrical resistance of different soil types depending on their salinity and moisture content at positive temperatures; dependences between p and negative temperature values for different types of soil and different moisture contents; fuzzy climatic soil models for a specific location in order to determine soil p deep into the earth throughout the year. Uncertainty of soil properties is a reason why calculations of grounding devices are very liberal with plenty to spare, and calculations of touch voltages and step voltages in the period of active snow melting or abundant rains give erroneous results. Since soil parameters are described by linguistic variables: humidity, salinity, and temperature etc., soil layers should be considered as overlapping fuzzy sets. Membership functions of the main components of the "soil conductivity" concept were formed on this basis. The new scientific result we obtained makes it possible to formalize natural language of the experts in the field of climatic factors influence on soil resistivity in order to represent this information with the help of computer algorithms for optimum design of energy supply security and step voltage and touch voltage calculation. Keywords: electrical safety specific, soil resistivity, modeling, fuzzy logic, membership functions, multi-parameter dependence, touch voltage, step-by-step pressure.
Analiza conditiilor de securitate electrica folosind logica fuzzy pentru a dezvolta modele
de conductivitate a solului
Manusov V.Z. 1, Zaytseva N.M. 2 1Universitatea Tehnica de Stat Novosibirsky, Novosibirsk, Federatia Rusa 2Universitatea Inovatoare din Eurasia, Pavlodar, Kazahstan Rezumat: Scopul lucrarii consta in demonstrarea posibilitatii de precizare a parametrilor de sol care afecteaza conductivitatea electrica, care este baza de calcul alegarii la pamant, a tensiunii de pas §i tensiunile de atingere. Aceasta, la randul sau, conduce la cre§terea nivelului de siguranta in functionarea echipamentelor electrice. Articolul este dedicat dezvoltarii de noi modele inovatoare de conductivitate electrica a solului, folosind teoria multimilor fuzzy §i a logicii fuzzy. Solutia problemei rezolvate include urmatoarele componente: modele fuzzy a rezistentii electrice specifice a solului p de diferite tipuri in functie de salinitatea solului §i de umiditate la temperaturi pozitive;dependenta de marimea p pentru valori negative ale temperaturii pentru diferite tipuri de sol cu umiditate variabila; modele climatice fuzzy pentru zone concrete utilizate pentru determinarea marimii p in functie de adancimea pe parcursul anului. in conditii reale de valori de temperatura §i umiditate ale straturilor de sol variaza in functie de adancimea aparitiei acestora, nivelul apei subterane, temperatura la suprafata, precipitatii §i alti factori climatici. Ace§ti factori determina o schimbare in sol a marimii p intr-o gama foarte larga, de la unitati pana la cateva zeci de mii de Ohm*m. Pe aceasta baza, am format functiile de atributie a principalelor componente pentru notiunea „conductanta electrica a solului". Rezultatul ¡jtiintific nou obtinut permite formalizarea limbajului natural a expertilor privind domeniile de influenta a factorilor climatici asupra conductantii specifice a solului pentru prezentarea acestei informatii cu ajutorul algoritmilor de calculator in cazul proiectarii optimale a obiectelor de siguranta electrica, tensiunii de pas a tensiunilor de atingere. Cuvinte-cheie: siguranta electrica, rezistenta specifica a solului, modelare, functii logice fuzzy,functii de membru, dependenta multiparametrica, tensiune de pas de tensiune atingere.
Анализ условий электробезопасности с учетом удельной электропроводности грунта на основе
нечеткой логики Манусов В.З.1, Зайцева Н.М.2
Новосибирский государственный технический университет, Новосибирск, Российская Федерация ^Инновационный Евразийский университет, Павлодар, Казахстан Аннотация. Цель работы - доказательство возможности уточнения параметров грунта, влияющих на его удельную электропроводность, которая является основой расчета заземлений, шагового напряжения и напряжения прикосновения. Это в свою очередь повышает уровень безопасности при эксплуатации электрооборудования. Статья посвящена разработке новых, нетрадиционных моделей электропроводности грунта с использованием теории нечетких множеств и нечеткой логики. Решение задачи представлено следующими разделами: нечеткие модели удельного электрического сопротивления р различного вида грунта в зависимости от засоленности и влажности при положительной температуре; зависимости р от отрицательных значений температуры для различных видов грунта и разной увлажненности; нечеткие климатические модели грунта конкретной местности для определения р грунта вглубь земли в течение года. В реальных условиях значения температуры и влажности слоев грунта меняются в зависимости от глубины их залегания, уровня грунтовых вод, температуры на поверхности земли, количества осадков и других климатических факторов. Эти факторы вызывают изменение р грунта в очень больших пределах: от единиц до десятков тысяч Ом^м. Такая неопределенность приводит к тому, что расчет устройств заземления выполняется с большим запасом, а расчет напряжения прикосновения и шагового напряжения в период активного снеготаянья или обильных дождей и вовсе приводит к ошибочным результатам. Поскольку параметры грунта описываются лингвистическими переменными: влажность, соленость, температура и т.д., то слои грунта следует рассматривать как нечеткие пересекающиеся множества. На этой основе сформированы функции принадлежности основных компонентов понятия «электропроводность грунта». Полученный новый научный результат позволяет формализовать естественный язык экспертов в области влияния климатических факторов на удельное сопротивление грунта для представления этой информации с помощью компьютерных алгоритмов при оптимальном конструировании объектов энергобезопасности и расчете шагового напряжения и напряжения прикосновения.
Ключевые слова: электробезопасность, удельное электрическое сопротивление грунта, моделирование, нечеткая логика, функции принадлежности, многопараметрическая зависимость, напряжение прикосновения, шаговое напряжение.
Введение
Электротравматизм по сравнению с другими видами производственного травматизма составляет небольшой процент, однако по числу травм с тяжелым и смертельным исходом занимает одно из первых мест. Наибольшее число электротравм (60-70%) происходит при работе с электрооборудованием и электроустановками напряжением до 1000В.
К факторам, увеличивающим
безопасность при работе с электрооборудованием, безусловно,
относятся средства, предупреждающие поражение электрическим током. Так, например, устройства заземления, предназначенные для обеспечения принятых режимов работы
энергоустановок, имеют функцию защиты от поражения электрическим током путем отвода в землю потенциала, появившегося на корпусе электрооборудования в случае его аварийной работы. Еще одним средством, увеличивающим уровень
безопасности, являются различные системы предупреждения о наличии опасности поражения электрическим током. Современные интеллектуальные комплексы
автоматизированных систем управления могут иметь такую функцию и при возникновении аварии оповещать персонал, работающий в непосредственной близости к опасным зонам.
Однако при решении такой задачи возникает ряд сложностей, который обусловлен неопределенностью значения удельного электрического сопротивления земли в конкретный момент времени, влияющего на величину отводимого потенциала.
Так все формулы расчета заземлений и определения потенциала при стекание тока в землю содержат р — удельное электрическое сопротивление грунта, которое зависит от ряда параметров и меняется в очень широких пределах[1-10].
Существующие методы и рекомендации не позволяют учесть сезонное изменение удельного электрического сопротивления, например, в сезон снеготаянья или в период обильных дождей, а также в период засухи или сильных морозов, что ведет к ошибочному
расчету напряжений прикосновения и шага и может способствовать травмированию обслуживающего персонала.
I. Постановка задачи
В последнее время для решения задач, содержащих различного вида
неопределенности, стали применять методы нечеткой логики [11-22]. Данная работа выполнена на этой же основе. В работе предпринята попытка получения моделей, способных решить описанную выше проблему.
Удельное электрическое сопротивление грунта колеблется в широких пределах - от единиц и десятков до тысяч Ом на метр, так как зависит от множества факторов, основные из которых: вид грунта, его влажность, температура и плотность. Грунты в природе залегают слоями и не являются исключительно «чистыми», то есть они содержат примеси различных видов.
Все это вносит элемент неопределенности при расчете параметров любого вида грунта. Поэтому в качестве основного математического аппарата моделирования удельного сопротивления грунта был использован аппарат нечеткой логики[1] и регрессионного моделирования.
II. Методы и алгоритмы решения
Разработка модели велась следующими этапами:
1) поиск зависимостей удельного электрического сопротивления различного вида грунта от влажности при температуре 20оС;
2) поиск зависимостей р от отрицательных значений температуры для различного вида грунта и разной увлажненности;
3) построение климатической модели грунта конкретной местности для определения р вглубь земли в течение года;
4) разработка методики расчета шагового напряжения на основе моделей грунта.
Результаты экспериментов по трем видам грунта показали, что на интервале влажности от 0 до 2% для песка, супеси от 0 до 5% и глин от 0 и 6% характер изменения
p(v) резко отличается от характера изменения р на интервале от 2, 5, 6% для песка, супеси и глин соответственно и до насыщения грунтов влагой. Поэтому построение эмпирических зависимостей p(v) было решено выполнять на основе регрессионного моделирования для каждого такого интервала в отдельности, объединяя в одно математическое выражение с помощью инструмента нечеткой логики: системы размытых правил.
В данном случае нечетко определяется влажность, которая характеризуется такими понятиями, как: - "сухо" - примерно от 0 % до 3-4%, "умеренно влажно" - примерно от 4 % до 12%, "влажно" - примерно от 12 % до насыщения грунта влагой.
Для определения функций принадлежности использовался прямой метод, в котором для каждого xeE (x - влажность грунта) задается значение /uA(x)(функция принадлежности). Это правомерно, поскольку метод применяется для измеримых понятий, таких как давление, температура и т.д., а влажность, несомненно, относится к этим понятиям.
Суть модели заключается в следующем. Для каждого вида грунта разрабатывается система размытых правил. При этом она может быть записана следующим образом [11,12]
IF (xx е Ли) И (x 2 е Л 2 J И ... И (xk е Aki ) THEN y = "л(x),i = 1,. ,N
(1)
где Ajj -нечеткое подмножество, то есть нечеткий интервал, для переменной xj с функцией принадлежности ц Ajj (x), N - число
правил (число интервалов); y = "i(x) -функция, определяющая локальное решение модели от набора: x = (xx,...,xk).
Для одномерной зависимости (речь идет о параметре влажность грунта) система размытых правил (1) приобретает вид IF x е Ai THEN y = л ( x ), i = 1,... ,N (2) где Ai имеют функцию принадлежности цA, (x) . Для описания нечеткой
принадлежности к интервалам приняты экспоненциальные функции:
Цл (x) = e(-(x-d)2) и
РКОБЬЕМЕЬЕ ЕКЕКвБТЮП КЕОЮМТЛЬЕ 2(34)2017 ЕЬЕСТКОЕКЕКОЕТЮЛ
ц м (х) = 1" е'
(-а(х-й)2 )
(3)
где коэффициенты с и й определяют степень нечеткости и положение на числовой оси параметра влажность.
ц,+1 = 1 -ц ■■•
Регрессионную кривую строили с использованием показательной функции:
Л <(х ) = Ь о • Ь 1х
(4)
Для данной модели функции принадлежности имеют вид, показанный на рисунке 1, где V -влажность.
Рис.1 Семейство функций принадлежности для нечетких интервалов влажности.
Значения интервала влажности ("сухо") и v2 ("влажно") определены приблизительно, первый интервал -0 < V < 2 %, второй - V > 6 %, область переключения 2 < V < 6 %, функции принадлежности, причем в области, где соседние ц 1 и ц 2 имеют ненулевые
значения, и всегда выполняется равенство:
2
^ ц ■ = 1. Ниже представлены выражения
1=1
(5), определяющие их значения. На рис.1 сплошной линией изображена функция принадлежности к нечеткому интервалу «сухо», а штриховой линией - функция принадлежности - «влажно».
ц1 =1
Ц2 =1
1,0 < v< 2
-0.8^-2)2 2 < < 6 е ,2 < V < 6
0, V > 6
'0,0 < V < 2
, -0.8^-2)2 2 < < 6 1 - е 4 ,2 < V < 6
1, V > 6
В результате получены уравнения удельного сопротивления для песка рП, супеси рС и
глины Ргл ,[23]:
р п = 6 • 0^ц/ V ) + 1 • 0Л v•ц 2 (V )
рс = 90 • 0Л^ц/V) + 0^3 • 0•8V • ц2(V) , (6)
п
с Р ГЛ
= 100 • 0^25v •ц1(V) + 3 • 0^ •ц2(V)
где V - влажность грунта в процентах.
Сравнение значений, полученных с помощью модели, и справочных данных показало, что модель хорошо описывает зависимость удельного электрического сопротивления от влажности.
Например, для песка сильно увлажненного грунтовыми водами, разброс значений удельного сопротивления составляет 10-60 Ом-м[24]. По расчетам с помощью (6) для такой же увлажненности получается 10 Ом-м. Для умерено увлажненного песка - 60-130 Омм, в то же время по модели - 70 Омм, для слегка влажного песка - 400-1500 Омм, по модели - 350 Омм при влажности 2%, т.е на интервале от 5% до насыщения совпадают с экспериментальными точками , а при уменьшении влаги от 5% до нуля отклоняются от них незначительно.
Известны градации засоленности грунтов: незасоленные (0,25 - 0,30% содержания солей), слабозасоленные (0,30 - 0,50%), среднезасоленные (0,50 - 1,0%), сильнозасоленные (1,0 - 2,0%), солончаки (2,0 - 4,0%). Для учета степени засоленности введен коэффициент, величина которого определяется с помощью нечеткой логики.
Рис. 2 Семейство функций принадлежности для нечетких интервалов засолённости.
Для определения величины этого коэффициента введено два нечетких интервала «слабозасоленные» и «засоленные» и определены функции принадлежности
где р20- удельное сопротивление при
температуре 20оС, 1 - температура грунта.
Для выявления закономерности изменения удельного электрического сопротивления при замерзании влаги в грунте был произведен ряд экспериментов.
Эксперименты показали, что при сухом грунте с уменьшением температуры р уменьшается по экспоненте, а при влажном - в районе от 00 C до -10 С происходит скачкообразное увеличение р, что объясняется кристаллическим строением исследуемых грунтов (пески, супеси).
Для расчета величины этого скачка для песка и супеси были выведены следующие формулы[23]:
н1 =
Ц 2 =
1, 0 < sol < 0.4
e"08(s°l-2)2 , 0.4 < sol < 2
0, sol > 2
0, 0 < sol < 0.4
1 _ e-0.8-(soi-2)2 , 0.4 < sol < 2
1, sol > 3
(7)
На рис.2 сплошной линией изображена функция принадлежности к нечеткому интервалу «слабозасоленные», а штриховой линией - функция принадлежности -«засоленные».
Значения коэффициентов засоленности вычисляются с помощью следующих функций:
K
= S "Л i( so1 ) • н i( so1 )
где " 1 (sol ) = 1 л 2 (sol ) = 0. 001 .
?
Для получения удельного электрического сопротивления грунта при конкретной влажности и засоленности необходимо значения, полученные с помощью выражения 6 умножить на Ksol.
Известна экспериментальная
зависимость расчета удельного
электрического сопротивления
влагосодержащих грунтов [24] от температуры
Р = Р2
-0.022(t-20 )
Р П = Р П ( - 0,03 •v2 + 0.86 •v - 1.9 ) р ^ = р С ( 0.024 •v2 - 0.022 • v + 0.2 )
(8)
где V- влажность, выраженная в долях, рПи рС - удельные сопротивления песка и супеси
при 0 0С , а рП и рС при -10С.
Расчет удельного электрического сопротивления песка и супеси при температуре ниже — 10 C может быть выполнен с помощью зависимостей:
рп =рП • 0.87(t+1), Рс =рС • 0.88'
(t+1)
(9)
Для глин и суглинков скачкообразного
0 р
увеличения р в районе 0 с не происходит ввиду аморфного строения данного вида грунта.
Поэтому для расчета удельного электрического сопротивления глин в области отрицательных температур предлагается следующая зависимость:
Р гл = Р гл • 0 88
(t+1)
(10)
где р^ - удельное сопротивление глины при 00 С.
Анализ значений р, полученных с помощью разработанных зависимостей, показывает, что они входят в диапазоны, указанные в [24].
Таким образом, полученная модель может быть применена для определения удельного электрического сопротивления любого типа
sol
i = 1
e
грунта от влажности и температуры с достаточной для практики точностью.
Климатическая модель параметров грунта. Для построения климатической модели грунта был выполнен анализ годового хода температуры на поверхности земли и для различных глубин по разным районам СНГ.
Данные были получены от метеостанций г.Павлодара и г.Новосибирска и взяты из литературы.
Для моделирования колебания t на заданной глубине, на основе вышесказанного и анализа данных метеостанций, была найдена следующая зависимость [25]:
/ = 1П - Ль • cos
2л 365
(g - 20 • к)
(11)
где Ак - амплитуда годового колебания t (°С) на к(м) от поверхности, вычисляемая
как:
А,
Лп • чк
а
- величина,
характеризующая убывание температуры вглубь грунта: д = ехр( 1п( Лпт / Ап )/ ) ,
/ к ПТ
АП - амплитуда годового колебания 1;(°С) температуры поверхности грунта;, ИПТ -глубина постоянной годовой температуры для средних широт И=15 метров; коэффициент «20» - учитывает запаздывание (в сутках) максимума (или минимума) годового колебания t на 1м глубины относительно t поверхности грунта (например, если максимум t поверхности наблюдается в июле, то на глубине в 5 м он будет через 100 суток); g - дни, 1;П.(температура слоя постоянной годовой температуры) для средних широт приблизительно рана 8 оС,И - глубина от поверхности земли в метрах; АПТ -амплитуда колебания t на глубинах с постоянной годовой температурой (приблизительно 0,1оС),
Л п =
(Т - т
— ' п тах п тт
Л
здесь Т
Птах
максимальная среднемесячная температура, ТПтт - минимальная среднемесячная температура за данный год.
Основной электрический проводник в грунте - влага. Именно она оказывает самое существенное влияние на сопротивление грунта. Поэтому при построении модели удельного сопротивления грунта в
реальных условиях[26] необходимо включить определение влажности его слоев в данном конкретном месте в определенный период времени.
Моделирование зависимости V по глубине до 1м было выполнено с помощью нечеткой логики. Нечетко определялся интервал значения «глубина».
Для получения общего выражения моделирования для любого вида грунта использовалось два терма: первый - «верхний» с ФП ц1(к), второй - «нижний» с ц2 (к) . Соответственно этому, введено два интервала: первый - 0 < к < 20 см (выбрано 20см и 100см, т. к. метеостанции измеряют влажность грунта на этих глубинах), второй - к > 100см, область переключения между интервалами:
20 < к < 100см, где И - глубина грунта. Вид функций принадлежности ц1(к) и ц2 (к) представлен на рис.3. Для терма «верхний» на первом интервале ц1(к) = 1, на втором ц1 (к) = 0, на переходном интервале
ц 1 (к) = 1 -
к - 20 80
Если ц1 (И) убывающая функция, то цх(к) = 1 -ц2(к), при этом
ц2 (к Л =
к - 20
х - а
100 - 20 Ь - а
Рис. 3. Семейство ФП для определения влажности по глубинам.
ФП ц 1 (к ) и ц 2 (к ) имеют вид (12).
Определение функций Л\(к) и л2(к) выполнено на основе анализа данных метеостанций за период в 4 года, в результате которого были построены многофакторные модели влажности грунта на глубинах 0.2м и 1 м.
ц/к) =
Ц 2 (Ь) = <
1, 0 < к < 20
к — 20 1 — , 20 < к < 100 80
0, к >100
0, 0 < к < 20
к — 20
80
1, к > 100
20 < к < 100
(12)
В соответствии с теорией нечетких множеств выражение для определения влажности представляется так
2
г =1
(13)
Определение функций Ц1(к) и |2 (к) выполнено на основе анализа данных метеостанций за период в 4 года, в результате которого были построены многофакторные модели влажности грунта на глубинах 0.2м и 1м. С помощью метода наименьших квадратов была получена трехфакторная модель:
у20 = 8.06+0.03-ос — 0.20 • tnг — 0.81-Ч; у100 = —4. 10+0. 03-ос — 0. 007 • * +1. 3-и,
Для построения модели определения влажности в любое время года на глубине до 1м была введена лингвистическая переменная «месяц» с двумя термами «холодный» и «теплый». В первый временной интервал моделирования входят месяцы, в которых температура воздуха устойчиво держится ниже нуля градусов Цельсия (для рассматриваемой климатической зоны с ноября по февраль), во второй - с апреля по октябрь. Таким образом, можно считать, что функция принадлежности ц 1 (т ) = 1 и ц2(т ) = 0 с ноября по февраль для первого терма «холодный», а с апреля по октябрь - ц(т)=0 и ц2(т)=1 для второго терма «теплый», где т - номер месяца. Интервал переключения - март (месяц активного снеготаяния) и ноябрь (месяц, когда температура грунта приближается к отрицательным значениям). Причем, в месяц активного снеготаяния к количеству выпавших за данный месяц осадков необходимо добавить осадки, накопленные в виде снежного покрова за зимние месяцы. Функции принадлежности имеют вид (см. рис.4):
Ц 1(т) =
0,
( —3• (т — 4)2 )
,( —3• (т —102 ))
т < 3, т > 11
3 < т < 4
4 < т < 10
, 10 < т < 11
и
где Ос - количество осадков, мм; 1пг -температура на поверхности грунта, 0С, ив -
средняя скорость ветра, м/с.
Для оценки модели и определения процента разброса функции V относительно ее среднего значения в зависимости от количества выпавших осадков,
температуры на поверхности грунта и скорости ветра был вычислен коэффициент детерминированности Бух. Для глубины
20см Вх =09, для глубины 100см Бух=0.96,
что говорит об адекватности модели реальным данным.
Функции Г|/к) = у20 и Г[2(И) = Ую0 из выражения (14). Исходя из (14) зависимость определения влажности V для глубины до 1м имеет вид:
Ц 2 (т ) = 1 — Ц1 (т )
(16)
Рис.4. Семейство функций принадлежности для определения влажности по месяцам.
В результате для глубины 1м влажность в зависимости от ос, *пг,ив и времени года (по месяцам) вычисляется по формуле:
V = Х|1(т) • №(т) =11(т) • Ц1(т) +12(т) • Ц2(т),(17) 1=1
V к = V
цх(И) + v100 •Цг(к) (15) где 11(т) = Vkm=10; — влажность за октябрь,
20
1
рассчитанная по (14), так как для рассматриваемого региона (средние широты) в зимние месяцы (т=11,12,1,2) V не изменяется; л2(т) = vk - влажность за
теплые месяцы (т=4^10).
Влажность грунтов, расположенных ниже 1м от поверхности, практически не зависит от климатических условий, а определяется уровнем грунтовых вод. Для разработки модели определения влажности грунта была использована высота ИКП капиллярного поднятия (водоподъемная способность): для песка эта величина составляет 0,5м, для супеси - 1м, для суглинка - 1,5м, для глины - 2,5м, причем высота и скорость подъема зависят от структуры грунта. В песчаных грунтах вода поднимается невысоко, но достаточно быстро: за 1-5 дней, а в глинистых -медленно: в течение 5-14 дней. Для моделирования влажности введена нечеткая переменная, «близость грунтовых вод», и определено два терма: первый - «рядом» с функцией принадлежности ц 1 (кГВ ) и
второй «далеко» с ц2(кГВ) . Следовательно, будет два интервала: первый на промежутке 0 < кГВ < 0. 8ккп м; второй - кГВ > кш м, а
интервал
переключения
0.8ккп < кГВ < ккп, м. Значение 0.8ккп выбрано по той причине, что с этого расстояния от уровня воды грунт начинает «высыхать». Для терма «рядом» на первом интервале функция принадлежности ц 1 (кГВ ) =1, на втором ц1(кГВ) =0, на интервале переключения
к -0 8• к
ц1 (кГВ) = 1- 08 ккп
0.. 2 • кк
Для
терма
«далеко» этих же
ц2 (кГВ ) = 1 -ц/кд ) на интервалах. Функции принадлежности ц 1 (кГВ ) и ц2(кГВ ) имеют вид
ИОгв ) =
1, при 0 < ИГВ < 0.8ИКП ИГВ -0.8 • И,
1-
0.2 • ИК
^, 0.8Икп < ИГВ < ИКП
11кп
0, при игв > икп ц2(игв ) = ^ИОгв )
(18)
Формула для определения V запишется следующим образом
2
v = Е Л ¡( к ГВ ) • ц ■( к ГВ )
1 = 1
Зависимость л 1 (И ГВ ) построена по данным геологоразведочного центра в виде кривой, полученной при использовании метода наименьших квадратов, например для песка ц1(кТВ ) = 18 - 27 • кГВ , а Л2(кГВ ) =4.5 (4.5 -влажность естественного залегания песка в процентах). В результате модель определения V грунтов выглядит следующим образом [26]
^есок = (18 - 27 • кГВ ) •ц1(кГВ ) + 4 5 ^2 (кГВ ) ^супесь = (23 - 19 • кГВ ) •ц1(кГВ ) + 3 8 ^2 (ккГВ ) (19) ^глина (50 12. 8 кГВ ) •ц/кз ) +18 • ц2 (кГВ )
ГВ
= (45 - 25 • кГВ ) •ц/ка ) + 7.5 •ц (кГВ )
Здесь коэффициенты 18, 23, 50, 45 - это максимальная влагоёмкость песка, супеси, глины и суглинка соответственно. Адекватность модели подтверждается с помощью критерия Фишера, который равен 3.34 для сухого грунта и 3.59 для влажного, что выше табличного 3.26 (при уровне значимости ^=0.05).
Таким образом, зная количество слоев, их вид и глубину залегания грунтовых вод, с помощью модели (19) можно определить влажность грунта в любое время года на глубинах ниже 1м с достаточной для практики точностью. Погрешность моделирования влажности не превышает 20% для сухих грунтов и 10% для влажных.
Анализ электробезопасности с учетом нечетких моделей р грунта в зависимости от метеофакторов. На основании разработанных климатических моделей удельного
электрического сопротивления грунта была создана программа, с помощью которой определяется удельное электрическое сопротивление слоев грунта при указании их мощности (толщины) и вида. Такие данные определяют в обязательном порядке при изыскательских работах во время проектирования объектов энергобезопасности. Для учета климатического влияния на удельное
V
суглинок
электрическое сопротивление грунта требуется ввести усредненные помесячные климатические данные за последние 12 месяцев (температуру, количество осадков и силу ветра), а также уровень грунтовых
вод и мощность водоносного слоя, если грунтовые воды находятся в зоне установки устройств заземления.
р Ом м Удельное электрическое сопротивление грунта
300
200
100
о 1
1 2 3 4 5 6 7 8 9 10 11 12 13 14
Рис.5. Изменение удельного электрического сопротивления грунта по глубине залегания.
На рис. 5 представлен расчет удельного электрического сопротивления р четырех слоев грунта на глубинах от 0м до 14м. Верхний слой мощностью 3м - супесь, следующий слой мощностью 3м - песок, следующий слой -водоносный песок (толщина водоносного слоя 2м), ниже идет слой глины (толщина слоя-6м). На глубине 5.5м значение удельного электрического сопротивления песка начинает резко падать, так как далее идет водоносный слой, р которого близко к нулю, точнее 0.50м-м.
На рис.5 приведены кривые, характеризующие удельное электрическое сопротивление слоев грунта в начале декабря (линия с маркерами) и в марте, в период активного снеготаянья (сплошная линия). Рисунок хорошо иллюстрирует падение р верхнего слоя земли практически до нуля в начале весны (сплошная линия) и возрастание р верхнего слоя земли в декабре, когда температура устойчиво держится ниже нуля (линия с маркерами).
Глубина промерзания грунта в данной климатической зоне на уровне 1.4м. Поэтому р ниже этого уровня резко падает до величины, соответствующей
Во время активного снеготаянья или обильного дождя удельное электрическое сопротивление верхнего слоя земли резко падает,
и потенциал, появившийся в случае аварии, может распространиться на значительное расстояние от места входа его в землю и привести к электротравматизму работающего персонала Современные континентальные осадки имеют удельное сопротивление 1-3 Ом-м[24]. положительным температурам.
В промышленных зонах удельное сопротивление талых вод и осадков меньше 1 Омм. Исходя из этих величин, был выполнен расчет распространения потенциала по поверхности земли от места его входа (см. рис.6).
Рис. 6. Изменение потенциала на поверхности земли в случае замыкания фазы на землю в период активного снеготаянья (кривая с крупными маркерами) и в период без обильной влаги.
Для расчета потенциала поверхности земли, для удаленной на расстояние х метров точки от заземлителя используется известная формула:
1 3 P , ф = —3— In
2 к l
vx
+ l2 + l
(20)
где p - удельное сопротивление земли (1 Омм при обильных дождях и активном снеготаянье), x - расстояние от места входа потенциала, l -величина слоя, имеющего сопротивление p.
Опасным напряжением считается 36в, которое может быть снижено при ряде условий. Таким образом, с использованием предлагаемых моделей и выражений 20 определяется опасное расстояние для работающего персонала.
Из учета, что длина шага 0.8м опасным расстоянием от заземления будет 2.5 м, о чем должно быть сообщено по средствам общей связи.
Опасным для прикосновения напряжением в данном случае будут обладать все проводящие электричество предметы на расстоянии в 9м от точки входа тока в землю. При обычном увлажнении опасного потенциала не возникает (рис.6, кривая с мелкими маркерами).
Выводы
1. Показано, что моделирование удельного электрического сопротивления p и других оказывающих влияние на p показателей грунта с использованием нечеткой логики позволяет определять значение этих величин в зависимости от вида грунта и климатических условий с погрешностью не более 20%.
2. На основе экспериментов получены функции принадлежности и нечеткие интервалы для таких показателей как: вид грунта, влажность, засоленность, а также изменение температуры и влажности слоев грунта под воздействием климатических факторов конкретной местности, что позволило построить модели p грунта.
3. Предложенные нечеткие модели позволяют определить параметры грунта в любое время года и на любой глубине с достаточной для практики точностью.
Литература (References)
[1] Gillies D. A. et al. Current North American assessment and refurbishment practices of substation grounding systems //IEEE Transactions on Power Delivery. - 2005. -V. 20. - №. 3. - pp. 1886-1889.
[2] Averbukh M. A., Zabusov V. V. Calculating the Parameters of Grounding Devices Using Fuzzy-Sets and
Fuzzy Logic Theory.// Electrical Technology Russia. 2009. №7- pp. 2-4.
[3] Guemes J. A., Hernando F. E. Method for calculating the ground resistance of grounding grids using FEM //IEEE Transactions on Power Delivery. - 2004. - V. 19. -№. 2. - P. 595-600.
[4] Nahman J., Paunovic I. Resistance to earth of earthing grids buried in multi-layer soil //Electrical Engineering. -2006. - V. 88. - №. 4. - pp. 281-287.
[5] Khan Y. et al. Efficient use of low resistivity material for grounding resistance reduction in high soil resistivity areas //TENCON 2010-2010 IEEE Region 10 Conference.
- IEEE, 2010. - pp. 620-624.
[6] Guemes-Alonso J. A. et al. A practical approach for determining the ground resistance of grounding grids //IEEE transactions on Power Delivery. - 2006. - V. 21. -№. 3. - pp. 1261-1266.
[7] He J. et al. Decreasing grounding resistance of substation by deep-ground-well method //IEEE transactions on power delivery. - 2005. -V. 20. - №. 2. -pp. 738-744.
[8] Nahman J., Paunovic I. Effects of the local soil nonuniformity upon performances of ground grids //IEEE Transactions on Power Delivery. - 2007. - V. 22. - №. 4.
- pp. 2180-2184.
[9] Trifunovic J., Kostic M. Analysis of influence of imperfect contact between grounding electrodes and surrounding soil on electrical properties of grounding loops //Electrical Engineering. - 2014. - V. 96. - №. 3. -pp. 255-265.
[10] Salam M. A. Grounding resistance measurement by grid electrode in Brunei Darussalam //International Journal of Energy Technology and Policy. - 2012. - V. 8. - №. 2.
- pp. 196-208
[11] Zadeh L. Fuzzy Sets // Information and Control. -1965. - № 8. - pp. 338-353.
[12] Zadeh L. A. Fuzzy sets as a basis for a theory of possibility //Fuzzy sets and systems. - 1999. - V. 100. -pp. 9-34.
[13] Mendel J. M. Uncertain rule-based fuzzy logic systems: introduction and new directions. - Upper Saddle River : Prentice Hall PTR, 2001. - pp. 131-184.
[14] Ross T. J. Fuzzy logic with engineering applications.
- John Wiley & Sons, 2009. 628 P
[15] Coupland S., John R. Geometric type-1 and type-2 fuzzy logic systems //IEEE Transactions on Fuzzy Systems. - 2007. - V. 15. - №. 1. - pp. 3-15.
[16] Greenfield S. et al. The collapsing method of defuzzification for discretised interval type-2 fuzzy sets //Information Sciences. - 2009. - V. 179. - №. 13. - pp. 2055-2069.
[17] Nie M., Tan W. W. Towards an efficient type-reduction method for interval type-2 fuzzy logic systems //Fuzzy Systems.IEEE World Congress on Computational Intelligence. IEEE International Conference on. - IEEE, 2008. - pp. 1425-1432.
2
x
[18] Naz M., Shabir M. Fuzzy Soft Sets and their Algebraic Structures //World Applied Sciences Journal. -2013. - V. 22. - pp. 45-61.
[19] Rutkovskaya D., Pilin'skii M., Rutkovskii L. Neyironnye sety, geneticheskie algoritmy i nechetkie sistemy. [Neural networks, genetic algorithms and indistinct systems.] Moscow. - Goryachaya liniya-Telekom, 2008. P 383.
[20] Kruglov V.V., Dli M.I.,Golunov R.Yu. Nechetkaya logika i iskusstvennye neyronnye sety. [Fuzzy logic and artificial neural networks]. Moscow, Izdatelstvo Fiziko-matematicheskoy literatury Publ., 2001, 224 P.
[21] Manusov V.Z., Ahyoev D.S. Technical diagnostics of electric equipment with the use of fuzzy logic models//Applied Mechanics and Materials. 2015. Vol. 792: Energy Systems, Materials and Designing in Mechanical Engineering.- pp.324-329.
[22] Andrzej Piegat. Fuzzy Modelling and Control. Phisica-Verlag. Heidelbeng. 2001. 798 P.
[23] Zaytseva, N.M., D.S. Zaytsev and M. Ya. Kletsel. Dependences of Soil Electrical Resistivity on Moisture Content and Temperature. // Electrical Technology Russia, 2008, Issue 9, pp 30-34.
[24] Karyakin R.N. Zazemlyayuchshie ustroystva elektroustanovok. Spravochnik. [Grounding Devices of Electroinstallations. Directory] Moscow, Energoservis Publ., 2002, 218 P.
[25] Mark Yakovlevich Kletsel, Natalia Mikhaelovna Zaytseva and Bibigul Beysembayevna Issabekova. Automation of Calculation of Resistances of Grounds// World Applied Sciences Journal, 2013, Vol. 25, Issue 6, pp.921-926.
[26] N.M.Zaytseva, B.B.Isabekova & M.Ya.Kletsel'. Determination of soil parameters to calculate soil resistivity.//Russian Electrical Engineering, 2015, Vol. 86, Issue 5, pp 275-281.
Сведения об авторах.
Манусов Вадим Зиновьевич. Новосибирский государственный технический университет, кафедра системы электроснабжения предприятий, д-р техн. наук, профессор. Основное направление исследований: применение методов искусственного интеллекта для планирования и оптимизации режимов электроэнергетических систем. E-mail: [email protected]
Зайцева Наталья Михайловна, Инновационный Евразийский университет, кафедра энергетики, металлургии и информационных технологий, кандидат техн. наук, доцент. Основное направление исследований: применение методов искусственного интеллекта для планирования и оптимизации функционирования электротехнических комплексов. E-mail: [email protected]