Tổng quan nghiên cứu

Năng lượng tái tạo đang trải qua giai đoạn tăng trưởng vượt bậc trên quy mô toàn cầu. Theo các báo cáo năng lượng quốc tế, năm 2016 đã chứng kiến công suất năng lượng tái tạo lắp đặt mới đạt mức kỷ lục 161 GW với tổng vốn đầu tư lên đến 249,8 tỷ USD. Trong đó, điện mặt trời quang điện (PV) khẳng định vị thế dẫn đầu khi bổ sung hơn 98 GW công suất trong năm 2017, nâng tổng công suất tích lũy toàn cầu vượt mốc 402 GW, tương ứng mức tăng trưởng ấn tượng 33%. Tại Việt Nam, tiềm năng bức xạ rất dồi dào với số giờ nắng trung bình tại miền Bắc đạt khoảng 1.100 giờ/năm và miền Nam đạt trên 2.100 giờ/năm, mở ra triển vọng to lớn để hiện thực hóa mục tiêu công suất điện mặt trời đạt 250 MWp vào năm 2025.

Mặc dù có tốc độ phát triển mạnh mẽ, việc khai thác tối đa công suất từ các hệ thống pin mặt trời vẫn đối mặt với thách thức kỹ thuật lớn do tính phi tuyến sâu sắc của đặc tuyến dòng điện - điện áp (I-V) và công suất - điện áp (P-V). Trong thực tế vận hành, hiệu suất phát điện không chỉ phụ thuộc vào công nghệ chế tạo bán dẫn mà còn biến thiên liên tục theo cường độ bức xạ và nhiệt độ môi trường. Các mô hình toán học truyền thống thường giả định các tham số cấu trúc, đặc biệt là điện trở nối tiếp và điện trở song song, là các đại lượng cố định ở điều kiện thử nghiệm tiêu chuẩn. Điều này dẫn đến sai số mô phỏng tăng cao khi điều kiện thời tiết thay đổi.

Luận văn thạc sĩ chuyên ngành Kỹ thuật điện của tác giả Võ Thị Ánh Tuyết, thực hiện dưới sự hướng dẫn khoa học của Tiến sĩ Nguyễn Trung Nhân tại Trường Đại học Công nghiệp Thành phố Hồ Chí Minh, tập trung giải quyết triệt để hạn chế trên. Mục tiêu cốt lõi của nghiên cứu là xây dựng mô hình toán học động cho hệ thống pin mặt trời, có khả năng tự động điều chỉnh linh hoạt các tham số điện trở ký sinh theo điều kiện vận hành thực tế thông qua giải thuật khớp ba điểm đặc trưng: ngắn mạch, hở mạch và điểm công suất cực đại. Nghiên cứu mang lại ý nghĩa thực tiễn to lớn khi giúp giảm sai số dự báo công suất xuống dưới mức 1,2%, hỗ trợ tối ưu hóa thuật toán bám điểm công suất cực đại (MPPT) và nâng cao hiệu quả vận hành kinh tế cho các dự án điện mặt trời.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên nền tảng lý thuyết quang điện bán dẫn và mô hình mạch điện tương đương phi tuyến. Hiệu ứng quang điện trong pin mặt trời là quá trình biến đổi trực tiếp năng lượng photon ánh sáng thành năng lượng điện thông qua sự giải phóng các hạt mang điện tại lớp tiếp giáp p-n. Mô hình toán học tổng quát của pin quang điện được phát triển từ phương trình diode bán dẫn Shockley kết hợp định luật dòng điện Kirchhoff.

Khung lý thuyết của đề tài tập trung phân tích và so sánh các mô hình mạch tương đương từ cấp độ đơn giản đến nâng cao:

  • Mô hình lý tưởng một diode: Bao gồm một nguồn dòng quang điện tỷ lệ thuận với bức xạ mặt trời mắc song song với một diode lý tưởng mô tả hiện tượng phân cực.
  • Mô hình thực tế có tính đến tổn thất thuần trở: Bổ sung điện trở nối tiếp đại diện cho điện trở khối bán dẫn, điện trở tiếp xúc kim loại, cùng điện trở song song phản ánh dòng rò qua lớp tiếp giáp p-n.
  • Mô hình hai diode: Mở rộng để xem xét tổn thất tái hợp vùng điện tích không gian, phục vụ việc đối chiếu độ chính xác.

Nghiên cứu tập trung vào năm khái niệm chuyên ngành then chốt bao gồm dòng quang điện, dòng điện bão hòa ngược của diode, điện áp hở mạch, dòng điện ngắn mạch và điểm công suất cực đại. Phương trình đặc tuyến dòng điện - điện áp của mô hình một diode cải tiến được thiết lập dưới dạng phương trình siêu việt phi tuyến, trong đó dòng điện đầu ra phụ thuộc chặt chẽ vào điện áp đầu cực, nhiệt độ tiếp giáp và hệ số lý tưởng của diode.

Phương pháp nghiên cứu

Nghiên cứu sử dụng phương pháp mô phỏng số học kết hợp kiểm chứng thực nghiệm dựa trên cơ sở dữ liệu kỹ thuật từ nhà sản xuất thương mại:

  • Cỡ mẫu và đối tượng nghiên cứu: Dữ liệu kiểm chứng được thu thập từ thông số kỹ thuật tiêu chuẩn của 2 dòng tấm pin thương mại phổ biến trên thị trường là Kyocera KC200GT (công suất định danh 200W, gồm 54 tế bào quang điện đa tinh thể mắc nối tiếp) và Solarex MSX60 (công suất định danh 60W, gồm 36 tế bào quang điện). Ngoài ra, mô hình còn sử dụng tập dữ liệu đối chiếu từ tấm pin Siemens SM110-24.
  • Phương pháp chọn mẫu: Phương pháp lấy mẫu có chủ đích được áp dụng nhằm lựa chọn các chủng loại pin mặt trời đại diện cho công nghệ silic tinh thể chiếm trên 85% thị phần thương mại, đảm bảo tính phổ quát cho giải thuật.
  • Phương pháp phân tích: Tác giả phát triển giải thuật lặp số học trên môi trường phần mềm MATLAB/Simulink. Lý do lựa chọn phương pháp lặp giải tích kết hợp mô phỏng số là vì phương trình đặc tính I-V là phương trình ẩn phi tuyến, không thể giải trực tiếp bằng phương pháp đại số thông thường. Giải thuật lặp cho phép liên tục tinh chỉnh giá trị điện trở nối tiếp và điện trở song song để đồ thị công suất đạt cực đại chính xác tại điểm công suất cực đại danh định từ bảng dữ liệu.
  • Kế hoạch triển khai: Toàn bộ quá trình nghiên cứu, xây dựng thuật toán và mô phỏng kiểm chứng được thực hiện hoàn chỉnh trong thời gian 6 tháng.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình mô phỏng và kiểm chứng trên phần mềm MATLAB đã mang lại những phát hiện then chốt về tính chính xác và độ ổn định của mô hình động:

Thứ nhất, ở điều kiện đo tiêu chuẩn với cường độ bức xạ 1000 W/m2 và nhiệt độ cell 25 độ C, mô hình đề xuất đạt độ khớp gần như tuyệt đối với thông số thực tế của tấm pin Kyocera KC200GT. Công suất cực đại mô phỏng đạt 200,14 W so với giá trị danh định 200,00 W, sai số công suất chỉ ở mức 0,07%. Trong khi đó, các mô hình tham số tĩnh truyền thống có mức sai số dao động từ 2,5% đến 3,8%.

Thứ hai, khi nhiệt độ làm việc của tấm pin tăng cao từ 25 độ C lên 75 độ C dưới bức xạ 1000 W/m2, điện áp hở mạch suy giảm mạnh theo hệ số nhiệt độ âm khoảng -0,35%/K. Mô hình động đề xuất đã tự động điều chỉnh các thông số cấu trúc, duy trì sai số tương đối của đường cong I-V và P-V dưới 1,15%. Ngược lại, mô hình cố định thông số ghi nhận mức sai lệch vượt quá 4,80% tại vùng lân cận điểm cực đại.

Thứ ba, dưới điều kiện bức xạ mặt trời thay đổi từ 200 W/m2 đến 1000 W/m2 ở nhiệt độ không đổi 25 độ C, dòng điện ngắn mạch biến thiên tỷ lệ thuận với cường độ sáng. Giải thuật lặp đã xác lập giá trị điện trở song song tăng phi tuyến khi bức xạ giảm, giúp hạn chế hiện tượng sụt áp và kiểm soát sai số tính toán công suất dưới 0,95%.

Thứ tư, khi kiểm chứng trên mô-đun Solarex MSX60 công suất 60W, mô hình động tiếp tục thể hiện tính vượt trội với sai số trung bình giảm hơn 62% so với các mô hình mạch chuẩn trước đây, khẳng định khả năng thích ứng linh hoạt trên nhiều dải công suất khác nhau.

Thảo luận kết quả

Các kết quả nghiên cứu đã phản ánh chính xác bản chất vật lý của vật liệu bán dẫn trong điều kiện môi trường thực tế. Điện trở nối tiếp chủ yếu đại diện cho điện trở tiếp xúc của các điện cực kim loại và điện trở nội của lớp bán dẫn, vốn chịu ảnh hưởng mạnh bởi nhiệt độ tiếp giáp. Trong khi đó, điện trở song song phản ánh dòng điện rò rỉ qua các mép của tế bào quang điện và sự tái hợp hạt mang điện, đại lượng này biến đổi rõ rệt khi mật độ bức xạ photon thay đổi.

Việc giả định hai tham số này là hằng số như trong các nghiên cứu trước đây là nguyên nhân chính gây ra sự sai lệch đáng kể ở vùng điện áp làm việc cao và mức chiếu xạ thấp. Giải thuật lặp do tác giả đề xuất đã giải quyết triệt để nhược điểm này bằng cách sử dụng ba điểm dữ liệu chuẩn từ nhà sản xuất để liên tục cập nhật giá trị điện trở ký sinh phù hợp với trạng thái nhiệt - quang thực tế.

Toàn bộ dữ liệu mô phỏng và thực nghiệm được biểu diễn trực quan qua hệ thống đồ thị đặc tuyến I-V và P-V đa biến, cùng các bảng biểu phân tích sai số chi tiết tại các mức nhiệt độ 25 độ C, 50 độ C và 75 độ C. Kết quả cho thấy đường cong mô phỏng từ mô hình đề xuất bám sát hoàn toàn các điểm đo thực nghiệm, cung cấp công cụ toán học tin cậy để nghiên cứu và thiết kế các hệ thống chuyển đổi năng lượng mặt trời.

Đề xuất và khuyến nghị

Dựa trên kết quả nghiên cứu đạt được, tác giả đưa ra 4 nhóm giải pháp kỹ thuật có tính ứng dụng cao:

  1. Tích hợp giải thuật thông số động vào bộ điều khiển MPPT: Các kỹ sư thiết kế phần cứng và lập trình nhúng cần tiến hành tích hợp thuật toán lặp xác định điện trở ký sinh vào firmware của các bộ biến tần hòa lưới và biến tần độc lập. Mục tiêu là nâng cao hiệu suất bám điểm công suất cực đại thêm từ 3% đến 5% trong điều kiện bức xạ thay đổi đột ngột, hoàn thành thử nghiệm trong vòng 6 tháng tới.
  2. Xây dựng khối thư viện mô phỏng chuẩn hóa trên MATLAB/Simulink: Các trường đại học kỹ thuật và viện nghiên cứu năng lượng cần phát triển bộ công cụ mô phỏng mã nguồn mở dựa trên mô hình động này. Giải pháp này giúp các kỹ sư hệ thống rút ngắn khoảng 40% thời gian tính toán, phân tích và tối ưu hóa chuỗi pin quang điện trước khi lắp đặt thực tế.
  3. Ban hành quy trình chuẩn đoán và đo kiểm suy hao tham số định kỳ: Doanh nghiệp vận hành nhà máy điện mặt trời cần xây dựng quy trình đo quét đường cong đặc tuyến I-V định kỳ 6 tháng một lần. Việc đối chiếu các tham số điện trở nối tiếp và song song thực tế với mô hình chuẩn sẽ giúp phát hiện sớm hiện tượng thoái hóa tấm pin, kiểm soát tỷ lệ suy giảm công suất dưới mức 0,7%/năm.
  4. Cải tiến cấu trúc thông gió và tản nhiệt cho dàn pin mặt trời: Đơn vị thi công cơ điện EPC cần tối ưu khoảng cách lắp đặt giữa tấm pin và mái công trình để tăng cường lưu thông không khí, phấn đấu duy trì nhiệt độ cell dưới 45 độ C vào mùa nắng nóng, hạn chế mức sụt giảm công suất do hiệu ứng nhiệt độ từ 8% đến 12%.

Đối tượng nên tham khảo luận văn

Nội dung và giải thuật trong luận văn mang lại giá trị học thuật và ứng dụng thực tiễn sâu sắc cho 4 nhóm đối tượng:

  • Kỹ sư thiết kế và vận hành hệ thống điện mặt trời (EPC Solar Engineers): Tài liệu cung cấp cơ sở toán học chuẩn xác để mô phỏng sản lượng phát điện hàng năm, tính toán tổn thất chuỗi pin và tối ưu hóa cấu hình ghép nối biến tần trong các dự án áp mái và trang trại điện mặt trời.
  • Giảng viên và nhà nghiên cứu chuyên ngành Kỹ thuật điện, Năng lượng tái tạo: Luận văn là tài liệu tham khảo chuyên sâu về mô hình hóa phi tuyến phần tử bán dẫn quang điện, phục vụ công tác giảng dạy các môn học về biến đổi năng lượng tái tạo và phát triển đề tài nghiên cứu liên quan.
  • Học viên cao học và sinh viên kỹ thuật: Cung cấp phương pháp luận mẫu mực về quy trình giải thuật toán lặp số học trên MATLAB/Simulink, kỹ năng phân tích đặc tuyến phi tuyến và cách tiếp cận giải quyết bài toán kỹ thuật từ bảng dữ liệu thương mại.
  • Nhà phát triển thuật toán điều khiển công suất: Tham khảo mô hình toán học chính xác để kiểm thử, đánh giá và nâng cấp các thuật toán bám điểm công suất cực đại như thuật toán xáo trộn và quan sát (P&O), điện dẫn gia tăng (Incremental Conductance) trong môi trường mô phỏng thời gian thực.

Câu hỏi thường gặp

Tại sao cần xây dựng mô hình động thay vì sử dụng các mô hình tham số cố định?

Trong thực tế, điện trở nối tiếp và điện trở song song của tấm pin biến thiên liên tục theo nhiệt độ môi trường và cường độ bức xạ. Các mô hình cố định từ nhà sản xuất chỉ chính xác ở điều kiện tiêu chuẩn 25 độ C và 1000 W/m2, dẫn đến sai số mô phỏng công suất tăng lên từ 3% đến 5% khi thời tiết thay đổi. Mô hình động giúp triệt tiêu sai số này, duy trì độ chính xác dưới 1,2%.

Hai tham số điện trở nối tiếp và điện trở song song ảnh hưởng như thế nào đến đặc tuyến pin?

Điện trở nối tiếp làm giảm độ dốc của đường cong đặc tính I-V tại vùng điện áp hở mạch, trực tiếp gây sụt giảm điện áp đầu ra và công suất cực đại khi dòng điện tăng. Trong khi đó, điện trở song song nhỏ làm tăng độ dốc tại vùng dòng ngắn mạch, gây rò rỉ dòng điện qua lớp bán dẫn và làm suy giảm nghiêm trọng hệ số lấp đầy của pin mặt trời.

Giải thuật đề xuất yêu cầu những thông số đầu vào nào từ bảng dữ liệu thương mại?

Giải thuật chỉ yêu cầu các thông số kỹ thuật tiêu chuẩn sẵn có trên bảng dữ liệu của nhà sản xuất bao gồm: dòng điện ngắn mạch, điện áp hở mạch, dòng điện tại điểm công suất cực đại, điện áp tại điểm công suất cực đại, cùng các hệ số suy giảm nhiệt độ của dòng điện và điện áp mà không cần tiến hành đo đạc phức tạp trong phòng thí nghiệm.

Mô hình động này có thể áp dụng cho các công nghệ pin quang điện mới không?

Mô hình toán học hoàn toàn có thể áp dụng cho nhiều công nghệ pin mặt trời khác nhau như Silic đơn tinh thể, đa tinh thể và pin màng mỏng. Bằng cách cập nhật các giá trị đặc trưng hở mạch, ngắn mạch và điểm cực đại tương ứng từ bảng thông số của từng công nghệ, giải thuật lặp sẽ tự động tìm ra cặp giá trị điện trở ký sinh tối ưu nhất.

Mô hình toán học này hỗ trợ tối ưu hóa thuật toán MPPT như thế nào?

Mô hình động cung cấp đặc tuyến công suất P-V chính xác tức thời theo điều kiện bức xạ và nhiệt độ thực tế. Điều này cho phép các kỹ sư mô phỏng chính xác hành vi bám điểm cực đại của bộ điều khiển MPPT trong môi trường số, hạn chế hiện tượng dao động quanh điểm cực đại và giảm tổn thất năng lượng chuyển đổi từ 2% đến 4%.

Kết luận

  • Luận văn đã xây dựng thành công mô hình toán học động cho hệ thống pin mặt trời với khả năng tự điều chỉnh điện trở ký sinh theo điều kiện vận hành.
  • Giải thuật lặp số học đề xuất đã triệt tiêu sự phụ thuộc vào các tham số cố định, đảm bảo công suất cực đại mô phỏng luôn trùng khớp với dữ liệu thực nghiệm.
  • Kết quả kiểm chứng trên các tấm pin Kyocera KC200GT và Solarex MSX60 khẳng định mô hình duy trì sai số công suất dưới 1,2% trên toàn bộ dải nhiệt độ từ 25 độ C đến 75 độ C.
  • Đóng góp chính của công trình là cung cấp một giải pháp mô hình hóa chính xác cao, dễ dàng tích hợp vào phần mềm mô phỏng mạch điện và hệ thống điều khiển công nghiệp.
  • Trong giai đoạn tiếp theo, nghiên cứu cần được mở rộng để đánh giá hiện tượng bóng che cục bộ và phát triển thuật toán bám điểm cực đại toàn cục trong thời gian thực.

Quý độc giả, kỹ sư và nhà nghiên cứu quan tâm có thể ứng dụng ngay mô hình động này vào các dự án thiết kế và mô phỏng hệ thống điện mặt trời để nâng cao độ chính xác và hiệu quả kinh tế cho công trình.