Nghiên cứu đáp ứng quá độ của hệ thống nối đất đối với dòng điện sét trên đường dây truyền tải việt na

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI - NGUYỄN XUÂN PHÚC NGHIÊN CỨU ĐÁP ỨNG QUÁ ĐỘ CỦA HỆ THỐNG NỐI ĐẤT ĐỐI VỚI DÒNG ĐIỆN SÉT TRÊN ĐƯỜNG DÂY TRUYỀN TẢI VIỆT NAM LUẬN ÁN TIẾN SĨ KỸ THUẬT ĐIỆN Hà Nội – 2020 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI - NGUYỄN XUÂN PHÚC NGHIÊN CỨU ĐÁP ỨNG QUÁ ĐỘ CỦA HỆ THỐNG NỐI ĐẤT ĐỐI VỚI DÒNG ĐIỆN SÉT TRÊN ĐƯỜNG DÂY TRUYỀN TẢI VIỆT NAM Ngành: Kỹ thuật điện Mã số: 9520201 LUẬN ÁN TIẾN SĨ KỸ THUẬT ĐIỆN NGƯỜI HƯỚNG DẪN KHOA HỌC: TS PHẠM HỒNG THỊNH PGS.TS TRẦN VĂN TỚP Hà Nội – 2020 LỜI CAM ĐOAN Tôi xin cam đoan cơng trình nghiên cứu thân dựa hướng dẫn tập thể hướng dẫn khoa học tài liệu tham khảo trích dẫn Các kết đạt luận án xác, trung thực chưa cơng bố cơng trình khác Hà Nội, ngày … tháng … năm 2020 XÁC NHẬN CỦA TẬP THỂ HƯỚNG DẪN GV HƯỚNG DẪN GV HƯỚNG DẪN TÁC GIẢ LUẬN ÁN TS Phạm Hồng Thịnh PGS TS Trần Văn Tớp Nguyễn Xuân Phúc i LỜI CẢM ƠN Trước hết, tơi xin bày tỏ lịng biết ơn sâu sắc đến tập thể hướng dẫn khoa học, TS Phạm Hồng Thịnh PGS TS Trần Văn Tớp ln tận tình giúp đỡ, hướng dẫn, bảo tơi suốt q trình nghiên cứu Hai thầy dành nhiều thời gian tâm huyết, hỗ trợ mặt để tơi hồn thành luận án Tơi xin trân trọng cảm ơn Lãnh đạo trường Đại học Bách Khoa Hà Nội, Phòng đào tạo – Bộ phận đào tạo Sau đại học, Viện điện Bộ môn Hệ thống điện tạo điều kiện thuận lợi cho nghiên cứu sinh suốt trình học tập nghiên cứu Tôi xin chân thành cảm ơn Giảng viên cán Bộ môn Hệ thống điện trường Đại học Bách Khoa Hà Nội tận tình hỗ trợ giúp đỡ trình thực luận án Tơi xin trân trọng cảm ơn Lãnh đạo Viện Năng lượng, Trung tâm Tư vấn Năng lượng Chuyển giao công nghệ tạo điều kiện thuận lợi cho suốt thời gian qua Xin chân thành cảm ơn quan tâm, giúp đỡ động viên đồng nghiệp Viện Năng lượng Cuối cùng, xin chân thành cảm ơn bạn bè gia đình ln bên cạnh động viên, giúp đỡ tơi hồn thành luận án Hà Nội, ngày … tháng … năm 2020 Tác giả luận án Nguyễn Xuân Phúc ii MỤC LỤC LỜI CAM ĐOAN i LỜI CẢM ƠN ii MỤC LỤC iii DANH MỤC TỪ VIẾT TẮT vi DANH MỤC HÌNH VẼ viii DANH MỤC BẢNG BIỂU xiii MỞ ĐẦU 1 Tính cấp thiết đề tài Mục đích nghiên cứu luận án .2 Đối tượng phạm vi nghiên cứu .2 Phương pháp nghiên cứu Ý nghĩa khoa học thực tiễn đề tài Những đóng góp luận án Cấu trúc luận án Chương TỔNG QUAN 1.1 TỔNG QUAN VỀ LƯỚI ĐIỆN TRUYỀN TẢI VIỆT NAM 1.1.1 Tình hình phát triển lưới điện truyền tải .6 1.1.2 Tình hình cố lưới điện truyền tải 1.1.3 Các hệ thống nối đất đường dây truyền tải Việt Nam 1.2 NỐI ĐẤT CHO ĐƯỜNG DÂY TRUYỀN TẢI TRÊN THẾ GIỚI 13 1.2.1 Quy phạm cho hệ thống nối đất 13 1.2.2 Hình dạng nối đất .15 1.3 NGHIÊN CỨU ĐÁP ỨNG QUÁ ĐỘ CỦA HỆ THỐNG NỐI ĐẤT CỦA ĐƯỜNG DÂY TRUYỀN TẢI 15 1.3.1 Nghiên cứu nước .15 1.3.2 Nghiên cứu giới 16 1.4 KẾT LUẬN 21 Chương 23 TÍNH TỐN NỐI ĐẤT VÀ PHƯƠNG PHÁP PHẦN TỬ HỮU HẠN TRONG TÍNH TỐN NỐI ĐẤT 23 iii 2.1 PHƯƠNG PHÁP PHẦN TỬ HỮU HẠN (FEM) TRONG TÍNH TỐN HỆ THỐNG NỐI ĐẤT 23 2.1.1 Giới thiệu chung 23 2.1.2 Các bước giải toán sử dụng phương pháp FEM 24 2.2 TÍNH TỐN KIỂM CHỨNG GIỮA PHƯƠNG PHÁP FEM VỚI THỰC NGHIỆM 26 2.2.1 Tính toán điện trở chiều cọc đất phân tầng 26 2.2.2 Tính tốn điện áp độ tia nối đất 28 2.2.3 Tính tốn điện áp q độ cột điện 31 2.2.4 Tính tốn điện áp bước điện áp tiếp xúc 32 2.3 KẾT LUẬN 34 Chương 35 ĐÁP ỨNG CỦA HỆ THỐNG NỐI ĐẤT TRONG MIỀN TẦN SỐ 35 3.1 ĐẶT VẤN ĐỀ 35 3.2 THÔNG SỐ MÔ PHỎNG 37 3.3 KẾT QUẢ MÔ PHỎNG 38 3.3.1 Đáp ứng miền tần số hệ thống dạng cọc 38 3.3.2 Đáp ứng hệ thống dạng tia 41 3.3.3 Đáp ứng hệ thống dạng cọc - tia 44 3.3.4 Hệ thống dạng tia quấn vòng 47 3.3.5 Ảnh hưởng hình dạng hệ thống nối đất .48 3.3.6 Ảnh hưởng điện trở suất đất .51 3.3.7 Ảnh hưởng số điện môi đất 52 3.3.8 Tính tốn cho đất phân tầng 53 3.4 KẾT LUẬN 57 Chương 59 ĐÁP ỨNG QUÁ ĐỘ CỦA HỆ THỐNG NỐI ĐẤT TRONG MIỀN THỜI GIAN 59 4.1 ĐẶT VẤN ĐỀ 59 4.2 MƠ HÌNH MƠ PHỎNG 61 4.3 ĐÁP ỨNG QUÁ ĐỘ CỦA CÁC HỆ THỐNG NỐI ĐẤT 61 4.3.1 Hệ thống tia .61 4.3.2 Ảnh hưởng chiều dài tia .63 4.3.3 Hệ thống nhiều tia .66 iv 4.3.4 Ảnh hưởng số lượng tia 68 4.3.5 Hệ thống dạng cọc – tia 71 4.3.6 Ảnh hưởng số lượng cọc 72 4.3.7 Ảnh hưởng hình dạng hệ thống nối đất .74 4.4 HIỆN TƯỢNG PHÓNG ĐIỆN TRONG ĐẤT 75 4.4.1 Các giá trị đặc trưng cho tượng phóng điện đất 75 4.4.2 Sơ đồ tính toán 77 4.4.3 Kết nghiên cứu 78 4.5 KẾT LUẬN 89 Chương 91 ĐÁP ỨNG CỦA HỆ THỐNG NỐI ĐẤT TRONG HỆ ĐẦY ĐỦ CỘT VÀ NỐI ĐẤT 91 5.1 ĐẶT VẤN ĐỀ 91 5.2 TÍNH TỐN TỔNG TRỞ SĨNG CỦA CỘT KHI BỎ QUA HỆ THỐNG NỐI ĐẤT 92 5.2.1 Lý thuyết tổng trở sóng cột 92 5.2.2 Tính tốn tổng trở sóng cột phương pháp giải tích 93 5.2.3 Tính tốn tổng trở sóng cột theo phương pháp số 94 5.2.4 Đánh giá phương pháp mô .95 5.2.5 Tổng trở sóng cột điển hình đường dây truyền tải 96 5.3 ĐIỆN ÁP ĐỈNH CỘT TRONG HỆ CỘT – HỆ THỐNG NỐI ĐẤT 96 5.3.1 Mơ hình nghiên cứu 96 5.3.2 Hệ cột điện 500 kV tia nối đất 98 5.3.3 Ảnh hưởng số lượng tia đến điện áp đỉnh cột 99 5.3.4 Ảnh hưởng hình dạng nối đất 101 5.3.5 Ảnh hưởng kích thước cột 103 5.3.6 Ảnh hưởng hình dạng xung 105 5.3.7 So sánh với mơ hình nối đất tập trung 107 5.4 KẾT LUẬN 109 KẾT LUẬN VÀ KIẾN NGHỊ 111 KẾT LUẬN 111 KIẾN NGHỊ 112 DANH MỤC CÁC CƠNG TRÌNH CÔNG BỐ CỦA LUẬN ÁN 113 TÀI LIỆU THAM KHẢO 114 PHỤ LỤC 123 v DANH MỤC TỪ VIẾT TẮT Từ viết tắt Nguyên nghĩa HTĐ Hệ thống điện EVN Tập đoàn điện lực Việt Nam EVNNPT Tổng công ty truyền tải điện Quốc gia FEM Phần tử hữu hạn MoM Phương pháp Mô men FDTD Sai phân hữu hạn miền thời gian GPR Điện dâng đất (Ground Potential Rise) (V) Z Tổng trở () 10 ρ Điện trở suất đất (.m) 11 ε0 Hằng số điện môi chân không (ε0 = 8,854.10-12 F/m) 12 μ0 Độ từ thẩm chân không (μ = 1,256.10-6 H/m) 13 εr Hằng số điện môi tương đối 14 μr Độ từ thẩm tương đối 15 re Điện trở đơn vị (/m) 16 L Điện cảm đơn vị (H/m) 17 G Điện dẫn đơn vị (S/m) 18 C Điện dung đơn vị (F/m) 19 J Mật độ dòng điện (A/m2) 20 E Cường độ điện trường (V/m) 21 D Mật độ điện cảm (C/m2) 22 H Cường độ từ trường (A/m) 23 B Cảm ứng từ (T/m2) 24 A Vec tơ từ (A) STT vi Từ viết tắt Nguyên nghĩa 25 V Điện (V) 26 i(t) Xung dòng biến đổi theo thời gian (A) 27 v(t) Điện áp độ biến đổi theo thời gian (V) 28 Zp Tổng trở xung (impulse impedance) hệ thống nối đất (Ω) 29 Vp Điện cực đại (V) 30 Ip Dòng điện cực đại (A) 31  Thời gian đầu sóng (μs) 32 ZT Tổng trở sóng cột điện (Ω) 33 Zg Tổng trở sóng dây chống sét (Ω) 34 αT Hệ số phản xạ đỉnh cột 35 T Thời gian truyền sóng cột 36 Zw Tổng trở dạng sóng cột (tower wave impedance) (Ω) 37  Hệ số suy giảm (damping factor) (Ω) 38 c Vận tốc ánh sáng (299.792 km/s) STT vii DANH MỤC HÌNH VẼ Hình 1.1 Chiều dài đường dây 220 kV 500 kV giai đoạn 2011 – 2019 [2, 3] Hình 1.2 Chiều dài đường dây 220 kV 500 kV giai đoạn 2020 – 2030 [4] Hình 1.3 Thống kê cố lưới truyền tải Việt Nam [2, 3] Hình 1.4 Nối đất dạng tia thẳng Hình 1.5 Cách bố trí tia thơng thường Hình 1.6 Cách bố trí tia đặc biệt Hình1.7 Cấu trúc nối đất kiểu cọc chơn sâu 10 Hình 1.8 Nối đất dạng cọc – tia 10 Hình 1.9 Thống kê khoảng cách cọc hệ thống nối đất 123 vị trí cột đường dây 220 kV Hịa Bình – Tây Hà Nội [6] 11 Hình 1.10 Cấu trúc điển hình nối đất kiểu quấn vòng 12 Hình 1.11 Thống kê hình dạng hệ thống nối đất sử dụng lưới điện truyền tải PTC1 quản lý [7] 13 Hình 1.12 Các cách bố trí tia nối đất [22] 15 Hình 1.13 Nối đất dạng vịng kín [13] 15 Hình 1.14 Sơ đồ mạch thay tia nối đất [29] 17 Hình 1.15 Mơ hình thay phương pháp tính phân bố dòng áp dọc theo chiều dài tia nối đất Lorentzou et.al [43] 17 Hình 1.12 Các phần từ lưới [ 83] 25 Hình 2.1 Cách thiết lập mơ hình tính tốn điện trở nối đất cọc nối đất lớp đất hai tầng so sánh đối chứng với thí nghiệm Ahmeda 26 Hình 2.2 Kết mơ cho cọc nối đất thí nghiệm [88] 27 Hình 2.3 Mơ hình mơ cho tia nối đất thí nghiệm Harid [89] 28 Hình 2.4 Biên ”lumped port” modul RF (Lumped port boundary) [86] 29 Hình 2.5 Thiết lập biên “Lumped port” mơ hình hệ thống nối đất 29 Hình 2.6 Dòng điện 4,9/13 s biên độ 2,4 A vào điện cực mơ hình COMSOL 30 Hình 2.7 Phân bố điện trường mặt phẳng song song với mặt đất qua tia 30 Hình 2.8 So sánh kết tính tốn điện điện cực phương pháp FEM (đường màu đỏ) với thí nghiệm Harid [89] (đường màu xanh) 30 Hình 2.9 Dịng xung ramp 1/70 s biên độ 1A vào đỉnh cột 31 Hình 2.10 Mơ hình mơ thí nghiệm Yamada [90] 31 viii [101] L Grcev (2008), Impulse efficiency of ground electrodes, IEEE Transactions on Power Delivery, Vol 24, No.1, pp 441-451 [102] J Wang, A C Liew, M Darveniza (2005), Extension of dynamic model of impulse behavior of concentrated grounds at high currents IEEE Transactions on Power Delivery, Vol 20, No 3, pp 2160-2165 [103] E E Oettle (1988), A new general estimation curve for predicting the impulse impedance of concentrated earth electrodes, IEEE Transactions on Power Delivery, Vol 3, No 4, pp 2020-2029 [104] Võ Viết Đạn (1972), Giáo Trình Kỹ Thuật Điện Cao Áp, Khoa đại học chức trường đại học Bách Khoa Hà Nội [105] Nguyễn Minh Chước (2002), Hướng dẫn thiết kế đồ án tốt nghiệp kỹ thuật điện cao áp, Nhà xuất Hà Nội [106] CIGRE (1991), CIGRE guide to procedure for estimating the lightning performances of transmission line, CIGRE SC33-WG01 [107] B R Gupta (1981), Impulse characteristics of grounding electrodes, Inst Eng (India), vol 64, no 4, pp 178-182 [108] Jinliang He et al (2005), Effective Length of counterpoise wire under lightning current, IEEE Transaction on Power Delivery, Vol 20, No 2, pp.1585-1591 [109] E E Oettle (1988), A new general estimation curve for predicting the impulse impedance of concentrated earth electrodes, IEEE Transactions on Power Delivery, Vol 3, No 4, pp 2020-2029 [110] M I Lorentzou, N D Hatziargyriou, B C Papadias, S Sekioka (2001), Modeling of Soil Ionization for Calculation of the Transient Response of Grounding Electrodes, 2001 International Conference on Power Systems Transients [111] Y Liu, N Theethayi, R Thottappillil, R M Gonzalez, M Zitnik (2004), An improved model for soil ionization around grounding system and its application to stratified soil, Journal of Electrostatics, Vol 60, No 2-4, pp.203-209 [112] A C Liew, M Darveniza (1974), Dynamic model of impulse characteristics of concentrated earths, Proceedings of the Institution of electrical Engineers, Vol 121, No 2, pp 123-135 [113] Y Gao, J He, J Zou, R Zeng, X Liang (2004), Fractal simulation of soil breakdown under lightning current, Journal of Electrostatics, Vol 61, No 3-4, pp 197-207 [114] B R Gupta, B Thapar (1980), Impulse impedance of grounding grids, IEEE Transactions on Power Apparatus and Systems, No 6, pp 2357-2362 [115] Y Baba, M Ishii (1999), Numerical electromagnetic field analysis on measuring methods of tower surge impedance, IEEE Transactions on Power Delivery, Vol 14, No 2, pp 630-635 121 [116] W A Chisholm, Y L Chow, K D Srivastava (1983), Lightning surge response of transmission towers, IEEE Transactions on Power Apparatus and Systems, No 9, pp 3232-3242 [117] J J LaForest (1981), Transmission-line reference book 345 kV and above, p.570 (No EPRI-EL-2500), General Electric Co., Pittsfield, MA (USA), Large Transformer Division General Electric Co., Schenectady, NY (USA) Electric Utility Systems Engineering Dept [118] C F Wagener, A R Hileman (1960), A new approach to the calculation of the lighting performance of transmission lines III, Electrical Engineering, Vol 78, No 4, pp 996-1020 [119] M A Sargent, M Darveniza (1969), Tower surge impedance IEEE Transactions on Power Apparatus and Systems, No.5, pp 680-687 [120] C A Jordan (1934), Lightning computations for transmission lines with overhead ground wires, Part II General Electric Review, No 4, pp 180-186 [121] IEEE Working Group Report (1993), Estimating Lighnting Performance of Transmission Lines II- Updates to Analytical Models [122] M Ishii, T Kawamura, T Kouno, E Ohsaki, K Shiokawa, K Murotani, and T Higuchi (1991), Multistory transmission tower model for lightning surge analysis, IEEE Transactions on Power Delivery, Vol 6, No 3, pp 1327–1335 [123] T Ito, T Ueda, H Watanabe, T Funabashi, and A Ametani (2003), Lightning flashover on 77-kV systems: Observed voltage bias effects and analysis, IEEE Transactions on Power Delivery, Vol 18, No 2, pp 545–550 [124] P C A Mota, M L R Chaves, J R Camacho (2012), Power line tower lightning surge impedance computation, a comparison of analytical and finite element methods, Internaltional Conference on Renewable Energies and Power Quality, Spain [125] T X Cao, T Pham and S Boggs (2013), Computation of tower surge impedance in transmission line, 2013 IEEE Electrical Insulation Conference (EIC), Ottawa, ON, pp 7780 [126] F P Dawalibi, W Ruan, S Fortin, J Ma, W K Daily (2001), Computation of power line structure surge impedances using the electromagnetic field method, IEEE/PES Transmission and Distribution Conference and Exposition Developing New Perspectives, IEEE, Vol 2, pp 663-668 [127] Motoyama, Hideki (1996), Experimental study and analysis of breakdown characteristics of long air gaps with short tail lightning impulse, IEEE Transactions on Power Delivery, Vol 11, No 2, pp 972-979 [128] Ninh Văn Nam (2020), Nghiên cứu số giải pháp giảm cố sét cho đường dây truyền tải điện không, Luận án tiến sĩ, Đại học Bách Khoa Hà Nội 122 PHỤ LỤC 123 PHỤ LỤC CHƯƠNG Chương mô đáp ứng độ miền tần số với dải tần từ Hz đến MHz cho điện cực cọc, tia, cọc – tia tia quấn vòng Các kịch mô sau: - Kịch 1: mô cho cọc nối đất (PL 3.1) PL 3.1 Mô hình COMSOL cọc nối đất có chiều dài 10 m, 20 m 35 m, đặt đất có điện trở suất ρ = 1000 .m số điện môi tương quan r = Riêng trường hợp cọc dài 35 m, có thêm trường hợp điện trở suất 2000 .m 5000 .m (r = cho trường hợp); số điện môi 20, 40, 60, 80 (ρ = 1000 .m cho trường hợp) Kịch 2: mô cho tia nối đất (PL 3.2) PL 3.2 Mơ hình COMSOL tia nối đất có chiều dài 10 m, 20 m, 40 m 50 m, đặt độ sâu 0,8 m đất có điện trở suất ρ = 1000 .m số điện môi tương quan r = 124 Trường hợp tia, chiều dài tia 50 m khoảng cách tia m, m 10 m - Kịch 3: mô cho hệ thống cọc - tia (PL 3.3) PL 3.3 Mơ hình COMSOL hệ thống cọc - tia, đặt độ sâu 0,8 m đất có điện trở suất ρ = 1000 .m số điện môi tương quan r = Với chung giá trị điện trở chiều 33 Ω, có cấu hình cọc tia: + Tia dài 36 m kết hợp với 12 cọc, cọc dài m, khoảng cách cọc m + Tia dài 40 m kết hợp cọc, cọc dài m, khoảng cách cọc m + Tia dài 44 m kết hợp cọc, cọc dài m, khoảng cách cọc 20 m - Kịch 4: mô cho hệ thống quấn vịng (PL 3.4) PL 3.4 Mơ hình COMSOL hệ thống tia quấn vòng, đặt độ sâu 0,8 m đất có điện trở suất ρ = 1000 .m số điện môi tương quan r = Trong trường hợp quấn vòng, vòng nằm mặt phẳng với khoảng cách vòng m cạnh hình chữ nhật lớn 10 m x 12 m Trường hợp vịng quấn đầu tia, 125 vịng quấn hình chữ nhật với chiều dài cạnh m x 12 m, đoạn tia thẳng 15 m PHỤ LỤC CHƯƠNG Chương mô đáp ứng độ miền thời gian với dịng xung có hình dạng 1,2/50 s, biên độ kA không xét đến ion hóa biên độ 30 kA xét đến ion hóa cho điện cực tia, cọc – tia, cọc tia quấn vịng Các kịch mơ sau: - Kịch 1: mô cho tia nối đất PL 4.1 Mơ hình COMSOL tia nối đất đặt độ sâu 0,8 m có chiều dài 20 m, 30 m, 40 m, 50 m đặt đất có điện trở suất ρ = 1000 .m số điện môi tương quan r = - Kịch 2: mô cho tia song song (PL 4.2) 126 PL4.2 Mơ hình COMSOL tia nối đất có chiều dài 50 m đặt cách m, 10 m 20 m độ sâu 0,8 m đất có điện trở suất ρ = 1000 .m số điện môi tương quan r = Trường hợp tia dài 20 m, khoảng cách tia m - Kịch 3: mô cho tia song song (PL 4.3) PL4.3 Mơ hình COMSOL tia nối đất có tổng chiều dài 40 m đặt cách m, độ sâu 0,8 m đất có điện trở suất ρ = 1000 .m số điện môi tương quan r = - Kịch 4: mô cho cọc – tia độ sâu 0,8 m đất có điện trở suất ρ = 1000 .m số điện môi tương quan r = 5, mô hình tương tự PL 3.3 - Kịch 5: mô cho cọc dài 35 m đặt đất có điện trở suất ρ = 1000 .m số điện mơi tương quan r = 5, mơ hình tương tự PL 3.1 - Kịch 6: mô cho tia quấn vòng vòng, đặt đất có điện trở suất ρ = 1000 .m số điện môi tương quan r = 5, mơ hình tương tự PL 3.4 - Kịch 7: mơ hình phóng điện đất Mơ hình COMSOL kết hợp Matlab với script sau: clear; clc; model=mphload('ion'); fid=fopen(filename,'wt'); model.hist.disable; ro = 1000; Ec = 241000*(ro)^0.215 127 Imax = 30; r1 = 0.01; rion1 = 0.1; rion2 = 0.1; rion3 = 0.1; rion4 = 0.1; rion5 = 0.1; ro1 = ro; ro2 = ro; ro3 = ro; ro4 = ro; ro5 = ro; l=10; model.param.set('ro',ro); model.param.set('ro1',ro1); model.param.set('ro2',ro2); model.param.set('ro3',ro3); model.param.set('ro4',ro4); model.param.set('ro5',ro5); model.param.set('rion1',rion1); model.param.set('rion2',rion2); model.param.set('rion3',rion3); model.param.set('rion4',rion4); model.param.set('rion5',rion5); model.param.set('Imax',Imax); for tm = [0:0.1:5] 128 model.param.set('tm',tm); model.sol('sol1').run; I11 = mphint2(model,'(temw.Jx*nx+temw.Jy*ny+temw.Jz*nz)','surface','selection',[ 694,695,699,702,708,709,713,716,723,724,728,731,738,739,743,746,753,754,758,761,768,76 9,773,776,779,780,784,787,798,801,808,809,813,816,823,824,828,831,838,839,843,846,853, 854,858,861]) I12 = mphint2(model,'(temw.Jx*nx+temw.Jy*ny+temw.Jz*nz)','surface','selection',[ 524,525,529,532,538,539,543,546,553,554,558,561,568,569,573,576,583,584,588,591,598,59 9,603,606,609,610,614,617,623,624,628,631,638,639,643,646,653,654,658,661,668,669,673, 676,683,684,688,691]) I13 = mphint2(model,'(temw.Jx*nx+temw.Jy*ny+temw.Jz*nz)','surface','selection',[ 354,355,359,362,368,369,373,376,383,384,388,391,398,399,403,406,413,414,418,421,428,42 9,433,436,439,440,444,447,453,454,458,461,468,469,473,476,483,484,488,491,498,499,503, 506,513,514,518,521]) I14 = mphint2(model,'(temw.Jx*nx+temw.Jy*ny+temw.Jz*nz)','surface','selection',[ 184,185,189,192,198,199,203,206,213,214,218,221,228,229,233,236,243,244,248,251,258,25 9,263,266,269,270,274,277,283,284,288,291,298,299,303,306,313,314,318,321,328,329,333, 336,343,344,348,351]) I15 = mphint2(model,'(temw.Jx*nx+temw.Jy*ny+temw.Jz*nz)','surface','selection',[ 10,11,16,21,27,29,33,36,43,44,48,51,58,59,63,66,73,74,78,81,88,89,93,96,99,100,104,107,11 3,114,118,121,128,129,133,136,143,144,148,151,158,159,163,166,173,174,178,181]) U0 = mphglobal(model,'temw.Vport_1'); I0 = mphglobal(model,'temw.Iport_1'); Z0 = U0(2)/I0(2) rion1 = (ro)*I11(2)/2/3.14/l1/Ec rion2 = (ro)*I12(2)/2/3.14/l1/Ec rion3 = (ro)*I13(2)/2/3.14/l1/Ec rion4 = (ro)*I14(2)/2/3.14/l1/Ec rion5 = (ro)*I15(2)/2/3.14/l1/Ec 129 rss = max(rion1, rion2); rss = max(rss, rion3); rss = max(rss, rion4); rss = max(rss, rion5); if rss r1 ro1 = 0.07*ro; model.param.set('rion1',rion1); model.param.set('ro1',ro1); else ro1 = ro; rion1 = 0.1; model.param.set('rion1',rion1); model.param.set('ro1',ro1); end if rion2 > r1 ro2 = 0.07*ro; model.param.set('rion2',rion2); model.param.set('ro2',ro2); 130 else ro2 = ro; rion2 = 0.1; model.param.set('rion2',rion2); model.param.set('ro2',ro2); end if rion3 > r1 ro3 = 0.07*ro; model.param.set('rion3',rion3); model.param.set('ro3',ro3); else ro3 = ro; rion3 = 0.1; model.param.set('rion3',rion3); model.param.set('ro3',ro3); end if rion4 > r1 ro4 = 0.07*ro; model.param.set('rion4',rion4); model.param.set('ro4',ro4); else ro4 = ro; rion4 = 0.1; model.param.set('rion4',rion4); model.param.set('ro4',ro4); end if rion5 > r1 ro5 = 0.07*ro; model.param.set('rion5',rion5); model.param.set('ro5',ro5); 131 else ro5 = ro; rion5 = 0.1; model.param.set('rion5',rion5); model.param.set('ro5',ro5); end model.sol('sol1').run; I21 = mphint2(model,'(temw.Jx*nx+temw.Jy*ny+temw.Jz*nz)','surface','selection',[ 694,695,699,702,708,709,713,716,723,724,728,731,738,739,743,746,753,754,758,761,768,76 9,773,776,779,780,784,787,798,801,808,809,813,816,823,824,828,831,838,839,843,846,853, 854,858,861]) I22 = mphint2(model,'(temw.Jx*nx+temw.Jy*ny+temw.Jz*nz)','surface','selection',[ 524,525,529,532,538,539,543,546,553,554,558,561,568,569,573,576,583,584,588,591,598,59 9,603,606,609,610,614,617,623,624,628,631,638,639,643,646,653,654,658,661,668,669,673, 676,683,684,688,691]) I23 = mphint2(model,'(temw.Jx*nx+temw.Jy*ny+temw.Jz*nz)','surface','selection',[ 354,355,359,362,368,369,373,376,383,384,388,391,398,399,403,406,413,414,418,421,428,42 9,433,436,439,440,444,447,453,454,458,461,468,469,473,476,483,484,488,491,498,499,503, 506,513,514,518,521]) I24 = mphint2(model,'(temw.Jx*nx+temw.Jy*ny+temw.Jz*nz)','surface','selection',[ 184,185,189,192,198,199,203,206,213,214,218,221,228,229,233,236,243,244,248,251,258,25 9,263,266,269,270,274,277,283,284,288,291,298,299,303,306,313,314,318,321,328,329,333, 336,343,344,348,351]) I25 = mphint2(model,'(temw.Jx*nx+temw.Jy*ny+temw.Jz*nz)','surface','selection',[ 10,11,16,21,27,29,33,36,43,44,48,51,58,59,63,66,73,74,78,81,88,89,93,96,99,100,104,107,11 3,114,118,121,128,129,133,136,143,144,148,151,158,159,163,166,173,174,178,181]) U02 = mphglobal(model,'temw.Vport_1') I02 = mphglobal(model,'temw.Iport_1') 132 Z02 = U02(2)/I02(2) end end 133 PHỤ LỤC CHƯƠNG PL 5.1 Sơ đồ cột 500 kV, Tỷ lệ 1/300 134 PL 5.2 Sơ đồ cột 220 kV, Tỷ lệ 1/200 135 ... truyền tải Việt Nam  Nghiên cứu đáp ứng độ miền tần số hệ thống nối đất điển hình cho đường dây truyền tải  Nghiên cứu đáp ứng độ miền miền thời gian hệ thống nối đất điển hình cho đường dây truyền. .. 1.2 NỐI ĐẤT CHO ĐƯỜNG DÂY TRUYỀN TẢI TRÊN THẾ GIỚI 13 1.2.1 Quy phạm cho hệ thống nối đất 13 1.2.2 Hình dạng nối đất .15 1.3 NGHIÊN CỨU ĐÁP ỨNG QUÁ ĐỘ CỦA HỆ THỐNG NỐI ĐẤT CỦA... NGUYỄN XUÂN PHÚC NGHIÊN CỨU ĐÁP ỨNG QUÁ ĐỘ CỦA HỆ THỐNG NỐI ĐẤT ĐỐI VỚI DÒNG ĐIỆN SÉT TRÊN ĐƯỜNG DÂY TRUYỀN TẢI VIỆT NAM Ngành: Kỹ thuật điện Mã số: 9520201 LUẬN ÁN TIẾN SĨ KỸ THUẬT ĐIỆN NGƯỜI HƯỚNG