KR20070080913A - 순간 전압 강하 평가를 위한 취약 지역 결정 시스템, 방법,및 상기 방법을 실행시키기 위한 컴퓨터 판독 가능한프로그램을 기록한 매체 - Google Patents

순간 전압 강하 평가를 위한 취약 지역 결정 시스템, 방법,및 상기 방법을 실행시키기 위한 컴퓨터 판독 가능한프로그램을 기록한 매체 Download PDF

Info

Publication number
KR20070080913A
KR20070080913A KR1020060012419A KR20060012419A KR20070080913A KR 20070080913 A KR20070080913 A KR 20070080913A KR 1020060012419 A KR1020060012419 A KR 1020060012419A KR 20060012419 A KR20060012419 A KR 20060012419A KR 20070080913 A KR20070080913 A KR 20070080913A
Authority
KR
South Korea
Prior art keywords
voltage
voltage drop
point
line
accident
Prior art date
Application number
KR1020060012419A
Other languages
English (en)
Other versions
KR100782354B1 (ko
Inventor
장길수
박창현
Original Assignee
고려대학교 산학협력단
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 고려대학교 산학협력단 filed Critical 고려대학교 산학협력단
Priority to KR1020060012419A priority Critical patent/KR100782354B1/ko
Publication of KR20070080913A publication Critical patent/KR20070080913A/ko
Application granted granted Critical
Publication of KR100782354B1 publication Critical patent/KR100782354B1/ko

Links

Images

Classifications

    • E—FIXED CONSTRUCTIONS
    • E04—BUILDING
    • E04B—GENERAL BUILDING CONSTRUCTIONS; WALLS, e.g. PARTITIONS; ROOFS; FLOORS; CEILINGS; INSULATION OR OTHER PROTECTION OF BUILDINGS
    • E04B1/00—Constructions in general; Structures which are not restricted either to walls, e.g. partitions, or floors or ceilings or roofs
    • E04B1/62—Insulation or other protection; Elements or use of specified material therefor
    • E04B1/74—Heat, sound or noise insulation, absorption, or reflection; Other building methods affording favourable thermal or acoustical conditions, e.g. accumulating of heat within walls
    • E04B1/88—Insulating elements for both heat and sound
    • E04B1/90—Insulating elements for both heat and sound slab-shaped
    • E—FIXED CONSTRUCTIONS
    • E04—BUILDING
    • E04C—STRUCTURAL ELEMENTS; BUILDING MATERIALS
    • E04C2/00—Building elements of relatively thin form for the construction of parts of buildings, e.g. sheet materials, slabs, or panels
    • E04C2/02—Building elements of relatively thin form for the construction of parts of buildings, e.g. sheet materials, slabs, or panels characterised by specified materials
    • E04C2/10—Building elements of relatively thin form for the construction of parts of buildings, e.g. sheet materials, slabs, or panels characterised by specified materials of wood, fibres, chips, vegetable stems, or the like; of plastics; of foamed products
    • E04C2/16—Building elements of relatively thin form for the construction of parts of buildings, e.g. sheet materials, slabs, or panels characterised by specified materials of wood, fibres, chips, vegetable stems, or the like; of plastics; of foamed products of fibres, chips, vegetable stems, or the like
    • D—TEXTILES; PAPER
    • D21—PAPER-MAKING; PRODUCTION OF CELLULOSE
    • D21J—FIBREBOARD; MANUFACTURE OF ARTICLES FROM CELLULOSIC FIBROUS SUSPENSIONS OR FROM PAPIER-MACHE
    • D21J3/00—Manufacture of articles by pressing wet fibre pulp, or papier-mâché, between moulds
    • D21J3/12—Manufacture of articles by pressing wet fibre pulp, or papier-mâché, between moulds of sheets; of diaphragms
    • E—FIXED CONSTRUCTIONS
    • E04—BUILDING
    • E04B—GENERAL BUILDING CONSTRUCTIONS; WALLS, e.g. PARTITIONS; ROOFS; FLOORS; CEILINGS; INSULATION OR OTHER PROTECTION OF BUILDINGS
    • E04B1/00—Constructions in general; Structures which are not restricted either to walls, e.g. partitions, or floors or ceilings or roofs
    • E04B1/62—Insulation or other protection; Elements or use of specified material therefor
    • E04B1/74—Heat, sound or noise insulation, absorption, or reflection; Other building methods affording favourable thermal or acoustical conditions, e.g. accumulating of heat within walls
    • E04B2001/742—Use of special materials; Materials having special structures or shape
    • E04B2001/746—Recycled materials, e.g. made of used tires, bumpers or newspapers
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02W—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO WASTEWATER TREATMENT OR WASTE MANAGEMENT
    • Y02W30/00—Technologies for solid waste management
    • Y02W30/50—Reuse, recycling or recovery technologies
    • Y02W30/64—Paper recycling

Abstract

순간 전압 강하 평가를 위한 취약 지역 결정 시스템은 시스템 정보 독출부, 취약 지역 결정부, 결과 출력부를 포함한다. 시스템 정보 독출부는 전력 시스템의 정보를 독출한다. 취약 지역 결정부는 상기 시스템 정보, 및 상기 전력 시스템 내의 사고가 모의되는 사고 선로상의 지점과 상기 사고에 의해 전압 변화가 발생하는 목적 모선상의 전압 관계인 제 1 전압 강하식을 이용하여 사고 선로상의 양단 지점에 대응하는 목적 모선상의 전압을 산출하고, 양단 지점에 대응하는 목적 모선상의 전압이 모두 소정의 한계 전압보다 큰 경우가 아닌 경우 양단 지점 사이의 지점에 대응하는 목적 모선상의 전압을 산출하고, 사고 선로상의 양단 지점, 및 양단 지점 사이의 지점에 대응하는 목적 모선상의 전압들을 이용하여 제 1 전압 강하식을 단순화한 제 2 전압 강하식을 도출하고, 제 2 전압 강하식을 이용하여 사고 선로 상의 지점에 대응하는 상기 목적 모선상의 전압을 산출하며, 소정의 한계 전압 이하인 목적 모선상의 전압에 대응하는 사고 선로의 지점을 취약 지역으로 결정한다. 결과 출력부는 취약 지역 결정 결과를 출력한다.
순간 전압 강하, 취약 지역, 2차 보간, 할선법

Description

순간 전압 강하 평가를 위한 취약 지역 결정 시스템, 방법, 및 상기 방법을 실행시키기 위한 컴퓨터 판독 가능한 프로그램을 기록한 매체{System and method for determining areas of vulnerability for voltage sags assessment, and a medium having computer readable program for executing the method}
도 1은 본 발명에 따른 순간 전압 강하 평가를 위한 취약 지역 결정 시스템의 개략적인 블록도.
도 2는 사고 선로와 목적 모선을 포함하는 전력 시스템의 도면.
도 3은 본 발명에 따른 순간 전압 강하 평가를 위한 취약 지역 결정 방법의 개략적인 흐름도.
도 4는 도 3의 취약 지역 결정 단계의 더욱 상세한 흐름도.
도 5는 간단한 전력 시스템의 도면.
도 6은 도 5의 전력 시스템의 조류 계산 결과의 표
도 7은 모선 1과 2 사이 선로에서의 사고에 의한 모선 5에서의 전압 강하와 보간 곡선.
도 8은 모선 1과 3 사이 선로에서의 사고에 의한 모선 7에서의 전압 강하와 보간 곡선.
도 9는 IEEE-30 모선 15번 모선에 대한 취약 지역 계산 결과의 표
도 10은 도 9의 계산 결과를 전력 시스템에 도시한 도면.
본 발명은 전력 시스템에 관한 것으로서, 더욱 상세하게는 전력 품질을 개선하기 위한 시스템, 및 방법에 관한 것이다.
순간 전압 강하(voltage sag)는 가장 중요한 전력품질 문제 중 하나로 짧은 시간 동안의 전압강하 현상을 말한다. IEEE Std 1159-1995에 따르면 순간 전압 강하의 크기 범위는 0.1에서 0.9[p.u.]이고 지속시간은 0.5cycle에서 1분 범위로 정의된다. 일반적으로 순간 전압 강하는 전력 계통에서의 사고(fault)에 의해 많이 발생한다.
사고의 발생으로 계통의 전압은 떨어지게 되고 이로 인해 전압 변동에 민감한 부하들은 오동작을 하거나 동작을 멈추게 된다. 전압 변동에 민감한 부하의 사용은 계속 증가하고 있으며 순간 전압 강하 발생으로 인한 경제적 피해 또한 증가하고 있는 추세이다.
순간 전압 강하는 전압강하의 크기(magnitude)와 그 지속 시간(duration)이라는 중요한 두 가지 특성을 가진다. 순간 전압 강하의 전압 크기는 사고 위치와 전력 계통의 구성 형태에 영향을 받는다. 그 외에도 사고의 종류, 사고전 전압, 변압기 결선, 고장 임피던스 등의 영향도 받는다. 전압강하의 크기는 백분율이나 per-unit의 rms 전압으로 표시되며, 일반적으로 고장 계산에 의해서 구해진다.
순간 전압 강하의 지속 시간은 사고가 발생한 계통에서 고장 전류의 지속 시간과 같다고 볼 수 있다. 그러므로 실제적인 순간 전압 강하 지속시간은 브레이커나 릴레이, 퓨즈 등의 계통 보호 기기의 동작시간과 보호기기간 상호 작용을 고려한 약간의 시간 지연(time delay)의 합과 같다.
순간 전압 강하 평가를 위해서는 민감 부하에 대한 정확한 취약 지역(Area of Vulnerability)의 계산이 가장 중요하다. 취약 지역이라 함은 해당 부하 모선에 순간 전압 강하를 발생시키는 고장 위치를 말한다. 취약 지역을 계산하는 것은 이러한 고장 위치들을 모두 찾는 것이다. 정확한 취약 지역을 계산함으로써 특정 목적 모선에서의 연간 순간 전압 발생 횟수를 예측할 수 있다.
예를 들어 어떠한 민감 부하가 0.8[p.u.] 이하의 전압강하에 대해 오동작을 하거나 동작을 멈추게 된다면 그 부하의 한계 전압은 0.8[p.u.]이며 해당 부하에 대한 취약 지역은 부하가 연결되어 있는 모선의 전압을 0.8[p.u.] 이하로 만드는 계통 내 사고 위치들을 모두 찾는 것이다.
취약 취약지역 결정을 위한 종래의 방법들은 대표적으로 두 가지가 있다.
첫 번째는 critical distance법이 있다. 이 방법은 방사상(radial) 계통에서의 취약지역을 결정하는 방법으로 전압 분배 법칙에 의해 정확한 취약지역을 계산한다. 그러나 이 방법은 일반적인 망(mesh) 구조의 전력계통에는 적용되지 못하는 한계가 있다.
두 번째는 고장 위치(fault positions)법 이 있다. 이 방법은 모든 형태의 계통에 적용가능하나 그 정확도와 효율이 떨어지는 문제가 있다. 이 방법은 계통의 임의의 위치들에 사고를 모의하여 특정 부하 모선의 전압을 한계 값 이하로 만드는 위치를 찾는 방법으로 랜덤하게 위치를 선정하거나 선로의 길이에 비례하여 일정 구간별 위치를 선정하여 계산한다. 이 방법은 사고 모의 수와 모의 위치에 의해 그 정확도가 결정된다. 높은 정확도를 위해서는 많은 수의 사고 모의가 필요하여 비효율적인 방법이다.
본 발명은 모든 형태의 계통에 적용 가능하면서도 빠르고 정확하게 취약 지역 계산을 수행할 수 있는 시스템 및 방법을 제공하는 것을 목적으로 한다.
상기 문제점을 해결하기 위해 본 발명에 따른 순간 전압 강하 평가를 위한 취약 지역 결정 시스템은 시스템 정보 독출부, 취약 지역 결정부, 결과 출력부를 포함한다.
시스템 정보 독출부는 전력 시스템의 정보를 독출한다. 취약 지역 결정부는 상기 시스템 정보, 및 상기 전력 시스템 내의 사고가 모의되는 사고 선로상의 지점과 사고에 의해 전압 변화가 발생하는 목적 모선상의 전압 관계인 제 1 전압 강하식을 이용하여 사고 선로상의 양단 지점에 대응하는 목적 모선상의 전압을 산출하고, 사고 선로상의 양단 지점, 및 양단 지점 사이의 지점에 대응하는 목적 모선상의 전압들을 이용하여 제 1 전압 강하식을 단순화한 제 2 전압 강하식을 도출하고, 제 2 전압 강하식을 이용하여 사고 선로 상의 지점에 대응하는 목적 모선상의 전압을 산출하며, 소정의 한계 전압 이하인 목적 모선상의 전압에 대응하는 사고 선로 의 지점을 취약 지역으로 결정한다. 결과 출력부는 취약 지역 결정 결과를 출력한다.
민감 부하가 연결된 민감 부하 모선인 목적 모선상의 전압 분포와 한계 전압 사이의 관계를 이용하여 취약 지역을 결정함으로써 모든 형태의 계통에 적용 가능하면서도 빠르고 정확하게 취약 지역 계산을 수행할 수 있게 된다.
제 2 전압 강하식은 제 1 전압 강하식의 2차 보간식일 수 있고, 사고 선로상의 양단 지점 사이의 지점은 양단 지점의 중간 지점일 수 있다.
목적 모선상의 전압이 한계 전압과 동일한 경우의 사고 선로상의 지점을 산출하기 위해 할선법을 이용할 수 있다. 할선법을 이용함으로써, 취약 지역을 빠르고 더욱 정확하게 결정할 수 있게 된다.
아울러, 상기 시스템의 방법 형태의 발명과 상기 방법을 실행시키기 위한 컴퓨터 판독 가능한 프로그램을 기록한 매체가 함께 청구된다.
이하, 첨부된 도면을 참조하여 본 발명의 바람직한 실시예를 설명한다.
도 1은 본 발명에 따른 순간 전압 강하 평가를 위한 취약 지역 결정 시스템의 개략적인 블록도이다. 도 1에서 취약 지역 결정 시스템(100)은 전력 시스템(200)과 사용자 단말(300)과 연결되어 있으며, 시스템 정보 독출부(110), 취약 지역 결정부(120), 및 결과 출력부(130)를 포함하고 있다.
시스템 정보 독출부(110)는 전력 시스템(200)의 정보를 독출하고, 취약 지역 결정부(120)는 시스템 정보를 이용하여 소정의 선로에 대한 취약 정보를 결정한다. 결과 출력부(130)는 취약 지역 결정부(120)가 수행한 취약 지역 결정 결과를 사용 자 단말(300)로 출력한다.
취약 지역 결정부(120)가 취약 지역 결정하는 원리는 다음과 같다.
도 2는 사고 선로와 목적 모선을 포함하는 전력 시스템의 도면이다. 도 2와 같이 어떤 선로 F-T 사이의 임의의 사고 위치 K 가 있을 때 K에서의 사고에 의한 계통 내 임의의 모선 m의 전압 강하식을 유도할 수 있다.
불평형 사고에 의한 전압 강하 크기는 각 사고 형태 및 각 상에 따라 그 값을 달리하기 때문에 불평형 사고에 대한 취약지역 또한 각 사고 형태 및 각 상에 따라 따로따로 결정된다.
3상 평형 사고 및 불평형 사고 즉, 1선 지락 사고 (Single line-to-ground fault), 선간 단락 사고 (Line-to-line fault), 2선 지락 사고 (Double line-to-ground fault)에 대한 전압 강하식을 유도할 수 있다.
임의의 선로상의 사고로 인한 임의의 민감 부하 모선의 각 상 전압을 계산하기 위해서는 사고위치와 민감 부하 모선 사이의 sequence transfer impedance, 사고 위치에서의 sequence driving point impedance, 사고위치에서의 사고 전 전압을 알아야 한다.
사고위치 K와 모선m사이의 정상, 역상, 영상의 sequence transfer impedance는 다음과 같다. 위 첨자 0은 영상, 1은 정상, 2는 역상을 나타낸다.
Figure 112006009526068-PAT00001
Figure 112006009526068-PAT00002
: 사고위치 K와 모선F사이의 각 sequence transfer impedance
Figure 112006009526068-PAT00003
: 사고위치 K와 모선T사이의 각 sequence transfer impedance
p: 사고 위치와 전체 선로 길이에 대한 비 (0<=p<=1)
사고위치 K에서의 정상, 역상, 영상의 sequence driving point impedance는 다음과 같다.
Figure 112006009526068-PAT00004
Figure 112006009526068-PAT00005
: 모선 F에서의 각 sequence driving point impedance
Figure 112006009526068-PAT00006
: 모선 T에서의 각 sequence driving point impedance
Figure 112006009526068-PAT00007
: 선로의 각 sequence impedance
사고위치 K에서의 사고전 전압은 다음과 같다.
Figure 112006009526068-PAT00008
Figure 112006009526068-PAT00009
: 모선 F에서의 사고전 전압
Figure 112006009526068-PAT00010
: 모선 T에서의 사고전 전압
이상의 각 sequence impedance들과 사고전 전압식을 이용한 불평형 및 평형 고장에 대한 전압 강하식은 다음과 같다.
1) Single line-to-ground fault (SLGF)
임의의 선로 F-T 사이에서 A상에 SLGF가 발생했을 때 임의의 민감 부하 모선 m의 각 상전압은 다음과 같다.
Figure 112006009526068-PAT00011
식(1)
Figure 112006009526068-PAT00012
식(2)
Figure 112006009526068-PAT00013
식(3)
Figure 112006009526068-PAT00014
: 모선 m의 사고전 전압
a: complex number operator,
Figure 112006009526068-PAT00015
2) Line-to-line fault (LLF)
임의의 선로 F-T 사이에서 B와 C 상에서 LLF가 발생했을 때 임의의 목적 모선 m의 각 상전압은 다음과 같다. LLF에 의한 전압 강하 계산은 정상과 역상 성분만 필요하다.
Figure 112006009526068-PAT00016
식(4)
Figure 112006009526068-PAT00017
식(5)
Figure 112006009526068-PAT00018
식(6)
3) Double line-to-ground fault (DLGF)
임의의 선로 F-T 사이에서 B와 C 상에서 DLGF가 발생했을 때 임의의 목적 모선 m의 각 상전압은 다음과 같다.
Figure 112006009526068-PAT00019
식(7)
Figure 112006009526068-PAT00020
식(8)
Figure 112006009526068-PAT00021
식(9)
4) Three phase fault (3PF)
임의의 선로 F-T 사이에서 3PF가 발생했을 때 임의의 목적 모선 m의 각 상전압은 다음과 같다. 3PF에 의한 전압 강하 계산은 정상 성분만 필요하다.
Figure 112006009526068-PAT00022
식(10)
위의 전압 강하식들은 모두 p에 대한 식으로서, 유도된 식들을 이용해 특정 목적 모선 m의 전압의 크기를 얼마 이하로 만드는 임의의 선로 F-T에서의 사고위치(critical point) p를 계산할 수 있다.
이때의 p를 critical point라고 한다. 불평형 사고에 있어 사고가 발생한 상(faulted phase)의 전압은 급격히 떨어지게 되고 사고가 나지 않은 상(non-faulted phase)의 전압은 일반적으로 사고전보다 올라가거나 변화가 없다. 물론 변압기 결선이라던지 각 sequence impedance 간의 차이로 인해 non-faulted phase의 전압도 떨어질 수 있다.
그러나 순간 전압 강하 해석에서는 전압이 떨어지는 상만이 주 관심사이며 일반적으로 임의의 선로 F-T를 따라 발생하는 사고에 의한 특정 목적 모선 m에서의 전압 강하는
Figure 112006009526068-PAT00023
에 대해 2차 곡선과 유사한 형태를 가진다.
이에 빠른 계산을 위해 한 선로당 critical point는 최대 2개까지 존재하며 선로 양단 모선에서의 사고로 인한 특정 목적 모선에서의 전압 강하는 선로상의 사고로 인한 전압 강하보다 크다라고 가정한다. 대부분의 경우 이 조건을 만족한다. 이러한 조건으로부터 다음을 유추할 수 있다.
1) 민감 부하의 전압 민감도(한계값)가 p=0과 p=1에서의 전압강하 크기보다 작다면 해당 선로는 취약지역에 포함되지 않는다.
2) 민감 부하의 전압 민감도(한계값)가 해당선로에 대한 전압 강하식의 최대값보다 큰 값이면 해당 선로 전체가 취약지역에 포함된다.
3) 민감 부하의 전압 민감도(한계값)가 p=0에서의 전압강하 크기와 p=1에서 의 전압강하 크기 사이의 값이면 선로에는 하나의 critical point가 존재하며 해당선로의 일정한 부분만이 취약지역에 포함되게 된다.
4) 민감 부하의 전압 민감도(한계값)가 해당선로에 대한 전압 강하식의 최대값보다 작고 p=0과 p=1에서의 전압강하 크기보다 크면 선로에는 두 개의 critical point가 존재하며 해당선로의 두 부분이 취약지역에 포함되게 된다.
위 사실로부터 우리는 계통 내 모든 선로에 대해 전압 강하식과 위의 4가지 사실을 적용해 정확한 취약지역을 얻을 수 있다.
선로의 critical point를 찾는 방법은 여러 가지가 있을 수 있다. 본 발명에서는 2차 보간법(quadratic interpolation)과 할선법(secant method)을 이용한다. 이 방법은 방사 및 망구조의 모든 계통에 다 적용가능하며 빠르고 정확하게 취약지역 계산을 할 수 있는 장점이 있다.
할선법은 비선형 방정식의 해를 찾는 수치해석방법 중 하나이며 뉴턴-랩슨 법에 비해 수렴 속도는 떨어지나 미분 계산이 필요 없어 전압 강하식과 같이 미분식을 유도하기 어려운 경우 이점이 많다.
비선형 방정식 f(x)=0의 해를 구하기 위한 일반적인 할선법은 다음과 같다.
1) 반복 루프 (f(x): 비선형 방정식, x from , x end : 2개의 초기값, x new : 새로운 추정값)
Figure 112006009526068-PAT00024
Figure 112006009526068-PAT00025
2) 수렴 판단
Figure 112006009526068-PAT00026
위의 일반적인 할선법을 전압 강하식에 적용하면 다음과 같다. 취약지역 계산에는 전압의 크기만을 고려하므로 f(p)의 크기 즉
Figure 112006009526068-PAT00027
가 목적 함수가 된다.
반복 루프 (f(p): 전압 강하식, p from , p end : 2개의 초기값, p new : 새로운 추정값)
Figure 112006009526068-PAT00028
Figure 112006009526068-PAT00029
수렴 판단
Figure 112006009526068-PAT00030
이것은
Figure 112006009526068-PAT00031
형태의 전압 강하식의 해를 구하는 꼴로 일반적인 할선법이다. 이것은 전압 강하식의 크기를 0으로 만드는 p를 찾는 것이다.
그러나 특정 한계 전압값에 대한 전압 강하식을
Figure 112006009526068-PAT00032
형태로 변환하지 않고 바로 해를 구하기 위해 다음과 같이 수정할 수 있다. 즉,
Figure 112006009526068-PAT00033
형태의 전압 강하식에 바로 적용할 수 있게 간단히 수정할 수 있다. 이의 장점은 다양한 한계 전압값에 대해서
Figure 112006009526068-PAT00034
형태로의 변환 없이 바로 적용할 수 있다는 것이다.
반복 루프 (f(p): 전압 강하식, p from , p end : 2개의 초기값, p new : 새로운 추정값)
Figure 112006009526068-PAT00035
Figure 112006009526068-PAT00036
수렴 판단
Figure 112006009526068-PAT00037
이와 같은 할선법 사용의 장점은 다음과 같다. 해를 찾기 위해 전압 강하식의 미분 계산이 필요 없다. 또한, 전압 강하식을
Figure 112006009526068-PAT00038
형태로 변환할 필요가 없기 때문에 다양한 전압 민감도 즉, 전압 한계값에 대해 식의 변형 과정 없이 바로 계산할 수 있는 이점이 있다.
할선법을 이용해 해를 찾기 위해서는 2개의 (초기값)starting value가 필요하다. 더욱 빠르게 해를 찾기 위해서는 해에 가까운 초기값의 선정이 필요하다. 이를 위해서는 제안하는 방법에서는 전압 강하식에 대한 2차 보간을 이용한다.
2차 보간법에는 대표적으로 Newton 보간법과 Lagrange 보간법이 있다. 두 방법 모두 나름의 장점이 있으나 2차와 같은 저차 보간에서는 그 성능의 차이가 없으며 어느 것을 사용해도 무방하다. 2차 보간의 경우 프로그램이 간단하여 계산부담 이 없다. 2차 보간식을
Figure 112006009526068-PAT00039
라고 하면 근의 공식을 이용해 쉽게 해를 찾을 수 있다.
이상의 방법들은 3상 평형 사고 및 불평형 사고인 1선 지락, 2선 지락, 선간 단락 사고에도 적용 가능하다.
도 3은 본 발명에 따른 순간 전압 강하 평가를 위한 취약 지역 결정 방법의 개략적인 흐름도이다.
먼저, 순간 전압 강하 평가를 위한 취약 지역 결정 시스템이 목적하는 전력 시스템의 정보를 독출한다(S110). 다음으로, 독출한 시스템 정보를 이용하여 전력 시스템 내의 목적하는 선로에 대한 취약 지역을 결정한다(S120). 마지막으로, 취약 지역 결정 결과를 소정의 출력 장치로 출력한다(S130).
도 4는 도 3의 취약 지역 결정 단계의 더욱 상세한 흐름도이다. 도 4에서 민감 부하 모선의 한계 전압에 대한 취약지역 계산 순서는 다음과 같다.
1) 계통의 사고전 전압을 알기 위해 조류계산을 수행한다.
2) 각 고장 형태에 따라 계산에 필요한positive, negative, zero sequence 어드미턴스 행렬을 구성하고, 각각에 대해 LU 분해(삼각 분해)를 수행한다.
3) 선로의 두 연결 모선에 대한 sequence 임피던스 행렬의 각 요소를 LU 행렬에 대한 forward 와 back substitution을 통해 계산한다.
4) 1)과 3)의 계산 결과를 이용해 전압 강하식을 형성한다.
5) 구성된 전압 강하식을 이용해 p=0과 p=1에서의 전압강하 크기를 계산한 다. 주어진 한계 전압값이 p=0와 p=1에서의 전압강하 값보다 작으면 해당 선로는 취약지역에 포함되지 않으며 다음 선로에 대해 단계3)부터 다시 적용해 나간다. 주어진 한계 전압값이 p=0 또는 p=1에서의 전압강하 값보다 크면 다음 단계로 넘어간다.
6) 구성된 전압 강하식을 이용해 p=0.5에서의 전압강하 값을 계산한다.
7) 단계5)와 6)에서 계산한 p=0, 0.5, 1에서의 전압강하 값 세 개를 이용해 2차 보간을 실시한다. 주어진 한계 전압값이 p=0와 p=1에서의 전압강하 값의 사이 값이면 단계9)로 넘어간다. 주어진 한계 전압값이 p=0와 p=1에서의 전압강하 값보다 크면 다음 단계로 넘어간다.
8) 보간 2차식에서 최대일 때의 p값
Figure 112006009526068-PAT00040
를 계산한다.
Figure 112006009526068-PAT00041
는 0과 1 사이의 값이어야 한다. 이
Figure 112006009526068-PAT00042
를 구성된 전압 강하식에 대입하여 계산한다. 이때의 값을 전압 강하식의 최대값으로 가정한다. 주어진 한계 전압값이 이 최대값보다 크거나 2차 보간 방정식의 판별식이 0보다 작거나 같으면 해당선로 전체를 취약지역에 포함시키고 다음 선로에 대해 단계3)부터 다시 적용해 나간다. 주어진 한계 전압값이 이 최대값보다 작고 p=0와 p=1에서의 전압강하 값보다 크면 해당선로에는 2개의 critical point 가 존재하게 된다.
9) 근의 공식을 이용해 2차 보간식의 해를 계산한다. 보간식의 해가
Figure 112006009526068-PAT00043
라면 두 초기값은
Figure 112006009526068-PAT00044
와
Figure 112006009526068-PAT00045
또는
Figure 112006009526068-PAT00046
와
Figure 112006009526068-PAT00047
로 결정한다.
10) 단계9)에서 결정된 초기값과 할선법을 이용해 critical point를 계산한다.
11) 다음 선로에 대해 단계3)부터 다시 적용해 나간다. 모든 선로에 대해 다 적용될 때까지 행한다.
12) 계산된 critical point들을 이용해 취약지역을 결정한다.
도 5는 본 발명에 따른 순간 전압 강하 평가를 위한 취약 지역 결정 방법을 적용하기 위한 간단한 7 모선 계통의 도면이고, 도 6은 도 5의 전력 시스템의 조류 계산 결과의 표이다.
본 실시예에서는 두 선로에 대해서 본 발명을 적용해 3상 평형 사고에 대한 취약 지역을 결정해 본다. 먼저, 모선 1과 2 사이의 선로에서의 취약 지역을 결정한다.
민감 부하 모선은 모선 5로 가정하며 한계 전압값은 0.25[p.u.]로 가정한다. 해당선로의 정상 임피던스는
Figure 112006009526068-PAT00048
이다.
3상 평형 사고에 대해서는 정상 성분만 고려하면 된다. 모선 1과 2에 대한 driving-point impedance와 transfer impedance를 LU 분해를 통해 계산한다. 결과는 다음과 같다.
Figure 112006009526068-PAT00049
위의 결과들을 이용해 식(10)을 구성한다.
Figure 112006009526068-PAT00050
위의 식을 이용해 p=0과 p=1일 때의 전압강하 크기를 계산하면 각각 0.3704와 0.1533[p.u.]이다. 주어진 한계전압이 0.25이고 이 두 값의 사이 값이므로 해당선로에는 하나의 critical point가 존재한다. p=0.5에서의 전압강하 크기를 계산하면 0.2903이다. 위의 세 값을 이용해 2차 보간을 실시한다. 보간 방정식은 다음과 같다.
Figure 112006009526068-PAT00051
근의 공식을 이용해
Figure 112006009526068-PAT00052
를 만족하는 해
Figure 112006009526068-PAT00053
를 구하면 0.6703 하나뿐이다.
도 7은 모선 1과 2 사이 선로에서의 3상 평형 사고에 의한 모선 5에서의 전압강하 크기와 2차 보간식의 그래프를 나타낸다.
더욱 정확한 해를 찾기 위해 할선법을 적용한다. 두 초기값은
Figure 112006009526068-PAT00054
와
Figure 112006009526068-PAT00055
로 한다.
할선법에 의해 계산된 정확한 critical point
Figure 112006009526068-PAT00056
는 0.6771이다. 그러므로 해당선로에서의 취약지역은
Figure 112006009526068-PAT00057
의 범위이다. 즉, 전체 선로 길이의 약 67%에서 끝까지가 취약지역에 포함된다. 본 계산에서 사용된 수렴 오차는 0.0001이었으며 한 번의 반복으로 해에 수렴해 아주 빠르게 계산됨을 확인하였다.
다음으로, 모선 1과 3 사이의 선로에서의 취약지역을 결정한다.
민감 부하 모선은 모선 7로 가정하며 한계 전압값은 0.5[p.u.]로 가정한다. 해당선로의 정상 임피던스는
Figure 112006009526068-PAT00058
이다. 모선 1과 3에 대한 driving-point impedance와 transfer impedance를 LU 분해를 통해 계산한다. 결과는 다음과 같다.
Figure 112006009526068-PAT00059
Figure 112006009526068-PAT00060
위의 결과들을 이용해 식(10)을 구성한다.
Figure 112006009526068-PAT00061
위의 식을 이용해 p=0과 p=1일 때의 전압강하 크기를 계산하면 각각 0.4507과 0.4758[p.u.]이다. p=0.5에서의 전압강하 크기를 계산하면 0.5956이다. 위의 세 값을 이용해 2차 보간을 실시한다. 보간 방정식은 다음과 같다.
Figure 112006009526068-PAT00062
주어진 한계 전압이 0.5이므로 p=0 와 1에서의 전압강하 크기보다 크다. 보 간식의 최대값일 때의
Figure 112006009526068-PAT00063
는 0.5238이고 이때의 전압강하 크기는 0.5960 [p.u.] 이다. 그러므로 해당 선로는 두 개의 critical point를 가지며 선로의 두 부분이 취약지역에 포함되게 된다.
근의 공식을 이용해
Figure 112006009526068-PAT00064
을 만족하는 두 개의 해
Figure 112006009526068-PAT00065
과
Figure 112006009526068-PAT00066
를 구하면 각각 0.0981과 0.9496이다.
도 8은 모선 1과 3 사이 선로에서의 3상 평형 사고에 의한 모선 7에서의 전압강하 크기와 2차 보간식의 그래프를 나타낸다.
더욱 정확한 해를 찾기 위해 할선법을 적용한다. 두 개의
Figure 112006009526068-PAT00067
,
Figure 112006009526068-PAT00068
에 대한 각각의 초기값은
Figure 112006009526068-PAT00069
과
Figure 112006009526068-PAT00070
,
Figure 112006009526068-PAT00071
과
Figure 112006009526068-PAT00072
로 한다.
할선법에 의해 계산된 정확한 critical point
Figure 112006009526068-PAT00073
은 0.0841,
Figure 112006009526068-PAT00074
는 0.9607이다. 그러므로 해당선로에서의 취약지역은
Figure 112006009526068-PAT00075
과
Figure 112006009526068-PAT00076
의 두 부분이다. 본 계산에서 사용된 수렴 오차는 0.0001이었으며 각 두 번의 반복으로 해에 수렴해 아주 빠르게 계산됨을 확인하였다.
마지막으로, IEEE-30 모선 계통에 대해 본 발명을 적용해 본다.
IEEE-30 모선 계통은 30개의 모선과 37개의 선로, 4개의 변압기로 연결되어있다. 수렴 오차는 0.0001로 하여 계산하였다. 목적 모선은 15번 모선으로 가정하였으며 한계 전압은 0.7[p.u.]로 가정하였다. 모든 발전기의 정상, 역상, 영상 impedance는 각 j0.3, j0.2, j0.05로 가정하였고 모든 변압기 결선은 grounded wye-grounded wye로 가정하였다. 사고는 B와 C 상에서의 LLF를 가정하였고 각 B 상과 C상에 대한 취약지역을 계산하였다. B상에 대한 취약지역은 전압 강하식 (5)를 C상에 대한 취약지역은 전압 강하식 (6)을 구성하여 계산된다.
도 9는 IEEE-30 모선 계통의 15번 모선에 대한 취약 지역 계산 결과의 표이다. 도 9의 표에는 B와 C 상에 대한 각 취약지역과 할선법 반복횟수를 같이 나타내었다.
도 10은 도 9의 계산 결과를 전력 시스템에 도시한 도면이다. 도 10의 실선과 점선 영역에서의 LLF는 15번 모선 부하에 순간 전압 강하를 발생시킨다는 것을 보여 준다.
이상에서, 취약지역을 빠르고 정확하게 계산할 수 있는 방법이 설명되었다. 본 발명은 3상 평형 사고뿐 아니라 1선 지락, 2선 지락, 선간 단락과 같은 불평형 사고에 대해서도 모두 적용 가능하다.
본 발명에 의해 모든 형태의 계통에 적용 가능하면서도 빠르고 정확하게 순간 전압 강하 평가를 위한 취약 지역 계산을 수행할 수 있게 된다.
본 발명은 2차 보간 다항식과 전압 강하식만을 이용하여 다양한 한계 전압값에 대한 취약지역을 빠르게 계산할 수 있는 장점이 있다.
본 발명이 비록 일부 바람직한 실시예에 의해 설명되었지만, 본 발명의 범위는 이에 의해 제한되어서는 아니 되고, 첨부된 특허청구범위에 의해 뒷받침되는 상기 실시예의 변형이나 개량에도 미쳐야 할 것이다.

Claims (9)

  1. 전력 시스템의 정보를 독출하는 시스템 정보 독출부;
    상기 시스템 정보, 및 상기 전력 시스템 내의 사고가 모의되는 사고 선로상의 지점과 상기 사고에 의해 전압 변화가 발생하는 목적 모선상의 전압 관계인 제 1 전압 강하식을 이용하여 상기 사고 선로상의 양단 지점에 대응하는 상기 목적 모선상의 전압을 산출하고,
    상기 양단 지점에 대응하는 목적 모선상의 전압이 모두 소정의 한계 전압보다 큰 경우가 아닌 경우 상기 양단 지점 사이의 지점에 대응하는 목적 모선상의 전압을 산출하고,
    상기 사고 선로상의 양단 지점, 및 양단 지점 사이의 지점에 대응하는 상기 목적 모선상의 전압들을 이용하여 상기 제 1 전압 강하식을 단순화한 제 2 전압 강하식을 도출하고,
    상기 제 2 전압 강하식을 이용하여 상기 사고 선로 상의 지점에 대응하는 상기 목적 모선상의 전압을 산출하며,
    상기 소정의 한계 전압 이하인 목적 모선상의 전압에 대응하는 사고 선로의 지점을 취약 지역으로 결정하는 취약 지역 결정부: 및
    상기 취약 지역 결정 결과를 출력하는 결과 출력부를 포함하는 것을 특징으로 하는 순간 전압 강하 평가를 위한 취약 지역 결정 시스템.
  2. 제 1항에 있어서,
    상기 양단 지점 사이의 지점은 상기 양단 지점의 중간 지점인 것을 특징으로 하는 순간 전압 강하 평가를 위한 취약 지역 결정 시스템.
  3. 제 1항에 있어서,
    상기 제 2 전압 강하식은 상기 사고 선로상의 지점에 대한 상기 제 1 전압 강하식의 2차 보간식인 것을 특징으로 하는 순간 전압 강하 평가를 위한 취약 지역 결정 시스템.
  4. 제 3항에 있어서,
    상기 목적 모선상의 전압이 상기 한계 전압과 동일한 경우의 상기 사고 선로상의 지점을 산출하기 위해 할선법을 이용하는 것을 특징으로 하는 순간 전압 강하 평가를 위한 취약 지역 결정 시스템.
  5. 순간 전압 강하 평가를 위한 취약 지역 결정 시스템이 목적하는 전력 시스템의 정보를 독출하는 단계;
    상기 시스템 정보, 및 상기 전력 시스템 내의 사고가 모의되는 사고 선로상의 지점과 상기 사고에 의해 전압 변화가 발생하는 목적 모선상의 전압 관계인 제 1 전압 강하식을 이용하여 상기 사고 선로상의 양단 지점에 대응하는 상기 목적 모선상의 전압을 산출하는 단계;
    상기 양단 지점에 대응하는 목적 모선상의 전압이 모두 소정의 한계 전압보다 큰 경우가 아닌 경우 상기 양단 지점 사이의 지점에 대응하는 목적 모선상의 전압을 산출하는 단계;
    상기 사고 선로상의 양단 지점, 및 양단 지점 사이의 지점에 대응하는 상기 목적 모선상의 전압들을 이용하여 상기 제 1 전압 강하식을 단순화한 제 2 전압 강하식을 도출하는 단계;
    상기 제 2 전압 강하식을 이용하여 상기 사고 선로 상의 지점에 대응하는 상기 목적 모선상의 전압을 산출하는 단계;
    상기 소정의 한계 전압 이하인 목적 모선상의 전압에 대응하는 사고 선로의 지점을 취약 지역으로 결정하는 단계; 및
    상기 취약 지역 결정 결과를 출력하는 단계를 포함하는 것을 특징으로 하는 순간 전압 강하 평가를 위한 취약 지역 결정 방법.
  6. 제 5항에 있어서,
    상기 양단 지점 사이의 지점은 상기 양단 지점의 중간 지점인 것을 특징으로 하는 순간 전압 강하 평가를 위한 취약 지역 결정 방법.
  7. 제 5항에 있어서,
    상기 제 2 전압 강하식은 상기 제 1 전압 강하식의 2차 보간식인 것을 특징으로 하는 순간 전압 강하 평가를 위한 취약 지역 결정 방법.
  8. 제 7항에 있어서,
    상기 목적 모선상의 전압이 상기 한계 전압과 동일한 경우의 상기 사고 선로상의 지점을 산출하기 위해 할선법을 이용하는 것을 특징으로 하는 순간 전압 강하 평가를 위한 취약 지역 결정 방법.
  9. 제 5 내지 8항 중 어느 한 항의 방법을 실행시키기 위한 컴퓨터 판독 가능한 프로그램을 기록한 매체.
KR1020060012419A 2006-02-09 2006-02-09 순간 전압 강하 평가를 위한 취약 지역 결정 시스템, 방법,및 상기 방법을 실행시키기 위한 컴퓨터 판독 가능한프로그램을 기록한 매체 KR100782354B1 (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020060012419A KR100782354B1 (ko) 2006-02-09 2006-02-09 순간 전압 강하 평가를 위한 취약 지역 결정 시스템, 방법,및 상기 방법을 실행시키기 위한 컴퓨터 판독 가능한프로그램을 기록한 매체

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020060012419A KR100782354B1 (ko) 2006-02-09 2006-02-09 순간 전압 강하 평가를 위한 취약 지역 결정 시스템, 방법,및 상기 방법을 실행시키기 위한 컴퓨터 판독 가능한프로그램을 기록한 매체

Publications (2)

Publication Number Publication Date
KR20070080913A true KR20070080913A (ko) 2007-08-14
KR100782354B1 KR100782354B1 (ko) 2007-12-11

Family

ID=38601142

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020060012419A KR100782354B1 (ko) 2006-02-09 2006-02-09 순간 전압 강하 평가를 위한 취약 지역 결정 시스템, 방법,및 상기 방법을 실행시키기 위한 컴퓨터 판독 가능한프로그램을 기록한 매체

Country Status (1)

Country Link
KR (1) KR100782354B1 (ko)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101106120B1 (ko) * 2010-06-04 2012-01-20 고려대학교 산학협력단 신호 교정 장치 및 방법
CN107436396A (zh) * 2017-09-05 2017-12-05 南方电网科学研究院有限责任公司 一种考虑电压暂降的电网薄弱环节的辨识方法
CN109828184A (zh) * 2019-01-18 2019-05-31 东南大学 一种基于互近似熵的电压暂降源识别方法
CN111157837A (zh) * 2018-11-08 2020-05-15 中国石油化工股份有限公司 一种油田配电网的电压暂降凹陷域的分析方法
CN111157836A (zh) * 2018-11-08 2020-05-15 中国石油化工股份有限公司 一种适用于油田配电网确定故障线路范围的方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100425417B1 (ko) * 2002-05-06 2004-03-30 명지대학교 전력 계통의 선간 단락 고장 지점 검출 방법
KR100538563B1 (ko) * 2004-11-22 2005-12-27 대광전기공업(주) 전력계통 모의시험 시스템 및 그 방법

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101106120B1 (ko) * 2010-06-04 2012-01-20 고려대학교 산학협력단 신호 교정 장치 및 방법
CN107436396A (zh) * 2017-09-05 2017-12-05 南方电网科学研究院有限责任公司 一种考虑电压暂降的电网薄弱环节的辨识方法
CN107436396B (zh) * 2017-09-05 2019-08-02 南方电网科学研究院有限责任公司 一种考虑电压暂降的电网薄弱环节的辨识方法
CN111157837A (zh) * 2018-11-08 2020-05-15 中国石油化工股份有限公司 一种油田配电网的电压暂降凹陷域的分析方法
CN111157836A (zh) * 2018-11-08 2020-05-15 中国石油化工股份有限公司 一种适用于油田配电网确定故障线路范围的方法
CN111157837B (zh) * 2018-11-08 2021-11-05 中国石油化工股份有限公司 一种油田配电网的电压暂降凹陷域的分析方法
CN111157836B (zh) * 2018-11-08 2021-11-05 中国石油化工股份有限公司 一种适用于油田配电网确定故障线路范围的方法
CN109828184A (zh) * 2019-01-18 2019-05-31 东南大学 一种基于互近似熵的电压暂降源识别方法
CN109828184B (zh) * 2019-01-18 2020-11-03 东南大学 一种基于互近似熵的电压暂降源识别方法

Also Published As

Publication number Publication date
KR100782354B1 (ko) 2007-12-11

Similar Documents

Publication Publication Date Title
KR20090010536A (ko) 순간 전압 강하 평가를 위한 취약 지역 결정 시스템, 방법,및 상기 방법을 실행시키기 위한 컴퓨터 판독 가능한프로그램을 기록한 매체
Su et al. Estimating the voltage stability margin using PMU measurements
McLaren et al. Software models for relays
Ahsaee et al. A novel fault-location algorithm for long transmission lines compensated by series FACTS devices
CN108603911B (zh) 用于确定三相串联补偿输电线中的故障位置的系统和方法
Yu et al. Wide-area backup protection and protection performance analysis scheme using PMU data
Alzaareer et al. A new sensitivity approach for preventive control selection in real-time voltage stability assessment
KR100782354B1 (ko) 순간 전압 강하 평가를 위한 취약 지역 결정 시스템, 방법,및 상기 방법을 실행시키기 위한 컴퓨터 판독 가능한프로그램을 기록한 매체
CN111141995B (zh) 一种基于比幅原理的线路双端稳态量测距方法和系统
JP2006204069A (ja) 単独運転検出方法および単独運転検出装置
Khan et al. A phase-shifting transformer protection technique based on directional comparison approach
Mrehel et al. Implementation and Evaluation a SIMULINK Model of a Distance Relay in MATLAB/SIMULINK
CN111937264A (zh) 用于在多终端电力传输系统中进行保护的方法和装置
Das et al. Distribution fault location using short-circuit fault current profile approach
Alwash et al. New impedance‐based fault location method for unbalanced power distribution systems
Ghaedi et al. Modified WLS three-phase state estimation formulation for fault analysis considering measurement and parameter errors
KR20190091859A (ko) 보호계전기 설치적합성 검사방법 및 검사장치
Leite et al. State estimation of distribution networks through the real-time measurements of the smart meters
Altonen et al. Novel algorithm for earth-fault location in compensated MV-networks
Bidadfar et al. Dynamic Thévenin equivalent and reduced network models for PMU-based power system voltage stability analysis
CN110879332A (zh) 一种适用于小电流接地系统的单相接地故障选相方法
Polymeneas et al. Margin-based framework for online contingency selection in unbalanced networks
Kader et al. A novel directional comparison pilot protection scheme based on neutral-point zero-sequence current of UHV auto-transformer
Enriquez et al. Enhanced time overcurrent coordination
Dam et al. A breaker‐oriented fault analysis methodology

Legal Events

Date Code Title Description
A201 Request for examination
E902 Notification of reason for refusal
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20110914

Year of fee payment: 5

FPAY Annual fee payment

Payment date: 20121004

Year of fee payment: 6

LAPS Lapse due to unpaid annual fee