Phương pháp nhận dạng và hiệu chỉnh các mạch vòng điều khiển nhà máy nhiệt điện đốt than phun

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI PHƯƠNG PHÁP NHẬN DẠNG VÀ HIỆU CHỈNH CÁC MẠCH VÒNG ĐIỀU KHIỂN NHÀ MÁY NHIỆT ĐIỆN ĐỐT THAN PHUN LUẬN ÁN TIẾN SĨ KỸ THUẬT NHIỆT Hà Nội - 2022 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI PHƯƠNG PHÁP NHẬN DẠNG VÀ HIỆU CHỈNH CÁC MẠCH VÒNG ĐIỀU KHIỂN NHÀ MÁY NHIỆT ĐIỆN ĐỐT THAN PHUN Ngành: Kỹ thuật nhiệt Mã số: 9520115 LUẬN ÁN TIẾN SĨ KỸ THUẬT NHIỆT NGƯỜI HƯỚNG DẪN KHOA HỌC: TS Lê Đức Dũng PGS TS Bùi Quốc Khánh Hà Nội - 2022 LỜI CAM ĐOAN Tôi xin cam đoan luận án “Phương pháp nhận dạng hiệu chỉnh mạch vòng điều khiển nhà máy nhiệt điện đốt than phun” nghiên cứu hoàn thành hướng dẫn Các số liệu nêu luận án hoàn toàn trung thực Những kết luận án chưa cơng bố cơng trình khác Hà Nội, ngày 04 tháng 01 năm 2022 Tập thể hướng dẫn Nghiên cứu sinh i LỜI CẢM ƠN Trước hết, tơi xin bày tỏ lịng biết ơn chân thành sâu sắc tới tập thể hướng dẫn tận tâm, tận lực hỗ trợ, động viên, hướng dẫn mặt chuyên môn suốt trình tơi thực luận án Tơi xin chân thành cảm ơn giúp đỡ, góp ý, thời gian thầy cô Viện Khoa học Công nghệ Nhiệt-Lạnh, Trường Đại học Bách khoa Hà Nội hướng dẫn báo cho vấn đề có liên quan đến luận án Tơi xin gửi lời cảm ơn tới Ban Giám hiệu Trường Đại học Bách khoa Hà Nội, lãnh đạo Viện liên quan hỗ trợ tơi hồn thành chương trình học nghiên cứu sinh Tôi chân thành cảm ơn Ban lãnh đạo Tổng công ty Phát điện hỗ trợ nhiều điều kiện trình tìm hiểu thực tế bố trí cơng việc hợp lý để tơi hồn thành chương trình học tập nghiên cứu sinh Cuối cùng, tơi xin cảm ơn gia đình, bạn bè người thân bên cạnh quan tâm, động viên tạo điều kiện thuận lợi để giành thời gian thực luận án, đồng thời ln sẻ chia khó khăn động viên tơi q trình nghiên cứu Hà Nội, ngày 04 tháng 01 năm 2022 Nghiên cứu sinh ii MỤC LỤC LỜI CAM ĐOAN i LỜI CẢM ƠN ii MỤC LỤC iii DANH MỤC CÁC KÝ HIỆU VÀ CHỮ VIẾT TẮT vi DANH MỤC CÁC BẢNG BIỂU viii DANH MỤC CÁC HÌNH VẼ, ĐỒ THỊ ix MỞ ĐẦU CHƯƠNG KHÁI QUÁT CHUNG VỀ NHẬN DẠNG VÀ CHỈNH ĐỊNH CÁC BỘ ĐIỀU KHIỂN TRONG NHÀ MÁY NHIỆT ĐIỆN 11 NGUYÊN LÝ LÀM VIỆC CƠ BẢN NHÀ MÁY NHIỆT ĐIỆN ĐỐT THAN PHUN 111 Lò nhà máy nhiệt điện đốt than phun 112 Tuabin nhà máy nhiệt điện 12 NGUYÊN NHÂN DẪN ĐẾN CẦN PHẢI HIỆU CHỈNH VÀ CHỈNH ĐỊNH CÁC BỘ ĐIỀU KHIỂN TRONG NHÀ MÁY NHIỆT ĐIỆN ĐỐT THAN PHUN CÔNG NGHỆ CẬN TỚI HẠN 121 Đặc trưng trình nhiệt nhà máy nhiệt điện 122 Công tác sửa chữa bảo dưỡng hệ thống thiết bị 123 Thành phần nhiên liệu than 124 Những ví dụ thực trạng vận hành NMNĐ Duyên Hải 10 13 KHÁI QUÁT CHUNG VỀ PHƯƠNG PHÁP NHẬN DẠNG THAM SỐ CỦA CÁC MẠCH VỊNG ĐIỀU KHIỂN TRONG CƠNG NGHIỆP 14 131 Phương pháp mơ hình đối tượng 14 132 Ngun tắc nhận dạng đối tượng điều chỉnh công nghiệp 16 14 LÝ THUYẾT CHUNG VỀ PHƯƠNG PHÁP HIỆU CHỈNH VÀ CHỈNH ĐỊNH CÁC MẠCH VÒNG ĐIỀU KHIỂN 17 141 Chỉnh định hệ thống điều khiển 18 142 Đánh giá mạch vòng cần chỉnh PID 23 143 Hiệu chỉnh lượng đặt, thơng số khởi tạo tồn q trình nhiệt nhà máy điện theo thời gian 24 15 KHÁI QT CÁC CƠNG TRÌNH NGHIÊN CỨU VỀ NHẬN DẠNG VÀ CHỈNH ĐỊNH CÁC BỘ ĐIỀU KHIỂN TRONG CÔNG NGHIỆP 25 151 Khái quát nhận dạng đối tượng điều khiển công nghiệp 25 152 Khái quát chỉnh định PID NMNĐ 26 153 Nghiên cứu ứng dụng mạng nơron nhận dạng điều khiển nhà máy nhiệt điện đốt than 31 iii 16 ĐỊNH HƯỚNG NGHIÊN CỨU VỀ NHẬN DẠNG VÀ CHỈNH ĐỊNH HỆ ĐIỀU KHIỂN TRONG NHÀ MÁY NHIỆT ĐIỆN ĐỐT THAN PHUN CÔNG NGHỆ CẬN TỚI HẠN 34 1KẾT LUẬN CHƯƠNG 36 CHƯƠNG 37 PHÂN TÍCH YẾU TỐ ẢNH HƯỞNG TỚI CÁC MẠCH VÒNG ĐIỀU KHIỂN CƠ BẢN TRONG NHÀ MÁY NHIỆT ĐIỆN ĐỐT THAN PHUN 37 CHƯƠNG 37 2SƠ ĐỒ CÁC MẠCH VÒNG ĐIỀU KHIỂN CƠ BẢN CỦA NHÀ MÁY NHIỆT ĐIỆN THAN 37 2PHÂN TÍCH MẠCH VỊNG CẤP NƯỚC 40 2PHÂN TÍCH MẠCH VỊNG ĐIỀU KHIỂN KHĨI GIÓ 44 3 Thành phần vai trị hệ thống khói gió 44 Mạch vòng điều khiển nồng độ oxy dư khói thải 47 3 Mạch vịng điều khiển áp suất âm buồng đốt 49 PHÂ G ĐIỀU KHIỂN N NHIỆT ĐỘ HƠI TÍCH 51 MẠC HỆ ĐIỀU KHIỂN H PHỤ TẢI NHIỆT VÒN 54 25 Ph ân tích cấu trú c điề u ển phụ tải nhi ệt 56 25 Ph ân tích mạ ch vịn g điề u ển phụ tải nhi ệt 61 2ĐÁNH GIÁ TÁC ĐỘNG CỦA CÁC THỐNG SỐ QUÁ TRÌNH 68 Mô hệ điều khiển 68 Kiểm chứng kết hồn thiện mơ hình theo số liệu NMNĐ Duyên Hải 70 Đánh giá chất lượng điều khiển 72 Đánh giá tác động thơng số q trình tới đáp ứng hệ thống 73 1KẾT LUẬN CHƯƠNG 77 CHƯƠNG 78 NGHIÊN CỨU SỬ DỤNG PHẦN TỬ DỰ BÁO KẾT HỢP BỘ ĐIỀU KHIỂN ĐỂ CHỈNH ĐỊNH CÁC MẠCH VÒNG ĐIỀU KHIỂN 78 3PHƯƠNG PHÁP KẾT HỢP KHÂU DỰ BÁO VỚI BỘ ĐIỀU KHIỂN 78 1 Phương pháp luận 78 Thiết kế hệ thống điều khiển kết hợp dự báo theo tín hiệu số 80 3PHẦN MỀM HỖ TRỢ CHỈNH ĐỊNH THAM SỐ BỘ ĐIỀU KHIỂN 85 Chức phần mềm 85 2 Các cửa sổ phần mềm 87 3ỨNG DỤNG PHẦN MỀM CHỈNH THAM SỐ BỘ ĐIỀU KHIỂN 89 3 Thiết kế điều khiển PI 89 3 Thiết kế điều khiển PID dạng thực 92 3 Thiết kế mạch vòng điều khiển nồng độ oxy dư tổ máy S1 NMNĐ Duyên Hải 94 iv 3 Thiết kế điều khiển Fuel Master Boiler Master tổ máy S1 NMNĐ Duyên Hải 96 34 CHUYỂN ĐỔI THAM SỐ BỘ ĐIỀU KHIỂN VÀO THỰC TẾ 99 Cấu trúc điều khiển thực tế NMNĐ Duyên Hải 99 Phương pháp chuyển đổi tham số điều khiển 100 35 MỘT SỐ NỘI DUNG CẦN PHÁT TRIỂN TRONG TƯƠNG LAI 102 36 KẾT LUẬN CHƯƠNG 103 CHƯƠNG 104 ỨNG DỤNG MẠNG NƠRON NHÂN TẠO NHẬN DẠNG VÀ HIỆU CHỈNH HỆ THỐNG ĐIỀU KHIỂN TRONG NHÀ MÁY NHIỆT ĐIỆN ĐỐT THAN 104 41 VẤN ĐỀ NGHIÊN CỨU ỨNG DỤNG MẠNG NƠRON TRONG NHÀ MÁY NHIỆT ĐIỆN ĐỐT THAN 104 1 Nguyên lý ứng dụng mạng nơron nhân tạo nhận dạng hiệu chỉnh hệ điều khiển nhà máy nhiệt điện đốt than 104 Phương án triển khai nghiên cứu thực tế 106 42 ỨNG DỤNG MẠNG NƠRON NHẬN DẠNG TRẠNG THÁI VẬN HÀNH HỆ ĐIỀU KHIỂN NHÀ MÁY NHIỆT ĐIỆN ĐỐT THAN 106 Thiết kế mạng nơron nhận dạng hệ điều khiển 106 2 Kết mô 109 43 ỨNG DỤNG MẠNG NƠRON ĐỂ HIỆU CHỈNH LƯỢNG ĐẶT HỆ ĐIỀU KHIỂN PHỤ TẢI NHIỆT NHÀ MÁY NHIỆT ĐIỆN ĐỐT THAN 112 Điều chỉnh lưu lượng than nhiệt trị thay đổi 112 Phương pháp luận sử dụng mạng nơron 114 3 Nhận dạng thời gian thực 116 4 Tác động theo thời gian dự báo 119 44 ỨNG DỤNG MẠNG NƠRON ĐÁNH GIÁ CẤU TRÚC HỆ ĐIỀU KHIỂN PHỤ TẢI NHIỆT 120 4 Kịch mô 120 4 Kết Matlab/Simulink 122 45 KẾT LUẬN CHƯƠNG 126 KẾT LUẬN VÀ KIẾN NGHỊ 127 DANH MỤC CÁC CƠNG TRÌNH ĐÃ CƠNG BỐ CỦA LUẬN ÁN 128 TÀI LIỆU THAM KHẢO 129 v DANH MỤC CÁC KÝ HIỆU VÀ CHỮ VIẾT TẮT vi STT Ký hiệu Ý nghĩa a0 ,a1 , ,an α A(ω) BF CCS C CFB Lị tầng sơi (Circulating Fluidized Boiler) c,u,ω Vận tốc tuyệt đối, vận tốc quay, vận tương tương đối tuabin DCS Hệ thống điều khiển phân tán (Distributed control system) 10 ε Sai lệch điều khiển 11 i Biến đếm, chạy thuật toán 12 IDF 13 j 14 F (ω ) Hàm hỗ trợ xây dựng tham số điều khiển 15 FDF Quạt gió cấp 2, gió (Forced Draft Fan) 16 γ Tần số phi thứ nguyên 17 G Hàm truyền trình 18 Gc ( s) Hàm truyền điều khiển (controller) 19 Gp ( s) Hàm truyền đối tượng trình (process) 20 GD ( s ) Hàm truyền nhiễu trình 21 H(s) 22 Hf Nhiệt trị làm việc nhiên liệu 23 Im Phần ảo số thực đặc tính tần số 24 L Cơng có ích 25 m0 Hằng số tỉ lệ, tích phân, vi phân điều khiển PID 26 µ Tín hiệu điều khiển 27 m0 Chỉ số dao động cứng Hằng số quán tính mẫu thức Hệ số hàm truyền Biên độ dao động hệ kín Chế độ điều khiển lò theo tuabin Hệ thống điều khiển phối hợp lò – tuabin (Coordinated control system) Nhiệt dung riêng vật chất Quạt khói (Induced Draft Fan) Đơn vị ảo j = −1 Hàm truyền hệ hở this this this this this this tbxTd_F_S Text = TdR ToString(); tbxTf_F_S Text = TfR ToString(); nudKp_Y_S Value = Math Round(Convert ToDecimal(kпR), nudTi_Y_S Value = Math Round(Convert ToDecimal(TiR), nudTd_Y_S Value = Math Round(Convert ToDecimal(TdR), nudTf_Y_S Value = Math Round(Convert ToDecimal(TfR), 3); 3); 3); 3); //Kp_S, Ti_S, Td_S, Tf_S, MR_S; Kp_S = kпR; Ti_S = TiR; Td_S = TdR; Tf_S = TfR; MR_S = MR; toolTip SetToolTip(this tbxKp_F_S, "EMERSON: Kp(PGANIN) = " + kпR); toolTip SetToolTip(this tbxTi_F_S, "EMERSON: ti(INTG) = " + TiR / kпR); toolTip SetToolTip(this tbxTd_F_S, "EMERSON: Kd(INTG) = " + TdR * kпR); toolTip SetToolTip(this tbxTf_F_S, "EMERSON: td(INTG) = " + TfR); toolTip SetToolTip(this nudKp_Y_S, "EMERSON: Kp(PGANIN) = " + kпR); toolTip SetToolTip(this nudTi_Y_S, "EMERSON: ti(INTG) = " + TiR / kпR); toolTip SetToolTip(this nudTd_Y_S, "EMERSON: Kd(INTG) = " + TdR * kпR); toolTip SetToolTip(this nudTf_Y_S, "EMERSON: td(INTG) = " + TfR); toolTip SetToolTip(this nudKp2_Y_K, "EMERSON: Kp(PGANIN) = " + kпR); toolTip SetToolTip(this nudTi2_Y_K, "EMERSON: ti(INTG) = " + TiR / kпR); toolTip SetToolTip(this nudTd2_Y_K, "EMERSON: Kd(INTG) = " + TdR * kпR); toolTip SetToolTip(this nudTf2_Y_K, "EMERSON: td(INTG) = " + TfR); } private void CodeGraph_F_S() // Hàm đặc tính tần số mạch hở ПИ и ПИД одноконтурного регулятора { // Vẽ vòng tròn trữ ổn định uint число_mr = 100; double MR = Convert ToDouble(this nudM_S Text); double[] MRRe = new double[число_mr]; for (int i = 0; i < число_mr; i++) { double a = / Convert ToDouble(число_mr); MRRe[i] = MRRe_Rotach(i * a, MR); } double[] MRIm = new double[число_mr]; for (int i = 0; i < число_mr; i++) { double a = / Convert ToDouble(число_mr); MRIm[i] = MRIm_Rotach(i * a, MR); } // Vẽ đặc tính cứng hệ hở double OMEGAstart = Convert ToDouble(this nudStartwFW_S Text); double OMEGAend = * Convert ToDouble(this nudEndwFW_S Text); UInt32 nw = * Convert ToUInt32(this nudPointFW_S Text); double DeltaOMEGA = (OMEGAend - OMEGAstart) / Convert ToDouble(nw); double OMEGAr = Convert ToDouble(this nudwFW_S Text); // Vẽ đặc tính cứng hệ hở double[] HRe_S = new double[nw]; // Mảng giá trị đặc tính cứng trục thực for (int i = 0; i < nw; i++) { double w = OMEGAstart + i * DeltaOMEGA; HRe_S[i] = GraphHRe_S(w); } double[] HIm_S = new double[nw]; // Mảng giá trị đặc tính cứng trục ảo for (int i = 0; i < nw; i++) { double w = OMEGAstart + i * DeltaOMEGA; HIm_S[i] = GraphHIm_S(w); 36 } ///////////////////////////////////////////////////////////// chart_F_S Series Clear(); // Vẽ vòng tròn trữ ổn định graph(chart_F_S, "Контур Ротач", MRRe, MRIm, colorLevels[23], typeLevels[0]); // Vẽ điểm (-1j,0) graphPoint(chart_F_S, "Точка-1j0", -1, 0, colorLevels[23]); // Vẽ đặc tính cứng hệ hở kiểm tra việc có tiếp xúc đường tròn dự trữ graph(chart_F_S, "Разомкнутой КЧХ", HRe_S, HIm_S, colorLevels[0], typeLevels[0]); } private void CodeGraph_Y_S() // Đáp ứng độ đầu { if (this ch_Y_S Checked == true) { double I2 = 0; double tstart = Convert ToDouble(this nudSt_t_S Text); double tend = Convert ToDouble(this nudEn_t_S Text); UInt32 nt = Convert ToUInt32(this nudP_y_S Text); double Deltat = tend / Convert ToDouble(nt); double[] YOt = new double[nt]; for (int i = 0; i < nt; i++) { double tO = tstart + i * Deltat; YOt[i] = tO; } double OMEGAend = Convert ToDouble(this nud_W_S Text); UInt32 nOMEGA = Convert ToUInt32(this nudP_W_S Text); double DeltaOMEGA = OMEGAend / Convert ToDouble(nOMEGA); double[] YOy = new double[nt]; for (int i = 0; i < nt; i++) { double y = 0; // Tính giá trị đáp ứng đầu hệ thống for (int ik = 0; ik < nOMEGA; ik++) { double ikk; if (ik == 0) ikk = 00001; else ikk = Convert ToDouble(ik); double a = 0; a = / Math PI * (MYRe_S(ikk * DeltaOMEGA) / ikk) * Math Sin(ikk * DeltaOMEGA * (Deltat * i)); y = y + a; } // Tính tiêu bình phương tích phân sai số điều khiển if (this cbx_S_TrP SelectedIndex == 0) { YOy[i] = y; I2 = I2 + YOy[i] * YOy[i] * Deltat; } else { YOy[i] = y; I2 = I2 + (YOy[i] - 1) * (YOy[i] - 1) * Deltat; } } if (this ch_MY_S Checked == true) { число_graficYO = число_graficYO + 1; graph(chart_Y_S, "YO" + число_graficYO, YOt, YOy, colorLevels[число_colorYO], typeLevels[0]); } 37 else { chart_Y_S Series Clear(); число_colorYO = 0; graph(chart_Y_S, "YO", YOt, YOy, colorLevels[число_colorYO], typeLevels[0]); } число_colorYO = число_colorYO + 1; this tbxI_S Text = String Format("{0:0 00000}", I2); } } PL4 HƯỚNG DẪN SỬ DỤNG PHẦN MỀM TÍNH TỐN THAM SỐ BỘ ĐIỀU KHIỂN ĐỐI TƯỢNG Q TRÌNH NHIỆT Thiết kế điều khiển PID theo phương pháp Rotach V Y Giả sử đối tượng cần điều khiển có hàm truyền dạng qn tính bậc hai có trễ sau [24]: −τ µ s 9,6e−18 s kµ e = Gp ( s ) = (T1s + 12s)(T+ 1) (54s + 1)2 Bằng phương pháp Mạnh N V , điều khiển thiết kế buộc phải dạng PID, đó, với phương pháp Rotach V Y thiết kế điều khiển PI Sử dụng phần mềm, thu hai điều khiển tương ứng theo hai phương pháp: G pM     + Td s  = 0, 642  +  Ti s   108s  ( s ) = k p 1 + + 27s     G pR ( s ) = k p 1 +  = 0,1511 + 74, 647s   Ti s   Trường hợp cần thiết kế điều khiển PID theo phương pháp Rotach V Y , bước thực cụ thể sau: - Tại cửa sổ “Begin” chọn: Phương pháp turning “Ротач Виталий Яковлевич”; mạch thiết kế hệ thống “Single Outline”; điều khiển “PID”; chọn mơ hình đối tượng “First order plus dead time”; nhập tham số hàm truyền đổi tượng Tích chuột vào nút “Next” nhấn “Enter” bàn phím để cửa sổ thiết kế vẽ đặc tính cứng mạch vịng hở - Tại cửa sổ “Method professor Ротач Виталий Яковлевич” nhập: tiêu trự tần số M Nhấn chọn “Build a support function” để vẽ hàm bổ sung Hiệu chỉnh tham số trục tọa độ chọn “Clear” để thu giao điểm đồ thị với trục tần số (trục hồnh) - Sử dụng chuột tính “Current position” để xác định giá trị tần số giao điểm - Nhập giá trị cộng hưởng vào ô “Resonance frequency” tích chọn “Build a support function” nút “Clear” để thu tham số điều khiển PID đặc tính cứng hệ hở tiếp xúc với đường tròn dự trữ dao động Gc ( s )p =1k +     + + Td s  = 0, 432   Ti s   51,572s 38 + 26, 017  Hình 20 Hàm bổ sung đặc tính cứng mạch hở Hình 21 Đáp ứng hệ thống theo điều khiển PID Rotach V Y Nhận xét: - Phương pháp Rotach V Y cho phép thiết kế linh hoạt điều khiển PI PID hàm truyền đối tượng dạng qn tính bậc hai có trễ - Với đối tượng quán tính bậc 2, điều khiển Rotach V Y đưa lượng đặt xác Thiết kế điều khiển mạch tầng theo Mạnh N V Thiết kế điều khiển mạch hai vòng điều khiển nhiệt độ nhiệt NMNĐ ng Bí MR Hàm truyền đối tượng vịng lấy theo [8] Gp2 ( s ) = G p1 ( s) = 0,9424 (1 + 0,05244s) 1,923 (0,5671s +1)(2,9056s +1) e−0,4219 s e−0,8826 s Trình tự thiết kế điều khiển theo phương pháp Mạnh N V sau: 39 - Tại cửa sổ “Begin” chọn dạng mạch vòng “Cascade”, điều khiển vòng vịng ngồi dạng “PID”, chọn mơ hình đối tượng vịng vịng ngồi dạng “Second order plus dead time” Sau đó, nhập tham số hàm truyền đối tượng vịng trong, vịng ngồi trước nhấn “Next” - Tại cửa sổ “Method professor NGUYEN VAN MANH” nhập giá trị tần số cắt mong muốn trước nhấn “Enter” để thu tham số điều khiển theo đặc tính xấu hình dáng đặc tính mềm hệ hở vịng trong, vịng ngồi Tham số hai điều khiển: ( )2     0,196cG s =1  0,026s+ ( )1+ 1,518cG s =   0,105s   3,473s  1 0,474s+ + Hai điều khiển tương tự trích dẫn [18] Tuy nhiên, phần mềm cung cấp khả tính tốn tối ưu hóa nghiệm đặc tính mềm vịng trục thực cách nhấn chọn “Optimization”, từ đưa giá trị hiệu chỉnh hệ số khuếch đại điều khiển vịng ngồi giúp hai đặc tính mềm hệ hở qua điểm ( −1 j,0) Hai điều khiển thu đặc tính mềm hệ hở dẫn Hình 22 ( )1   0,026s+ 0,196cG s =   +  0,105s  ( )1   0,8088cG s =1  0,474s+ +  3,473s  Hình 22 Đặc tính mềm hệ hở vịng ngồi Nhận xét: Tham số điều khiển vịng ngồi trùng khớp với kết [8] 40 Phương pháp cho phép thiết kế điều mạch hai vòng dễ dàng với thao tác, nhiên, độ tin cậy độ trễ đối tượng vòng nhỏ nhiều đối tượng vịng ngồi Hình 23 Đáp ứng bước đầu theo lượng đặt nhiễu đối tượng dao động Nhận xét: Đáp ứng đầu hệ thống tương đồng nói khớp với phần mềm trước Thời gian dự báo 1s phát huy tác dụng 3 Hiệu chỉnh điều khiển vịng ngồi theo Mạnh N V Tiến hành thiết kế điều khiển mạch điều khiển mức nước bao NMNĐ ng Bí MR với mơ hình hàm truyền đối tượng vịng vịng ngồi sau [8]: Gp2 ( s ) = G p1 ( s) = 1, 411 + 6,956s + 120,362s 0,055 s ( 4,608s + 1)( 4,608 + 1) e−3,488s e−34,671s Với phương pháp Mạnh N V , cần xấp xỉ hàm truyền đối tượng vòng sang dạng chuẩn, thu kết sau: G p1 ( s) = 0,055 s (5,517s + 1) e−43,086 s Thực thao tác theo hướng dẫn tính tốn với phương pháp Mạnh N V thu kết sau:   + + 17,303s Gc ( s) = 1,048   6,956s  Gc1 (s) = 313(1+ 517s) 41 Hình 24 Đặc tính mềm hệ hở sau kéo dãn điều khiển vịng ngồi Tiến hành “Optimization” thu điều khiển vịng ngồi: Gc1 (s) = 2743(1+ 517s) Nhận xét: - Phần mềm tính tốn điều khiển bền vững mơ phịng đáp ứng đầu hệ thống tương tự tài liệu [8]; - Phương pháp Mạnh N V chưa thể áp dụng phần mềm với tham số điều khiển thiếu thuật tốn tối ưu hóa “Vượt khe” nên hàm truyền tương đương vịng ngồi khơng thể tối ưu hóa Vì thế, cần sử dụng phương áp thiết kế tham số điều khiển khác (ở phương án Rotach V Y ) Thiết kế điều khiển hệ hai tầng theo Rotach V Y Sử dụng phần phềm thiết kế điều khiển PI PID mạch hai vịng điều khiển phụ tải nhiệt NMNĐ Mơng Dương biến đổi mạnh [9] theo phương pháp [24] với hàm truyền đối tượng phần sở1 sau: Gp2 ( s ) = 7,832 4,873 e−0,786 s G p1 ( s) = + 1,685s (1 + 6,934s) e−3,512 s Trình tự thiết kế điều khiển vịng ngồi theo phương pháp [24]: - Thiết kế điều khiển vịng theo trình tự 3 thu tham số điều khiển PI đặc tính cứng sau:  1   1 +  1, 23654s  Gc ( s)p= k += 0, 28966   Ti s  Khi xây dựng theo phương pháp Nguyễn Văn Mạnh cần có thành phần bất định mơ hình hàm truyền 42 Hình 25 Đặc tính cứng hàm bổ sung đối tượng vòng - Chuyển sang trang “Configure PID Two-ACP”, chọn “Build a support function” để thu đồ thị hàm bổ sung với đối tượng tương đương sau có tham số điều khiển vòng - Hiệu chỉnh tham số trục đồ thị (ω, F (ω )) , giá trị “Resonance frequency” để thu điểm cực đại (đối với điều khiển PI nghiệm giao trục hoành điều khiển PID) tham số điều khiển vịng ngồi ()  1   + c pG s k= 1+  0,12699 =  Ti s   9,75482s    Hình 26 Đặc tính cứng hàm bổ sung đối tượng tương đương 43 Hình 27 Đáp ứng bước đầu hệ thống theo nhiễu Nhận xét: - Với tham số hai điều khiển nhận từ phần mềm, đặc tính cứng hệ hở hệ thống không vi phạm vùng cấm, đảm bảo độ dự trữ ổn định biên độ dao động không vượt giá trị 1,55 cho trước - Đáp ứng đầu hệ thống mạch vòng kín hệ hai tầng dễ dàng xây dựng với tiêu bình phương sai số điều khiển đạt 51,4135 đơn vị PL5 CODE CHƯƠNG TRÌNH TẠO VÀ HUẤN LUYỆN CÁC MẠNG NƠRON MISO NHẬN DẠNG QUÁ TRÌNH NHIỆT NMNĐ MISO Áp suất MISO Lưu lượng MISO Công suất điện MISO Nồng độ Oxi dư MISO Lưu lượng khói Nhiên liệu Khơng khí Nước cấp Hình 28 Mơ hình tạo mạng nơron MISO net = network (1,2, [1; 0], [1 ;0 ],[0 0; 0], [0 1] ); net layers{1} size = 200; net layers{2} size = 1; 44 % hidden layer transfer function net layers{1} transferFcn = 'logsig'; % liệu a (ví dụ x đầu vào, y đầu ra) dạng cột thêm câu lệnh x = MW'; y = Nnc' để chuyển vị x = INPUT'; y = Pbl'; view(net); %để xem cấu hình mạng net = configure(net,x,y); initial_output = net(x); net trainFcn = 'trainlm'; net performFcn = 'mse'; net trainParam epochs = 10000; net = train(net,x,y); final_output = net(y); net inputconnect = [2 ; 0]; gensim(net, 1); %tạo mạng noron simulink NNout = sim(net,[MWS1';BTHAN';]) % Tính đầu Nơron theo đầu vào PL6 CODE CHƯƠNG TRÌNH TẠO VÀ HUẤN LUYỆN CÁC MẠNG NƠRON NHẬN DẠNG VÀ HIỆU CHỈNH LƯỢNG ĐẶT HỆ THỐNG ĐIỀU KHIỂN CẤP NHIÊN LIỆU Ne(t) Wf(t) NN αNN(t) Pq n(t) Hình 29 Mơ hình tạo mạng nơ ron đầu vào đầu net = network (1,2, [1; 0], [1 ;0], [0 0; 0], [0 1]); net layers{1} size = 50; net layers{2} size = 50; % hidden layer transfer function net layers{1} transferFcn = 'logsig'; x = BTUIMPUT'; y = Alpha'; view(net); %để xem cấu hình mạng net = configure(net,x,y); initial_output = net(x); net trainFcn = 'trainlm'; 45 net performFcn = 'mse'; net trainParam epochs = 100000; net = train(net,x,y); final_output = net(y); net inputconnect = [2 ; 0]; gensim(net, 250); %tạo mạng noron simulink NNout = sim(net,[MWS1';BTHAN';Pqn’;]) % Tính đầu Nơron theo đầu vào PL7 ỨNG DỤNG MẠNG NƠRON NHẬN DẠNG ẢNH HƯỞNG DAMPER TỚI ÁP SUẤT BUỒNG LỬA CỦA HỆ ĐIỀU KHIỂN NMNĐ ĐỐT THAN Hiện nay, lò vận hành thường khơng có thiết bị đo, truyền dẫn trực tuyến nhiệt độ tâm buồng đốt hệ thống DCS Tâm lửa thường đánh giá thông qua dụng cụ cầm tay lần sửa chữa Trong trình vận hành, đại lượng sử dụng để nhận định tương đối tâm lửa nhiệt độ cuối buồng đốt, đầu vào nhiệt cấp Do để đánh giá ảnh hưởng damper quạt khói tới áp suất buồng đốt, luận án nghiên cứu ba mơ hình nhận dạng minh họa Hình 30 %DID_A(t) NN Pbl(t) NN Pbl(t) a) %DID_A(t) b) %DID_B(t) %DID_A(t) %DID_B(t) c) NN Pbl(t) Wair2(t) Hình 30 Mơ hình nhận dạng áp suất buồng đốt a) Một tín hiệu đầu vào; b) Hai tín hiệu đầu vào; c)Ba tín hiệu đầu vào Cả ba mơ hình thực có điểm chung có tín hiệu đầu áp suất buồng đốt Điểm khác nằm số tín hiệu đầu vào: 46 - Phương án 1: Mơ hình thứ với đại lượng đầu vào độ mở damper cánh tĩnh quạt khói IDF A; - Phương án 2: Mơ hình thứ hai với hai tín hiệu độ mở quạt khói IDF A B; - Phương án 3: Mơ hình thứ bổ sung tín hiệu lưu lượng gió cấp hai với vai trị nhiễu q trình a) b) c) Hình 31 Dữ liệu phục vụ huấn luyện mạng nơron nhận dạng áp suất buồng lửa a) Độ mở damper quạt khói A,B; b) Lưu lượng gió 2; c) Áp suất buồng lửa 47 Kịch mô phỏng, phương pháp số liệu sử dụng tương thích với mục Sử dụng mạng nơron có hai lớp ẩn (lớp thứ gồm 100 nơron, lớp thứ có nơron) Dữ liệu phục vụ huấn luyện mạng nơron trích suất từ tập liệu vận hành tổ máy số Duyên Hải ngày 30/7/2018 (Phụ lục 1) mơ Hình 32 Bao gồm: lưu lượng gió cấp cấp vào lị hơi, độ mở damper cánh tính đầu hút quạt khói số 1,2 áp suất buồng lửa Kết trả mạng nơron, giá trị áp suất chân không buồng lửa thực tế minh họa Hình 32 a) b) c) Hình 32 Kết mơ hình mạng nơron theo phương án a) Phương án 1; b) Phương án 2; c) Phương án 48 Bảng tổng hợp kết huấn luyện mạng nơron: STT Nơron Performance Gradient Mu Phương án 874 8,46 e-6 0,01 Phương án 871 2,20 e-7 100 Phương án 863 12,8 e-7 1000 Kết so sánh sai lệch đầu mạng nơron giá trị lượng đặt -50Pa mơ tả Hình 33 Hình 33 Kết sai lệch đầu mạng nơron lượng đặt Nhận xét: - Khi thay đổi bổ sung tín hiệu đầu vào, sai lệch đầu mạng nơron với giá trị đặt ngày lớn Đầu mạng nơron dao động qua giá trị đặt mạch vòng điều khiển áp suất buồng lửa -50 Pa Xu hướng trả mơ hình mạng nơron tuân theo xu hướng trình cháy - Với tác động bổ sung nhiễu trình lưu lượng gió cấp 2, áp suất buồng đốt ngày lệch khỏi giá trị đặt mạch vòng Đặt yêu cầu cho điều khiển tác động lại độ mở damper, dẫn tới tín hiệu áp suất buồng đốt có dao động lớn đảm bảo bao quanh giá trị đặt - Sử dụng mạng nơron đánh giá tình trạng vận hành thiết bị thực nguyên tắc sau: + So sánh có sai lệch độ mở damper quạt khói A, B dẫn tới áp suất buồng đốt trả mạng nơron thay đổi bất thường so với vận hành đưa cảnh báo hai thiết bị xảy cố + Sử dụng tập liệu mẫu từ lịch sử vận hành phục vụ huấn luyện mạng nơron để nhận dạng đại lượng đặc trưng phản ánh tính trạng thiết bị thơng qua đại lượng đầu vào tác động tới chế độ vận hành thiết bị, gọi giai 49 đoạn học Thời gian học lâu liệu phục vụ huấn luyện nhiều thể tính xác mơ hình mạng nơron, đồng thời giảm sai số nhiễu không mong muốn + Giai đoạn đánh giá & dự báo xác định theo nguyên tắc sử dụng liệu đầu vào thực tế để kiểm chứng đầu mạng, sau đó, so sánh với nguồn sở liệu mẫu (được xác định thời điểm vận hành ổn định) để nhận biết tình trạng thiết bị vận hành dự kiến thời gian xảy cố 50 ... BẢN NHÀ MÁY NHIỆT ĐIỆN ĐỐT THAN PHUN 111 Lò nhà máy nhiệt điện đốt than phun 112 Tuabin nhà máy nhiệt điện 12 NGUYÊN NHÂN DẪN ĐẾN CẦN PHẢI HIỆU CHỈNH VÀ CHỈNH ĐỊNH CÁC BỘ ĐIỀU KHIỂN TRONG NHÀ MÁY...BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI PHƯƠNG PHÁP NHẬN DẠNG VÀ HIỆU CHỈNH CÁC MẠCH VÒNG ĐIỀU KHIỂN NHÀ MÁY NHIỆT ĐIỆN ĐỐT THAN PHUN Ngành: Kỹ thuật nhiệt Mã số: 9520115... phương án nhận dạng, hiệu chỉnh giúp trì ổn định chất lượng điều khiển nhà máy nhiệt điện Do vậy, nghiên cứu sinh lựa chọn đề tài nghiên cứu ? ?Phương pháp nhận dạng hiệu chỉnh mạch vòng điều khiển