- Trang Chủ
- Năng lượng
- Xây dựng công cụ mô phỏng vị trí mặt trời và cường độ bức xạ tới trái đất theo thời gian và vị trí địa lý
Xem mẫu
- TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ NĂNG LƯỢNG - TRƯỜNG ĐẠI HỌC ĐIỆN LỰC
(ISSN: 1859 - 4557)
XÂY DỰNG CÔNG CỤ MÔ PHỎNG VỊ TRÍ MẶT TRỜI VÀ CƯỜNG ĐỘ
BỨC XẠ TỚI TRÁI ĐẤT THEO THỜI GIAN VÀ VỊ TRÍ ĐỊA LÝ
DEVELOPE A PROGRAM TO SIMULATE THE SUN MOTION AND RADIATION
RELATED TO DAY'S TIME AND GEOGRAPHIC LOCATION ON THE EARTH
Phạm Anh Tuân
Trường Đại học Điện lực
Ngày nhận bài: 04/01/2019, Ngày chấp nhận đăng: 28/03/2019, Phản biện: TS. Phạm Mạnh Hải
Tóm tắt:
Các ứng dụng năng lượng mặt trời ngày càng chiếm một vị trí quan trọng đối với cuộc sống, vì đây
là nguồn năng lượng sạch, bao phủ khắp Trái Đất và vô tận. Tính đến năm 2050 nhu cầu sử dụng
năng lượng từ điện mặt trời sẽ chiếm khoảng 20 % tổng nhu cầu điện năng trên toàn thế giới [1,2].
Năng lượng mặt trời không đồng nhất trên Trái Đất; nó phụ thuộc vào tọa độ, thời gian trong ngày,
ngày trong năm, dữ liệu khí tượng... Trong bài toán dự báo, bức xạ mặt trời có thể được tính trên cơ
sở là các phương trình mô tả chuyển động tương đối của Trái Đất - Mặt Trời, kết hợp với các thông
tin khác như vị trí địa lý, ngày trong năm, sự hấp thụ và phản xạ của các tầng mây... Trong bài báo
này, chúng tôi trình bày các phương trình thiên văn học mô tả chuyển động của Trái Đất - Mặt Trời,
qua đó lập trình để tính toán và mô phỏng nhằm đưa ra các đồ thị mô tả: toạ độ và góc phương vị
của Mặt Trời theo thời gian, năng lượng bức xạ chiếu đến Trái Đất trên nền matlab... Kết quả tính
toán và mô phỏng được ứng dụng để tính và đưa ra một số đồ thị kỹ thuật nhằm mô tả chuyển
động của Mặt Trời và lượng bức xạ tại các toạ độ khác nhau tại vị trí toạ độ bất kỳ trên Trái Đất và
là cơ sở cho việc chọn hướng lắp đặt dàn pin mặt trời hoặc các thiết bị thu năng lượng mặt trời và
dự báo tổng năng lượng sẽ thu được.
Từ khóa:
Bức xạ mặt trời, mô phỏng năng lượng mặt trời, điện mặt trời, Matlab/GUI.
Abstract:
Solar energy applications increasingly occupy an important position for life, as this is a clean energy,
renewable resource and covering all over the earth. By 2050, the demand for solar energy will
account for about 20% of the total electricity demand worldwide [1,2]. Solar energy is not uniform
on earth; it depends on coordinates, time of day, day of year, meteorological data, etc. For the
energy prediction, solar radiation can be calculated on the basis of equations describing relative
motion of the earth-sun, combined with other information such as geographic location, day of the
year, absorption and reflection of cloud layers, etc. In this paper, we present astronomical equations
describing the motion of the earth-sun, thereby programming to calculate and simulate to produce
descriptive graphs: coordinates and angles the position of the sun over time, radiant energy to the
earth under Matlab/GUI framework, etc. The simulation results are applied to calculate and give a
number of technical graphs to describe the sun's movement and the amount of radiation at different
coordinates at any position on earth and is the basis for choosing the direction of installing solar
panels or solar collectors and forecasting the total energy to be collected.
48 Số 19
- TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ NĂNG LƯỢNG - TRƯỜNG ĐẠI HỌC ĐIỆN LỰC
(ISSN: 1859 - 4557)
Keywords:
Solar radiation, solar energy simulation, solar electric system, Matlab/GUI.
1. GIỚI THIỆU CHUNG 8766 giờ. Quĩ đạo của Trái Đất quanh
Năng lượng mặt trời có thể chia thành 2 Mặt Trời tại các thời điểm trong năm
dạng gồm năng lượng bức xạ nhiệt và được mô tả như hình 1 [4-6].
năng lượng ánh sáng. Các ứng dụng của
năng lượng mặt trời bao gồm cả nhiệt
năng và điện năng ngày càng trở nên quan
trọng với cuộc sống. Điểm mạnh năng
lượng mặt trời so với nguồn năng lượng
khác là miễn phí, sạch và phong phú. Hơn
nữa, năng lượng mặt trời có thể khai thác
ở hầu hết các nơi điểm trên bề mặt trái
đất. Nó càng quan trọng hơn khi chi phí
Hình 1. Vị trí tương đối của Trái Đất và Mặt Trời
nhiên liệu hóa thạch cao và những tác
động làm ô nhiễm môi trường từ việc Đối với một điểm quan sát cụ thể trên
đốt chúng. Năng lượng mặt trời được Trái Đất, vị trí mặt trời có thể được xác
truyền đến trái đất thông qua bức xạ, định bởi hai góc, bao gồm: góc cao độ (α)
bao gồm hai phần: phần bức xạ nằm trên và góc phương vị (θs); các góc này được
bầu khí quyển và phần bức xạ dưới bầu mô tả như trong hình 2, trong đó:
khí quyển. Phần năng lượng bức xạ dưới
bầu khí quyển có thể được đo bằng Góc cao độ α (hay còn gọi là góc nhập xạ
các thiết bị đo năng lượng mặt trời như: - Altitude) là góc tạo bởi chùm tia sáng
nhật xạ kế (pyranometer, solarimeter), mặt trời so với mặt phẳng xích đạo có thể
trực xạ kế (pyrheliometer), nhật xạ ký được xác định bởi phương trình 1:
(actinography) [3]. Ngoài ra, nó cũng có sin sin . sin cos . cos . cos
thể được tính toán dựa trên các nguyên lý
thiên văn học và hình học, thông qua đó (1)
sẽ mô hình hóa nguồn mặt trời chiếu tới Với: L là vĩ độ tại vị trí quan sát, δ là góc
trái đất. Trong nghiên cứu này chúng tôi nghiêng tại vị trí quan sát (được tính theo
sẽ xây dựng công cụ để tính toán mô phương trình 2), và ω là góc giờ tại vị trí
phỏng vị trí và cường độ bức xạ năng quan sát (được tính theo phương trình 8).
lượng mặt trời theo thời gian và vị trí địa
lý trên nền Matlab/GUI. 23,45°. sin (2)
2. CƠ SỞ LÝ THUYẾT VÀ MÔ HÌNH ở đây N là ngày qui đổi trong năm, N được
TOÁN HỌC tính qui ước từ ngày 1 tháng 1 (N=1); các
Trái Đất xoay quanh Mặt Trời trong một ngày tiếp theo sẽ được tính tròn theo ngày,
quỹ đạo hình elip, thời gian của mỗi vòng ví dụ ngày 2 tháng 1 thì N=2, ngày 3 tháng
Trái Đất quay quanh Mặt Trời là khoảng 1 thì N=3…
Số 19 49
- TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ NĂNG LƯỢNG - TRƯỜNG ĐẠI HỌC ĐIỆN LỰC
(ISSN: 1859 - 4557)
Trái Đất nghiêng so với mặt phẳng quĩ 15°. (7)
đạo của nó quanh Mặt Trời; phương trình
Trong phương trình 7, hệ số 15° chính là
thực nghiệm để điều chỉnh độ lệch tâm
góc quay được của Trái Đất mỗi giờ (một
(EoT) do độ nghiêng trục này được tính
như sau (theo đơn vị phút): ngày Trái Đất quay được 360° trong vòng
24 giờ); LMST được hiểu là thời gian qui
đổi theo phút của múi giờ (giờ GMT).
9,87. sin 2 7,53. cos 1,5. sin
Trên cơ sở đó có thể tính được góc giờ tại
(3) vị trí quan sát:
ớ 81 (4) 15°. 12 (8)
Một số nhận xét: vào ngày 21/3 hoặc 4°.
23/9, Mặt Trời sẽ lên thiên đỉnh tại xích (9)
đạo và vào ngày 22/6 hoặc ngày 22/12,
Mặt Trời lên thiên đỉnh tại chí tuyến Bắc Theo qui ước này, vào buổi sáng góc giờ
và chí tuyến Nam. sẽ mang giá trị qui ước là âm và vào buổi
chiều giờ góc sẽ mang giá trị qui ước là
Góc phương vị (Azimuth) là góc giữa
dương.
vector bắc và mặt phẳng chiếu vuông góc
của Mặt Trời xuống mặt phẳng xích đạo; Năng lượng ánh sáng và cường độ bức xạ
góc phương vị có thể được tính theo công phát ra từ Mặt Trời sẽ tỉ lệ với bình
thức: phương khoảng cách tới điểm nhận bức
xạ. Cường độ bức xạ chiếu tới Trái Đất sẽ
sin (5)
tỉ lệ với bình phương khoảng cách thực
Hệ số hiệu chỉnh thời gian (tính bằng giữa Trái Đất và Mặt Trời (Gex) được tính
phút) tại vị trí khảo sát trong múi giờ (giờ bởi phương trình sau [5]:
GMT) của vị trí đó khi xét đến ảnh hưởng
của độ lệch tâm EoT được tính theo . (10)
phương trình 6 (với LOD là kinh độ của
trong ®ã Rav lμ kho¶ng c¸ch trung b×nh
vị trí khảo sát):
gi÷a MÆt Trêi vμ Tr¸i §Êt vμ R lμ kho¶ng
LMT = 4.(LOD – LSMT) + EoT (6) c¸ch thùc gi÷a MÆt Trêi vμ Tr¸i §Êt øng
víi thêi ®iÓm cÇn kh¶o s¸t. Theo c¸c sè
liÖu tÝnh to¸n thèng kª th× tØ lÖ kho¶ng
c¸ch thùc so víi kho¶ng c¸ch trung b×nh
sÏ phô thuéc vμo ngμy qui ®æi trong n¨m
vμ cã gi¸ trÞ kho¶ng:
1 0,0333. cos (11)
Thay thế phương trình 11 vào phương
trình 10 ta có:
Hình 2. Mô tả góc cao độ
và góc phương vị của Mặt Trời
. 1 0,0333. cos (12)
50 Số 19
- TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ NĂNG LƯỢNG - TRƯỜNG ĐẠI HỌC ĐIỆN LỰC
(ISSN: 1859 - 4557)
Gex chÝnh lμ cường ®é bøc x¹ n¨ng lượng ViÖc tÝnh to¸n phÇn n¨ng lượng t¸n x¹
mÆt trêi chiÕu xuèng Tr¸i §Êt. thùc tÕ rÊt khã so víi tÝnh to¸n bøc x¹ mÆt
trêi trùc tiÕp. PhÇn t¸n x¹ phô thuéc vμo
Cường ®é bøc x¹ cña MÆt Trêi trªn mÆt
líp khÝ quyÓn vμ tÇng m©y. Mét phÇn
®Êt (GT) được hiÓu lμ bøc x¹ cña MÆt Trêi
n¨ng lượng chiÕu xuèng Tr¸i §Êt sÏ bÞ
xuyªn qua bÇu khÝ quyÓn xuèng tíi ®é cao
ph¶n x¹ tõ bÒ mÆt trë l¹i bÇu trêi vμ mét
mùc nước biÓn. Bøc x¹ mÆt trêi trªn mÆt
phÇn t¸n x¹ xuèng mÆt ®Êt. Theo m« h×nh
®Êt bao gåm hai thμnh phÇn: phÇn chiÕu
®¼ng hướng (coi lượng t¸n x¹ ®Õn mét ®Þa
xuèng trùc tiÕp (GB) vμ phÇn t¸n x¹ bëi
®iÓm cã cường ®é b»ng nhau tõ mäi
c¸c ®¸m m©y vμ tÇng khÝ quyÓn (GD).
hướng), n¨ng lượng t¸n x¹ cã thÓ xÊp xØ
(GD).GT=GB+GD (13) b»ng:
360
N 100 GB (18)
Khi chïm bøc x¹ mÆt trêi ngoμi Tr¸i §Êt
GD 0,095 0,04sin
®i qua bÇu khÝ quyÓn, mét phÇn cña chïm 365
tia nμy bÞ hÊp thô bëi bÇu khÝ quyÓn.
Th«ng thường ®èi víi mét ngμy trêi 3. XÂY DỰNG CÔNG CỤ MÔ PHỎNG VỊ
quang, 70% bøc x¹ mÆt trêi xuèng Tr¸i TRÍ MẶT TRỜI VÀ CƯỜNG ĐỘ SÁNG
§Êt lμ bøc x¹ mÆt trêi trùc tiÕp. Theo m« TỚI TRÁI ĐẤT THEO THỜI GIAN VÀ VỊ
TRÍ ĐỊA LÝ
h×nh ASHRAE, bøc x¹ mÆt trêi trùc tiÕp
tíi bÒ mÆt tr¸i ®Êt (GB, nom) cã thÓ được Trên cơ sở lý thuyết đã trình bày ở trên,
tÝnh như sau: chúng tôi xây dựng thuật toán và lập trình
K xây dựng công cụ mô phỏng cường độ
GB ,nom Ae sin
(14) bức xạ năng lượng mặt trời theo thời gian
và vị trí địa lý trên Matlab/GUI. Công cụ
Trong đó A là một thông lượng biểu kiến này cho phép tính toán và mô phỏng vị trí
chiếu tới bề mặt trái đất và K là độ sâu tương đối của Trái Đất và Mặt Trời,
quang học (đại lượng này không thứ cường độ bức xạ của mặt trời trên bề mặt
nguyên). A và K phụ thuộc vào ngày qui trái đất... Sơ đồ khối của thuật toán này
đổi trong năm và có thể được tính theo được mô tả trong hình 3.
công thức thực nghiệm như sau:
Hình 4 là giao diện chính của công cụ mô
360
A 1160 75sin N 275 (15) phỏng vị trí và cường độ bức xạ của Mặt
365 Trời tới Trái Đất.
360
K 0,174 0, 035sin N 100 (16) D÷ liÖu ®Çu vμo gåm: täa ®é ®Þa lý, ngμy
365 th¸ng, mói giê, sè ®iÓm tÝnh trªn ®å thÞ...
Ngoμi ra ®Ó thuËn tiÖn trong sö dông,
Theo đó, cường độ bức xạ mặt trời trực
c«ng cô nμy cho phÐp người dïng chän
tiếp trên độ cao mực nước biển GB có thể
nhËp vÞ trÝ cÇn kh¶o s¸t theo mét sè ®Þa
được tính theo công thức 17.
danh được t¹o s½n trªn nÒn file excel, tiÕp
GB=GB,nomsin (17) ®ã d÷ liÖu sÏ tù ®éng cËp nhËp c¸c th«ng
Số 19 51
- TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ NĂNG LƯỢNG - TRƯỜNG ĐẠI HỌC ĐIỆN LỰC
(ISSN: 1859 - 4557)
tin vÒ kinh ®é vμ vÜ ®é, mói giê cña ®iÓm tÝnh được cμi ®Æt, sau ®ã khi tÝnh to¸n
cÇn kh¶o s¸t. c«ng cô cho phÐp xuÊt ®å thÞ m« pháng
(cöa sæ bªn tr¸i giao diÖn chÝnh) hoÆc xuÊt
D÷ liÖu ®Çu ra gåm: gãc cao ®é, gãc
d÷ liÖu dưới d¹ng sè.
phương vÞ, cường ®é bøc x¹ mÆt trêi t¹i
c¸c ®iÓm trªn bÒ mÆt khÝ quyÓn, t¹i c¸c Người dùng có thể khảo sát các dữ liệu
®iÓm ®o trªn Tr¸i §Êt. D÷ liÖu nμy được đầu ra hoặc dựng đồ thị trên cùng một
tÝnh theo tõng ngμy víi kho¶ng thêi gian trục trong nhiều ngày.
Hình 3. Sơ đồ khối thuật toán mô phỏng vị trí của Mặt Trời và bức xạ tới Trái Đất theo thời gian
Hình 4. Giao diện của phần mềm mô phỏng vị trí của Mặt Trời và bức xạ tới Trái Đất theo thời gian
52 Số 19
- TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ NĂNG LƯỢNG - TRƯỜNG ĐẠI HỌC ĐIỆN LỰC
(ISSN: 1859 - 4557)
4. KẾT QUẢ Quan s¸t trªn ®å thÞ thÊy gãc cao ®é lín
nhÊt (®iÓm Al-2) gÇn trïng víi gãc
phương vÞ b»ng 0 (®iÓm Az-2, MÆt Trêi
trïng víi phương nh×n hướng vÒ phÝa b¾c)
vμo thêi ®iÓm kho¶ng 11 giê 45 phót
trong ngμy (tøc lμ vμo kho¶ng gi¸ trÞ 705
phót qui ®æi øng víi trôc tung cña ®å thÞ).
Ngoμi ra đường phương vÞ cũng cho thÊy
thêi ®iÓm r¹ng s¸ng (®iÓm Al-, kho¶ng
thêi gian 380 phót quy ®æi - gãc cao ®é
vượt qua 0) còng lμ thêi ®iÓm gãc phương
vÞ øng víi gi¸ trÞ trªn trôc tung vμo
kho¶ng -65o (®iÓm Az-1); ®iÒu nμy cho
biÕt MÆt Trêi kh«ng ë hướng chÝnh ®«ng
Hình 5A. mô phỏng góc cao độ (Alitude angle)
và góc phương vị (Azimuth angle) tại thành mμ ë gãc kho¶ng 65o ®«ng nam. Gãc cao
phố Hà Nội ngày 1 tháng 2 ®é còng cho thÊy MÆt Trêi lÆn (®iÓm
A1-3, gãc cao ®é tiÕn tíi 0) vμo kho¶ng
Mét kÕt qu¶ m« pháng gãc cao ®é vμ gãc
1035 thêi ®iÓm phót quy ®æi, tøc kho¶ng
phương vÞ cña thμnh phè Hμ Néi (vÜ ®é
17 giê 15 phót. H×nh 5B minh häa vÞ trÝ
21o 02’ 00” b¾c, kinh ®é 105o 51’ 00”
cña MÆt Trêi nh×n tõ ®iÓm kh¶o s¸t.
®«ng) ngμy 1 th¸ng 12 như trong h×nh 5A;
trong ®ã trôc hoμnh lμ trôc thêi gian trong Một kết quả mô phỏng để so sánh về
ngμy qui ®æi theo phót (0 giê 00 phót sÏ cường độ bức xạ của Mặt Trời tại mặt đất
tương øng víi gi¸ trÞ phót 0), trôc tung lμ cho khu vực thành phố Hà Nội và khu vực
gi¸ trÞ gãc, tÝnh theo ®é (tÝnh chung cho c¶ thành phố Hồ Chí Minh (vĩ độ 10° 50’
gãc cao ®é vμ gãc phương vÞ). 00” bắc, kinh độ 106° 37’ 58” đông) vào
ngày 1 tháng 5 như hình 6. Kết quả này
được thực hiện mô phỏng trong điều kiện
bỏ qua sự suy giảm bức xạ do hấp thụ và
phản xạ của mây. Các kết quả trên hình
thể hiện mặt trời xuất hiện vào khoảng 6
giờ sáng (ứng với điểm qui đổi khoảng
360 phút), và lặn vào khoảng 18 giờ (ứng
với điểm qui đổi khoảng 1080 phút) và
cường độ bức xạ lớn nhất tại thời điểm
khoảng 12 h, trưa khi đó cường độ bức xạ
vào khoảng 1050 W/m2 đối với Hà Nội và
Hình 5B. Minh họa vị trí của Mặt Trời
1100 W/m2 đối với thành phố Hồ Chí
tại thành phố Hà Nội ngày 1 tháng 12 Minh).
Số 19 53
- TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ NĂNG LƯỢNG - TRƯỜNG ĐẠI HỌC ĐIỆN LỰC
(ISSN: 1859 - 4557)
Một kết quả khác tính toán mô phỏng bức
xạ mặt trời trong ngày mồng 1 tại các
tháng khác nhau từ tháng 1 đến tháng 6
tại thành phố Hà Nội như hình 8 (các
đường được đánh số từ 1 đến 6 ứng với
các tháng). Kết quả cho thấy cường độ
bức xạ thay đổi rõ rệt từ mức cường độ
bức xạ khoảng 740 W/m2 trong tháng 1 đã
tăng lên đến khoảng 1130 W/m2 trong
tháng 6. Ngoài ra thông qua các đường
Hình 6. Mô phỏng so sánh cường độ sáng bức xạ này ta cũng có thể dễ dàng quan
tại Hà Nội và Thành phố Hồ Chí Minh sát thấy thời gian xuất hiện cường độ bức
xạ (ban ngày) của các tháng mùa hè tăng
lên nhiều so với các tháng mùa đông.
Hình 7. Mô phỏng bức xạ mặt trời ngày 18 trên
tại bề mặt khí quyển và Trái Đất tại Hà Nội
Hình 8. Mô phỏng bức xạ mặt trời
KÕt qu¶ tÝnh to¸n m« pháng kÕt qu¶ so trong các ngày mồng 1 của các tháng 1
s¸nh cường ®é bøc x¹ t¹i bÒ mÆt khÝ đến tháng 6 tại Hà Nội
quyÓn vμ bÒ mÆt tr¸i ®Êt t¹i Hμ Néi ngμy 1 5. KẾT LUẬN
th¸ng 8 như trong h×nh 7. M« pháng nμy
cho thÊy t¹i thêi ®iÓm cường ®é bøc x¹ Các thiết bị đo bức xạ mặt trời như
lín nhÊt vμo kho¶ng 12 giê. Đường phÝa nhật xạ kế (pyranometer, solarimeter),
trªn (đường xanh) lμ cường ®é bøc x¹ t¹i trực xạ kế (pyrheliometer), nhật xạ ký
bÒ mÆt khÝ quyÓn; gi¸ trÞ lín nhÊt vμo (actinography) thường được lắp đặt tại
kho¶ng 1330 W/m2. Đường phÝa dưới các địa điểm được lựa chọn ở các vùng cụ
(đường mμu vμng) lμ cường ®é bøc x¹ ë thể. Trên thực tế, để thực hiện đo và tổng
bÒ mÆt tr¸i ®Êt, gi¸ trÞ lín nhÊt chØ cßn hợp dữ liệu đo này có chi phí khá lớn nên
kho¶ng 1030 W/m2 (do bÞ suy gi¶m bëi không thể đặt chúng tại mọi điểm trên
khÝ quyÓn). Trái Đất. Các tính toán này cũng được
54 Số 19
- TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ NĂNG LƯỢNG - TRƯỜNG ĐẠI HỌC ĐIỆN LỰC
(ISSN: 1859 - 4557)
dùng để dự báo năng lượng mặt trời bức nhiều ngày; mô phỏng này giúp người
xạ ở mọi địa điểm trên Trái Đất. Nó có ý dùng dễ dàng so sánh, phân tích các kết
nghĩa đặc biệt trong việc xác định dữ liệu quả. Kết quả này là cơ sở cho việc tính
tại những vị trí không có hoặc chưa thể chọn hướng lắp đặt các thiết bị khai thác
lắp đặt được thiết bị đo năng lượng mặt năng lượng mặt trời (dàn pin mặt trời,
trời. thiết bị thu nhiệt…) ở một vị trí bất kỳ để
đảm bảo năng lượng bức xạ chiếu tới nó
Các giá trị mà công cụ mô phỏng này tính là lớn nhất. Ngoài ra, các kết quả mô
được mặc dù không mới, tuy nhiên kết phỏng cũng là cơ sở cho việc tích phân để
quả mô phỏng giúp người dùng nhanh dự báo tổng năng lượng chiếu đến một
chóng có được đồ thị mô phỏng trong điểm bất kỳ.
TÀI LIỆU THAM KHẢO
[1] https://www.equinor.com/en/how-and-why/sustainability/energy-perspectives.html, 2018.
[2] IRENA, Global Energy transformation: A roadmap to 2050 (2019 edition). International
Renewable Energy Agency, 2019.
[3] S. Radiation and D. Models, Solar Radiation and Daylight Models. Elsevier, 2004.
[4] V. Badescu, Modeling Solar Radiation at the Earth’s Surface. Spinger, 2008.
[5] T. Khatib, Modeling of photovoltaic systems using MATLAB. Wiley, 2016.
[6] “https://www.pveducation.org/pvcdrom/terrestrial-solar-radiation,” 2018.
Giới thiệu tác giả:
Tác giả Phạm Anh Tuân nhận bằng Thạc sĩ ngành kỹ thuật điện năm 2006,
nhận bằng Tiến sĩ ngành khoa học vật liệu năm 2017 tại Trường Đại học
Bách khoa Hà Nội. Tác giả hiện là giảng viên Khoa Kỹ thuật điện, Trường
Đại học Điện lực.
Lĩnh vực nghiên cứu: vật liệu pin mặt trời và mô phỏng hệ thống điện mặt trời.
Số 19 55
- TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ NĂNG LƯỢNG - TRƯỜNG ĐẠI HỌC ĐIỆN LỰC
(ISSN: 1859 - 4557)
46 Số 19
nguon tai.lieu . vn