Dynamic evaluation of maximum power point tracking operation with PV array simulator

Dynamic evaluation of maximum power point tracking operation with PV array simulator

Solar Energy Materials & Solar Cells 75 (2003) 537–546 Dynamic evaluation of maximum power point tracking operation with PV array simulator Hiroshi M...

936KB Sizes 1 Downloads 61 Views

Solar Energy Materials & Solar Cells 75 (2003) 537–546

Dynamic evaluation of maximum power point tracking operation with PV array simulator Hiroshi Matsukawaa,*, Koukichi Koshiishib, Hirotaka Koizumib, Kosuke Kurokawab, Masayasu Hamadac, Liu Boc a

Resources Total System Co. Ltd., 5 Floor Kyoritsu Building, 2-3-11 Shinkawa Chuo-ku, Tokyo 140-0033, Japan b Tokyo University of Agriculture and Technology, 2-24-16 Naka-cho, Koganei, Tokyo 184-8588, Japan c Myway Labs Co., Ltd., 1-19-3 Shinyokohama, Kohoku-ku, Yokohama, Kanagawa 222-0033, Japan

Abstract This paper presents an innovative method to measure the dynamic control ability of maximum power point tracking for PV inverters under the condition of irradiance fluctuation. The PV array I–V curve simulator is a kind of indoor testing facility and easy to be adopted by industries. Basic functions are given by a specially designed PV array I–V curve simulator composed of the active power load. Most of the parameters are controllable by sophisticated software with capability of treating a lot of 1-s data for a very long period of time. In this paper, detailed structure of the equipment is described and test examples are also given by using a commercial PV inverter. r 2002 Elsevier Science B.V. All rights reserved. Keywords: Photovoltaic system; Inverter; MPPT; Array simulator; I–V curve; Irradiance

1. Introduction It is effective for the test of PV inverter using PV array I–V curve simulator which is able to reproduce variable conditions. Tokyo University of Agriculture and Technology (TUAT) and Myway Labs Co., Ltd. have developed a PV array I–V curve simulator. The simulator is based on the active power load (APL) of Myway and software for I–V curve simulation of TUAT. The basic functions of APL allow several modes such as constant voltage and/or current power supply, active load without power consumption. Another control loop was provided for allowing PV array I–V curve mode. The definition of I–V curve can be made through specially *Corresponding author. Tel.: +81-3-3551-6345; fax: +81-3-3553-8954. E-mail address: [email protected] (H. Matsukawa). 0927-0248/03/$ - see front matter r 2002 Elsevier Science B.V. All rights reserved. PII: S 0 9 2 7 - 0 2 4 8 ( 0 2 ) 0 0 1 8 5 - X

538

H. Matsukawa et al. / Solar Energy Materials & Solar Cells 75 (2003) 537–546

developed program on a PC. Matsukawa et al. [1,2] have reported that the I–V curve can be simulated by using this program. So, irregular pattern caused by shading etc., can also be simulated. PV array structure is arranged freely. Furthermore, when an I–V curve data set is formed, irradiance fluctuation is easily treated. In this paper, detailed structure of the equipment is described. Then, results of tests are also introduced. Effectiveness of the dynamic evaluation of maximum power point tracking (MPPT) operation with this simulator has been cleared.

2. PV array I–V curve simulator A basic structure of PV array I–V curve simulator including APL [3] is shown in Fig. 1. The APL connected PC by optical fiber to RS-232C. At first, an I–V curve is generated on the PC. Then the data of the I–V curve is transmitted to APL. Following the data, the APL performs as a PV array simulator. While its operation, parameters of APL are transmitted to PC and monitored every 1 s. The I–V, P–V characteristics, value of Pmax and output power are also indicated. Characteristics of the PV array I–V curve simulator are summarized as follows: * *

Simulating irradiance every 1 s for 9-h operation. Preinstalled standard patterns of irradiance fluctuation (additional files are installed).

Fig. 1. Basic structure of PV array I–V curve simulator.

H. Matsukawa et al. / Solar Energy Materials & Solar Cells 75 (2003) 537–546

539

Table 1 Specification of APL Heading

Specification

Utility side

Grid voltage (V) Grid frequency (Hz) Power factor of power grid Current harmonic distortion factor Rated capacity (kVA)

3f/200710% 50/60 Over 0.98 THD p5%, each p3% 710

DC side

DC input/output voltage (V) DC input/output current (A) Input/output capacity

DC 1–400 DC 0B725 77.5 kW continuous (max. 710 kW/10 min)

Other specification

Maximum efficiency Permissible temperature range (1C) Permissible humidity range Size (mm) Weight (kg)

89% (powering/regenerative) 0–40 10–80% RH W430, D550, H305 52

* * * * *

Faithfully simulated output characteristics of PV array (parameters are set). Size, weight, and cost reduction. Operation by a PC with easy handling. Simulating an enough scale of PV array for a standard housing. Further utilization as constant voltage and/or current power supply. Specification of this PV array I–V curve simulator was shown in Table 1.

3. Tests of APL response time The PV array I–V curve simulator (APL) was tested using an experimental circuit model shown in Fig. 2. 3.1. Tests for sudden irradiance fluctuation Fig. 3 shows the output voltage and current waveforms observed by PV array I–V curve simulator. In this test, irradiance was stepped down from 1.0 to 0.5 kW/m2. The load resistance was set to have Pmax at 0.5 kW/m2 of irradiance. The measured response time in this test was approximately 25 ms. 3.2. Tests for rapid load variation Fig. 4 shows output voltage and current waveforms observed by the simulator in the test for rapid load variation from 1 k O to 100 O (a), and from 100 O to 1 k O (b). The measured response time was 16 and 25 ms, respectively.

540

H. Matsukawa et al. / Solar Energy Materials & Solar Cells 75 (2003) 537–546 3 AC200V ±10% 3φ R

AC

S T

APL

200V : 200V

+

VDC IDC

Fig. 2. Circuit model for measurement.

Fig. 3. Result of sudden irradiance fluctuation test.

4. Tests of PV inverter for MPPT control with the simulator 4.1. Default data of irradiance fluctuation The commercial PV inverter was tested with the simulator for MPPT control. Fig. 5 shows the default data sets prepared for a typical clear sky pattern ‘‘Clear’’ and three kinds of fluctuating patterns, ‘‘Vary-1’’, ‘‘Vary-2’’ and ‘‘Vary-3’’. Table 2 shows proportion of the fluctuation for default irradiance patterns. These patterns can examine MPPT function of PV inverter very realistically. These are chosen from actual observation with 1 s sampling at Tsukuba around 4 years since 1996. Relatively, irradiance fluctuation frequency is defined in the following equations. Maximum irradiance fluctuation rate means maximum value of irradiance fluctuation in a second and average irradiance fluctuation means the average value. PN Gi AVR ¼ i¼1 ; ð1Þ N Fi ¼

m X Gi 1 Gi  ; AVR 2m þ 1 i¼m AVR

rffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 1 XN sðFi Þ ¼ ðF  Fi Þ2 ; i¼1 i N

ð2Þ

ð3Þ

H. Matsukawa et al. / Solar Energy Materials & Solar Cells 75 (2003) 537–546

Fig. 4. Result of rapid load variation test.

541

H. Matsukawa et al. / Solar Energy Materials & Solar Cells 75 (2003) 537–546

irradiance

542 1.2

1.2

1

1

0.8

0.8

0.6

0.6

0.4

0.4

0.2

0.2

0 7:00:00

0 7:00:00

9:00:00 11:00:00 13:00:00 15:00:00 17:00:00 (1) Clear

9:00:00 11:00:00 13:00:00 15:00:00 17:00:00 (2) Vary-1 time

time 1.2

irradiance

1.2 1

1

0.8

0.8

0.6

0.6

0.4

0.4

0.2

0.2

0 7:00:00

9:00:00 11:00:00 13:00:00 15:00:00 17:00:00 (3) Vary-2 time

0 7:00:00

9:00:00 11:00:00 13:00:00 15:00:00 17:00:00 (4) Vary-3 time

Fig. 5. Default irradiance patterns.

Table 2 Proportion of fluctuation for default irradiance patterns

Relatively irradiance fluctuation frequency (W/m2)

at 8 a.m.–4 p.m. at 11 a.m.–1 p.m. Maximum irradiance fluctuation rate (W/m2/s) Average irradiance fluctuation rate (W/m2/s) at 8 a.m.–4 p.m. Daily amount of irradiance (Wh/day) at 8 a.m.–4 p.m

Clear

Vary-1

Vary-2

Vary-3

100

363

173

138

100

397

182

169

10 0.55

681 19.94

249 8.31

220 2.59

6500

3930

4840

5150

where AVR is the average of irradiance, Gi the irradiance and Fi the irradiance fluctuating factor. We define a factor KPM for evaluating the inverter’s MPPT control capability P KPM ¼

PVArray Output Power P 100ð%Þ: Pmax

ð4Þ

H. Matsukawa et al. / Solar Energy Materials & Solar Cells 75 (2003) 537–546

543

4.2. Calculation of array temperature variation The resumption of array temperature usually requires many parameters and it is very difficult to simulate the temperature. However, we can easily calculate the temperature by the following equations: TCR ¼ TA þ DT;

ð5Þ

DT ¼ ð6:036 þ 0:274V þ 0:071V 2 Þ þ HAE ð45:63  5:91V þ 0:333V 2 Þ;

ð6Þ

where TCR ; TA ; V ; and HAE are, respectively, array temperature, ambient temperature, wind velocity, and irradiance. TA ¼ 251C and V ¼ 1:0 m/s are assumed. The developed simulator allows an operational mode of high-speed fluctuation of irradiance and corresponding temperature variation.

inverter

3 AC200V ±10% 3φ

1 AC100V 1φ

-

R

AC

S T

200V : 200V

APL

+



V I

4.0kW

I0

Fig. 6. Circuit model for measurement at interactive connected.

Table 3 Specification and rating of the inverter used for the test Rated capacity (kW) Rated input voltage (V) Rated AC output voltage (V) Rated frequency (Hz) Power conversion efficiency Output power factor at fundamental frequency Current harmonic distortion Permissible temperature range (1C) Permissible humidity range (%) Inverter mode

Control method Switching mode Isolation mode Output mode

4.0 DC240 AC202712 50/60 95% (max) Over 95% THDo5%, Eacho3% 10–40 25–85 Utility interactive operating: voltage stiff current control Independency operating: voltage stiff voltage control Maximum power point tracking control leading reactive power control PWM mode Transformerless mode (booster chopper mode) Single-phase two-wire system

H. Matsukawa et al. / Solar Energy Materials & Solar Cells 75 (2003) 537–546

544

4.3. Circuit model of the test for PV inverter’s MPPT control The commercial PV inverter was connected to APL, and the PV inverter was tested for MPPT control. Fig. 6 shows the experimental circuit model and Table 3 shows specification and rating of this PV inverter. 4.4. Experimental results

250

2500

200

2000

150

1500

100

1000

50

500 0 11:00:00

11:30:00 12:00:00

0 12:30:00 13:00:00

250

2500

200

2000

150

1500

100

1000

50

500 0 11:00:00

11:30:00 12:00:00

0 12:30:00 13:00:00

250

2500

200

2000

150

1500

100

1000

50

500 0 11:00:00

11:30:00 12:00:00

0 12:30:00 13:00:00

3500

300

3000

250

2500

200

2000

150

1500

100

1000

50

500 0 11:00:00

11:30:00 12:00:00

0 12:30:00 13:00:00

Voltage [V] Temp. [deg]

300

3000

Power[W] Deviation[W] iation[W]

3500

(c) 300

3000

250

2500

200

2000

150

1500 100

1000

50

500 0 11:00:00

11:30:00 12:00:00

0 12:30:00 13:00:00

(e (e)

3500

300

3000

250

2500

200

2000

150

1500

100

1000

50

500 0 11:00:00

11:30:00 12:00:00

0 12:30:00 13:00:00

Voltage [V] Temp. [deg]

3500

Power[W] Deviation[W] iation[W]

(d)

Voltage [V] Temp. [deg]

Power[W] Deviation[W] iation[W]

300

3000

(b)

Voltage [V] Temp. [deg]

Power[W] Deviation[W] iation[W]

(a (a)

3500

Voltage [V] Temp. [deg]

300

3000

Power[W] Deviation[W] iation[W]

3500

Voltage [V] Temp. [deg]

Power[W] Deviation[W] iation[W]

Using the irradiance data of array temperature, a commercial inverter has been tested. The results are shown in Fig. 7. The power Pmax and voltage VPmax at Pmax

(f (f) Pmax Vpmax

Parray Varray

Deviation Temp.

Fig. 7. Measured and calculated parameters of the tests as functions of time. For the irradiance data: (a) and (b) Vary-1 with 1 and 10 min, (c) and (d) Vary-2 with 1 and 10 min and (e) and (f) Vary-3 with 1 and 10 min.

H. Matsukawa et al. / Solar Energy Materials & Solar Cells 75 (2003) 537–546

545

MPPT Matching Factor [%]

98 Vary-1 Vary-1

97

Vary-2 Vary-2 Vary-3 Vary-3

96 95 94 93 92 0

2

4 6 Temperature Fluctuation Factor[%]

8

10

Fig. 8. An example of array temperature fluctuation effects on MPPT.

point, the output power Parray and voltage Varray of PV array, deviation of Parray from Pmax ; and temperature are shown as functions of time for the irradiance data, Vary-1, Vary-2, and Vary-3. For each case, two sets of data are shown by moving average of irradiance data with 1 and 10 min averaging spans. All the data are stored in the monitoring PC. From the results of these tests, using Eq. (4), KPM for evaluating the inverter’s MPPT control capability of the inverter is calculated and is shown in Fig. 8. Transient control mismatch can be clearly observed by the PV array I–V curve simulator. These results show that MPPT control capability decreases as irradiance of temperature increases.

5. Conclusions The detailed structure of the PV array I–V curve simulator including APL has been described. The results of tests have been shown. The dynamic evaluation of MPPT operation has become possible by means of using this simulator. That is very effective for testing commercial inverters with MPPT control function.

References [1] H. Matsukawa, M. Shioya, T. Yamada, T. Sugiura, K. Kurokawa, Investigation of photovoltaic array simulation method for architecture, Proceedings of the JSES/JWEA Joint Conference’99, November 1999, pp. 57–60 (in Japanese).

546

H. Matsukawa et al. / Solar Energy Materials & Solar Cells 75 (2003) 537–546

[2] H. Matsukawa, M. Shioya, K. Kurokawa, Study on simple assessment method of BIPV power generation for architects, Proceedings of the 28th IEEE PVSC-2000, September 2000, pp. 1648–1651. [3] K. Koshiishi, H. Matsukawa, K. Kurokawa, M. Hamada, Liu bo, Evaluation of maximum power point tracking operation with PV array I–V curve simulator, National convention record IEEJ, March 2001, p. 3042 (in Japanese).