(Luận văn) nghiên cứu ứng dụng lý thuyết hàm tương quan để xác định vị trí sự cố của đường dây tải điện

ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC KỸ THUẬT CÔNG NGHIỆP NGUYỄN SỸ HẢI lu an n va tn to gh NGHIÊN CỨU ỨNG DỤNG LÝ THUYẾT p ie HÀM TƯƠNG QUAN ĐỂ XÁC ĐỊNH VỊ TRÍ SỰ CỐ d oa nl w CỦA ĐƯỜNG DÂY TẢI ĐIỆN ll u nf va an lu oi m LUẬN VĂN THẠC SỸ KHOA HỌC z at nh KỸ THUẬT ĐIỆN z m co l gm @ an Lu Thái Nguyên - Năm 2020 n va ac th si ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC KỸ THUẬT CÔNG NGHIỆP NGUYỄN SỸ HẢI lu an NGHIÊN CỨU ỨNG DỤNG LÝ THUYẾT n va HÀM TƯƠNG QUAN ĐỂ XÁC ĐỊNH VỊ TRÍ SỰ CỐ gh tn to CỦA ĐƯỜNG DÂY TẢI ĐIỆN p ie NGÀNH: KỸ THUẬT ĐIỆN d oa nl w MÃ SỐ: 8.52.02.01 LUẬN VĂN THẠC SỸ KHOA HỌC ll u nf va an lu KỸ THUẬT ĐIỆN oi m z at nh NGƯỜI HƯỚNG DẪN KHOA HỌC: TS Đỗ Trung Hải z m co l gm @ an Lu Thái Nguyên – Năm 2020 n va ac th si CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc lập – Tự – Hạnh phúc BẢN XÁC NHẬN CHỈNH SỬA LUẬN VĂN THẠC SĨ Họ tên tác giả luận văn: Nguyễn Sỹ Hải Đề tài luận văn: Nghiên cứu ứng dụng lý thuyết hàm tương quan để xác định vị trí cố đường dây tải điện Chuyên ngành:Kỹ thuật điện Mã số: : 8.52.02.01 lu an Tác giả, Cán hướng dẫn khoa học Hội đồng chấm luận văn xác nhận tác n va giả sửa chữa, bổ sung luận văn theo biên họp Hội đồng ngày 04/10/2020 với - Sửa lỗi câu chữ, thuật ngữ luận văn phù hợp với chuyên ngành gh tn to nội dung sau: p ie - Bổ sung thêm danh mục hình vẽ bảng biểu theo quy định w - oa nl Thái Nguyên,ngày 26 tháng 10 năm 2020 Tác giả luận văn d Cán hướng dẫn u nf va an lu TS Đỗ Trung Hải Nguyễn Sỹ Hải ll oi m z at nh CHỦ TỊCH HỘI ĐỒNG z gm @ m co l PGS.TS Nguyễn Hữu Công an Lu n va ac th i si LỜI CAM ĐOAN Tơi xin cam đoan cơng trình nghiên cứu thân Các nghiên cứu kết trình bày luận văn trung thực chưa công bố công bố luận văn trước Thái Nguyên, ngày 20 tháng năm 2020 Tác giả luận văn lu an n va Nguễn Sỹ Hải p ie gh tn to d oa nl w ll u nf va an lu oi m z at nh z m co l gm @ an Lu n va ac th ii si LỜI CẢM ƠN Trong trình thực luận văn này, nhận nhiều ý kiến đóng góp, động viên từ thầy giáo, bạn đồng nghiệp người thân gia đình Lời đầu tiên, xin bày tỏ long biết ơn tới TS Đỗ Trung Hải tận tình hướng dấn, ln hỗ trợ khích lệ suốt thời gian làm luận văn để tơi hồn thành luận văn Tơi xin gửi lời cảm ơn chân thành tới tồn thể thầy giáo tham gia giảng dạy khóa học chuyên ngành Kỹ thuật điện cho ý kiến quý báu suốt q trình học tập lu an Tơi xin gửi lời cảm ơn tới Thầy giáo, Cô giáo khoa Điện Phòng Đào n va tạo Nhà trường tạo điều kiện thuận lợi mặt để tơi hồn thành nội tn to dung luận văn p ie gh Thái Nguyên, ngày 20 tháng năm 2020 d oa nl w HỌC VIÊN lu ll u nf va an Nguyễn Sỹ Hải oi m z at nh z m co l gm @ an Lu n va ac th iii si MỤC LỤC LỜI CAM ĐOAN II LỜI CẢM ƠN III MỤC LỤC IV DANH MỤC CÁC KÍ HIỆU VÀ CHỮ VIẾT TẮT VII DANH MỤC HÌNH VẼ, BẢNG BIỂU VIII MỞ ĐẦU 1.Tính cấp thiết Mục đích nghiên cứu .1 Nhiệm vụ nghiên cứu lu an Đối tượng phạm vi nghiên cứu n va Phương pháp nghiên cứu tn to Nội dung luận văn gh CHƯƠNG 1: TỔNG QUAN VỀ CÁC PHƯƠNG PHÁP ĐỊNH VỊ SỰ CỐ TRÊN p ie ĐƯỜNG DÂY TRUYỀN TẢI ĐIỆN .3 w 1.1 Đặt vấn đề oa nl 1.2 Một số phương pháp xác định vị trí cố d 1.2.1 Phương pháp tính tốn dựa trở kháng (đo lường phía) [9,10] an lu 1.2.1.1 Phương pháp điện kháng đơn u nf va 1.2.1.2 Phương pháp TAKAGI 1.2.1.3 Phương pháp TAKAGI cải tiến ll oi m 1.2.2 Phương pháp đo lường từ hai phía [11,12] cố z at nh 1.2.3 Phương pháp định vị cố dựa nguyên lý sóng lan truyền từ điểm z 1.2.4 Phương pháp định vị cố dựa nguyên lý sóng lan truyền từ đầu @ gm đường dây 10 l 1.3 Kết luận chương .10 m co CHƯƠNG 2: NGHIÊN CỨU QUÁ TRÌNH TRUYỀN SÓNG TRÊN ĐƯỜNG DÂY an Lu TẢI ĐIỆN 11 2.1 Mơ hình đường dây truyền tải điện [2,4] 11 n va ac th iv si 2.2 Nguyên lý lan truyền sóng đường dây [4,6] 13 2.2.1 Tổng trở sóng ZC 13 2.2.2 Hệ số truyền sóng  14 2.2.3 Vận tốc truyền sóng v 14 2.3 Sóng điện từ đường dây tải điện không cố [3,4] 15 2.3.1 Sóng lan truyền đường dây khơng có cố với tải cuối đường dây trở 17 2.3.2 Sóng lan truyền đường dây khơng có cố với tải cuối đường dây dạng (R nt L): 17 lu 2.3.3 Sóng lan truyền đường dây khơng có cố với tải cuối đường dây dạng an (R|| L): 18 va 2.3.4 Sóng lan truyền đường dây khơng có cố với tải cuối đường dây dạng n 2.3.5 Sóng lan truyền đường dây khơng có cố với tải cuối đường dây dạng ie gh tn to (R ||C): 19 p (R nt C): 19 w 2.4 Sóng điện từ đường dây tải điện có điểm cố .20 oa nl 2.5 Kết luận chương .21 d CHƯƠNG 3: ỨNG DỤNG LÝ THUYẾT HÀM TƯƠNG QUAN ĐỂ XÁC ĐỊNH VỊ lu an TRÍ SỰ CỐ TRÊN ĐƯỜNG DÂY TRUYỀN TẢI ĐIỆN 22 u nf va 3.1 Tín hiệu chirp [13,14] 22 3.2 Hàm tương quan [7] 25 ll oi m 3.3 Sử dụng hàm tương quan để xác định vị trí cố [1] 26 z at nh 3.4 Kết luận chương .30 CHƯƠNG 4: CÁC KẾT QUẢ TÍNH TỐN VÀ MƠ PHỎNG 31 z 4.1 Phần mềm Matlab-Simulink mô hệ thống điện [8] 31 @ gm 4.1.1 Thư viện Sources: gồm có khối nguồn constant, khối Ramp Step, l Khối Pulse Generator, Khối sine Wave hình 4.1 31 m co 4.1.2 Thư viện thị Sinks bao gồm khối Display, XY Graph, To File, To an Lu Workspace, Scope hình vẽ 32 4.1.3 Thư viện SimPowerSystems: Công cụ mô lưới điện 33 n va ac th v si 4.1.4 Thư viện nguồn: bao gồm nguồn chiều, xoay chiều pha, pha, nguồn dịng, nguồn áp, nguồn có điều khiển nguồn khơng có điều khiển hình vẽ 34 4.1.5 Thư viện Elements 34 4.1.6 Thư viện khối đo lường: 38 4.2 Mô sóng lan truyền đường dây dài để xác định vận tốc truyền sóng .38 4.2.1 Mơ hình mơ 38 4.2.2 Kết mô 41 4.3 Mơ sóng lan truyền đường dây dài không phân nhánh đường dây lu an có cố: 42 n va 4.3.1 Mơ hình mơ 42 4.3.2.1 Sự cố ngắn mạch ba pha chạm đất 43 4.3.2.2 Sự cố pha chạm đất 45 p ie gh tn to 4.3.2 Kết mô số loại cố khác nhau: 43 w 4.3.2.3 Sự cố ngắn mạch hai pha 45 oa nl 4.4 Kết luận chương .46 d KẾT LUẬN 47 an lu TÀI LIỆU THAM KHẢO .48 u nf va PHỤ LỤC 1: CHƯƠNG TRÌNH XÁC ĐỊNH VẬN TỐC TRUYỀN SÓNG 50 PHỤ LỤC 2: CHƯƠNG TRÌNH XÁC ĐỊNH VỊ TRÍ SỰ CỐ 52 ll oi m z at nh z m co l gm @ an Lu n va ac th vi si DANH MỤC CÁC KÍ HIỆU VÀ CHỮ VIẾT TẮT Ý nghĩa Ký hiệu Đơn vị lu an R0 Điện trở đơn vị chiều dài đường dây /km L0 Điện cảm đơn vị chiều dài đường dây H/km H0 Điện dung đơn vị chiều dài đường dây F/km G0 Điện dẫn đơn vị chiều dài đường dây S/km Vref Sóng tín hiệu điện áp phản hồi V Vinc Sóng tín hiệu điện áp chiều có biên độ Vinc (sóng tới) V Vận tốc truyền sóng đường dây truyền tải điện I Dòng điện A Chiều dài đường dây km Chiều dài từ đầu đường dây đến điểm cố km w  Km/s n va V gh tn to l p ie Lfault Điện trở cố IF Dòng điến cố ZL Tổng trở đường dây  hệ số khúc xạ  hệ số phản xạ oa nl Rf d A lu ll u nf va an  oi m z at nh z m co l gm @ an Lu n va ac th vii si DANH MỤC HÌNH VẼ Hình 1.1: Sơ đồ minh họa cố đường dây truyền tải sử dụng phương pháp điện kháng đơn Hình 1.2: Minh họa phương pháp TAKAGI mạch điện pha hai nguồn Hình 1.3: Sơ đồ nguyên lý đường dây bị cố với hai nguồn cấp Hình 1.4: Sơ đồ thay đường dây cố Hình 1.5: Sự lan truyền phản xạ sóng dịng điện đường dây 10 Hình 1: Mơ hình đường dây truyền tải hình PI pha 11 Hình 2.2: Mơ hình phân đoạn đường dây truyền tải hình PI ba pha 12 Hình 2.3: Mơ hình Petersen tương đương để giải tốn truyền sóng 16 lu an Hình 2.4: Mơ hình Petersen tương đương mạch có tải trở 17 n va Hình 2.5: Mơ hình Petersen tương đương mạch có tải R nối tiếp L .17 tn to Hình 2.6: Mơ hình Petersen tương đương mạch có tải R song song L 18 gh Hình 2.7: Mơ hình Petersen tương đương mạch R song song C 19 p ie Hình 2.8: Mơ hình Petersen tương đương mạch R nối tiếp C 19 w Hình 3.1: Sơ đồ khối tổng quan phương pháp xác định vị trí cố đường dây oa nl truyền tải điện 22 d Hình 3.2 tín hiệu chirp 25 an lu Hình 3.3: Tín hiệu đầu đường dây đo khơng có cố 27 u nf va Hình 3.4: Đồ thị hàm tương quan tín đầu tới tín hiệu phản hồi đo khơng có cố .27 ll oi m Hình 3.5: Lưu đồ thuật tốn xác định thời điểm sóng phản hồi 30 z at nh Hình 4.1 Thư viện khối nguồn 31 Hình 4.2 Thư viện khối hiển thị 33 z Hình 4.3 Thư viện cơng cụ mơ SimPowerSystems 33 @ gm Hình 4.4 Thư viện khối nguồn SimPowerSystems 34 l Hình 4.5 Thư viện Elements SimPowerSystems 35 m co Hình 4.6 Block cài đặt thơng số cho đường dây thông số dải .36 an Lu Hình 4.7 Block cài đặt thơng số cho máy cắt pha .36 Hình 4.8 Block cài đặt thông số cho cổng kết nối 37 n va ac th viii si lu an va n Hình 4.10 Mơ hình khối đo lường p ie gh tn to d oa nl w ll u nf va an lu m oi Hình 4.11 Mơ hình khối phát xung tín hiệu hình chirp vào đầu đường dây z at nh Để mơ tín hiệu phát xung vào đầu đường dây luận văn sử dụng khối tín hiệu hình chirp nối với khối tín hiệu Step khối tín hiệu Signal z gm @ Generator thơng qua khối nhân tín hiệu Trong mơ hình tín hiệu chirp với tần số bắt đầu 400 khz tần số kết thúc l m co 500khz với thời gian 50.10-6(s) Hình 4.12 Tín hiệu Signal Generator có tần số 100.102Hz hình 4.13 Với tín hiệu chirp biểu diễn theo hàm số: an Lu x(t)  cos(2 f ) cos  2f (t)t    n va ac th 39 si Trong tần số biến đổi xác định: f t  k t  f0 lu an n va gh tn to Hình 4.11 Mơ hình thơng số tín hiệu Hình 4.12 Mơ hình thơng số tín hiệu chirp ie chirp p Với mơ hình thiết bị phát xung, mơ hình khối đo lường mơ hình đường nl w dây trình bày luận văn xây dựng mơ hình mơ thành phần d 4.13 oa sóng lan truyền phản xạ đường dây pha với nguồn phát xung chirp hình ll u nf va an lu oi m z at nh z m co l gm @ Hình 4.13: Mơ hình mơ xác định thành phần sóng lan truyền phản xạ an Lu đường dây pha với nguồn phát xung chirp n va ac th 40 si 4.2.2 Kết mơ Sử dụng mơ Hình 4.13 ta gửi tín hiệu chirp vào đầu đường dây truyền tải, tín hiệu đo đầu đường dây khơng có cố hình 4.14 Phân tích hàm tương quan dựa sở tính giá trị tương quan tín hiệu tới (tín hiệu mẫu chuẩn) với đoạn tín hiệu phản xạ đo đồ thị hàm tương quan Hình Khi hàm tương quan tín hiệu phản xạ tín hiệu mẫu xấp xỉ nói thời gian tín hiệu phản xạ Biết thời gian sóng phản xạ chiều dài đường dây tính tốn vận tốc truyền sóng 0.8 0.6 lu 0.4 an 0.2 va n tn to -0.2 gh -0.4 p ie -0.6 -0.8 Time(s) -4 nl w x 10 d oa Hình 4.14: Tín hiệu đầu đường dây đo khơng có cố an lu va 0.8 u nf 0.6 ll 0.4 m 0.2 oi z at nh -0.2 -0.4 -0.6 z 1000 2000 3000 4000 Time(s) gm -1 @ -0.8 5000 6000 7000 8000 l khơng có cố m co Hình 4.15: Đồ thị hệ số tương quan tín hiệu tới tín hiệu phản xạ đo an Lu n va ac th 41 si Từ mơ hình đường dây chương trình so sánh hệ số tương quan phụ lục 1, ta xác định thời điểm sóng phản xạ từ cuối đường dây hệ số tương quan điện áp đầu đường dây tín hiệu chirp R  0,99 Khoảng thời gian từ ban đầu tới thời điểm là: T0  601,3(  s) nên ta có vận tốc sóng chạy đường dây xét là: v 2l  232828,8  km / s  T0 4.3 Mô sóng lan truyền đường dây dài khơng phân nhánh đường dây có cố: lu 4.3.1 Mơ hình mô an Với đường dây cố luận văn sử dụng phần tử Simulink để xây dựng n va khối mô ngắn mạch pha, ngắn mạch pha ngắn mạch pha chạm đất to p ie gh tn hình vẽ oa nl w d Hình 4.8: Sơ đồ khối mơ dạng ngắn mạch lu ll u nf va an Khi đường dây có cố luận văn sử dụng mơ hình mơ hình vẽ oi m z at nh z m co l gm @ an Lu Hình 4.17: Sơ đồ mơ đường dây cố n va ac th 42 si Các giá trị thông số sau cố ngắn mạch sử dụng để mô phỏng: - Điện trở cố Rfault: R = 20  - Điện cảm cố: Lfaulf: L= 1mH - Loại cố: - Chạm đất pha - Ngắn mạch pha - Ngắn mạch pha chạm đất 4.3.2 Kết mô số loại cố khác nhau: lu 4.3.2.1 Sự cố ngắn mạch ba pha chạm đất an Sử dụng mơ hình 4.17 với cố ngắn mạch pha chạm đất 40km ta va thu tín hiệu đo đầu đường dây hình 4.18 Sử dụng hàm tương quan n tn to thu đồ thị hệ số tương quan tín tới tín hiệu phản xạ đo cố ie gh pha chạm đất 40 km hình 4.20 p 0.8 0.6 nl w 0.4 oa d 0.2 an lu va u nf -0.2 ll -0.4 oi m -0.6 z at nh -0.8 Time(s) -4 x 10 z 40km m co l gm @ Hình 4.18: Tín hiệu đầu đường dây có cố ngắn mạch pha chạm đất an Lu n va ac th 43 si 0.8 0.6 0.4 0.2 -0.2 -0.4 -0.6 -0.8 0.5 1.5 Time(s) 2.5 3.5 -4 x 10 Hình 4.19: Tín hiệu phóng to đầu đường dây có cố ngắn mạch pha lu an chạm đất 40km va n tn to 0.8 0.6 gh 0.4 ie p 0.2 w oa nl -0.2 -0.4 d -0.6 lu 1000 2000 3000 4000 5000 6000 7000 8000 u nf va -1 an -0.8 ll Hình 4.20: Đồ thị hệ số tương quan tín tới tín hiệu phản xạ đo khơng m oi có cố pha chạm đất 40km z at nh Trên hình 4.18 hình phóng to 4.19 nhận thấy có tín hiệu hình chirp z phản hồi từ vị trí cố từ cuối đường dây Trên hình 4.20 đồ thị hàm tương gm @ quan tín hiệu tới tín hiệu phản hồi Trong có thời điểm sóng phản hồi từ vị l trí cố thời điểm sóng phản hồi từ cuối đường dây có hệ số tương quan cao m co Thực tương tự điểm cố: 20Km, 30 Km, 50 Km, 60 Km ta thu kết bảng 4.1 an Lu n va ac th 44 si Bảng 1: Bảng kết xác định vị trí cố pha chạm đất áp dụng phương pháp hàm tương quan Lfault Dạng cố Hệ số R L (km) Sai số(m) 20 ABCG 0.994 19.83 170 30 ABCG 0.993 29.815 185 40 ABCG 0.992 39.81 190 50 ABCG 0.993 49.78 220 60 ABCG 0.992 59.79 210 4.3.2.2 Sự cố pha chạm đất Thực tương tự với cố pha chạm đất 20Km, 30Km, 40Km, 50 Km, 60 lu Km ta bảng kết bảng 4.2 an n va Bảng 4.2: Bảng kết xác định vị trí cố pha chạm đất áp dụng phương pháp hàm tương quan Dạng cố Hệ số R L (km) Sai số(m) 20 AG 0.934 19.932 68 AG 0.94 29.96 40 AG 0.93 39.91 90 AG 0.93 49.912 88 AG 0.94 59.92 80 ie gh tn to Lfault p 30 40 d oa nl 60 w 50 an lu va 4.3.2.3 Sự cố ngắn mạch hai pha u nf Thực tương tự với cố ngắn mạch pha 20Km, 30Km, 40Km, 50 Km, ll 60 Km ta bảng kết bảng 4.3 m oi Bảng 4.3: Bảng kết xác định vị trí cố ngắn mạch pha áp dụng Lfault Dạng cố Hệ số R Sai số(m) 20 AB 0.93 19.98 20 30 AB 0.94 29.965 35 40 AB 0.93 45 50 AB 0.93 49.952 48 60 AB 0.92 59.93 70 z L (km) l z at nh phương pháp hàm tương quan gm @ 39.955 m co an Lu n va ac th 45 si 4.4 Kết luận chương Trong chương mơ q trình truyền sóng đường dây 110kV Sơn La đoạn từ E17.1 Mộc Châu E19.6 Mai Châu với chiều dài đường dây 71 km sử dụng dây dẫn AC185 Sử dụng phương pháp phân tích hệ số tương quan tín hiệu tới tín hiệu phản xạ để phát thời điểm sóng phản xạ từ cuối đường dây hay từ điểm cố Căn vào thời điểm phát xung thời điểm sóng phản xạ từ điểm cố để định vị cố đường dây không phân nhánh với thay đổi vị trí cố dạng cố lu an n va p ie gh tn to d oa nl w ll u nf va an lu oi m z at nh z m co l gm @ an Lu n va ac th 46 si KẾT LUẬN Luận văn nghiên cứu phương pháp định vị cố dựa nguyên lý sóng lan truyền từ đầu đường dây Từ mơ hình đường dây tải điện, nguyên lý lan truyền sóng thơng số q trình truyền sóng; phân tích sóng điện từ đường dây truyền tải điện trường hợp đường dây không cố đường dây cố Sử dụng phương pháp phân tích hệ số tương quan tín hiệu tới tín hiệu phản xạ để phát điểm cố Các kết mơ q trình truyền sóng, xác định vị trí cố thực đường dây 110kV Sơn La đoạn từ E17.1 Mộc Châu E19.6 Mai Châu lu an với chiều dài đường dây 71 km sử dụng dây dẫn AC185 n va Phương pháp hàm tương quan cần đo lường từ đầu đường dây khơng u tn to cầu đồng tín hiệu mơ xác định vị trí cố đường dây không phân nhánh gh với dạng cố: p ie - Chạm đất pha w - Ngắn mạch hai pha oa nl - Ngắn mạch ba pha chạm đất d Ở vị trí cố khác cho kết tốt, sai số khoảng cách cố thực ll u nf va an lu xác định phương pháp nghiên cứu chấp nhận oi m z at nh z m co l gm @ an Lu n va ac th 47 si TÀI LIỆU THAM KHẢO Dương Hòa An (2019), " Nghiên cứu ứng dụng phương pháp đại nhận dạng cố ngắn mạch đường dây truyền tải điện ", luận án tiến sĩ kỹ thuật, Trường Đại học Kỹ thuật Công nghiệp Trần Bách (2004), "Lưới điện Hệ thống điện" tập & 2, NXB Khoa học Kỹ thuật, Hà Nội Lại Khắc Lãi (2009), "Cơ sở lý thuyết mạch tập 2", NXB Đại học Thái Nguyên Nguyễn Bình Thành (1978), "Giáo trình Cơ sở kỹ thuật điện 1,2", NXB ĐHBK Hà Nội lu an Thành Lương (2013), “Công nghệ định vị cố”, Tạp chí Khoa học Cơng nghệ n va Điện, số tn to Trần Đình Long (2000),” Bảo vệ hệ thống điện”, NXB Khoa học Kỹ thuật gh Nguyễn Thanh Sơn Lê Khánh Luận, "Lý thuyết Xác suất Thống kê", NXB Đại học p ie Quốc gia TPHCM oa nl 2003 w Nguyễn Phùng Quang (2003), "Matlab & Simulink", NXB Khoa học Kỹ thuật, d Takagi (1982), "Development of a new fault locator using the one-terminal lu u nf va pp 2892–2898 an voltage and current data", IEEE Trans.Power App Syst , vol PAS-101, no 8, 10 Karl Zimmerman and David Costello, "Impedance-Based Fault Location ll oi m Experience ", Rural Electric Power Conference, 2006 IEEE z at nh 11 M Kezunovic and B Perunicic (1996), "A utomated transmission line fault analysisusing synchronized sampling at two ends", IEEE Trans Power Syst , vol z 11, no 1, pp 441–447 @ gm 12 Minambres J Zamora I, Mazon A, Alvarez-Isasi R, Lazaro J (1996), "Fault m co Generation, Transmission and distribution, 143(1) l location on two-terminal transmission lines based on voltages", IEE Proceedings: an Lu n va ac th 48 si 13 N G Paulter ( 2001), "An assessment on the accuracy of time-domain reflectometry for measuring the characteristic impedance of transmission line", IEEE Trans Instrum Meas., vol 50, pp 1381 –1388, 2001 14 M Ohmiya H Yamada, Y Ogawa, K Itoh (1991), "Super resolution techniques for time-domain measurements with a network analyzer", IEEE Trans Antennas Propag., vol 39, pp 177 –183, 1991 lu an n va p ie gh tn to d oa nl w ll u nf va an lu oi m z at nh z m co l gm @ an Lu n va ac th 49 si PHỤ LỤC 1: CHƯƠNG TRÌNH XÁC ĐỊNH VẬN TỐC TRUYỀN SĨNG tao file d? li?u xu?t h́ nh t=ScopeData7.time; y=ScopeData7.signals.values; %plot(ScopeData7.time,ScopeData7.signals.values) grid on; xlabel('Time(s)') print -dtiff -r300 khongsuco Coeff = zeros(1,length(x1)-length(mau)+1); for i=1:(length(x1)-length(mau)+1) lu an tmp=x1(i:(i+length(mau)-1)); va tmp_coeff=corrcoef(mau,tmp); n tn to Coeff(i) = tmp_coeff(1,2); gh end; p ie figure(2) w plot(Coeff); grid on d y1=Coeff; oa nl % xac dinh gia tri thoi diem t1 ll u nf va % figure(3) an mauy1=y1(1:1000); lu n=length(y1); z if y1(i)==maxy1 z at nh for i=1:1:1000 oi maxy1=max(mauy1); m %plot(mauy1) @ gm t1=i; l end m co end an Lu % xac dinh gia tri thoi diem t2 y2=Coeff; n va ac th 50 si n=length(y2); y2=y2(1000:n); %figure(4) %plot(y2) maxy2=max(y2); for i=1000:1:n if y1(i)==maxy2 t2=i; end end lu an % va detat2=(t2-t1)*1000/10000; % ns n tn to v0=2*46.7/detat2*10^6 p ie gh % -xac dinh van toc d oa nl w ll u nf va an lu oi m z at nh z m co l gm @ an Lu n va ac th 51 si PHỤ LỤC 2: CHƯƠNG TRÌNH XÁC ĐỊNH VỊ TRÍ SỰ CỐ tao file d? li?u xu?t h́ nh t=ScopeData7.time; y=ScopeData7.signals.values; %plot(ScopeData7.time,ScopeData7.signals.values) grid on; xlabel('Time(s)') print -dtiff -r300 khongsuco Coeff = zeros(1,length(x1)-length(mau)+1); for i=1:(length(x1)-length(mau)+1) lu an tmp=x1(i:(i+length(mau)-1)); va tmp_coeff=corrcoef(mau,tmp); n tn to Coeff(i) = tmp_coeff(1,2); gh end; p ie figure(2) w plot(Coeff); grid on d y1=Coeff; oa nl % xac dinh gia tri thoi diem t1 ll u nf va % figure(3) an mauy1=y1(1:1000); lu n=length(y1); z if y1(i)==maxy1 z at nh for i=1:1:1000 oi maxy1=max(mauy1); m %plot(mauy1) @ gm t1=i; l end m co end an Lu % xac dinh gia tri thoi diem t2 y2=Coeff; n va ac th 52 si n=length(y2); y2=y2(1000:n); %figure(4) %plot(y2) maxy2=max(y2); for i=1000:1:n if y1(i)==maxy2 t2=i; end end lu an % va detat=(t2-t1)*1000/10000; % ns n p ie gh tn to Lfault=detat*v0/2*10^(-6) d oa nl w ll u nf va an lu oi m z at nh z m co l gm @ an Lu n va ac th 53 si