CHƯƠNG 3 : MƠ HÌNH HÓA PV
3.1 Đặc tuyến và thông số cơ bản của PV
Như đã giải thích ở những phần trên, trong cấu trúc của PV có mối nối P-N. Do đó có thể mơ hình hóa PV bằng một diode mắc song song với một nguồn dịng như hình 3.1.
Hình 3. 1.Mơ hình PV
Dựa theo ký hiệu của hình vẽ, áp dụng định luật Kirchhoff cho nút:
I ISC Id (3.1) Áp dụng công thức Shockley cho diode, (3.1) trở thành:
I ISC I0 (e1) (3.2)
Khi khơng có ánh sáng ISC = 0, đặc tuyến I-V của PV giống như đặc tuyến của diode nhưng đảo ngược. Khi có ánh sáng, đặc tuyến đảo ngược đó sẽ được nâng lên một lượng đúng bằng ISC.
Hình 3. 2 Đặc tuyến PV
Trên đặc tuyến của PV có những thơng số cơ bản nhưng rất quan trọng như dòng ngắn mạch, áp hở mạch và điểm công suất cực đại (MPP). Chỉ khi vận hành ở
q V / A kT Hình 3. 3 Thơng số cơ bản PV
Một thông số quan trọng khác của PV là hệ số lấp đầy (Fill Factor-FF), như hình 3.4. Hệ số FF mơ tả chất lượng mối nối bán dẫn của PV và thể hiện khả năng hấp thu ánh sáng để tạo năng lượng.
(3.3)
Hình 3. 4 Hệ số lấp đầy
FF luôn < 1 và thay đổi tùy theo vật liệu PV. FF càng gần bằng 1, PV càng vận hành tốt. Đối với những cell pin chất lượng cao, FF có khi lên đến 0.85 hoặc hơn. Đối với các PV thương mại thông thường, FF chỉ ở mức 0.68.
q V / A kT 3.2 Phương trình tốn của PV
Có nhiều cách, tùy theo độ chính xác khác nhau, để mơ hình hóa PV như mơ hình 1 diode, mơ hình 2 diode. Trong bài này, mơ hình 1 diode và bỏ đi điện trở song song, hình 3.7, được áp dụng.
Hình 3. 5 Mơ hình 1 diode
Hình 3. 6 Mơ hình 2 diode
Hình 3. 7 Mơ hình 1 diode đơn giản
I =
Áp dụng cơng thức Shockley cho diode ta có:
I = (3.4) Trong đó: + q là điện tích electron (1,6 x C)
+ là dòng bão hòa của diode + là dòng ngắn mạch của PV + A là hệ số phẩm chất của diode + k là hằng số Boltzmann = 1,38x
thay đổi theo nhiệt độ: (3.7) ở nhiệt độ chuẩn: (3.8) với Trong đó:
ISC: dòng quang điện, xem như bằng với dòng ngắn mạch tương ứng. I0: dòng bão hòa của diode.
Rs: trở nối tiếp.
T: nhiệt độ tuyệt đối, Kelvin. G: bức xạ, W/m2.
q: điện tích electron = 1.6e-19 Coulomb. k: hằng số Boltzmann = 1.38e-23. A: hệ số phẩm chất của diode. Vg: điện áp bandgap.
Các phương trình trên sẽ được giải theo phương pháp Newton-Ralphson để tính toán (IPV,VPV) ở trường hợp bức xạ và nhiệt độ bất kỳ. Một số thơng số trong phương trình sẽ được chọn để đơn giản tính tốn. PV của hãng Solartech SPM080P với công suất cực đại 80 watt được chọn để mơ hình hóa.
Hình 3. 8 Datasheet solartech
Hệ số K0 được chọn là 0.065% dựa theo hệ số nhiệt của PV. Còn Rs được lấy là 0.005 Ohms.
Hình 4. 1 Điểm hoạt động của tải
Tải trở:
Phương trình mơ tả đơn giản cho tải trở: V = IR hay I =
Đặc tuyến I-V đơn giản là một đường thẳng có độ dốc là Khi R tăng độ dốc của đặc tuyến sẽ hạ thấp, ngược lại khi R giảm độ dốc của đặc tuyến sẽ nâng cao. Đường thẳng I-V của tải trở sẽ giao đường I-V của PV ở điểm làm việc.
Hình 4. 2 Đặc tuyến điện trở
Như hình 4.2, chỉ có trường hợp tải trở là Rm thì cơng suất vận hành của PV là cực đại.
Khi điện trở không thay đổi, điểm làm việc cũng sẽ thay đổi khi điều kiện mơi trường khác đi. Như hình 4.3, chỉ có một điều kiện bức xạ ứng với cơng suất PV là cực đại.
Hình 4. 3 Đặc tuyến trở khi thay đổi bức xạ
4.2 Điểm công suất cực đại 4.2.1 Mục đích của MPPT
Qua những phần trên, cho thấy đối với những tải và điều kiện môi trường khác nhau, điểm vận hành của PV sẽ thay đổi. Khi đó, PV có thể sẽ khơng được vận hành ở công suất cực đại. Để sử dụng tối ưu công suất PV, việc bắt công suất cực đại (mppt) là rất cần thiết.
Hình 4. 4 Buck-Boost
MPPT (Maximum Power Point Tracker) là phương pháp dị tìm điểm làm việc có cơng suất cực đại của hệ thống điện mặt trời thơng qua việc đóng mở khóa điện tử của bộ biến đổi DC/DC. Phương pháp MPPT được sử dụng rất phổ biến trong hệ thống PV làm việc độc lập và đang dần được áp dụng trong hệ quang điện làm việc với lưới.
được giới thiệu sơ ở phần dưới.
4.2.2 Bộ Buck-Boost
Cấu hình khơng cách ly của bộ Buck-Boost được minh họa ở hình 4.5. Bộ biến đổi bao gồm nguồn áp Vs, khóa điều khiển S, cuộn cảm L, diode D, tụ lọc đầu ra C và tải R. Khi khóa đóng, diode sẽ tắt (khơng dẫn), cuộn cảm nhận năng lượng từ nguồn, tụ sẽ cấp năng lượng cho tải. Khi khóa ngắt, diode sẽ dẫn, cuộn cảm sẽ cấp năng lượng cho tải và nạp năng lượng cho tụ.
Hình 4. 5 Mạch điện Buck-Boost
Tỷ số vào-ra của điện áp sẽ phụ thuộc vào tỷ số đóng cắt:
(4.1) Mơ tả sơ về dịng và áp của linh kiện trong bộ Buck-Boost khi khóa đóng cắt (lúc ổn định):
Hình 4. 6 Giản đồ hoạt động
Trong đó:
vL, iL: áp và dòng cuộn cảm vS, iS: áp và dòng nguồn vO, iO: áp và dòng ra
iC: dòng tụ
Khi thiết kế cần lưu ý giá trị cuộn cảm và tụ điện để đảm bảo dòng hoạt động liên tục và điện áp ra có độ nhấp nhơ bé:
(4.2) Vì sao bộ Buck-Boost có thể giữ được điểm cực đại của PV?
Thử xét tải trở ở đầu ra bộ Buck-Boost và PV ở đầu vào bộ Buck-Boost, khi hoạt động ổn định, phương trình cân bằng cơng suất vào ra:
(4.3) với Vpv và Ipv là áp và dòng PV, Vo và Io là áp và dòng tải.
đổi x. Dù trường hợp tải thay đổi hoặc điều kiện mơi trường thay đổi, cơng suất vẫn có thể giữ cực đại.
Hình 4. 7 Cách thay đổi đặc tuyến bộ Buck-Boost
4.3 Giới thiệu các giải thuật của MPPT
Trong khn khổ của luận văn này chỉ phân tích 2 phương pháp MPPT được ứng dụng rộng rãi và đã trở nên phổ biến, quen thuộc và cho được một số hiệu quả làm việc. Đó là Phương pháp nhiễu loạn và quan sát P&O, Phương pháp điện dẫn gia tăng INCond
4.3.1 Phương pháp nhiễu loạn và quan sát (Perturbation And Observe - P&O)
Trong thuật toán này điện áp hoạt động của pin mặt trời (PMT) bị nhiễu bởi một gia số nhỏ ΔV và kết quả làm thay đổi cơng suất, ΔP được quan sát. Hình 4.9 mơ tả ngun lý hoạt động của thuật tốn P&O, từ đó có thể suy ra cách thức hoạt động của thuật toán như sau:
- Nếu điểm hoạt động của hệ thống đang di chuyển theo hướng 1 (ΔP < 0 và ΔV < 0) thì cần tăng điện áp hoạt động lên để di chuyển điểm hoạt động tới điểm MPP.
- Nếu điểm hoạt động của hệ thống đang di chuyển theo hướng 2 (ΔP > 0 và ΔV > 0) thì cần tăng điện áp hoạt động lên để di chuyển điểm hoạt động tới điểm MPP.
- Nếu điểm hoạt động của hệ thống đang di chuyển theo hướng 3 (ΔP > 0 và ΔV < 0) thì cần giảm điện áp hoạt động để di chuyển điểm hoạt động tới điểm MPP. - Nếu điểm hoạt động của hệ thống đang di chuyển theo hướng 4 (ΔP < 0 và ΔV > 0) thì cần giảm điện áp hoạt động để di chuyển điểm hoạt động tới điểm MPP.
Hình 4. 8 Cách thực hiện P&O
Hình 4. 9 Ngun lý thực hiện thuật tốn P&O
Giải thích thuật tốn:
Bộ điều khiển MPPT sẽ đo các giá trị dòng điện I và điện áp V, sau đó tính tốn độ sai lệch ∆P, ∆V và kiểm tra:
Hình 4. 10 Lưu đồ P&O
Thuật tốn P&O có ưu điểm là đơn giản, dễ thực hiện. Khối lượng tính tốn ít. Trong tồn bộ giải thuật chỉ có phép tính cộng trừ và so sánh, khơng hề có phép tính chia.
4.3.2 Phương pháp điện dẫn gia tăng (Incremental Conductance-IncCond)
Phương pháp điện dẫn gia tăng (Incremental Conductance-IncCond) dựa trên việc phân tích đạo hàm của cơng suất theo điện áp. Ở bên trái điểm cực đại, đạo hàm có giá trị dương. Ở bên phải điểm cực đại đạo hàm có giá trị âm.
Giải thuật IncCond dựa trên việc so sánh những biểu thức trên để xác định điểm cực đại. Tỷ số độ thay đổi dòng điện và điện áp được so sánh với tỷ số dòng
điện và điện áp tức thời để biết được PV đang vận hành ở bên trái hay bên phải so với điểm cực đại.
(4.7)
Ta lại có:
(4.8) Vậy, cơng thức 4.1 trở thành:
(4.9)
Khi điểm vận hành của PV đến được điểm cực đại, giá trị đặt sẽ được giữ khơng đổi. Khi có sự thay đổi về điều kiện mơi trường và tải, giải thuật sẽ so sánh tiếp để điều chỉnh giá trị đặt nhằm bắt công suất cực đại mới.
Hình 4. 11 Ngun lý thực hiện thuật tốn InCond
CHƯƠNG 5: MƠ PHỎNG
5.1 Mơ hình mơ phỏng MPP trong Matlab Simulink
Hình 5. 1 Sơ đồ tổng quát hệ thống tìm MPP trong Matlab Simulink
Mơ hình được mơ phỏng trên thư viện simulink của phần mềm MATLAB 2014ra. Mô phỏng trên miền rời rạc, kiểu mô phỏng fix-step với bước tính Tstep = 2e-5 (s), phương pháp giải ode-4 Runge-Kutta.
Khối PV được xây dựng dựa trên S-Function.
Khối P&O và InCond xây dựng trong Embbeded Function. Giá trị các biến được đưa ra workspace để khảo sát và vẽ đồ thị.
Các giá trị đầu vào như nhiệt độ và bức xạ được thay đổi theo hàm nấc, tại các điểm thời gian 0.33s và 0.66s
Các thông số của PV: + Điện áp hở mạch: 30.8V + Dòng ngắn mạch: 8.7A + Số cell pin mắc nối tiếp: 50
Hình 5. 2 Sơ đồ chi tiết triển khai giải thuật P&O và InCond
Hình 5. 3 Mơ hình hóa PV BX1
Trong hình 5.4, ngõ ra của PI là tín hiệu PWM với tần số 10khz để khiển bộ Boost Converter. Các hệ số Kp, Ki được lựa chọn trong q trình mơ phỏng. Giá trị phù hợp là Kp = 1, Ki = 3.
Hình 5. 5 Mơ hình bộ boost converter
Hình 5. 6 Mơ hình MPPT
Trong hình 5.6, các giải thuật P&O và Incond sẽ được nhúng vào trong các khối Embbedded Function. Chi tiết được mô tả trong mục 5.2 và 5.3
Giải thuật P&O được thể hiện trong khối Embedded Function của Matlab Simulink:
function [Vref,I] = MPPT( V, P, Tsamp)
if isempty(Vref1) Vref1 = 0; end if isempty(V1) V1 = 0; end if isempty(P1) P1 = 0; end if isempty(N) N = 0; end
if (Tsamp > 0.001*N) %chu ky lay mau la 1ms
N = N + 1;
if((P - P1)> deltaP)%cong suat tang
if((V - V1) > deltaV) % dien ap tang --> ben trai mppt Vref = Vref1 + POV;
I = 1;
else %ben phai mppt Vref = Vref1 - POV; I = 2;
end
else %cong suat giam
Vref = Vref1 - POV; I = 3;
else %ben trai mppt Vref = Vref1 + POV; I = 4; end end P1 = P; V1 = V; Vref1 = Vref; else Vref = Vref1; I = 0; end end
Giải thuật IncCond được thể hiện trong khối Embedded Function của Matlab Simulink:
function [Vref,n] = MPPT( V, I, song)
%#eml persistent Vref1; persistent V1; persistent I1; Inc = 0.2 ; a = 0.2; if isempty(Vref1) Vref1 = 0; end if isempty(V1) V1 = 0; end if isempty(I1)
dV = 0; end if abs(dI)< 1e-3 dI = 0; end if song ==0 Vref = Vref1; n = 1; else
if dV == 0 %thay doi nho tren Vref, lam cho dien ap pin khong thay doi
%can duy tri do thay doi Vref cho den khi nao diep ap pin thay doi that %su
if dI == 0
Vref1 = Vref1; n = 2;
elseif dI > 0 %tren dac tuyen V-I, chieu search di len Vref1 = Vref1 + Inc;
n = 3;
else %tren dac tuyen V-I, chieu search di xuong Vref1 = Vref1 -Inc;
n = 4; end
else % dac tuyen P-I
Vref = Vref1; n = 5;
elseif dI/dV + I/V > a % ben trai mppt Vref1 = Vref1 + Inc;
n = 6; else
Vref1 = Vref1 - Inc; % ben phai mppt n = 7;
end
end
Vref = Vref1;
5.2 Mơ hình hóa và mơ phỏng PV
Mơ hình PV được xây dựng trên SimPowerSystem của Simulink để đảm bảo đặc tính điện. Mơ hình gồm có một khối (dùng Embedded function) để tính đầu ra IPV dựa theo VPV và điều kiện mơi trường. IPV tính được sẽ đưa vào một nguồn dòng điều khiển, giống như cấu trúc của PV theo cơng thức 3.4.
Hình 5. 7 Mơ hình trên Simulink
Trong đó Ns và Np chính là số panel PV nối tiếp và song song trong một mảng (array) PV.
Hình 5. 8 Đặc tuyến V-I khi nhiệt độ thay đổi, bức xạ khơng đổi
Hình 5. 10 Đặc tuyến P-I khi nhiệt độ thay đổi, bức xạ không đổi
Trong các hình 5.8 đến 5.10 mơ tả các mối quan hệ như V-I, P-V và P-I của tấm pin mặt trời khi nhiệt độ trên bề mặt tấm pin thay đổi ở các điểm như 25, 35, 45 và 55 độ C. Trong khảo sát này, bức xạ được giữ khơng đổi là 900W/m2.
Ở hình 5.8, ta thấy khi nhiệt độ thay đổi chỉ làm thay đổi áp hở mạch của pin, khơng ảnh hưởng tới dịng ngắn mạch. Chính vì vậy, điểm cơng suất cực đại có xu hướng dịch sang trái khi nhiệt độ tăng.
Ở hình 5.9, khi nhiệt độ tăng, cơng suất giảm. Có 4 điểm mà tại đó PV cho cơng suất cực đại tương ứng với 4 điều kiện hoạt động về nhiệt độ. Điện áp tại MPP giảm dần khi nhiệt độ tăng. Tại 55oC là 18.55V, công suất cực đại tương ứng là 132.7W. Tại 45oC là 19.54V, công suất cực đại tương ứng là 140.2W. Tại 35oC là 20.56V, công suất cực đại là 147.7W. Tại 25oC là 21.52V, công suất cực đại là 155.1W. Các giải thuật tìm kiếm MPP, phải tìm kiếm thành cơng các điểm này để ép PV luôn luôn cho công suất lớn nhất tùy thuộc vào từng điều kiện hoạt động.
Ở hình 5.10, nếu xét theo đồ thị P-I thì điểm cơng suất cực đại chỉ di chuyển theo trục tung khi nhiệt độ thay đổi. Thành phần dịng ngắn mạch ít thay đổi so với điện áp hở mạch.
Hình 5. 11 Đặc tuyến V-I khi bức xạ thay đổi, nhiệt độ khơng đổi
Hình 5. 13 Đặc tuyến P-I khi bức xạ thay đổi, nhiệt độ không đổi
Trong các hình 5.11 đến 5.13 mơ tả các mối quan hệ như V-I, P-V và P-I của tấm pin mặt trời khi bức xạ trên bề mặt tấm pin thay đổi ở các điểm như 600, 700, 800 và 900W/m2. Trong khảo sát này, nhiệt độ được giữ khơng đổi là 35 độ C.
Ở hình 5.11, ta thấy khi bức xạ thay đổi, nhiệt độ không đổi, áp hở mạch của các tấm pin không thay đổi nhiều nhưng dòng ngắn mạch thay đổi rất lớn. Điều này dẫn tới điểm cơng suất cực đại sẽ ít thay đổi theo hồnh độ.
Ở hình 5.12, ta thấy cơng suất thay đổi rất nhiều khi bức xạ thay đổi. Tại bức xạ 600W/m2, điện áp tại MPP là 20.94V, công suất cực đại là 101W. Tại bức xạ 700W/m2, điện áp tại MPP là 20.85V, công suất cực đại là 117.3W. Tại bức xạ 800W/m2, điện áp tại MPP là 20.7V, công suất cực đại là 132.7W. Tại bức xạ 900W/m2, điện áp tại MPP là 20.5V, công suất cực đại là 147.7W. Trong trường hợp bức xạ tăng dần, điện áp tại MPP không thay đổi nhiều. Các điểm MPP sẽ được gặp lại khi áp dụng các giải thuật tìm kiếm MPP.
Ở hình 5.13, quan sát trên biểu đồ P-I, ta thấy điểm công suất cực đại di