KR102238494B1 - Sensorless control System of DFIG - Google Patents

Sensorless control System of DFIG Download PDF

Info

Publication number
KR102238494B1
KR102238494B1 KR1020200005324A KR20200005324A KR102238494B1 KR 102238494 B1 KR102238494 B1 KR 102238494B1 KR 1020200005324 A KR1020200005324 A KR 1020200005324A KR 20200005324 A KR20200005324 A KR 20200005324A KR 102238494 B1 KR102238494 B1 KR 102238494B1
Authority
KR
South Korea
Prior art keywords
rotor
dfig
current
unit
axis
Prior art date
Application number
KR1020200005324A
Other languages
Korean (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 KR1020200005324A priority Critical patent/KR102238494B1/en
Application granted granted Critical
Publication of KR102238494B1 publication Critical patent/KR102238494B1/en

Links

Images

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03DWIND MOTORS
    • F03D7/00Controlling wind motors 
    • F03D7/02Controlling wind motors  the wind motors having rotation axis substantially parallel to the air flow entering the rotor
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05BINDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
    • F05B2270/00Control
    • F05B2270/30Control parameters, e.g. input parameters
    • F05B2270/337Electrical grid status parameters, e.g. voltage, frequency or power demand
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05BINDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
    • F05B2270/00Control
    • F05B2270/70Type of control algorithm
    • F05B2270/703Type of control algorithm integral
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05BINDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
    • F05B2270/00Control
    • F05B2270/70Type of control algorithm
    • F05B2270/705Type of control algorithm proportional-integral
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/72Wind turbines with rotation axis in wind direction

Abstract

The present invention relates to a sensor-less control system for a double-fed induction generator (DFIG) wind turbine. More particularly, the present invention is to provide a sensor-less control system for a DFIG wind turbine, which comprises: a first converter connected to rotor windings of the DFIG; and a control unit for controlling the first converter. The control unit uses an angular velocity estimator for estimating the angular velocity of the DFIG rotor based on the current flowing in the DFIG rotor, so that the turbine of the wind generator can be accurately controlled without using a sensor.

Description

DFIG 센서리스 제어 시스템{Sensorless control System of DFIG}DFIG sensorless control system {Sensorless control System of DFIG}

본 발명은 DFIG(Doubly-Fed Induction Generator) 센서리스 제어 시스템에 관한 것으로, 보다 상세하게는 풍력 발전기의 터빈(Turbine)이 센서를 사용하지 않고, 강인하고, 정확하게 제어될 수 있는, DFIG 센서리스 제어 시스템을 제공함에 있다.The present invention relates to a DFIG (Doubly-Fed Induction Generator) sensorless control system, and more particularly, a DFIG sensorless control that can be controlled robustly and accurately without the use of a sensor for the turbine of a wind power generator. It is in providing the system.

신재생 에너지는 최근 화석 연료의 사용으로 인한 환경 오염 문제, 그리고 점차 화석 연료가 고갈되는 문제로 인해 대체 에너지원으로 주목을 받고 있다. 신재생 에너지원으로는 태양광 에너지, 풍력 에너지, 바이오 에너지, 조력 에너지 등 다양한 종류가 있지만, 가장 연구가 활발한 분야는 태양광 발전, 풍력 발전 분야가 대표적이다. 그 중에서도 풍력 발전 분야는 모터, 발전기와 같은 전기기기 분야 및 컨버터, 인버터와 같은 전력변환장치 분야에서 지속적인 연구 개발로 많은 발전이 이루어졌다.Renewable energy has recently attracted attention as an alternative energy source due to the problem of environmental pollution due to the use of fossil fuels and the problem of gradually depleting fossil fuels. There are various types of renewable energy sources such as solar energy, wind energy, bio energy, and tidal energy, but the most active fields of research are solar power generation and wind power generation. Among them, a lot of development has been made through continuous research and development in the field of electric devices such as motors and generators, and power conversion devices such as converters and inverters.

풍력 발전기의 AC 모터 구동 시스템에서, AC 모터의 회전자 위치를 파악하는 것은 풍력 발전기의 터빈을 제어하는 것에 있어서 아주 중요한 기술적 요소 중 하나이다. AC 모터의 회전자 위치를 파악하기 위해서 종래에는 엔코더와 같은 위치 및 속도 검출 센서를 많이 사용하였다. 그러나 센서의 가격이 비싸기 때문에, 전체 시스템을 구성하기 위한 원가가 증가한다는 단점이 있고, 위치 및 속도 검출 센서는 온도, 습도 등 주위 환경과 시스템 내부 파라미터 등의 영향을 받을 수 있기 때문에, 검출된 값에 대한 신뢰성이 떨어질 수 있다. 또한, 발전기의 작동 중에 위치 검출 센서가 파손될 우려가 있기 때문에, 유지 보수에 대한 부담이 있다는 단점이 있다.In the AC motor drive system of a wind generator, knowing the position of the rotor of the AC motor is one of the very important technical factors in controlling the turbine of the wind generator. In order to grasp the rotor position of an AC motor, a position and speed detection sensor such as an encoder has been used a lot in the past. However, since the price of the sensor is expensive, there is a disadvantage that the cost for configuring the entire system increases, and since the position and speed detection sensor can be affected by the surrounding environment such as temperature and humidity, and the system internal parameters, the detected value The reliability of this may be reduced. In addition, there is a disadvantage in that there is a burden on maintenance because there is a possibility that the position detection sensor may be damaged during the operation of the generator.

상기와 같은 단점으로 인해, AC 모터의 회전자 위치를 파악하고, 풍력 발전 시스템을 제어함에 있어서, 위치 및 속도 검출 센서를 사용하지 않고, 센서를 대체할 수 있는 알고리즘을 바탕으로 한 소프트웨어를 설계함으로써, 비용을 절감하고, 신뢰성을 높일 수 있는 방법에 대한 많은 연구 결과가 제시되고 있다. 일례로 폐루프 제어 방식 중에서도 위상동기루프(PLL, Phase Locked Loop)를 포함하여 설계하는 방법이 대표적이다.Due to the above drawbacks, by designing a software based on an algorithm that can replace the sensor, without using a position and speed detection sensor, in determining the position of the rotor of the AC motor and controlling the wind power generation system. In addition, many research results have been suggested on ways to reduce cost and increase reliability. For example, among the closed loop control methods, a design method including a phase locked loop (PLL) is a representative method.

상기 위상동기루프를 사용하여 AC 모터의 회전자의 위치 및 속도를 검출하는 방식은 한국등록특허 제 10-1091970호(“비선형 관측기를 이용한 표면부착형 영구자석 동기전동기의 센서리스 제어 방법 및 시스템”, 선행기술 1)에서 개시하고 있다.The method of detecting the position and speed of the rotor of an AC motor using the phase-locked loop is described in Korean Patent No. 10-1091970 (“Sensorless control method and system of a surface-attached permanent magnet synchronous motor using a nonlinear observer”. , Disclosed in prior art 1).

도 1은 상기 선행기술 1의 일실시예에 따른 전체 상세 제어 시스템의 구성도를 도시하고 있다.1 is a block diagram of an overall detailed control system according to an embodiment of the prior art 1.

도 1에 도시된 바와 같이, 회전자의 위치 검출기(100)은 비선형 관측기(110)를 기반으로 설계되고, 비선형 관측기(110)는 PWM 인버터에서 출력되는 3상 중 2 상의 전류값 및 단자 전압을 입력받아 상기 비선형 관측기(110)에서 설계된 수학적 모델을 거쳐 추정 회전자 위치에 대한 sin값과 cos값을 출력하고, 상기 cos값과 sin값을 기반으로 역 탄젠트 함수(120)를 통해 추정 회전자 위치가 검출된다.As shown in Figure 1, the rotor position detector 100 is designed based on the nonlinear observer 110, the nonlinear observer 110 is the current value and the terminal voltage of the two phases output from the PWM inverter. It receives the input and outputs the sin value and the cos value for the estimated rotor position through the mathematical model designed by the nonlinear observer 110, and the estimated rotor position through the inverse tangent function 120 based on the cos value and the sin value. Is detected.

상기 검출된 회전자 위치는 위상동기루프(PLL)를 기반으로 설계된 속도 검출기(200)로 입력되고, 상기 속도 검출기의 비례적분제어기(210)에서 설계된 모델을 통해 회전자의 각속도가 출력된다. 또한, 상기 비례적분제어기(210)에서 출력된 회전자의 각속도를 피드백 하여 조정하는 폐루프를 형성하기 위해, 적분기(220)가 구비되어 있다.The detected rotor position is input to the speed detector 200 designed based on a phase-locked loop (PLL), and the angular speed of the rotor is output through a model designed by the proportional integral controller 210 of the speed detector. In addition, an integrator 220 is provided in order to form a closed loop that adjusts by feeding back the angular velocity of the rotor output from the proportional integral controller 210.

선행기술 1에서 개시되어 있는 비선형 관측기(110)를 이용한 회전자 위치 검출기(100)와 위상동기루프를 이용한 속도 검출기(200)를 사용하여 AC 모터 제어 시스템을 설계했을 때, 엔코더나 리졸버와 같은 센서를 사용하지 않고, AC 모터의 회전자 위치와 속도를 검출할 수 있다는 장점이 있다. 보다 더 상세하게는, 센서를 사용하지 않기 때문에, 모터의 원가를 줄일 수 있고, 센서가 파손될 우려가 없기 때문에, 유지 및 보수에 대한 부담을 없앨 수 있으며, 외부 요인에 의해 신뢰성이 저하되는 점이 개선될 수 있는 장점이 있다.When designing an AC motor control system using the rotor position detector 100 using the nonlinear observer 110 and the speed detector 200 using a phase-locked loop disclosed in Prior Art 1, sensors such as encoders and resolvers It has the advantage of being able to detect the rotor position and speed of an AC motor without using a. More specifically, since the sensor is not used, the cost of the motor can be reduced, and there is no fear of damage to the sensor, so the burden of maintenance and repair can be eliminated, and the reliability deterioration due to external factors is improved. There is an advantage that can be.

그러나, 상기 비선형 관측기(110)를 이용한 회전자 위치 추정기(100)와 위상동기루프를 이용한 속도 검출기(200) 모두 내부 파라미터 변동에 민감할 수 있기 때문에, 정상 상태의 운전 조건에서는 문제없이 작동될 수 있지만, 다양한 운전 조건에서는 상태 변수에 대해서 예측 불가능한 값이 출력될 수 있다는 단점이 있다. 보다 더 자세하게는, 온도 변화와 같은 이상 조건에 의해서 고정자 저항을 비롯한 상태변수에 영향을 줄 수 있는 다른 파라미터들이 변동됨으로써, 회전자 위치 및 속도 값에 대한 신뢰성이 떨어질 수 있으며, 이는 전체 모터 제어 시스템의 성능을 저하시킬 수 있다.However, since both the rotor position estimator 100 using the nonlinear observer 110 and the speed detector 200 using the phase-locked loop can be sensitive to fluctuations in internal parameters, they can be operated without problems under normal operating conditions. However, there is a disadvantage in that unpredictable values may be output for state variables under various driving conditions. More specifically, other parameters that may affect state variables including stator resistance are changed due to abnormal conditions such as temperature change, which may reduce the reliability of the rotor position and speed values, which is the result of the overall motor control system. It can degrade the performance of.

1. 한국등록특허 제 10-1780265호(등록일 20107년 09월 14일)1. Korean Patent Registration No. 10-1780265 (Registration date September 14, 20107)

따라서 본 발명은 상기한 바와 같은 종래 기술의 문제점을 해결하기 위해 안출된 것으로, 본 발명의 목적은 이중 여자 유도 발전기(DFIG)를 속도 검출 센서를 사용하지 않고 제어하기 위해서, 복수 차수의 적분 필터(IF : Integral Filter)와 회전자 위치 보정기(RPC : Rotor Position Corrector)를 적용하여, 회전자의 위치 및 속도가 정확하게 추정될 수 있는, DFIG 센서리스 제어 시스템을 제공함에 있다.Accordingly, the present invention was conceived to solve the problems of the prior art as described above, and an object of the present invention is to control a double excitation induction generator (DFIG) without using a speed detection sensor, in order to control a multi-order integral filter ( By applying IF: Integral Filter) and Rotor Position Corrector (RPC), it is to provide a DFIG sensorless control system that can accurately estimate the position and speed of the rotor.

상기한 바와 같은 문제점을 해결하기 위한 본 발명에 의한 DFIG 센서리스 제어 시스템에 있어서, 전원(100)에 고정자 권선이 연결된 DFIG(200), 상기 DFIG(200)의 회전자 권선과 연결된 제 1 컨버터(300), 상기 제 1 컨버터(300)를 제어하는 제어부(400)를 포함하고, 상기 제어부(1000)는 상기 DFIG(200) 회전자에 흐르는 전류를 기초로 상기 DFIG(200) 회전자의 각속도를 추정하는 각속도 추정부(1100)를 포함하되, 상기 각속도 추정부(1100)는 왜곡된 전류에 의해서 추정 오차가 커지지 않도록 복수 차수의 적분 필터를 포함하는 것을 특징으로 한다.In the DFIG sensorless control system according to the present invention for solving the above problems, the DFIG 200 having a stator winding connected to the power supply 100, a first converter connected to the rotor winding of the DFIG 200 ( 300), including a control unit 400 for controlling the first converter 300, the control unit 1000 based on the current flowing through the rotor of the DFIG 200, the angular velocity of the rotor of the DFIG 200 It includes an angular velocity estimating unit 1100 to estimate, and the angular velocity estimating unit 1100 is characterized in that it includes an integral filter of a plurality of orders so as not to increase an estimation error due to the distorted current.

더 나아가, 상기 각속도 추정부(1100)는 상기 DFIG(200)의 회전자 전류 (

Figure 112020004606316-pat00001
,
Figure 112020004606316-pat00002
)을 입력받아서 정지좌표계로 좌표 변환하여 정지좌표계상의 DFIG(200) 회전자 전류 (
Figure 112020004606316-pat00003
,
Figure 112020004606316-pat00004
)를 생성하는 제 1 좌표변환부(1110), 상기 제 1 좌표변환부(1110)로부터 상기 정지좌표계상의 DFIG(200) 회전자 전류 (
Figure 112020004606316-pat00005
,
Figure 112020004606316-pat00006
)를 입력받고, 상기 복수 차수의 적분 필터(1121, 1122)를 이용하여 제1차 적분 필터링 된 회전자 전류 (
Figure 112020004606316-pat00007
,
Figure 112020004606316-pat00008
)와, 상기 제1차 적분 필터링 된 회전자 전류 (
Figure 112020004606316-pat00009
,
Figure 112020004606316-pat00010
)에 비해서 90도 위상차를 갖는 제2차 적분 필터링 된 회전자 전류 (
Figure 112020004606316-pat00011
,
Figure 112020004606316-pat00012
)를 산출하는 필터부(1120), 상기 필터부(1120)로부터 상기 제1, 2차 적분 필터링 된 회전자 전류 (
Figure 112020004606316-pat00013
,
Figure 112020004606316-pat00014
,
Figure 112020004606316-pat00015
,
Figure 112020004606316-pat00016
) 및 상기 전원(100)의 각주파수
Figure 112020004606316-pat00017
를 입력받고, 입력 각주파수
Figure 112020004606316-pat00018
을 산출하는 입력 각주파수 추정부(1130) 및 상기 입력 각주파수
Figure 112020004606316-pat00019
과 상기 전원(100)의 각 주파수
Figure 112020004606316-pat00020
를 기초로 상기 DFIG(200) 회전자의 추정 각속도
Figure 112020004606316-pat00021
을 산출하는 각속도 산출부(1140)를 포함하는 것을 특징으로 한다.Furthermore, the angular velocity estimating unit 1100 is the rotor current of the DFIG 200 (
Figure 112020004606316-pat00001
,
Figure 112020004606316-pat00002
) Is received and the coordinates are converted to the stationary coordinate system, and the DFIG(200) rotor current (
Figure 112020004606316-pat00003
,
Figure 112020004606316-pat00004
) From the first coordinate converting unit 1110, the first coordinate converting unit 1110 to the DFIG 200 on the stationary coordinate system, the rotor current (
Figure 112020004606316-pat00005
,
Figure 112020004606316-pat00006
) Is input, and the first-order integral filtered rotor current (
Figure 112020004606316-pat00007
,
Figure 112020004606316-pat00008
) And the first-order integral filtered rotor current (
Figure 112020004606316-pat00009
,
Figure 112020004606316-pat00010
), the second-order integral filtered rotor current (
Figure 112020004606316-pat00011
,
Figure 112020004606316-pat00012
) To calculate the filter unit 1120, the first and second integral filtered rotor current from the filter unit 1120 (
Figure 112020004606316-pat00013
,
Figure 112020004606316-pat00014
,
Figure 112020004606316-pat00015
,
Figure 112020004606316-pat00016
) And each frequency of the power supply 100
Figure 112020004606316-pat00017
Is input, input angular frequency
Figure 112020004606316-pat00018
The input angular frequency estimating unit 1130 to calculate and the input angular frequency
Figure 112020004606316-pat00019
And each frequency of the power supply 100
Figure 112020004606316-pat00020
Based on the estimated angular velocity of the DFIG (200) rotor
Figure 112020004606316-pat00021
It characterized in that it comprises an angular velocity calculation unit 1140 for calculating.

더 나아가, 상기 필터부(1120)는 상기 정지좌표계상의 DFIG(200) 회전자 전류 (

Figure 112020004606316-pat00022
,
Figure 112020004606316-pat00023
)와 상기 제1차 적분 필터링 된 회전자 전류 (
Figure 112020004606316-pat00024
,
Figure 112020004606316-pat00025
)의 차이에 기초한 값에서 상기 제2차 적분 필터링 된 회전자 전류 (
Figure 112020004606316-pat00026
,
Figure 112020004606316-pat00027
)를 뺀 값에 상기 전원(100)의 각 주파수
Figure 112020004606316-pat00028
를 곱한 후, 적분하여 상기 제1차 적분 필터링 된 추정 회전자 전류 (
Figure 112020004606316-pat00029
,
Figure 112020004606316-pat00030
)를 산출하는 것을 특징으로 한다.Further, the filter unit 1120 is the DFIG 200 on the stationary coordinate system, the rotor current (
Figure 112020004606316-pat00022
,
Figure 112020004606316-pat00023
) And the first-order integral filtered rotor current (
Figure 112020004606316-pat00024
,
Figure 112020004606316-pat00025
) At the value based on the difference in the second-order integral filtered rotor current (
Figure 112020004606316-pat00026
,
Figure 112020004606316-pat00027
) Minus each frequency of the power supply 100
Figure 112020004606316-pat00028
After multiplying and integrating, the first-order integral filtered estimated rotor current (
Figure 112020004606316-pat00029
,
Figure 112020004606316-pat00030
) Is characterized in that it calculates.

더 나아가, 상기 필터부(1120)는 상기 제1차 적분 필터링 된 회전자 전류 (

Figure 112020004606316-pat00031
,
Figure 112020004606316-pat00032
)를 적분하고, 상기 전원의 각주파수
Figure 112020004606316-pat00033
를 곱하여 상기 제2차 적분 필터링 된 회전자 전류 (
Figure 112020004606316-pat00034
,
Figure 112020004606316-pat00035
)를 산출하는 것을 특징으로 한다.Furthermore, the filter unit 1120 includes the first integrally filtered rotor current (
Figure 112020004606316-pat00031
,
Figure 112020004606316-pat00032
) And the angular frequency of the power supply
Figure 112020004606316-pat00033
Multiply by the second-order integral filtered rotor current (
Figure 112020004606316-pat00034
,
Figure 112020004606316-pat00035
) Is characterized in that it calculates.

더 나아가, 상기 입력 각주파수 추정부(1130)는 하기 수학식 1에 기초하여 입력 각주파수

Figure 112020004606316-pat00036
을 산출하는 것을 특징으로 한다.Further, the input angular frequency estimation unit 1130 is based on the following equation (1), the input angular frequency
Figure 112020004606316-pat00036
It characterized in that it calculates.

[수학식 1][Equation 1]

Figure 112020004606316-pat00037
Figure 112020004606316-pat00037

(이 때,

Figure 112020004606316-pat00038
: 전원의 각주파수,
Figure 112020004606316-pat00039
,
Figure 112020004606316-pat00040
: 제1차 적분 필터링 된 d축 및 q축 추정 회전자 전류,
Figure 112020004606316-pat00041
,
Figure 112020004606316-pat00042
: 제2차 적분 필터링 된 d축 및 q축 추정 회전자 전류)(At this time,
Figure 112020004606316-pat00038
: Each frequency of power,
Figure 112020004606316-pat00039
,
Figure 112020004606316-pat00040
: First-order integral filtered d-axis and q-axis estimated rotor current,
Figure 112020004606316-pat00041
,
Figure 112020004606316-pat00042
: Second-order integral filtered d-axis and q-axis estimated rotor current)

더 나아가, 상기 각속도 산출부(1130)는 상기 입력 각주파수

Figure 112020004606316-pat00043
과 상기 전원(100)의 각주파수
Figure 112020004606316-pat00044
의 차이를 입력받아 상기 DFIG(200) 회전자의 추정 각속도
Figure 112020004606316-pat00045
을 출력하는 비례 적분 제어기(1141)를 포함하는 것을 특징으로 한다.Further, the angular velocity calculation unit 1130 is the input angular frequency
Figure 112020004606316-pat00043
And each frequency of the power supply 100
Figure 112020004606316-pat00044
The estimated angular velocity of the rotor of the DFIG 200 by receiving the difference
Figure 112020004606316-pat00045
It characterized in that it comprises a proportional integral controller 1141 for outputting.

더 나아가, 상기 제 1 좌표변환부(1110)는 상기 각속도 산출부(1130)로부터 입력받는 상기 DFIG 회전자의 추정 각속도

Figure 112020004606316-pat00046
을 적분하여 얻어진 추정 회전자 위치
Figure 112020004606316-pat00047
를 이용하여 정지좌표계로 좌표 변환하는 것을 특징으로 한다.Further, the first coordinate conversion unit 1110 is the estimated angular velocity of the DFIG rotor received from the angular velocity calculation unit 1130
Figure 112020004606316-pat00046
The estimated rotor position obtained by integrating
Figure 112020004606316-pat00047
It is characterized in that the coordinates are converted to a stationary coordinate system by using.

더 나아가, 상기 DFIG(200) 회전자의 추정 각속도

Figure 112020004606316-pat00048
을 적분하여 산출되는 추정 회전자 위치
Figure 112020004606316-pat00049
와 회전자 위치 보정부(1200)에서 산출되는 회전자 위치 보정값
Figure 112020004606316-pat00050
을 기초로 상기 DFIG(200)의 회전자 추정 위치를 산출하는 제 2 회전자 위치 추정부(1310)를 더 포함하는 것을 특징으로 한다.Furthermore, the estimated angular velocity of the rotor of the DFIG 200
Figure 112020004606316-pat00048
Estimated rotor position calculated by integrating
Figure 112020004606316-pat00049
And the rotor position correction value calculated by the rotor position correction unit 1200
Figure 112020004606316-pat00050
It characterized in that it further comprises a second rotor position estimating unit 1310 for calculating the estimated rotor position of the DFIG 200 on the basis of.

더 나아가, 상기 회전자 위치 보정부(1200)는 제 1차 적분 필터링 된 회전자 전류 (

Figure 112020004606316-pat00051
,
Figure 112020004606316-pat00052
)를 입력받아 정규화 하여 제 1차 적분 필터링 및 정규화 된 회전자 전류 (
Figure 112020004606316-pat00053
,
Figure 112020004606316-pat00054
)를 출력하는 정규화부(1210), 상기 제 1차 적분 필터링 및 정규화 된 회전자 전류 (
Figure 112020004606316-pat00055
,
Figure 112020004606316-pat00056
)를 좌표 변환하여 d상의 정규화 된 회전자 전류
Figure 112020004606316-pat00057
를 산출하는 제 2 좌표변환부(1220) 및 상기 좌표 변환되어 얻어진 d상의 정규화 된 회전자 전류
Figure 112020004606316-pat00058
를 기초로 상기 회전자 위치 보정값
Figure 112020004606316-pat00059
을 산출하는 보정량 산출부(1230)를 포함하는 것을 특징으로 한다.Further, the rotor position correction unit 1200 is the first integral filtered rotor current (
Figure 112020004606316-pat00051
,
Figure 112020004606316-pat00052
) Is received and normalized to filter the first-order integral and normalized rotor current (
Figure 112020004606316-pat00053
,
Figure 112020004606316-pat00054
A normalization unit 1210 outputting ), the first-order integral filtering and normalized rotor current (
Figure 112020004606316-pat00055
,
Figure 112020004606316-pat00056
) By converting the coordinates to the normalized rotor current in d phase
Figure 112020004606316-pat00057
A second coordinate conversion unit 1220 that calculates and the normalized rotor current of the d-phase obtained by the coordinate conversion
Figure 112020004606316-pat00058
Based on the rotor position correction value
Figure 112020004606316-pat00059
It characterized in that it comprises a correction amount calculation unit 1230 for calculating.

더 나아가, 상기 보정량 산출부(1230)는 상기 d상의 정규화 된 회전자 전류

Figure 112020004606316-pat00060
와 정지좌표계상의 DFIG(200)의 d축 회전자 전류
Figure 112020004606316-pat00061
이 정규화 된
Figure 112020004606316-pat00062
과의 차이에 정지좌표계상의 DFIG(200)의 q축 회전자 전류
Figure 112020004606316-pat00063
이 정규화 된
Figure 112020004606316-pat00064
을 곱한 값을 적분한 적분값을 기초로 상기 회전자 위치 보정값
Figure 112020004606316-pat00065
을 산출하는 것을 특징으로 한다.Furthermore, the correction amount calculation unit 1230 is the normalized rotor current of the d-phase
Figure 112020004606316-pat00060
And d-axis rotor current of DFIG(200) in stationary coordinate system
Figure 112020004606316-pat00061
Is normalized
Figure 112020004606316-pat00062
Q-axis rotor current of DFIG (200) on the stationary coordinate system at the difference with
Figure 112020004606316-pat00063
Is normalized
Figure 112020004606316-pat00064
The rotor position correction value based on the integral value multiplied by
Figure 112020004606316-pat00065
It characterized in that it calculates.

더 나아가, 상기 보정량 산출부(1220)는 비례적분제어기(1222)에 상기 적분값을 입력하여 상기 회전자 위치 보정값

Figure 112020004606316-pat00066
을 산출하는 것을 특징으로 한다.Furthermore, the correction amount calculation unit 1220 inputs the integral value to the proportional integral controller 1222 to determine the rotor position correction value.
Figure 112020004606316-pat00066
It characterized in that it calculates.

더 나아가, 상기 정지좌표계상의 DFIG(200) 회전자 전류

Figure 112020004606316-pat00067
,
Figure 112020004606316-pat00068
는 하기 수학식 2 내지 수학식 5로부터 산출되는 것을 특징으로 한다.Furthermore, the DFIG 200 rotor current on the stationary coordinate system
Figure 112020004606316-pat00067
,
Figure 112020004606316-pat00068
Is characterized in that it is calculated from the following Equation 2 to Equation 5.

[수학식 2][Equation 2]

Figure 112020004606316-pat00069
Figure 112020004606316-pat00069

[수학식 3][Equation 3]

Figure 112020004606316-pat00070
Figure 112020004606316-pat00070

[수학식 4][Equation 4]

Figure 112020004606316-pat00071
Figure 112020004606316-pat00071

[수학식 5][Equation 5]

Figure 112020004606316-pat00072
Figure 112020004606316-pat00072

(이 때,

Figure 112020004606316-pat00073
: 정지좌표계상의 회전자 d축 전류,
Figure 112020004606316-pat00074
: 정지좌표계상의 회전자 q축 전류,
Figure 112020004606316-pat00075
: 정지좌표계상의 회전자 d축 쇄교 자속,
Figure 112020004606316-pat00076
: 정지좌표계상의 회전자 q축 쇄교 자속,
Figure 112020004606316-pat00077
: 자화 인덕턴스,
Figure 112020004606316-pat00078
: 고정자 인덕턴스,
Figure 112020004606316-pat00079
: 정지좌표계상의 고정자 d축 전류,
Figure 112020004606316-pat00080
: 정지좌표계상의 고정자 q축 전류,
Figure 112020004606316-pat00081
: 정지좌표계상의 고정자 d축 전압,
Figure 112020004606316-pat00082
: 정지좌표계상의 고정자 q축 전압,
Figure 112020004606316-pat00083
: 전원의 각주파수,
Figure 112020004606316-pat00084
: 고정자 저항)(At this time,
Figure 112020004606316-pat00073
: Rotor d-axis current in the stationary coordinate system,
Figure 112020004606316-pat00074
: Rotor q-axis current in the stationary coordinate system,
Figure 112020004606316-pat00075
: The d-axis linkage flux of the rotor in the stationary coordinate system,
Figure 112020004606316-pat00076
: Rotor q-axis linkage flux in the stationary coordinate system,
Figure 112020004606316-pat00077
: Magnetizing inductance,
Figure 112020004606316-pat00078
: Stator inductance,
Figure 112020004606316-pat00079
: Stator d-axis current on the stationary coordinate system,
Figure 112020004606316-pat00080
: Stator q-axis current on the stationary coordinate system,
Figure 112020004606316-pat00081
: Stator d-axis voltage on the stationary coordinate system,
Figure 112020004606316-pat00082
: Stator q-axis voltage on the stationary coordinate system,
Figure 112020004606316-pat00083
: Each frequency of power,
Figure 112020004606316-pat00084
: Stator resistance)

본 발명의 바람직한 실시예에 따른 DFIG 센서리스 제어 시스템에 의하면, 풍력 발전기에 구비된 모터의 회전자 위치 및 속도를 검출하기 위해서 엔코더 또는 레졸버와 같은 센서를 사용하지 않기 때문에, 모터의 부피를 줄일 수 있고, 더 나아가서는 전체 터빈 시스템의 부피와 중량을 감소시킬 수 있다는 장점이 있다.According to the DFIG sensorless control system according to a preferred embodiment of the present invention, since a sensor such as an encoder or a resolver is not used to detect the rotor position and speed of a motor provided in a wind power generator, the volume of the motor is reduced. And, furthermore, it has the advantage of reducing the volume and weight of the entire turbine system.

보다 더 상세하게는, 종래에는 모터의 회전자 위치 및 속도를 검출하기 위해서 센서를 사용했기 때문에, 모터의 부피에 속도 검출 센서의 부피가 추가되어, 기어 박스 등 터빈 시스템 내부에서 각 부품을 배치하는 데 어려움이 있었고, 전체 부피와 중량 또한 증가한다는 단점이 있는 반면, 본 발명에서는 속도 검출 센서를 사용하지 않기 때문에 공간을 활용하여 많은 부품을 각각 배치하는 것이 더욱 수월하고, 이로 인해 부피와 중량을 감소시킬 수 있다는 장점이 있다.More specifically, since sensors were used to detect the rotor position and speed of a motor in the past, the volume of the speed detection sensor is added to the volume of the motor, so that each component is placed inside the turbine system such as a gearbox. While there is a disadvantage in that the overall volume and weight are also increased, in the present invention, since the speed detection sensor is not used, it is easier to place many parts separately using space, thereby reducing the volume and weight. There is an advantage of being able to do it.

더불어, 본 발명의 바람직한 실시예에 따른 DFIG 센서리스 제어 시스템에 의하면, 가격이 비싼 속도 검출 센서를 사용하지 않기 때문에, 종래 보다 전체 시스템을 구성하기 위한 원가를 절감할 수 있다는 장점이 있다.In addition, according to the DFIG sensorless control system according to the preferred embodiment of the present invention, since the expensive speed detection sensor is not used, there is an advantage in that the cost for configuring the entire system can be reduced compared to the prior art.

또한, 본 발명의 바람직한 실시예에 따른 DFIG 센서리스 제어 시스템에 의하면, 복수의 적분을 수행하는 필터부(1120)가 적용됨에 따라, 측정된 전류에서 리플(Ripple)이 존재하더라도 회전자의 속도가 정확하게 추정될 수 있다. 그리고 회전자 위치 보정부(RPC, 1200)를 적용함에 따라 추정된 회전자의 위치에 대한 오차를 크게 줄임으로써, 회전자의 위치가 보다 더 정확하게 확인될 수 있다는 장점이 있다. In addition, according to the DFIG sensorless control system according to a preferred embodiment of the present invention, as the filter unit 1120 performing a plurality of integrations is applied, the speed of the rotor is increased even if there is a ripple in the measured current. Can be estimated accurately. In addition, by applying the rotor position correction unit (RPC, 1200), there is an advantage that the position of the rotor can be more accurately identified by greatly reducing an error in the estimated position of the rotor.

더불어, 상기 필터부(1120)와 회전자 위치 보정부(1200)를 통해 풍력 발전 시스템을 제어함으로써, 외부 조건의 변동 및 이상 조건에도 회전자의 위치 및 속도가 강인하고, 정확하게 검출될 수 있다는 장점이 있다.In addition, by controlling the wind power generation system through the filter unit 1120 and the rotor position correction unit 1200, the position and speed of the rotor are robust and can be accurately detected even under fluctuations in external conditions and abnormal conditions. There is this.

도 1은 선행기술 1의 일실시예에 따른 전체 제어 시스템의 블록도이다.
도 2는 DFIG 풍력 발전 시스템의 일반적인 전체 블록도이다.
도 3의 (A)는 표준 적분 필터의 블록도 이고, 도 3의 (B)는 수정된 적분 필터의 블록도 이다.
도 4은 본 발명의 바람직한 실시예에 따른 DFIG 센서리스 제어 시스템의 적분 필터 기반 속도 추정기의 블록도이다.
도 5는 본 발명의 바람직한 실시예에 따른 DFIG 센서리스 제어 시스템의 필터부의 블록도이다.
도 6은 본 발명의 바람직한 실시예에 따른 DFIG 센서리스 제어 시스템의 회전자 위치 보정부의 블록도이다.
도 7은 본 발명의 바람직한 실시예에 따른 DFIG 센서리스 제어 시스템의 회전자 위치 보정값

Figure 112020004606316-pat00085
과 입력 각주파수
Figure 112020004606316-pat00086
의 조건에 따른 상기 신호
Figure 112020004606316-pat00087
의 변화 그래프이다.
도 8은 본 발명의 바람직한 실시예에 따른 DFIG 센서리스 제어 시스템의 슬립각 추정에 대한 블록도이다.
도 9는 본 발명의 바람직한 실시예에 따른 DFIG 센서리스 제어 시스템의 동기화 과정에서의 슬립 추정 방법에 대한 블록도이다.
도 10 내지 도 12는 본 발명의 바람직한 실시예에 따른 DFIG 센서리스 제어 시스템의 실험 결과 그래프이다.1 is a block diagram of an overall control system according to an embodiment of Prior Art 1;
2 is a general block diagram of a DFIG wind power generation system.
3(A) is a block diagram of a standard integral filter, and FIG. 3(B) is a block diagram of a modified integral filter.
4 is a block diagram of an integral filter-based velocity estimator of a DFIG sensorless control system according to a preferred embodiment of the present invention.
5 is a block diagram of a filter unit of a DFIG sensorless control system according to a preferred embodiment of the present invention.
6 is a block diagram of a rotor position correction unit of a DFIG sensorless control system according to a preferred embodiment of the present invention.
7 is a rotor position correction value of the DFIG sensorless control system according to a preferred embodiment of the present invention
Figure 112020004606316-pat00085
And input angular frequency
Figure 112020004606316-pat00086
The signal according to the condition of
Figure 112020004606316-pat00087
Is a graph of change.
8 is a block diagram for estimating a slip angle of a DFIG sensorless control system according to a preferred embodiment of the present invention.
9 is a block diagram of a slip estimation method in a synchronization process of a DFIG sensorless control system according to a preferred embodiment of the present invention.
10 to 12 are graphs of experimental results of a DFIG sensorless control system according to a preferred embodiment of the present invention.

이하, 상기한 바와 같은 구성을 가지는 본 발명의 바람직한 실시예에 따른 컨버터를 위한 센서 고장 진단 시스템의 도시된 도면들을 참조하여 상세히 설명한다. 다음에 소개되는 도면들은 당업자에게 본 발명의 사상이 충분히 전달될 수 있도록 하기 위해 실시예로서 제공되는 것이다. 따라서, 본 발명은 이하 제시되는 도면들에 한정되지 않고 다른 형태로 구체화될 수도 있다. 또한, 명세서 전반에 걸쳐서 동일한 참조번호들은 동일한 구성요소들을 나타낸다.Hereinafter, a sensor failure diagnosis system for a converter according to a preferred embodiment of the present invention having the above-described configuration will be described in detail with reference to the drawings. The drawings introduced below are provided as examples in order to sufficiently convey the spirit of the present invention to those skilled in the art. Accordingly, the present invention is not limited to the drawings presented below and may be embodied in other forms. In addition, the same reference numbers throughout the specification indicate the same elements.

이 때, 사용되는 기술 용어 및 과학 용어에 있어서 다른 정의가 없다면, 이 발명이 속하는 기술 분야에서 통상의 지식을 가진 자가 통상적으로 이해하고 있는 의미를 가지며, 하기의 설명 및 첨부 도면에서 본 발명의 요지를 불필요하게 흐릴 수 있는 공지 기능 및 구성에 대한 설명은 생략한다.In this case, unless there are other definitions in the technical terms and scientific terms used, they have the meanings commonly understood by those of ordinary skill in the art to which this invention belongs, and the gist of the present invention in the following description and accompanying drawings Descriptions of known functions and configurations that may unnecessarily obscure are omitted.

더불어, 시스템은 필요한 기능을 수행하기 위하여 조직화되고 규칙적으로 상호 작용하는 장치, 기구 및 수단 등을 포함하는 구성 요소들의 집합을 의미한다.In addition, a system refers to a set of components including devices, devices, and means that are organized and regularly interact to perform a required function.

DFIG 풍력 발전 시스템 모델DFIG wind power system model

도 2는 일반적인 DFIG 풍력 발전 시스템의 전체 블록도를 도시하고 있다. 도 2에서 도시된 바와 같이, DFIG(200) 풍력 발전 시스템에서 제 1 컨버터(GSC : Grid Side Converter, 400)와 제 2 컨버터(RSC : Rotor Side Converter, 300)는 각각의 회로에 대응하는 내부 전류, 즉 제 1 컨버터(400) 측의 3상 전류(

Figure 112020004606316-pat00088
,
Figure 112020004606316-pat00089
,
Figure 112020004606316-pat00090
), 제 2 컨버터(300) 측의 3상 전류 (
Figure 112020004606316-pat00091
,
Figure 112020004606316-pat00092
,
Figure 112020004606316-pat00093
)의 제어 루프와 외부 DC 링크 전압
Figure 112020004606316-pat00094
및 고정자 유효 전력
Figure 112020004606316-pat00095
, 무효 전력
Figure 112020004606316-pat00096
에 대한 제어 루프에 있어서 하나의 출력이 다음단의 입력이 되도록 하는 캐스케이스(Cascade) 제어 구조를 가진다. Figure 2 shows the overall block diagram of a typical DFIG wind power generation system. As shown in FIG. 2, in the DFIG 200 wind power generation system, the first converter (GSC: Grid Side Converter, 400) and the second converter (RSC: Rotor Side Converter, 300) have internal currents corresponding to respective circuits. That is, the three-phase current on the side of the first converter 400 (
Figure 112020004606316-pat00088
,
Figure 112020004606316-pat00089
,
Figure 112020004606316-pat00090
), the three-phase current of the second converter 300 side (
Figure 112020004606316-pat00091
,
Figure 112020004606316-pat00092
,
Figure 112020004606316-pat00093
) Of the control loop and external DC link voltage
Figure 112020004606316-pat00094
And stator active power
Figure 112020004606316-pat00095
, Reactive power
Figure 112020004606316-pat00096
It has a cascade control structure in which one output becomes the input of the next stage in the control loop for.

상기 DFIG(200) 풍력 발전 시스템의 모델링에 따른 전압 방정식은 하기 수학식 6 내지 수학식 9와 같이 표현될 수 있다.The voltage equation according to the modeling of the wind power generation system of the DFIG 200 may be expressed as Equations 6 to 9 below.

[수학식 6][Equation 6]

Figure 112020004606316-pat00097
Figure 112020004606316-pat00097

[수학식 7][Equation 7]

Figure 112020004606316-pat00098
Figure 112020004606316-pat00098

[수학식 8][Equation 8]

Figure 112020004606316-pat00099
Figure 112020004606316-pat00099

[수학식 9][Equation 9]

Figure 112020004606316-pat00100
Figure 112020004606316-pat00100

(이 때,

Figure 112020004606316-pat00101
: d축 고정자 전압,
Figure 112020004606316-pat00102
: q축 고정자 전압,
Figure 112020004606316-pat00103
: d축 회전자 전압,
Figure 112020004606316-pat00104
: q축 회전자 전압,
Figure 112020004606316-pat00105
: d축 고정자 전류,
Figure 112020004606316-pat00106
: q축 고정자 전류,
Figure 112020004606316-pat00107
: d축 회전자 전류,
Figure 112020004606316-pat00108
: q축 회전자 전류,
Figure 112020004606316-pat00109
: 고정자 저항,
Figure 112020004606316-pat00110
: 회전자 저항,
Figure 112020004606316-pat00111
: 고정자 각주파수,
Figure 112020004606316-pat00112
: 회전자 각속도,
Figure 112020004606316-pat00113
: d축 고정자 쇄교 자속,
Figure 112020004606316-pat00114
: q축 고정자 쇄교자속,
Figure 112020004606316-pat00115
: d축 회전자 쇄교 자속,
Figure 112020004606316-pat00116
: q축 회전자 쇄교 자속,
Figure 112020004606316-pat00117
: d축 고정자 쇄교 자속의 시간 t에 대한 미분,
Figure 112020004606316-pat00118
: q축 고정자 쇄교 자속의 시간 t에 대한 미분,
Figure 112020004606316-pat00119
: d축 회전자 쇄교 자속의 시간 t에 대한 미분,
Figure 112020004606316-pat00120
: q축 회전자 쇄교 자속의 시간 t에 대한 미분)(At this time,
Figure 112020004606316-pat00101
: d-axis stator voltage,
Figure 112020004606316-pat00102
: q-axis stator voltage,
Figure 112020004606316-pat00103
: d-axis rotor voltage,
Figure 112020004606316-pat00104
: q-axis rotor voltage,
Figure 112020004606316-pat00105
: d-axis stator current,
Figure 112020004606316-pat00106
: q-axis stator current,
Figure 112020004606316-pat00107
: d-axis rotor current,
Figure 112020004606316-pat00108
: q-axis rotor current,
Figure 112020004606316-pat00109
: Stator resistance,
Figure 112020004606316-pat00110
: Rotor resistance,
Figure 112020004606316-pat00111
: Stator angular frequency,
Figure 112020004606316-pat00112
: Rotor angular velocity,
Figure 112020004606316-pat00113
: d-axis stator flux linkage,
Figure 112020004606316-pat00114
: q-axis stator flux linkage,
Figure 112020004606316-pat00115
: d-axis rotor linkage flux,
Figure 112020004606316-pat00116
: q-axis rotor linkage flux,
Figure 112020004606316-pat00117
: Differentiation of the d-axis stator flux linkage with respect to time t,
Figure 112020004606316-pat00118
: The derivative of the q-axis stator flux linkage with respect to time t,
Figure 112020004606316-pat00119
: The derivative of the d-axis rotor linkage flux with respect to time t,
Figure 112020004606316-pat00120
: The derivative of the q-axis rotor linkage flux with respect to time t)

상기 수학식 6 내지 수학식 9 에서 고정자와 회전자의 쇄교 자속은 하기 수학식 10 내지 수학식 13과 같이 표현될 수 있다.In Equations 6 to 9, the flux linkage between the stator and the rotor may be expressed as Equations 10 to 13 below.

[수학식 10][Equation 10]

Figure 112020004606316-pat00121
Figure 112020004606316-pat00121

[수학식 11][Equation 11]

Figure 112020004606316-pat00122
Figure 112020004606316-pat00122

[수학식 12][Equation 12]

Figure 112020004606316-pat00123
Figure 112020004606316-pat00123

[수학식 13][Equation 13]

Figure 112020004606316-pat00124
Figure 112020004606316-pat00124

(이 때,

Figure 112020004606316-pat00125
: 자화 인덕턴스,
Figure 112020004606316-pat00126
: 고정자 인덕턴스,
Figure 112020004606316-pat00127
: 회전자 자기 인덕턴스,
Figure 112020004606316-pat00128
,
Figure 112020004606316-pat00129
: 고정자 누설 인덕턴스,
Figure 112020004606316-pat00130
: 회전자 누설 인덕턴스)(At this time,
Figure 112020004606316-pat00125
: Magnetizing inductance,
Figure 112020004606316-pat00126
: Stator inductance,
Figure 112020004606316-pat00127
: Rotor magnetic inductance,
Figure 112020004606316-pat00128
,
Figure 112020004606316-pat00129
: Stator leakage inductance,
Figure 112020004606316-pat00130
: Rotor leakage inductance)

상기 수학식 10 내지 수학식 13를 이용하여 고정자 유효 전력

Figure 112020004606316-pat00131
및 무효 전력
Figure 112020004606316-pat00132
는 하기 수학식 14 및 수학식 15와 같이 표현될 수 있다.Stator active power using Equations 10 to 13
Figure 112020004606316-pat00131
And reactive power
Figure 112020004606316-pat00132
May be expressed as in Equation 14 and Equation 15 below.

[수학식 14][Equation 14]

Figure 112020004606316-pat00133
Figure 112020004606316-pat00133

[수학식 15][Equation 15]

Figure 112020004606316-pat00134
Figure 112020004606316-pat00134

(이 때,

Figure 112020004606316-pat00135
: 고정자 자화 전류)(At this time,
Figure 112020004606316-pat00135
: Stator magnetizing current)

상기 DFIG(200) 풍력 발전 시스템 모델을 기초로 한 본 발명의 바람직한 실시예에 따른 DFIG 센서리스 제어 시스템의 각 요소에 대해 이하에서 상세히 설명한다.Each element of the DFIG sensorless control system according to a preferred embodiment of the present invention based on the DFIG 200 wind power generation system model will be described in detail below.

복수 차수의 적분 필터(Integral Filter)Multiple order integral filter

도 3의 (A)는 표준 복수 차수의 적분 필터의 블록도에 대해서 도시하고 있고, 도 3의 (B)는 수정된 복수 차수의 적분 필터의 블록도에 대해서 도시하고 있다. 도 3의 (A)에 도시된 바와 같이, 상기 복수 차수의 적분 필터는 점선의 박스가 그려진 부분이고, 적분 필터의 전달 함수는 하기 수학식 16 및 수학식 17과 같다.Fig. 3A shows a block diagram of a standard multi-order integration filter, and Fig. 3B shows a block diagram of a modified multi-order integration filter. As shown in (A) of FIG. 3, the multi-order integral filter is a portion drawn with a dotted box, and the transfer function of the integral filter is shown in Equations 16 and 17 below.

[수학식 16][Equation 16]

Figure 112020004606316-pat00136
Figure 112020004606316-pat00136

[수학식 17][Equation 17]

Figure 112020004606316-pat00137
Figure 112020004606316-pat00137

(이 때,

Figure 112020004606316-pat00138
: 입력 신호,
Figure 112020004606316-pat00139
,
Figure 112020004606316-pat00140
: 출력 신호,
Figure 112020004606316-pat00141
Figure 112020004606316-pat00142
는 직교 관계,
Figure 112020004606316-pat00143
: 라플라스 변환 상수,
Figure 112020004606316-pat00144
: 튜닝 주파수,
Figure 112020004606316-pat00145
: 감쇠 계수)(At this time,
Figure 112020004606316-pat00138
: Input signal,
Figure 112020004606316-pat00139
,
Figure 112020004606316-pat00140
: Output signal,
Figure 112020004606316-pat00141
and
Figure 112020004606316-pat00142
Is an orthogonal relationship,
Figure 112020004606316-pat00143
: Laplace transform constant,
Figure 112020004606316-pat00144
: Tuning frequency,
Figure 112020004606316-pat00145
: Attenuation factor)

만약, 입력 신호가

Figure 112020004606316-pat00146
라고 할 때, 두 출력 신호
Figure 112020004606316-pat00147
,
Figure 112020004606316-pat00148
은 하기 수학식 18 및 수학식 19와 같이 표현될 수 있다.If the input signal is
Figure 112020004606316-pat00146
When said, the two output signals
Figure 112020004606316-pat00147
,
Figure 112020004606316-pat00148
May be expressed as in Equation 18 and Equation 19 below.

[수학식 18][Equation 18]

Figure 112020004606316-pat00149
Figure 112020004606316-pat00149

[수학식 19][Equation 19]

Figure 112020004606316-pat00150
Figure 112020004606316-pat00150

(이 때,

Figure 112020004606316-pat00151
: 상수,
Figure 112020004606316-pat00152
: 주파수,
Figure 112020004606316-pat00153
: 시간,
Figure 112020004606316-pat00154
: 위상각,
Figure 112020004606316-pat00155
)(At this time,
Figure 112020004606316-pat00151
: a constant,
Figure 112020004606316-pat00152
: frequency,
Figure 112020004606316-pat00153
: time,
Figure 112020004606316-pat00154
: Phase angle,
Figure 112020004606316-pat00155
)

상기 복수 차수의 적분 필터의 적용은 응답 시간이 비교적 길기 때문에, 도 3의 (B)에 도시되어 있는 수정된 복수 차수의 적분 필터가 사용될 수 있다. 상기 수정된 복수 차수의 적분 필터에서는 이득 계수

Figure 112020004606316-pat00156
를 사용하고,
Figure 112020004606316-pat00157
가 더해지는 이득 회로가 추가되는데, 이는 응답 시간을 줄이면서 출력 신호에 대한 추정 속도를 향상 시킬 수 있다. 그러나, 고조파가 발생 할 수 있는 단점이 생기게 되고, 이러한 단점으로 인해 복수 차수의 적분 필터의 응답 시간을 줄이면서, 고조파가 제거 될 수 있는 절충점이 고려되어야 한다.Since the application of the multiple-order integration filter has a relatively long response time, the modified multiple-order integration filter shown in Fig. 3B can be used. In the modified multi-order integral filter, the gain factor
Figure 112020004606316-pat00156
And use
Figure 112020004606316-pat00157
A gain circuit to which is added is added, which can improve the estimation speed for the output signal while reducing the response time. However, there is a disadvantage in that harmonics may occur, and due to this disadvantage, a trade-off in which harmonics can be eliminated while reducing the response time of the multi-order integral filter must be considered.

상기 복수 차수의 적분 필터에 대해서는 이하에서 개시되는 본 발명의 바람직한 실시예에 따른 DFIG 센서리스 제어 시스템의 각속도 추정부(1100)에서 자세하게 설명한다.The multi-order integral filter will be described in detail in the angular velocity estimation unit 1100 of the DFIG sensorless control system according to a preferred embodiment of the present invention disclosed below.

복수 차수의 적분 필터 기반 각속도 추정부Angular velocity estimation unit based on multi-order integral filter

도 4는 본 발명의 바람직한 실시예에 따른 DFIG 센서리스 제어 시스템의 적분 필터 기반 각속도 추정부(1100)가 도시되어 있다. 4 shows an integral filter-based angular velocity estimation unit 1100 of a DFIG sensorless control system according to a preferred embodiment of the present invention.

도 4에 도시된 바와 같이, 상기 DFIG(200)의 회전자 전류 (

Figure 112020004606316-pat00158
,
Figure 112020004606316-pat00159
)와 기존에 추정된 회전자 속도
Figure 112020004606316-pat00160
을 적분하여 얻어진 추정 회전자 위치
Figure 112020004606316-pat00161
는 상기 속도 추정기의 제 1 좌표변환부(1110)로 입력되고, 상기 제 1 좌표변환부(1110)에서 정지좌표계상의 DFIG(200)의 회전자 전류 (
Figure 112020004606316-pat00162
,
Figure 112020004606316-pat00163
)가 출력된다. 4, the rotor current of the DFIG 200 (
Figure 112020004606316-pat00158
,
Figure 112020004606316-pat00159
) And the previously estimated rotor speed
Figure 112020004606316-pat00160
The estimated rotor position obtained by integrating
Figure 112020004606316-pat00161
Is input to the first coordinate conversion unit 1110 of the speed estimator, and the rotor current of the DFIG 200 on the stationary coordinate system in the first coordinate conversion unit 1110 (
Figure 112020004606316-pat00162
,
Figure 112020004606316-pat00163
) Is displayed.

상기 출력된 정지좌표계상의 DFIG(200) 회전자 전류 (

Figure 112020004606316-pat00164
,
Figure 112020004606316-pat00165
)와 전원(100) 주파수
Figure 112020004606316-pat00166
은 상기 복수 차수의 적분 필터를 기반으로 설계된 필터부(1120)로 입력되어 제1차 적분 필터링 된 회전자 전류 (
Figure 112020004606316-pat00167
,
Figure 112020004606316-pat00168
)와 상기 제1차 적분 필터링 된 회전자 전류 (
Figure 112020004606316-pat00169
,
Figure 112020004606316-pat00170
)에 비해서 90도 위상차를 갖는 제2차 적분 필터링 된 회전자 전류 (
Figure 112020004606316-pat00171
,
Figure 112020004606316-pat00172
)가 출력된다.DFIG (200) rotor current on the output stationary coordinate system (
Figure 112020004606316-pat00164
,
Figure 112020004606316-pat00165
) And power (100) frequency
Figure 112020004606316-pat00166
Is input to the filter unit 1120 designed based on the multi-order integral filter and the first-order integral filtered rotor current (
Figure 112020004606316-pat00167
,
Figure 112020004606316-pat00168
) And the first-order integral filtered rotor current (
Figure 112020004606316-pat00169
,
Figure 112020004606316-pat00170
), the second-order integral filtered rotor current (
Figure 112020004606316-pat00171
,
Figure 112020004606316-pat00172
) Is displayed.

상기 제 1, 2차 적분 필터링 된 회전자 전류 (

Figure 112020004606316-pat00173
,
Figure 112020004606316-pat00174
,
Figure 112020004606316-pat00175
,
Figure 112020004606316-pat00176
)는 입력 각주파수 추정부(Input frequency estimator, 1130)로 입력되고, 하기 수학식 1에 의해 입력 각주파수
Figure 112020004606316-pat00177
이 연산된다.The first and second integral filtered rotor currents (
Figure 112020004606316-pat00173
,
Figure 112020004606316-pat00174
,
Figure 112020004606316-pat00175
,
Figure 112020004606316-pat00176
) Is input to an input angular frequency estimator (1130), and input angular frequency by Equation 1 below
Figure 112020004606316-pat00177
Is calculated.

[수학식 1][Equation 1]

Figure 112020004606316-pat00178
Figure 112020004606316-pat00178

상기 연산된 입력 각주파수

Figure 112020004606316-pat00179
은 각속도 산출부(1140)로 입력되고, 상기 전원(100) 주파수
Figure 112020004606316-pat00180
에서 상기 입력 각주파수
Figure 112020004606316-pat00181
을 뺀 값이 상기 각속도 산출부(1140)의 비례적분제어기(1141)로 입력된다. 상기 비례적분제어기(1141)에 의해 추정 회전자 각속도
Figure 112020004606316-pat00182
이 출력되고, 폐루프의 피드백 회로로 구비되는 제 1 회전자 위치 추정부(1150)로 상기 추정 회전자 각속도
Figure 112020004606316-pat00183
이 입력되어 추정 회전자 위치
Figure 112020004606316-pat00184
가 출력된다.Input angular frequency calculated above
Figure 112020004606316-pat00179
Is input to the angular velocity calculation unit 1140, and the frequency of the power supply 100
Figure 112020004606316-pat00180
At the input angular frequency
Figure 112020004606316-pat00181
The value minus is input to the proportional integral controller 1141 of the angular velocity calculator 1140. Estimated rotor angular velocity by the proportional integral controller 1141
Figure 112020004606316-pat00182
Is output, and the estimated rotor angular velocity to the first rotor position estimating unit 1150 provided as a feedback circuit of the closed loop
Figure 112020004606316-pat00183
Is inputted and estimated rotor position
Figure 112020004606316-pat00184
Is displayed.

도 5는 본 발명의 바람직한 실시예에 따른 DFIG 센서리스 제어 시스템의 복수 차수의 적분 필터 기반 각속도 추정부에서 필터부(1120)를 도시하고 있다. 5 illustrates a filter unit 1120 in an angular velocity estimation unit based on a multi-order integral filter of a DFIG sensorless control system according to a preferred embodiment of the present invention.

도 5에 도시된 바와 같이, 필터부(1120)는 d축 추정 회전자 전류에 대하여 필터링 하는 제 1 적분 필터부(1121), q축 추정 회전자 전류에 대하여 필터링 하는 제 2 적분 필터부(1122)를 포함하며, 상기 필터부(1120)에는 도 3의 (A)에 도시되어 있는 복수 차수의 적분 필터가 포함되어 있다.As shown in FIG. 5, the filter unit 1120 includes a first integral filter unit 1121 that filters the d-axis estimated rotor current, and a second integral filter unit 1122 that filters the q-axis estimated rotor current. ), and the filter unit 1120 includes an integral filter of a plurality of orders shown in (A) of FIG. 3.

제 1 적분 필터부(1121)와 제 2 적분 필터부(1122)는 제1차 적분 필터링 된 회전자 전류 (

Figure 112020004606316-pat00185
,
Figure 112020004606316-pat00186
)와 상기 제1차 적분 필터링 된 회전자 전류 (
Figure 112020004606316-pat00187
,
Figure 112020004606316-pat00188
)에 비해서 90도 위상차를 갖는 제2차 적분 필터링 된 회전자 전류 (
Figure 112020004606316-pat00189
,
Figure 112020004606316-pat00190
)을 출력하게 된다.The first integral filter unit 1121 and the second integral filter unit 1122 have a first-order integral filtered rotor current (
Figure 112020004606316-pat00185
,
Figure 112020004606316-pat00186
) And the first-order integral filtered rotor current (
Figure 112020004606316-pat00187
,
Figure 112020004606316-pat00188
), the second-order integral filtered rotor current (
Figure 112020004606316-pat00189
,
Figure 112020004606316-pat00190
) Is displayed.

또한, 상기 감쇠계수

Figure 112020004606316-pat00191
의 값에 의해서 상기 제 1 적분 필터부(1121) 및 제 2 적분 필터부(1122)의 정착 시간(리플이 반복되면서 수렴되는 값의 ±1% 또는 ±2.5% 이내로 진입하는 시간)이 결정될 수 있다.In addition, the attenuation coefficient
Figure 112020004606316-pat00191
The settling time of the first integral filter unit 1121 and the second integral filter unit 1122 (the time to enter within ±1% or ±2.5% of the value converged while the ripple is repeated) may be determined by the value of .

회전자 위치 보정부Rotor position correction unit

도 6은 본 발명의 바람직한 실시예에 따른 DFIG 센서리스 제어 시스템의 회전자 위치 보정부(1200)의 블록도가 도시되어 있다. 6 is a block diagram of a rotor position correction unit 1200 of a DFIG sensorless control system according to a preferred embodiment of the present invention.

상기 회전자의 위치

Figure 112020004606316-pat00192
는 상기 각속도 추정부(1100)에서 개시한 바와 같이
Figure 112020004606316-pat00193
이 적분되는 것에 의해 추정될 수 있지만, 속도 추정 오차와 초기 회전자 위치 정보 부족 등 여러 요인들에 의해 잘못된 추정값이 생성될 수 있다. 이러한 문제점을 보완하기 위해 회전자 위치 보정값
Figure 112020004606316-pat00194
이 상기 추정 회전자 위치
Figure 112020004606316-pat00195
에 보상됨으로써, 상기 회전자의 위치가 정확하게 추정될 수 있다.The position of the rotor
Figure 112020004606316-pat00192
Is, as started in the angular velocity estimation unit 1100
Figure 112020004606316-pat00193
Although this can be estimated by integration, an incorrect estimate may be generated due to various factors such as a speed estimation error and a lack of initial rotor position information. To compensate for this problem, the rotor position correction value
Figure 112020004606316-pat00194
This said estimated rotor position
Figure 112020004606316-pat00195
By being compensated for, the position of the rotor can be accurately estimated.

도 6에 도시된 바와 같이, 회전자 위치 보정부(1200)는 상기 각속도 추정부(1100)의 상기 필터부(1120)로부터 제1차 적분 필터링 된 회전자 전류 (

Figure 112020004606316-pat00196
,
Figure 112020004606316-pat00197
)를 입력받아 정규화부(1210)에서 상기 제 1차 적분 필터링 된 회전자 전류를 정규화 하여 정규화 된 회전자 전류 (
Figure 112020004606316-pat00198
,
Figure 112020004606316-pat00199
)가 출력되고, 기존에 출력된 회전자 위치 보정값
Figure 112020004606316-pat00200
을 피드백 받아 좌표가 변환된 d상의 정규화 된 회전자 전류
Figure 112020004606316-pat00201
를 출력하는 제 2 좌표 변환부(1220) 및 상기 출력된 d상의 추정 회전자 전류
Figure 112020004606316-pat00202
를 입력받아 상기 회전자 위치 보정값
Figure 112020004606316-pat00203
을 출력하는 보정량 산출부(1230)가 포함된다.As shown in FIG. 6, the rotor position correction unit 1200 includes a first-order integrally filtered rotor current from the filter unit 1120 of the angular velocity estimating unit 1100 (
Figure 112020004606316-pat00196
,
Figure 112020004606316-pat00197
) Is received, and the first-order integral filtered rotor current is normalized in the normalization unit 1210, and the normalized rotor current (
Figure 112020004606316-pat00198
,
Figure 112020004606316-pat00199
) Is output, and the previously outputted rotor position correction value
Figure 112020004606316-pat00200
The normalized rotor current of the d-phase whose coordinates are converted by receiving the feedback
Figure 112020004606316-pat00201
The second coordinate converting unit 1220 outputting the output and the estimated rotor current of the d-phase
Figure 112020004606316-pat00202
Receive the input and the rotor position correction value
Figure 112020004606316-pat00203
A correction amount calculation unit 1230 for outputting is included.

정지좌표계상의 DFIG(200) 회전자 전류 및 쇄교 자속은 하기 수학식 2 내지 수학식 5로 표현될 수 있고, 상기 정지좌표계상의 회전자 전류는 상기 보정량 산출부(1230)에서 상기 회전자 위치 보정값

Figure 112020004606316-pat00204
이 출력되기 위한 연산에 사용된다.The rotor current and the flux linkage of the DFIG 200 on the stationary coordinate system can be expressed by the following equations 2 to 5, and the rotor current on the stationary coordinate system is the rotor position correction value in the correction amount calculation unit 1230
Figure 112020004606316-pat00204
Is used in the operation to be output.

[수학식 2][Equation 2]

Figure 112020004606316-pat00205
Figure 112020004606316-pat00205

[수학식 3][Equation 3]

Figure 112020004606316-pat00206
Figure 112020004606316-pat00206

[수학식 4][Equation 4]

Figure 112020004606316-pat00207
Figure 112020004606316-pat00207

[수학식 5][Equation 5]

Figure 112020004606316-pat00208
Figure 112020004606316-pat00208

(이 때,

Figure 112020004606316-pat00209
: 정지좌표계상의 회전자 d축 전류,
Figure 112020004606316-pat00210
: 정지좌표계상의 회전자 q축 전류,
Figure 112020004606316-pat00211
: 정지좌표계상의 회전자 d축 쇄교 자속,
Figure 112020004606316-pat00212
: 정지좌표계상의 회전자 q축 쇄교 자속,
Figure 112020004606316-pat00213
: 자화 인덕턴스,
Figure 112020004606316-pat00214
: 고정자 인덕턴스,
Figure 112020004606316-pat00215
: 정지좌표계상의 고정자 d축 전류,
Figure 112020004606316-pat00216
: 정지좌표계상의 고정자 q축 전류,
Figure 112020004606316-pat00217
: 정지좌표계상의 고정자 d축 전압,
Figure 112020004606316-pat00218
: 정지좌표계상의 고정자 q축 전압,
Figure 112020004606316-pat00219
: 전원의 각주파수,
Figure 112020004606316-pat00220
: 고정자 저항)(At this time,
Figure 112020004606316-pat00209
: Rotor d-axis current in the stationary coordinate system,
Figure 112020004606316-pat00210
: Rotor q-axis current in the stationary coordinate system,
Figure 112020004606316-pat00211
: The d-axis linkage flux of the rotor in the stationary coordinate system,
Figure 112020004606316-pat00212
: Rotor q-axis linkage flux in the stationary coordinate system,
Figure 112020004606316-pat00213
: Magnetizing inductance,
Figure 112020004606316-pat00214
: Stator inductance,
Figure 112020004606316-pat00215
: Stator d-axis current on the stationary coordinate system,
Figure 112020004606316-pat00216
: Stator q-axis current on the stationary coordinate system,
Figure 112020004606316-pat00217
: Stator d-axis voltage on the stationary coordinate system,
Figure 112020004606316-pat00218
: Stator q-axis voltage on the stationary coordinate system,
Figure 112020004606316-pat00219
: Each frequency of power,
Figure 112020004606316-pat00220
: Stator resistance)

상기 정지좌표계상의 회전자 전류 및 쇄교 자속은 고정자 인덕턴스

Figure 112020004606316-pat00221
에 의해 오류 값이 출력될 수 있으므로, 상기 회전자 전류에 대해 정규화 과정이 필요하다. 상기 제1차 적분 필터링 된 회전자 전류 (
Figure 112020004606316-pat00222
,
Figure 112020004606316-pat00223
)는 정규화부(1210)에서 정규화 된 회전자 전류(
Figure 112020004606316-pat00224
,
Figure 112020004606316-pat00225
)로 변환될 수 있고, 정지좌표계상의 DFIG(200) 회전자 전류 (
Figure 112020004606316-pat00226
,
Figure 112020004606316-pat00227
) 또한 정규화 되어 보정량 산출부로 입력된다.The rotor current and the flux linkage in the stationary coordinate system are the stator inductance
Figure 112020004606316-pat00221
Since an error value can be output by, a normalization process is required for the rotor current. The first order integral filtered rotor current (
Figure 112020004606316-pat00222
,
Figure 112020004606316-pat00223
) Is the rotor current normalized by the normalization unit 1210 (
Figure 112020004606316-pat00224
,
Figure 112020004606316-pat00225
) Can be converted into, and the DFIG (200) rotor current in the stationary coordinate system (
Figure 112020004606316-pat00226
,
Figure 112020004606316-pat00227
) Is also normalized and input to the correction amount calculation unit.

상기 정규화 된 회전자 전류를 기초로 하여 도 6의 피드백인 회전자 위치 보정값

Figure 112020004606316-pat00228
을 제외하고, 상기 회전자 위치 보정부(1200)의 입력을 고려하면 하기 수학식 20 내지 수학식 22와 같다.Based on the normalized rotor current, the rotor position correction value, which is the feedback of FIG. 6,
Figure 112020004606316-pat00228
Except for, when the input of the rotor position correction unit 1200 is considered, Equations 20 to 22 are shown below.

[수학식 20][Equation 20]

Figure 112020004606316-pat00229
Figure 112020004606316-pat00229

[수학식 21][Equation 21]

Figure 112020004606316-pat00230
Figure 112020004606316-pat00230

[수학식 22][Equation 22]

Figure 112020004606316-pat00231
Figure 112020004606316-pat00231

상기 수학식 20 내지 수학식 22 및 이득 상수

Figure 112020004606316-pat00232
에 기초하여, 신호
Figure 112020004606316-pat00233
는 하기 수학식 23과 같이 도출될 수 있다.Equations 20 to 22 and gain constant
Figure 112020004606316-pat00232
Based on the signal
Figure 112020004606316-pat00233
Can be derived as in Equation 23 below.

[수학식 23][Equation 23]

Figure 112020004606316-pat00234
Figure 112020004606316-pat00234

상기 신호

Figure 112020004606316-pat00235
는 적분기(1231)로 입력되어 신호
Figure 112020004606316-pat00236
가 출력되고, 비례적분제어기(1232)에 의해서 회전자 위치 보정값
Figure 112020004606316-pat00237
이 출력된다. 여러 사이클에 걸쳐 도 6의 폐루프가 수행되고, 상기 입력 각주파수
Figure 112020004606316-pat00238
이 상기 전원(100) 주파수
Figure 112020004606316-pat00239
에 근접하면 상기 수학식 23의
Figure 112020004606316-pat00240
이 DC 성분으로 간주되고, 이는 상기 적분기(1231)에 의해 신호
Figure 112020004606316-pat00241
에 영향을 줄 수 있다.Reminder signal
Figure 112020004606316-pat00235
Is input to the integrator 1231 and the signal
Figure 112020004606316-pat00236
Is output, and the rotor position correction value by the proportional integral controller 1232
Figure 112020004606316-pat00237
Is output. The closed loop of FIG. 6 is performed over several cycles, and the input angular frequency
Figure 112020004606316-pat00238
This power supply 100 frequency
Figure 112020004606316-pat00239
When it is close to, in Equation 23
Figure 112020004606316-pat00240
Is regarded as a DC component, which is signaled by the integrator 1231
Figure 112020004606316-pat00241
Can affect

도 7은 상기 회전자 위치 보정값

Figure 112020004606316-pat00242
과 상기 입력 각주파수
Figure 112020004606316-pat00243
의 조건에 따라 상기 신호
Figure 112020004606316-pat00244
의 변화 그래프를 도시하고 있다. 상기 회전자 위치 보정값
Figure 112020004606316-pat00245
과 상기 입력 각주파수
Figure 112020004606316-pat00246
의 조건은 표 1에 도시되어 있다.7 is the rotor position correction value
Figure 112020004606316-pat00242
And the input angular frequency
Figure 112020004606316-pat00243
The signal according to the condition of
Figure 112020004606316-pat00244
It shows the graph of change. The rotor position correction value
Figure 112020004606316-pat00245
And the input angular frequency
Figure 112020004606316-pat00246
The conditions of are shown in Table 1.

CaseCase

Figure 112020004606316-pat00247
Figure 112020004606316-pat00247
Figure 112020004606316-pat00248
Figure 112020004606316-pat00248
Case 1Case 1
Figure 112020004606316-pat00249
Figure 112020004606316-pat00249
Figure 112020004606316-pat00250
Figure 112020004606316-pat00250
Case 2Case 2
Figure 112020004606316-pat00251
Figure 112020004606316-pat00251
Figure 112020004606316-pat00252
Figure 112020004606316-pat00252
Case 3Case 3
Figure 112020004606316-pat00253
Figure 112020004606316-pat00253
Figure 112020004606316-pat00254
Figure 112020004606316-pat00254
Case 4Case 4
Figure 112020004606316-pat00255
Figure 112020004606316-pat00255
Figure 112020004606316-pat00256
Figure 112020004606316-pat00256
Case 5Case 5
Figure 112020004606316-pat00257
Figure 112020004606316-pat00257
Figure 112020004606316-pat00258
Figure 112020004606316-pat00258
Case 6Case 6
Figure 112020004606316-pat00259
Figure 112020004606316-pat00259
Figure 112020004606316-pat00260
Figure 112020004606316-pat00260

도 7과 표 1에 도시된 바와 같이 신호

Figure 112020004606316-pat00261
는 상기 제 1차 적분 필터링 및 정규화 된 d축 회전자 전류
Figure 112020004606316-pat00262
가 상기 정지좌표계상의 정규화 된 d축 회전자 전류
Figure 112020004606316-pat00263
에 뒤처지는 위상일 때, 양의 값을 나타내고, 상기 제1차 적분 필터링 및 정규화 된 d축 회전자 전류
Figure 112020004606316-pat00264
가 상기 정지좌표계상의 정규화 된 d축 회전자 전류
Figure 112020004606316-pat00265
에 앞서는 위상일 때, 음의 값을 나타낸다. 그렇기 때문에, 상기 제 1차 적분 필터링 된 d축 및 q축 회전자 전류
Figure 112020004606316-pat00266
와 상기 정지좌표계상의 d축 및 q축의 회전자 전류
Figure 112020004606316-pat00267
사이의 위상 차이를 얻기 위해서, 신호
Figure 112020004606316-pat00268
Figure 112020004606316-pat00269
의 피드백을 이용하여
Figure 112020004606316-pat00270
의 위상을 이동시킴으로써, 0으로 제어 되어야 한다.Signal as shown in Figure 7 and Table 1
Figure 112020004606316-pat00261
Is the first-order integral filtered and normalized d-axis rotor current
Figure 112020004606316-pat00262
Is the normalized d-axis rotor current in the stationary coordinate system
Figure 112020004606316-pat00263
When the phase lags behind, it represents a positive value, and the first order integral filtering and normalized d-axis rotor current
Figure 112020004606316-pat00264
Is the normalized d-axis rotor current in the stationary coordinate system
Figure 112020004606316-pat00265
When it is a phase preceding a, it represents a negative value. Therefore, the first-order integral filtered d-axis and q-axis rotor currents
Figure 112020004606316-pat00266
And rotor currents in the d-axis and q-axis on the stationary coordinate system
Figure 112020004606316-pat00267
In order to get the phase difference between the signal
Figure 112020004606316-pat00268
Is
Figure 112020004606316-pat00269
Using feedback from
Figure 112020004606316-pat00270
By shifting the phase of, it should be controlled to zero.

슬립각 추정Slip angle estimation

도 8은 본 발명의 바람직한 실시예에 따른 DFIG 센서리스 제어 시스템의 슬립각 추정에 대한 블록도를 도시하고 있다.8 is a block diagram of estimating a slip angle of a DFIG sensorless control system according to a preferred embodiment of the present invention.

도 8에서 도시된 바와 같이, 추정 회전자 위치

Figure 112020004606316-pat00271
은 상기 각속도 추정부(1100)로부터 추정된 회전자 위치
Figure 112020004606316-pat00272
와 상기 회전자 위치 보정부(1200)에서 출력된 회전자 위치 보정값
Figure 112020004606316-pat00273
이 더해짐으로써 산출된다.As shown in Fig. 8, the estimated rotor position
Figure 112020004606316-pat00271
Is the rotor position estimated from the angular velocity estimating unit 1100
Figure 112020004606316-pat00272
And the rotor position correction value output from the rotor position correction unit 1200
Figure 112020004606316-pat00273
Is calculated by adding

여기에서 추정 슬립각

Figure 112020004606316-pat00274
은 위상동기루프(PLL) 알고리즘에 의해 얻어진 전원(100) 위상각
Figure 112020004606316-pat00275
로부터 상기 추정 회전자 위치
Figure 112020004606316-pat00276
을 감산함으로써 산출될 수 있다.Where the estimated slip angle
Figure 112020004606316-pat00274
Is the phase angle of the power supply (100) obtained by the phase-locked loop (PLL) algorithm
Figure 112020004606316-pat00275
From the estimated rotor position
Figure 112020004606316-pat00276
It can be calculated by subtracting.

동기화 과정에서의 센서리스 제어 방식Sensorless control method in synchronization process

도 10은 본 발명의 일실시예에 따른 DFIG 센서리스 제어 시스템의 동기화 과정에서의 슬립 추정 방법에 대한 블록도를 도시하고 있다.10 is a block diagram of a slip estimation method in a synchronization process of a DFIG sensorless control system according to an embodiment of the present invention.

도 10에 도시된 바와 같이, 동기화 과정에서는 특정 슬립 추정 방법이 요구되는데, 전원(100)측 전압과 DFIG(200) 고정자 전압 사이의 위상차에 대한 사인값

Figure 112020004606316-pat00277
를 입력으로 하여 동기화 과정에서의 슬립
Figure 112020004606316-pat00278
이 산출될 수 있다. 상기 전원(100)측 전압과 고정자 전압 사이의 위상차에 대한 사인값
Figure 112020004606316-pat00279
는 하기 수학식 24와 같이 표현될 수 있다.As shown in FIG. 10, a specific slip estimation method is required in the synchronization process, and the sine value of the phase difference between the power supply 100 side voltage and the DFIG 200 stator voltage
Figure 112020004606316-pat00277
Slip in the synchronization process by entering
Figure 112020004606316-pat00278
Can be calculated. Sine value of the phase difference between the power supply 100 side voltage and the stator voltage
Figure 112020004606316-pat00279
May be expressed as in Equation 24 below.

[수학식 24][Equation 24]

Figure 112020004606316-pat00280
Figure 112020004606316-pat00280

(이 때,

Figure 112020004606316-pat00281
: 전원(100)측 위상각,
Figure 112020004606316-pat00282
: 고정자 전압 위상각,
Figure 112020004606316-pat00283
,
Figure 112020004606316-pat00284
: 고정자 전압과 전원(100)측 전압 사이의 위상차)(At this time,
Figure 112020004606316-pat00281
: Power supply 100 side phase angle,
Figure 112020004606316-pat00282
: Stator voltage phase angle,
Figure 112020004606316-pat00283
,
Figure 112020004606316-pat00284
: The phase difference between the stator voltage and the voltage on the power supply (100) side)

상기

Figure 112020004606316-pat00285
값이 비례적분제어기(1410)에 입력되어 추정 슬립 각속도
Figure 112020004606316-pat00286
이 출력되고, 상기 추정 슬립 각속도
Figure 112020004606316-pat00287
는 적분기(1420)에 입력되어 상기 동기화 과정 슬립
Figure 112020004606316-pat00288
이 산출된다.remind
Figure 112020004606316-pat00285
The value is input to the proportional integral controller 1410 to determine the estimated slip angular velocity.
Figure 112020004606316-pat00286
Is output, and the estimated slip angular velocity
Figure 112020004606316-pat00287
Is input to the integrator 1420 to slip the synchronization process
Figure 112020004606316-pat00288
Is calculated.

상기 고정자 전압 위상각

Figure 112020004606316-pat00289
는 상기 동기화 과정에서의 슬립
Figure 112020004606316-pat00290
을 제어함으로써 상기 전원(100)측 위상각
Figure 112020004606316-pat00291
와 매칭시킬 수 있다. 이와 동시에 회전자 속도
Figure 112020004606316-pat00292
은 하기 수학식 25와 같이 산출될 수 있다.The stator voltage phase angle
Figure 112020004606316-pat00289
Is a sleep in the synchronization process
Figure 112020004606316-pat00290
By controlling the phase angle of the power supply 100 side
Figure 112020004606316-pat00291
Can be matched with. At the same time, the rotor speed
Figure 112020004606316-pat00292
May be calculated as in Equation 25 below.

[수학식 25][Equation 25]

Figure 112020004606316-pat00293
Figure 112020004606316-pat00293

(이 때,

Figure 112020004606316-pat00294
: 전원측 주파수)(At this time,
Figure 112020004606316-pat00294
: Power side frequency)

동기화 과정은 고정자 전압이 진폭, 위상각 및 주파수의 관점에서 전원(100) 전압과 동기화되었을 때, 고정자 단자가 전원(100)측 단자와 전기적으로 연결됨으로써 완료될 수 있다.The synchronization process may be completed when the stator voltage is synchronized with the power supply 100 voltage in terms of amplitude, phase angle, and frequency, by electrically connecting the stator terminal to the power supply 100 side terminal.

제 1 실시예Embodiment 1

도 10 및 도 11는 본 발명의 바람직한 실시예에 따른 DFIG 센서리스 제어 시스템의 실험 결과 그래프를 도시하고 있다.10 and 11 are graphs of experimental results of a DFIG sensorless control system according to a preferred embodiment of the present invention.

하기 표 2는 DFIG 센서리스 제어 시스템의 실험에 사용된 DFIG(200)의 파라미터 값을 도시하고 있다.Table 2 below shows the parameter values of the DFIG 200 used in the experiment of the DFIG sensorless control system.

파라미터(Parameters)Parameters 값(Values)Values 고정자 주파수Stator frequency 60Hz60Hz 고정자 정격 전력Stator rated power 3kW3kW 고정자 정격 선간 전압Stator rated line voltage 220V220V 고정자 정격 전류Stator rated current 14.7A14.7A 극 쌍의 수Number of pole pairs 22 고정자 저항Stator resistance

Figure 112020004606316-pat00295
Figure 112020004606316-pat00295
회전자 저항Rotor resistance
Figure 112020004606316-pat00296
Figure 112020004606316-pat00296
자화 인덕턴스Magnetizing inductance
Figure 112020004606316-pat00297
H
Figure 112020004606316-pat00297
H
고정자 누설 인덕턴스Stator leakage inductance
Figure 112020004606316-pat00298
H
Figure 112020004606316-pat00298
H

도 10의 (a) 내지 도 10의 (f)는 초기 회전자 속도가 1600rpm인 동기화 모드에서 DFIG 센서리스 제어 시스템의 성능을 도시하고 있다. 도 10의 (a)는 고정자 전압의 진폭을 조정하기 위해 d축 회전자 전류가 제어되고 있음을 보여주고, 도 10의 (b)는 q축 회전자 전류가 0으로 제어되는 것에 대한 그래프가 도시되어 있다. d축 및 q축의 기준 회전자 전류 (

Figure 112020004606316-pat00299
,
Figure 112020004606316-pat00300
)에 비해 출력되는 회전자 전류 (
Figure 112020004606316-pat00301
,
Figure 112020004606316-pat00302
)는 약간의 오차와 리플이 있지만, 상기 기준 회전자 전류 값에 맞게 제어가 되는 것을 보여주고 있다.10A to 10F illustrate the performance of the DFIG sensorless control system in the synchronization mode in which the initial rotor speed is 1600rpm. FIG. 10(a) shows that the d-axis rotor current is controlled to adjust the amplitude of the stator voltage, and FIG. 10(b) shows a graph of the q-axis rotor current being controlled to zero. Has been. Reference rotor current in d-axis and q-axis (
Figure 112020004606316-pat00299
,
Figure 112020004606316-pat00300
) Compared to the output rotor current (
Figure 112020004606316-pat00301
,
Figure 112020004606316-pat00302
) Shows that there is some error and ripple, but it is controlled according to the reference rotor current value.

또한 도 10의 (c)와 (d)는 고정자와 전원에서의 전압과 위상각에 대한 그래프가 도시되어 있다. 위상이 다른 상기 고정자 전압과 상기 전원(100)측 전압이 진폭, 위상각 및 주파수 관점에서 동기화가 되는 그래프에 대해서 도시되어 있다. In addition, (c) and (d) of FIG. 10 are graphs of voltages and phase angles of the stator and the power source. A graph in which the stator voltage having different phases and the voltage on the power supply 100 side are synchronized in terms of amplitude, phase angle, and frequency is shown.

또한 도 10의 (e)는 기준 회전자 속도

Figure 112020004606316-pat00303
에 대한 추정 회전자 속도
Figure 112020004606316-pat00304
의 오차
Figure 112020004606316-pat00305
의 그래프가 도시되어 있으며, 도 10의 (f)는 추정된 슬립각
Figure 112020004606316-pat00306
그래프를 도시하고 있다. In addition, (e) of Figure 10 is the reference rotor speed
Figure 112020004606316-pat00303
Estimated rotor speed for
Figure 112020004606316-pat00304
Error of
Figure 112020004606316-pat00305
A graph of is shown, and (f) of FIG. 10 is an estimated slip angle
Figure 112020004606316-pat00306
It shows a graph.

도 10의 (a) 내지 도 10의 (f)에서 도시하고 있는 바와 같이, 동기화 프로세스는 약 2 주기 정도의 시간이 필요하며, 도달하고자 하는 기준 회전자 전류 및 기준 회전자 각속도에 대해서 약간의 리플과 오차가 발생하지만, 동기화에 걸리는 시간인 2 주기 정도 후에는 DFIG 고정자 전압과 상기 전원(100)측 전압이 동기화가 된다는 것을 그래프에서 도시하고 있다.As shown in FIGS. 10A to 10F, the synchronization process requires a time of about 2 cycles, and a slight ripple for the reference rotor current and the reference rotor angular velocity to be reached. The over error occurs, but the graph shows that the DFIG stator voltage and the power supply 100 side voltage become synchronized after about 2 cycles, which is the time required for synchronization.

도 11은 풍속 변화 조건에서 구동하는 DFIG의 센서리스 제어 성능을 도시하고 있다.11 shows the sensorless control performance of the DFIG driven under a wind speed change condition.

도 11의 (a)에 도시된 바와 같이 평균 풍속 8.4m/s에서 변동이 있는 풍속으로 실험 조건이 주어졌다. 상기 DFIG의 불리한 구동 조건에 대한 센서리스 제어 성능 그래프가 도 11의 (b) 내지 도 11의 (g)에 도시되어 있다.As shown in (a) of FIG. 11, the experimental conditions were given as the wind speed fluctuating at an average wind speed of 8.4 m/s. Sensorless control performance graphs for the unfavorable driving conditions of the DFIG are shown in FIGS. 11B to 11G.

도 11의 (b) 및 도 11의 (c)에 도시되어 있는 바와 같이, 동기화 이전의 제어에서는 d축 회전자 전류가 제어되었으나, 동기화 이후, MPPT(Maximum Power Point Tracking) 제어로써 풍력 발전 시스템이 제어될 때, q축 회전자 전류가 제어된다. 이에 따라 제어되는 유효전력, 무효전력, 추정 각속도, 각속도 오차는 도 10의 (d) 내지 도 10의 (g)에 도시되어 있다.11(b) and 11(c), in the control before synchronization, the d-axis rotor current was controlled, but after synchronization, the wind power generation system was controlled by MPPT (Maximum Power Point Tracking) control. When controlled, the q-axis rotor current is controlled. Active power, reactive power, estimated angular velocity, and angular velocity errors controlled accordingly are shown in FIGS. 10D to 10G.

도 10에 도시되어 있는 바와 같이, 고조파의 영향으로 인해 리플이 있지만, 아주 미비한 수준이고, 각각 기준 제어 값의 평균치에 맞게 제어 되고 있다. 도 10의 그래프를 통해 풍속의 변동이 있는 상황에서도 강인하고, 정확하게 제어될 수 있는 시스템이라는 장점이 있다.As shown in Fig. 10, there is ripple due to the influence of harmonics, but the level is very insignificant, and each is controlled to fit the average value of the reference control values. Through the graph of FIG. 10, there is an advantage of a system that is robust and can be accurately controlled even in a situation where wind speed fluctuates.

도 11의 (a) 내지 도 11의 (f)는 본 발명의 바람직한 실시예에 따른 DFIG 센서리스 제어 시스템의 이상 조건 적용시의 실험 그래프를 도시하고 있다.11(a) to 11(f) illustrate experimental graphs when applying abnormal conditions of the DFIG sensorless control system according to a preferred embodiment of the present invention.

도 11의 (a)에 도시된 바와 같이 전원 측의 전압에 대해서까지 교란된 전원(100)측 전압 조건이 적용되고, 두 점선 사이의 구간인 왜곡 및 불균형 조건이 함께 적용된 구간에서는 유효전력 및 정 각속도 오차가 변동이 생기는 형태의 그래프가 출력되었지만, 유효전력과 추정 각속도 모두 평균값에는 이상이 없을 정도로 제어 되었고, 각속도 오차는 최대 33rpm 으로 출력된 기준 각속도 및 추정 각속도의 1.71% 정도 수준으로 아주 미비한 수준인 것으로 그래프에서 도시하고 있다.As shown in (a) of FIG. 11, the voltage condition of the power supply 100, which is disturbed up to the voltage of the power supply side, is applied, and in the section in which the distortion and unbalance conditions, which are the section between the two dotted lines, are applied together, the active power and positive A graph in which the angular velocity error fluctuates was output, but both the active power and the estimated angular velocity were controlled so that there was no abnormality in the average value. Is shown in the graph.

본 발명의 바람직한 실시예에 따른 DFIG 센서리스 제어 시스템의 실험에서 시사하는 바와 같이 속도 검출 센서가 없이 강인하고, 정확하게 시스템이 제어될 수 있는 DFIG 센서리스 제어 시스템을 제공할 수 있다는 장점이 있다.As suggested in the experiment of the DFIG sensorless control system according to the preferred embodiment of the present invention, there is an advantage in that it is possible to provide a DFIG sensorless control system that is robust and capable of accurately controlling the system without a speed detection sensor.

이상과 같이 본 발명에서는 구체적인 구성 등과 같은 특정 사항들과 한정된 실시예 도면에 의해 설명되었으나 이는 본 발명의 보다 전반적인 이해를 돕기 위해서 제공된 것 일 뿐, 본 발명은 상기의 일 실시예에 한정되는 것이 아니며, 본 발명이 속하는 분야에서 통상의 지식을 가진 자라면 이러한 기재로부터 다양한 수정 및 변형이 가능하다.As described above, in the present invention, specific matters such as specific configurations and the like have been described with reference to the drawings of limited embodiments, but this is provided only to help a more general understanding of the present invention, and the present invention is not limited to the above-described embodiment. , If one of ordinary skill in the field to which the present invention belongs, various modifications and variations are possible from these descriptions.

따라서, 본 발명의 사상은 설명된 실시예에 국한되어 정해져서는 아니 되며, 후술하는 특허 청구 범위뿐 아니라 이 특허 청구 범위와 균등하거나 등가적 변형이 있는 모든 것들은 본 발명 사상의 범주에 속한다고 할 것이다.Therefore, the spirit of the present invention is limited to the described embodiments and should not be determined, and all things that are equivalent or equivalent to the scope of the claims as well as the claims to be described later belong to the scope of the spirit of the present invention. .

100 : 전원(Grid)
200 : 이중 여자 유도 발전기(DFIG)
300 : 제 1 컨버터(RSC)
400 : 제 2 컨버터(GSC)
500 : 프로펠러(Propeller)
600 : 기어박스(Gearbox)
1000 : 제어부(Controller)
1100 : 각속도 추정부
1110 : 제 1 좌표변환부
1120 : 필터부 1121 : 제 1 적분 필터
1122 : 제 2 적분 필터
1130 : 입력 각주파수 추정부
1140 : 각속도 산출부 1141 : 각속도 산출부 비례적분 제어기
1150 : 제 1 회전자 위치 추정부
1200 : 회전자 위치 보정부
1210 : 정규화부
1220 : 제 2 좌표변환부
1230 : 보정량 산출부 1231 : 보정량 산출부 적분기
1232 : 보정량 산출부 비례적분제어기
1300 : 슬립각 추정부
1310 : 제 2 회전자 위치 추정부
1320 : 슬립각 추정기
1400 : 동기화 과정 슬립각 추정기
1410 : 동기화 과정 슬립각 추정기 비례적분 제어기
1420 : 동기화 과정 슬립각 추정기 적분기
100: Power (Grid)
200: Double Excitation Induction Generator (DFIG)
300: first converter (RSC)
400: second converter (GSC)
500: Propeller
600: Gearbox
1000: Controller
1100: angular velocity estimation unit
1110: first coordinate conversion unit
1120: filter unit 1121: first integral filter
1122: second integral filter
1130: input angular frequency estimation unit
1140: angular velocity calculation unit 1141: angular velocity calculation unit proportional integral controller
1150: first rotor position estimation unit
1200: rotor position correction unit
1210: normalization unit
1220: second coordinate conversion unit
1230: correction amount calculation unit 1231: correction amount calculation unit integrator
1232: Correction amount calculation part proportional integral controller
1300: slip angle estimation unit
1310: second rotor position estimation unit
1320: slip angle estimator
1400: Synchronization process slip angle estimator
1410: Synchronization process slip angle estimator proportional integral controller
1420: Synchronization process slip angle estimator integrator

Claims (12)

DFIG(Doubly Fed Induction Generator) 센서리스 제어 시스템에 있어서,
전원에 고정자 권선이 연결된 DFIG;
상기 DFIG 의 회전자 권선과 연결된 제 1 컨버터; 및
상기 제 1 컨버터를 제어하는 제어부;를 포함하고,
상기 제어부는,
상기 DFIG 회전자에 흐르는 전류를 기초로 상기 DFIG 회전자의 각속도를 추정하는 각속도 추정부를 포함하되,
상기 각속도 추정부는 왜곡된 전류에 의해서 추정 오차가 커지지 않도록 복수 차수의 적분 필터를 포함하는 것;
을 특징으로 하는 DFIG 센서리스 제어 시스템.
In the Doubly Fed Induction Generator (DFIG) sensorless control system,
DFIG with the stator winding connected to the power source;
A first converter connected to the rotor winding of the DFIG; And
Includes; a control unit for controlling the first converter,
The control unit,
Including an angular velocity estimation unit for estimating the angular velocity of the DFIG rotor based on the current flowing through the DFIG rotor,
The angular velocity estimating unit includes a multi-order integration filter so that an estimation error is not increased by a distorted current;
DFIG sensorless control system, characterized in that.
제 1항에 있어서,
상기 각속도 추정부는,
상기 DFIG의 회전자 전류 (
Figure 112020004606316-pat00307
,
Figure 112020004606316-pat00308
)을 입력받아서 정지좌표계로 좌표 변환하여 정지좌표계상의 DFIG 회전자 전류 (
Figure 112020004606316-pat00309
,
Figure 112020004606316-pat00310
)를 생성하는 제 1 좌표변환부;
상기 제 1 좌표변환부로부터 상기 정지좌표계상의 DFIG 회전자 전류 (
Figure 112020004606316-pat00311
,
Figure 112020004606316-pat00312
)를 입력받고, 상기 복수 차수의 적분 필터를 이용하여 제 1차 적분 필터링 된 회전자 전류 (
Figure 112020004606316-pat00313
,
Figure 112020004606316-pat00314
)와, 상기 제1차 적분 필터링 된 회전자 전류 (
Figure 112020004606316-pat00315
,
Figure 112020004606316-pat00316
)에 비해서 90도 위상차를 갖는 제 2차 적분 필터링 된 회전자 전류 (
Figure 112020004606316-pat00317
,
Figure 112020004606316-pat00318
)를 산출하는 필터부;
상기 필터부로부터 상기 제 1, 2차 적분 필터링 된 회전자 전류 (
Figure 112020004606316-pat00319
,
Figure 112020004606316-pat00320
,
Figure 112020004606316-pat00321
,
Figure 112020004606316-pat00322
) 및 상기 전원의 각주파수
Figure 112020004606316-pat00323
를 입력받고, 입력 각주파수
Figure 112020004606316-pat00324
을 산출하는 입력 각주파수 추정부;
상기 입력 각주파수
Figure 112020004606316-pat00325
과 상기 전원의 각 주파수
Figure 112020004606316-pat00326
를 기초로 상기 DFIG 회전자의 추정 각속도
Figure 112020004606316-pat00327
을 산출하는 각속도 산출부;를 포함하는 것
을 특징으로 하는 DFIG 센서리스 제어 시스템.
The method of claim 1,
The angular velocity estimation unit,
The rotor current of the DFIG (
Figure 112020004606316-pat00307
,
Figure 112020004606316-pat00308
) Is input and the coordinates are converted to the stationary coordinate system, and the DFIG rotor current in the stationary coordinate system (
Figure 112020004606316-pat00309
,
Figure 112020004606316-pat00310
A first coordinate conversion unit that generates );
DFIG rotor current in the stationary coordinate system from the first coordinate conversion unit (
Figure 112020004606316-pat00311
,
Figure 112020004606316-pat00312
) Is input, and the first-order integral filtered rotor current (
Figure 112020004606316-pat00313
,
Figure 112020004606316-pat00314
) And the first-order integral filtered rotor current (
Figure 112020004606316-pat00315
,
Figure 112020004606316-pat00316
The second-order integral filtered rotor current (
Figure 112020004606316-pat00317
,
Figure 112020004606316-pat00318
A filter unit that calculates );
The first and second integrally filtered rotor currents from the filter unit (
Figure 112020004606316-pat00319
,
Figure 112020004606316-pat00320
,
Figure 112020004606316-pat00321
,
Figure 112020004606316-pat00322
) And each frequency of the power supply
Figure 112020004606316-pat00323
Is input, input angular frequency
Figure 112020004606316-pat00324
An input angular frequency estimating unit that calculates;
The input angular frequency
Figure 112020004606316-pat00325
And each frequency of the power supply
Figure 112020004606316-pat00326
Based on the estimated angular velocity of the DFIG rotor
Figure 112020004606316-pat00327
Including; an angular velocity calculation unit that calculates
DFIG sensorless control system, characterized in that.
제 2항에 있어서,
상기 필터부는,
상기 정지좌표계상의 DFIG 회전자 전류 (
Figure 112020004606316-pat00328
,
Figure 112020004606316-pat00329
)와 상기 제 1차 적분 필터링 된 회전자 전류 (
Figure 112020004606316-pat00330
,
Figure 112020004606316-pat00331
)의 차이에 기초한 값에서 상기 제 2차 적분 필터링 된 회전자 전류 (
Figure 112020004606316-pat00332
,
Figure 112020004606316-pat00333
)를 뺀 값에 상기 전원의 각 주파수
Figure 112020004606316-pat00334
를 곱한 후, 적분하여 상기 제 1차 적분 필터링 된 추정 회전자 전류 (
Figure 112020004606316-pat00335
,
Figure 112020004606316-pat00336
)를 산출하는 것
을 특징으로 하는 DFIG 센서리스 제어 시스템.
The method of claim 2,
The filter unit,
DFIG rotor current in the stationary coordinate system (
Figure 112020004606316-pat00328
,
Figure 112020004606316-pat00329
) And the first-order integral filtered rotor current (
Figure 112020004606316-pat00330
,
Figure 112020004606316-pat00331
) At the value based on the difference in the second-order integral filtered rotor current (
Figure 112020004606316-pat00332
,
Figure 112020004606316-pat00333
) Minus each frequency of the power supply
Figure 112020004606316-pat00334
After multiplying by and integrating, the first-order integral filtered estimated rotor current (
Figure 112020004606316-pat00335
,
Figure 112020004606316-pat00336
To yield)
DFIG sensorless control system, characterized in that.
제 3항에 있어서,
상기 필터부는,
상기 제 1차 적분 필터링 된 회전자 전류 (
Figure 112020004606316-pat00337
,
Figure 112020004606316-pat00338
)를 적분하고, 상기 전원의 각주파수
Figure 112020004606316-pat00339
를 곱하여 상기 제 2차 적분 필터링 된 회전자 전류 (
Figure 112020004606316-pat00340
,
Figure 112020004606316-pat00341
)를 산출하는 것
을 특징으로 하는 DFIG 센서리스 제어 시스템.
The method of claim 3,
The filter unit,
The first-order integral filtered rotor current (
Figure 112020004606316-pat00337
,
Figure 112020004606316-pat00338
) And the angular frequency of the power supply
Figure 112020004606316-pat00339
Multiply by the second-order integral filtered rotor current (
Figure 112020004606316-pat00340
,
Figure 112020004606316-pat00341
To yield)
DFIG sensorless control system, characterized in that.
제 2항에 있어서,
상기 입력 각주파수 추정부는,
하기 수학식 1에 기초하여 입력 각주파수
Figure 112020004606316-pat00342
을 산출하는 것;
을 특징으로 하는 DFIG 센서리스 제어 시스템.
[수학식 1]
Figure 112020004606316-pat00343

(이 때,
Figure 112020004606316-pat00344
: 전원의 각주파수,
Figure 112020004606316-pat00345
,
Figure 112020004606316-pat00346
: 제1차 적분 필터링 된 d축 및 q축 추정 회전자 전류,
Figure 112020004606316-pat00347
,
Figure 112020004606316-pat00348
: 제2차 적분 필터링 된 d축 및 q축 추정 회전자 전류)
The method of claim 2,
The input angular frequency estimation unit,
Input angular frequency based on Equation 1 below
Figure 112020004606316-pat00342
To yield;
DFIG sensorless control system, characterized in that.
[Equation 1]
Figure 112020004606316-pat00343

(At this time,
Figure 112020004606316-pat00344
: Each frequency of power,
Figure 112020004606316-pat00345
,
Figure 112020004606316-pat00346
: First-order integral filtered d-axis and q-axis estimated rotor current,
Figure 112020004606316-pat00347
,
Figure 112020004606316-pat00348
: Second-order integral filtered d-axis and q-axis estimated rotor current)
제 2항에 있어서,
상기 각속도 산출부는,
상기 입력 각주파수
Figure 112021016407915-pat00349
과 상기 전원의 각주파수
Figure 112021016407915-pat00350
의 차이를 입력받아 상기 DFIG 회전자의 추정 각속도
Figure 112021016407915-pat00351
을 출력하는 비례 적분 제어기를 포함하는 것;
을 특징으로 하는 DFIG 센서리스 제어 시스템.
The method of claim 2,
The angular velocity calculation unit,
The input angular frequency
Figure 112021016407915-pat00349
And each frequency of the power supply
Figure 112021016407915-pat00350
Receive the difference of the estimated angular velocity of the DFIG rotor
Figure 112021016407915-pat00351
Including a proportional integral controller for outputting;
DFIG sensorless control system, characterized in that.
제 2항에 있어서,
상기 제 1 좌표변환부는,
상기 각속도 산출부로부터 입력받는 상기 DFIG 회전자의 추정 각속도
Figure 112020004606316-pat00352
을 적분하여 얻어진 추정 회전자 위치
Figure 112020004606316-pat00353
를 이용하여 정지좌표계로 좌표 변환하는 것;
을 특징으로 하는 DFIG 센서리스 제어 시스템.
The method of claim 2,
The first coordinate conversion unit,
Estimated angular velocity of the DFIG rotor received from the angular velocity calculation unit
Figure 112020004606316-pat00352
The estimated rotor position obtained by integrating
Figure 112020004606316-pat00353
Converting the coordinates to the stationary coordinate system by using;
DFIG sensorless control system, characterized in that.
제 2항에 있어서,
상기 DFIG 회전자의 추정 각속도
Figure 112021016407915-pat00354
을 적분하여 산출되는 추정 회전자 위치
Figure 112021016407915-pat00355
와 회전자 위치 보정부에서 산출되는 회전자 위치 보정값
Figure 112021016407915-pat00356
을 기초로 상기 DFIG의 회전자 추정 위치를 산출하는 제 2 회전자 위치 추정부;를 더 포함하는 것
을 특징으로 하는 DFIG 센서리스 제어 시스템.
The method of claim 2,
Estimated angular velocity of the DFIG rotor
Figure 112021016407915-pat00354
Estimated rotor position calculated by integrating
Figure 112021016407915-pat00355
And the rotor position correction value calculated by the rotor position correction unit
Figure 112021016407915-pat00356
A second rotor position estimating unit for calculating the estimated rotor position of the DFIG on the basis of; further comprising
DFIG sensorless control system, characterized in that.
제 8항에 있어서,
상기 회전자 위치 보정부는,
상기 제 1차 적분 필터링 된 DFIG 회전자 전류 (
Figure 112020004606316-pat00357
,
Figure 112020004606316-pat00358
)를 입력받아서 정규화 하여 제 1차 적분 필터링 및 정규화 된 회전자 전류 (
Figure 112020004606316-pat00359
,
Figure 112020004606316-pat00360
)를 출력하는 정규화부;
상기 제 1차 적분 필터링 및 정규화 된 회전자 전류 (
Figure 112020004606316-pat00361
,
Figure 112020004606316-pat00362
)를 좌표 변환하여 d상의 정규화 된 회전자 전류
Figure 112020004606316-pat00363
를 산출하는 제 2 좌표변환부; 및
상기 좌표 변환되어 얻어진 d상의 정규화 된 회전자 전류
Figure 112020004606316-pat00364
를 기초로 상기 회전자 위치 보정값
Figure 112020004606316-pat00365
을 산출하는 보정량 산출부; 를 포함하는 것
을 특징으로 하는 DFIG 센서리스 제어 시스템.
The method of claim 8,
The rotor position correction unit,
The first-order integral filtered DFIG rotor current (
Figure 112020004606316-pat00357
,
Figure 112020004606316-pat00358
) Received and normalized to filter the first-order integral and normalized rotor current (
Figure 112020004606316-pat00359
,
Figure 112020004606316-pat00360
A normalization unit that outputs );
The first order integral filtered and normalized rotor current (
Figure 112020004606316-pat00361
,
Figure 112020004606316-pat00362
) By converting the coordinates to the normalized rotor current in d phase
Figure 112020004606316-pat00363
A second coordinate conversion unit that calculates; And
The normalized rotor current of the d phase obtained by the above coordinate transformation
Figure 112020004606316-pat00364
Based on the rotor position correction value
Figure 112020004606316-pat00365
A correction amount calculating unit that calculates a; Including
DFIG sensorless control system, characterized in that.
제 9항에 있어서,
상기 보정량 산출부는,
상기 d상의 정규화 된 회전자 전류
Figure 112020004606316-pat00366
와 정지좌표계상의 DFIG의 d축 회전자 전류
Figure 112020004606316-pat00367
이 정규화 된
Figure 112020004606316-pat00368
과의 차이에 정지좌표계상의 DFIG의 q축 회전자 전류
Figure 112020004606316-pat00369
이 정규화 된
Figure 112020004606316-pat00370
을 곱한 값을 적분한 적분값을 기초로 상기 회전자 위치 보정값
Figure 112020004606316-pat00371
을 산출하는 것
을 특징으로 하는 DFIG 센서리스 제어 시스템.
The method of claim 9,
The correction amount calculation unit,
The normalized rotor current of the d phase
Figure 112020004606316-pat00366
And d-axis rotor current of DFIG in stationary coordinate system
Figure 112020004606316-pat00367
Is normalized
Figure 112020004606316-pat00368
DFIG's q-axis rotor current in the stationary coordinate system at the difference with
Figure 112020004606316-pat00369
Is normalized
Figure 112020004606316-pat00370
The rotor position correction value based on the integral value multiplied by
Figure 112020004606316-pat00371
To yield
DFIG sensorless control system, characterized in that.
제 10항에 있어서,
상기 보정량 산출부는,
비례적분제어기에 상기 적분값을 입력하여 상기 회전자 위치 보정값
Figure 112020004606316-pat00372
을 산출하는 것
을 특징으로 하는 DFIG 센서리스 제어 시스템.
The method of claim 10,
The correction amount calculation unit,
The rotor position correction value by inputting the integral value to a proportional integral controller
Figure 112020004606316-pat00372
To yield
DFIG sensorless control system, characterized in that.
제 10항에 있어서,
상기 정지좌표계상의 DFIG 회전자 전류
Figure 112020004606316-pat00373
,
Figure 112020004606316-pat00374

하기 수학식 2 내지 수학식 5로부터 산출되는 것
을 특징으로 하는 DFIG 센서리스 제어 시스템.
[수학식 2]
Figure 112020004606316-pat00375

[수학식 3]
Figure 112020004606316-pat00376

[수학식 4]
Figure 112020004606316-pat00377

[수학식 5]
Figure 112020004606316-pat00378

(이 때,
Figure 112020004606316-pat00379
: 정지좌표계상의 회전자 d축 전류,
Figure 112020004606316-pat00380
: 정지좌표계상의 회전자 q축 전류,
Figure 112020004606316-pat00381
: 정지좌표계상의 회전자 d축 쇄교 자속,
Figure 112020004606316-pat00382
: 정지좌표계상의 회전자 q축 쇄교 자속,
Figure 112020004606316-pat00383
: 자화 인덕턴스,
Figure 112020004606316-pat00384
: 고정자 인덕턴스,
Figure 112020004606316-pat00385
: 정지좌표계상의 고정자 d축 전류,
Figure 112020004606316-pat00386
: 정지좌표계상의 고정자 q축 전류,
Figure 112020004606316-pat00387
: 정지좌표계상의 고정자 d축 전압,
Figure 112020004606316-pat00388
: 정지좌표계상의 고정자 q축 전압,
Figure 112020004606316-pat00389
: 전원의 각주파수,
Figure 112020004606316-pat00390
: 고정자 저항)
The method of claim 10,
DFIG rotor current in the stationary coordinate system
Figure 112020004606316-pat00373
,
Figure 112020004606316-pat00374
Is
What is calculated from the following equations 2 to 5
DFIG sensorless control system, characterized in that.
[Equation 2]
Figure 112020004606316-pat00375

[Equation 3]
Figure 112020004606316-pat00376

[Equation 4]
Figure 112020004606316-pat00377

[Equation 5]
Figure 112020004606316-pat00378

(At this time,
Figure 112020004606316-pat00379
: Rotor d-axis current in the stationary coordinate system,
Figure 112020004606316-pat00380
: Rotor q-axis current in the stationary coordinate system,
Figure 112020004606316-pat00381
: The d-axis linkage flux of the rotor in the stationary coordinate system,
Figure 112020004606316-pat00382
: Rotor q-axis linkage flux in the stationary coordinate system,
Figure 112020004606316-pat00383
: Magnetizing inductance,
Figure 112020004606316-pat00384
: Stator inductance,
Figure 112020004606316-pat00385
: Stator d-axis current on the stationary coordinate system,
Figure 112020004606316-pat00386
: Stator q-axis current on the stationary coordinate system,
Figure 112020004606316-pat00387
: Stator d-axis voltage on the stationary coordinate system,
Figure 112020004606316-pat00388
: Stator q-axis voltage on the stationary coordinate system,
Figure 112020004606316-pat00389
: Each frequency of power,
Figure 112020004606316-pat00390
: Stator resistance)
KR1020200005324A 2020-01-15 2020-01-15 Sensorless control System of DFIG KR102238494B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020200005324A KR102238494B1 (en) 2020-01-15 2020-01-15 Sensorless control System of DFIG

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020200005324A KR102238494B1 (en) 2020-01-15 2020-01-15 Sensorless control System of DFIG

Publications (1)

Publication Number Publication Date
KR102238494B1 true KR102238494B1 (en) 2021-04-09

Family

ID=75444223

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020200005324A KR102238494B1 (en) 2020-01-15 2020-01-15 Sensorless control System of DFIG

Country Status (1)

Country Link
KR (1) KR102238494B1 (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101423212B1 (en) * 2014-05-22 2014-07-24 전북대학교산학협력단 Voltage control system and method at the point of common coupling of wind power plant
KR20140112929A (en) * 2013-03-15 2014-09-24 현대중공업 주식회사 Doubly fed induction generator wind power generating system
JP2014231785A (en) * 2013-05-29 2014-12-11 エムエイチアイ ヴェスタス オフショア ウィンド エー/エス Floating body type wind power generation device
KR101780265B1 (en) 2016-12-23 2017-09-21 (주)파워닉스 Apparatus for bypassing a submodule of a modular multilevel converter

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20140112929A (en) * 2013-03-15 2014-09-24 현대중공업 주식회사 Doubly fed induction generator wind power generating system
JP2014231785A (en) * 2013-05-29 2014-12-11 エムエイチアイ ヴェスタス オフショア ウィンド エー/エス Floating body type wind power generation device
KR101423212B1 (en) * 2014-05-22 2014-07-24 전북대학교산학협력단 Voltage control system and method at the point of common coupling of wind power plant
KR101780265B1 (en) 2016-12-23 2017-09-21 (주)파워닉스 Apparatus for bypassing a submodule of a modular multilevel converter

Similar Documents

Publication Publication Date Title
Abdelrahem et al. Robust predictive control for direct-driven surface-mounted permanent-magnet synchronous generators without mechanical sensors
Liu et al. Control design and experimental verification of the brushless doubly‐fed machine for stand‐alone power generation applications
Amrane et al. Design and implementation of high performance field oriented control for grid-connected doubly fed induction generator via hysteresis rotor current controller
CN106655947B (en) A kind of permanent magnet synchronous motor control algolithm improving low capacity dc-link capacitance voltage transient stability
Monfared et al. Direct active and reactive power control of single-phase grid-tie converters
JP4022630B2 (en) Power conversion control device, power conversion control method, and program for power conversion control
USRE45880E1 (en) Method and system for multiphase current sensing
CA2740404C (en) Sensorless optimum torque control for high efficiency ironless permanent magnet machine
Nian et al. Improved control strategy for stand-alone distributed generation system under unbalanced and non-linear loads
US10491146B2 (en) System and method for compensating for generator-induced flicker in a wind turbine
US9753067B2 (en) Determining a phase and a frequency of an electric quantity of an operating electrical device
US20130027000A1 (en) Reduction of noise and vibrations of an electromechanical transducer by using a modified stator coil drive signal comprising harmonic components
CN110350551B (en) Direct current amplitude-frequency control method and system for voltage source type grid-connected conversion device
JP2008048593A (en) Method and system for controlling permanent magnet motor drive system
Nguyen et al. Sensorless control of DFIG wind turbine systems based on SOGI and rotor position correction
Solís-Chaves et al. A direct power control for DFIG under a three phase symmetrical voltage sag condition
CN107204727B (en) A kind of low capacity thin-film capacitor permanent magnet synchronous motor direct-axis current given controled method
Ramos et al. Rotor current controller with voltage harmonics compensation for a DFIG operating under unbalanced and distorted stator voltage
CN101499665A (en) Idle grid connection control method for speed variant frequency constant dual feedback asynchronous wind power generator set
KR102238494B1 (en) Sensorless control System of DFIG
Nguyen et al. Sensorless control scheme of DFIG wind energy conversion systems based on SOGIs and FLL
JP6592232B2 (en) Control circuit for controlling power conversion circuit, power conversion apparatus including the control circuit, and method
Dinesh et al. Independent operation of DFIG-based WECS using resonant feedback compensators under unbalanced grid voltage conditions
Kumar et al. A Robust Rotor Flux Based S-PLL for Permanent Magnet Synchronous Generator
KR20140039060A (en) Method of regulating the power of an energy conversion installation and energy conversion installation driven by such a method

Legal Events

Date Code Title Description
E701 Decision to grant or registration of patent right
GRNT Written decision to grant