Tính toán bù tối ưu cho lưới điện phân phối có xét đến sự phát triển của phụ tải

Nội dung text: Tính toán bù tối ưu cho lưới điện phân phối có xét đến sự phát triển của phụ tải

1. TÍNH TOÁN BÙ TỐI ƯU CHO LƯỚI ĐIỆN PHÂN PHỐI CÓ XÉT ĐẾN SỰ PHÁT TRIỂN CỦA PHỤ TẢI DISTRIBUTION OF ELECTRICAL NETWORK WITH THE OPTIMUM CACULATION OF COMPENSATION WILL BE CONSIDERED THE DEVELOPMENT OF LOAD POWER Lê Kỷ(1); Đinh Công Thành(2) (1): Trường Đại Học Bách Khoa Tp. Hồ Chí Minh. (2): Trường Đại Học Sư Phạm Kỹ Thuật Tp. Hồ Chí Minh. TÓM TẮT Lưới điện phân phối (LĐPP) liên tục phát triển theo không gian và thời gian, ngoài việc đáp ứng đủ nhu cầu công suất tác dụng cho phụ tải, cần chú trọng đáp ứng nhu cầu công suất phản kháng (CSPK). Việc bù CSPK sẽ cho phép nâng cao chất lượng điện năng cũng như hiệu quả kinh tế của LĐPP. Luận văn trình bày kết quả nghiên cứu xây dựng cơ sở dữ liệu và các chỉ số về điều kiện khí hậu LĐPP cho chương trình PSS/ADEPT. Từ đó sử dụng chương trình để tính toán, đánh giá hiệu quả các phương án bù tối ưu phía trung áp và phía hạ áp, nhằm đề xuất phương án bù tối ưu cho LĐPP. Từ khóa: Lưới điện phân phối, dự báo phụ tải điện, dự báo phụ tải ngắn hạn, ABSTRACT The distribution network developed continuously in space and time, beside sufficiently meeting the active power need for additional load, it’s also very important to meet reactive power need. Reactive power compensation will enhance electricity quality as well as economic efficiency of the distribution network. In this article, we present the study result of establishing data bases and economic indexes of the distribution network to PSS/ADEPT. Then, we can use this program to calculate and evaluate the efficiency of medium and low voltage optimal compensation with the aim of proposing the most optimal alternative for the distribution network. Key: distribution line, electric load forecasting, short-term load forecasting models, 1. Giới thiệu điểm bù công suất lưới điện phân phối có sự phát Ngày nay, với sự phát triển nhanh chóng của nền triển của phụ tải trong tương lai. kinh tế và đời sống nhân dân được nâng cao, do đó 2. Lợi ích bù công suất phản kháng nguồn điện cũng phải đáp ứng được những đòi hỏi về Việc bù Công suất phản kháng cho lưới điện sẽ công suất và chất lượng được đảm bảo. Vấn đề công có những tích cực như sau: suất phát ra phải được truyền tải và tận dụng một cách a. Giảm được tổn thất công suất trong mạng hiệu quả nhất, không để lãng phí quá nhiều ảnh điện hưởng đến kinh tế là một bài toán được rất nhiều đề Ta có tổn thất công suất trên đường dây trước tài nghiên cứu. Tổn hao công suất là vấn đề ảnh khi bù được xác định: hưởng đến chất lượng điện năng, do đó cần phải có 푃2+푄2 ∆푃 = 2 ∙ 푅 (1) các biện pháp đưa ra để giảm lượng tổn thất công suất 푈đ trên lưới. Cùng với đó là sự phát triển của phụ tải điện Tổn thất công suất trên đường dây sau khi bù: 푃2 + (푄−푄 )2 trong tương lại, việc dự báo sự tăng trưởng đúng đắng ∆푃 = ∙ 푅 (2) 푈2 sự phát triển của phụ tải sẽ góp phần vào việc giảm đ Khi giảm Q truyền tải trên đường dây, ta giảm chi phí đầu tư ban đầu và giảm được chi phí bán điện được thành phần tổn thất công suất Q gây ra. cho khách hàng. Vì vậy, cần có một biện pháp là tìm
2. b. Giảm được tổn thất điện áp trong mạng Bảng 3.1: Thời gian vận hành các chế độ vận hành điện Thời gian tồn tại Thời Hệ số phụ Tổn thất điện áp trước khi bù được xác định theo (p.u*) Chế độ gian tải công thức: (Relative 푃.푅+푄. (giờ) (p.u ) ∆푈 = (3) duration) 푈 Thấp 21h – Tổn thất điện áp sau khi bù: 9/24 = 0.3750 0.68 điểm 5h 푃.푅 + (푄−푄 .) ∆푈 = (4) Bình 06h – 푈 11/24 = 0.4583 0.8 Khi ta giảm Q truyền tải trên đường dây sẽ giảm thường 16h 17h – được tổn thất điện áp trên đường dây. Cao điểm 4/24 = 0.1667 1 c. Tăng khả năng truyền tải của đường dây 20h và máy biến áp * 1 p.u = 24h, 1p.u = 100% (full load) Khả năng truyền tải của đường dây và máy biến Khi đó việc tối ưu dung lượng bù cần xem xét áp phụ thuộc vào điều kiện phát nóng, tức phụ thuộc cả 3 chế độ vận hành của xuất tuyến 477AX. Giả sử vào dòng điện cho phép của chúng. Dòng điện chạy ta không giới hạn về số lượng, dung lượng tụ bù trung trên dây dẫn và máy biến áp được tính như sau: thế ứng động và tụ trung thế cố định (trường hợp mỗi √푃2+푄2 = (5) bộ tụ có dung lượng 300 kVar) ta dùng PSS tính toán √3푈 bù tối ưu ở trung thế, tìm tổng P, Q tiết kiệm được Từ công thức (5) cho thấy với cùng một tình sau bù và tính ra lợi nhuận sau khi bù cho 3 chế độ trạng phát nóng nhất định của đường dây và máy biến vận hành tải khác nhau. áp (tức I = const) chúng ta có thể tăng khả năng Như vậy, khi xem xét lưới điện vận hành ở chế truyền tải Công suất tác dụng bằng cách giảm Công độ khác nhau (Cao điểm, bình thường và thấp điểm) suất phản kháng mà chúng phải tải đi. Vì thế khi vẫn việc tối ưu vị trí và dung lượng tụ bù trên xuất tuyến giữ nguyên đường dây và máy biến áp, nếu giảm 477AX với kết quả là: tổng chi phí đầu tư lắp đặt thiết lượng Q phải truyền tải thì khả năng truyền tải của bị bù ban đầu là 142,075.57 $, sau 8 năm chi phí tiết chúng sẽ được tăng lên, góp phần làm ổn định điện kiệm được trên hệ thống sau khi khấu trừ các chi phí áp, tăng khả năng phát điện của máy phát điện lạm phát, trượt giá thì số tiền tiết kiệm được là Việc bù Công suất phản kháng ngoài việc nâng 324,615.60 $. cao hệ số công suất cosφ còn đưa đến hiệu quả là Trường hợp 2: Xem xét lắp đặt thiết bị bù giảm được chi phí kim loại màu tức giảm được tiết trong điều kiện vị trí bù phù hợp trong thực tế diện dây dẫn nên tiết kiệm được chi phí đầu tư xây Dựa trên lần chạy thứ 1 tìm vị trí và dung lượng dựng lưới điện. Giảm được chi phí điện năng bù tối ưu có xét đến 3 trường hợp phụ tải cao điểm, 3. Tính toán bù CSPK bằng chương trình bình thường và thấp điểm của một ngày trung bình. PSS/ADEPT cho phát tuyến 477AX Điện lực Cà Mau đã khảo sát vị trí bù theo các điểm Thực hiện mô phỏng các xuất tuyến 477AX vào bù của lần chạy thứ 1. Qua thực tế cho thấy tại một chương trình PSS/ADEPT, tiến hành tính phân bố số vị trí như nút 13,14,19,47 bị hạn chế cho việc lắp công suất, tính toán các phương án bù cho xuất tuyến đặt tụ bù. 477AX theo các trường hợp sau: Như vậy theo như bảng kết quả như trên xuất Trường hợp 1: Không giới hạn dung lượng tụ tuyến 477AX có thêm 2 vị trí bù mới tại nút 12,32 bù và vị trí bù trên xuất tuyến 477AX với dung lượng tụ bù là 1800 kVar và 300 kVar. Ta Giả sử trong 1 ngày xuất tuyến 477AX vận hành thấy tổn hao công suất và chi phí tiết kiệm là tương theo các chế độ phụ tải cao điểm, thấp điểm và bình đồng với trường hợp1. thường như bảng 3.1 Trường hợp 3: Khảo sát vị trí bù khi xét đến số lượng tụ cần bù Theo tiêu chí giảm lượng tổn thất công suất năm 2014 là 8,6% của tổng công ty điện lực việt nam
3. EVN. Đồng thời theo khảo sát bằng phần mềm suất 300 kVar của hãng SAMWHA để thực hiện bù Pss/Adept trong 2 trường hợp nêu trên, điện lực Cà trung thế cho xuất tuyến 477AX. Mau lên dự án mua 15 bộ tụ bù 3 pha cố định có công Ta thực hiện mô phỏng và đạt được kết quả sau: suất 300 kVar và 05 bộ tụ bù 3 pha ứng động có công Bảng 3.2. Bù trung áp tại nhiều vị trí có giới hạn thiết bị bù Tụ cố định Tụ điều chỉnh Chi phí tổn thất Delta P (kW) Delta Q (kVAr) kVAR (nút) kVAR (nút) năng lượng (USD) Ban đầu - - 216.21 427.06 1,651,843.68 21(900); 18(1200); Tối ưu 20(2100); 174.57 336.46 1,313,610.23 44(300) 18(1500) Tiết kiệm 41.64 90.60 338,233.45 Chi phí đặt tụ 56,082.46 cố định Chi phí đặt tụ 45,719.52 điều chỉnh Chi phí lắp 101,801.98 đặt tụ Tổng Tiết 236,431.47 kiệm Kết luận: Qua trường hợp 1 và trường hợp 2 ta Bài toán dự báo phụ tải điện được thiết kế và xây thấy tổng số tiền tiết kiệm được sau khi khấu trừ các dựng dựa vào dữ liệu đầu vào của mạng là Nhiệt độ chi phí lạm phát, trượt giá thì số tiền tiết kiệm được Tmin và nhiệt độ cao Tmax và công suất phụ tải P (Kw) là 324,615.60 $ là lớn nhất. Tuy nhiên trong trường tương ứng cùng ngày, giờ trong vài năm trước. Mạng hợp 2 cũng với số tiền tiết kiệm được tương đương noron trong bài này có 1 lớp ẩn, các mạng có nhiều TH1 nhưng ta đã khảo sát và loại bỏ một số vị trí đặt lớp ẩn sẽ phức tạp hơn và việc huấn luyện mạng tốn tụ không phù hợp với điều kiện vận hành thực tế. nhiều thời gian. Việc chọn số nơron trong lớp ẩn phải Trường hợp 3 với bài toán sử dụng nguồn vốn đầu tư được cân nhắc; nếu quá nhiều thì làm cho mạng quá định mức, từ đó ta nghiên cứu bài toán kinh tế lựa chuyên biệt dẫn đến làm mất đi khả năng phổ dụng chọn các vị trí bù tối ưu và số tiền tiết kiệm được là của nó. Nếu số nơron không đủ thì mạng khó học. Vì 236,431.47 $ vậy các thông số sau được chọn cho đầu vào: 4. Ứng dụng mạng neural cho bài toán dự Có 8 thông số đầu vào và một thông số đầu ra. báo phụ tải điện ngắn hạn.Những yếu tố ảnh - Nhiệt độ Tmax trong một giờ. hưởng đến dự báo phụ tải - Nhiệt độ Tmin trong một giờ. Để tìm ra được các phương pháp tối ưu đối với - Giờ làm việc trong một ngày. bài toán dự báo phụ tải điện với sai số nhỏ thì ta cần - Số ngày làm việc trong một tuần xác định được các yếu tố ảnh hưởng đến phụ tải, ví - Số ngày nghỉ trong tuần. dụ: thứ của ngày trong tuần, ngày lễ, các đặc điểm - Công suất tải trung bình của ngày trước. đặc trưng của ngày cần dự báo, nhiệt độ trong ngày, - Công suất tải trung bình trong một giờ cùng công suất phụ tải trong ngày theo giờ, . Chính vì ngày hôm trước. vậy để nâng cao độ chính xác của quá trình dự báo - Công suất tải trung bình trong cùng một giờ và phụ tải ta cần phải xét đến các yếu tố ảnh hưởng đến cùng một ngày so với tuần trước. phụ tải. - Số neural của lớp ẩn là 20 b. Dữ liệu đầu vào cho bài toán dự báo phụ - Một đầu ra là công suất dự báo P. tải.
4. 01/01/2013, ta có kết quả so sánh trong cùng ngày như sau: Hình 1: Cấu trúc mạng neural Trong quá trình nghiên cứu và xây dựng mạng neural dự báo phụ tải điện các dữ liệu phụ tải P và nhiệt độ Tmax ,Tmin trong ngày của các ngày làm việc từ ngày 01/01/2011 đến 31/12/2012 làm các ví dụ để học và dự báo. Từ số liệu trên ta xây dựng tập Hình 2. Kết quả dự báo so với kết quả thực tế. mẫu dùng để huấn luyện và kiểm tra khả năng dự báo mạng của neural đã qua huấn luyện. Một mẫu cụ thể trong tập hợp mẫu được xây dựng từ tám thông số đầu vào đã nêu ở trên. Kết quả dự báo cho năm 2013 được kiểm tra lại so với kết quả thực tế. Dựa trên hai bảng số liệu của hai năm 2011 và năm 2012, ta nhập vào hàm dự báo để tính giá trị dự báo của năm 2013. Từ đó, thông qua kết quả của hàm dự báo và bảng số liệu công suất thực tế của Hình 3. Sai số dự báo. Tương tự như trên, ta lập bảng kết quả dự báo cho các tổng tất cả các giờ trong năm và lấy giá trị trung bình của tổng tất cả các giờ trong năm Bảng 4.4: Kết quả dự báo cho các năm tiếp theo Năm 2013 Dự báo Dự báo Dự báo Dự báo Giờ Năm 2011 Năm 2012 Thực tế Dự báo Năm 2014 Năm 2015 Năm 2016 Năm 2017 1 7491 7615 7756 7296 7369 7766 9275 9815 2 7143 7304 7460 6995 6998 7435 9151 9650 3 6895 7112 7271 6768 6723 7189 8916 9472 4 6765 7009 7203 6630 6553 7043 8782 9121 5 6758 7032 7300 6636 6546 7064 8812 9385 6 6892 7136 7318 6879 6808 7360 9127 9686 7 7130 7374 7504 7591 7567 8159 9663 9979 8 7279 7503 7690 8418 8691 9284 9920 10239 8754 8023 8124 8287 8350 8710 9118 9580 9918 8755 8111 8278 8365 8557 8891 9234 9625 10149 8756 8168 8296 8484 8669 8976 9336 9802 10365 8757 8202 8412 8582 8790 9007 9403 9978 10553 8758 8163 8313 8491 8532 8843 9257 9645 10283 8759 8015 8201 8414 8380 8659 9015 9348 10094 8760 7818 7972 8097 8187 8449 8931 9212 9982 8761 7566 7764 7911 8045 8230 8630 9088 9825 Trung 7781 7910 8072 8293 8689 9062 9485 9913 bình
5. Từ dữ liệu là tám thông số đầu vào với số liệu công suất ban đầu của phát tuyến 477AX, ta có được kết quả dự báo cho các năm tiếp theo như sau: Tổng công suất Sai số dự Tỷ lệ gia tăng so với công suất tải tăng lên(kw) báo % thực của năm 2013 (kw) Năm 2013 8293 2.73 Năm 2014 8689 - 8689/8072=1.076 Năm 2015 9062 - 9062/8072=1.122 Năm 2016 9485 - 9485/8072=1.175 Năm 2017 9913 - 9913/8072=1.228 Trong bảng 4.4 kết quả dự báo được tính từ ngày 01/01/2013 đến 31/12/2017 công suất dự báo trong năm 2013 với tổng công suất dự báo là 8293 kW và sai số dự báo là 2.73%. 5. Tính dung lượng bù tối ưu cho phát tuyến kiệm được sau khi gắn tụ bù vào lưới điện phân phối khi có dự báo phụ tải điện. phát tuyến 477AX của lưới điện thành phố Cà Mau. Khi công suất của lưới điện tăng lên theo kết quả Giá trị Tỷ lệ gia tăng so với công suất tải thực dự báo của lưới điện hiện hữu thì chúng ta sẽ tính của năm 2013 (kW) dùng để đặt vào giá trị tỷ lệ tải được dung lượng bù công suất phản kháng vào lưới dự báo cho các năm tiếp theo, từ đó ta dùng phần điện, chi phí lắp đặt tụ bù vào lưới điện và chi phí tiết mềm CAPO để tính dung lượng bù tối ưu cho xuất tuyến 477AX được thể hiện ở bảng 4.5. Bảng 4.5: Kết quả tính toán bù trung thế theo dự báo. 2013 2014 2015 2016 2017 Tổng tải (pu) 1 1.076 1.122 1.175 1.228 Cố định Ứng Cố định Ứng Cố định Ứng Cố định Ứng Cố định Ứng (15x300 động (15x300 động (15x300 động (15x300 động (15x300 động kVar) (5x300 kVar) (5x300 kVar) (5x300 kVar) (5x300 kVar) (5x300 kVar) kVar) kVar) kVar) kVar) 21(900) 18(1200) 21(900) 18(1200) 21(900) 18(1200) 21(900) 18(1200) 21(900) 18(1200) Vị trí tụ bù 20(2100) 44(300) 20(2100) 44(300) 20(2100) 44(300) 20(2100) 44(300) 20(2100) 44(300) 18(1500) 18(1500) 18(1500) 18(1500) 18(1500) 12(300) 12(300) 12(300) 12(300) 7(300) 7(300) 7(300) 6(300) 6(300) 10(300) Giá tiền tiết 236,431 5,244 3,029 1,440 86 kiệm ($) ($) ($) ($) ($) 6. Kết Luận Tiếng Việt - Giai đoạn 1: Thực hiện bù trong năm 2013 [1]. Lê Thị Thanh Hà (2000), Trường Đại Học cho cấu hình hiện hữu với 15x300 kVar dung lượng Bách Khoa Hà Nội, “Luận văn tốt nghiệp cao học”. cố định tại các nút 21, 20, 18 và 5x300 kVar tại các [2]. LM Markovits, Trần Đình Long, Bùi Ngọc nút 18 và 44. Thư, Bùi Thiên Dụ, Hà Học Trạc (dịch), Các chế độ - Giai đoạn 2: Thực hiện bù thêm trong năm của hệ thống năng lượng, NXB KH-KT 2014 với tụ cố định 1x300 kVar tại nút số 12 tương [3]. Sách “Truyền tải và Phân phối Hệ Thống ứng với sự phát triển của phụ tải trong năm 2014. Điện” của tác giả Hồ Văn Hiến, Nhà Xuất Bản Đại - Giai đoạn 3: Thực hiện bù thêm trong năm Học Quốc Gia Tp. Hồ Chí Minh. 2015 với tụ cố định 1x300 kVar tại nút số 7 tương [4] “Áp dụng PSS/ADEPT 5.0 trong Lưới điện ứng với sự phát triển của phụ tải trong năm 2015. phân phối” của PGS.TS. Nguyễn Hữu Phúc của Đai Học Điện Lực. TÀI LIỆU THAM KHẢO
6. [5]. Đinh Thành Việt - Trường Đại học Bách Wu, Senior Member, IEEE, “Electric Load khoa, Đại học Đà Nẵng, “Nghiên cứu lựa chọn Forecasting Based on Locally Weighted Support phương thức vận hành tối ưu cho lưới điện phân phối Vector Regression” IEEE transactions on systems, tỉnh phú yên.” man, and cybernetics—part c: applications and [6]. Trần Kỳ Phúc (12/2006), “Nghiên cứu ứng reviews, vol. 40, no. 4, july 2010. dụng mạng Nơ ron nhân tạo trong dự báo ngắn hạn [15]. A.J. Al-Shareef and M.F. Abbod phụ tải điện (giai đoạn I)”, Viện năng lượng. “Intelligent Modelling Techniques of Power Load [7]. Nguyễn Hoàng Việt1 , Trần Anh Dũng1 , Forecasting for the Western Area of Saudi Arabia”, Nguyễn Quang Thi2, 1Trường Đại học Bách Khoa, JKAU: Eng. Sci., Vol. 21 No.1 pp: 3-18 ĐHQG-HCM, 2Công ty Điện lực Tp.HCM, “Mạng (1431A.H./2010 A.D.) DOI: 10.4197 / Eng. 21-1.1. wavelet cho bài toán dự báo phụ tải ngắn hạn trong [16]. Y. Rui A.A. El-Keib, “A Review of ANN- các ngày đặc biệt”, Tạp chí phát triển kh&cn, tập 10, based Short-Term Load Forecasting Models”, số 06 – 2007”. Department of Electrical Engineering University of [8]. Nguyễn Sĩ Dũng 1 , Lê Hòai Quốc 2, Trường Alabama, Tuscaloosa, AL 35487. Đại học Công nghiệp thành phố Hồ Chí Min, 2 trường [17]. Nima Amjady and Farshid Keynia, “A New Đại học Bách khoa thành phố Hồ Chí Minh, “Một Neural Network Approach to Short Term Load thuật toán về huấn luyện mạng neural networks trên Forecasting of Electrical Power Systems”, cơ sở phương pháp conjugate gradient”, Tạp chí phát Department of Electrical Engineering, Semnan triển kh&cn, tập 10, số 06 – 2007”. University, Semnan, Iran. [9]. Trần Kỳ Phúc , “Nghiên cứu ứng dụng mạng [18]. Mark Hudson Beale, Martin T. Hagan, nơ ron nhân tạo trong dự báo ngắn hạn phụ tải điện Howard B. Demuth, “Neural Network Toolbox™ 7 (giai đoạn 2)”, Báo cáo tổng kết đề tài nghiên cứu User’s Guide”. khoa học cấp bộ mã số: i- 145, hà nội, 12-2008. [19]. Farid Naimi, “Java Based Load [10].“Bù tối ưu công suất phản kháng lưới điện Forecasting for a Small Number of Public phân phối”, Tạp chí khoa học và công nghệ, Đại học Properties” Day of delivery March 20th 2008. Đà nẵng - số 2(25).2008 [20]. Jayant D. Sawarkar 1, Umesh L. Kulkarni Tiếng Anh 2, Dr. Sudhir kumar Sawarkar, “Prediction of Short [11]. Howard Demuth Mark eale, “Neural Term Electric Load Using Artificial Neural Network Network Toolbox For Use with MATLAB®”. ”, Available Online at www.ijecse.org ISSN- 2277- [11]. Hesham K. Alfares* and Mohammad 1956, International Journal of Electronics and Nazeeruddin “Electric load forecasting: literature Computer Science Engineering 992. survey and classifition cation of methods”, [21]. Ehab E. Elattar, Student Member, IEEE, John International Journal of Systems Science, 2002, (Yannis) Goulermas, Member, IEEE, and Q. H. Wu, volume 33, number 1, pages 23 ±34. Senior Member, IEEE, “Electric Load Forecasting [13]. Saahil Shenoy and Dimitry Gorinevsky Based on Locally Weighted Support Vector “Risk Adjusted Forecasting of Electric Power Loa”, Regression” IEEE Transactions on systems, man, 2014 American Control Conference (ACC) June 4-6, and cybernetics—part c: applications and reviews, 2014. Portland, Oregon, USA. vol. 40, no. 4, july 2010 [14]. Ehab E. Elattar, Student Member, IEEE, John (Yannis) Goulermas, Member, IEEE, and Q. H. GV Hướng Dẫn: Ts. Lê Kỹ HVTH: Đinh Công Thành Ký Tên và Ghi rõ họ tên Ký Tên và Ghi rõ họ tên
7. BÀI BÁO KHOA HỌC THỰC HIỆN CÔNG BỐ THEO QUY CHẾ ĐÀO TẠO THẠC SỸ Bài báo khoa học của học viên có xác nhận và đề xuất cho đăng của Giảng viên hướng dẫn Bản tiếng Việt ©, TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT TP. HỒ CHÍ MINH và TÁC GIẢ Bản quyền tác phẩm đã được bảo hộ bởi Luật xuất bản và Luật Sở hữu trí tuệ Việt Nam. Nghiêm cấm mọi hình thức xuất bản, sao chụp, phát tán nội dung khi chưa có sự đồng ý của tác giả và Trường Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh. ĐỂ CÓ BÀI BÁO KHOA HỌC TỐT, CẦN CHUNG TAY BẢO VỆ TÁC QUYỀN! Thực hiện theo MTCL & KHTHMTCL Năm học 2016-2017 của Thư viện Trường Đại học Sư phạm Kỹ thuật Tp. Hồ Chí Minh.